海洋科学  2023, Vol. 47 Issue (1): 13-24   PDF    
http://dx.doi.org/10.11759/hykx20220527001

文章信息

林纪江, 王平, 余世维, 栾虹, 范哲南, 牛建伟, 郑兆勇. 2023.
LIN Ji-jiang, WANG Ping, YU Shi-wei, LUAN Hong, FAN Zhe-nan, NIU Jian-wei, ZHENG Zhao-yong. 2023.
闽粤交界处大埕湾海域表层沉积物分布运移特征及沉积动力环境分区
Distribution and migration characteristics of surficial settlements and division of dynamic sedimentary environment in Dacheng Bay between Fujian and Guangdong
海洋科学, 47(1): 13-24
Marina Sciences, 47(1): 13-24.
http://dx.doi.org/10.11759/hykx20220527001

文章历史

收稿日期:2022-05-27
修回日期:2022-09-10
闽粤交界处大埕湾海域表层沉积物分布运移特征及沉积动力环境分区
林纪江1,2, 王平1,2, 余世维3, 栾虹4, 范哲南1,2, 牛建伟5, 郑兆勇2,6     
1. 国家海洋局南海信息中心, 广东 广州 510310;
2. 自然资源部海洋环境探测技术与应用重点实验室, 广东 广州 510310;
3. 广东省海洋发展规划研究中心, 广东 广州 510220;
4. 国家海洋局南海标准计量中心, 广东 广州 510310;
5. 中国科学院南海海洋研究所, 广州 510301;
6. 国家海洋局南海规划与环境研究院, 广东 广州 510310
摘要:大埕湾位于闽粤交界处, 是典型的岬角弧形开放型海湾, 本文基于实测的海流、表层沉积物数据, 分析了大埕湾海域表层沉积物的分布特征, 计算了再悬浮通量与输运通量, 结合粒径趋势分析方法探讨了表层沉积物运移趋势, 利用Flemming三角图示法对沉积动力环境进行分区并分析了泥沙输运的主要控制因素。结果表明, 研究区表层沉积物有砂质粉砂、粉砂、含砾泥和泥4种类型, 沉积物颗粒较细, 整体以粉砂和泥为主, 砂质粉砂和含砾泥呈圆斑状零星分布; 粒级组分含量和粒度参数平面分布分区明显, 黏土组分、分选系数和偏态值呈现“东北小西南大”的特点, 粉砂组分和峰态值则表现出相反趋势。观测期间表层沉积物日再悬浮量为41.86 kg/m2, 日净输运通量40.36 kg/m, 逐时再悬浮通量值和单宽输沙率大小与潮流流速变化相关性较好, 且涨潮段大于落潮段, 净输沙指向涨潮方向, 由西南向东北方向输沙。研究区沉积动力环境可划分为DⅡ、EⅡ和EⅢ 3个分区, 沉积物搬运以悬移为主, 潮流和波浪是研究区泥沙起动及输运的主要动力, 泥沙来源主要有近岸侵蚀供沙、东西两侧海域输沙以及少量短源季节性陆源碎屑补给3方面。
关键词分布特征    输运通量    运移趋势    沉积动力分区    大埕湾    
Distribution and migration characteristics of surficial settlements and division of dynamic sedimentary environment in Dacheng Bay between Fujian and Guangdong
LIN Ji-jiang1,2, WANG Ping1,2, YU Shi-wei3, LUAN Hong4, FAN Zhe-nan1,2, NIU Jian-wei5, ZHENG Zhao-yong2,6     
1. South China Sea Information Center, SOA, Guangzhou 510310, China;
2. Key Laboratory of Marine Environmental Detection Technology and Application, MNR, Guangzhou 510310, China;
3. Research Centre for the Marine Development and Planning of Guangdong Province, Guangzhou 510220, China;
4. South China Sea Standard Metrology Center, SOA, Guangzhou 510310, China;
5. South China Sea Institute of Oceanology, Chinese Academy of Sciences, Guangzhou 510301, China;
6. South China Sea Institute of Planning and Environmental Research, SOA, Guangzhou 510310, China
Abstract: Dacheng Bay is a typical cape arc open bay located between Fujian and Guangdong. Based on the measured current and surficial settlement data, this study analyzes the distribution characteristics of the surficial settlement in Dacheng Bay. Moreover, it discusses the transport trend of the surficial settlement by calculating the resuspension and transport flux and combined with the particle size trend analysis method. Furthermore, it partitions the dynamic sedimentary environment using Flemming's triangle graphical method and analyzes the main control factors of sediment transport. The results show four types of surficial settlements in the study area: sandy silt, silt, gravelly mud, and mud. The sediment particles are fine and mainly consist of silt and mud; moreover, sandy silt and gravelly mud are distributed sporadically in round spots. The plane distribution zoning of the component content and parameters of the particle size is obvious. The clay component, sorting coefficient, and skewness values are small and large in the northeast and southwest, respectively, whereas the silt component and peak value show the opposite trend. During the observation period, the daily resuspension amount of surficial settlement was 41.86 kg·m−2 and the daily net transport flux was 40.36 kg·m−1. The hourly resuspension flux and sediment transport rate per unit width correlate well with the variation in tidal current velocity, and the rising tide section exceeds the falling tide section. The net sediment transport points to the rising tide direction and transports sediment from the southwest to the northeast. The dynamic sedimentary environment of the study area can be divided into three zones: DⅡ, EⅡ, and EⅢ. The sediment transport is mainly suspended. Tidal currents and waves are the main driving forces for sediment initiation and transport in the area under study. The sediment sources mainly include nearshore erosion and sediment supply, sediment transport from the East and West seas, and a small amount of short source seasonal land-based debris supply.
Key words: distribution characteristics    transport flux    migration trend    sedimentary dynamic division    Dacheng Bay    

表层沉积物既是水体中营养物质的载体, 又是污染物循环扩散的搬运工, 其分布运移特征对海洋沉积、生态环境、岸线冲淤以及物化循环等具有重要指示意义[1], 在海岸带环境研究中具有重要意义。水动力条件、搬运介质及水深地形等不同的沉积环境下, 沉积物及其粒度参数的分布格局也不尽相同[2], 受沉积物搬运“源—汇”过程的制约, 其粒度参数是分析研究沉积物与沉积动力环境关系的重要纽带[3], 可用以反演水动力条件及沉积物物质来源[4-5], 判别沉积物输运方式与趋势[6-7], 指示沉积动力环境和海平面变化过程[8], 分析区域沉积动力特征[9-11]

大埕湾位于闽粤交界处, 是一个开放的弧形海湾, 东侧是诏安湾, 西侧是拓林湾, 西南侧是南澳岛, 岛-峡道-海湾地形相互作用对研究区的海洋动力过程产生较大影响, 也影响到沉积物的再悬浮、输运及沉积。前人对大埕湾近岸海域的研究主要集中在岸滩稳定及防护对策[12-13]、表层沉积物中重金属分布及地化环境[14-15]、台风波浪模拟分析[16]等方面, 受调查资料所限, 对该海域的沉积动力过程与环境、沉积物输运等方面的研究较少。本文基于实测的海流、表层沉积物数据, 分析了大埕湾海域表层沉积物的底质类型、粒级组分及粒度参数分布特征, 计算了表层沉积物再悬浮通量与输运通量, 结合粒径趋势分析方法探讨了表层沉积物运移趋势, 对沉积动力环境进行分区并分析了泥沙输运的主要控制因素。研究成果不仅有助于深入了解大埕湾陆、海相互作用的沉积环境过程和海岸演变趋势, 丰富对开放性海湾沉积动力特征的认识, 还可为大埕湾砂质岸线保护、岸滩整治修复、海岸工程合理开发建设提供参考。

1 研究区概况与资料方法 1.1 研究区概况

大埕湾位于闽粤交界处, 呈半月状的岬角弧形开放型砂质海湾, 湾口朝向东南, 东起宫口港, 西至鸡草角。在地貌上自成一个较完整的沙坝泻湖海岸, 分布有宫口港(泻湖)、渡西沙嘴(拦湾沙坝)、潮汐通道及落潮沙洲等次级地貌单元, 海滩物质主要为细砂[12]。有季节性小溪流—港仔从海湾出海, 局部地貌受径流注入影响, 底质以细颗粒的粉砂和黏土质粉砂为主[11, 14], 水深在20 m以浅(图 1)。波浪以风浪为主, 波向和频率季节变化明显, 6—8月常浪向为SW向和SE向, 其余月份常浪向为E向和NE向, 强浪向为NE向、ENE向和E向, 夏季受台风影响波浪强但持续时间短, 冬季受东北季风影响波浪弱于台风浪但持续时间长, 平均波高冬季大于夏季[12, 13, 16], 潮汐属不正规半日潮, 平均潮差1.64 m, 受地形影响, 潮流主要为呈NE-SW(S)的往复流, 落潮流向SW(S), 涨潮流向NE, 落潮历时小于涨潮历时, 余流主要与涨潮方向一致, 在较强东北季风影响下, 也可出现西南向余流[12, 13]。大埕湾降水充沛, 据南澳海洋站1986年—2015年统计, 雨季(4—9月)平均年降水量共1 025.7 mm, 占全年降水量的81%, 受季风和热带气旋影响, 个别月份累年月平均降水量在155 mm以上。

图 1 研究区调查站位分布图 Fig. 1 Distribution map of survey stations in the study area
1.2 资料来源与处理

交通运输部天津水运工程科学研究所于2015年9月在大埕湾近岸海域非等间距布设49个底质取样站位(图 1), 最大采样间距1.5 km, 最小采样间距0.5 km, 用蚌式采泥法采样(曙光采泥器), 现场密封保存带回实验室, 取0~2 cm的表层样品, 用Wi­n­ner2008A激光粒度仪测试粒度, 测试范围: 0.05~2 000 μm。

交通运输部天津水运工程科学研究所于2015年9月在大埕湾近岸海域布设1个海流站位, 采用SLC 9-2直读式海流计测量流速流向, 流速测量范围: 3~ 350 cm/s, 准确度: ≤±1.5%。流向测量范围: 0~360°, 准确度: ≤±4°。海流观测采用六点法(表层、0.2H、0.4H、0.6H、0.8H、底层)连续整点观测26 h。观测期间海况良好, 风和浪均不大。

1.3 研究方法与公式

本文利用FOLK等[17]的公式方法来计算粒度参数, 采用FOLK[18]分类命名方法对研究区沉积物进行定名。粒度参数分级的定性描述参考贾建军等[19]2002年提出的分级定性描述(表 1)。采用GAO等[20]粒径趋势分析沉积物净输运趋势, 该方法根据不同类型的粒径出现概率的差异化趋势, 可以较好地反映沉积物的净输运趋势; 利用FIE­MMING[21]三角图示法对研究区沉积动力环境进行分析, 该方法根据表层沉积物中砂、粉砂、黏土的百分比进行划分区域, 可指示沉积动力环境的强弱。

表 1 粒度参数定性描述[19] Tab. 1 Qualitative description of grain size parameters
分选 偏态 峰态
分选程度 分选系数 偏态 偏态值 峰态 峰态值
分选极好 < 0.35 极负偏 < –1.50 非常窄 < 0.72
分选好 0.35~0.5 负偏 –1.50~–0.33 很窄 0.72~1.03
分选较好 0.5~0.71 近对称 –0.33~+0.33 中等 1.03~1.42
分选中等 0.71~1.00 正偏 +0.33~+1.50 1.42~2.75
分选较差 1.00~2.00 极正偏 > +1.50 很宽 2.75~4.50
分选差 2.00~4.00 非常宽 > 4.50
分选极差 > 4.00

再悬浮通量采用PARTHENIADES[22]提出的公式:

$\frac{\mathrm{d} M_E}{\mathrm{~d} t}=\mathrm{E}\left(\frac{\tau_0}{\tau_{c r}}-1\right), $ (1)

式(1)中, E是再悬浮常数, 本文取1.25×10–4进行计算; $\tau_0$为底切应力; $\tau_{c r} $为临界切应力。

表层沉积物的单宽输沙率能较好地反映其输运速率, 本文参考BARTHOLDY等[23]提出的半经验公式来计算单宽输沙率以及涨落潮段输沙通量。具体公式见式(2)~(6):

$q_{s b}=0.5 u_f^{\partial} V_D^2 d^{-1}, $ (2)
$u_f=\frac{u}{\left(8.5+2.5 \ln \left(\frac{z}{k}\right)\right)}, $ (3)
$k=\frac{d}{12}, $ (4)
$ V_D=\left(6+2.5 \ln \left(\frac{D}{k}\right)\right) u_f, $ (5)
$Q_{s b}=q_{s b} \times T, $ (6)

式中, qsb表示单宽输沙率(kg/(m/s)), VD表示垂线平均流速, uf表示摩阻流速(m/s), d表示中值粒径, 0.5表示经验系数(单位: kg·s4·m–5), u表示实测流速, k表示底床粗糙度(m), z表示潮流测量高度(m), D表示水深(m), Qsb表示平均单宽通量(kg/m), T表示时间。

2 海流与底质特征 2.1 海流特征

C1站位实测资料显示, 流速最大值为52 cm/s, 出现在表层, 最小流速值为8 cm/s, 出现在表层和0.2H层, 各层平均流速值自表层至底层先增大后减小(表 2, 图 2), 垂向平均流速值介于14~47 cm/s, 底层流速值介于12~40 cm/s, 落潮流向SW, 涨潮流向NE, 涨潮历时大于落潮历时。

表 2 各层流速流向极值统计 Tab. 2 Statistics of the extreme value of velocity and flow direction of each layer
站位
C1
表层 0.2 H 0.4 H 0.6 H 0.8 H 底层 垂向平均
流速/(cm/s) 流向/(°) 流速/(cm/s) 流向/(°) 流速/(cm/s) 流向/(°) 流速/(cm/s) 流向/(°) 流速/(cm/s) 流向/(°) 流速/(cm/s) 流向/(°) 流速/(cm/s) 流向/(°)
最大值 52 304 48 334 50 330 46 342 44 356 40 320 47 331
最小值 8 28 8 30 12 26 14 14 12 20 12 8 14 21
平均值 23 158 24 156 27 147 28 139 27 140 24 131 27 145

图 2 实测底层流速流向过程曲线 Fig. 2 Process curve of measured bottom flow velocity and direction
2.2 表层沉积物特征

表层沉积物FOLK分类命名投点结果(图 3)显示, 研究区表层沉积物有砂质粉砂、粉砂、含砾泥和泥4种类型, 沉积物颗粒较细(图 4), 整体以粉砂和泥为主(图 5表 3), 分别占样品总数的55.1%、34.7%, 较粗颗粒的砂质粉砂和含砾泥在研究区海域内呈圆斑状零星分布, 二者总共占比10.2%。各类沉积物平面分布特征如下:

图 3 研究区表层沉积物的FOLK分类命名分布 Fig. 3 Sediment type of surficial settlement according to FOLK's classification in the study area

图 4 各类型沉积物粒级频率曲线 Fig. 4 Grain size frequency curve of various sediments 注: N为样品数, 不同颜色曲线表示不同站位

图 5 研究区表层沉积物类型分布图 Fig. 5 Distribution of surficial settlement types in the area under study

表 3 表层沉积物类型及粒度参数 Tab. 3 Statistics of the surficial settlement type and grain size parameters
沉积物类型 样品数 取值 砾/% 砂/% 粉砂/% 黏土/% Mz/φ 偏态 峰态 分选系数
粉砂 27 最大值 0 5 76.12 30.35 7.27 0.64 1.19 2.27
最小值 0 0 65.8 18.92 6.51 –0.15 0.73 1.29
平均值 0 1.24 72.13 23.33 7.08 0.03 1.02 1.60
17 最大值 0 4.23 65.51 42.46 7.58 0.37 0.75 2.33
最小值 0 0.94 53.93 32.45 6.97 0.19 0.65 2.1
平均值 0 2.24 60.37 37.39 7.36 0.27 0.69 2.21
砂质粉砂 2 最大值 0 11.73 68.15 19.6 6.72 0.11 1.21 2.03
最小值 0 8.9 64.73 19.17 6.56 0.1 1.05 1.9
平均值 0 10.32 66.44 19.39 6.64 0.11 1.13 1.97
含砾泥 3 最大值 0.03 3.46 71.84 37.05 7.3 0.53 1.07 2.29
最小值 0.03 0 59.46 23.69 6.94 –0.14 0.7 1.35
平均值 0.03 1.54 67.50 29.44 7.11 0.24 0.86 1.94

粉砂: 在研究区的东北侧广泛分布并零星分布于西南侧海域, 砂含量相对较低仅占0~5%, 粉砂含量占65.8%~76.12%, 黏土含量占18.92%~30.35%, 偏态在–0.15~0.64, 属于近对称至正偏, 峰态值在0.73~ 1.19之间, 属于很窄至中等, 分选在1.29~2.27, 属于分选较差至分选差, 平均粒径在6.51~7.27φ。

泥: 广泛分布于研究区域西南侧海域, 砂含量相对较低仅占0.94%~4.23%, 粉砂含量占53.93%~ 65.51%, 黏土含量占32.45%~42.46%, 偏态在0.19~ 0.37, 属于近对称至正偏, 峰态值在0.65~0.75, 属于非常窄至很窄, 分选在2.1~2.33, 属于分选差, 平均粒径在6.97~ 7.58φ。

砂质粉砂: 砂质粉砂小面积分布于研究区中部近岸, 砂含量占8.9%~11.73%, 粉砂含量占64.73%~68.15%, 黏土含量占19.17%~19.6%, 偏态在0.1~0.11, 属于近对称, 峰态值在1.05~1.21, 属于中等, 分选在1.9~2.03, 属于分选较差至分选差, 平均粒径在6.56~6.72φ。

含砾泥: 呈东西向零星分布于研究区域, 砾含量为0.03%, 砂含量相对较低仅占0~3.46%, 粉砂含量占59.46%~71.84%, 黏土含量占23.69%~37.05%, 偏态在–0.14~0.53, 属于近对称至正偏, 峰态值在0.7~1.07, 属于非常窄至中等, 分选在1.35~2.29, 属于分选较差至分选差, 平均粒径在6.94~7.3φ。

粒级组分含量和粒度参数平面分布分区明显(图 6), 黏土组分、分选系数和偏态值呈现“东北小西南大”的特点, 粉砂组分和峰态值则表现出相反趋势。平均粒径φ值表现为东北侧海域分布均匀, 西南侧海域差值较大。整体上东北侧海域泥沙颗粒略粗于西南侧海域, 分选好于西南侧海域。

图 6 表层沉积物组分含量与粒度参数分布图 Fig. 6 Distribution of component content and particle size parameters in the surficial settlement
3 沉积物再悬浮与输运特征 3.1 再悬浮特征

根据前述公式计算了潮周期内C1站位的再悬浮通量(图 7), 结果显示, 潮周期内C1站位再悬浮通量值在4.50×10–5~1.26×10–3 kg/(m2/s), 日再悬浮量为41.86 kg/m2, 涨潮段再悬浮通量大于落潮段, 再悬浮通量变化与潮流流速变化相关性较好。研究区海底表层沉积物主要为粉砂和泥, 颗粒较细, 在研究区较强水动力作用下, 泥沙较易起动, 流速增大, 再悬浮通量随之增大, 反之, 再悬浮通量逐渐减小; 当流速降至最低时, 因底部切应力小于临界值, 悬浮泥沙发生沉降。

图 7 潮周期内观测站位再悬浮通量与流速的变化 Fig. 7 Variation between resuspension and velocity in the tidal cycle at the observed station
3.2 表层沉积物输运通量

根据公式(2)~(6)计算得到单日内研究区涨、落潮段平均单宽通量和平均单宽净通量, 结果见表 4图 8, 表中涨、落潮单宽输沙率为涨、落潮段各自的平均值, 本文根据研究区潮位曲线, 先分出涨、落潮段, 计算出每个小时的单宽输沙率, 然后求其平均, 作为涨、落潮段的单宽输沙率。

表 4 表层沉积物输运通量 Tab. 4 Transport fluxes in surficial settlements
站位 涨潮单宽输沙率(kg/(m/s)) 落潮单宽输沙率(g/(m/s)) 涨潮输沙(kg/(m/d)) 落潮输沙(kg/(m/d)) 净输沙(kg/(m/d))
1 1.28 0.67 64.59 24.23 –40.36
注: “–”表示净输沙方向指向涨潮方向

图 8 表层沉积物单日逐时单宽输沙率 Fig. 8 Sediment transport rate per hour in a day of surficial settlement

结果显示, 研究区表层沉积物单日逐时单宽输沙率在0.024~4.648 g/(m/s), 潮流流速变化对逐时单宽输沙率的影响显著, 单宽输沙率最小值出现在3时, 对应的实测流速值较低, 仅14 cm/s, 最大值出现在14时, 对应的实测流速值较高, 为40 cm/s, 水深变化对逐时单宽输沙率的影响不明显。涨潮段单宽输沙率1.28 g/(m/s), 落潮段单宽输沙率0.67 g/(m/s), 日净输沙40.36 kg/m, 净输沙指向涨潮方向, 总体呈NE向, 即由西南向东北方向输沙, 与前人学者[13]根据水下地形的变化特征分析大埕湾岸段泥沙运动是由西南向东北处运动的结论相一致。

3.3 表层沉积物运移趋势

Gao-Collins粒径趋势分析方法分析表层沉积物净运移趋势的关键是特征距离的选择, 前有学者[19, 24]认为最佳特征距离为最大采样距离或$\sqrt{2} $倍的最大采样距离, 本文选取了1.5和2 km的特征距离, 分别绘制了表层沉积物运移趋势图(图 9), 图中矢量的方向代表了该采样点沉积物的净搬运方向, 而该矢量的长度并不代表沉积物搬运率的大小, 仅表示该粒径趋势的显著性。结果显示, 1.5和2 km的特征距离下, 大埕湾表层沉积物运移趋势基本一致, 但2 km的特征距离下趋势规律更显著, 由西南向东北运移, 指示研究区沉积物净输运方向单一, 指向东北, 与前文通过计算沉积物净输运通量及阎新兴等[12]根据水下地形的变化特征分析的近岸泥沙运移趋势结论一致。

图 9 大埕湾表层沉积物运移趋势 Fig. 9 Grain size net trend of the surficial settlements in Dacheng Bay
4 沉积动力环境分区与泥沙输运控制因素 4.1 沉积动力环境分区

Flemming三角图解法以砂、粉砂、黏土为三端元, 根据三者两两比值将三角图分为25个分区。本文49个样品中, 只有3个样品含有砾组分, 且含量仅为0.03%, 可忽略砾组分的影响, 结果(图 10)显示研究区沉积物主要分布于EⅡ和EⅢ区, DⅡ区有少量分布, 就各沉积物类型来看, 砂质粉砂分布于DⅡ区, 泥分布于EⅢ区, 含砾泥和粉砂在EⅡ和EⅢ区均有分布且粉砂主要分布于EⅡ区。样品主要集中于E区, 说明研究区沉积物粒径较集中, 以粉砂和黏土为主, 样品主要分布于Ⅱ和Ⅲ区, 粉砂黏土比值较高, 指示研究海域递变悬移组分较多, 水体扰动环境较大, 水动力环境较强。

图 10 研究区表层沉积物Flemming三角分区图 Fig. 10 Flemming triangle division map of surficial settlements in the study area

砂泥比值、粉砂黏土比值可较好地反映沉积动力条件强弱和水体扰动强弱环境, 其比值等值线分布见图 11。由图 11可见研究区内砂泥比值较低, 均低于0.14, 大部分区域低于0.04, 表明研究区内沉积物颗粒整体偏细,高于0.04的区域仅在研究区中部、南部部分海域出现, 指示该处沉积物粗颗粒物质较多, 东北侧边缘处砂泥比值为0, 表明此处海域沉积物颗粒最细, 水动力条件相对较弱。砂泥比值呈现“西南大东北小”的特征, 指示水动力由西南向东北减弱。粉砂黏土比值在1.27~3.85, 分区明显, 东北侧海域大于2.5, 颗粒整体较粗, 西南侧海域小于2.5, 指示研究区东北侧海域水体中递变悬移组分较多, 水体扰动性强于西南侧海域。

图 11 研究区沉积物砂/泥比值、粉砂/黏土比值分布 Fig. 11 Distribution map of sand/mud and silt/clay in the study area

根据Flemming三角图分区结果, 结合FOLK沉积物分类及研究区砂泥比值、粉砂黏土比值分布情况, 对研究区进行沉积环境分区, 分区结果见图 12

图 12 研究区沉积动力环境分区图 Fig. 12 Dynamic sediment environment divisions in the study area

DⅡ区: 主要分布于研究海域的中部, 靠近季节性小溪流-港仔入海口的下方, 底质以砂质粉砂为主, 颗粒平均粒径较粗, 水动力强, 水体扰动性一般, 泥沙运动受潮流和波浪的控制, 也会受季节溪流的影响, 沉积物搬运以跃移和悬移为主。

EⅡ区: 主要集中于研究区东北侧海域, 底质对应为粉砂和含砾泥, 颗粒平均粒径较细, 砂泥比值低, 粉砂黏土比值高, 递变悬浮组分含量较高, 水动力较强, 水体扰动环境强, 泥沙运动受潮流和波浪的控制, 沉积物搬运以悬移为主。

EⅢ区: 主要集中于研究区西南侧海域, 底质对应为泥、含砾泥和粉砂, 颗粒平均粒径最细, 砂泥比值较低, 粉砂黏土比值较高, 递变悬浮组分含量较高, 水动力较强, 水体扰动环境较强, 泥沙运动受潮流和波浪的控制, 沉积物搬运以悬移为主。

4.2 泥沙输运控制因素分析

通常来讲, 物质来源、水动力条件和海底地形等是决定海洋沉积物分布格局及输运特征的主要因素, 物质来源决定了沉积物的最初样貌, 波浪、海流、风暴潮等水动力在偶发或周期性作用下扰动水体, 改变并重新塑造沉积物的样貌, 同时控制沉积物的时空分布, 海底地形则制约沉积物的输运过程。

在相近的水动力和海底地形条件下, 沉积物类型(图 5)与粉砂、黏土粒级组分和分选系数、偏态、峰态等粒径参数(图 6)分布呈现的“东北-西南”明显分区格局, 应该是泥沙来源不同所致, 大埕湾海域无大江大河入海, 只有短源季节溪流或冲沟入海, 但研究区降水充沛且集中, 根据前文较粗颗粒的砂质组分高含量集中在港仔入海口下方的平面分布特征, 初步认为在洪水期溪流或冲沟可能会携带部分陆源碎屑入海。大埕湾西侧有黄岗河流入拓林湾, 更西侧的韩江入海携沙量较为客观, 东侧有黄岗河支流流入宫口港以及存在较大的诏安湾, 在沿岸潮流的作用下, 东西两侧细颗粒的悬浮泥沙将输运至大埕湾, 并在潮流流速减小时发生沉积。亦即研究区西侧海域有拓林湾(黄岗河)和韩江的泥沙补给, 研究区东侧海域有宫口港(黄岗河东支)和诏安湾泥沙输入, 但由于近年频繁修建围堤工程, 使得湾内纳潮量减少, 潮流携带输运泥沙能力减弱, 两侧输沙量也有下降。阎新兴等[13]和袁国明等[14]分别从水下地形变化和沉积物主控因子等角度也得出相近结论。蔡锋等[12]从岸滩地貌演变方面认为大埕湾中部和西部岸滩呈侵蚀状态, 尽管修筑了人工防护设施, 但由于直接受偏E向优势浪的影响, 加之沿岸裂流、海岸活动、人工采沙等因素, 近岸侵蚀供沙现象依然存在。

潮流是大埕湾海域泥沙起动及输运的主要动力, 研究区海域表层沉积物在单纯潮流作用下即可起动, 再悬浮通量与潮流具有较好的对应关系, 且日再悬浮量较为可观。波浪对研究区沉积物分布与输运也有较大影响, 波浪在近岸掀动底质泥沙进入水体, 在潮流作用下往复搬运, 彭程等[16]基于SWAN模型模拟的大埕湾波浪分布显示, ESE~SE向波浪可直接入射湾内, 由于波浪较强, 直接控制湾内波浪方向, SSE~SW向波浪需绕射进入湾内和入射波浪共同控制湾内波浪方向。尤其研究区位于闽粤交界处, 属于热带季风气候区, 受台风和东北季风的双重影响, 且具有年季重复性, 波浪对泥沙分布和运移的影响不可忽视。

综上, 大埕湾泥沙来源主要有近岸侵蚀供沙、东西两侧海域输沙以及少量短源季节性陆源碎屑补给3方面, 潮流和波浪是大埕湾海域泥沙起动及输运的主要动力。

5 结论

研究区表层沉积物有砂质粉砂、粉砂、含砾泥和泥4种类型, 沉积物颗粒较细, 整体以粉砂和泥为主, 砂质粉砂和含砾泥呈圆斑状零星分布; 粒级组分含量和粒度参数平面分布分区明显, 黏土组分、分选系数和偏态值呈现“东北小西南大”的特点, 粉砂组分和峰态值则表现出相反趋势。

观测期间表层沉积物日再悬浮量为41.86 kg/m2, 日净输运通量40.36 kg/m, 逐时再悬浮通量值和单宽输沙率大小与潮流流速变化相关性较好, 且涨潮段大于落潮段。

研究区水体扰动环境较强, 潮流和波浪是研究区泥沙起动及输运的主要动力, 研究区表层沉积物净输运方向单一, 指向涨潮方向, 由西南向东北运移。

根据Flemming图示法和沉积物粒度分布情况, 研究区沉积动力环境可划分为DⅡ、EⅡ和EⅢ 3个分区, 沉积物搬运以悬移为主; 泥沙来源主要有近岸侵蚀供沙、东西两侧海域输沙以及少量短源季节性陆源碎屑补给3方面。

参考文献
[1]
林纪江, 胡日军, 王平, 等. 庙岛海峡周边海域表层沉积物再悬浮及悬浮泥沙输运机制[J]. 海洋地质与第四纪地质, 2022, 42(3): 9-24.
LIN Jijiang, HU Rijun, WANG Ping, et al. Surface sediment resuspension and suspended sediment transportation mechanism in the waters around Miaodao Strait[J]. Marine Geology & Quaternary Geology, 2022, 42(3): 9-24. DOI:10.16562/j.cnki.0256-1492.2021122002
[2]
瞿洪宝, 苟鹏飞, 孙龙飞, 等. 海南岛崖州湾表层沉积物空间分布特征及其受控机制[J]. 海洋学报, 2021, 43(12): 70-81.
QU Hongbao, GOU Pengfei, SUN Longfei, et al. Spatial distribution characteristics and control mechanism of surficial settlements in Yazhou Bay, Hainan Island[J]. Acta Oceanologica Sinica, 2021, 43(12): 70-81. DOI:10.12284/j.issn.0253-4193.2021.12.hyxb202112007
[3]
杨阳, 高抒, 周亮, 等. 海南新村港潟湖表层沉积物粒度特征及其沉积环境[J]. 海洋学报, 2016, 38(1): 94-105.
YANG Yang, GAO Shu, ZHOU Liang, et al. Grain size characteristics and sedimentary environment of surficial settlements in Xincungang lagoon, Hainan[J]. Acta Oceanologica Sinica, 2016, 38(1): 94-105. DOI:10.3969/j.issn.0253-4193.2016.01.009
[4]
袁萍. 渤海表层沉积物的空间分布及其与物源和沉积动力环境的关系[D]. 青岛: 中国海洋大学, 2015.
YUAN Ping. Spatial distribution of surficial settleme nts in the Bohai Sea and their relationship with provenance and sedimentary dynamic environment[D]. Qing-dao: Ocean University of China, 2015.
[5]
王伟伟, 付元宾, 李树同, 等. 渤海中部表层沉积物分布特征与粒度分区[J]. 沉积学报, 2013, 31(3): 478-485.
WANG Weiwei, FU Yuanbin, LI Shutong, et al. Distribution characteristics and grain size division of surficial settlements in the central Bohai Sea[J]. Acta Sedimentology, 2013, 31(3): 478-485. DOI:10.14027/j.cnki.cjxb.2013.03.011
[6]
胡日军, 吴建政, 朱龙海, 等. 东海舟山群岛海域表层沉积物运移特性[J]. 中国海洋大学学报(自然科学版), 2009, 39(3): 495-500, 442.
HU Rijun, WU Jianzheng, ZHU Longhai, et al. Characteristics of surficial settlement transport in the Zhoushan Islands in the East China Sea[J]. Periodical of Ocean University of China, 2009, 39(3): 495-500, 442. DOI:10.16441/j.cnki.hdxb.2009.03.022
[7]
吕纪轩, 胡日军, 李毅, 等. 烟台北部近岸海域表层沉积物粒度分布及沉积动力环境特征[J]. 海洋地质前沿, 2020, 36(4): 27-36.
LV Jixuan, HU Rijun, LI Yi, et al. Grain size distribution and sedimentary dynamic environment characteristics of surficial settlements in the coastal waters of northern Yantai[J]. Marine Geology Frontiers, 2020, 36(4): 27-36.
[8]
刘成, 胡日军, 朱龙海, 等. 庙岛群岛海域沉积动力环境分区及沉积物输运趋势[J]. 海洋地质前沿, 2018, 34(8): 24-33.
LIU Cheng, HU Rijun, ZHU Longhai, et al. Sedimentary dynamic environment division and sediment transport trend in Miaodao Islands sea area[J]. Frontiers in Marine Geology, 2018, 34(8): 24-33.
[9]
杨旭辉, 冯秀丽, 褚忠信, 等. 中国东部陆架表层沉积物粒度特征及其沉积环境浅析[J]. 中国海洋大学学报(自然科学版), 2012, 42(Z2): 132-140.
YANG Xuhui, FENG Xiuli, CHU Zhongxin, et al. Analysis of grain size characteristics and sedimentary environment of surficial settlements on the continental shelf in eastern China[J]. Periodical of Ocean University of China, 2012, 42(Z2): 132-140.
[10]
肖晓, 冯秀丽, 石要红, 等. 广西沿岸海域表层沉积物粒度分布与输运特征[J]. 中国海洋大学学报(自然科学版), 2018, 48(7): 65-72.
XIAO Xiao, FENG Xiuli, SHI Yaohong, et al. Grain size distribution and transport characteristics of surficial settlements in coastal waters of Guangxi[J]. Periodical of Ocean University of China, 2018, 48(7): 65-72.
[11]
陈翰, 陈忠, 颜文, 等. 汕头近岸海域表层沉积物粒度特征及其输运趋势[J]. 沉积学报, 2014, 32(2): 314-324.
CHEN Han, CHEN Zhong, YAN Wen, et al. The grain size characteristics and transport trend of surficial settlements in the coastal area of Shantou[J]. Journal of sedimentation, 2014, 32(2): 314-324.
[12]
蔡锋, 苏贤泽, 高智勇, 等. 闽粤交界的大埕湾岸滩稳定分析及岸滩防护对策[J]. 台湾海峡, 2003, 4: 518-525.
CAI Feng, SU Xianze, GAO Zhiyong, et al. Stability analysis and beach protection countermeasures of Dacheng Bay at the junction of Fujian and Guangdong[J]. Taiwan Strait, 2003, 4: 518-525.
[13]
阎新兴, 刘国亭, 蔡嘉熙. 潮州三百门电厂煤港工程海岸动力地貌特征分析[J]. 水道港口, 2004, 25(1): 38-44.
YAN Xinxing, LIU Guoting, CAI Jiaxi. Analysis on coastal dynamic geomorphology of coal harbor project for Sanbaimen Power Plant[J]. Waterway Port, 2004, 25(1): 38-44.
[14]
袁国明, 何桂芳. 广东柘林湾至大埕湾海域表层沉积物重金属分布与环境地球化学特征[J]. 地球与环境, 2015, 43(2): 190-197.
YUAN Guoming, HE Guifang. Distribution and environmental geochemical characteristics of heavy metals in surficial settlements from Zhelin Bay to Dacheng Bay, Guangdong[J]. Earth and Environment, 2015, 43(2): 190-197.
[15]
张小华, 钟立峰, 苗莉, 等. 粤东典型海湾表层沉积物Cd的分布及其潜在生态危害评价[J]. 热带海洋学报, 2013, 32(2): 118-127.
HANG Xiaohua, ZHONG Lifeng, MIAO Li, et al. Distribution and potential ecological hazard assessment of Cd in surficial settlements of typical bays in eastern Guangdong[J]. Journal of Tropical Oceanography, 2013, 32(2): 118-127.
[16]
彭程, 张慈珩, 管宁. 粤东大埕湾海域台风浪模拟与分析[J]. 中国水运(下半月), 2017, 17(6): 111-115.
PENG Cheng, ZHANG Ciheng, GUAN Ning. Simula tion and analysis of typhoon waves in Dacheng Bay, eastern Guangdong[J]. China Water Transportation (Seco nd Half of the Month), 2017, 17(6): 111-115.
[17]
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.
[18]
FOLK R L, ANDRESWS P B, LEWIS D W. Detrital sedimentary rock classification and nomenclature for use in New Zealand[J]. New Zealand Journal of Geology and Geophysics, 1970, 13: 937-968.
[19]
贾建军, 高抒, 薛允传. 图解法与矩法沉积物粒度参数的对比[J]. 海洋与湖沼, 2002, 33(6): 577-582.
JIA Jianjun, GAO Shu, XUE Yunchuan. Comparison of sediment grain size parameters between graphical method and moment method[J]. Oceanologia et Limno logia Sinica, 2002, 33(6): 577-582.
[20]
GAO S, COLLINS M. Net sediment transport patterns inferred from grain-size trends, based upon definition of "transport vectors"[J]. Sedimentary Geology, 1992, 81(1/2): 47-60.
[21]
FLEMMING B W. A revised textural classification of gravel free muddy sediments on the basis of ternary diagrams[J]. Continental Shelf Research, 2000, 20: 1125-1137.
[22]
PARTHENIADES E. Erosion and deposition of cohesive soil[J]. Journal of the Hydraulics Division, 1965, 2(4): 150-139.
[23]
BARTHOLDY J, BARTHOLOMAE A, FLEMMING B W. Grain-size control of large compound flow-transverse bedforms in a tidal inlet of the Danish Wadden Sea[J]. Marine Geology, 2002, 188(3/4): 391-413.
[24]
高抒. 沉积物粒径趋势分析: 原理与应用条件[J]. 沉积学报, 2009, 27(5): 826-836.
GAO Shu. Trend analysis of sediment particle size: principle and application conditions[J]. Chinese Journal of Sedimentology, 2009, 27(5): 826-836.