长江口无机氮和活性磷酸盐水质基准研究*

武亚男 张昊飞 王 璐 李正炎，

(1. 中国海洋大学环境科学与工程学院 青岛 266100； 2. 国家海洋局东海环境监测中心 上海 201206； 3. 自然资源部海洋生态监测与修复技术重点实验室 上海 201206； 4. 中国环境科学研究院 北京 100012； 5. 中国海洋大学海洋环境与生态教育部重点实验室 青岛 266100)

受淡水径流和海洋潮汐共同作用影响， 河口作为咸淡水的“过滤器”， 盐度变化剧烈， 生态系统复杂敏感(沈焕庭等， 2001； 陆健健， 2003)。氮、磷等营养物质大量输入容易造成河口及其邻近海域富营养化， 藻类大量繁殖， 破坏生态系统完整性(李俊龙等， 2015)。作为我国第一大河口， 长江口紧邻上海、浙江和江苏三个重要经济区， 氮、磷等营养盐的污染持续加重， 近50 年来， 长江口门区域可溶性无机氮(dissolved inorganic nitrogen， DIN)的浓度升高了四倍以上， 可溶性活性磷酸盐(soluble reactive phosphorus， SRP)升高了两倍左右， 尤其在2000 年以后， 赤潮事件暴发频率和面积显著上升， 对生态安全、海水养殖业和人类健康都构成了严重威胁(刘录三等， 2011a； 张璇， 2012； 于仁成等， 2017)。但目前我国尚没有适用于河口环境的水质标准， 现行的《地表水环境质量标准》(GB 3838-2002)和《海水水质标准》(GB 3097-1997)难以满足河口区的管理和保护(孟伟等， 2006)。因此， 制定符合我国区域特点的河口营养物基准十分必要， 这不仅是营养物监测、评价和管理的科学依据， 而且可以为营养物标准的制定提供理论和方法(Wuet al， 2010)。

河口营养物基准(nutrient criteria in estuarine waters)是指河口营养状态参数对生态环境不产生不良或有害影响的最大阈值(USEPA， 1998)。美国系统地研究了河流、湖泊和水库、河口和近岸海域以及湿地等生态系统营养物基准的理论方法， 相继颁布了不同环境类型的营养物基准技术指南， 并对佛罗里达州19 个河口的营养物基准进行了初步研究， 在生态学分区的基础上， 主要运用频数分布法和压力-响应模型法对营养物基准值进行推导和验证， 目前已经形成比较成熟的营养物基准体系(USEPA， 2000a， 2000b， 2001， 2008)。

在这些理论方法的基础上， 我国学者也在探索河口营养基准和标准制定技术。我国河口众多， 河口的分类和分区是营养物基准制定首先需要考虑的因素， 也是河口近岸海域环境管理的重要手段之一(郑丙辉等， 2013)。营养物基准制定过程中， 频数分布法和压力-响应模型是最常用的方法， 前者利用收集到的河口历史和现状数据进行统计分析确定营养物基准值， 后者通过建立营养指标与响应指标的关系来确定营养指标阈值(霍守亮等， 2017)。本文以我国长江口为研究对象， 根据河口内部环境特征开展生态学分区， 在此基础上， 运用频数分布法和压力-响应模型对不同分区的营养指标基准值进行推导， 将赤潮优势藻培养实验得到的生态响应值作为参考。

1 材料与方法

1.1 研究区域和数据来源

长江口是中国的第一大河口， 位处东海和黄海的分界， 整体呈现三级分汊、四口入海的格局， 窄口端江面宽度5.8 公里， 宽口江面宽度90 公里， 其邻近海域相对开阔， 水深变化较为明显， 在123°E 以西的区域内， 水深在60m 以内， 而在123°E 以东的海域， 水深可达到90m (中国海湾志编纂委员会， 1998)。长江口水域是长江径流和近岸海域多种暖流、水团的交汇区， 气象地形条件和生物化学特征都比较复杂(苏纪兰， 2001； Weiet al， 2007)。

我国学者将长江口的范围划定为自江阴(120.26°E)到水下三角洲前缘(123°E)， 按普理查德理论来看， 长江口上缘应为徐六泾(121°E)， 所以本研究的研究区域定为30°30′—32°00′N， 121°00′—123°00′E(沈焕庭等， 2003)。研究数据主要来源于原国家海洋局海洋环境调查数据。长江口盐度空间分布分析主要基于2016和2017 年33 个采样点的表底层数据(图3)； 分区检验运用2018 年3 月、5 月、8 月和10 月96 个站位的调查数据， 调查指标包括水温、pH、溶解氧(DO)、化学需氧量(COD)、总氮(TN)、总有机碳(TOC)等； 营养物基准指标研究运用了2016—2017 年DIN、SRP、pH、DO、悬浮物(SS)、叶绿素a(chla)和盐度(S)数据； 基准推导采用2015—2018 年春夏秋各季节的监测数据。监测站位覆盖了长江口的咸淡水混合区， 所有监测指标的样品采集、预处理、分析与鉴定均按《海洋监测规范》(GB 17378.4-2007)等相关标准方法进行。

1.2 河口生态分区

河口生态分区的出发点是河口生态系统对营养物的敏感性， 河口区域的盐度和地理形态特征是影响其敏感性的主要自然因素(Liuet al， 2018)。地理形态特征是流体动力学和混合因素的产物， 盐度是生物地球化学过程的主要控制参数， 因此， 这些因素被认为是河口生态分区的主要因素(刘录三等， 2011b)。另外， 基于每个采样点的盐度值， 运用Surfer 中的内插法对离散数据格网化获得盐度变化的空间模式。在此基础上， 运用单因素方差分析(One-way ANOVA)， 对水体多参数指标进行分区之间的差异显著性检验， 以验证分区的合理性。

1.3 营养物基准指标研究

结合调查数据， 运用相关性分析(correlation analysis， CA)、主成分分析(principal component analysis， PCA)以及线性回归分析(linear regression analysis， LRA)等方法进行分析。相关性分析是研究两两指标的相关性， 通常用Pearson 相关系数评价相关性； 主成分分析是通过正交变换将一组可能存在相关性的变量转换为一组线性不相关的变量， 转换后的这组变量叫主成分， 以此来实现降维的效果； 线性回归是利用线性回归方程的最小平方函数对一个或多个自变量和因变量之间关系进行建模的一种回归分析。通过以上方法来分析与营养状况显著相关的因子， 进而确定长江口水体营养物基准的关键指标。

1.4 基准推导方法

频数分布法是对营养物数据进行统计分析， 基准值的得出是基于区域内生境状况的判断， 生境状况有完好、部分退化和完全退化三类。参照点一般选取区域内受人类活动损害最小且对该水体或邻近水体的生物学完整性具有代表性的点位， 以此来衡量区域内相对未受干扰的营养状态(USEPA， 1998)。将营养物数据按从小到大的顺序分别排列， 对于生境状况完好的区域， 选择所有参照点数据分布曲线的中值作为基准值； 对于生境部分退化， 但参照点数量充足(大于等于所有观测点数量的10%)的区域， 取参照点营养指标频数分布曲线的上25%数值作为基准值； 对于生境部分退化， 参照点不足(小于所有观测点数量的10%)的区域， 取所有观测点营养指标频数分布曲线的下25%数值作为基准值， 需要指出的是， 对于原始数据大多处于退化状态， 可选取5%作为基准值； 而生境完全退化的区域， 不能使用频数分布法(图1)。

压力-响应模型通过构建原因变量(TN、TP、DIN、SRP 等)与响应变量(chla、DO、TOC、透明度等)之间的关系模型来确定营养物基准值。Chla是反映藻类生物量最常用的指标， DO 是生态健康和栖息地功能评价的重要指标， 受水体扰动、光强和温度、浮游植物及大型动植物的种类和生物量等因素影响， 本 研究将两者作为响应指标以建立压力-响应模型(USEPA， 2010； 李薇， 2014)。

Chla和DO 含量与营养物浓度之间的关系比较细微， 而分类回归树(classification and regression tree， CART)是解决响应变量和原因变量之间非线性、分层及高阶关系的有效方法， 并能测定导致生态发生变化的数字化数值(De'Athet al， 2000)。分类回归树是一种非参数模拟方法， 其通过构建一个二元递归函数， 将由自变量定义的因变量划分为尽可能同质的类别， 每一次划分都由自变量的一次最佳划分值来完成， 分组后响应指标对应的最大降幅或变异性被分成两个独立的节点， 重复此过程， 直到形成一组按照某种标准使节点无法再分的终节点。本研究采用CART 模型来揭示原因指标与响应变量之间复杂的相互作用关系， 从而确定营养物基准值。长江口3—10 月份有发生赤潮的历史记录， 所以对2016—2018 年3—10月的数据进行CART 分析(刘录三等， 2011a)。

以上的营养物基准指标、频数分布统计法和分区检验(One-way ANOVA)的分析和作图在SPSS 24.0 和Origin 9.2 软件中进行， 运用R 语言(R x64 3.5.1)实现CART 的建立。

1.5 生态响应值推导方法

1.5.1 研究方法 本校验试验以中肋骨条藻(Skeletonema costatum)和东海原甲藻(Prorocentrum donghaiense)为研究藻种。中助骨条藻属于温带和热带海洋硅藻类， 有广温广盐的特征， 虽然无毒， 但增殖速率很快， 春夏秋季均是其繁殖高峰期(王江涛等， 2012； 中华人民共和国自然资源部， 1989—2019)。中肋骨条藻是长江口的赤潮种之一， 温度25°C 左右和 盐度14—20 条件下最常见， 是口门区赤潮的主要优势藻(李照， 2018)。东海原甲藻属温带甲藻类， 最适温度为20°C 左右， 最适盐度为25—34， 是我国长江口及其邻近海域发生赤潮事件的关键种， 与中肋骨条藻相比， 东海原甲藻赤潮面积更大， 主要在口外区出现， 持续时间更长， 会导致爆发海域的大面积低氧现象， 对其海洋生物和区域内生态系统产生极其不良的影响(陈翰林等， 2006； 张清春， 2011； 戴鑫烽等， 2014； 黄备等， 2014 )。

根据长江口现场调查监测数据， 设置合适的氮磷浓度， 以S-logistic 1 模型拟合浮游藻生长曲线， 具体形式如下：

式中，Bt代表t时浮游藻的生物量(cells/mL)，Bf代表终止生物量， 即浮游藻生长达到稳定时的生物量(cells/mL)， k 为S-logistic 1 模型参数，t为培养时间(d)，t0是拐点处的培养时间， 即Bt=Bf/2 时的培养时间。

营养物敏感性是指水生态系统对营养物浓度和结构变化响应的敏感程度， 用生态学指标的数值随营养物浓度变化而变化的规律来确定(张远等， 2006)。本研究中浮游藻终止生物量随营养盐浓度的变化规律基本呈“S”形， 也采用 S-logistic 1 来拟合， 具体形式见公式(2)， 生态响应值的计算方法如图2。

图2 生态响应值的计算方法 Fig.2 Calculation method of ecological response value

1.5.2 浮游藻的培养 浮游藻的培养采用人工海水(石岩峻， 2004)， 将其置于2L 的锥形瓶中， 中肋骨条藻的培养海水盐度调至20， 东海原甲藻盐度为30， 以分别接近两分区的实际海水盐度和最适生长盐度， 高压蒸汽灭菌(121°C， 0.1MPa) 30min 后自然冷却， 并充分振荡， 以恢复原有的溶解气体量。两藻种的实验均分两组进行， 分别加入NaNO3和NaH2PO4以设置氮磷的浓度， 氮组固定磷的浓度， 磷组固定氮的浓度， 每个浓度水平设置三个平行样。培养实验初始将生长良好的指数期藻液置于 1 0 0 m L 离心管中， 于1700r/min 的速度下离心10min， 去除上清液保留下 层藻液， 以尽量减少藻液中残留的营养盐。中肋骨条藻和东海原甲藻的相应培养浓度如表1， 其他的营养液均按照f/2 配方加入(Rytheret al， 1962； Guillard， 1975)。培养瓶置于GXZ 恒温培养箱， 设置中肋骨条藻 培 养 温 度 为 25°C， 东 海 原 甲 藻 20°C， 光 强3000—5000lx， 光暗比12︰12， 每天定时摇动三次， 保持气体的溶解量， 防止浮游藻聚集。每天同一时间， 用BD Accuri C6 Plus 流式细胞仪测定浮游藻量， 指标下降时结束。

2 结果

2.1 长江口生态分区

利用插值法， 得出表层和底层的盐度空间格局(图3)。从盐度分布来看， 与底层相比， 表层冲淡水向外海扩展较远， 且不同年份盐度梯度变化差异较大；而底层盐度在不同年份的空间格局非常相似， 河口拦门沙附近梯度变化最剧烈。根据沉积物特点来看， 口门拦门沙附近10m 等深线以内的区域， 冲积滩特征明显， 主要为粗颗粒沙， 而口外区10—60m 等深线的区域， 沉积物较细(陈沈良等， 2009)。而沉积物分布特征可以反映水动力条件、地貌类型以及泥沙来源等因素的影响， 所以是长江口生态分区的重要参考因素(陈沈良等， 2009)。在分析以上因素的基础上， 同时考虑行政管理的便利性， 前哨海岸、佘山和鸡骨礁为几个标志性位置点， 所以将其的弧形连线作为分区的界限。

表1 浮游藻氮组和磷组的培养浓度(μmol/L) Tab.1 Culture concentration of phytoplankton for nitrogen group and phosphorus group

图3 长江口不同年份盐度分布 Fig.3 Saltine distribution of Changjiang (Yangtze) River estuary in different years

图4 长江口分区结果 Fig.4 Preliminary results of classifications on Changjiang River estuary

根据以上结果， 将长江口分为两个区域： I 区(口门区)和II 区(口外区)(图4)。I 区西起徐六泾， 此为冬季盐水入侵锋所及位置， 河口分汊的重要节点； 东至 江苏启东前哨海岸、佘山(31°25.3′N， 122°14.4′E)和鸡骨礁(31°12′N， 122°19′E)大致呈弧形的连接分界线。I区属于咸淡水过渡区， 盐度为7.78±7.88(平均值±标准差)， 90%的采样点盐度值在20 以内， 此区域分布于长江口最大浑浊带， 透明度较低， 强烈受淡水径流的影响， 陆源污染物随长江流域直接排入此区域， 营养物污染严重。II 区西起分界线， 东至123°E 河口锋前沿附近， 盐度为28.76±4.44， 90%的采样点盐度值在33 以内， 盐度变化相对较小， 此区域表层会受冲淡水的影响小， 而底层主要为盐度大于30 的咸水， 受海水冲刷和自净作用， 营养物浓度较低。

首先对两个分区水体表层的水温、pH、DO、COD、TOC 和TN 等进行方差齐性检验， 筛选出方差齐性的数据， 然后进行单因素方差分析， 各指标具体统计结果详见表2。结果显示， 两个分区之间各个参数的差异显著性水平均小于0.01， 差异性显著， 说明该分区基本合理。

2.2 营养物基准指标筛选

通过主成分分析， 得出指标间的相关系数矩阵(表3)、主成分特征值及累积贡献率(表4)。根据主成分分析的基本要求， 发现前三个主成分的特征值大于1， 且累积贡献率达到了74.703%， 所以提取三个主成分， 从而达到数据降维的目的。为了对各主成分相关指标进行分析， 将每个主成分与各指标相关系数进行旋转， 得出旋转后的成分与原始指标变量的因子载荷矩阵(表5)。

表2 2018 年长江口不同分区水体表层各指标基本统计值 Tab.2 Basic statistics values of surface indicators in different segments of Changjiang River estuary in 2018

表3 长江口各指标因子之间的相关系数及检验结果 Tab.3 Correlation coefficients and significance test results of various index factors in Changjiang River estuary

其中第一个主成分的特征值和贡献率最大， 其包含的信息最多。从因子载荷矩阵来看， 与第一主成分关系密切的监测指标有三个， 分别是DIN、SRP 和盐度。它反映了两方面的信息， 一方面是营养指标， 河口区营养盐含量是富营养化评价的重要指标； 另一方面是盐度， 盐度是河口生态分区的关键指标之一， 它在一定程度上反映了河口区域的本底状况， 是生物地球化学过程的主要控制参数(刘录三等， 2011b)。所以将第一主成分定义为营养物基准分区指标和营养指标。

表4 特征值及主成分贡献率与累积贡献率 Tab.4 Eigenvalue， principal component contribution rate and cumulative contribution rate

表5 旋转后的成分矩阵 Tab.5 Rotated component matrix

对长江口各指标因子间进行两两相关性分析， 并对其进行双尾显著性检验(表3)。相关关系表明， DIN、SRP、DO 和SS 与盐度呈显著负相关， 且与DO和SS 相比， DIN 和SRP 与盐度的相关性较强， 而pH与盐度呈显著正相关， chla与盐度相关性不显著。从线性回归图中可以看出， 随着盐度的增加， DIN 和SRP 浓度均逐渐降低， 与SRP 相比， DIN 表现出更好的保守性， 在该地区的其他研究中也观察到了类似的相关性(张平， 2001； 辛明， 2014)(图5)。

图5 长江口营养盐浓度与盐度的相关关系 Fig.5 Relationship between nutrient concentrations and salinity in Changjiang River estuary

2.3 基准值推导

2.3.1 频数分布法 本研究中生境状况的判断运用生态系统健康评价的结果， 健康及以上等级为完好， 亚健康等级为部分退化， 不健康及以下等级为严重退化， 研究表明， 长江口生态系统处于亚健康状态， 即生境部分退化(周晓蔚等， 2011； 中华人民共和国生态环境部， 2019)。参照点的选择参考了生态环境脆弱性评价结果和历史赤潮记录。口门区接近海岸线， 工业、渔业和旅游业等产业发达， 较20 世纪90 年代初期， 区域内生物多样性明显降低， 浮游生物物种密度年际间变幅较大， 生态系统结构不稳定， 且历史记录中佘山附近的水域发生过赤潮事件， 参照点难寻， 所以选取所有监测数据营养指标频数分布曲线的下5%对应值作为口门区的营养物基准值(周晓蔚等， 2011； 何彦龙等， 2019)。口外区虽离海岸较远， 但也受到来自长江径流污染物的影响， 营养物结构和水平发生了显著变化， 此外还受水体层化、上升流加强、黑潮入侵等多种因素影响， 区域内有大面积的赤潮现象， 选取没有赤潮记录以及处于微度和轻度脆弱等级的点位作为参照点， 将参照点营养指标频数分布曲线的上25%对应值作为口外区营养物基准值(刘录三等， 2011a； 于仁成等， 2017)。利用SPSS 24 进行频数分析拟合， 绘制长江口DIN 和SRP 的频数分布统计表(表6)和频数分布直方图(图6， 图7)， 得到口门区DIN 和SRP 基准值分别为0.614 和0.029mg/L； 口外区DIN和SRP 基准值分别为0.300 和0.022mg/L。

2.3.2 压力-响应模型 运用CART 法建立压力-响应模型， 将分类回归树的根节点营养指标浓度值作为基准值。SRP 和DIN 浓度分别是影响区域内chla和DO 浓度变化的首要原因变量， SRP 浓度值0.05705 和0.00845mg/L 分别为图8 口门区和口外区分类回归树的根节点， DIN 浓度值0.6273 和0.2951mg/L 分别为图9 口门区和口外区分类回归树的根节点， 保留三位小数， 得到口门区DIN 和SRP 的基准值分别为0.627和0.057mg/L， 口外区DIN 和SRP 基准值分别为0.295 和0.008mg/L。

表6 长江口2015—2018 年DIN 和SRP 频数分布法统计结果 Tab.6 Statistical results of DIN and SRP by frequency distribution method in Changjiang River estuary during 2015—2018

图6 口门区DIN 和SRP 频数分布直方图 Fig.6 Histogram of frequency distribution of DIN and SRP in the mouth area

图7 口外区DIN 和SRP 频数分布直方图 Fig.7 Histogram of frequency distribution of DIN and SRP in the outside-mouth area

图8 口门区(a)和口外区(b) DIN 和SRP 浓度值对chl a 浓度值划分的分类回归树模型 Fig.8 The CART of values of chl a divided by DIN and SRP concentrations in the mouth area (a) and outside-mouth area (b)

2.3.3 基准推荐值 为了尽量实现对河口生态系统的保护， 选取频数分布法和压力-响应模型结果的低值作为基准推荐值， 最终得出口门区DIN 和SRP的基准推荐值分别为0.614 和0.029 mg/L， 口外区DIN 和SRP 的基准推荐值分别为0.295 和0.008mg/L。与长江口的其他研究相比， 本研究口门区和口外区的基准推荐值均小于其结果， 这一方面是所选用数据不同的原因， 另一方面可能是由于频数分布法的 运用有差别， 这就需要建立更加完整的河口基准数据库， 探索系统的分类分区体系， 完善营养物基准推导方法(表7)( Yanget al， 2019)。

图9 口门区(a)和口外区(b)DIN 和SRP 浓度值对DO 浓度值划分的分类回归树模型 Fig.9 The CART of values of DO divided by DIN and SRP concentrations in the mouth area (a) and outside-mouth area (b)

2.4 生态响应值

根据公式(1)曲线拟合， 得到两种赤潮敏感藻的生长曲线图(图10， 11)。对于中肋骨条藻， 当DIN 浓度固定为120μmol/L 时， 随着SRP 浓度的增大， 浮游藻的最终生物量逐渐增加， 但增加变缓， 固定SRP 为3.0μmol/L 时， 随着DIN 浓度的增大浮游藻的最终生物量也呈现同样的规律。而对于东海原甲藻， 当DIN 固定时为120μmol/L 时， 随着SRP 浓度的增大， 浮游藻的最终生物量逐渐增加， 而当固定SRP为2.0μmol/L时， DIN浓度在0—80μmol/L 的范围内变化时， 浮游藻的最终生物量逐渐增加， 但超过80μmol/L 时， 最终生物量开始下降， 这可能是由于东海原甲藻对磷酸盐的吸收效率高， 氮磷比超过40 则生长受限(张璇， 2012)。

表7 不同方法推导的营养物基准值与其他研究结果的比较 Tab.7 Nutrient criteria derived from various methods in comparison with other research results

图10 中肋骨条藻氮组(a)和磷组(b)的生长曲线拟合图 Fig.10 Growth curves of Skeletonema costatum of nitrogen group (a) and phosphorus group (b)

运用1.5.1 小节中生态响应值的计算方法， 得到中肋骨条藻对DIN 和SRP 的生态响应值分别为0.535和0.025mg/L， 东海原甲藻对DIN 和SRP 的生态响应值分别为0.254 和0.043mg/L。分别将中肋骨条藻和东海原甲藻的响应值作为I 区和II 区的生态响应值， 得到长江口I 区DIN 和SRP 的生态响应值分别为0.535 和0.025mg/L， II 区DIN 和SRP 的生态响应值分别为0.254 和0.043mg/L (表8)。

图11 东海原甲藻氮组(a)和磷组(b)的生长曲线拟合图 Fig.11 Growth curves of Prorocentrum donghaiensis of nitrogen group (a) and phosphorus group (b)

表8 浮游藻对DIN 和SRP 的生态响应值(单位： mg/L) Tab.8 Ecological response of phytoplankton to DIN and SRP (unit： mg/L)

3 讨论

本研究对长江口的生态分区结果符合连续性和便利性原则(Perillo， 1995)。口门区淡水和海水的混合过程非常强烈， 当地的水产养殖业集中在这一区域， 且来自地表径流的污染物在此处排放入海， 所以， 对此区域水文和营养物污染等方面的评估和管理非常关键(林晓娟， 2018)。口外区主要以海水为主， 虽然营养物浓度较低， 但发生赤潮的风险仍然较大， 所以对营养指标的总量控制也很重要(刘录三等， 2011a)。两个分区的动力学特征和地理特征有明显不同， 同时与水行政管理部门对河海的划界非常接近(张华明， 2009)。此外， 单因素方差分析的结果也表明， 两个分区在多项指标中都有显著性差异， 这与河口中的生物、化学和物理过程密切相关， 由此说明此分区的合理性(Teleshet al， 2010)。然而， 河口区域的水质和生物群落还会受河口的气候、水域特征(例如土地利用)和人类活动等多种因素的影响， 此研究初步实现了分区管理， 但对一些重点管控区域缺乏保护， 例如自然保护区、饮用水水源保护区、特别保护海岛、重要渔业海域等， 这就需要在此基础上根据其功能和海洋生态红线等因素， 通过进一步的生态分区或水体功能分级的标准制定来有效评估污染状况并实施污染控制， 实现对河口的管理和保护(Elliottet al， 2002； 席北斗等， 2013)。

本研究主要关注DIN 和SRP， 而未对TN 和TP做相关分析， 这是由我国水环境管理的实际情况决定的， 现行的《海水水质标准》(GB 3097-1997)中， 水质监测与评价中只要求DIN 和SRP， 造成TN 和TP数据缺失。本研究运用了近几年的监测数据对营养物基准值进行推导， 还可以参考历史数据， 将未被开发的状态作为参照状态。20 世纪80 年代长江口基本无富营养化现象， 1984—2015 年间， 长江口及其邻近水域DIN 和SRP 浓度有线性上升趋势， 1988 年是营养盐含量快速增长的起始年， 因此， 将1984—1987 年营养盐要素的平均值作为参考值， DIN 和SRP 的平均浓度值分别为0.0705 和0.0008mg/L (范海梅等， 2015， 2019)。可以发现， 本研究得到的基准推荐值远远大于此参考值， 说明近30 年来长江口营养物的排放已严重超标， 要使水质恢复到接近未受干扰的状态需要大量的时间和努力， 如何把营养物浓度控制在一个恰当的范围是最关键的问题， 需既能保证初级生产量以及渔业的正常发展， 又不发生富营养化造成生态系统的失衡。另外一个值得关注的问题是营养盐结构， 本研究中得到口门区和口外区的基准推荐值氮磷比分别为46.82 和29.65， 与Redfield 比率16 相比， 明显偏大， 特别是口门区， 将近是16 的3 倍(Redfield， 1960)。这与长江口的氮通量大有关， 区域内有较大程度的磷限制， 这与之前的研究一致(李磊等， 2010； 陈慧敏等， 2011)。

为了有效控制富营养化， 需要考虑河口生态系统(尤其是浮游藻)对营养物的敏感性， 而长江口区域中肋骨条藻和东海原甲藻是对营养物最敏感的赤潮藻类(张远等， 2006； 中华人民共和国自然资源部， 1989—2019)。本研究中， 生态响应值是指区域内赤潮优势藻对营养盐浓度产生大量繁殖响应的临界值， 通过浮游藻生物量随营养指标浓度变化的“S”型曲线拐点值得到， 以防止河口区浮游藻的大量暴发。然而， 虽然单物种藻类培养可以初步得出生态响应值， 但未考虑种间竞争、浮游动物等的影响， 所以需要进一步进行多种浮游藻的微宇宙群落试验， 以校验营养物基准值。

我国对长江口的研究相对完整， 本研究以长江口为例探讨了不同生态分区的营养指标基准值， 然而， 为了实现河口及近岸的环境管理， 我们必须对河口区域进行长期的监测和研究， 建立完整的河口基准数据库， 探索科学合理的河口分类和分区方法， 对多项营养物指标进行分析和推导， 以建立完整的营养物基准体系， 从而实现对赤潮等有害事件的防控， 保证河口生态系统的可持续性。

4 结论

(1) 根据盐度分布、地形特点和行政管理的便利性， 将长江口分为两个生态区： I 区(口门区)和II 区(口外区)， 分区检验发现两生态区各指标差异性显著， 分区基本合理。

(2) 基于主成分分析、相关性分析和线性回归分析， 得出DIN 和SRP 是制定长江口营养物基准的关键指标， 而盐度是其重要影响因素之一。DIN、SRP、DO 和SS 与盐度呈显著负相关， 与SRP 相比， DIN 表现出更好的保守性。

(3) 采用频数分布法和压力-响应模型法得到口门区DIN 和SRP 基准推荐值分别为0.614 和0.029mg/L， 口外区基准推荐值分别为0.295 和0.008mg/L。

(4) 通过中肋骨条藻和东海原甲藻室内培养实验， 得到口门区DIN 和SRP 的生态响应值分别为0.535 和0.025mg/L， 口外区DIN 和SRP 的生态响应值分别为0.254 和0.043mg/L。