Concentration of suspended particulate matter and magnetic minerals from surface seawater in Russian Arctic Seas: Distribution and influencing factors
-
摘要: 为了解俄罗斯极地海悬浮体颗粒物（Suspended Particulate Matter，SPM）的分布特征及其沉积学意义，对2019年中俄北极联合考察（AMK78航次）在喀拉海-拉普捷夫海-东西伯利亚海表层海水采集的SPM样品开展了颗粒浓度、组成和岩石磁学分析。研究发现，SPM主要由陆源碎屑和硅质浮游生物碎屑组成。悬浮体浓度在德米特里拉普捷夫海峡及其东侧海域最高，鄂毕河和叶尼塞河河口外侧次之，其他海域整体较低。悬浮体中携磁矿物为单畴、多畴磁铁矿。陆源碎屑集中分布在近岸和河流入海口附近海域,离海岸和河口较远海域悬浮体中硅质生物碎屑占比升高。SPM浓度主要受到河流搬运入海的陆源碎屑通量和海岸侵蚀的影响，而SPM中磁性矿物与流域内岩石类型有关，并受西伯利亚沿岸流的影响，磁性矿物集中在洋流流速缓慢的区域，粒径粗的磁性矿物分布在沿岸地区，可能与海岸侵蚀有关。Abstract: To understand the distribution characteristics and sedimentological significance of suspended particulate matter (SPM) in the Russian Arctic Seas, the particle concentration, composition, and rock magnetism of SPM samples collected in the surface seawater of the Kara Sea-Laptev Sea-Eastern Siberia Sea during the Sino-Russian joint Arctic expedition (AMK78 voyage) were analyzed. Results show that the SPM was mainly composed of terrigenous detritus and siliceous plankton debris. The SPM concentration was the highest in the Dmitry Laptev Strait and its eastern sea areas, followed by the outer sides of the Ob and Yenisei estuaries, and were generally lower in the other sea areas. The magnetic minerals in the SPM were single-domain and multi-domain magnetite. The terrigenous detritus was concentrated in the waters near the coast and the estuary of the river, and the proportion of siliceous bioclasts in the SPM in the sea farther from the coast and the estuary increases. The SPM concentration was mainly affected by the flux of terrigenous debris transported into the sea by rivers and coastal erosion, while the magnetic minerals in SPM were related to the rock types in the watershed and were affected by the Siberian Coastal Current. Coarse magnetic minerals were distributed in the coastal areas, which might be related to coastal erosion.
-
海湾具有独特的自然环境和优越的地理位置，对人类生存及社会发展具有极其重要的意义[1]，然而海湾是一个相对封闭的海洋环境，其水交换能力较差，生态环境较为脆弱，对人工构筑物建设等人类活动的干扰较为敏感[2]。人工构筑物的建设可以有效增加自然岸线、提高海洋资源利用效率，但人工构筑物的建设会改变原有海域的地形条件，直接或间接地影响海洋生态环境[3]，导致所在海域潮流、潮位、波浪等水动力条件发生变化[4-5]，甚至经常引起更大程度的极端事件以及海平面上升现象[6]。因此，为了保护海湾环境实现可持续发展，研究人工岛建设等人类活动对海湾环境的影响尤为重要[1, 7]。
围填海工程建设等人类活动对水动力环境的影响，一直以来吸引着国内外学者的广泛关注。近岸潮汐的变化与围填海工程的建设密切相关[8]，Byun等[9]认为在潮汐占主导地位、壅塞的河口或港湾环境中建造堤坝或海堤等围填海工程，可能导致潮汐振幅显著增加和水流速度下降。围填海工程的建设不会导致潮流特征大范围改变，仅在工程区域附近有较大变化[10]，流向受工程影响会发生偏转，工程防波堤等会引起挑流作用的区域流速增大，其余区域流速普遍减小[11-12]。法国兰斯河口世界第二大潮汐电站[13]以及加拿大芬迪湾潮汐动力潟湖工程[14]对附近潮流也产生了类似的影响。海岸的形态对于波浪能的耗散具有关键作用[15]，围填海工程造成的海岸形态的改变往往会中断与波浪条件相关联的沿岸泥沙输运，同时会破坏海岸抵抗风暴的能力[16]，但设计良好的防波堤会有效减小入港波浪的有效波高，从而达到保护港口的目的[17-18]。围填海造成的海域面积减小会直接导致纳潮量减少[19-20]，例如马尔代夫国际机场扩建工程[21]、雅加达湾巨型海堤[22]以及胶州湾围填海工程[23]等。水交换也会受到围填海工程的影响[24]，岸线大规模的改变会导致水交换能力减弱，通过修复岸线对水交换能力具有促进作用[25]，岸线变迁引起的水动力环境的改变是影响水交换率的主要因素[26]。此外，还有学者在水动力环境变化的基础上探讨了冲淤环境以及地貌演变对海岸工程建设的响应[27-29]，同时对于围填海工程造成的水质污染[30-31]、海洋生物受损[32-33]等方面的研究也有很多学者涉猎。然而，前人的研究多集中于河口或者较为开阔的海湾或海域中围填海工程的建设，对于小海湾中大规模离岸人工岛群建设对海洋环境影响的相关研究却鲜有报道，因此，本文针对龙口湾水动力特征及其对人工岛群建设的响应进行较为全面的研究。
龙口湾面积约84.13 km2，龙口人工岛面积约35.23 km2，约占龙口湾面积的41%，因此，在龙口湾这类小海湾中建设大规模的围填海工程，势必对其海洋生态环境造成明显影响。刘星池[34]等通过建立龙口湾水沙数值模型来预测人工岛建设对海洋环境的影响，人工岛建设会导致龙口湾潮流流速流向均发生不同程度的变化，冲淤特征表现为人工岛北侧和西南侧海域以冲淤为主，西侧海域以侵蚀为主[35]，人工岛外悬浮泥沙浓度和悬沙通量大于岛内水道[36]，同时表层沉积物组合特征也产生了重要变化[37]。但是，目前对于人工岛群建设对龙口湾水环境的影响研究主要集中在潮流场变化，缺乏其对水动力环境的系统研究；同时，龙口人工岛的实际建设方案与设计方案存在明显差异，导致前人对于该区域的研究存在误差，不能完全反映人工岛群建设对龙口湾水动力特征的影响。因此，本文利用龙口湾以及附近海域的水深地形、潮汐潮流等实测资料，利用Mike21数学模型，分析了人工岛建设前后龙口湾的水动力特征，并在此基础上探讨了人工岛群建设对龙口湾潮流、波浪、潮位、纳潮量以及水交换的影响。
1. 研究区概况
龙口湾为莱州湾的一个附属海湾，湾廓呈半圆形，是典型的连岛坝成因的次生海湾[38]。20世纪90年代以来龙口湾海岸开发活动不断增多，极大地改变了龙口湾的海岸形态和地形地貌[39]，其中以龙口人工岛规模最大。除航道外，湾内水深不足10 m，湾外水深一般为10～20 m（图1）。屺坶岛以南的连岛海岸以基岩海岸或人工岸线为主，人工岛群以南主要为砂质海岸。表层沉积物主要以粉砂和砂为主，湾外沉积物粒径较粗[40]。潮汐性质为不规则半日潮，潮流性质以不规则半日潮流为主，潮流运动形式主要为往复流[41]。波浪以风浪为主，湾内常浪向为SW向，强浪向为WSW向，湾外常浪和强浪向均为NE向[42]。
2. 研究方法
2.1 潮流场数值模拟
本文潮流场模拟采用Mike21中的Flow Model FM HD模块[43]进行，模拟采用非结构三角网格剖分计算域，计算域范围如图2所示，坐标范围为36°59′15.743″～40°59′21.417″N、117°32′22.881″～122°39′30.992″E，覆盖整个渤海海域及部分北黄海海域[44]。人工岛群建设前后除工程区域外其他位置的网格一致，增强了前后对比的准确性和可靠性。
对位于潮滩区的干、湿网格采用动边界的方法进行处理。为能清楚了解研究区所在海域潮流场特征，对该海域网格作局部加密处理。水深资料采用中国人民解放军海军航海保证部制作的海图1∶100万（10011号）、1∶15万（11370号、11570号、11710号、11840号、11910号、11932号）水深资料及研究区附近海域最新实测水深资料。
2.2 波浪场数值模拟
波浪场模拟采用Mike21中的SW浅水波浪数值模块，该模型广泛适用于大范围或者局部区域的波浪预报和分析以及不同历史条件下近海、海岸和港口结构物设计过程中的波浪情况预报[45]。
波浪数值模拟的计算域及网格、岸界和水深资料与本文潮流场数值模拟设置一致，分别模拟了研究区N向和SW向六级风（12 m/s）作用24小时下的波浪场状况，以此来探讨人工岛群建设对龙口湾波浪场的影响。
2.3 纳潮量计算
纳潮量是一个水域可以接纳潮水的体积，海湾的纳潮量不仅是衡量海湾开发价值的一个水文指标，而且也是反映湾内外水交换的一个重要参数[46]，本文对纳潮量的计算采用叶海桃[47]等对三沙湾纳潮量的算法，其公式如下：
$$ P = \Delta H{A_0} + \mathop \sum \nolimits_{i = 1}^n \Delta H_i'{A_i}$$ (1) 式中：
$ P $ 为纳潮量，单位：m3；$ \Delta H $ 为潮差，单位：m；$ {A}_{0} $ 为最低潮位下水域面积，单位：m2；$ \Delta {H}_{i}' $ 为潮滩上第$ i $ 个网格高潮位时的水深，单位：m；$ {A}_{i} $ 为第$ i $ 个网格上最低潮位时的水域面积，单位：m2。由于龙口人工岛最西侧填海区域超出了龙口湾界线（37°32′15.982″～37°40′12.993″N、120°13′26.819″～120°13′43.580″E），无法采用原龙口湾的界线来计算纳潮量，因此本文选取从屺坶岛至石虎咀断面和海岸线封闭的区域为本文计算区域（图1），以下简称“计算区域”。
2.4 水交换计算
水交换是海洋环境科学研究的一个基本命题，水交换率的计算是研究海湾自净能力的重要指标和手段[48]。本文对水交换率的计算采用保守污染物浓度扩散的方法[49]，在Mike21水动力数值模型中，初始时刻将计算区域内示踪剂浓度设为1，外海域设为0，其他水动力条件保持不变，得到水交换率等于（1–浓度值）×100%。水交换计算区域与纳潮量计算区域一致，计算时长30 d。
3. 结果
3.1 模型验证
本文收集了大连、旅顺、鲅鱼圈、曹妃甸、大口河、潍坊港、北隍城、八角、烟台港、龙口港等10个潮位站历史观测资料并采用傅氏分析方法[50]进行调和分析，选用M2、S2、K1、O1四个分潮的调和常数预报出大潮期的潮位变化，同时结合中国海洋大学2017年5月10—11日在工程附近进行的2个站位的潮位观测资料进行潮位验证；采用2017年5月10— 11日（大潮）6个站位27小时海流同步连续观测资料进行流向流速验证，验证结果表明，模拟结果和实测值吻合较好。限于篇幅限制，本文仅列出研究区周边的潮位验证曲线（C4站位和龙口港）和潮流验证曲线（C4、C5和C6站位），见图3和图4。
3.2 潮流场特征
3.2.1 人工岛建设前潮流场
潮流场数值模拟结果表明，在涨急时刻，湾外潮流流向整体上为西南向，受地形影响较小，流向变化较小，湾内潮流方向主要为西南向和南向，在屺坶岛西南侧潮流转为东南向，进入湾内则变为东向，向南逐渐转为西南向；流速在屺坶岛西侧最大，可达0.68 m/s，龙口湾内流速较小，整体小于0.2 m/s，且越靠近岸线流速越小。落急时刻的整体规律与涨急时刻相反，潮流带为明显的往复流形式，湾外潮流流向整体上为东北向，湾内潮流流向在龙口湾内受地形影响发生逆时针旋转，转为西北向至屺坶岛流向湾外；流速在屺坶岛西北侧最大，可达0.8 m/s，龙口湾内流速分布与涨急时刻类似（图5）。
![]()
图 5 人工岛建设前涨急和落急时刻潮流场分布图（大潮期）Figure 5. Current fields during flood and ebb of tide before construction of artificial island (spring tide)3.2.2 人工岛建设后潮流场
人工岛建设后，龙口湾潮流场特征发生了明显的改变。涨急时刻，湾外潮流流向整体上为西南向，流向变化较小，湾内潮流流向较为复杂，整体上围绕人工岛外侧流出湾内，小部分进入岛内水道，在人工岛防波堤内存在一个漩涡；湾外流速较大，为0.2～0.8 m/s，流速由北向南逐渐减小，湾内流速整体小于0.1 m/s。落急时刻潮流流向整体上与涨急时刻相反，湾外潮流流向整体上为东北向，湾内潮流在人工岛南部逆时针转为西北向，在人工岛防波堤处一部分向北流出湾内，一部分转为东方向进入龙口港围绕岸线流向湾外；湾外流速为0.2～0.8 m/s，湾内流速整体上小于0.1 m/s（图6）。
![]()
图 6 人工岛建设后涨急和落急时刻潮流场分布图（大潮期）Figure 6. Current fields during flood and ebb of tide after construction of artificial island (spring tide)3.3 波浪场特征
3.3.1 SW向6级风作用下波浪场特征
人工岛建设前，SW向6级风（风速12 m/s）作用12小时情况下，湾外有效波高较大，大部分区域有效波高大于1.2 m，受SW向风和地形的影响，波高等值线形态向南凸出。湾内有效波高整体小于1.2 m，近岸人工构筑物的掩蔽作用[51]较为明显，在龙口港内以及屺坶岛南侧部分半封闭区域形成波影区，有效波高小于0.4 m；受湾内水深影响，波高等值线变化率较快，且逐渐趋于与岸线平行（图7）。
![]()
图 7 SW向6级风下人工岛建设前后波浪场分布图Figure 7. Wave fields before and after construction of artificial island under southwest wind人工岛建设后，SW向6级风作用12小时情况下，湾外的波浪特征并未发生明显的变化，有效波高整体上大于1.2 m，波高等值线形态向南凸出。湾内有效波高整体减小，屺坶岛南侧局部区域有效波高能达到1.3 m左右，其他区域有效波高小于1.1 m。人工岛的建设使得龙口湾的地形变得更为复杂，除龙口港内和屺坶岛附近形成波影区外，人工岛水道内也形成波影区，波影区有效波高小于0.4 m。
3.3.2 N向6级风作用下波浪场特征
人工岛建设前，N向6级风（12 m/s）作用12小时情况下，湾外有效波高整体较大，大部分区域有效波高大于1.65 m，受N向风和地形的影响，波高等值线形态向南凸出。湾内有效波高整体小于1.2 m，波高等值线形态向东凸出，受湾内水深影响，波高等值线变化率较快，且逐渐趋于与岸线平行（图8）。
![]()
图 8 N向6级风下人工岛建设前后波浪场分布图Figure 8. Wave fields before and after construction of artificial island under north wind人工岛建设后，N向6级风作用12小时情况下，湾外的波浪特征并未发生明显的变化，有效波高整体上大于1.65 m。湾内有效波高整体减小，在人工岛防波堤处波高骤减，形成波影区，湾内有效波高整体小于1.35 m，人工岛建设区域有效波高形态变化较大，人工岛北侧区域有效波高形态变化较小。
3.4 纳潮量
本文分别计算了人工岛建设前后大潮期和小潮期的纳潮量，结果如表1所示，人工岛建设前大潮期纳潮量1.3620×108 m3，小潮期纳潮量9.1227×107 m3，平均纳潮量1.1371×108 m3；人工岛建设后大潮期纳潮量1.1749×108 m3，小潮期纳潮量7.8660×107 m3，平均纳潮量9.8075×107 m3。相比工程建设前后，大潮期纳潮量减小了13.74%，小潮期纳潮量减小了13.78%，平均减小了13.75%，可见人工岛的建设对计算区域纳潮量的影响很大。
表 1 人工岛建设前后纳潮量Table 1. Tidal prism before and after construction of artificial island潮况 建设前纳潮量/m3 建设后纳潮量/m3 变化量/m3 变化率/% 大潮 1.3620×108 1.1749×108 −1.8710×107 −13.74 小潮 9.1227×107 7.8660×107 −1.2567×107 −13.78 平均 1.1371×108 9.8075×107 −1.5635×107 −13.75 3.5 水交换
根据计算区域30 d的水交换率计算结果（图9），人工岛建设前，30 d平均水交换率为62.58%，在计算边界附近水交换能力最强，水交换率为80%～90%，水交换率由边界处向湾内依次递减，越靠近岸线减小速率越快，在龙口港附近最弱，水交换率小于20%。人工岛建设后，30 d平均水交换率59.82%，在计算边界附近水交换率在80%左右，在龙口港附近及人工岛内部水域水交换率迅速减小，龙口湾内水交换率小于20%，局部区域小于5%，人工岛内部水道水交换率为0～80%，南部和西部区域可达80%，中间水道局部区域水交换率小于5%。
![]()
图 9 人工岛建设前后30 d水交换率分布图Figure 9. 30 days water exchange rate before and after construction of artificial island4. 讨论
4.1 人工岛的建设对潮流场的影响
（1）潮流流向变化
人工构筑物的建设会引起潮流流向的偏转[52]，本文中将人工岛建设前后的流场进行了叠加并将变化明显的区域进行了局部放大（图10、图11）。可见人工岛建设后，在涨急时刻人工岛西北侧（A区域）流向顺时针偏转，人工岛北侧区域（B区域）由原来的开放海域变成半封闭海域，潮流流向在此变化较为复杂，整体上顺时针偏转，人工岛西侧（C区域）由于人工岛防波堤的挑流作用形成了一个逆时针旋转的漩涡，人工岛南侧（D区域）潮流流向逆时针偏转（图10）。在落急时刻人工岛西北侧（A区域）流向顺时针偏转，人工岛北侧（B区域）流向整体上逆时针偏转，人工岛西侧（C区域）和人工岛南侧（D区域）流向逆时针偏转（图11）。
![]()
图 10 人工岛建设前后涨急时刻流向对比图（大潮期）Figure 10. Current direction during flood of tide before and after construction of artificial island (spring tide)![]()
图 11 人工岛建设前后落急时刻流向对比图（大潮期）Figure 11. Current direction during ebb of tide before and after construction of artificial island (spring tide)（2）潮流流速变化
本文分别选取了工程建设前后大潮期涨急时刻和落急时刻的潮流场进行了流速对比（图12）。通过对比分析，龙口人工岛建设前后，研究区的潮流场特征发生了较为明显的变化，尤其是在人工岛附近区域变化更为明显，具体表现在：人工岛西部靠近防波堤的区域由于挑流作用[53]流速增大，最大流速变化超过0.4 m/s，变化率可达60%，涨急时刻流速变化相比落急时刻变化较大，距离防波堤越远，流速变化越小。李池鸿[54]等基于Mike21模型对新建的洋浦神头港区南港区码头工程进行了研究，结果表明工程后防波堤口门附近涨、落潮都会因挑流的影响局部流速会变大。人工岛的建设阻挡了部分区域原有潮流的流动，致使潮流流速减小，在人工岛北部，流速整体减小，涨急时刻减小范围大多为0.08～0.16 m/s（变化率30%～70%），极少数区域超过0.16 m/s；落急时刻流速变化较大，流速减小超过0.16 m/s（变化率约70%）的区域大幅增加；在人工岛内部水道以及人工岛南部区域，流速减小，涨落急时刻变化较为相似，流速减小范围整体上为0.08～0.16 m/s（变化率30%～70%），在界河口的一小部分区域流速减小范围超过0.16 m/s。潮流流速的变化可能是引起人工岛北侧、西南侧和内部水道产生淤积以及人工岛西侧发生冲刷[35]的主要原因。
![]()
图 12 人工岛建设前后涨急和落急时刻流速对比图（大潮期）Figure 12. Current velocity during flood and ebb of tide before and after construction of artificial island (spring tide)（3）余流场变化
为了分析人工岛建设对龙口湾余流场的影响，本文选取了潮流数模结果中一个完整大小潮周期（15天）的流速流向计算了人工岛建设前后研究区的欧拉余流（图13），并进行了差值对比（图14）。人工岛建设前，龙口湾湾外除屺坶岛附近区域外余流流向整体上为北向，流速普遍在0.05 m/s左右，在屺坶岛西北侧和西南侧分别形成一个顺时针漩涡和逆时针漩涡，流速为0.05～0.15 m/s，龙口湾西侧存在一个顺时针漩涡，其流速在0.05 m/s左右；湾内余流流向较为复杂，余流从龙口湾西北侧进入后分为两支，一支整体呈逆时针运移在龙口湾北侧流出，另一支呈顺时针运移在龙口湾南侧流出。
![]()
图 13 人工岛建设前后余流特征Figure 13. Residual current characteristics before and after the construction of artificial island![]()
图 14 人工岛建设前后余流流速对比图Figure 14. Comparison of residual flow velocity before and after the construction of artificial island人工岛建设后，湾外余流场未发生明显变化，湾内余流场变化较为明显，具体表现在人工岛西侧形成一个逆时针旋转的新漩涡，且流速增大，增大范围为0.03～0.4 m/s，人工岛北部余流整体上分为两支，一支逆时针运移在屺坶岛附近流出龙口湾，流速较人工岛建设前减小，另一支顺时针运移在人工岛防波堤处流出，流速略有增大，人工岛西南侧区域余流整体上向西南方向运移，流速较人工岛建设前减小。
4.2 人工岛建设对波浪场的影响
海域波浪场的分布特征与水深、地形和风速风向等要素密切相关，在岸线走向、海底地形、风速风向的影响下，波浪在传播过程中会发生一定的折减、绕射等衰减现象[55-56]。本文为了更好地分析人工岛的建设对波浪场的影响，将人工岛建设前后的有效波高做了差值对比（图15），在SW向6级风作用下，湾外有效波高未发生明显变化；湾内变化主要集中在人工岛北部及内部水道区域，其中人工岛内部区域有效波高变化最为明显，人工岛对该区域的掩蔽作用较强，导致有效波高在此明显减小，有效波高减小范围为0.3～1.2 m（变化率24%～96%）；人工岛北部海域人工岛的掩蔽作用较小，该区域有效波高减小，减小范围为0.15～0.6 m（变换率12%～48%）。在N向6级风作用下，湾外有效波高未发生明显变化，湾内变化主要集中在人工岛水道及人工岛西南侧区域，人工岛的掩蔽作用导致该区域的有效波高相对人工岛建设前整体减小，在人工岛防波堤处有效波高减小可达1.2 m，变化率可达96%。
![]()
图 15 人工岛建设前后波浪场对比图Figure 15. Comparison of wave field before and after the construction of artificial island波浪对比结果表明，龙口人工岛的建设导致工程附近区域的波浪有效波高整体减小，其影响范围主要集中在人工岛周边。顾杰[57]等研究结果也表明工程实施前后的波高变化仅集中在工程区域附近，工程的建设使该区域波高显著减小。
4.3 人工岛的建设对潮位的影响
围填海的建设会导致潮位产生一定的变化[58]，为此，我们在人工岛北部、西部、南部以及人工岛内部共选取了12个代表点（图1）分别计算了大潮期一个潮周期内的最大潮差变化，计算结果如表2所示。通过对比分析发现，在人工岛北部潮差减小，最大潮差变化在–0.017 m左右，这主要是由于人工岛的建设使得该区域的潮汐动力减弱，致使潮差变小[59]；在人工岛西侧潮差变化较小，呈现出离人工岛越远变化越小的趋势；人工岛南部潮差增大，距人工岛1000 m的位置潮差增大约0.014 m，远离人工岛变化逐渐变小；人工岛内部水道由于壅水作用[60]潮差变化最为明显，最大潮差变化可达0.047 m。胶州湾跨海大桥的建设对胶州湾潮差的影响普遍小于0.01 m[61]，相比龙口人工岛建设引起的潮差变化较小，其主要原因是龙口人工岛占据的过水面积比例更大，对水流特征的影响更明显。
表 2 人工岛建设前后代表点潮位变化Table 2. Tide changes before and after construction of artificial island (spring tide)位置 站号 工程前最大潮差/m 工程后最大潮差/m 最大潮差变化/m 人工岛北 1 1.080 1.063 −0.017 2 1.076 1.063 −0.013 人工岛西 3 1.095 1.087 −0.008 4 1.100 1.094 −0.006 5 1.109 1.104 −0.005 人工岛南 6 1.107 1.121 0.014 7 1.118 1.126 0.008 8 1.138 1.139 0.001 人工岛内 9 1.085 1.066 −0.019 10 1.092 1.139 0.046 11 1.091 1.070 −0.021 12 1.097 1.144 0.047 4.4 人工岛的建设对纳潮量的影响
人工岛群的建设直接占用了计算区域的海域面积，改变了原有水动力环境，龙口人工岛建设后，本文计算区域的海域面积减少20.68%，大潮期纳潮量减少13.74%，小潮期纳潮量减少13.78%，平均纳潮量减少13.75%。本文整理了前人对于不同海湾纳潮量的研究成果（表3），通过对比发现，莱州湾和罗源湾围填海造成的海域面积减少对纳潮量的影响相对较小[62-63]（纳潮量变化率小于海域面积变化率），锦州湾、芝罘湾、湛江湾以及象山港海域面积减小引起的纳潮量变化较大[64-67]（纳潮量变化率大于海域面积变化率），可见，本文计算区域纳潮量的变化与莱州湾和罗源湾较为相似。总体而言，人工岛群的建设是造成计算区域纳潮量减少的主要原因，这与诸多海湾受围填海影响导致纳潮量减小的结果是相符的。
表 3 不同海湾纳潮量变化对比Table 3. Variation of tide prism in different bays4.5 人工岛的建设对水交换的影响
人工岛群建设后，本文研究区域的地形岸线变得更为复杂，对原有的水交换能力产生了一定的影响，尤其是在人工岛群附近及其内部水道，水交换率有了明显的改变，我们将人工岛建设后的30 d水交换率做了差值对比（图16），在人工岛建设后，计算区域30 d平均水交换率减小2.76%，具体表现在人工岛北部靠近屺坶岛的区域和人工岛南部界河附近区域水交换率明显变大，在人工岛内部水道及北部附近区域水交换率明显减小。本文收集了其他海湾水交换率的变化情况[62, 66-68]（表4），通过对比发现，水交换率的变化与海域面积的变化并无明显相关关系，尤其是湛江湾和锦州湾，在减少更多面积的情况下，平均水交换率变化幅度反而减小。因此，围填海造成的海域面积减小对水交换率的直接影响不大，其引起的水动力环境的改变是影响水交换率变化的主要原因[66]。
![]()
图 16 人工岛建设前后30 d水交换率对比图Figure 16. Variation of 30 days water exchange rate before and after construction of artificial island表 4 不同海湾水交换率变化对比Table 4. Variation of water exchange rate in different bays5. 结论
（1）人工岛建设显著改变了龙口湾潮流场特征，湾内受到人工岛的阻挡，流速普遍减小，局部区域潮流运动形式由往复流变为旋转流，流向变化较大，余流形成多个涡旋；湾外由于堤头挑流作用导致局部区域流速增大和余流流速增大的特征，潮流运动形式未发生明显改变。
（2）受人工岛的掩蔽作用，人工岛及附近区域的波浪有效波高普遍减小，减小幅度主要为0.3～1.2 m，其中在西南向6级风作用下，人工岛北部有效波高减小明显，在北向6级风作用下，人工岛西南部有效波高减小明显。
（3）龙口湾潮位出现北部最大潮差变小、南部最大潮差增大的格局，壅水作用导致人工岛内部水道潮差变化最为明显，最大潮差变化可达0.047 m。
（4）人工岛建设直接占据了龙口湾海域面积，导致其纳潮量明显减小。人工岛建设导致龙口湾水交换率整体减小，呈现出南部和北部增大、人工岛北侧以及内部水道减小的特征，人工岛造成的水动力环境的改变是影响水交换率变化的主要原因。
（5）人工岛建设显著改变了龙口湾水深地形及海湾形态，导致龙口湾内纳潮量减小、潮流、波浪以及水交换等水动力条件普遍减弱，是引起龙口湾水动力条件变化的根本因素。
-
![]()
![]()
图 3 典型SPM颗粒组分扫描电镜照片
a、 b. P1站位，c、 d. P10站位，e. P15站位，f、 g. P39站位，h、i. P42站位。
Figure 3. The SEM images of representative SPM compositions in sites P1(a, b), P10 (c, d), P15 (e), P39 (f, g), and P42 (h, i)
![]()
图 4 P17站位低温κ-T曲线
蓝色线：磁化率随温度变化曲线，橙色线：求导曲线。
Figure 4. Low temperature κ-T curve at site P17
Blue line: magnetic susceptibility curve with temperature, Orange line: derivative curve.
![]()
图 5 SPM代表性样品Loop曲线
蓝色线：顺磁矫正前，粉色线：顺磁矫正后。
Figure 5. Loop curves of representative samples of SPM
Blue line: before paramagnetic correction, pink line: after paramagnetic correction.
![]()
图 6 磁性矿物Day图
SD：单畴，MD：多畴，SP：超顺磁，PSD：假单畴。
Figure 6. The Day plot of magnetic minerals
SD: single domain, MD: multidomain,SP: superparamagnetic, PSD: pseudo-single domain.
![]()
图 7 SPM磁化率空间分布特征
Figure 7. Spatial distribution characteristics of SPM magnetic susceptibility
![]()
图 8 SPM频率磁化率空间分布特征
Figure 8. Spatial distribution characteristics of SPM frequency magnetic susceptibility
-
[1] Kravchishina M D, Lisitsyn A P, Klyuvitkin A A, et al. Suspended particulate matter as a main source and proxy of the sedimentation processes[M]//Lisitsyn A P, Demina L L. Sedimentation Processes in the White Sea. Cham: Springer, 2018: 13-48.
[2] 胡吉连, 杜晓琴. 舟山海域悬浮体的特征及输运机制[J]. 海洋地质与第四纪地质, 2020, 40(6):39-48 doi: 10.16562/j.cnki.0256-1492.2019111304 HU Jilian, DU Xiaoqin. Characteristics and transport mechanism of suspended particles in offshore area of Zhoushan Islands [J]. Marine Geology & Quaternary Geology, 2020, 40(6): 39-48. doi: 10.16562/j.cnki.0256-1492.2019111304
[3] 李文建, 王珍岩, 黄海军. 夏季南黄海悬浮体粒度分布及其影响因素[J]. 海洋地质与第四纪地质, 2020, 40(6):49-60 LI Wenjian, WANG Zhenyan, HUANG Haijun. Grain size distribution pattern and influencing factors of suspended matters in the southern Yellow Sea during summer season [J]. Marine Geology & Quaternary Geology, 2020, 40(6): 49-60.
[4] Stein R. The Late Mesozoic-Cenozoic arctic ocean climate and sea ice history: a challenge for past and future scientific ocean drilling [J]. Paleoceanography and Paleoclimatology, 2019, 34(12): 1851-1894. doi: 10.1029/2018PA003433
[5] Larkin C S, Piotrowski A M, Hindshaw R S, et al. Constraints on the source of reactive phases in sediment from a major Arctic river using neodymium isotopes [J]. Earth and Planetary Science Letters, 2021, 565: 116933. doi: 10.1016/j.jpgl.2021.116933
[6] 汪卫国, 方建勇, 陈莉莉, 等. 楚科奇海悬浮体含量分布及其颗粒组分特征[J]. 极地研究, 2014, 26(1):79-88 WANG Weiguo, FANG Jianyong, CHEN Lili, et al. The Distribution and composition of suspended particles in the Chukchi Sea [J]. Chinese Journal of Polar Research, 2014, 26(1): 79-88.
[7] Kravchishina M, Lein A, Burenkov V, et al. Distribution and sources of suspended particulate matter in the Kara Sea[C]//Complex Interfaces Under Change: Sea-River-Groundwater-Lake. Gothenburg: IAHS, 2014: 42-48.
[8] Wegner C, Hölemann J A, Dmitrenko I, et al. Suspended particulate matter on the Laptev Sea shelf (Siberian Arctic) during ice-free conditions [J]. Estuarine, Coastal and Shelf Science, 2003, 57(1-2): 55-64. doi: 10.1016/S0272-7714(02)00328-1
[9] Rachold V, Eisenhauer A, Hubberten H W, et al. Sr isotopic composition of suspended particulate material (SPM) of east Siberian rivers: sediment transport to the arctic ocean [J]. Arctic and Alpine Research, 1997, 29(4): 422-429. doi: 10.2307/1551990
[10] Stein R, Dittmers K, Fahl K, et al. Arctic (palaeo) river discharge and environmental change: evidence from the Holocene Kara Sea sedimentary record [J]. Quaternary Science Reviews, 2004, 23(11-13): 1485-1511. doi: 10.1016/j.quascirev.2003.12.004
[11] Matul A G, Khusid T A, Mukhina V V, et al. Recent and Late Holocene environments on the southeastern shelf of the Laptev Sea as inferred from microfossil data [J]. Oceanology, 2007, 47(1): 80-90. doi: 10.1134/S0001437007010110
[12] 李秋玲, 乔淑卿, 石学法, 等. 北极东西伯利亚陆架沉积物物源: 来自黏土矿物和化学元素的证据[J]. 海洋学报, 2021, 43(3):76-89 LI Qiuling, QIAO Shuqing, SHI Xuefa, et al. Sediment provenance of the East Siberian arctic shelf: evidence from clay minerals and chemical elements [J]. Haiyang Xuebao, 2021, 43(3): 76-89.
[13] Lien V S, Trofimov A G. Formation of Barents sea branch water in the north-eastern Barents Sea [J]. Polar Research, 2013, 32(1): 18905. doi: 10.3402/polar.v32i0.18905
[14] Schauer U, Loeng H, Rudels B, et al. Atlantic water flow through the Barents and Kara seas [J]. Deep Sea Research Part I:Oceanographic Research Papers, 2002, 49(12): 2281-2298. doi: 10.1016/S0967-0637(02)00125-5
[15] Rozhkova A Y, Dmitrenko I A, Baukh D, et al. Variations in characteristics of the Barents branch of the Atlantic Water in the Nansen Basin under the influence of atmospheric circulation over the Barents Sea [J]. Doklady Earth Sciences, 2008, 418(1): 149-154. doi: 10.1134/S1028334X08010339
[16] Weingartner T J, Danielson S, Sasaki Y, et al. The Siberian coastal current: a wind- and buoyancy-forced arctic coastal current [J]. Journal of Geophysical Research:Oceans, 1999, 104(C12): 29697-29713. doi: 10.1029/1999JC900161
[17] 贾福福, 沙龙滨, 李冬玲, 等. 西伯利亚极地海域第四纪以来古海洋环境研究进展[J]. 极地研究, 2020, 32(2):250-263 doi: 10.13679/j.jdyj.20190074 JIA Fufu, SHA Longbin, LI Dongling, et al. Review of research on Quaternary paleoceanography of the Siberian arctic seas [J]. Chinese Journal of Polar Research, 2020, 32(2): 250-263. doi: 10.13679/j.jdyj.20190074
[18] 田引, 白学志, 黄颖祺. 北冰洋穿极流强度和源头位置变动机制分析[J]. 极地研究, 2021, 33(4):529-544 doi: 10.13679/j.jdyj.20210034 TIAN Yin, BAI Xuezhi, HUANG Yingqi. Analysis of the variation in intensity and source region of the arctic transpolar drift [J]. Chinese Journal of Polar Research, 2021, 33(4): 529-544. doi: 10.13679/j.jdyj.20210034
[19] Deng C, Zhu R, Jackson M J, et al. Variability of the temperature-dependent susceptibility of the Holocene Eolian deposits in the Chinese loess plateau: a Pedogenesis indicator [J]. Physics and Chemistry of the Earth, Part A:Solid Earth and Geodesy, 2001, 26(11-12): 873-878. doi: 10.1016/S1464-1895(01)00135-1
[20] Verwey E J W. Electronic conduction of magnetite (Fe3O4) and its transition point at low temperatures [J]. Nature, 1939, 144(3642): 327-328.
[21] Tauxe L, Bertram H N, Seberino C. Physical interpretation of hysteresis loops: micromagnetic modeling of fine particle magnetite [J]. Geochemistry, Geophysics, Geosystems, 2002, 3(10): 1-22.
[22] Roberts A P, Tauxe L, Heslop D, et al. A critical appraisal of the “day” diagram [J]. Journal of Geophysical Research:Solid Earth, 2018, 123(4): 2618-2644. doi: 10.1002/2017JB015247
[23] Dunlop D J. Theory and application of the Day plot (Mrs/Ms versus Hcr/Hc) 1. Theoretical curves and tests using titanomagnetite data [J]. Journal of Geophysical Research:Solid Earth, 2002, 107(B3): 2056. doi: 10.1029/2001JB000486
[24] Dunlop D J. Theory and application of the day plot (Mrs/Ms versus Hcr/Hc) 2. Application to data for rocks, sediments, and soils [J]. Journal of Geophysical Research:Solid Earth, 2002, 107(B3): 2057. doi: 10.1029/2001JB000487
[25] Thompson R, Oldfield F. Environmental Magnetism[M]. London: Allen & Unwin, 1986.
[26] 敖红, 邓成龙. 磁性矿物的磁学鉴别方法回顾[J]. 地球物理学进展, 2007, 22(2):432-442 doi: 10.3969/j.issn.1004-2903.2007.02.015 AO Hong, DENG Chenglong. Review in the identification of magnetic minerals [J]. Progress in Geophysics, 2007, 22(2): 432-442. doi: 10.3969/j.issn.1004-2903.2007.02.015
[27] Rudenko О, Taldenkova Е, Ovsepyan Y, et al. A multiproxy-based reconstruction of the mid- to Late Holocene paleoenvironment in the Laptev Sea off the Lena River Delta (Siberian Arctic) [J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2020, 540: 109502. doi: 10.1016/j.palaeo.2019.109502
-
期刊类型引用(9)
1. 佀同浩，褚宏宪，孙菲菲，白大鹏，朱龙海，胡日军，林超然. 曹妃甸南部海域表层沉积物侵蚀特征. 海洋地质前沿. 2025(02): 92-104 . 
百度学术
2. 王超，姜胜辉，王世臣，李春川. 山东裕龙岛跨河道工程对沉积动力环境影响研究. 海洋地质与第四纪地质. 2024(01): 191-202 . 本站查看
3. 纪超，姜奇，马殿光，吕彪，周俊伟. 人工岛与多工程影响下的潮流动力响应特征. 科学技术与工程. 2024(05): 1837-1847 . 百度学术
4. 韩慧慧，高飞，丁咚，马乐. 围填海工程对三门湾纳潮量和水交换的影响. 海洋地质前沿. 2024(05): 40-50 . 百度学术
5. 朱金龙，孙伟，陈少伟，冯银银，魏潇，刘宁，徐艳东，王万冠，李姿墨. 岸线变迁对莱州湾纳潮量和水交换的影响研究. 海洋湖沼通报. 2024(03): 12-18 . 百度学术
6. 杨超，朱龙海，张晓东. 基于大量卫星图像研究莱州湾东北部界河海滩近40年的海岸线演变. 海洋学报. 2024(12): 13-25 . 百度学术
7. 聂超辉，吴园园，朱金龙，徐艳东，杜鹏，王菡. 莱州湾西南部岸线变迁对海域水动力影响的数值研究. 海洋环境科学. 2023(02): 228-236 . 百度学术
8. 郭杨子，王世宁，曹晓航. 围堤工程中潮流对抛石漂移的影响研究. 中国港湾建设. 2023(07): 23-29 . 百度学术
9. 张一鸣. 人类活动对沉积环境的影响研究综述. 城市地质. 2022(02): 184-192 . 百度学术
其他类型引用(2)
下载:
计量
- 文章访问数: 828
- HTML全文浏览量: 279
- PDF下载量: 71
- 被引次数: 11
