黄河三角洲周边海域表层沉积物地球化学特征分布及影响因素

雷志宇，刘荣波，胡日军,3，仇建东，朱龙海,3，朱峰，张晓东,3

1.中国海洋大学海洋地球科学学院，青岛 266100

2.潍坊市地下水及地质环境保护重点实验室，潍坊 261021

3.中国海洋大学海底科学与探测技术教育部重点实验室，青岛 266100

4.中国地质调查局青岛海洋地质研究所，青岛 266237

5.潍坊市土地储备中心，潍坊 261045

研究海洋沉积物的地球化学元素含量及其分布规律不但可以指示沉积物的物质来源，还可以为物质搬运以及沉积过程提供依据[1]。黄河三角洲是中国的大型河口三角洲之一，河口及其邻近海域是陆源物质向海输运的关键区域[2]。大量研究表明，物质来源[3]、水动力环境[4]、沉积物粒度等是黄河三角洲周边海域表层沉积物地球化学特征重要的影响因素[5]。蓝先洪等[6]对整个渤海西部区域的地球化学研究发现，黄河为该区域主要的物质来源，水动力环境以及矿物分布对表层沉积物的元素分布具有重要影响。也有研究表明，研究区不同海域的影响因素类型及其作用强度差异明显[7-8]。黄河三角洲北部海域的黄河入海口水道在1976年由刁口河改道为清水沟后，黄河输入的泥沙只有少量向西北方向流入渤海湾南岸[9]。由于缺乏物质来源，废弃老黄河口海域在水动力环境的作用下受到强烈的侵蚀[10]，但该海域依然具有黄河物质特性[11]。位于黄河三角洲东部的现行黄河口海域具有沉积动力条件复杂、冲淤速度快等显著特征，黄河泥沙输入是该海域重要的物质来源[12]，入海泥沙量显著影响着元素的含量与分布[13]，这些物质同时也受潮流、余流等海流作用的影响[14]。也有学者认为，由于黄河输沙量巨大，现行黄河口海域中人类活动的影响并不强烈[15]，因此研究区表层沉积物地球化学研究也常忽略人类活动的影响；然而，在埕北老黄河口海域，陆域工业农业[16]、海上养殖[17]、海上石油开采等人类活动所造成的污染也明显影响着该海域的底质环境[18]。

综上所述，研究区沉积动力环境复杂多变，人类活动明显增强，地球化学特征影响因素众多且不同区域差异明显。关于该区地球化学特征分布规律研究近年来取得了一定的成果，但对自然因素和人类活动共同作用的影响考虑不足，同时对不同区域影响因素的差异也缺乏系统性探讨。因此本文在前人工作的基础上研究并完善了黄河三角洲周边海域表层沉积物常微量元素分布规律及其影响因素。

1 材料与方法

1.1 资料来源

表层沉积物由山东省第四地质矿产勘查院2015年9月和青岛海洋地质研究所2012年7月所采集。海底表层沉积物采用蚌式取样器进行采集，所有站位样品的重量不少于2 kg，少于2 kg时在原地附近重新采样，连续两次采集仍不达标时，移位再取。采集完毕后样品置于聚乙烯袋，做好标记并妥善放置以备实验室的粒度分析以及元素测定，共采集站位180个，具体站位见图1。

1.2 粒度分析

粒度测试在砂粒级组分含量较高的条件下使用筛析法；在粉砂和黏土组分含量较高的条件下使用激光法。筛析法用1/4 Φ间隔分级；激光法采用Mastersizer-2000型激光粒度分析仪（测量范围为0.02～2 000 μm，偏差＜1%，重现性Φ50＜1%）。本文对于沉积物的粒级划分采用乌顿[20]-温德华[21]（Udden-Wentworth）等比制Φ粒级标准，粒度参数的计算采用Folk和Ward图解法计算[22]，采用Folk分类法进行沉积物的分类和定名。

1.3 元素分析

1.3.1 沉积物样品预处理

将50.00 mg烘干后的沉积物放置于聚四氟乙烯溶样内胆中，并加入少量高纯水以润湿样品。在内胆中加入1.50 mL的高纯HNO3和1.50 mL的高纯HF后摇匀并烘干。冷却后取出内胆并置于电热板上蒸至湿盐状，再加入1 mL 的HNO3蒸干。加入 3 mL高纯HNO3（1∶1），0.5 mL Rh内标溶液（1 000×10-9），加盖及钢套密闭，放入 150℃ 的烘箱中保持24 h。冷却后，将提取液转移至干净的聚酯瓶中，使用去离子水后定容待测[23]。

1.3.2 沉积物样品测试以及质量控制

常量元素以及Cr、Zn、Zr的测试使用ICPOES；Ni、Cu、Cd、Pb 的测试使用ICP-MS。仪器采用iCAP6300（电感耦合等离子体发射光谱仪）和X series2（电感耦合等离子体质谱仪）。空白对照使用硝酸溶液（2%），将标准物质（GBW07311和GBW07313）作为标准溶液以绘制标准工作曲线。监控样品使用GBW07309，用于质量控制的监控样品标准值与仪器所获得测定值之间的差值均不大于5%，并且元素测试重复实验的分析结果也都在10%以内。

1.4 数值模拟

1.4.1 模型设置

数值模拟采用MIKE21中的Flow Model FM HD模块，使用标准Galerkin有限元法进行水平空间离散和非结构网格剖分计算域，时间上使用显示迎风差分格式离散动量方程与运输方程。

1.4.2 计算区域设置

本次数值模拟的计算范围如图2所示，计算坐标 范 围 为 36°59 ′15.743 ″~40°59 ′21.417 ″N、177°32′22.881″~122°39′30.992″E。模型采用三角网格进行计算，动边界的处理使用干、湿网格进行。模拟计算中将本次研究区及其周边海域的计算网格进行加密处理，计算三角单元共计37 217个，节点共计20 827个。

图2 研究区计算网格图Fig.2 Gridding of the study area

1.4.3 模型验证

潮流潮位验证点概况见表1，潮流验证引用中国海洋大学于2019年5月31日8时至2019年6月1日9时在研究区附近海域的6个站位进行的26 h同步测流资料中底层流速流向对数值模拟计算结果进行验证。潮位验证基于大连、大口河、北隍城、八角、塘沽、曹妃甸6个潮位站的M2、S2、K1、O14个分潮的调和常数进行调和分析，预测结果用于数值模拟结果的验证。

表1 潮流潮位验证点概况Table 1 Information of tidal current and tidal level verified on site

1.5 数据分析

本文所进行的多元统计分析采用SPSS22.0软件，包括相关性分析、聚类分析以及因子分析。聚类分析采用组间距离法，距离系数的计算拟采用Euclidean距离平方；因子分析采用Kaiser-Meyer-Olkin和Bartlett球度检验，旋转采用最大方差法，最后提取特征值大于1的因子。

2 结果

2.1 表层沉积物粒度分布特征

黄河三角洲周边海域的平均粒径（Mz）范围为3.14～8.64 Φ（图3a），均值为 5.91 Φ，沉积物类型主要以砂质粉砂和粉砂为主（图3b）。埕北老黄河口海域中，近岸的沾利河、草桥沟、挑河和刁口河河口及其沿岸海域的平均粒径Φ值较小，远岸海域的平均粒径Φ值较大。该海域由近岸到远岸沉积物类型变化特征大致为粉砂质砂、砂质粉砂、粉砂、泥、粉砂。这些不同类型的沉积物总体来说在该海域以大致平行于岸线的条状带分布。现行黄河口海域中，现黄河口、老清水沟河口及其沿岸海域的平均粒径Φ值最小，东部远岸海域的平均粒径Φ值相对较高。莱州湾西部海域存在明显的Φ值较高的区域，研究区东南部的平均粒径Φ值低于莱州湾西部海域。粉砂主要分布在现黄河口以东的15 km处，以及莱州湾西部细粒沉积物富集的海域，其他海域大部分属于砂质粉砂。

图3 研究区表层沉积物分布图Fig.3 Distribution of surface sediments in the study area

2.2 表层沉积物常微量元素分布特征

研究区内常量元素含量以SiO2和Al2O3为主，平均含量分别为59.38%和12.30%，其次是CaO和TFe2O3，平均含量分别为6.77% 和4.46%；MnO、K2O、MgO的平均含量较低，分别为2.71%、2.52%、2.37%；TiO2和P2O5的平均含量最低，分别为0.66%和0.09%。微量元素中Zr的平均含量最高，为 273.10 mg/kg；其次为Zn和Cr，分别为79.60 mg/kg和67.45 mg/kg；Ni、Cu和Pb的平均含量分别为27.58、23.73和22.60 mg/kg；Cd的平均含量最低，为0.13 mg/kg。常微量元素含量见表2和表3，元素空间分布特征见图4。

图4 表层沉积物常微量元素含量分布图Fig.4 Distribution in concentration of major and minor elements of the surface sediments

表2 研究区表层沉积物常量元素含量Table 2 Contents of major elements of surface sediments in the study area %

表3 研究区表层沉积物微量元素含量Table 3 Contents of trace elements of surface sediments in the study area mg/kg

2.3 地球化学分区

研究区不同海域的自然环境以及人类活动差异显著[7]，其中物质输入[8]、水动力环境以及人为因素不仅会控制研究区沉积物的理化性质[3]，而且也会直接影响多种元素的分布规律[24]。为了探讨黄河三角洲周边海域表层沉积物地球化学特征的区域性差异，本文利用系统聚类法对研究区进行地球化学分区。聚类分析是常用于识别元素空间分布主要影响因素的一种多元统计方法，其在地球化学的研究中能够指示各常微量元素相互的联系[25]。R型聚类分析包括了本次研究的16种元素（图5）。通过结果可知当选定的判别距离为15时，各化学元素能够被归为4种元素组合。组合1包含Al2O3、K2O 、TFe2O3、MnO 、Cu、MgO 、Cr、Ni、CaO、Pb；组合2 包含Zn、Cd；组合3包含P2O5、TiO2；组合4包含SiO2、Zr。对于选择进行Q型聚类分析的元素，结合R型聚类分析结果以及各元素组分代表意义后选择 TFe2O3、P2O5、Al2O3、CaO、Cd、Ni、SiO2作为地球化学分区（Q型聚类分析）的基础，结果见图6。

图5 R型聚类分析结果图Fig.5 R-type cluster analysis results

图6 表层沉积物地球化学分区结果图Fig.6 Geochemical division of surface sediments

本次研究在地球化学分区结果的基础上将黄河三角洲北部的I-1、I-2、I-3三个区域归入埕北老黄河口海域，黄河三角洲东部以及莱州湾内的II-1、II-2、II-3三个区域归入现行黄河口海域。埕北老黄河口海域中，Al2O3、K2O、TFe2O3、MnO、Cu、MgO、Cr、Ni、CaO在I-1区含量相对较低，在I-2区和I-3区含量相对较高；SiO2、Zr的分布规律与上述元素相反，主要富集于I-1区。现行黄河口海域中，Al2O3、K2O、TFe2O3、MnO、Cu、MgO、Cr、Ni、Pb、Cu、Zn高值区主要位于II-1区东部远岸海域以及II-2区；SiO2高值区主要位于II-3区。

2.4 潮余流分布特征

研究区位于渤海西部，虽然渤海内余流强度较弱，但是渤海与黄海水交换的能力较差，故潮流的往复运动是该区域物质输运的重要动力过程。潮流的非线性运动会产生潮致余流，潮致余流的特征对研究区的物质运输、生态环境等有着重要影响[26]。

本文通过数值模拟提取了2015年9月研究区两个潮周期的潮余流平均值。根据研究区潮余流模拟结果分析可知（图7），在埕北老黄河口远岸海域的潮余流主要指向东北向，近岸海域的余流较为复杂，在沾利河、草桥沟、挑河河口处存在一个逆时针环流。现行黄河口海域中，现黄河口以东的海域潮余流方向也表现为东北向，莱州湾西部域内主要存在一个顺时针环流，并且在老清水沟河口处的流速较大，这与前人的研究结果一致[27]。

3 讨论

3.1 粒度对地球化学元素分布的影响

沉积物粒度是影响各元素分布的重要影响因素，不同粒级沉积物对多种常微量元素的富集作用差异明显，这种现象被称作“元素粒度控制规律”[5]。为了探讨研究区表层沉积物粒度对各元素组分的影响程度，本文通过各元素与平均粒径的皮尔逊相关性系数来进行分析（表4和表5）。研究区的TFe2O3、Al2O3、MnO、MgO、K2O、Cr、Ni、Cu与平均粒径表现出较强的正相关性，沉积物粒度是控制它们分布的重要因素。其分布规律主要表现为I-3区、II-2区以及II-1区以东海域含量较高，这些海域均具有细粒沉积物富集的特征。埕北老黄河口海域的CaO、TiO2和现行黄河口海域的Zn、Cd、Pb与上述元素具有相似的分布规律，也受沉积粒度影响明显。SiO2与上述元素的分布规律相反，在I-1、II-3以及II-1区的河口沿岸海域含量较高，这些区域主要富集粗粒沉积物。埕北老黄河口海域的Zr与SiO2的分布规律也大致相同，受沉积物粒度控制明显。P2O5在两个区域与平均粒径相关性都不明显，表明研究区沉积物粒度并不是P2O5主要的影响因素。该结果表明研究区中大部分元素受沉积物粒度的影响明显。

表4 埕北老黄河口海域表层沉积物元素相关性分析Table 4 Correlation coefficients between values of elements and average particle size in Chengbei Old Yellow River estuary

表5 现行黄河口海域表层沉积物元素相关性分析Table 5 Correlation coefficients between values of elements and average particle size in modern Yellow River estuary

3.2 水动力环境对地球化学元素分布的影响

前人研究表明，研究区表层沉积物地球化学特征受水动力作用的控制明显。埕北老黄河口海域在1974年黄河改道后，沉积物的主要物质来源消失，在潮流的作用下海流会不断淘蚀废弃叶瓣，因此I-1区主要规律表现为沉积物细粒组分被运移，粗粒组分就地沉积[28]。在欧拉余流主要运输方向中，悬浮物质的运输路径和余流方向大致一致[29]。通过本文潮致欧拉余流的模拟结果得出埕北老黄河口海域潮余流主要汇聚区为东北远岸海域（图7）。结合表层沉积物粒度分布特征以及各元素分布规律可以发现，研究区内受沉积物粒度控制较为明显的TFe2O3、Al2O3、MgO、MnO、K2O、Cr、Ni、Cu均在I-2区和I-3区有明显富集。除此之外，在余流流向较复杂的I-1区河口处也存在着逆时针余环流，这些余流漩涡通常也会导致悬浮物质和其他物质的沉积[30]，因此细粒沉积物以及上述元素在该区域也存在不同程度的富集。

图7 欧拉潮余流场分布Fig.7 Distribution of tidal residual flow

现行黄河口海域中，现黄河口处河流感潮段较短、河流潮差较小，黄河水入海时仅在较小的范围与海水混合，黄河所携带的泥沙大量地集中于河口及其附近海域。此外，受河口沙坝的影响，沉积物难以直接输送到远离河口的海域[31]，入海泥沙粗粒级组分易堆积在 河口附近[32]，导致II-1区河口沿岸海域大部分元素含量较低。莱州湾海域沉积物的空间分布受潮流影响，该海域涨潮的潮流流向在莱州湾南岸主要向西，当潮流到达莱州湾西北岸时流向大致变为北向，这将会与现黄河口南部的南西向涨潮潮流相顶托，两种潮流相互影响导致部分细粒组分在II-2区落淤。同时莱州湾海域内落潮的潮流流向大致为东北向，部分细粒沉积物会被落潮潮流往东北向运移[33]。结合潮致余流场的分布可以发现，II-2区周边存在着顺时针余环流，较好地指示了细粒物质运输规律。细粒沉积物富含亲黏土矿物且具有较强的元素吸附能力，导致II-2区内表层沉积物中大部分元素含量较高。因此，本文认为水动力环境主要通过改变研究区沉积物粒度的空间格局进而控制表层沉积物地球化学特征。

3.3 物源输入对地球化学元素分布的影响

元素比值可以一定程度上消除沉积物粒度对各元素分布所造成的影响[34-35]。沉积物中Al含量随粒度的减小而增加，Al的溶解度在天然水的环境下相对较低，且迁移能力相对较弱，因此，Al常被用于元素比值的研究。本分析所选择的主要是受粒度影响较弱以及受人类活动影响通常较为明显的元素（P2O5[36]、CaO、Cd、Zn、Pb[13]）。P2O5/Al2O3、CaO/Al2O3、Cd/Al2O3、Pb/Al2O3和Zn/Al2O3的比值结果如图8所示。

3.3.1 黄河物质的影响

前人研究表明，研究区内黄河是沉积物的主要来源[37]。现黄河口入海泥沙主要来自黄土高原，由于方解石等碳酸盐矿物的含量较高，黄土具有富Ca的特点[38]。因此本文根据CaO/Al2O3结果（图8a）来探讨黄河物质输入对研究区表层沉积物中元素分布的影响程度[39]。通过结果可以看到，越靠近黄河三角洲CaO/Al2O3比值越高，具体表现为由黄河入海口为中心向南北扩散，再向东北逐渐减弱。对于埕北老黄河口海域来说，黄河三角洲河道从刁口河向清水沟河道迁移后黄河沉积物供应中断[40-41]，水动力不断对沉积物进行改造，导致废弃三角洲叶瓣遭受广泛的侵蚀[42]。虽然I-1区长期在水动力作用下被侵蚀，但其表层沉积物依然具有黄河沉积物的物质特性[11]。这表明沉积物的分布特征受水动力和黄河物质输入的共同影响。另外，CaO的分布特征表明黄河物质输入对现行黄河口海域的影响明显大于埕北老黄河口海域。

图8 元素比值分布图Fig.8 Element ratio distribution

3.3.2 人类活动的影响

黄河三角洲地处中国环渤海经济带之中，近年来工农业的发展、海上油气开采以及海上养殖等活动都对该海域的底质环境影响明显[41]。P2O5/Al2O3比值结果（图8b）显示，现行黄河口海域莱州湾内与埕北老黄河海域西部的比值明显高于其他区域。有研究认为，渤海沉积物中自生钙结合磷与原生碎屑结合磷占有较大比重，这些含磷矿无法为生物所利用[43]。而黄河作为该区域主要物质来源，这些矿物主要通过黄河携带的大量泥沙进入莱州湾，导致莱州湾湾内表层沉积物中磷的含量较高[44]。埕北老黄河口海域受黄河物质的直接影响相对较弱，入海河流输入的污水中携带大量的活性磷酸盐通过初级生产力转化为颗粒态磷，最后不断向沉积物中沉积和矿化，造成I-1区表层沉积物中磷的含量较高[45]。同时有研究认为，养殖区的贝类、鱼类等生物摄入的大部分含磷物质大都通过排泄的方式沉积于底质环境中，因此沿岸广泛分布的养殖区域（图9）会导致表层沉积物磷含量偏高[46]。

图9 海上主要开发利用活动分布图Fig.9 Distribution of major offshore development and utilization activities

Pb/Al2O3、Cd/Al2O3、Zn/Al2O3的比值（图8c——e）在埕北老黄河口海域中主要在I-1区较高，而东北部的I-2区和I-3区比值较低。I-1区沿岸的挑河、沾利河、神仙沟等河流及其周边排污口是工业废水及生活废水的重要的排污途径。有研究认为，该区域的Pb、Cd的废水主要来自医药、农药等化工企业及电子、电镀等工业[47]。根据调查，挑河、沾利河两侧分布着多家大型畜禽养殖基地，动物养殖业向饲料中添加含重金属的添加剂会造成养殖业废水中含有一定量的Zn等重金属，这些污染物随着河流运输入海并在河口沉积，导致I-1区表层沉积物中重金属含量偏高[48]。也有研究认为，埕北老黄河口海域中的海上油气开采所排放的污染物也会影响沉积物重金属的分布[49]。该区域海上分布大量钻井平台，油气开采所产生的废弃泥浆中常含有Pb、Zn、Cu等重金属元素[50]，结合埕北老黄河口海域的Cd/Al2O3、Zn/Al2O3分布规律，它们在主要油气开采区也表现出高值，并且出现了部分点源污染现象。对于现行黄河口海域，除了来自黄河直接物质输入之外，莱州湾内部分元素也受湾内其他河流以及人类活动的共同影响[51]。相对于黄河的高输沙量，排放到该区域的污染物会被稀释导致人类活动的影响较弱[52]。

3.4 常微量元素分布主要的影响因素

因子分析方法可以从原始的高维数据中将有相同规律的参数进行整合并提取公因子[53]，在探讨表层沉积物地球化学特征分布的影响因素中有广泛的应用[54]。为了系统地划分研究区表层沉积物各化学元素分布的主要影响因素，本文基于SPSS22.0对各元素组分进行因子分析。埕北老黄口海域中提取了3个主要影响因子，累积方差贡献为80.25%；现行黄河口海域中提取了4个主要影响因子，累积方差为86.70%。本文认为这些因子反映了研究区大部分元素含量信息（表6）。

表6 因子分析结果Table 6 Results of factor analysis

3.4.1 埕北老黄河口海域常微量元素分布影响因素

埕北老黄河口海域的F1中，方差占比为57.68%，其 中Al2O3、TFe2O3、MgO、K2O、MnO、Cu、Ni、Cr表现出了较高的正载荷，CaO表现出了一定的正载荷，SiO2和Zr表现出较高的负载荷。通过Al2O3、TFe2O3、MgO 、TiO2、K2O 、MnO 、Cu 、Ni、Cr、SiO2的相关性分析发现，这些元素之间相关系数较高且与平均粒径的相关性显著，说明这些元素分布的影响因素相似[55]。有研究表明，Al2O3、TFe2O3、MgO、K2O、MnO易赋存于细颗粒黏土矿物中[56]，SiO2主要富集于石英和长石含量较高的矿物。对于Cu、Ni、Cr来说，由于细粒沉积物具有表面积大、有机质含量高、黏粒结合能力强等特点，因此金属离子通常优先与细粒沉积物颗粒相结合。通过上述分析可以得出，F1因子中正载荷代表细粒沉积，负载荷表示粗粒沉积。

F2中方差占比为13.53%，其中Zn、Cd、Pb表现出了较强的正载荷，其他元素组分的载荷不明显。通过相关性分析可知，在埕北老黄河口海域粒度颗粒大小对于Zn、Cd、Pb的影响并不明显。通过Zn/Al2O3、Cd/Al2O3、Pb/Al2O3的分析结果得出，陆域工农业污水[16]、生活废水以及海上油气开采等人类活动对其分布的影响明显[49]，因此本文将该因子划分为人类活动的影响。

F3中方差占比为9.04%，其中P2O5表现出较强的正载荷。前文已经通过P2O5/Al2O3的分布特征分析得到该海域受人类污染物质的输入以及生物作用的影响明显，因此本文将该因子归类于人类活动与生物作用的共同影响。

综上所述，埕北老黄河口海域元素组分分布主要受沉积物粒度的影响，水动力环境是改变沉积物粒度的主要动力因素，沉积物粒度的改变同时又伴随着元素含量的变化，人类活动则进一步改变了部分元素组分的空间分布规律。

3.4.2 现行黄河口海域常微量元素分布影响因素

现行黄河口海域F1中，方差占比45.48%，Al2O3、TFe2O3、CaO、MgO、K2O、MnO、Cu、Ni、Zn、Pb表现出了较高正载荷，SiO2表现为较高负载荷。这些元素大多与粒度相关性较为显著，并且各元素组分相关性较高，说明其来源和影响因素相似[57]。结合前文分析，受黄河物质中含量较高的方解石等碳酸盐矿物的影响，CaO的分布规律受黄河物质的影响明显。因此本文根据上述结果将该因子的归类为受沉积物粒度和黄河物质输入的共同影响。

F2中，方差占比19.86%，Cd和Ni具有较高的正载荷。通过元素的相关性分析可知，该区域内Cd、Ni与粒度的相关性显著，因此本文将这些微量重金属元素归类于沉积物粒度的影响。Cr在F1、F2、F4中均有一定的正载荷，这说明其影响构成复杂，结合其与平均粒径以及其他元素之间的相关性可知，其受沉积物粒度影响显著。

F3中，方差占比11.40%，其中P2O5和TiO2的正载荷最高，其他化学组分的载荷均不明显。通过P2O5/Al2O3的分布特征分析也得出，该海域中表层沉积物的磷主要以黄河输入的碎屑为主。同时，TiO2多存在于陆源碎屑重矿物中，其中Ti以离子的形式迁移极少，主要以重矿物钛铁矿、金红石的形式进行搬运[58]。也有研究表明，磷在地壳中主要以P2O5的形式存在，副矿物以榍石（CaTiSiO4O）、磷灰石和锆石等为主[59]。根据前文的元素相关性分析可知，现行黄河口海域中P2O5与TiO2的相关性显著，说明它们具有相似的来源或行为。由此本文将该因子划分为主要受黄河输入的碎屑物质和重矿物分布的共同影响[6]。

F4中，方差占比为9.96%，载荷较高的主要是Zr。Zr主要赋存于重矿物锆石中，其具有极难溶解、迁移能力较弱的特征。此外，Zr在现行黄河口海域中与P2O5具有一定的正相关性，说明P2O5也可能会与锆石（ZrSiO4）形成含磷重矿物[59]，因此本文认为Zr的分布受锆石等重矿物分布的影响。

现行黄河口海域大多数元素分布主要受沉积物粒度和黄河物质输入的共同控制，其次部分元素还受重矿物等因素的影响。与埕北老黄河口海域相比，现行黄河口海域受黄河物质的影响更为明显，而人类活动的影响相对较弱。

4 结论

（1）研究区可分为6个地球化学区域，大部分元素（TFe2O3、Al2O3、MgO、MnO、K2O、Cr、Ni、Cu）高值区主要位于埕北老黄河口远岸海域的I-2区、I-3区以及现行黄河口海域南部的II-2区、东部的II-1区远岸海域。SiO2高值区主要位于埕北老黄河口近岸海域的I-1区和现行黄河口海域东南部的II-3区。

（2）大部分元素含量与平均粒径的相关性显著，粒度是影响元素分布的直接因素。潮流控制着研究区沉积物的起动及运移，潮余流的汇聚区和环流区既是细粒沉积物的聚集中心，是大部分元素的高值区，水动力环境主要通过改变表层沉积物粒度的空间格局进而控制元素的分布。

（3）埕北老黄河口海域主要受沉积物粒度的影响，其次为人类活动，人类污染物的输入对Cd、Zn、Pb、P2O5的影响最为显著；现行黄河口海域主要受沉积物粒度和黄河物质输入的共同控制，黄河物质对CaO、P2O5的分布影响显著。自然因素是研究区表层沉积物地球化学特征分布的主要原因，人类活动则进一步改变了部分元素组分的空间分布规律。