基于珠江虎门河口2016年枯季(1月3~10日)水文水质多要素同步观测数据,采用Pearson相关分析和主成分分析等方法研究河口动力及氮磷营养盐的潮周期变化,并探讨虎门河口营养盐的环境影响机制。结果显示,大、小潮观测期间氨氮浓度变化范围为0.59~2.06 mg/l,硝酸盐氮浓度变化范围为1.20~3.54 mg/l,亚硝酸盐氮浓度变化范围为0.14~0.49 mg/l,磷酸盐浓度变化范围为0.025~0.063 mg/l,且垂向平均浓度分别为1.20 mg/l,1.84 mg/l,0.32 mg/l和0.04 mg/l,揭示虎门河口动力驱动下氮磷营养盐的潮周期差异显著。水文要素盐度和含沙量亦随潮汐涨落变化较大,而水温无明显潮周期变化,且各因子对氮磷营养盐的影响机制不同。外海水团对氮营养盐的稀释混合作用强于磷营养盐,而悬浮泥沙对氮营养盐以吸附作用为主,对磷营养盐吸附与解吸附作用共存。研究结果将为珠江口的污染治理及生态修复提供科学依据。 This study aims to investigate the tidal variations of hydrology and dissolved nutrients, based on extensive field measurements during January 3~10, 2016 in the Humen outlet of the Pearl River estuary, and to indicate the nutrient dynamics. The concentration of NH4 varied from 0.59~2.06 mg/l, 1.20~3.54 mg/l for NO3, 0.14~0.49 mg/l for NO2, and 0.025~0.063 mg/l for PO4 during the sampling time, which demonstrates a significant difference in tidal changes. Salinity and suspended sediment fluctuated with tides, and no big tidal difference was found in water temperature. The role of marine waters on the diluting-mixing process of nitrogen was stronger than that on phosphorus. Compared with nitrogen, phosphorus was scarce in this study, and was a limited factor. The results enable us to develop effective coastal management and to benefit from more focused studies in the functioning of ecological health.
任明吉
中国水利水电第七工程局有限公司,四川 成都
收稿日期:2019年6月26日;录用日期:2019年7月16日;发布日期:2019年7月26日
基于珠江虎门河口2016年枯季(1月3~10日)水文水质多要素同步观测数据,采用Pearson相关分析和主成分分析等方法研究河口动力及氮磷营养盐的潮周期变化,并探讨虎门河口营养盐的环境影响机制。结果显示,大、小潮观测期间氨氮浓度变化范围为0.59~2.06 mg/l,硝酸盐氮浓度变化范围为1.20~3.54 mg/l,亚硝酸盐氮浓度变化范围为0.14~0.49 mg/l,磷酸盐浓度变化范围为0.025~0.063 mg/l,且垂向平均浓度分别为1.20 mg/l,1.84 mg/l,0.32 mg/l和0.04 mg/l,揭示虎门河口动力驱动下氮磷营养盐的潮周期差异显著。水文要素盐度和含沙量亦随潮汐涨落变化较大,而水温无明显潮周期变化,且各因子对氮磷营养盐的影响机制不同。外海水团对氮营养盐的稀释混合作用强于磷营养盐,而悬浮泥沙对氮营养盐以吸附作用为主,对磷营养盐吸附与解吸附作用共存。研究结果将为珠江口的污染治理及生态修复提供科学依据。
关键词 :河口动力,营养盐,主成分分析,虎门
Copyright © 2019 by author(s) and Wuhan University.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
河口区既是陆地与海洋的耦合作用带,又是物质由河流流向海洋的主要运输通道及源/汇区,位置独特,生物地球化学过程复杂 [
珠江河口具有“三江汇流、八口入海”的自然属性,大量氮、磷和有机污染物通过八大口门(图1,虎门、蕉门、洪奇沥、横门、磨刀门、虎跳门、鸡啼门和崖门)进入珠江口并输送至南海。不同河口类型(潮优型、河优型)其动力结构存在显著差异,导致河口生态过程不同 [
珠江虎门河口是东江干流的主要入海口,平均年径流量为6.03 × 1010 m3,年均输沙量为6.58 × 106 t,为珠江八大口门中潮汐动力最强的潮优型河口,潮汐吞吐量为八大口门之首,潮汐为不正规半日混合潮型,潮流为往复流,年均潮差为1.61~1.66 m,平均涨潮流速为1.3 m/s,落潮流速1.4 m/s,多年平均进潮总量为2.29 × 1011 m3,是其下泄径流总量的3.6倍 [
图1. 珠江口位置及虎门采样点布设
水沿虎门上溯至狮子洋,可达距口门42 km的黄埔。涨潮时,潮水由口外向湾涌入,地形束窄,能量聚集,因此潮差增大,流束集中,水流对河床冲刷作用强烈;落潮时,水流退泄,流速加大,使河槽刷深加剧。
在珠江虎门河口进行连续26 h定点观测及样品采集,观测时间为2016年1月3~4日(小潮)和9~10日(大潮)。采用有机玻璃采水器,逐时采集观测期间垂线表(0.2 H,H代表观测时水深)、中(0.6 H)、底(0.8 H)三层水样。水文要素温度、盐度和水深使用optical backscatter point sensor (OBS-3A)现场测定。用于含沙量和营养盐分析的水样,现场抽滤,滤液与滤膜带回实验室分析。所有样品的采集、运输及处理均按照《海洋监测规范》和《海洋调查规范》执行。
本研究采用Pearson相关分析和主成分分析(principal component analysis, PCA)等方法研究珠江虎门河口动力对氮磷营养盐的影响。Pearson相关分析和主成分分析使用SPSS软件完成。
本研究中虎门河口的动力结构主要关注大、小潮涨落过程中流速、水温、盐度和含沙量的潮周期变化,逐时分布如图2所示。
连续观测期间,大潮平均流速为0.49 m/s,小潮平均流速为0.21 m/s (文中实测数据用正值表示落潮流,用负值表示涨潮流)。大潮平均流速明显大于小潮,流速变化越大,水体紊动作用越强,对营养盐的影响也越显著。大潮期间涨潮流速范围为0.06~0.68 m/s,落潮流速范围为0.1~0.91 m/s,落潮历时(16 h)大于涨潮历时(10 h)。小潮涨潮流速范围为0.09~0.39 m/s,落潮流速范围为0.02~0.41 m/s,且涨潮历时(6 h)明显小于落潮历时(20 h)。与大潮相比,垂向流速分布出现明显差异,主要表现在两次涨急时刻。落急时刻,流速出现明显的垂向梯度,表层流速大于底层流速。
图2. 珠江虎门河口水文要素各水层及垂向平均的潮周期变化(小潮:3~4日,大潮:9~10日)
虎门河口观测期间大、小潮水温差异不显著。大潮平均温度为18.76℃,略高于小潮的18.31℃。表层水体温度易受气温影响,本次观测时间为枯季,温度变化范围小,气温低,大潮期间温度垂向分布均匀,小潮期间表层温度略低于底层。
采用Hansen和Rattray所提出的分层系数来分析珠江虎门河口盐淡水混合情况 [
连续观测大潮期间,虎门河口表层悬沙浓度为0.09~0.76 kg/m3,中层为0.11~1.44 kg/m3,底层为0.19~2.49 kg/m3。潮汐涨落过程中,水体悬浮颗粒物含量随流速的增大,水体中含沙量浓度增大,在涨急和落急时刻,因泥沙再悬浮,悬沙浓度出现峰值,中、底层的悬沙浓度明显增大,表层浓度增大不明显。小潮涨落潮过程中表、中和底层悬沙浓度范围分别为0.08~0.26 kg/m3、0.05~0.24 kg/m3和0.09~0.41 kg/m3,表、中层差异较小,底层悬沙浓度较大,变化趋势一致,且变化范围较小,较大悬沙浓度出现在憩流时刻近底层附近。
虎门河口枯季大、小潮期间各项营养盐浓度涨落潮逐时分布剖面图如图3所示。大潮时,PO4、NH4、NO2和NO3浓度的变化范围分别为0.025~0.063 mg/L、0.67~2.06 mg/L、0.14~0.30 mg/L和1.25~3.54 mg/L,垂线平均浓度分别为0.046 mg/L、1.41 mg/L、0.24 mg/L和1.82 mg/L。垂向上,PO4、NH4各层平均浓度均表现为表层 > 中层 > 底层,表层水体PO4、NH4的平均浓度为底层的121%、130%;NO2垂向分布相对比较均匀;NO3表层浓度略低,表层平均浓度比中、底层低10%左右。
小潮期间,PO4、NH4、NO2和NO3浓度的变化范围分别为0.029~0.062 mg/L、0.59~1.74 mg/L、0.26~0.49 mg/L和1.20~2.42 mg/L,垂线平均浓度分别为0.039 mg/L、1.00 mg/L、0.39 mg/L和1.85 mg/L。垂向上,PO4、NH4、NO2和NO3各层平均浓度表现为表层 > 中层 > 底层,表层水体PO4、NH4、NO2和NO3的平均浓度分别为底层的109%、136%、126%和117%。
对比大小潮各项营养盐的浓度,大潮PO4、NH4和NO3垂线平均浓度和各水层的平均浓度均大于小潮,垂线平均浓度分别为小潮的1.24、1.52、1.05倍;NO2垂线平均浓度和各水层的平均浓度则是小潮更大,大潮NO2垂线平均浓度仅为小潮的67%。
海洋浮游生物N/P的正常摩尔比值为12~22,平均值为16,与大洋水体比例相似。水环境可利用的氮与磷在这一比例范围被认为是最适宜浮游生物生长。虎门河口可溶态的氮磷比远超出适宜范围(表1),大潮垂线平均DIN/DIP达177.5 (DIN = NO3 + NO2 + NH4),小潮垂线平均DIN/DIP达188.6。结果显示,该水域氮营养盐充足,其中浮游植物生长可利用的溶解态氮非常丰富,而可利用的溶解态磷显得十分贫乏,表明磷是该水域水环境的主要限制因素,即磷限制,与已有研究吻合 [
图3. 虎门河口大、小潮期间营养盐的时间剖面图(左:大潮,右:小潮)
统计值 | 表层 | 中层 | 底层 | 垂线平均 | |
---|---|---|---|---|---|
大潮 | 最小值 | 92.2 | 113.0 | 138.4 | 120.2 |
最大值 | 265.6 | 536.9 | 272.7 | 352.7 | |
平均值 | 164.4 | 187.8 | 180.3 | 177.5 | |
小潮 | 最小值 | 80.5 | 144.4 | 84.3 | 121.1 |
最大值 | 322.4 | 309.1 | 286.3 | 290.5 | |
平均值 | 199.0 | 199.2 | 167.5 | 188.6 |
表1. 虎门河口枯季可溶性氮磷比
枯水期,河口区以高盐度低营养盐的冲淡水为主体。受径–潮动力影响,水体分层或强烈混合导致营养盐分布及结构发生剧烈的变化 [
虎门河口大、小潮各项营养盐与水文要素间的Pearson相关系数如表2所示。结果显示,小潮PO4与流速、水深、盐度和悬沙浓度的相关性均不显著,而大潮PO4与流速显著相关;小潮NO3以及大潮NO2均与悬沙浓度呈显著负相关,说明悬浮泥沙对以上营养盐以吸附作用为主,悬沙浓度越高,吸附的营养盐越多,导致水体中溶解态营养盐浓度下降;大潮NO3与悬沙浓度呈正相关,说明悬浮泥沙对NO3以解吸附作用为主,悬沙浓度越大,解吸的NO3越多,NO3浓度提升;大、小潮氮营养盐(NH4、NO2、NO3和DIN)均与盐度呈负相关,这是因为高盐度低营养盐的外海水团对其具有稀释作用,盐淡水混合使水体盐度升高而营养盐浓度下降。
流速 | 悬沙浓度 | 盐度 | 温度 | 水深 | ||
---|---|---|---|---|---|---|
小潮 | NH4 | −0.055 | 0.029 | −0.206 | 0.28 | −0.171 |
NO2 | −0.538* | −0.284 | −0.402* | −0.305* | −0.369* | |
NO3 | −0.289 | −0.505* | −0.491* | −0.628* | −0.422* | |
DIN | −0.347* | −0.413* | −0.597* | −0.322* | −0.512* | |
PO4 | 0.156 | 0.19 | 0.146 | 0.577* | 0.121 | |
大潮 | NH4 | 0.488* | −0.392* | −0.09 | 0.165 | −0.46* |
NO2 | −0.012 | −0.585* | −0.466* | −0.036 | −0.32* | |
NO3 | −0.022 | 0.206 | −0.331* | −0.442* | −0.171 | |
DIN | 0.222 | −0.065 | −0.337* | −0.273 | −0.383* | |
PO4 | 0.515* | −0.27 | 0.274 | 0.415* | −0.224 |
表2. 虎门河口水文要素与各项营养盐的相关关系
表中,*代表显著相关。
分别绘制大、小潮悬沙浓度、盐度与DIN的分层散点分布图(图4、图5),探讨虎门河口大、小潮各水层悬沙浓度和盐度对无机氮的影响。可看出,虎门河口小潮期间,各层DIN浓度均随悬沙浓度的增大而减小,表明小潮期间悬浮泥沙对DIN以吸附作用为主;而大潮期间,各层DIN浓度均随悬沙浓度的增大而增大,表明大潮期间悬浮泥沙对DIN主要以解吸附作用为主。虎门河口大、小潮各层DIN浓度均随盐度的增加而降低,体现了盐淡水混合作用下,高盐度低营养盐的外海水团对水体氮营养盐的稀释作用。
为进一步探究虎门河口枯季营养盐的主要驱动因子,分别对大、小潮各项营养盐的影响因素进行主成分分析,结果如表3、表4所示。结果显示,小潮前三个主成分和大潮前四个主成分累计贡献率分别达83.2%和89.0%,可代表变量的大部分信息。
虎门河口大、小潮前两个主成分旋转因子载荷矩阵如图6所示,虎门河口小潮第一主成分主要由水文要素主导(特征值为3.25,占总方差的36.15%);第二主成分主要由PO4和NH4贡献较多(特征值为2.33,占总方差的25.87%)。各个主成分相互独立,揭示了不同水文因子对营养盐的贡献。小潮时水温和含沙量对磷酸盐的影响较
图4. 虎门河口枯季小潮期间悬沙浓度、盐度与DIN分层散点分布
图5. 虎门河口枯季大潮期间悬沙浓度、盐度与DIN分层散点分布
大,而流速对硝酸盐和亚硝酸盐作用较显著。同时显示盐度与磷营养盐的正相关关系,与氮营养盐呈负相关,与Pearson相关分析结构一致。小潮观测期间氮磷营养盐的主要影响因素为盐度、水深、含沙量和水温。
虎门河口大潮第一主成分中PO4和NH4为主要组成部分(特征值为2.50,占总方差的27.78%);第二主成分中盐度贡献较多(特征值为2.10,占总方差的23.29%)。结果显示,盐度和水温对磷酸盐的影响较大,而流速对硝酸盐和亚硝酸盐的影响较突出。大潮观测期间氮磷营养盐的主要影响因素为流速和盐度。综上所述,虎门河口枯季大、小潮动力结构不同,其对氮磷营养盐的影响机制也不同。
主成分 | 初始特征值 | 旋转后的平方和 | ||||
---|---|---|---|---|---|---|
合计 | 方差% | 累计% | 合计 | 方差% | 累计% | |
1 | 4.005 | 44.497 | 44.497 | 3.254 | 36.155 | 36.155 |
2 | 1.996 | 22.181 | 66.679 | 2.328 | 25.869 | 62.024 |
3 | 1.485 | 16.497 | 83.175 | 1.904 | 21.151 | 83.175 |
4 | 0.598 | 6.643 | 89.818 | |||
5 | 0.444 | 4.933 | 94.75 | |||
6 | 0.234 | 2.601 | 97.351 | |||
7 | 0.149 | 1.653 | 99.005 | |||
8 | 0.064 | 0.713 | 99.717 | |||
9 | 0.025 | 0.283 | 100 |
表3. 虎门河口枯季小潮环境因子主成分分析结果
主成分 | 初始特征值 | 旋转后的平方和 | ||||
---|---|---|---|---|---|---|
合计 | 方差的% | 累计% | 合计 | 方差的% | 累计% | |
1 | 3.198 | 35.535 | 35.535 | 2.501 | 27.785 | 27.785 |
2 | 2.471 | 27.455 | 62.99 | 2.096 | 23.293 | 51.078 |
3 | 1.265 | 14.053 | 77.044 | 1.956 | 21.73 | 72.808 |
4 | 1.076 | 11.96 | 89.004 | 1.458 | 16.196 | 89.004 |
5 | 0.47 | 5.219 | 94.223 | |||
6 | 0.284 | 3.156 | 97.379 | |||
7 | 0.138 | 1.535 | 98.913 | |||
8 | 0.061 | 0.683 | 99.596 | |||
9 | 0.036 | 0.404 | 100 |
表4. 虎门河口枯季大潮环境因子主成分分析结果
图6. 虎门河口枯季环境因子第一、第二主成分旋转因子载荷矩阵:(a) 小潮;(b) 大潮
本研究基于珠江虎门河口2016年枯季连续观测数据,分析河口动力和氮磷营养盐的潮周期变化,探讨河口动力对营养盐的控制作用,主要结论如下:
1) 氮营养盐与盐度呈显著负相关,说明虎门河口氮营养盐主要是河流入海的贡献,其分布主要受到沿岸入海河流径流以及外海洋流强烈混合作用及其相互消长影响。
2) 氮营养盐与悬沙浓度总体上呈显著负相关,悬浮泥沙对氮营养盐以吸附作用为主;磷营养盐与悬沙浓度关系不显著,悬沙对磷营养盐吸附与解吸作用并存。
3) 枯季氮磷比远大于海洋浮游生物正常摩尔比值,该水域整体表现为磷限制。
衷心感谢中山大学河口海岸研究所罗向欣、牛丽霞和罗珮珍提供河口基础数据,在论文写作中给予建议和指导并对论文进行修改。
国家自然科学基金青年基金项目(编号:51709289,41506104)。
任明吉. 珠江虎门河口营养盐的分布特征及其环境影响Distribution Characteristics of Dissolved Nutrients and Their Environmental Effect in the Pearl River Estuary[J]. 水资源研究, 2019, 08(04): 371-381. https://doi.org/10.12677/JWRR.2019.84043