Impacts of morphological evolution of the Huanghe River mouth by artificial regulation on deltaic sedimentation
-
摘要:
自2002年实施调水调沙以来，由于入海径流、沉积物的通量和组成发生变化，黄河现行河口三角洲叶瓣不断向海淤积造陆，水下三角洲坡度变陡。地形地貌变化对调水调沙期间入海泥沙沉积格局的影响成为黄河口研究的重要科学问题。本文基于Delft3D模型系统，利用岸线、水深和河流水沙数据构建三维水沙数值模型，对2002年、2008年、2014年和2019年调水调沙期间现行河口近岸海域泥沙的输运和沉积过程进行模拟。结果表明，随着水深、岸线的变化，黄河口近岸海域动力环境增强，泥沙的横向输运增强，纵向输运相应减弱；进而导致黄河入海泥沙堆积体的横向长度增加约30%，纵向长度减小约27%，厚度、形态也相应变化。本研究揭示了地形地貌变化条件下，黄河调水调沙期间入海泥沙在河口的沉积格局及动力机制，对深入理解黄河口近岸海域水动力-地貌耦合系统有重要参考价值。
Abstract:Since the water and sediment regulation scheme (WSRS) was implemented in 2002, the present active Huanghe (Yellow) River delta lobe has continuously prograded seaward and the slope of delta has become steeper due to the changes of the river runoff, and riverine sediment flux and components. The impact of morphological evolution of the river mouth on the sedimentation pattern during the WSRS has become a crucial scientific issue. A Delft3D-based three-dimensional hydro-sediment coupling numerical model was established to simulate the transport and sedimentation of riverine sediment in the river mouth during the WSRS conducted in 2002, 2008, 2014, and 2019. Results show that the hydrodynamics in the area were enhanced and the river mouth progressed. Meanwhile, the along-shore transport of sediment was increased while the cross-shore transport was weakened correspondingly. The along-shore extent of the deposition was increased by ～30% while the cross-shore extent was reduced by ～27%, and the thickness and shape of deposition center were changed significantly. This study provided a reference for better understanding the hydrodynamic-morphology coupling system off the Huanghe river mouth.
-
河口是河流与海洋相互作用最强烈的区域，每年接受全球大河携带的约120亿t入海泥沙[1]，其中超过90%的泥沙以悬浮状态向海运输，在径流-潮汐动力相互作用下在河口输运、扩散及沉降，是地球系统中物质循环的关键环节[2]。大量的陆源物质在河口及其近岸海域快速堆积，形成大型三角洲沉积体[3]。在人类活动加剧和全球气候变化的背景下，河口三角洲作为世界上人口最集中、经济最发达的区域之一，入海泥沙的沉积格局和冲淤演化机理已经受到国内外学者的广泛关注[4-8]，是“海岸带陆地-海洋相互作用”[9]、“未来地球海岸”[10]、“海陆边缘科学计划——从源到汇”[11]等重大国际研究计划的重点科学问题。
黄河是中国第二大河，也是全球海洋最重要的陆源物质供应端源之一。黄河以水少沙多、水沙异源著称，年均向渤海输送沉积物量达1.08×109 t[12]，约占全球河流入海沉积物通量的6%[13]，形成了全球造陆最快的三角洲。自20世纪50年代以来，黄河入海水沙通量显著降低[14-15]，尤其是自2002年以来实施的调水调沙工程，显著改变了黄河的入海径流量、沉积物的通量和组成，入海水沙的时空分配也发生了变化[16]。黄河口的冲淤格局产生快速的响应，现行黄河三角洲2002年起由缓慢蚀退转变为快速造陆[17]。同时，调水调沙期间产生的人造洪峰冲刷下游河床，入海粗颗粒泥沙比例增加；入海泥沙粒径增大、沉降速率增加，导致大量的泥沙在河口地区迅速沉降，水下三角洲坡度变陡[18-20]。
随着现行黄河三角洲陆上和水下部分地形地貌的变化，黄河口近岸海域的潮汐、波浪等海洋动力发生了明显的变化[21-24]，河口潮流切变锋的形成位置、传播方式和强度等也发生了调整[25-27]，对入海泥沙在河口的输运和沉积过程产生了影响，深刻改变了黄河口的沉积格局。此外，黄河口泥沙的快速堆积有利于河口拦门沙的形成，河口拦门沙对黄河入海径流起到阻水和分汊的作用，致使尾闾河道易出汊或摆动，使得河口近岸区域的沉积动力环境更为复杂[28-30]，对黄河口近岸海域的沉积演化也有重要影响。因此，探讨调水调沙工程实施以来的新入海水沙情势下，地形地貌变化对黄河口沉积格局的影响有重要的研究价值，对实现陆海统筹、流域-海洋协调发展具有重要现实意义。
本文应用Delft3D模型系统建立的黄河口近岸海域三维水沙动力模型，基于2002、2008、2014和2019年的水深、岸线数据，设置对照实验，模拟调水调沙期间不同地形地貌边界下黄河现行河口近岸海域的沉积动力环境和入海泥沙沉积格局，探究地形地貌变化对黄河口沉积格局的影响。
1. 区域概况
黄河自1855年铜瓦厢决口改道以来，携带大量泥沙在渤海西南部堆积形成现代黄河三角洲，并在自然和人为因素影响下，发生了多次尾闾河道的改道[31]。1996年黄河在清8断面人工出汊改道入海，后逐渐形成了现行黄河三角洲叶瓣（图1）。黄河尾闾摆动频繁，2006年汛期洪水漫滩向北出汊入海，之后于2008年并汊消亡，2013年向东出汊，形成双槽入海的格局延续至今[32]，分汊河道在径流冲刷下不断下蚀加深，流量不对称性显著降低[30]。
![]()
图 1 黄河三角洲遥感影像及2019年近岸海域定点连续观测站位Figure 1. Remote sensing images of the Huanghe River delta and the deployment of in-situ observation sites off the Huanghe river mouth in 2019黄河口属于弱潮河口，外海潮波进入渤海后受地球自转偏向力和渤海地形的影响，一支向北传入辽东湾，一支向南传入莱州湾[33]，形成两大潮波系统，经渤海湾湾顶反射在黄河神仙沟口五号桩和秦皇岛外海形成M2分潮无潮点[34]。潮流以流向平行海岸的往复流为主，涨潮流东南流向，落潮流西北流向，与等深线平行的往复流与沿河口方向的径流共同作用，将黄河口的悬浮泥沙沿西北和东南方向输移[35]。潮流受到无潮点驻波加强，在神仙沟口M2无潮点附近形成高速流区[36]；受到凸出地形的影响，在清水沟老河口附近形成了另一个高速流区[37]。黄河口位于半封闭的渤海内，受到外海传入波浪的影响较小，波浪类型以风浪为主。黄河三角洲夏季盛行偏南风，冬季多为偏北风，主要大浪为秋季开始寒潮南侵造成的偏北向大浪[38]。余流以风生余流为主，春夏季表层余流以东北向为主，秋冬季表层余流以南向为主。底层余流受风的影响很小，表现为由海向陆的上升补偿流。另外，还有与岸线和浅水地形密切相关的潮致余流[39]和汛期黄河入海导致的径流性余流。
2. 数据与方法
2.1 数据资料
本研究使用的岸线基于Landsat卫星遥感影像，通过ENVI计算NDWI提取水边线得到。渤海模型水深来自于海图资料，黄河口近岸水深数据来源于黄河水利委员会历年的实测水深数据。模型验证使用的水文定点观测资料（流速、流向、悬浮泥沙浓度等）来自于课题组2019年在黄河三角洲近岸海域现场观测数据（M1、M2、A30、B30站，图1），流速使用固定于海底三脚架的声学多普勒剖面仪（TRDI Workhorse 600 kHz型）测得，悬浮泥沙浓度通过采集水样并应用抽滤法获得。
2.2 模型设置
Delft3D模型系统是由荷兰Deltares开发，集合水动力、波浪、泥沙输运、地貌、水质、生态、粒子追踪7个开源模块的综合性模拟软件，可以实现河口近岸区域不同尺度的冲淤演变过程模拟和机制分析[40-45]。本研究利用已经完成相关设置和验证的Delft3D模型[30]，对不同年份调水调沙期间黄河口的泥沙输运、沉积过程进行数值模拟。
为提高模型分辨率并保证计算效率，使用渤海模型和黄河口近岸模型进行嵌套。渤海模型采用矩形网格，黄河口近岸嵌套模型采用曲线正交网格，垂向上均采用σ坐标系。模型模拟时间自2019年6月2日至2019年8月21日，根据CFL收敛条件最终确定渤海模型计算步长5 min，黄河口近岸嵌套模型0.5 min。
黄河口近岸海域嵌套模型的河流边界设置在北汊和东汊分汊点上游11 km处，使用利津站实测的日均水沙数据驱动（图2）。其中，6月21日至8月8日期间悬浮沉积物组成由利津水文站采集的水样经粒度测试获得。对于缺少实测数据的时间段，分别使用调水调沙前、后沉积物组成实测值的平均值给定。模型中泥沙组分选择非黏性的砂和黏性的粉砂、黏土三种。非黏性组分沙的干密度设置为
1600 kg/m3，黏性组分干密度设置为500 kg/m3，泥沙的静水沉降速率等参数参考前人设置[41-42]并进行微调，按表1进行设置。河口区域会发生海水和淡水的混合，使得泥沙因絮凝发生加速沉降，细颗粒泥沙的絮凝作用明显[46]。在模型中，泥沙最大絮凝盐度设置为18‰，粉砂的絮凝沉降速率设置为0.5 mm/s，黏土的絮凝沉降速率设置为0.22 mm/s。模型初始床面的沉积物层设置为混合层，层厚设置为统一值20 m，底质沉积物组分含量根据前人[47]研究成果数字化得到。对课题组2015年在黄河三角洲近岸海域采集的表层沉积物样品进行粒度测试得到中值粒径[48]，采用窦国仁公式[49]计算初始床面的临界启动应力，结果如图3所示。![]()
图 2 利津水文站实测黄河日均水沙数据（a）和悬浮沉积物组成（b）Figure 2. Daily river runoff and suspended sediment concentration at gauge station Lijin (a) and components of suspended sediment (b)表 1 模型黏性泥沙和非黏性泥沙参数设置Table 1. The physical parameter settings for cohesive and non-cohesive sediment in the model泥沙类型 泥沙类型 中值粒径/μm 沉降速率/
(mm·s−1)侵蚀速率/
(kg·m−2·s−1)非黏性泥沙 砂 85 − 5.0×10−5 黏性泥沙 粉砂 16 0.12 黏土 11 0.03 ![]()
图 3 黄河口近岸海域表层沉积物中值粒径分布（a）和临界起动应力（b）Figure 3. Distribution of median grain size (a) and critical shear stress (b) of surface sediment off the Huanghe River Mouth2.3 模型验证
经过验证，渤海模型的潮位、振幅和迟角模拟效果良好，可以为嵌套模型提供准确的水动力边界条件，黄河口近岸海域模型的水位、流速、流向模拟效果良好，表明模型可以准确刻画研究区水动力过程[30]。利用M1、M2站位连续观测数据对黄河口近岸海域模型的悬浮泥沙浓度进行验证（图4），模拟值与实测值具有较好的一致性；模拟结果计算得到的相关系数和均方根误差（表2）表明，模拟的悬浮泥沙整体上与实际情况相符。对模型水动力和泥沙的验证结果表明，模型可以模拟黄河口近岸海域的水沙动力过程，满足研究地形地貌变化对河口沉积动力格局影响的要求。
![]()
图 4 M1、M2站位悬浮泥沙观测值（红色点）与模拟值对比（实线）Figure 4. Comparison between observed SSC (suspended sediment concentration) (red point) and modeled SSC (solid line) at Stations M1 and M2表 2 M1、M2站位悬浮泥沙浓度验证结果Table 2. The validation on the suspended sediment concentration at Stations M1 and M2站位 分层 相关系数 均方根误差/(kg/m3) M1 表层 0.90 2.66 中层 0.90 1.73 底层 0.67 1.85 M2 表层 0.76 0.12 中层 0.72 0.18 底层 0.71 0.13 2.4 研究方法
自2002年实施调水调沙工程以来，黄河尾闾河道和河口海岸地区的地貌冲淤格局发生了显著变化[32,50-51]。本研究中采用控制变量法，分别利用2002、2008、2014和2019年实测水深、岸线和河流水沙边界数据构建模型（图5a），对调水调沙期间黄河口动力环境、物质输运和泥沙沉积进行模拟，诊断不同地形地貌条件对黄河口沉积格局的影响。为了研究现行河口垂向上的物质输运和沉积过程，沿现行河口北汊轴线方向设置了断面CQ（图5a），断面CQ水下斜坡及坡度如图5b所示。
![]()
图 5 2002、2008、2014和2019年黄河三角洲岸线（a）和断面CQ水下斜坡坡度变化（b）Figure 5. The shoreline evolution of the Huanghe River delta (a) and slope gradient changes of subaqueous delta (b) along the section CQ in 2002, 2008, 2014, and 20193. 结果
岸线、水深和水下三角洲斜坡的坡度是悬浮泥沙在河口输运、沉积的重要地形边界。为了分析地形地貌变化对黄河口沉积格局的影响，模拟在2002、2008、2014和2019年不同的地形地貌条件下，黄河入海泥沙形成的堆积体在平面、垂向上空间分布的差异。
在4个年份的模拟中，入海泥沙在河口附近沉积形成堆积体的厚度、范围和形态有显著的区别（图6）。2002年，调水调沙期的入海泥沙主要堆积在河口15 m等深线以浅的区域，沉积中心位于5～10 m等深线范围内，最大堆积厚度约6 m；河口右岸堆积体范围明显大于左岸，形成呈北窄南宽的“葫芦型”堆积体（图6a）。2008年，入海泥沙的堆积范围与2002年大致相同，最大堆积厚度约5 m，减少了约1 m；堆积体沿西北-东南向平直分布，长轴大致与5 m等深线走向一致，堆积体整体上呈“扁足状”，堆积体尾部在最南侧发生偏转（图6b）。2014年，随着清8河口向海淤进，入海泥沙沉积的位置向海推进，沉积中心位于10 m等深线附近，堆积体边缘位于15 m等深线上，最大堆积厚度约6 m，呈“水滴状”（图6c）。2019年，入海泥沙在东汊口门附近形成新的沉积中心，同时北汊口门处的沉积中心相应削弱、堆积体厚度减小，两沉积中心均位于5～10 m等深线范围内。在两沉积中心的影响下，入海泥沙在北汊、东汊口门间形成了连续的“弯月状”堆积体，堆积厚度约4～5 m，减小了约1～2 m，边缘位于15 m等深线以浅（图6d）。相较于2002年、2008年和2014年的模拟结果，可以明显观察到2019年的堆积体长轴增长、短轴缩短。同时，观察到河口“马鞍形”分流沙岛不断发育增大，导致北汊河道向北偏转，东汊河道向南偏转。
![]()
图 6 2002、2008、2014及2019年黄河口近岸海域泥沙堆积体Figure 6. Riverine sediment deposition patterns off the Huanghe river mouth in 2002, 2008, 2014, and 2019调水调沙期间黄河入海泥沙主要在现行河口附近发生沉积（图6），在河流输入水沙量相同的前提下，通过统计北汊口门外各深度范围内的泥沙堆积体积和总堆积体积，计算得到各深度范围内泥沙堆积的体积百分比（图7），以此定量分析泥沙堆积体分布范围对河口地形地貌变化的响应。如图7所示，调水调沙期的入海泥沙仅有极少部分沉积于2 m等深线以浅和15 m等深线以深的区域：在4个年份的模拟中，在2 m以浅区域发生沉积的泥沙占比分别为3.6%、1.3%、2.0%和3.1%，15 m以深的区域发生沉积的泥沙占比分别为0.4%、0.1%、1.1%和0.5%，2014年在15 m以深区域发生沉积的入海泥沙量相对较大。2～15 m水深区域是入海悬浮泥沙集中沉积的区域，全部年份均有超过96%的泥沙在该区域堆积（图7）。在4个年份的模拟中，在2～5 m水深区域内发生沉积的泥沙比例分别为22.1%、27.0%、8.8%和11.4%，其中2002年和2008年堆积体积是2014年和2019年的2～3倍；5～10 m水深区域发生沉积的泥沙比例分别为48.6%、60.7%、47.3%、54.1%，是入海泥沙最主要的沉积区域；10～15 m水深区域发生沉积的泥沙比例分别为25.4%、10.9%、40.8%、30.9%，其中2014年体积百分数分别为2002和2008年的1.6倍和3.7倍。
![]()
图 7 黄河口近岸海域不同深度范围内泥沙堆积体积百分比Figure 7. The volume percentage of riverine sediment accumulation at different water depth off the Huanghe river mouth除了在平面上存在明显的空间差异，河口断面CQ的模拟结果显示，随着河口水下三角洲斜坡坡度的增加（图5b），入海泥沙堆积体的位置、厚度和形态，以及堆积体向海推进的距离，在垂向上也存在着明显的差异（图8）。调水调沙期间河口水下斜坡沉积地层具有明显的斜层理，并在受水动力的侵蚀、破坏下，在沉积体顶部形成了不整合的薄层堆积体（图8）。2002年，河口水下斜坡向海推进最远距离约4 km，最大淤积厚度约5.6 m（图8a）；2008年，河口水下斜坡向海推进的最远距离约2.4 km，最大淤积厚度约6 m（图8b）；这两个年份入海泥沙堆积在12 m以浅的区域，堆积体向海延伸至距河口8～10 km处尖灭（图8a-b）。随着水下斜坡坡度的增加，水下斜坡向海推进距离和堆积厚度均减小。2014年，河口水下斜坡向海最远推进距离约1.3 km，最大沉积厚度约为5 m（图8c）；2019年，河口水下斜坡向海推进最远距离约0.8 km，最大沉积厚度约为5 m（图8d），这两个年份入海泥沙堆积在15 m以浅的区域，堆积体向海延伸至距河口6～8 km处尖灭（图8c-d）。随着黄河东汊河道分流作用增强，北汊口门处的沉积中心削弱（图6），2019年水下斜坡最大推进距离相比于2014年减小了约0.3 km。
![]()
图 8 2002、2008、2014和2019年CQ断面的泥沙堆积体Figure 8. Riverine sediment deposition buildup at section CQ in 2002, 2008, 2014, and 20194. 讨论
4.1 河口物质输运变化对入海泥沙沉积格局的影响
余流是影响河口物质输运的重要因素，黄河三角洲近岸海域的余流以风生余流、潮致余流以及调水调沙期间入海径流引起的余流为主，包括入海径流射流效应引起的扩散流和淡水斜压效应引起的斜压流[52]。黄河口三角洲是弱潮河控型三角洲，调水调沙期间黄河口近岸海域的风速通常很小（2～4 m/s），大风天气频率很低，径流和潮汐是塑造河口地形及控制泥沙运移的主要动力[42,53] ，因此模型中没有考虑风和波浪[53-54] 。本文重点研究了径流与潮汐相互作用下河口附近海域余流场的变化时，河口泥沙输运与泥沙堆积体的关系。
径流引起的余流主要发生在河口近岸海域，是汛期的重要余流类型[55-57]。调水调沙期间，黄河在大约20天中输送了全年14%～50%的径流[16]。受入海径流的驱动，河口附近表层余流方向由陆向海，底层余流方向由海向陆[52]。河流冲淡水在口门附近向两河口侧输运，一支北上受渤海中部海域环流挤压向东偏转[55]，另一支在科氏力作用下发生偏转[58-59]，向南进入莱州湾，引起莱州湾盐度降低。
在黄河口近岸海域，径流引起的余流主要受控于汛期的入海径流，而潮致余流主要受控于地形地貌的变化[60-61]。随着地形地貌的变化，在潮汐动力环境变化和径流的共同作用下，黄河口调水调沙期间的物质输运格局变化显著，最终控制了入海泥沙的沉积格局。根据模拟结果，按照下式计算单宽泥沙通量F，得到泥沙单宽余通量变化分布（图9），以研究地形地貌变化条件下入海泥沙输运格局的响应：
![]()
图 9 黄河口近岸海域表层泥沙单宽余通量变化黑色箭头指示泥沙的输运方向。Figure 9. Variation in the net transport rate of riverine sediment in surface layer off the Huanghe river mouthBlack arrows indicate the direction of sediment transport.$$ F=\frac{1}{T}{\int }_{0}^{T}\overline{V}(x,y,z,t)\cdot C(x,y,z,t)\mathrm{d}t $$ (1) 式中，$ T $为选择的计算潮周期，选择的时间为7月15—31日（调沙期），$ \overline{V}(x,y,z,t) $表示点$ (x,y,z) $在$ t $时刻的矢量速度，$ C(x,y,z,t) $表示$ t $时刻的泥沙浓度。
现行河口表层存在单宽泥沙余通量大于1 kg·m−1·s−1的区域，代表泥沙输运高值区域，区域形态近似椭圆形，最大悬浮泥沙输运量约为20～32 kg·m−1·s−1（图9）。在河口轴线方向上，表层入海泥沙的单宽余通量由陆向海快速减小，到15 m等深线附近减小至0.001 kg·m−1·s−1；平面上，受河口附近余流的影响，表层悬浮泥沙在河口以北以西北向输运，在河口以南以东南向输运（图9）。随着调水调沙以来黄河入海悬浮泥沙浓度的降低，在4个年份的模拟中，河口泥沙扩散方式均以羽状流为主[20]。在现行河口底层，余流以补偿流性质的上升流为主[52]，入海泥沙受到由海向陆余流的顶托，底层悬浮泥沙输运量很小。入海泥沙在河口表层以羽状流形式输运易受到潮流切变锋的阻隔与捕集[26,62-63]，而在河口底层余流以补偿流性质的上升流为主[52]，泥沙在余流的顶托下输运量很小，导致黄河水沙向深水区扩散的范围有限。2002年、2008年时泥沙输运高值区域（单宽泥沙余通量＞1 kg·m−1·s−1的近“椭圆形”区域）的外海边缘大致位于10 m等深线处（图9a-b），泥沙输运低值区域（单宽泥沙余通量＞0.001 kg·m−1·s−1的近“椭圆形”区域）的外海边缘大致位于15 m等深线处。而随着现行河口向海淤进，现行河口表层的泥沙更容易向深水区输运并沉积，2014年、2019年时泥沙通量高值区域的外海边缘已经移动至10 m等深线以深的区域，泥沙输运低值区域的外海边缘已经超出15 m等深线（图9c-d）。
对模拟结果中泥沙输运高值区域、泥沙输运低值区域和入海泥沙堆积体的长短轴长度进行统计（图10），以阐明河口物质输运与沉积格局的关系。统计结果表明，泥沙输运高值和低值区域的长轴整体上均呈增大趋势（图10a-b），而短轴持续减小。表明随着河口的淤进和水深岸线的变化，河口表层悬浮泥沙沿等深线向河口两侧的输运增强，而沿河口轴线方向的输运相对减弱。2019年，随着东汊河道的发育，分流分沙能力增强，东汊口门附近的单宽泥沙通量由2014年的0.2 kg·m−1·s−1增加到2019年的1.0 kg·m−1·s−1。汊道流量不对称性显著降低，北汊口门附近入海泥沙量减小，导致2019年泥沙输运高值和低值区域的长、短轴长同步减小（图10）。
![]()
图 10 泥沙输运高值区域（＞1 kg·m−1·s−1）(a)、泥沙输运低值区域（＞0.001 kg·m−1·s−1) (b)和泥沙堆积体(c)的长、短轴变化Figure 10. Changes of major and minor axes of area with F＞1 kg·m−1·s−1 (a), F＞0.001 kg·m−1·s−1 (b), and riverine sediment deposition accumulation (c)入海泥沙的沉积过程受控于河口泥沙的输运过程，随着泥沙在河口横向输运的增强和纵向的输运减弱（图10a-b），河口泥沙堆积体的长轴呈增大趋势，短轴出现持续的减小（图10c）。与2002年相比，2019年时河口泥沙堆积体的横向长度增加约30%，纵向长度减少约27%，河口泥沙输运和沉积的变化呈现较好的一致性（图10）。在河口泥沙输运格局变化的影响下，河口泥沙堆积体的范围、厚度和形态对河口泥沙输运格局的变化产生了明显的响应。2002年，由于现行河口地形的挑流作用，入海泥沙在河口以北向西北输运，在河口以南向东南输运，在河口发生转向，最终形成了“U”型包围河口的“葫芦形”堆积体（图6a）。2008年，随着2006—2008年河口附近相继发生的出汊、并汊归股，河口向北偏转，入海沉积物可被搬运的距离增大，导致河口堆积体范围扩大，堆积厚度减小了约1 m，长轴长度减小（图6b）。同时，由于三角洲叶瓣东侧向海淤积凸出造陆，现行河口附近岸线变得平直，入海泥沙在河口附近形成了西北-东南向的平直分布的“扁足状”堆积体，堆积体尾部在潮流作用下向南发生偏转（图6b）。2014年时，随着现行河口向海淤进，河口至15 m等深线的距离缩短，水下三角洲斜坡坡度增大；且东汊河道分流分沙微弱，北汊是黄河水沙入海的主要通道[30]，入海泥沙更容易扩散到深水区。因此，2014年时10～15 m水深区域沉积的泥沙体积百分比达到40.8%，分别为2002和2008年的1.6倍和3.7倍，15～20 m水深沉积的泥沙体积百分比分别为2002和2008年的2.7倍和11.0倍。泥沙通量大于0.001 kg·m−1·s−1的区域明显超出15 m等深线（图9c），泥沙堆积体堆积范围扩大，边缘大致位于15 m等深线上，整体呈“水滴状”（图6c）。2019年时，汊道流量不对称性显著降低，现行河口北汊和东汊口门附近均形成了泥沙通量高值区（图9d）和沉积中心（图6d），随着河口往复流流速增强，两口门间形成了连续的“弯月状”堆积体。堆积体厚度约4 m，平行于岸线呈直线分布，长轴方向明显增长，短轴方向明显缩短（图6d）。由于北汊的行水能力相对减弱，泥沙扩散范围减小，北汊口门外10～20 m水深范围内沉积的泥沙体积百分比较2014年降低了约10%。
4.2 河口水下斜坡变化对入海泥沙沉积格局的影响
前人通过陆上三角洲和水下三角洲两个部分对河口三角洲的地形地貌变化进行了研究，其中水下三角洲的冲淤对三角洲地貌演化和海洋动力作用响应更加敏感[64-68]。随着河口不断向海淤进，在模拟的4个年份中，断面CQ水下斜坡坡度存在显著差异。2002、2008年河口水下斜坡的坡度较缓，分别为约1.2‰、1.0‰，容沙空间充足，水下斜坡为“宽缓”型；2014、2019年河口水下斜坡的坡度较陡，分别为约1.3‰、1.4‰，容沙空间较小，水下斜坡为“窄陡”型（图5）。不同类型水下斜坡的河口环流和高速流场的空间分布存在差异，进而显著影响了入海泥沙在河口的沉积格局。
当河口盐淡水产生纵向盐度梯度与河口纵向坡降相互作用，上层水体被驱动向海，而底层水体向陆便产生了河口环流[69]，其在盐度、悬沙和营养盐等物质输运过程中起到非常重要的作用[70]。根据模拟结果绘制了4个年份的河口轴线余流、平均悬浮泥沙浓度和盐度分布图（图11）。在现行河口附近，表层余流方向由陆向海，底层余流方向由海向陆，在一定水深处发育垂向环流。在不同类型水下斜坡的影响下，河口环流作用深度存在差异。在2002年、2008年“宽缓”型的水下斜坡影响下，河口垂向环流作用于8～10 m深度水域，距离口门约4～6 km（图11a, b）。而在2014年、2019年“窄陡”型的水下斜坡影响下，河口垂向环流作用于10～12 m深度水域，距离口门约2～4 km（图11c-d）。河口环流导致水体的层化结构增强，加快了泥沙沉降，阻隔泥沙向深水区输运，从而控制泥沙的沉积格局[71]。以2002年和2008年为例，2002年垂向环流陆侧的平均悬浮泥沙浓度约为2.0 kg/m3，海侧的平均悬浮泥沙浓度约为0.5 kg/m3，垂向环流阻隔了75%的悬浮泥沙；2008年垂向环流陆侧的平均悬浮泥沙浓度约为3.6 kg/m3，海侧的平均悬浮泥沙浓度约为1.2 kg/m3，垂向环流阻隔了67%的悬浮泥沙（图11a-b）。
![]()
图 11 河口CQ断面的余流、盐度和平均悬浮泥沙浓度黑色和绿色箭头使用不同比例尺表示余流大小，红色箭头示意河口环流的方向。Figure 11. Residual current, salinity, and suspended sediment concentration at section CQThe black and green arrows indicate the value of the residual flow on different scales,the red arrows indicate the direction of the estuarine circulation.研究表明，口门地貌形态改变会影响河口近岸海域的潮流场分布[72]。本研究在不同的水下斜坡类型条件下的模拟显示，在河口轴线方向上，高速流场的空间分布存在差异，反映了河口不同的水动力环境。2002、2008年“宽缓”型的河口水下斜坡水深坡度小，高速流场距口门远，在涨急时刻距离河口超过2 km，落急时刻距离河口超过4 km，最大涨急流速可达约0.8 m/s（图12），表明河口附近水动力较弱。而2014、2019年“窄陡”型的河口水下斜坡水深、坡度较大，高速流场紧邻口门，涨急时刻距河口2 km以内，落急时刻分别距河口约2 km、3.5 km，最大涨急流速可达约1.0 m/s（图12），表明河口附近水动力较强。
![]()
图 12 2002、2008、2014和2019年断面CQ(a)涨急和(b)落急流速分布Figure 12. The maximum flood current velocity (a) and maximum ebb current velocity (b) distributions along section CQ in 2002, 2008, 2014, and 2019高速流场与河口环流两个沉积动力要素在河口有较好的协同性。“宽缓”型的河口水下斜坡的坡度小、水深小，高速流场距离口门远，近岸浅水缓坡区域的水动力较弱（图12），河口环流向陆移动、作用水深小（图11），其对泥沙的捕获、阻隔作用有利于入海泥沙在河口附近的浅水区域发生沉积，有利于河口向海推进。在“宽缓”型水下斜坡的影响下，2002年和2008年在2～5 m水深区域发生沉积的泥沙占比可达22.1%和27%（图7）。“窄陡”型的河口水下斜坡坡度陡、水深大，高速流场靠近口门，近岸水动力较强（图12）；随着斜坡坡度增加，河口环流向海移动，作用水深大（图11），有利于泥沙向深水区扩散、沉积，并在较强的往复流作用下平行于岸线发生横向输运，减缓了河口向海推进速度，有利于流路的稳定。在“窄陡”型水下斜坡的影响下，2014年和2019年同区域发生沉积的泥沙百分比分别为8.8%和11.4%（图7），仅为2002年和2008年在2～5 m水深区域沉积泥沙的50%左右。
5. 结论
本文基于Delft3D模型系统建立的黄河口近岸海域的三维水沙数值模型，对2002、2008、2014和2019年地形地貌变化条件下，调水调沙期间黄河口的物质输运和沉积过程进行模拟，主要结论如下：
（1）黄河调水调沙期间的入海泥沙仅极少部分在2 m以浅和15 m以深区域发生沉积，96%以上的泥沙沉积于2～15 m水深区域，其中5～10 m水深区域是发生沉积的最主要区域。随着现行河口向海淤进，河口近岸区域地形地貌变化明显，入海泥沙发生沉积的区域向深水区移动。在不同的地形边界下，入海泥沙形成堆积体的堆积范围、厚度和形态存在显著差异，表明地形地貌影响了入海泥沙的沉积格局。
（2）随着地形地貌的变化，调水调沙期间黄河口的物质输送格局亦发生显著改变，最终控制了入海泥沙的沉积格局。在潮致余流和入海径流引起的扩散流、斜压流的共同作用下，调水调沙期间黄河口近岸海域泥沙的横向输运增强、纵向输运减弱，导致入海泥沙堆积体横向长度增加约30%，纵向长度减少约27%，堆积体的堆积范围、厚度和形态明显改变。此外，随着分汊河道的发育，流量不对称系数减小，东汊、北汊口门附近均形成了泥沙输运和沉积中心，两口门间形成了连续的“弯月状”沉积体。
（3）2002、2008年河口水下斜坡的坡度小、水深浅、容沙空间小，属于“宽缓”型水下斜坡，河口水动力较弱，垂向环流靠近口门，泥沙易在近岸缓坡上堆积，有利于河口沙嘴向海推进；随着河口地形地貌的变化，2014、2019年河口水下斜坡的坡度大、水深深、容沙空间大，属于“窄陡”型水下斜坡，近岸区域水动力强，垂向环流向海移动，泥沙不易在近岸缓坡发生沉积。随着河口往复流增强，泥沙发生横向输运和堆积，有利于流路的稳定。
致谢：感谢中国海洋大学海洋高等研究院海洋大数据中心提供的数据计算、存储服务和董玉杰提供的专业技术支持。
-
![]()
图 1 黄河三角洲遥感影像及2019年近岸海域定点连续观测站位
Figure 1. Remote sensing images of the Huanghe River delta and the deployment of in-situ observation sites off the Huanghe river mouth in 2019
![]()
图 2 利津水文站实测黄河日均水沙数据（a）和悬浮沉积物组成（b）
Figure 2. Daily river runoff and suspended sediment concentration at gauge station Lijin (a) and components of suspended sediment (b)
![]()
图 3 黄河口近岸海域表层沉积物中值粒径分布（a）和临界起动应力（b）
Figure 3. Distribution of median grain size (a) and critical shear stress (b) of surface sediment off the Huanghe River Mouth
![]()
图 4 M1、M2站位悬浮泥沙观测值（红色点）与模拟值对比（实线）
Figure 4. Comparison between observed SSC (suspended sediment concentration) (red point) and modeled SSC (solid line) at Stations M1 and M2
![]()
图 5 2002、2008、2014和2019年黄河三角洲岸线（a）和断面CQ水下斜坡坡度变化（b）
Figure 5. The shoreline evolution of the Huanghe River delta (a) and slope gradient changes of subaqueous delta (b) along the section CQ in 2002, 2008, 2014, and 2019
![]()
图 6 2002、2008、2014及2019年黄河口近岸海域泥沙堆积体
Figure 6. Riverine sediment deposition patterns off the Huanghe river mouth in 2002, 2008, 2014, and 2019
![]()
图 7 黄河口近岸海域不同深度范围内泥沙堆积体积百分比
Figure 7. The volume percentage of riverine sediment accumulation at different water depth off the Huanghe river mouth
![]()
图 8 2002、2008、2014和2019年CQ断面的泥沙堆积体
Figure 8. Riverine sediment deposition buildup at section CQ in 2002, 2008, 2014, and 2019
![]()
图 9 黄河口近岸海域表层泥沙单宽余通量变化
黑色箭头指示泥沙的输运方向。
Figure 9. Variation in the net transport rate of riverine sediment in surface layer off the Huanghe river mouth
Black arrows indicate the direction of sediment transport.
![]()
图 10 泥沙输运高值区域（＞1 kg·m−1·s−1）(a)、泥沙输运低值区域（＞0.001 kg·m−1·s−1) (b)和泥沙堆积体(c)的长、短轴变化
Figure 10. Changes of major and minor axes of area with F＞1 kg·m−1·s−1 (a), F＞0.001 kg·m−1·s−1 (b), and riverine sediment deposition accumulation (c)
![]()
图 11 河口CQ断面的余流、盐度和平均悬浮泥沙浓度
黑色和绿色箭头使用不同比例尺表示余流大小，红色箭头示意河口环流的方向。
Figure 11. Residual current, salinity, and suspended sediment concentration at section CQ
The black and green arrows indicate the value of the residual flow on different scales,the red arrows indicate the direction of the estuarine circulation.
![]()
图 12 2002、2008、2014和2019年断面CQ(a)涨急和(b)落急流速分布
Figure 12. The maximum flood current velocity (a) and maximum ebb current velocity (b) distributions along section CQ in 2002, 2008, 2014, and 2019
表 1 模型黏性泥沙和非黏性泥沙参数设置
Table 1 The physical parameter settings for cohesive and non-cohesive sediment in the model
泥沙类型 泥沙类型 中值粒径/μm 沉降速率/
(mm·s−1)侵蚀速率/
(kg·m−2·s−1)非黏性泥沙 砂 85 − 5.0×10−5 黏性泥沙 粉砂 16 0.12 黏土 11 0.03 
下载: 导出CSV
表 2 M1、M2站位悬浮泥沙浓度验证结果
Table 2 The validation on the suspended sediment concentration at Stations M1 and M2
站位 分层 相关系数 均方根误差/(kg/m3) M1 表层 0.90 2.66 中层 0.90 1.73 底层 0.67 1.85 M2 表层 0.76 0.12 中层 0.72 0.18 底层 0.71 0.13
下载: 导出CSV
-
[1] Syvitski J P M, Voeroesmarty C J, Kettner A J, et al. Impact of Humans on the Flux of Terrestrial Sediment to the Global Coastal Ocean[J]. Science, 2005, 308(5720):376-380. doi: 10.1126/science.1109454
[2] Syvitski J P M, Saito Y. Morphodynamics of deltas under the influence of humans[J]. Global and Planetary Change, 2007, 57(3):261-282.
[3] Wright L D. Sediment transport and deposition at river mouths: A synthesis[J]. Geological Society of America Bulletin, 1977, 88(6):857-868. doi: 10.1130/0016-7606(1977)88<857:STADAR>2.0.CO;2
[4] Bi N S, Wang H J, Yang Z S. Recent changes in the erosion-accretion patterns of the active Huanghe (Yellow River) delta lobe caused by human activities[J]. Continental Shelf Research, 2014, 90:70-78.
[5] Edmonds D A, Slingerland R L. Significant effect of sediment cohesion on delta morphology[J]. Nature Geoscience, 2009, 3(2):105-109.
[6] Nittrouer J A, Best J L, Brantley C, et al. Mitigating land loss in coastal Louisiana by controlled diversion of Mississippi River sand[J]. Nature Geoscience, 2012, 5(8):534-537. doi: 10.1038/ngeo1525
[7] Nittrouer J A, Viparelli E. Sand as a stable and sustainable resource for nourishing the Mississippi River delta[J]. Nature Geoscience, 2014, 7(5):350-354. doi: 10.1038/ngeo2142
[8] Temmerman B S, Kirwan M L, Karina H. Building land with a rising sea[J]. Science, 2015, 349(6248):9-11.
[9] Loicz IPO[Z]. Land-ocean interactions in the coastal zone: Science plan and implementation strategy, 2005.
[10] Future Eearth Coasts IPO[Z]. Strategy for research 2018–2028, 2018.
[11] Margins Office[Z]. NSF margins program science plans, 2003.
[12] Milliman J D, Syvitski J P M. Geomorphic/Tectonic Control of Sediment Discharge to the Ocean: The Importance of Small Mountainous Rivers[J]. Journal of Geology, 1992, 100(5):525-544. doi: 10.1086/629606
[13] Milliman J D, Meade R H. World-Wide Delivery of River Sediment to the Oceans[J]. Journal of Geology, 1983, 91(1):1-21. doi: 10.1086/628741
[14] Wang H J, Yang Z S, Saito Y, et al. Interannual and seasonal variation of the Huanghe (Yellow River) water discharge over the past 50 years: Connections to impacts from ENSO events and dams[J]. Global and Planetary Change, 2006, 50:212-225. doi: 10.1016/j.gloplacha.2006.01.005
[15] Wang H J, Yang Z S, Saito Y, et al. Stepwise decreases of the Huanghe (Yellow River) sediment load (1950-2005): Impacts of climate change and human activities[J]. Global and Planetary Change, 2007, 57:331-354. doi: 10.1016/j.gloplacha.2007.01.003
[16] Wang H J, Wu X, Bi N S, et al. Impacts of the dam-orientated water-sediment regulation scheme on the lower reaches and delta of the Yellow River (Huanghe): A review[J]. Global and Planetary Change, 2017, 157:93-113. doi: 10.1016/j.gloplacha.2017.08.005
[17] Wu X, Wang H J, Bi N S, et al. Impact of artificial floods on the quantity and grain size of river-borne sediment: A case study of a dam regulation scheme in the Yellow River catchment[J]. Water Resources Research, 2021, 57:e2021WR029581.
[18] Bi N S, Sun Z Q, Wang H J, et al. Response of channel scouring and deposition to the regulation of large reservoirs: A case study of the lower reaches of the Yellow River (Huanghe)[J]. Journal of Hydrology, 2019, 568:972-984. doi: 10.1016/j.jhydrol.2018.11.039
[19] Wu X, Bi N S, Kanai Y, et al. Sedimentary records off the modern Huanghe (Yellow River) delta and their response to deltaic river channel shifts over the last 200 years[J]. Journal of Asian Earth Sciences, 2015, 108:68-80. doi: 10.1016/j.jseaes.2015.04.028
[20] Wang H J, Bi N S, Saito Y, et al. Recent changes in sediment delivery by the Huanghe (Yellow River) to the sea: Causes and environmental implications in its estuary[J]. Journal of Hydrology, 2010, 391(3-4):302-313. doi: 10.1016/j.jhydrol.2010.07.030
[21] Wang N, Li G X, Xu J S, et al. The marine dynamics and changing trend off the modern Yellow River Mouth[J]. Journal of Ocean University of China, 2015, 14(3):433-445. doi: 10.1007/s11802-015-2764-0
[22] 王永刚, 魏泽勋, 方国洪, 等. 黄河口及其邻近海域水深和岸线变化对M2分潮影响的数值研究[J]. 海洋科学进展, 2014, 32(2):141-147 doi: 10.3969/j.issn.1671-6647.2014.02.003 WANG Yonggang, WEI Zexun, Fang Guohong, et al. A numerical study on the effect of changes in water depth and coastline on M2 tidal component near the Yellow River Estuary[J]. Advances in Marine Science, 2014, 32(2):141-147.] doi: 10.3969/j.issn.1671-6647.2014.02.003
[23] Zhu L H, Hu R J, Zhu H J, et al. Modeling studies of tidal dynamics and the associated responses to coastline changes in the Bohai Sea, China[J]. Ocean Dynamics, 2018, 68:1625-1648. doi: 10.1007/s10236-018-1212-2
[24] Rafael J B, Alejandro L, Miguel O. Implications of delta retreat on wave propagation and longshore sediment transport - Guadalfeo case study (southern Spain)[J]. Marine Geology, 2016, 382:1-16. doi: 10.1016/j.margeo.2016.09.011
[25] 卢昱岑, 沈永明, 张明. 地形演变对黄河口切变锋位置及盐度分布的影响[J]. 水动力学研究与进展, 2012, 27(3):348-358 LU Yucen, SHEN Yongming, ZHANG Ming. Influence of topography evolution on position of tidal shear front and distribution of salinity around Yellow River estuary[J]. Chinese Journal of Hydrodynamics, 2012, 27(3):348-358.]
[26] Wang N, Li G X, Qiao L L, et al. Long-term evolution in the location, propagation, and magnitude of the tidal shear front off the Yellow River Mouth[J]. Continental Shelf Research, 2017, 137:1-12. doi: 10.1016/j.csr.2017.01.020
[27] Wang N, Li K, Song D H, et al. Impact of tidal shear fronts on terrigenous sediment transport in the Yellow River Mouth: Observations and a synthesis[J]. Marine Geology, 2024, 469:107222. doi: 10.1016/j.margeo.2024.107222
[28] 王厚杰, 杨作升, 毕乃双, 等. 2005年黄河调水调沙期间河口入海主流的快速摆动[J]. 科学通报, 2005, 50(23):2656-2662 doi: 10.3321/j.issn:0023-074X.2005.23.016 WANG Houjie, YANG Zuosheng, BI Naishuang, et al. Dispersal pattern of suspended sediment in the shear frontal zone off the Huanghe (Yellow River) mouth[J]. Chinese Science Bulletin, 2005, 50(23):2656-2662.] doi: 10.3321/j.issn:0023-074X.2005.23.016
[29] 徐丛亮, 谷硕, 刘喆, 等. 黄河调水调沙14a来河口拦门沙形态变化特征[J]. 人民黄河, 2016, 38(10):69-73 doi: 10.3969/j.issn.1000-1379.2016.10.014 XU Congliang, GU Shuo, LIU Zhe, et al. Characteristic of the river mouth bar in the past 14 years of the Yellow River Water-Sediment Regulation[J]. Yellow River, 2016, 38(10):69-73.] doi: 10.3969/j.issn.1000-1379.2016.10.014
[30] 葛春海, 范勇勇, 巴旗, 等. 现行黄河口分汊河道的分流特征及其影响机制[J]. 海洋地质与第四纪地质, 2024, 44(2):131-145 GE Chunhai, FAN Yongyong, BA Qi, et al. Diversion characteristics of the branching channels in the Yellow River mouth and its influencing mechanisms[J]. Marine Geology & Quaternary Geology, 2024, 44(2):131-145.]
[31] Fan H, Huang H J, Zeng T Q, el al. River mouth bar formation, riverbed aggradation and channel migration in the modern Huanghe (Yellow River) delta, China[J]. Geomorphology, 2006, 74(1):124-136.
[32] 杨卓媛, 夏军强, 周美蓉, 等. 黄河口尾闾河道近期自然出汊过程及其机理探讨[J]. 泥沙研究, 2022, 47(1):65-72 YANG Zhuoyuan, XIA Junqiang, ZHOU Meirong, et al. Study on channel avulsion in the recent tail reach of the Yellow River Estuary[J]. Journal of Sediment Research, 2022, 47(1):65-72.]
[33] 胡春宏, 曹文洪. 黄河口水沙变异与调控Ⅰ——黄河口水沙运动与演变基本规律[J]. 泥沙研究, 2003(5):1-8 doi: 10.3321/j.issn:0468-155X.2003.05.001 HU Chunhong, CAO Wenhong. Variation, Regulation and Control of Flow and Sediment in the Yellow River Estuary I: Mechanism of Flow-Sediment Transport and Evolution[J]. Journal of Sediment Research, 2003(5):1-8.] doi: 10.3321/j.issn:0468-155X.2003.05.001
[34] 侍茂崇, 赵进平. 黄河三角洲半日潮无潮区位置及水文特征分析[J]. 山东海洋学院学报, 1985, 15(1):127-136 SHI Maochong, ZHAO Jinping. The analysis of hydrographical characteristics in the nontidal region M2 near the delta of the Huanghe River[J]. Journal of Shandong College of Oceanology, 1985, 15(1):127-136.]
[35] 姬泓宇. 新入海水沙情势下黄河三角洲地貌动态变化与演变机制[D]. 华东师范大学博士学位论文, 2021 JI Hongyu. Morphological variability of the Yellow River Delta and its dynamic mechanism under the new regime of river delivery[D]. Doctor dissertation of East China Normal University, 2021.]
[36] 王宝灿, 黄仰松. 海岸动力地貌[M]. 上海: 华东师范大学出版社, 1989 WANG Baocan, HUANG Yangsong. Coastal Dynamics[M]. Shanghai: East China Normal University Press, 1989.]
[37] Wang H J, Yang Z S, Li G X, et al. Wave climate modeling on the abandoned Huanghe (Yellow River) delta lobe and related deltaic erosion[J]. Journal of Coastal Research, 2006, 224(4):906-918.
[38] 藏启运. 黄河三角洲近岸泥沙[M]. 北京: 海洋出版社, 1996: 34-42 ZANG Qiyun, Nearshore Sediment of the Huanghe River and Its Delta[M]. Beijing: China Ocean Press, 1996: 34-42.]
[39] Qiao L L, Bao X W, Wu D X, et al. Numerical study of generation of the tidal shear front off the Yellow River mouth[J]. Continental Shelf Research, 2008, 28(14):1782-1790. doi: 10.1016/j.csr.2008.04.007
[40] Fagherazzi S, Edmonds D A, Nardin W, et al. Dynamics of river mouth deposits[J]. Reviews of Geophysics, 53(3): 642-672.
[41] 刘猛, 毕乃双, 纪金龙, 等. 现行黄河三角洲叶瓣蚀积演化对动力环境的影响[J]. 海洋地质前沿, 2018, 34(6):8-18 LIU Meng, BI Naishuang, JI Jinlong, et al. Evolution of the active deltaic lobe of Huanghe River and its response to hydrodynamics[J]. Marine Geology Frontiers, 2018, 34(6):8-18.]
[42] 凡姚申. 黄河三角洲近岸海床侵蚀过程及其动力机制[D]. 华东师范大学博士学位论文, 2019 FAN Yaoshen. Seabed erosion and its mechanism in the littoral area of Yellow River Delta[D]. Doctor dissertation of East China Normal University, 2019.]
[43] Ji H Y, Pan S Q, Chen S L, et al. Impact of river discharge on hydrodynamics and sedimentary processes at Yellow River Delta[J]. Marine Geology, 2020, 425:106210. doi: 10.1016/j.margeo.2020.106210
[44] 乔璐璐, Le D, 李珏, 等. 超强台风“威马逊”作用下红河三角洲海域水动力环境变化的数值研究[J]. 海洋科学, 2021, 45(4):64-74 doi: 10.11759/hykx20190508002 QIAO Lulu, LE DUC Cuong, LI Yu, et al. Numerical modeling of hydrodynamic changes due to super Typhoon Rammasun in the Red River Delta coastal area[J]. Marine Sciences, 2021, 45(4):64-74.] doi: 10.11759/hykx20190508002
[45] Umgiesser G, Ferrarin C, Bajo M, et al. Hydrodynamic modelling in marginal and coastal seas — The case of the Adriatic Sea as a permanent laboratory for numerical approach[J]. Ocean Modelling, 2022, 179:102123. doi: 10.1016/j.ocemod.2022.102123
[46] Gratiot N, Bildstein A, Anh T T, et al. Sediment flocculation in the Mekong River estuary, Vietnam, an important driver of geomorphological changes[J]. Comptes Rendus Géoscience, 2017, 349(6-7):260-268.
[47] 季有俊. 渤海海域泥沙输运对季节性因素及地形变化响应的数值模拟研究[D]. 中国海洋大学博士学位论文, 2010 JI Youjun. Numerical study on response of suspended sediment transport to the changes of seasonal factors and topography in the Bohai Sea[D]. Doctor dissertation of Ocean University of China, 2010.]
[48] Liu L, Wang H J, Yang Z S, et al. Coarsening of sediments from the Huanghe (Yellow River) delta-coast and its environmental implications[J]. Geomorphology, 401: 108105.
[49] 窦国仁. 河口海岸全沙模型相似理论[J]. 水利水运工程学报, 2001(1):1-12 doi: 10.3969/j.issn.1009-640X.2001.01.001 DOU Guoren. Similarity theory of total sediment transport modeling for estuarine and coastal regions[J]. Hydro-Science and Engineering, 2001(1):1-12.] doi: 10.3969/j.issn.1009-640X.2001.01.001
[50] Zheng S, Wu B S, Wang K R, et al. Evolution of the Yellow River delta, China: Impacts of channel avulsion and progradation[J]. International Journal of Sediment Research, 2017, 32(1):34-44. doi: 10.1016/j.ijsrc.2016.10.001
[51] 刘清兰, 陈俊卿, 陈沈良. 调水调沙以来黄河尾闾河道冲淤演变及其影响因素[J]. 地理学报, 2021, 76(1):139-152 doi: 10.11821/dlxb202101011 LIU Qinglan, CHEN Junqing, CHEN Shenliang. Spatiotemporal evolution of Yellow River estuarine channel and its influencing factors since the water-sediment regulation scheme[J]. Acta Geographica Sinica, 2021, 76(1):139-152.] doi: 10.11821/dlxb202101011
[52] 董年虎, 王广月. 渤海湾黄河入海口区余流特性分析[J]. 黄渤海海洋, 1997(1):64-69 DONG Nianhu, WANG Guangyue. Residual current analysis of the Yellow River mouth area in Bohai gulf[J]. Journal of Oceanography of Huanghai & Bohai Seas, 1997(1):64-69.]
[53] 王楠. 现代黄河口沉积动力过程与地形演化[D]. 中国海洋大学博士学位论文, 2014 WANG Nan. Sedimentary dynamics process and topographic evolution in the modern Yellow River Mouth[D]. Doctor dissertation of Ocean University of China, 2014.]
[54] Wu G X, Wang K M, Liang B C, et al. Modeling the Morphological Responses of the Yellow River Delta to the Water-Sediment Regulation Scheme: The Role of Impulsive River Floods and Density-Driven Flows. Water Resources Research, 2023, 59(7): e2022WR033003.
[55] 寿玮玮, 宗海波, 丁平兴. 夏季黄河入海径流对黄河口及附近海域环流影响的数值研究[J]. 海洋学报, 2016, 38(7):1-13 doi: 10.3969/j.issn.0253-4193.2016.07.001 SHOU Weiwei, ZONG Haibo, DING Xingping. Numerical study of the circulation influenced by runoff input in the Huanghe (Yellow) River estuary and adjacent waters in summer[J]. Haiyang Xuebao, 2016, 38(7):1-13.] doi: 10.3969/j.issn.0253-4193.2016.07.001
[56] LeBlond P H, Emery W J, Nicol T. A climatic model of runoff-driven coastal circulation[J]. Estuarine, Coastal and Shelf Science, 1986, 23(1):59-79. doi: 10.1016/0272-7714(86)90085-5
[57] Hill D F, Ciavola S J, Etherington L, et al. Estimation of freshwater runoff into Glacier Bay, Alaska and incorporation into a tidal circulation model[J]. Estuarine, Coastal and Shelf Science, 2009, 82(1):95-107. doi: 10.1016/j.ecss.2008.12.019
[58] Chao Shenn-Yu. River-forced estuarine plume[J]. Journal of Physical Oceanography, 1987, 18(1):72-88.
[59] Cléa D, Budgell W P, Toumi R. The Congo River plume: impact of the forcing on the far-field dynamics[J]. Journal of Geographic Research: Oceans, 2013, 118(C2):964-989.
[60] 梁书秀, 孙昭晨, Nakatsuji Keiji, 等. 渤海典型余环流及其影响因素研究[J]. 大连理工大学学报, 2006(1):103-110 doi: 10.3321/j.issn:1000-8608.2006.01.021 LIANG Shuxiu, SUN Zhaochen, NAKATSUJI Keiji, et al. Research on typical residual circulation and its driving factors in the Bohai Sea[J]. Journal of Dalian University of Technology, 2006(1):103-110.] doi: 10.3321/j.issn:1000-8608.2006.01.021
[61] 王悦, 林霄沛. 地形变化下渤海湾M2分潮潮致余流的相应变化及其对污染物输运的影响[J]. 中国海洋大学学报: 自然科学版, 2006, 36(1):1-6 WANG Yue, LIN Xiaopei. The variation of M2 constituent corresponding to the change of topography in Bohai Bay and its effects on the transport of pollutants[J]. Journal of Ocean University of China, 2006, 36(1):1-6.]
[62] Wang H J, Yang Z S, Li Y H, et al. Dispersal pattern of suspended sediment in the shear frontal zone off the Huanghe (Yellow River) mouth[J]. Continental Shelf Research, 2007, 27:854-871. doi: 10.1016/j.csr.2006.12.002
[63] 王厚杰, 杨作升, 毕乃双. 黄河口泥沙输运三维数值模拟 Ⅰ——黄河口切变锋[J]. 泥沙研究, 2006, 2:1-9 doi: 10.3321/j.issn:0468-155X.2006.03.001 WANG Houjie, YANG Zuosheng, BI Naishuang. 3-D simulation of the suspended sediment transport in the Yellow River mouth Ⅰ: Shear front off the Yellow River mouth[J]. Journal of Sediment Research, 2006, 2:1-9.] doi: 10.3321/j.issn:0468-155X.2006.03.001
[64] Jiang C, Pan S Q, Chen S L. Recent morphological changes of the Yellow River (Huanghe) submerged delta: causes and environmental implications. Geomorphology, 2017, 293: 93-107.
[65] Jiang C, Chen S L, Pan S Q, et al. Geomorphic evolution of the Yellow River Delta: quantification of basin-scale natural and anthropogenic impacts. Catena, 2018, 163: 361-377.
[66] Wu X, Bi N S, Xu J P, et al. Stepwise morphological evolution of the active Yellow River (Huanghe) delta lobe (1976-2013): Dominant roles of riverine discharge and sediment grain size[J]. Geomorphology, 2017, 292:115-127. doi: 10.1016/j.geomorph.2017.04.042
[67] Xing G P, Wang H J, Yang Z S, et al. Spatial and temporal variation in erosion and accumulation of the subaqueous Yellow River Delta (1976-2004)[J]. Journal of coastal research, 2016, 74(10074):32-47.
[68] Bi N S, Wang H J, Wu X, et al. Phase change in evolution of the modern Huanghe (Yellow River) Delta: Process, pattern, and mechanisms[J]. Marine Geology, 2021, 437:106516. doi: 10.1016/j.margeo.2021.106516
[69] Pritchard D W. Estuarine Hydrography[J]. Advances in Geophysics, 1952, 1:243-280.
[70] Geyer W R, MacCready P. The Estuarine Circulation[J]. Annual Review of Fluid Mechanics, 2014, 46(46):175-197.
[71] 谢荣耀, 刘锋, 罗向欣, 等. 河控型河口盐度层化对悬沙的捕集机制——以洪季磨刀门河口为例[J]. 海洋学报, 2021, 43(5):38-49 XIE Rongyao, LIU Feng, LUO Xiangxin, et al. Sediment trapping mechanism by salinity stratification in a river-dominated estuary: A case study of the Modaomen Estuary in flood season[J]. Haiyang Xuebao, 2021, 43(5):38-49.]
[72] 龚雪雷, 姬泓宇, 李鹏, 等. 黄河三角洲近岸潮汐动力对地貌演变的响应及其沉积效应[J]. 海洋学报, 2024, 46(2): 64-78 Response of tidal dynamics to geomorphic evolution and depositional effects in the Huanghe River Delta[J]. Haiyang Xuebao, 2024, 46(2): 64-78.]
计量
- 文章访问数: 186
- HTML全文浏览量: 45
- PDF下载量: 126
