Search
E-mail
RSS

Online submission Editor-in-chief Peer review Staff Login

Journal of Tropical Oceanography >

2021 , Vol. 40 >Issue 4: 110 - 121

DOI: https://doi.org/10.11978/2020093

Marine Geology

Interpreting the origin of coastal boulders on a coral reef flat at Xiaodonghai of Hainan Island based on storm wave energy analysis

  • WANG Yuanqi , 1 ,
  • YANG Yang 2 ,
  • ZHOU Liang 3 ,
  • WANG Yaping 1 ,
  • GAO Shu , 1
Expand
  • 1. State Key Laboratory of Estuarine and Coastal Research, East China Normal University, Shanghai 200241, China
  • 2. School of Marine Science and Engineering, Nanjing Normal University, Nanjing 210046, China
  • 3. School of Geography, Geomatics and Planning, Jiangsu Normal University, Xuzhou 221116, China
GAO Shu. email:

Copy editor: YAO Yantao

Received date: 2020-08-22

  Revised date: 2020-11-17

  Online published: 2020-12-01

Supported by

National Natural Science Foundation of China(41530962)

National Natural Science Foundation of China(41706095)

National Natural Science Foundation of China(41706096)

Copyright

Copyright reserved © 2021. Office of Acta Agronomica Sinica All articles published represent the opinions of the authors, and do not reflect the official policy of the Chinese Medical Association or the Editorial Board, unless this is clearly specified.
Fold

Abstract

Storm and tsunamis boulders may be distinguished by analyzing wave energy dissipation patterns. Here, we use the principles of storm energy conservation and the transport by storm-induced flow, to establish a theoretical distribution pattern of storm boulders on a coral reef flat. The movements of coral boulders with different sizes during a storm event are quantified, to derive their transport distances. The method is applied to the Xiaodonghai coral reef flat of Hainan Island, in combination with the information on historical typhoon events. The result is generally in agreement with the in situ records of the boulders, indicating that the deposition of the boulders is controlled by the storm energy level and the physical properties of the boulders, i.e., the boulder deposits at Xiaodonghai are caused by storm events. In the future, a combined analytical procedure may be established if the storm- and tsunami-induced transport theories are further improved, with improved assumptions for modelling.

Key words： storm boulders; wave energy distribution; storm surge; material transport; criteria for boulder deposit; Hainan Island

Cite this article

WANG Yuanqi , YANG Yang , ZHOU Liang , WANG Yaping , GAO Shu . Interpreting the origin of coastal boulders on a coral reef flat at Xiaodonghai of Hainan Island based on storm wave energy analysis[J]. Journal of Tropical Oceanography, 2021 , 40(4) : 110 -121 . DOI: 10.11978/2020093

*感谢华东师范大学河口海岸学国家重点实验室钱伟伟工程师、南京大学地理与海洋科学学院王丹丹、许眸莹、杨保明、樊一阳、祁雅莉参加了野外调查工作; 感谢刘桢峤对样品进行了室内实验分析。
海岸珊瑚巨砾是高能极端水文事件古记录最有价值的标志之一(Hayne et al, 2001; Scheffers et al, 2008; Yu et al, 2009; Terry et al, 2014, 2016)。海岸巨砾有风暴和海啸两种成因, 20世纪中期的研究普遍认为海岸珊瑚巨砾是由风暴潮及其伴随的大浪造成的(Mastronuzzi et al, 2004)。2004年印度洋海啸事件之后,人们开始更多地关注于其海啸成因(Scicchitano et al, 2007; Goff et al, 2010; Costa et al, 2011)。
风暴与海啸巨砾沉积具有很大的相似性。通常, 在平行岸线方向形成脊状分布的巨砾堆积一般考虑为风暴成因(Etienne et al, 2010), 在垂直岸线方向不规则分布的巨砾堆积则大多考虑为海啸成因(Paris et al, 2007; Etienne et al, 2011; Goto et al, 2012)。然而, 由于海岸珊瑚巨砾来源的不确定性以及海岸地形环境的复杂性(Nott, 1997; Goto et al, 2010b; Weiss, 2012), 对这一判据存在争议。在同时受到海啸和风暴潮事件影响的海岸, 如澳大利亚新南威尔士海岸(Bryant et al, 2001; Felton et al, 2003; Saintilan et al, 2005; Goff et al, 2009)、意大利西西里岛东南海岸(Scicchitano et al, 2007; Barbano et al, 2010)、地中海马耳他岛海岸(Biolchi et al, 2016; Causon et al, 2017), 巨砾的成因区分难题在近年来的海岸地貌研究中占据了重要位置。
目前对于风暴、海啸巨砾的区分研究大多从巨砾堆积年龄以及沉积动力角度入手。通过对珊瑚巨砾进行C14或U238测年, 并与已知的极端事件发生时间进行对比, 可大致判定巨砾的成因(Hayne et al, 2001; Terry et al, 2014)。但礁坪上的珊瑚巨砾因风化、波浪作用和生物侵蚀过程而发生剥蚀, 珊瑚死亡年龄也可能存在不确定性(Goto et al, 2007), 故该方法只能提供极端事件的时间上限(Yu et al, 2009)。对于巨砾搬运的水动力过程, Nott(2003)提出了不同初始状态下巨砾绕其长轴发生转动时的最小波高, 以推算极端事件的能量, 但该公式对巨砾运动所需起动波高的计算值与真实值相差较大(Mastronuzzi et al, 2004)。Lorang(2000, 2011)借鉴防波堤护面抛石的起动阈值, 推算了与海岸巨砾沉积位置相关的波浪周期, 并计算了巨砾堆积的极端事件能量, 但该方法只适用于海滩坡面上的巨砾堆积, 对于坡面之外的珊瑚礁平台环境并不适用。Scheffers等(2003)提出了“移动指数”概念对两类巨砾进行区分, 根据巨砾质量、沉积位置与岸线之间的距离、相对于平均海面的垂直高度三者的乘积近似地判断极端事件的能量。同样, 该方法也不适用于珊瑚礁平台上的巨砾沉积。Goto等(2010a, 2011)通过对冲绳岛、石垣岛已知风暴和海啸成因的巨砾进行研究, 得出西太平洋地区风暴、海啸巨砾的极限搬运距离, 获得了区分珊瑚巨砾归属成因的大致距离标准。
本文选取我国海南岛三亚小东海为研究区, 该地兼受风暴潮和海啸影响, 是珊瑚礁地形风暴、海啸巨砾成因区分研究的良好环境。从动力学角度出发, 考虑海岸珊瑚巨砾在风暴潮作用下的运动过程, 基于波浪能量分析, 对小东海珊瑚礁海岸的珊瑚巨砾沉积的理论空间分布格局进行刻画。通过对比野外实际观测的珊瑚巨砾分布与风暴巨砾的理论空间分布, 从而对巨砾成因进行判断。本项研究的目的是将巨砾输运过程与堆积位置分布格局相结合, 提出区分风暴、海啸巨砾的新思路。

1 研究区概况

海南岛三亚鹿回头半岛三面环海, 其南端及西海岸一带发育有众多岸礁, 小东海位置如图1所示。小东海西南沿岸珊瑚礁平台宽度约为200~300m, 珊瑚礁地形以-1.3m等深线为界, 内侧为礁坪带, 面积约0.27km2, 外侧为礁缘带, 面积约0.15km2(毛龙江 等, 2006)。礁坪内侧在低潮时暴露, 礁坪外侧通常处于淹没状态。
图1 海南岛南部小东海研究区地理位置

图a底图源自中华人民共和国自然资源部海南测绘地理信息局标准地图服务网站[审图号: 琼S(2019)063], 图b底图来自Google Earth

Fig. 1 Study area of Xiaodonghai on the southern coast of Hainan Island. The base map is from Google Earth

本区珊瑚礁主体发育于7000~6300 a BP全新世高海面时期(黄德银 等, 2005)。珊瑚礁坪内侧以蔷薇珊瑚、鹿角珊瑚为主; 礁坪外侧以角蜂巢珊瑚、扁脑珊瑚和鹿角珊瑚为主; 礁缘带斜坡方向活体珊瑚分布密集, 水下2~5m范围内发育有大面积的角蜂巢珊瑚、扁脑珊瑚、鹿角珊瑚、杯形珊瑚和沙珊瑚等, 表面覆盖率可达到50%~65%(王月 等, 2011)。
小东海位于热带, 受海洋性季风气候影响, 夏季盛行西南风。海域受潮汐和风浪影响, 潮汐潮差小, 日潮不等现象明显。全年平均波高为1.3m, 平均波周期为5.1s, 常浪向为E向。夏季平均波高较低, 月平均最小波高为0.9m; 冬季平均波高较大, 月平均最大波高1.9m。浪向季节变化显著, 9月底至次年3月中旬常浪向为E向, 6月中旬至8月中旬常浪向为SSE向(黄德银 等, 2005)。
海南岛是我国遭受风暴潮影响最严重的地区之一。冬季受冷空气影响, 灾害性海浪发生次数较多, 而夏季由于热带气旋影响, 灾害性海浪程度更为严重(李淑江 等, 2016)。近年来对海南省产生严重影响的灾害性海浪大多都是由热带气旋登陆造成。海南岛遭受海啸侵袭的概率较小, 但历史上海南岛确有因近海地震而引起的海啸事件(Sun et al, 2013)。

2 材料与方法

2.1 样品采集与实验室分析

2018年1月31至2月3日对小东海西南沿岸珊瑚礁坪开展野外调查(图2)。挑选礁平台上18块典型的海岸珊瑚巨砾, 使用测量精度为15mm的Riegl VZ4000地面3D激光扫描仪进行扫描, 获取其实体表面信息。采集完成后通过 3D 激光扫描仪处理软件 RiSCAN PRO 构建“点云”模型, 并导入软件Cloud Compare中, 计算出巨砾的三轴长度与体积。采用 GPS-RTK技术准确定位这些海岸珊瑚巨砾的位置, 将巨砾位置坐标输入到Google Earth中,测算巨砾与珊瑚礁坪前缘的水平距离。根据阿基米德排水法的巨砾密度测量结果(刘桢峤 等, 2019), 取平均值作为研究区珊瑚巨砾的密度。
图2 小东海西南沿岸珊瑚礁坪

a. 18块典型巨砾(XD01~XD18)的分布, 图中各巨砾的编号以01~18标示; b. 小东海西南沿岸珊瑚礁坪卫星图, 底图来自Google Earth; c. 编号为XD01的巨砾特征

Fig. 2 Distribution of 18 typical boulders on the coral reef flat at Xiaodonghai (a); satellite image of the coral reef flat (with the base image from Google Earth) (b); and (c) characteristics of No. XD01 boulder

2.2 风暴潮作用下的珊瑚巨砾移动距离计算

在前人对珊瑚巨砾起动流速和波高计算(Nott, 1997, 2003; Nandasena et al, 2011)的基础上, 可进一步拓展巨砾整体搬运距离的计算方法。当巨砾在风暴潮水流作用下以推移质方式搬运时, 其整体移动过程可分为两部分。当水流速度达到巨砾运动的起动流速后, 巨砾在水流拖曳力Fd和动摩擦力f的合力作用下从静止开始向礁坪内侧推移, 如图3a所示的A点位置。此后, 由于珊瑚礁坪上能量不断耗散, 作用于巨砾上的拖曳力Fd不断减小, 巨砾在此阶段作变加速运动。巨砾运动速度不断增大, 当运动至B点位置时达到最大值, 巨砾所受合力为0。在AB段运动过程中, 巨砾水平方向受力如图3b所示, 其中水流拖曳力Fd以及与床面之间的动摩擦力f可表达为(Nott, 2003):
${{F}_{d}}=\frac{1}{2}\rho {{C}_{d}}(ac)u_{i}^{2}$
${{f}_{{}}}=\mu \left( \frac{{{\rho }_{s}}-\rho }{{{\rho }_{s}}} \right)mg$
式中: ρ为海水密度, 取1.02g·cm-3; Cd为拖拽系数, 根据Nott (2003)的研究取Cd=2; a为巨砾长轴长度, 单位为m; c为巨砾短轴长度, 单位为m; μi为水流近底流速, 单位为m·s-1; ρs为巨砾密度, 单位为g·cm-3; m为巨砾质量, 单位为kg; g为重力加速度, 取9.81N·kg-1; µ为摩擦系数, 根据Bagnold(1963)的研究, 摩擦系数与巨砾大小和床面粗糙程度相关, 取值范围在0.32~0.75之间, 本文根据Hongo等(2018)的实验结果取0.73。
图3 风暴潮作用下珊瑚巨砾在礁坪上的运动过程

a. 珊瑚巨砾运动过程示意图; b. AB段运动过程中巨砾在水平方向上的受力; c. BC段运动过程中巨砾在水平方向上的受力

Fig. 3 The movements of coral boulders on the platform under the action of storm surge: (a) the movement from A to B and finally to C; (b) the horizontal forces exerted on the boulder in the segment AB; and (c) the horizontal forces in the segment BC

破波带内水流近底流速μi与沿程各点的波高H的关系近似为(Nott, 2003):
${{\mu }_{i}}=\sqrt{\text{g}H}$
根据牛顿第二定律推导, 可得巨砾在AB段运动的加速度a(单位: m·s-2)为:
$a=\frac{{{F}_{d}}-{{f}_{{}}}}{m}$
加速度的方向与巨砾运动方向一致, 指向礁坪内侧。
假定经过时间T1后, 巨砾可由A点运动至B点, 根据物体运动速度公式可推导巨砾在B点的运动速度vb为:
${{v}_{b}}=\int_{{{T}_{1}}}^{0}{a\text{d}t=\bar{a}{{T}_{1}}}$
$\bar{a}=\frac{{{F}_{d}}-2{{f}_{{}}}}{2m}$
两公式中的$\bar{a}$为巨砾在AB段的平均加速度(单位: m·s-2), 方向与巨砾运动方向一致, 指向礁坪内侧。
根据物体运动的速度位移公式可推导出巨砾在T1时间内的运动距离L1(即AB段的长度)为:
${{L}_{1}}=\frac{v_{\text{b}}^{2}}{2\bar{a}}$
从功能转换关系角度考虑, 波浪破碎后能量释放, 其中部分能量对巨砾做功。波成水流对巨砾做功, 不仅使得巨砾从A点搬运至B点, 还使得巨砾从静止状态加速至运动速度vb。假定波浪做功的大小为E1, 根据能量守恒定律, 巨砾在运动过程中的能量关系如下:
${{E}_{1}}=m\bar{a}{{L}_{1}}+\frac{1}{2}mv_{\text{b}}^{2}$
Power等(2013)根据能量通量平衡得出破波带内单位面积的波浪能量大小与波高的平方以及波群速度成正比。根据Nielsen (1992)的研究, 在波浪破碎所消耗的波浪能量中, 转化为泥沙运动的效率系数eB的量级为10-2, 故本文推导波浪对巨砾所做的功大小系数E1可近似地由下式确定:
${{E}_{1}}=\frac{1}{8}{{e}_{B}}\rho \text{g}\sqrt{\text{g}h}{{H}^{2}}S$
式中: 效率系数eB取3.5×10-2; H为波浪波高(单位: m); h为水深(单位: m); S为波浪做功的截面面积, 取巨砾长轴与短轴截面面积的10倍(单位: m2)。根据上述运动学方程及能量关系可计算出巨砾在AB段做变加速运动的距离L1。
当巨砾运动至B点时, 其速度不为零, 且将继续朝礁坪内侧推移, 直至最后静止而沉积在C点(图3a)。巨砾在BC段的水平方向受力如图3c所示, 其所受合力主要以床面动摩擦力f为主, 根据牛顿第二定律可推导巨砾在BC段的加速度为:
$a=\frac{{{f}_{{}}}}{m}\text{=}\mu \left( \frac{{{\rho }_{s}}-\rho }{{{\rho }_{s}}} \right)g$
加速度方向与巨砾运动方向相反, 指向礁坪外侧。从功能转换关系角度考虑, 巨砾从B点运动至C点过程中, 其动能将全部转化为巨砾位移L2:
$\frac{1}{2}mv_{\text{b}}^{2}=m\bar{a}{{L}_{2}}$
$\bar{a}=\frac{{{f}_{{}}}}{2m}$
两公式中的$\bar{a}$为巨砾在BC段的平均加速度(单位: m·s-2), 方向与巨砾运动方向相反, 指向礁坪外侧。
由上式推导可得, 巨砾在BC段的运动距离L2为:
${{L}_{2}}=\frac{mv_{\text{b}}^{2}}{{{f}_{{}}}}$
由此可得巨砾在风暴潮作用下的总搬运距离L, 即为AB段距离L1与BC段距离L2之和。

3 结果

3.1 野外实测巨砾参数

研究区礁坪有大量海岸珊瑚巨砾分布, 堆积形态呈现出叠瓦状, 向海倾斜。巨砾形状以近似长方体和椭球体为主。构成珊瑚巨砾沉积的珊瑚碎屑来源于枇杷珊瑚科、褶叶珊瑚科、滨珊瑚科和蜂巢珊瑚科的珊瑚石(练健生 等, 2010)。巨砾表面在高潮水位线以上附着有海洋生物, 可推断礁坪上的珊瑚巨砾沉积是在高能极端事件中从礁坪前缘搬运形成。
根据实验室测定结果(刘桢峤 等, 2019), 小东海西南沿岸礁平台分布的珊瑚巨砾的平均密度为1.82g·cm-3。所测量的珊瑚巨砾长轴(a)长度在1.63~3.47m范围内, 中轴(b)长度在1.34~3.02m范围内, 短轴(c)长度在1.17~2.11m范围内。根据Corey形状因子$S=c/\sqrt{ab}$计算野外观测巨砾的形状参数(Corey, 1949), 其形状因子在0.4~0.8之间, 平均为0.6。18块巨砾中, 最大的珊瑚巨砾体积为10.89m3, 最小为2.19m3。由此估算, 巨砾最大质量可达19820kg, 最小质量为4350kg, 平均质量为8275kg。珊瑚巨砾与礁坪前缘距离在3.18~116.0m范围内, 平均距离为53.4m (表1)。
表1 小东海西南沿岸礁平台巨砾参数

Tab. 1 Boulder parameters of a reef flat in the southwest coast of Xiaodonghai

巨砾编号 长轴长度/m 中轴长度/m 短轴长度/m Corey形状因子 质量/kg 巨砾位置与礁坪前缘的距离/m
XD01 2.94 2.58 2.11 0.8 19820 64.6
XD02 1.83 1.82 1.18 0.6 4914 53.8
XD03 2.15 1.72 1.19 0.6 6370 46.5
XD04 1.63 1.34 1.21 0.8 3986 47.6
XD05 1.88 1.67 1.22 0.7 5587 40.2
XD06 2.53 2.44 1.67 0.7 8627 31.8
XD07 2.29 1.85 1.19 0.6 6625 49.2
XD08 2.48 2.30 1.51 0.6 6334 43.1
XD09 3.47 1.93 1.39 0.5 6479 57.1
XD010 2.15 1.91 1.27 0.6 4350 70.6
XD011 2.83 2.46 1.64 0.6 9335 32.6
XD012 2.48 2.46 1.61 0.6 8427 37.1
XD013 3.12 2.47 1.17 0.4 9682 39.4
XD014 2.61 2.25 1.54 0.6 7225 57.1
XD015 3.34 2.23 1.48 0.5 9045 71.9
XD016 2.72 2.61 1.61 0.6 10683 116.0
XD017 3.10 2.22 1.31 0.5 7680 84.2
XD018 3.16 3.02 1.62 0.5 13777 72.0

3.2 小东海风暴巨砾的理论搬运距离

目前的海岸巨砾成因区分研究主要是通过巨砾的形状大小来推断其起动时所需的极端事件能量, 从而进行风暴巨砾或海啸巨砾的成因判定。本文聚焦于风暴潮波浪能量作用下珊瑚巨砾的整体搬运动力过程, 对小东海巨砾在不同风暴条件下的理论搬运距离L进行计算, 以确定小东海风暴巨砾的理论空间分布格局。
基于SWAN模型对1949—2013年台风级别的471次热带气旋的台风浪模拟结果, 可确定风暴潮的能量。南海北部地区1949—2013年间台风风暴潮引起的海浪最大有效波高为13.0m(尹洪强 等, 2014), 本文将其作为小东海风暴巨砾搬运波高的上限。根据莺歌海和东方所1964—2005年的观测资料, 海南岛由热带气旋引起的风暴潮海浪平均波高为5.9m, 故将风暴巨砾搬运所需海浪波高的下限定为6.0m。
不同形状的巨砾在风暴潮条件下的理论搬运距离有着较大差异。假定巨砾为椭球体, 其体积按照标准的椭球体计算公式$V=\frac{1}{6}\pi abc$确定。海岸珊瑚巨砾的质量范围介于0.5×103~2.0×104kg, Corey形状因子分别为0.4、0.6、0.8。将不同质量和形状的巨砾参数代入公式(1)、(2)及(6)中, 可得巨砾在AB段运动的平均加速度$\bar{a}$, 结合公式(7)、(8)、(9)可推得巨砾在AB段做变加速运动的距离L1。根据公式(12)和(13)可计算出巨砾在BC段的运动距离L2, 进而获得风暴潮波浪能量范围内巨砾从礁坪起动后的理论搬运距离L(表2)。将巨砾参数、风暴潮波浪波高和相应的理论搬运距离L绘制成图(图4)。
表2 风暴巨砾在礁坪上起动后的理论搬运距离

Tab. 2 Theoretical transport distances of storm boulders after mobilization on the reef flat

Corey形状因子 巨砾质量/(×103kg) 波高H=13m时的理论搬运距离/m 波高H=6m时的理论搬运距离/m
0.4 0.1 277.8 60.0
1 127.6 28.1
5 76.4 17.4
10 61.7 14.3
15 54.4 12.8
20 49.9 12.0
0.6 0.1 285.6 61.7
1 136.2 30.0
5 81.1 18.3
10 65.5 15.1
15 57.6 13.5
20 52.8 12.5
0.8 0.1 299.2 64.6
1 139.4 30.6
5 83.7 18.9
10 67.2 15.5
15 59.3 13.8
20 54.3 12.8
图4 风暴巨砾在礁坪上起动后的理论搬运距离

图中阴影区为小东海风暴巨砾从礁坪上起动的理论分布空间, 黑色标记点为Corey形状因子分别为0.4、0.6、0.8的野外观测巨砾的分布位置

Fig. 4 Theoretical transport distances of storm boulders: (a) Corey shape factor S=0.4; (b) S=0.6; and (c) S=0.8

由图4和表2可知, Corey形状因子S=0.4的风暴巨砾的理论搬运距离为12.0~277.8m, S=0.6的风暴巨砾为12.5~285.6m, 而S=0.8的风暴巨砾为12.8~299.2m; 小东海风暴巨砾理论搬运距离的最大平均值为287.5m, 最小平均值为12.4m。

4 讨论

4.1 风暴巨砾理论分布格局

台风和海啸都具有强大的搬运能力, 能使巨砾在海岸线附近形成堆积, 但两者的作用机制存在差异。海啸一般只经历一次来回或只有向岸输移的过程, 而台风过程影响时间较长, 往往存在反复多次的向岸运移和回流过程(高抒 等, 2019)。海啸可将重量达0.1×103~0.2×103kg的巨砾搬运500m以上的距离(Nandasena et al, 2013), 而台风引起的风暴潮将巨砾搬运至岸线前缘的距离一般不超过300m (Goto et al, 2010a)。由此可知, 风暴潮所蕴含的能量与搬运的巨砾质量及移动距离存在相应的关系。
当巨砾在礁坪上起动时, 风暴潮的波高大小对小东海巨砾整体的搬运距离有着重要影响。Corey形状因子为0.4、质量为0.1t的巨砾在波高分别为13m和6m的风暴潮作用下的理论搬运距离相差217.8m, 而质量为20t的巨砾则相差37.9m; Corey形状因子为0.8、质量为0.1t的巨砾在波高分别为13m和6m的风暴潮作用下的理论搬运距离相差234.6m, 而质量为20t的巨砾则相差41.5m。这表明质量越小的巨砾, 在风暴潮作用下的搬运距离越远, 且风暴潮波高的变化对其理论搬运距离存在明显的影响, 形成的风暴巨砾空间分布越广。随着巨砾质量的增大, 风暴潮波高的变化对理论搬运距离的影响减弱, 巨砾的理论空间分布范围相对缩小。
对于相同质量的巨砾, 形状的不同对其搬运距离也有明显的影响。由图5可知, Corey形状因子越大, 同等极端事件影响下, 巨砾搬运距离也会越远。当波高为13m时, 质量为0.1t、Corey形状因子分别为0.4与0.8的巨砾在搬运距离上相差21.4m, 而质量为20t的巨砾则相差4.4m; 当波高为6m时, 质量为0.1t、Corey形状因子分别为0.4与0.8的巨砾在搬运距离上相差4.6m, 而质量为20t的则相差0.8m。Corey形状因子为0.8的风暴巨砾在理论空间分布上较Corey形状因子为0.4的巨砾更远离珊瑚礁坪前缘, 且分布范围更广。
图5 不同形状的风暴巨砾在礁坪上起动后的理论搬运距离差异

Fig. 5 Differences of theoretical transport distance of storm boulders due to different shape factors

4.2 小东海巨砾的风暴成因

在目前的观测研究中, 对不同质量的风暴巨砾在礁坪上的理论空间分布格局进行系统描述的较少。根据珊瑚礁海岸上风暴巨砾的理论空间分布格局, 若野外实测巨砾在其理论空间分布范围内, 可判定该巨砾为风暴巨砾。本文以小东海珊瑚礁海岸为例, 将现场观测的巨砾质量、形状参数以及巨砾至珊瑚礁前缘的实际距离与小东海风暴巨砾的理论空间分布格局进行对比, 并结合区域极端事件记录, 对小东海的巨砾进行成因判别。现场实测巨砾在风暴潮情况下的理论分布范围用理论搬运距离区间表示(表3)。
表3 小东海野外实测巨砾在风暴潮作用下的理论空间分布范围与实测距离

Tab. 3 Theoretical spatial distribution and measured distances of boulders at Xiaodonghai

巨砾编号 最小理论搬运距离/m 最大理论搬运距离/m 巨砾与礁坪前缘的实测距离/m
XD01 12.9 54.4 64.6
XD02 18.5 82.3 53.8
XD03 17.2 75.5 46.5
XD04 20.2 89.9 47.6
XD05 18.0 78.7 40.2
XD06 15.9 69.5 31.8
XD07 17.0 74.8 49.2
XD08 17.2 75.8 43.1
XD09 16.6 72.8 57.1
XD010 19.1 84.8 70.6
XD011 15.4 66.9 32.6
XD012 16.0 70.0 37.1
XD013 14.5 62.3 39.4
XD014 16.5 74.5 57.1
XD015 15.1 65.6 71.9
XD016 14.8 64.1 116.0
XD017 16.2 68.8 84.2
XD018 13.5 57.5 72.0
由表3结果可知, XD01、XD015、XD016、XD017、XD018这5块巨砾的理论搬运距离与实测距离稍有偏离, 其余巨砾均位于小东海风暴巨砾的理论分布范围内。对于有所偏离的巨砾, 其沉积位置至礁坪前缘的实际距离超出了13m波高作用下所计算得出的理论搬运距离。
考虑到海南岛风暴潮事件多发, 珊瑚礁巨砾与礁坪前缘的实测距离可能是巨砾在风暴潮作用下多次移动的结果。根据遥感影像对比, 在2013年超强台风“海燕”引起的风暴潮作用下, 原本位于菲律宾萨马尔岛东部和帕劳岛的珊瑚礁坪上的部分巨砾发生了移动, 使巨砾与礁坪前缘之间的距离增大(May et al, 2015; Hongo et al, 2018)。因此, 推测上述5块巨砾目前所处的位置可能是经过多次风暴潮搬运的结果。
另一原因可能是巨砾曾受到了波高大于13m的极端事件搬运。通过波浪能量分析, 反推XD01、XD015、XD016、XD017、XD018被搬运至实测距离所需的风暴潮海浪波高, 结果如表4所示。
表4 巨砾被搬运至实测距离所需的风暴潮海浪波高

Tab. 4 Storm surge wave height required for boulder transport to measured distance

巨砾编号 巨砾与礁坪前缘的实测距离/m 风暴潮海浪波高/m
XD01 64.6 14.2
XD015 71.9 13.5
XD016 116.0 17.8
XD017 84.2 14.3
XD018 72.0 14.6
根据统计资料显示, 海南岛近70年来出现大于2m增水的风暴潮共有7次, 其中台风“威马逊”和“海鸥”引发了近年来少有的两次特大风暴潮灾害(陈哲, 2017)。根据莺歌海和东方所1964—2005年的观测资料, 海南岛由热带气旋引起的风暴潮海浪最大波高可达9m以上。而利用WAVEWATCH III海浪数值模式对南海1986—2005年发生的275次台风的数值模拟结果表明, 台风作用下海南岛东北部海域百年一遇的极值波高可达16m(苗庆生, 2011)。2005—2008年内, 海南岛因热带风暴引起的灾害性海浪实测波高均小于8m, 2008年“黑格比”强台风引起的海浪实测最大波高则在10m以上(李文欢 等, 2013)。
除了历史风暴浪影响之外, 海南岛也曾受到海啸的侵袭。根据Sun等(2013)的研究, 中国南海在公元1076年曾发生过海啸, 西沙群岛东岛以及广东南澳岛等均受到海啸影响。历史文献中海南省遭受海啸影响的记录有: 1906年琼山发生7.5级地震并引发海啸(Mak et al, 2007); 1992年在海南三亚也发生过海啸事件, 但波能较小, 最大波高仅为0.8m(叶琳, 2005); 2004年12月和2005年3月因受印度洋地震影响, 海南三亚也观测到了海啸波浪, 但海啸引起的最大波高在0.1m以下(Mak et al, 2007)。近百年来, 海南岛海啸事件发生的频率远低于风暴频率。
结合小东海极端事件历史记录可知, 巨砾XD01、XD015、XD017、XD018由风暴潮搬运而成的可能性较大, 而搬运XD016所需的风暴潮海浪波高应大于16m, 因而推测其成因可能是历史强风暴或海啸。在后期研究中还需对珊瑚巨砾进行年代测定, 并与已知台风或海啸记录作对比, 以进一步确定小东海珊瑚巨砾的成因, 同时判断台风或者海啸的发生频率以及强度变化。
此外, 小东海野外观测巨砾中处在风暴巨砾理论分布空间外的5块巨砾的质量普遍在1.0×104kg以上, 属于大型的珊瑚巨砾块体。Nandasena等 (2011)认为大型珊瑚巨砾与一般的珊瑚巨砾在起动流速上存在较大差异, 大型珊瑚巨砾需要更大的起动流速; 同时, 大型珊瑚巨砾在搬运过程中与礁坪间的接触面积更大, 单位面积所受摩擦力与一般巨砾相比也可能存在差异。因此, 对于质量大于1.0×104kg的珊瑚巨砾的成因判断方法, 仍需进一步研究。

4.3 风暴巨砾理论空间分布作为成因区分方法的适用性

小东海沿岸分布的珊瑚巨砾在台风作用下起动搬运所需的海浪波高应在7.2~15.8m之间(刘桢峤 等, 2019)。本文应用风暴巨砾的理论空间分布对小东海沿岸的巨砾成因进行判定。根据前述研究(表4)可知, 本文所测巨砾中除XD016外, 其余巨砾被搬运至现今位置所需的风暴潮波高均小于或等于14.6m, 这与刘桢峤等(2019)的研究结果相符。同时, 研究区大部分野外实测巨砾位于风暴巨砾理论分布空间的中部偏上(图4), 这些巨砾是在波高为9m左右的风暴潮作用下搬运至现今位置的, 产生如此波高的台风强度可能与帕劳群岛的超强台风“宝霞”相当(Hongo et al, 2018)。故本文提出的基于理论空间分布的风暴巨砾成因判别方法具有一定的合理性。但若将此方法推广至更广的研究区, 则还需要作进一步完善。
首先, 本文试图从风暴水流的输运机制出发, 确定不同大小的巨砾在礁平台上的平面分布格局, 同样的方法也可运用于海啸的情形。将两种分布格局与实际情况相对照, 以判别成因类型。由于小东海环境条件的制约, 目前尚未就海啸成因的巨砾分布给出相应的模拟结果。今后可在进一步完善风暴水流输运理论的基础上, 结合海啸模拟方法, 构建出区分风暴巨砾与海啸巨砾的判别标准。
其次, 本文所刻画的风暴巨砾理论空间分布格局中, 巨砾起动位置假设为礁坪前缘。然而, 巨砾的初始位置可能处于水下礁缘斜坡, 也可能处于珊瑚礁坪上, 即其来源存在不确定性。若巨砾从水下礁缘斜坡开始搬运, 其整体搬运动力过程还应考虑斜坡段的动力过程及能量转换。同一风暴潮作用情况下, 从水下礁缘斜坡起动的巨砾的最终沉积位置与礁坪前缘的距离应小于从礁坪起动的巨砾的搬运距离, 故从水下礁缘斜坡起动的巨砾在风暴潮作用下形成的风暴巨砾理论空间分布应更靠近海岸线。本文研究因不能确定斜坡起动的条件而无法考虑这一因素, 今后如能就该问题进行补充和完善, 可进一步提高对巨砾成因的判定以及风暴潮能量计算的准确性。
最后, 本文的模拟只考虑了一种珊瑚的密度。不同种类的珊瑚所构成的巨砾密度范围可在0.8~ 2.7g·cm-3间不等, 差异较大(Spiske, 2008), 今后可进一步考虑不同种类和密度的珊瑚巨砾。此外, 不同大小、不同质量的巨砾与礁坪之间的摩擦系数可能也存在差异。尤其对于质量在1.0×104kg以上的珊瑚巨砾, 目前的方法只能提供初步成因判断。

5 结论

本文选择兼受风暴潮和海啸影响的小东海珊瑚礁海岸为研究区, 通过波浪能量分析推导出巨砾在风暴潮作用下的搬运动力过程, 进而对小东海风暴巨砾的理论搬运距离进行定量估算。将水动力学过程与风暴巨砾空间分布格局相结合, 提出巨砾成因区分的新方法。本文获取的主要结论有:
1) 小东海风暴巨砾的理论空间分布与风暴潮引起的海浪波高、巨砾质量和形状等因素密切相关。
2) 运用风暴巨砾的理论空间分布对小东海实测巨砾进行成因区分, 结果表明小东海巨砾由风暴潮搬运形成的可能性较大。尽管部分巨砾的沉积位置不在本文模拟得到的风暴巨砾理论分布范围内, 但对于其是否为海啸成因仍需进一步验证。
3) 基于风暴水流的输运机制可建立风暴巨砾和海啸巨砾在礁平台上的平面分布格局, 进而形成区分这两种巨砾的判别标准。为实现该目标, 今后需进一步研究的问题包括模拟假设条件的完善、巨砾密度和质量尺度效应的界定、风暴和海啸水流输运理论的提升。

References

Publishing order | Descend order by publishing year | Descend order by cited within
[1]
陈哲, 2017. 浅谈近年影响海南岛风暴潮的因素探讨[J]. 科技风, (16): 126-128, 130 (in Chinese).

[2]
高抒, 贾建军, 杨阳, 等, 2019. 陆架海岸台风沉积记录及信息提取[J]. 海洋学报, 41(10):141-160.

GAO SHU, JIA JIANJUN, YANG YANG, et al, 2019. Obtaining typhoon information from sedimentary records in coastal-shelf waters[J]. Haiyang Xuebao, 41(10):141-160 (in Chinese with English abstract).

[3]
黄德银, 施祺, 张叶春, 2005. 海南岛鹿回头与全新世高海平面[J]. 海洋地质与第四纪地质, 25(4):5-11.

HUANG DEYIN, SHI QI, ZHANG YECHUN, 2005. The coral reef and high sea level in Luhuitou, Hainan Island during Holocene[J]. Marine Geology & Quaternary Geology, 25(4):5-11 (in Chinese with English abstract).

[4]
李淑江, 李泽文, 范斌, 等, 2016. 海南岛东南近岸海浪观测及统计特征[J]. 海洋科学进展, 34(1):1-9.

LI SHUJIANG, LI ZEWEN, FAN BIN, et al, 2016. Wave observation and statistical analysis in the southeast coast of Hainan Island[J]. Advances in Marine Science, 34(1):1-9 (in Chinese with English abstract).

[5]
李文欢, 石海莹, 2013. 海南省风暴潮灾害预报及防范系统研究[M]. 北京: 海洋出版社(in Chinese).

[6]
练健生, 黄晖, 黄良民, 等, 2010. 三亚珊瑚礁及其生物多样性[M]. 北京: 海洋出版社: 106 (in Chinese).

[7]
刘桢峤, 周亮, 高抒, 2019. 基于地面3D激光扫描技术的海南岛南部海岸巨砾沉积研究[J]. 海洋学报, 41(11):127-141.

LIU ZHENQIAO, ZHOU LIANG, GAO SHU, 2019. Application of the terrestrial laser scanner to the coastal boulders on the southern coast of Hainan Island[J]. Hiayang Xuubao, 41(11):127-141 (in Chinese with English abstract).

[8]
毛龙江, 张永战, 魏灵, 等, 2006. 海南岛三亚湾海滩研究[J]. 第四纪研究, 26(3):477-484.

MAO LONGJIANG, ZHANG YONGZHAN, WEI LING, et al, 2006. Study on beach characteristics in Sanya area of Hainan Island[J]. Quaternary Sciences, 26(3):477-484 (in Chinese with English abstract).

[9]
苗庆生, 2011. 南海台风浪场特征的计算与分析[D]. 青岛: 中国海洋大学.

MIAO QINGSHENG, 2011. Calculation and analysis of typhoon waves in South China Sea[D]. Qingdao: Ocean University of China (in Chinese with English abstract).

[10]
王月, 沈建伟, 龙江平, 2011. 海南岛三亚小东海珊瑚礁坪生态沉积分带和碳酸盐沉积作用[J]. 中国科学 D辑: 地球科学, 41(3):362-374.

WANG YUE, SHEN JIANWEI, LONG JIANGPING, 2011. Ecological-sedimentary zonations and carbonate deposition, Xiaodonghai Reef Flat, Sanya, Hainan Island, China[J]. Science China Earth Sciences, 54(3):359-371.

[11]
叶琳, 于福江, 吴玮, 2005. 我国海啸灾害及预警现状与建议[J]. 海洋预报, 22(S1):147-157.

YE LIN, YU FUJIANG, WU WEI, 2005. The disaster and warning of tsunami in China and the suggestion in future[J]. Marine Forecasts, 22(S1):147-157 (in Chinese with English abstract).

[12]
尹洪强, 梁书秀, 2014. 基于SWAN模式对台风海燕风浪的推算[J]. 中国水运, 14(6):105-107, 111 (in Chinese).

[13]
BAGNOLD R A, 1963. Beach and nearshore processes. Part I, mechanics of marine sedimentation[M]// THORNE C R, MACARTHUR R C, BRADLEY J B. The physics of sediment transport by wind and water. New York: American Society of Civil Engineers: 507-528.

[14]
BARBANO M S, PIRROTTA C, GERARDI F, 2010. Large boulders along the South-eastern Ionian coast of Sicily: storm or tsunami deposits?[J]. Marine Geology, 275(1-4):140-154.

DOI

[15]
BIOLCHI S, FURLANI S, ANTONIOLI F, et al, 2016. Boulder accumulations related to extreme wave events on the eastern coast of Malta[J]. Natural Hazards and Earth System Sciences, 16(3):737-756.

[16]
BRYANT E A, NOTT J, 2001. Geological indicators of large tsunami in Australia[J]. Natural Hazards, 24(3):231-249.

DOI

[17]
CAUSON D J, GAUCI R, 2017. Evidence of extreme wave events from boulder deposits on the South-East Coast of Malta (Central Mediterranean)[J]. Natural Hazards, 86(S2):543-568.

DOI

[18]
COREY A T, 1949. Influence of shape on the fall velocity of sand grains[D]. Fort Collins: Colorado Agricultural and Mechanical College.

[19]
COSTA P J M, ANDRADE C, FREITAS M C, et al, 2011. Boulder deposition during major tsunami events[J]. Earth Surface Processes and Landforms, 36(15):2054-2068.

DOI

[20]
ETIENNE S, PARIS R, 2010. Boulder accumulations related to storms on the south coast of the Reykjanes Peninsula (Iceland)[J]. Geomorphology, 114(1-3):55-70.

DOI

[21]
ETIENNE S, BUCKLEY M, PARIS R, et al, 2011. The use of boulders for characterising past tsunamis: lessons from the 2004 Indian Ocean and 2009 South Pacific tsunamis[J]. Earth-Science Reviews, 107(1-2):76-90.

DOI

[22]
FELTON E A, CROOK K A W, 2003. Evaluating the impacts of huge waves on rocky shorelines: An essay review of the book ‘Tsunami - the underrated hazard’[J]. Marine Geology, 197(1-4):1-12.

DOI

[23]
GOFF J, DOMINEY-HOWES D, 2009. Australasian palaeotsunamis — Do Australia and New Zealand have a shared trans-Tasman prehistory?[J]. Earth Science Reviews, 97(1-4):147-154.

DOI

[24]
GOFF J, WEISS R, COURTNEY C, et al, 2010. Testing the hypojournal for tsunami boulder deposition from suspension[J]. Marine Geology, 277(1-4):73-77.

DOI

[25]
GOTO K, CHAVANICH S A, IMAMURA F, et al, 2007. Distribution, origin and transport process of boulders deposited by the 2004 Indian ocean tsunami at Pakarang cape, Thailand[J]. Sedimentary Geology, 202(4):821-837.

DOI

[26]
GOTO K, MIYAGI K, KAWAMATA H, et al, 2010a. Discrimination of boulders deposited by tsunamis and storm waves at Ishigaki Island, Japan[J]. Marine Geology, 269(1-2):34-45.

DOI

[27]
GOTO K, KAWANA T, IMAMURA F, 2010b. Historical and geological evidence of boulders deposited by tsunamis, southern Ryukyu Islands, Japan[J]. Earth-Science Reviews, 102(1-2):77-99.

DOI

[28]
GOTO K, MIYAGI K, KAWANA T, et al, 2011. Emplacement and movement of boulders by known storm waves-field evidence from the Okinawa Islands, Japan[J]. Marine Geology, 283(1-4):66-78.

DOI

[29]
GOTO K, SUGAWARA D, IKEMA S, et al, 2012. Sedimentary processes associated with sand and boulder deposits formed by the 2011 Tohoku-Oki tsunami at Sabusawa Island, Japan[J]. Sedimentary Geology, 282:188-198.

DOI

[30]
HAYNE M, CHAPPELL J, 2001. Cyclone frequency during the last 5000 years at Curacoa Island, north Queensland, Australia[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 168(3-4):207-219.

DOI

[31]
HONGO C, KURIHARA H, GOLBUU Y, 2018. Coral boulders on Melekeok reef in the Palau Islands: an indicator of wave activity associated with tropical cyclones[J]. Marine Geology, 399:14-22.

DOI

[32]
LORANG M S, 2000. Predicting threshold entrainment mass for a boulder beach[J]. Journal of Coastal Research, 16(2):432-445.

[33]
LORANG M S, 2011. A wave-competence approach to distinguish between boulder and megaclast deposits due to storm waves versus tsunamis[J]. Marine Geology, 283(1-4):90-97.

DOI

[34]
MAK S, CHAN L S, 2007. Historical tsunamis in South China[J]. Natural Hazards, 43(1):147-164.

DOI

[35]
MASTRONUZZI G, SANSÒ P, 2004. Large boulder accumulations by extreme waves along the Adriatic coast of southern Apulia (Italy)[J]. Quaternary International, 120(1):173-184.

DOI

[36]
MAY S M, ENGEL M, BRILL D, et al, 2015. Block and boulder transport in Eastern Samar (Philippines) during Supertyphoon Haiyan[J]. Earth Surface Dynamics, 3(4):543-558.

DOI

[37]
NANDASENA N A K, PARIS R, TANAKA N, 2011. Reassessment of hydrodynamic equations: Minimum flow velocity to initiate boulder transport by high energy events (storms, tsunamis)[J]. Marine Geology, 281(1-4):70-84.

DOI

[38]
NANDASENA N A K, TANAKA N, SASAKI Y, et al, 2013. Boulder transport by the 2011 Great East Japan tsunami: comprehensive field observations and whither model predictions?[J]. Marine Geology, 346:292-309.

DOI

[39]
NIELSEN P, 1992. Coastal bottom boundary layers and sediment transport[M]. Singapore: World Scientific Publishing.

[40]
NOTT J, 1997. Extremely high-energy wave deposits inside the Great Barrier Reef, Australia: determining the cause—tsunami or tropical cyclone[J]. Marine Geology, 141(1-4):193-207.

DOI

[41]
NOTT J, 2003. Waves, coastal boulder deposits and the importance of the pre-transport setting[J]. Earth and Planetary Science Letters, 210(1-2):269-276.

DOI

[42]
PARIS R, LAVIGNE F, WASSMER P, et al, 2007. Coastal sedimentation associated with the December 26, 2004 tsunami in Lhok Nga, west Banda Aceh (Sumatra, Indonesia)[J]. Marine Geology, 238(1-4):93-106.

DOI

[43]
POWER H E, BALDOCK T E, CALLAGHAN D P, et al, 2013. Surf zone states and energy dissipation regimes - A similarity model[J]. Coastal Engineering Journal, 55(1): 1350003-1-1350003-18.

[44]
SAINTILAN N, ROGERS K, 2005. Recent storm boulder deposits on the Beecroft Peninsula, New South Wales, Australia[J]. Geographical Research, 43(4):429-432.

DOI

[45]
SCHEFFERS A, KELLETAT D, 2003. Sedimentologic and geomorphologic tsunami imprints worldwide—a review[J]. Earth-Science Reviews, 63(1-2):83-92.

DOI

[46]
SCHEFFERS A, KELLETAT D, VÖTT A, et al, 2008. Late Holocene tsunami traces on the western and southern coastlines of the Peloponnesus (Greece)[J]. Earth and Planetary Science Letters, 269(1-2):271-279.

DOI

[47]
SCICCHITANO G, MONACO C, TORTORICI L, 2007. Large boulder deposits by tsunami waves along the Ionian coast of South-eastern Sicily (Italy)[J]. Marine Geology, 238(1-4):75-91.

DOI

[48]
SPISKE M, BÖRÖCZ Z, BAHLBURG H, 2008. The role of porosity in discriminating between tsunami and hurricane emplacement of boulders — A case study from the lesser Antilles, Southern Caribbean. Earth and Planetary Science Letters, 268(3):384-396.

DOI

[49]
SUN LIGUANG, ZHOU XIN, HUANG WEN, et al, 2013. Preliminary evidence for a 1000-year-old tsunami in the South China Sea[J]. Scientific Reports, 3(1):1655.

DOI

[50]
TERRY J P, LAU A Y A, ETIENNE S, 2013. Reef-platform coral boulders: evidence for high-energy marine inundation events on tropical coastlines[M]. Singapore: Springer: 63-65.

[51]
TERRY J P, ETIENNE S, 2014. Potential for timing high-energy marine inundation events in the recent geological past through age-dating of reef boulders in Fiji[J]. Geoscience Letters, 1(1):14.

DOI

[52]
TERRY J P, OLIVER G J H, FRIESS D A, 2016. Ancient high-energy storm boulder deposits on Ko Samui, Thailand, and their significance for identifying coastal hazard risk[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 454:282-293.

DOI

[53]
WEISS R, 2012. The mystery of boulders moved by tsunamis and storms[J]. Marine Geology, 295-298:28-33.

DOI

[54]
YU KEFU, ZHAO JIANXIN, SHI QI, et al, 2009. Reconstruction of storm/tsunami records over the last 4000 years using transported coral blocks and lagoon sediments in the southern South China Sea[J]. Quaternary International, 195(1-2):128-137.

DOI

Options
Outlines

/

〈 〉