海洋学报

剖面编号	礁坪宽度/m	礁缘水深/m	礁后海滩坡度	礁后海滩宽度/m
P1	1 146.00	1.27	0.13	30.31
P2	288.59	0.69	0.05	66.35
P3	731.65	2.21	0.12	40.75

剖面编号	礁坪宽度/m	礁缘水深/m	礁后海滩坡度	礁后海滩宽度/m
P1	1 146.00	1.27	0.13	30.31
P2	288.59	0.69	0.05	66.35
P3	731.65	2.21	0.12	40.75

剖面编号	E_ss0/(10⁴m²)	E_ig0/(10⁴m²)	[(E_ss0−E_ssλ)/E_ssλ]/%	[(E_ig0−E_igλ)/E_igλ]/%	(E_ss0/E₀)/%	(E_ig0/E₀)/%
P1	0.41	0.18	–99.59	105.21	69.58	30.42
P2	0.14	0.95	–99.85	3 729.29	13.11	86.89
P3	5.11	3.99	–95.19	9 776.25	56.14	43.86

剖面编号	E_ss0/(10⁴m²)	E_ig0/(10⁴m²)	[(E_ss0−E_ssλ)/E_ssλ]/%	[(E_ig0−E_igλ)/E_igλ]/%	(E_ss0/E₀)/%	(E_ig0/E₀)/%
P1	0.41	0.18	–99.59	105.21	69.58	30.42
P2	0.14	0.95	–99.85	3 729.29	13.11	86.89
P3	5.11	3.99	–95.19	9 776.25	56.14	43.86

剖面编号	$ {A}_{rp}/{m}^{1/3} $	$\varepsilon \left({A}_{rp}\right)/\text{%}$
P1	0.26	3.71
P2	0.17	9.56
P3	0.28	1.06

剖面编号	$ {A}_{rp}/{m}^{1/3} $	$\varepsilon \left({A}_{rp}\right)/\text{%}$
P1	0.26	3.71
P2	0.17	9.56
P3	0.28	1.06

不同受损情况下珊瑚礁海岸动力地貌差异性分析

PDF下载

陈燕珊 ¹^,² , 戚洪帅 ²^,⁴^,⁵^,^* , 杨清书 ¹^,³^,⁴ , 蔡锋 ²^,⁴ , 刘根 ²^,⁴^,⁵ , 朱君 ²^,⁴ , 赵绍华 ²

海洋学报 | 论文 2022,44(3): 61-69

收起

海洋学报 | 论文 2022, 44(3): 61-69

不同受损情况下珊瑚礁海岸动力地貌差异性分析

全屏

陈燕珊¹^,², 戚洪帅²^,⁴^,⁵^,^*, 杨清书¹^,³^,⁴, 蔡锋²^,⁴, 刘根²^,⁴^,⁵, 朱君²^,⁴, 赵绍华²

作者信息

1.中山大学海洋工程与技术学院中山大学河口海岸研究所，广东广州 510275

2.自然资源部第三海洋研究所海洋与海岸地质研究室，福建厦门 361005

3.河口水利技术国家地方联合工程实验室，广东广州 510275

4.广东省海岸与岛礁工程技术研究中心，广东珠海 510275

5.自然资源部北部湾滨海湿地生态系统野外科学观测研究站，广西北海 536015

陈燕珊（1995-），女，广东省深圳市人，主要从事海岸动力地貌研究。E-mail：chenysh76@mail2.sysu.edu.cn

通讯作者:

戚洪帅，男，研究员，主要从事海岸与海滩过程以及海滩养护技术研究。E-mail: qihongshuai@tio.org.cn

Morphodynamic analysis on fringing reef coasts under different damage conditions

Yanshan Chen¹^,², Hongshuai Qi²^,⁴^,⁵^,^*, Qingshu Yang¹^,³^,⁴, Feng Cai²^,⁴, Gen Liu²^,⁴^,⁵, Jun Zhu²^,⁴, Shaohua Zhao²

Affiliations

1. Institute of Estuarine and Coastal Research, School of Marine Engineering and Technology, Sun Yat-sen University, Guangzhou 510275, China

2. Laboratory of Ocean and Coast Geology, Third Institute of Oceanography, Ministry of Natural Resources, Xiamen 361005, China

3. State and Local Joint Engineering Laboratory of Estuarine Hydraulic Technology, Guangzhou 510275, China

4. Guangdong Provincial Engineering Research Center of Coasts, Islands and Reefs, Zhuhai 510275, China

5. Observation and Research Station of Coastal Wetland Ecosystem in Beibu Gulf, Minstry of Natural Resources, Beihai 536015, China

出版时间: 2022-03-01 doi: 10.12284/hyxb2022025

文章导航

摘要

收起

探索珊瑚礁与海滩地貌之间动力地貌联系是认识珊瑚礁海岸变化的重要一环。本文以雷州半岛徐闻西落港珊瑚礁海岸为研究对象，应用RTK-GPS和无人船开展岸滩剖面和近岸水下地形的测量、结合海滩沉积物分析，基于FUNWAVE-TVD数值模型模拟并分析不同珊瑚礁地形地貌条件下波浪动力传播过程。结果显示，研究区珊瑚礁水下地形是影响礁后海滩地貌的主要因素。礁体形态不同，导致其礁后海滩在珊瑚礁地形控制下，短波波能和次重力波波能沿程呈现不同变化规律，最终导致礁后海滩的近岸波能主控频段的差异，在较窄的珊瑚礁海岸，次重力波占比较大。在不同主控频段波浪驱动下，礁后海滩平衡剖面呈现差异性特征，Muñóz-Pérez提出的珊瑚礁后海滩平衡剖面拟合中未考虑该因素的影响，需要在此认识基础上进一步改进。

关键词

珊瑚礁海岸 / FUNWAVE-TVD / 次重力波 / 动力地貌 / 海滩平衡剖面

Abstract

收起

Morphodynamic relationships between fringing reefs and back-reef beach play key roles in understanding the evolution of coral coasts. The fringing reef coasts of Leizhou Peninsula, in Xiluo Port, Xuwen, are taken as the research object in this paper. RTK-GPS and unmanned-surface-vessel are used to measure and analyze the topography of the beach profile and reef as well as the sediment sampling in beach. FUNWAVE-TVD is used to simulate and analyze the hydrodynamic process of short wave and infragravity wave cross the different reef. The results show that the reef topography in the study area has significant effects on the morphology of back-reef beach. Whether the nearshore wave energy dominated by the short wave band or the infragravity band is controlled by the topography. In the narrow reef, there is a predominance of infragravity energy near the back-reef beach. The dominated wave band in the nearshore wave energy is the main factor shape the back-reef beach equilibrium profile. However, the beach equilibrium profile model of coral coast proposed by Muñóz-Pérez did not take this factor into account, it need to further research.

Key words

fringing reef coasts / FUNWAVE-TVD / infragravity waves / morphodynamic / beach equilibrium profile

引用本文

陈燕珊, 戚洪帅, 杨清书, 蔡锋, 刘根, 朱君, 赵绍华. 不同受损情况下珊瑚礁海岸动力地貌差异性分析. 海洋学报, 2022 , 44 (3) : 61 -69 . DOI: 10.12284/hyxb2022025

Yanshan Chen, Hongshuai Qi, Qingshu Yang, Feng Cai, Gen Liu, Jun Zhu, Shaohua Zhao. Morphodynamic analysis on fringing reef coasts under different damage conditions[J]. Haiyang Xuebao, 2022 , 44 (3) : 61 -69 . DOI: 10.12284/hyxb2022025

正文

收起

1 引言

收起

在全球气候变化和人类活动频繁背景下，珊瑚礁海岸最常见的地貌组合——珊瑚礁–海滩系统对环境变化响应极其敏感，探索珊瑚礁与海滩地貌之间的动力地貌联系是认识珊瑚礁海岸变化的重要一环，可为海岸防护和生态修复提供科学依据。

波浪过程是珊瑚礁海岸研究中的重要科学问题，礁体能有效削弱外海传至近岸波能，对岸滩起到一定保护作用^[1]，目前多数研究集中于对影响礁体波浪传播过程的因素进行分析。研究表明，礁坪水深、礁坪宽度、岸滩坡度、礁体形态通过影响波浪传播变化如波浪破碎、摩阻损失、礁坪共振等过程，引起不同频段波能的重新分配，导致近岸波能主控频段分布呈现差异性特征^[2-8]。不同地貌对动力、沉积输运过程差异导致海滩地貌特征不同，尽管前人有所研究^[7,9-14]，但针对珊瑚礁海岸区域，近岸不同主控频段波浪影响下的海滩地貌特征分析研究存在不足，因此，有必要开展珊瑚礁海岸对区域波浪能量变化以及礁后海滩地貌的研究，探讨礁体地形与礁后海滩动力地貌之间的联系。数值模型是研究珊瑚礁海岸波浪过程常用的方法之一，相位平均模型如SWAN、Delft-3D、MIKE21 SW等模型对波浪过程的求解一般基于能量守恒原理，但此类模型无法精细刻画珊瑚礁海岸波浪传播过程。相位解析模型如非静压模型、Boussinesq类模型可模拟高色散完全非线性波浪过程，计算精度和效率相对更高。基于激波捕捉的完全非线性Boussinesq方程模型FUNWAVE可避免因引入人工率定经验参数带来的不确定性，被成功应用于珊瑚礁海岸波浪传播的模拟^[4,15]。本文通过对华南雷州半岛徐闻西落港珊瑚礁海岸沉积以及地形特征分析，运用FUNWAVE-TVD模拟分析不同礁体地形上波能传播过程，并结合礁后海滩平衡剖面特征，探讨珊瑚礁海岸礁后海滩动力地貌差异性及其机理。

2 材料和方法

收起

2.1 野外调查

研究区为徐闻西落港珊瑚礁海岸，地处广东雷州半岛西南海岸、灯楼角北部区域，与海南岛北部隔海相望（图1a），是我国珊瑚礁海岸（岸礁）的主要分布区域之一，该区域珊瑚礁因受到如炸鱼、礁坪采掘、养殖引起水质污染等人为活动的影响而发生一定损毁以及退化^[16]。该海区为热带季风气候，属全日潮，潮差约为1.5 m，波浪主要由外海传入，具有明显季节性变化，波高约为0.6 m^[17]。

为了研究其地形地貌特征，对该区域开展了野外调查。调查断面根据历史资料选择礁坪存在明显差异的剖面进行布设，剖面调查时间是在2020年8月10–11日（夏季），该处主要受夏季西南向浪影响，调查时间可以适当反映海滩状态。研究区域从南向北布设P1、P2和P3 3个珊瑚礁海岸剖面（图1b），3个剖面礁体受损程度不同，P2剖面海滩上珊瑚礁碎屑数量较多且颗粒较大，其次是P3剖面，而P1剖面海滩上珊瑚礁碎屑相对更少（图1c至图1e）。岸滩地形测量采用STONEX S9Ⅱ型的RTK-GPS，水平和高程误差分别为±1 cm和±2 cm，从后方植被线向海测至涉水深度（图1c至图1e虚线）。水下地形测量采用云洲SE40型的无人船搭载单波束完成，测量范围水深0.5～10 m（图1b虚线）水平误差和水深精度分别为±10 cm和±5 cm。研究区海滩沉积物差异不大，选择P1剖面进行表层沉积物采样，采样高程间距为0.5 m，共获得10个表层沉积物样品。

2.2 沉积物粒度分析

沉积物粒度分析采用激光粒度分析法和筛分法相结合的方式，对于相对较粗的沉积物样品（含砾石），取适量烘干，加入去离子水充分分散，使用筛分法对粒度大于2 mm的颗粒进行粒度分级。而对于其中粒度小于2 mm的部分沉积物，取适量样品溶液，加入0.05 mol/L六偏硫酸钠溶液约5 mL，使用马尔文公司生产的MS2000激光粒度仪（0.02～2 000 μm）进行粒度分析。粒度参数计算采用福克–沃德粒度参数公式^[18]，粒级标准采用尤登–温德华等比制Φ值粒级标准^[19]，沉积物类型采用谢帕德的沉积物粒度三角图解法^[20]。

2.3 数值模拟与验证

数值模型采用基于完全非线性Boussinesq方程的FUNWAVE-TVD 3.5模型，该模型适用于模拟近岸波浪浅化、反射、折射、破碎以及波浪在岸滩区域的爬高与下冲等过程^[21]。模拟区域徐闻剖面总长为2 000 m，礁体底床摩擦系数C_d=0.05，入射波谱型采用研究中常用的JONSWAP不规则波经验谱，谱峰升高因子γ=3.3，在礁缘外入射。徐闻月平均波高为0.6 m，月平均波周期为3.5～4.9 s^[16-17,22]，故正常波况波浪设置为谱峰周期T_p=4.6 s，入射波高H_i=0.6 m。采样频率为1.0 s，模拟时长以水位时间序列达到稳定状态并持续一段时间为止。模型区域设置如图2。模型输出水位时间序列最后相对稳定的200 s的数据用于谱分析，谱分析参数设置取默认值，采用Hanning窗口，50%重叠率，频域分辨率为0.009 8 Hz。短波频率一般介于0.04～1.00 Hz之间，次重力波指频率低于短波的海表面流体波动，频率一般介于0.005～0.04 Hz之间。由水位时间序列的谱分析可得波能密度–频率（S（f）-f）分布关系，进而可计算全谱波波能（式（1））、短波波能（式（2））以及次重力波波能（式（3））（f为频率，单位为Hz）。

(1)

$ E = \int_{{\text{ - }}\infty }^\infty {S\left( f \right)} {\text{d}}f, $

(2)

$ {E_{{\text{ss}}}} = \int_{0.04}^\infty {S\left( f \right)} {\text{d}}f, $

(3)

$ {E_{{\text{ig}}}} = \int_{0.005}^{0.04} {S\left( f \right)} {\text{d}}f . $

为对FUNWAVE-TVD数值模型模拟珊瑚礁海岸波浪传播过程的适用性进行验证，本文将模型模拟结果与密歇根大学进行的基于概化珊瑚礁地形的水槽模型实验^[23]结果进行误差分析（式（4）^[24]，图3），结果表明数值模型结果与水槽模型结果吻合较好，有效波高模拟误差Skill≈0.93，次重力波波高模拟误差Skill≈0.79，故该数值模型模拟精确度较高，适用于珊瑚礁海岸波浪过程的模拟。在对模型进行验证后，接下来对模型数值稳定性进行探究：数值稳定性与模型使用网格大小有关，一般以波长作为网格大小选取的参考依据。本文数值模型以2倍关系加密计算网格，计算区域网格大小ρΔx分别为L_p/21、L_p/41、L_p/82、L_p/165（L_p为谱峰波长），选取离岸50 m、100 m处波面变化稳定性进行了分析（图4），可知模型在Δx=L_p/165，即单位网格大小约为0.2 m时，结果已收敛。

(4)

$ S kill = 1 - \sqrt {\frac{{{{(1} /{N)\displaystyle\sum\limits_{n = 1}^N {{{({X_{\rm{p}}} - {X_{\rm{o}}})}^2}} }}}}{{{{(1} /{N)\displaystyle\sum\limits_{n = 1}^N {X_{\rm{o}}^2} }}}}}, $

式中，N是实验点总数；X_p是模型模拟值；X_o是实测值；Skill值越接近1，表明误差越小。

3 结果

收起

3.1 珊瑚礁海岸地形地貌与沉积特征

徐闻西落港珊瑚礁海岸P1、P2、P3 3个剖面礁坪基本淹没于平均海平面之下，且主要分布于潮下带，3个剖面珊瑚礁和礁后海滩地貌呈现明显差异性特征。P1剖面珊瑚礁发育较完善，礁坪宽度是3条剖面中最宽广的，礁缘水深为1.27 m，礁后海滩宽度小而坡度大，后滨无滩肩发育，整体处于稳定状态；P2剖面珊瑚礁受破坏较严重，礁坪窄短，礁缘水深为0.69 m，礁后海滩较宽缓，处于侵蚀状态；P3剖面活体珊瑚覆盖率高，礁坪宽度较P1小，但礁缘水深较大，礁后海滩滩肩发育，海滩坡度与P1近似，稳定性较好（图5，表1）。比较来看，P1和P3剖面海滩的反射性较强，而P2剖面滩面耗散性较强^[24]。

海滩表层沉积物类型以砂为主，中值粒径在0.24 Φ～1.95 Φ范围内变动，由岸向海，表层沉积物粒径逐渐变细，平均高潮面和平均海平面间的区域岸滩表层沉积物粒径变动较大，平均海平面以下，表层沉积物粒径基本不变（图6）。

3.2 不同剖面的波能差异

基于实测剖面地形，对3个不同特征的珊瑚礁海岸进行数值试验。模拟结果表明，受岸滩差异性反射作用影响，各剖面礁后海滩附近波浪呈不同程度振荡趋势分布，研究区不同剖面礁后海滩（0 m等深线）短波波能E_ss0、次重力波波能E_ig0、短波波能占全谱波波能比重（E_ss0/E₀）、次重力波波能占全谱波波能比重（E_ig0/E₀）、短波波能相对入射短波波能变化率（(E_ss0−E_ssλ)/E_ssλ）、次重力波波能相对入射次重力波波能变化率（(E_ig0−E_igλ)/E_igλ）均呈现明显差异性特征（图7，表2）。

各剖面礁后海滩短波相对外海都是衰减的，P2剖面近岸短波波能相对外海短波波能衰减率最大，短波波能最小；P1剖面近岸短波波能衰减率略小于P2剖面；P3剖面短波波能衰减率最小，到达礁后海滩短波剩余波能最大。各剖面礁后海滩次重力波波能相对外海都显著增大，其中P3剖面次重力波波能增长率最大，P1剖面次重力波波能增长率最小。由于各剖面短波波能和次重力波波能变化率不同，礁后海滩各频段波能占全谱波波能比重出现差异。P1和P3剖面礁后海滩短波波能占全谱波波能比重都大于对应剖面次重力波波能占全谱波波能比重，其中次重力波波能占全谱波波能比例分别约为30.42%和43.86%；P2剖面则与之不同，礁后海滩次重力波波能占全谱波波能的86.89%。总体来看，耗散性强的P2剖面礁后海滩波能以次重力波频段波浪为主，而对于反射性较强的P1和P3剖面礁后海滩波能主要以短波波浪为主。

4 讨论

收起

4.1 剖面波能差异性原因分析

研究区各剖面波能沿程变化规律基本与前人研究结果一致^[25-27]，模型试验可以较好的体现其波浪过程差异。不同剖面波能分布呈现明显差异，因波浪传播过程受礁体水深、礁坪宽度、地形坡度等因素影响显著。

波浪由外海传播至近岸时，波浪在礁缘附近发生破碎，波浪破碎引起水体强烈的紊动和漩涡，导致入射短波波能显著衰减。在波浪继续向岸传播过程中，礁坪水深越小，因礁坪底床对波浪水流的摩阻力而引起的波能损失更大，而礁坪宽度越大意味着波浪传播需经过更长的路程，导致波能衰减量越多^[5-6]，对于岸滩坡度相近的P1和P3剖面来说，P1剖面由于礁体水深较小且礁坪宽广，短波波能损失较大，导致到达近岸短波波能较小，相反，P3剖面由于礁体水深较大且礁坪宽度较小，其短波波能损失较小，到达近岸短波波能更大。礁冠处由于水深突变，进而引起破波点处更大的湍动能，波浪破碎效应更强烈，相对无礁冠结构礁体，存在礁冠结构的礁体对入射波能削弱性更大^[7-8]。故导致礁坪宽度最窄的P2剖面近岸短波剩余波能是3个剖面中最小，就因其礁缘处存在礁冠结构。

波浪破碎同样也会导致次重力波波能衰减，但与此同时，基于移动破波点机制的驱动，波能由短波频段向次重力波频段转移，对波浪的落后与超前时间相关分析或谱分析都发现了礁缘附近次重力波生成的现象^[27-30]。礁坪宽度、礁体水深和岸滩坡度对次重力波波能损耗影响与对短波类似^[3-5]。P1剖面破波带外岸坡反射性强，但其礁坪浅且宽广，到达近岸时短波和次重力波波能都衰减显著，而P3剖面礁体水深较大，进而导致其近岸短波和次重力波波能较大。P2剖面礁坪宽度小，次重力波波能衰减量相对较少。由于各剖面地形的差异影响，对应剖面短波波能和次重力波波能沿程呈现不同变化规律，最终导致礁后海滩的近岸波能主控频段的差异。P1和P3剖面近岸短波波能均大于次重力波波能，导致这两个剖面近岸波能以短波频段为主，相反，P2剖面近岸短波波能小于次重力波波能，导致该剖面近岸波能以次重力波频段为主。

4.2 礁后海滩动力地貌特征

徐闻西落港珊瑚礁海岸3条典型剖面珊瑚礁地形和礁后海滩地貌表现出明显的差异性，从动力地貌学角度两者是相互联系的。海滩地貌形态受近岸波浪、沉积物以及地质等多种因素影响，是区域波能在滩面耗散并驱动沉积物分配的结果，趋于形成海滩平衡剖面^[31]。海滩平衡剖面形态模式有助于理解岸滩平衡剖面形态与近岸波能之间的关系。前人对海滩平衡剖面的研究众多^[32-38]，其中最具代表性的是Dean基于美国大西洋海岸和墨西哥海岸504条海滩剖面的拟合结果，其提出的海滩平衡剖面模式，被广泛应用于浪控海滩。然而对于珊瑚礁海岸而言，礁体地形对礁后海滩水动力环境影响显著，该平衡剖面模式存在较大局限性，Muñóz-Pérez等^[9]基于Dean的海滩平衡剖面模型，提出了适用于珊瑚礁海岸礁后海滩平衡剖面模型

(5)

$ h = {A_{rp}}{x^m}, $

上式拟合起点位于坐标系原点，根据本文实测海滩地形，对海滩平衡剖面拟合起点进行相应的平移转化

(6)

$ h\text{′}={A}_{rp}(x-{x}_{0}{)}^{m}-{h}_{0}, $

式中，h为拟合起点为坐标系原点的剖面高程（单位：m）；h′为拟合起点平移转化后的剖面高程（单位：m）；h₀为岸线所在位置高程（平均大潮高潮线，单位：m）；x为离岸距离（单位：m）；x₀为岸线位置（单位：m）；A_rp、m为拟合常数；Muñóz-Pérez等^[9]取m=2/3。虽然该剖面模型形式与Dean海滩平衡剖面模型一致，但其A_rp值不同于Dean模型中的A值。A_rp与A满足

(7)

$ {A_{rp}}/A = \left[ {\beta /\varGamma } \right]_{{h_r}}^{4/3} \text{，} $

式中，Γ和β分别为受珊瑚礁保护和不受珊瑚礁保护剖面破波指数；h_r为礁坪水深（单位：m）。

对于一个宽广的礁体（宽度→∞），Γ值位于0.35～0.55之间^[39]，而β值位于0.65～1.1之间^[40]，故A_rp=1.250～4.604A。基于Muñóz-Pérez等^[9]礁后海滩平衡剖面形式对研究区礁后海滩进行拟合，计算误差ε（式（6））并取误差最小的拟合值A_rp（图8，表3）。拟合结果表明P1和P3剖面拟合误差较小，而P2剖面拟合误差达36.6%，差异性较大。

(8)

$ \varepsilon = \frac{\displaystyle\sum \left( h'_i - h'_{pi} \right)^2}{\displaystyle\sum {h'}^2_i} \times 100{\text{%}}, $

式中，$h_i^{\text{′}} $和$h_{pi}^{\text{′}} $分别为剖面实测高程和拟合高程值（单位：m）。

从海滩平衡剖面模型（式（6））数学形式上可知，拟合参数m的取值决定海滩剖面曲线的反射性强度，由前文分析可知P2剖面反射性弱，而P1和P3剖面反射性较强，故P2剖面与P1和P3剖面应采用不同的m值进行拟合。究其原因，Muñóz-Pérez等^[9]提出的珊瑚礁海滩平衡剖面模式是基于对西班牙南部海岸的珊瑚礁地形的研究而提出的，该区域礁坪水深都较大且礁体形态较平整，模型虽然考虑了珊瑚礁地形对波浪的作用结果，但未能体现不同礁体地形下的波浪传播差异以及该差异影响下波能主控频段的分布变化。不同礁体地形影响下的近岸波能主控频段是不同的，P2剖面礁后海滩处波能主控频段为次重力波频段，而P1和P3剖面礁后海滩处波能主控频段为短波频段。因此，结合实测地形地貌结果、平衡海滩剖面拟合结果和礁后海滩波浪分布结果可知，不同频段波浪影响下，礁后海滩地貌特征呈现差异性特征，次重力波主控海岸滩面耗散性较强，呈侵蚀地貌类型；短波主控海岸滩面反射性较强，海滩相对稳定。

5 结论

收起

基于对徐闻西落港珊瑚礁海岸典型断面的调查和基于FUNWAVE-TVD模型的波浪传播过程进行模拟分析，分析了不同珊瑚礁地形地貌条件下波浪动力传播过程，进而揭示了其对礁后海滩的动力地貌特征影响机理。

（1）在人为损毁下，珊瑚礁海岸地形呈现明显的差异性，导致其礁后海滩也呈现出不同的地貌特征，在礁坪发育好的的岸段，礁后海滩较窄且陡，滩面反射性较强；而礁坪受损的岸段，礁后海滩更宽缓，滩面耗散性较强。

（2）FUNWAVE-TVD模型分析表明，礁坪和礁冠地形是珊瑚礁波浪能量变化的主要因素，不同珊瑚礁水下地形对近岸水动力影响明显，由于短波波能和次重力波波能沿程呈现不同变化规律，最终导致礁后海滩的近岸波能主控频段的差异，在较窄的珊瑚礁海岸，次重力波占比较大。

（3）近岸水动力差异导致礁后海滩不同的动力地貌特征，塑造不同的平衡剖面形态。不同礁体地形影响下的近岸波能主控频段是不同的，进而导致礁后海滩平衡剖面呈现差异性特征，现有平衡剖面拟合中皆未考虑该因素的影响。

致谢　感谢自然资源部第三海洋研究所海洋与海岸地质研究室海滩组的何岩雨、肖哲宇、宋嘉诚以及胥海涛为野外调查、数据处理以及论文修改提供的帮助。

基金

收起

自然科学基金重点项目(41930538)
自然科学基金面上项目(42076058)
自然资源部第三海洋研究所基科项目(2019026)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

Young I R. Wave transformation over coral reefs[J]. Journal of Geophysical Research: Oceans, 1989, 94(C7): 9779−9789.

Beetham E, Kench P S, O’Callaghan J, et al. Wave transformation and shoreline water level on Funafuti Atoll, Tuvalu[J]. Journal of Geophysical Research: Oceans, 2016, 121(1): 311−326.

Nwogu O, Demirbilek Z. Infragravity wave motions and runup over shallow fringing reefs[J]. Journal of Waterway, Port, Coastal, and Ocean Engineering, 2010, 136(6): 295−305.

Su S F, Ma Gangfeng, Hsu T W. Boussinesq modeling of spatial variability of infragravity waves on fringing reefs[J]. Ocean Engineering, 2015, 101: 78−92.

Quataert E, Storlazzi C, Rooijen A V, et al. The influence of coral reefs and climate change on wave-driven flooding of tropical coastlines[J]. Geophysical Research Letters, 2015, 42(15): 6407−6415.

Shimozono T, Tajima Y, Kennedy A B, et al. Combined infragravity wave and sea-swell runup over fringing reefs by super Typhoon Haiyan[J]. Journal of Geophysical Research: Oceans, 2015, 120(6): 4463−4486.

De Alegria-Arzaburu A R, Mariño-Tapia I, Enriquez C, et al. The role of fringing coral reefs on beach morphodynamics[J]. Geomorphology, 2013, 198: 69−83.

Shi Jian, Zhang Chi, Zheng Jinhai, et al. Modelling wave breaking across coral reefs using a non-hydrostatic model[J]. Journal of Coastal Research, 2018, 85(S1): 501−505.

Muñóz-Pérez J J, Tejedor L, Medina R. Equilibrium beach profile model for reef-protected beaches[J]. Journal of Coastal Research, 1999, 15(4): 950−957.

Muñóz-Pérez J J, Medina R. Comparison of long-, medium- and short-term variations of beach profiles with and without submerged geological control[J]. Coastal Engineering, 2010, 57(3): 241−251.

Taaouati M, Parisi P, Passoni G, et al. Influence of a reef flat on beach profiles along the Atlantic Coast of Morocco[J]. Water, 2020, 12(3): 790.

Beach R A, Sternberg R W. Suspended-sediment transport in the surf zone: response to breaking waves[J]. Continental Shelf Research, 1996, 16(15): 1989−2003.

Hughes M G, Masselink G, Brander R W. Flow velocity and sediment transport in the swash zone of a steep beach[J]. Marine Geology, 1997, 138(1/2): 91−103.

Butt T, Russell P. Suspended sediment transport mechanisms in high-energy swash[J]. Marine Geology, 1999, 161(2/4): 361−375.

Ning Yue, Liu Weijie, Sun Zhilin, et al. Parametric study of solitary wave propagation and runup over fringing reefs based on a Boussinesq wave model[J]. Journal of Marine Science and Technology, 2019, 24(2): 512−525.

邵超. 海滩−珊瑚礁海岸侵蚀及适应性管理研究[D]. 厦门: 国家海洋局第三海洋研究所, 2016.

Shao Chao. Study on coastal erosion of beach-coral system and its adaptive management[D]. Xiamen: Third Institute of Oceanography, State Oceanic Administration, 2016.

王宝灿, 陈沈良, 龚文平, 等. 海南岛港湾海岸的形成与演变[M]. 北京: 海洋出版社, 2006.

Wang Baocan, Chen Shenliang, Gong Wenping, et al. Formation and Evolution of Embayment Coast of Hainan Island[M]. Beijing: China Ocean Press, 2006.

Folk R L, Ward W C. Brazos River Bar: a study in the significance of grain size parameters[J]. Journal of Sedimentary Petrology, 1957, 27(1): 3−26.

Udden J A. Mechanical composition of clastic sediments[J]. GSA Bulletin, 1914, 25(1): 655−744.

Shepard F P. Nomenclature based on sand-silt-clay ratios[J]. Journal of Sedimentary Petrology, 1954, 24(3): 151−158.

Shi Fengyan, Kirby J T, Harris J C, et al. A high-order adaptive time-stepping TVD solver for Boussinesq modeling of breaking waves and coastal inundation[J]. Ocean Modelling, 2012, 43−44: 36−51.

邵超, 戚洪帅, 蔡锋, 等. 海滩−珊瑚礁系统风暴响应特征研究——以1409号台风“威马逊”对清澜港海岸影响为例[J]. 海洋学报, 2016, 38(2): 121−130.

Shao Chao, Qi Hongshuai, Cai Feng, et al. Study on storm-effects on beach-coral reef system—Taking the response of Qinglangang Coast on No. 1409 Typhoon Rammasun as an example[J]. Haiyang Xuebao, 2016, 38(2): 121−130.

Demirbilek Z, Nwogu O G, Ward D L. Laboratory study of wind effect on Runup over fringing reefs. Report 1. Data report[R]. Vicksburg: Coastal and Hydraulics Laboratory, 2007.

Reniers A J H M, MacMahan J H, Thornton E B, et al. Modelling infragravity motions on a rip-channel beach[J]. Coastal Engineering, 2006, 53(2/3): 209−222.

da Fontoura Klein A H, de Menezes J T. Beach morphodynamics and profile sequence for a headland bay coast[J]. Journal of Coastal Research, 2001, 17(4): 812−835.

Brander R W, Kench P S, Hart D. Spatial and temporal variations in wave characteristics across a reef platform, Warraber Island, Torres Strait, Australia[J]. Marine Geology, 2004, 207(1/4): 169−184.

Filipot J F, Cheung K F. Spectral wave modeling in fringing reef environments[J]. Coastal Engineering, 2012, 67: 67−79.

Pomeroy A, Lowe R, Symonds G, et al. The dynamics of infragravity wave transformation over a fringing reef[J]. Journal of Geophysical Research: Oceans, 2012, 117(C11): C11022.

Van Dongeren A, Lowe R, Pomeroy A, et al. Numerical modeling of low-frequency wave dynamics over a fringing coral reef[J]. Coastal Engineering, 2013, 73: 178−190.

Ma Gangfeng, Su S F, Liu Shugang, et al. Numerical simulation of infragravity waves in fringing reefs using a shock-capturing non-hydrostatic model[J]. Ocean Engineering, 2014, 85: 54−64.

Buckley M L, Lowe R J, Hansen J E, et al. Mechanisms of wave-driven water level variability on reef-fringed coastlines[J]. Journal of Geophysical Research: Oceans, 2018, 123(5): 3811−3831.

邹志利. 海岸动力学[M]. 4版. 北京: 人民交通出版社, 2009.

Zou Zhili. Coastal Hydrodynamics[M]. 4th ed. Beijing: China Communications Press, 2009.

Larson M, Kraus N C. SBEACH: Numerical model to simulate storm-induced beach change[R]. Army Corps of Engineers, Waterway Experiment, Technical Report CERC-89-9, 1989.

Bodge K R. Representing equilibrium beach profiles with an exponential expression[J]. Journal of Coastal Research, 1992, 8(1): 47−55.

Inman D L, Elwany M H S, Jenkins S A. Shorerise and bar-berm profiles on ocean beaches[J]. Journal of Geophysical Research: Oceans, 1993, 98(C10): 18181−18199.

Lee P Z F. The submarine equilibrium profile: a physical model[J]. Journal of Coastal Research, 1994, 10(1): 1−17.

Kit E, Pelinovsky E. Dynamical models for cross-shore transport and equilibrium bottom profiles[J]. Journal of Waterway, Port, Coastal, and Ocean Engineering, 1998, 124(3): 138−146.

Wang Ping, Davis Jr R A. A beach profile model for a barred coast: case study from sand key, west-central Florida[J]. Journal of Coastal Research, 1998, 14(3): 981−991.

Dubois R N. An inverse relationship between the A and m coefficients in the Bruun/Dean equilibrium profile equation[J]. Journal of Coastal Research, 1999, 15(1): 186−197.

Nelson R C. Depth limited design wave heights in very flat regions[J]. Coastal Engineering, 1994, 23(1/2): 43−59.

2022年第44卷第3期

PDF下载

109

引用本文

BibTeX

文章信息

doi: 10.12284/hyxb2022025

接收时间：2021-03-24
首发时间：2026-02-01
出版时间：2022-03-01

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2021-03-24
修回日期：2021-08-15

基金

自然科学基金重点项目(41930538)

自然科学基金面上项目(42076058)

自然资源部第三海洋研究所基科项目(2019026)

作者信息

1.中山大学海洋工程与技术学院中山大学河口海岸研究所，广东广州 510275

2.自然资源部第三海洋研究所海洋与海岸地质研究室，福建厦门 361005

3.河口水利技术国家地方联合工程实验室，广东广州 510275

4.广东省海岸与岛礁工程技术研究中心，广东珠海 510275

5.自然资源部北部湾滨海湿地生态系统野外科学观测研究站，广西北海 536015

通讯作者:

戚洪帅，男，研究员，主要从事海岸与海滩过程以及海滩养护技术研究。E-mail: qihongshuai@tio.org.cn

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/hyxb/CN/10.12284/hyxb2022025

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT