大气氮沉降对南黄海初级生产过程影响的数值研究*

郭 伟， 史 洁，2** ， 朱君莹， 赵 亮

(1.中国海洋大学海洋环境与生态教育部重点实验室，山东 青岛 266100；2.青岛海洋科学与技术国家实验室海洋生态与环境科学功能实验室，山东 青岛 266071；3.天津科技大学海洋与环境学院，天津 300457)

大气沉降和河流输入是海洋生态系统氮元素的主要外部来源[1]，模式评估结果显示氮、磷等元素的大气沉降量与河流输入量相当，在某些海域甚至大于河流输入量[2-4]。由于经济的快速发展，人为氮氧化物的排放从1980—2000年增加了2.4倍[5]。已有研究表明，西北太平洋边缘海的大气硝酸盐沉降从1980—2003年增加了1.8～4.2倍[6]。Duce等[7]通过生态模型预测，到2030年时，中国近海以及西北太平洋的总大气氮沉降量会在2000年的基础上增加1.0～1.5倍。因此，大气沉降对海区氮循环及初级生产过程的贡献会越来越重要。

黄海位于中国大陆和朝鲜半岛之间，是典型的半封闭陆架海，平均水深44 m。黄海东、西沿岸大约有30多条河流注入，是黄海陆源物质的主要来源，但是对于远离陆地的黄海中部地区，河流的影响较弱，大气沉降是陆源物质向海洋输送的主要途径。Zhang等[8-9]通过观测发现，在西北太平洋近岸海区，包括黄海，尤其是在真光层，由于垂直输运较弱，与河流输送相比，大气沉降已经是黄海营养盐的主要来源之一。此外，Zhang等[10]通过在千里岩岛和嵊泗列岛的观测，定量估算了西北太平洋边缘海和陆架海的大气沉降量，认为大气沉降对海洋溶解无机氮(DIN)的贡献与河流输入同样重要；并且指出大气沉降对西北太平洋近岸海域中的营养物质的生物地球化学循环有着重要作用。也有研究[11]估算黄海全年大气氮沉降量可达494 GgN/a，约是河流输入量的5倍。

大气沉降作为海洋生态系统中氮元素的主要外源[2-4]，对海洋氮循环以及初级生产过程都具有重要影响。大气氮沉降会直接增加海洋氮储库，影响海洋中各营养盐(氮、磷等)含量和比例，以及氮循环路径和速率，改变浮游植物的营养盐限制状况，进而影响初级生产过程[1]。Qi等[11]基于海上站点的观测数据估算黄海的大气氮沉降通量，进而计算得出大气氮沉降可支持黄海0.3%～6.7%的初级生产力；也有研究[12]表明，在全年来看，来自大气沉降的氮支持了黄海17%～37%的新生产力；亚洲沙尘事件带来的营养盐沉降与黄海春季水华的强弱有明显的相关关系[13-15]。然而，由于海区自身营养盐状况的不同，且温度和光照也对浮游植物生长具有限制作用，导致大气沉降的氮不一定完全被浮游植物吸收，进而支持海区的初级生产，因而以往研究直接利用沉降量去估算对初级生产力的贡献往往存在高估。因此，本文利用三维物理-生物耦合模型来模拟大气氮沉降对南黄海主要初级生产过程的影响，并通过数值实验区分不同季节的大气氮沉降的贡献。

1 方法

1.1 物理-生态耦合模型

本研究所运用的三维物理-生物耦合模型主要包括水动力模块和生物地球化学模块。水动力模块基于Princeton Ocean Model[16-17](POM)，模拟区域范围为117.5°E～131.5°E，24.0°N～41.0°N，包括渤、黄、东海，以南黄海为主要研究区域(见图1)，为考虑海海界面和大河影响，将东海和渤海也包括在内。模型的空间分辨率均为1/18(°)，垂直分为21层，对上混合层进行加密[18]。模型的生态模块(见图2)基于NORWECOM[19-20]，生态过程经过修改使得更适合中国海生态过程的特点[21]。生态模块变量包括：大型硅藻(DIA)和小型鞭毛藻(FLA)两种藻类，无机氮(DIN)、无机磷(DIP)和硅酸盐(SIL)三种主要营养盐，以及两种碎屑(碎屑DET和生物硅SIS)。水体中的生物过程包括：浮游植物光合作用生长、呼吸和死亡以及碎屑的再矿化，底栖层中的生物过程包括碎屑的矿化溶出、再悬浮等过程，此外还考虑浮游植物和碎屑从海表到海底的沉降过程以及营养盐的沉积物-水界面交换通量、河流输入和大气沉降。水动力模块为生态模块提供流场、垂直和水平湍流混合系数，用来计算各生态变量的对流输运和湍流混合；提供温度用来计算浮游植物生长、呼吸等过程的温度限制因子。

((a)不同背景色表示不同大气氮沉降数据,沿岸黑色点为河流入海口；(b)研究区域水深分布，单位为m，其中黑色虚线表示36°N断面。(a) Different background colors represent different atmospheric nitrogen deposition data the estuaries are denoted by black dots around the coastline the bathymetry of the research area. (b) The black dotted line is the 36° N section.)

图1 模型模拟区域(a)和研究区域(b)
Fig.1 Model domain (a) and research area (b)

1.2 模型强迫条件及大气输入

物理模块的初始条件和开边界条件均来自大区的模拟结果[18]，包括日均的环流、水温、盐度和海面动力高度等外洋信号和潮汐强迫；大气强迫为2015年日均的NCEP/DOE再分析资料(https://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.html)，包括海表温度、降水、蒸发、气压、热通量和风应力；生态变量的初始条件和开边界条件来自 WOA2005(World Ocean Atlas)数据集和已发表的多年观测资料[22]；河流(主要的10条河流)对营养盐的输入来自已发表的观测资料(2015年河流泥沙公报)与已发表的观测结果[23-28]；计算水中光强所需的悬浮颗粒物(SPM)资料来自SeaWiFS卫星数据(https://oceandata.sci.gsfc.nasa.gov/SeaWiFS/)。

图2 生态动力学模型概念图[21]Fig.2 Schematic illustration of the biophysical model[21]

季节Season渤海Bohai Sea黄海Yellow Sea东海East China Sea干沉降①湿沉降②干沉降①湿沉降②干沉降①湿沉降②冬(1—3月)Winter (Jan～Mar)10.94.34.023.04.34.8春(4—6月)Spring (Apr～Jun)11.034.65.160.14.018.5夏(7—9月)Summer (Jul～Sep)11.113.13.733.12.66.4秋(10—12月)Autumn (Oct～Dec)11.11.54.318.43.53.2

Note:①Dry deposition/mmol·m-2;②Wet deposition/mmol·m-2

1.3 模型实验

本研究共包含5个模型实验，实验初始场、大气动力强迫、河流输入等均相同。模型在不加大气氮沉降的情况下运行两年，第二年末的结果作为所有数值实验的初始场。实验1为在此基础上继续运行一年，模拟得到不加大气氮沉降下2015年南黄海DIN和叶绿素浓度分布的季节变化；实验2为第三年添加大气氮沉降；实验3、4、5分别为第三年的冬、春、夏三个季节分别添加大气氮沉降(月均)。

2 结果

2.1 南黄海无机氮和叶绿素浓度季节变化特征

模型模拟得到东中国海物理场、DIN及叶绿素浓度分布的季节变化。其中，物理场的模拟结果在已发表文献中与实测结果进行了比较验证[32]。图3为模拟得到的南黄海DIN浓度水平(以表层为例)和垂直(以36°N断面为例)分布的季节变化与观测资料[33-36]的对比。可以看出，受近岸河流的排放影响，近岸区域营养盐浓度四季都较高。春季，表层DIN沿岸海域浓度最大可以达到25.0 mmol/m3，由沿岸向黄海中部迅速减小，黄海中部表层由于春季浮游植物旺发消耗了大量营养盐，DIN浓度较低，小于5.0 mmol/m3(见图3(b)、(f))；断面上，黄海中部真光层内DIN浓度较低，下层水水体中DIN浓度高(见图3(j))。夏季，由于水体层化加强，断面上可以看出营养盐跃层较春季有所加强，下层高DIN浓度海水难以向上补充，DIN浓度最小值出现在次表层，浓度小于1.0 mmol/m3(见图3(k))；表层DIN低浓度区域范围进一步扩大，浓度值进一步减小，小于2.0 mmol/m3(见图3(c)、(g))。秋季，由于气温降低、海表风加强，促进了海洋垂直混合，底层富营养盐海水可以向上层海水补充，加之温度降低限制了浮游植物生长的吸收，表层DIN浓度与夏季相比明显增加(见图3(d)、(h))；断面上也可以看出营养盐跃层减弱，次表层DIN低值区基本消失(见图3(l))。冬季，表层DIN浓度则达到了全年最大，黄海中部最大值可以达到5.0 mmol/m3，且低DIN区域范围是四季中最小(见图3(a)、(e))；断面上DIN浓度垂直分布趋向均匀(见图3(i))。与实测对比可以看出，模型可以成功量化DIN的主要生消机制，模拟得到的季节变化特征与实测结果相吻合。近岸水深较浅区域，高DIN浓度范围要大于实际观测的结果，由于本研究主要分析黄海中部区域的两个重要初级生产过程(春季水华和次表层叶绿素最大值)，模型对黄海中部的模拟结果准确性较高，因此以下基于模型结果的分析具有较高的可靠性。

((a)、(b)、(c)、(d)分别表示2、5、8和11月的表层DIN实测数据；(e)、(f)、(g)、(h)分别表示2、5、8和11月的表层DIN模拟结果；(i)、(j)、(k)、(l)分别表示2、5、8和11月的36°N断面DIN模拟结果。(a)～(d) show the observations of surface DIN in February, May, August, and November, respectively; (e)～(h) show the horizontal model results of surface DIN in February, May, August, and November, respectively; (i)～(l) show the vertical model results of DIN in the 36°N section in February, May, August, and November, respectively.)

图3 南黄海表层以及36°N断面DIN季节变化(单位：mmol/m3)
Fig.3 The seasonal variations of DIN concentrations in the surface of the Southern Yellow Sea and in the 36 ° N section(Unit: mmol/m3)

南黄海中部春季表层浮游植物水华和夏季次表层叶绿素最大值现象是两个重要初级生产过程。图4为模拟得到的黄海中部(34°N～36°N，123°E～125°E)表层叶绿素浓度的季节变化图(黑线)，全年看来，叶绿素浓度峰值出现在4—5月，达1.5 mg/m3，其余月份浓度较低，模型成功再现黄海中部春季水华现象。MODIS反演海表叶绿素浓度数据(灰线)也明显看出叶绿素浓度在春季的峰值。但在秋冬季节，MODIS反演结果明显大于模拟结果，这是由于卫星反演的叶绿素浓度受水体中碎屑和有色溶解有机物质等的影响[37]，反演值偏高。尤其是在秋冬季节，由于海表降温和大风影响，水体垂直混合增强，水体中悬浮物浓度增高，导致卫星反演叶绿素浓度很高，而此时低水温限制了浮游植物，实际水中叶绿素浓度很低，因此卫星反演结果在一定程度上存在高估[38-40]。图4中不同标记分别代表已发表的黄海中部现场观测的叶绿素浓度数据[41-45]，可以发现，模拟结果与现场观测基本一致。夏季(见图5)，模型模拟得到的叶绿素垂直分布可以看出(见图5(a))，除近岸区域由于陆源营养盐排放一直支持较高营养盐浓度(见图3(k))和较高的生产力外，黄海中部区域上层水体由于营养盐在春季水华期间耗尽，限制了浮游植物生长，导致叶绿素浓度很低，小于0.4 mg/m3。黄海中部叶绿素浓度最大值出现见次表层，大概40 m水层，叶绿素浓度可达1.8 mg/m3，模拟结果与现场观测结果[46](见图5(b))较为吻合。

图4 黄海中部表层叶绿素浓度季节变化与卫星及观测数据对比Fig. 4 The seasonal variations of the surface chlorophyll concentrations in the central Yellow Sea and the comparison with satellite and observation data

((a)，(b)分别表示模式模拟结果和实测结果，单位：mg/m3。(a) and (b) show the model results and the observations, respectively. Unit: mg/m3.)

2.2 南黄海氮营养盐响应

大气氮沉降的直接作用就是会增加海洋氮储库，其作用在表层尤为明显。图6为考虑大气氮沉降下(实验2)表层DIN浓度与不考虑时(实验1)的年平均差值分布图。从图中可以看出，大气氮沉降导致的黄海表层DIN浓度变化量空间差异显著，近岸海域增加量较大，可以达到3.0 mmol/m3，由近岸海域到黄海中部海域呈明显的递减趋势。这一变化主要决定于原黄海营养盐分布。近岸海域，由于受人类活动影响，河流等营养盐来源众多，导致本身DIN浓度较高，相对于浮游植物生长吸收来说存有富余，故大气沉降的氮在此处不能被浮游植物生长吸收，其效果是显著增加此处DIN浓度。黄海中部区域，本身营养盐不足，尤其是在春季水华期间，大气沉降的DIN是此处营养盐的重要来源，能够及时被浮游植物吸收利用，几乎没有剩余，故从营养盐分布来看，黄海中部差值较小，小于0.5 mmol/m3。

图6 添加大气氮沉降后黄海表层年均DIN浓度变化量水平分布(单位： mmol/m3)Fig.6 The horizontal distribution of the variation in yearly mean surface DIN concentration after adding the atmospheric nitrogen deposition in the Yellow Sea (Unit: mmol/m3)

大气氮沉降进入海洋水体后，有一部分在海洋表层被浮游植物利用，剩余部分则会沉降影响到下层水体的DIN浓度。从两个实验的差值在36°N断面分布来看(见图7)，大气氮沉降引起36°N断面DIN浓度变化量具有较明显的季节差异，冬季(2月：见图7(a))DIN浓度增加量水平和垂直方向都较为均匀，这是由于，浮游植物生长受到冬季低温的限制，因此不能吸收大气沉降的DIN，水体中DIN浓度的增加量基本与大气氮沉降量一致。春季(5月：见图7(b))DIN浓度在近岸海域增加明显，且由于冬季大气沉降的累积作用，春季近岸的增加量较冬季增大；而黄海中部区域，由于春季水华过程消耗了大量DIN，使得上层水体中冬季和春季沉降的DIN都被利用，图7(b)中显示两实验DIN浓度的变化并不明显。夏季(8月：见图7(c))DIN浓度在表层增加多(最大1.0 mmol/m3)，而次表层增加少，这是由于夏季表层水体中限制浮游植物生长的营养盐为无机磷，DIN相对富余，模型结果显示此时表层氮磷比约为40，大气沉降下来的DIN在表层不被浮游植物所消耗而存留；在次表层，氮磷比较小，约为10，大气沉降的氮到达此处被吸收，支持了次表层叶绿素最大值现象，故两个实验DIN浓度在次表层几乎没差别。秋季(11月：见图7(d))，由于水体垂直混合增强，DIN浓度的变化量垂向分布均匀，温度降低使得浮游植物生长较弱，大气沉降的氮累积下来，并被强混合带到下层水体；近岸水域由于全年沉降的氮都得不到利用，累积最为明显，增加量可达2.0 mmol/m3。

((a)、(b)、(c)、(d)分别表示2、5、8和11月的模拟结果。(a),(b),(c),and (d) show the model results in February,May,August,and November,respectively.)

图7 添加大气沉降后36°N断面DIN变化量垂向分布(单位： mmol/m3)
Fig. 7 The vertical distributions of the variations in DIN concentrations after adding the
atmospheric nitrogen deposition in the 36°N section (Unit: mmol/m3)

2.3 南黄海叶绿素响应

南黄海中部主要初级生产过程包括春季表层浮游植物水华和夏季次表层叶绿素最大值现象，图8为大气氮沉降对两个重要初级生产过程的影响。添加大气氮沉降后黄海中部海域春季表层(5月)叶绿素浓度增加的较为明显，增加量最大达到0.20 mg/m3(见图8(a))，可达水华期间峰值浓度的10%，而近岸区域叶绿素浓度增加并不明显。这是由于沿岸区域，由于受河流输入以及人类活动影响较大，海洋表层DIN较为丰富，大气氮沉降量与本底值相比较小，对浮游植物生长影响也就很小；黄海中部浮游植物水华过程消耗了大量营养盐，此时大气氮沉降是营养盐的有效补充，进而使水华强度增大。夏季(8月)，36°N断面上叶绿素浓度的变化量分布可以看出(见图8(b))，次表层的增幅明显，而表层没有明显的增加。海洋表层经历了春季水华过程之后，各种营养盐消耗殆尽，到了夏季，大气氮沉降有利地缓解了DIN的缺乏，但表层磷酸盐由于没有额外补充而成为限制因子，因此大气氮沉降不会促进浮游植物生长，叶绿素浓度增长并不明显；而在次表层，从表层进入的氮沉降到达次表层，缓解了此处的氮限制情况，促进此处初级生产过程，叶绿素浓度增加显著。

图8 添加大气沉降后春季表层(a)与夏季36°N断面叶绿素浓度变化量(b)的空间分布(单位： mg/m3)Fig.8 The spatial distributions of the variations in the chlorophyll concentrations in the surface layer during spring (a) and in the 36°N section during summer (b) after adding the atmospheric nitrogen deposition (Unit: mg/m3)

3 讨论

3.1 冬、春季大气氮沉降对黄海中部春季水华现象的影响

冬季，由于温度和光照限制了浮游植物对水中营养盐的吸收，此时大气沉降的营养盐就会存留在水中，与春季大气氮沉降一起支持浮游植物的旺发。为了区分冬季和春季的大气沉降对春季水华过程的影响，本研究设计了两个数值实验(实验3和4)。图9为分别只添加冬季或春季氮沉降下南黄海春季水华期间表层叶绿素浓度增加量的分布。冬季存留下的氮到了春季促进了春季水华强度(见图9(a))，与不添加大气沉降的结果相比，最大增加量为0.05 mg/m3以上，黄海中部区域平均增加量约为0.03 mg/m3。与冬季相比，春季的大气氮沉降对春季水华的促进作用更加明显，由图9(b)可以看出，增加量大于0.05 mg/m3的区域范围明显大于冬季添加大气沉降的结果，与不添加大气沉降的结果相比，最大值增加量可以达到0.10 mg/m3以上，黄海中部区域平均增加量约为0.05 mg/m3。结果表明，虽然冬季南黄海氮沉降量并不少，但由于冬季大风混合作用，通过海表进入水体中的氮被带到整个水层，到了春季仅有部分留在真光层中支持了水华过程。然而，春季的大气沉降可以直接被浮游植物生长所利用，直接缓解了水华期间营养盐缺乏的情况，对春季水华的促进作用更加明显。

3.2 冬、春、夏季大气氮沉降对黄海中部夏季次表层叶绿素最大值现象的影响

实验3、4和5分别为只在冬季、春季和夏季添加大气氮沉降的模拟，将所得结果与全年不添加大气氮沉降结果对比(见图10)，研究不同季节的氮沉降对夏季次表层叶绿素最大值现象的影响。冬季大气氮沉降一部分在真光层内支持春季水华期间的浮游植物生长，一部分则沉降到达海洋次表层，到了夏季时被浮游植物生长吸收，促进了此时次表层叶绿素最大值现象，与不添加大气沉降的结果相比，次表层叶绿素浓度的最大增加量约为0.03 mg/m3(见图10(a))。春季添加的大气氮沉降则大部分被春季海洋真光层内的水华过程所消耗，沉降到次表层的量较少，所以对夏季次表层的影响比冬季沉降要小，最大增加量约为0.02 mg/m3(见图10(b))。夏季，由于水体层化较强，此时的大气沉降的氮都在上层水体中存留，少部分到达较深的水层，对次表层叶绿素浓度的影响与春季沉降相当，小于冬季沉降(见图10(c))。

图9 在不同季节添加大气氮沉降后春季黄海表层的叶绿素浓度变化量的水平分布(单位： mg/m3)Fig.9 The horizontal distributions of the variations in the surface chlorophyll concentrations during spring after adding the atmospheric nitrogen deposition in different seasons in the Yellow Sea (Unit: mg/m3)

图10 在不同季节添加大气氮沉降后夏季36°N断面叶绿素浓度变化量的垂向分布(单位： mg/m3)Fig.10 The vertical distributions of the variations in the chlorophyll concentrations during summer after adding the atmospheric nitrogen deposition in different seasons in the 36°N section (Unit: mg/m3)

3.3 大气氮沉降对海洋中氮循环过程的影响

大气沉降作为海洋生态系统中氮元素的主要外源，会直接增加海洋氮贮库，影响海洋中营养盐浓度与比例，从而改变海洋氮循环过程及速率[1]。基于加与不加大气氮沉降的模型结果(实验1和实验2)，本研究定量计算了模型中所涉氮循环中的三个过程的改变，即浮游植物生长吸收、呼吸释放以及碎屑矿化。图11为三个过程全年平均水柱积分的变化量的分布图。可以看出，添加大气氮沉降之后，浮游植物生长吸收的无机氮量有明显增加，变化量最大可达2.0 mg·m-2·d-1以上，整个海域平均增大约10%。相应地，浮游植物呼吸释放无机氮量也有相应的增长，最大可达1.5 mg·m-2·d-1以上。浮游植物生长吸收的增加量大于呼吸释放的变化量，说明大气氮沉降对南黄海初级生产过程的促进作用显著。同时，碎屑矿化量也有增加，矿化过程是有机碎屑产生营养盐的过程，而有机碎屑的主要来源是浮游植物死亡，因此大气氮沉降在促进浮游植物生物量的同时，也相应促进了矿化过程，使得浮游植物生物量增长明显的南黄海东部海域矿化量也显著增加，变化量最大可达2.0 mg·m-2·d-1。

((a)，(b)，(c)分别表示浮游植物生长消耗的氮，浮游植物呼吸释放的氮，矿化过程产生的氮。(a), (b), and (c) show the DIN consumption due to growth of phytoplankton, the DIN release due to breathing of phytoplankton, the DIN release due to mineralization, respectively.)

图11 添加大气沉降后年均海洋水柱积分氮通量变化量的水平分布(单位： mg·m-2·d-1)
Fig.11 The horizontal distributions of the variations in the yearly mean vertically integrated nitrogen
fluxes after adding the atmospheric deposition (Unit: mg·m-2·d-1)

4 结语

春季水华和夏季次表层叶绿素最大值现象是南黄海两个重要的初级生产过程，各营养盐外部来源均会对其产生影响，其中大气沉降的作用不容忽视。本文利用三维物理-生态耦合模型定量估算了大气氮沉降的影响。大气氮沉降最明显的作用是增加了黄海氮储量，近岸区域最为明显，向中部海区逐渐递减，这主要是由于近岸海域无机氮来源众多，本身水体中DIN浓度较高，氮不是限制浮游植物生长的主要因子，因此大气氮沉降不会被浮游植物生长吸收，容易出现氮累积。春季水华期间，由大气氮沉降提供的叶绿素浓度的增加达到了峰值浓度的10%，数值实验显示冬季的氮沉降在水体中的存留可以促进接下来的春季浮游植物旺发，而春季的氮沉降影响更加显著。大气氮沉降对夏季次表层叶绿素最大值的贡献为6%，且冬季沉降的存留影响最大，夏季的沉降在水体强层化的情况下难以到达次表层，只留在表层提高了此处的氮磷比。

本研究仅考虑了大气无机氮沉降对海洋初级生产力的影响，同时有研究表明大气沉降的其他物质也对海洋生态系统产生一定影响，比如大气有机氮沉降可能会影响海洋生态系统的结构[47]，大气沉降中的无机磷也会促进特定海区的初级生产力[48]，并且大气沉降中的铁和铜对浮游植物生长有一定的抑制作用[48-50]。除此之外，本研究中模型使用的大气沉降数据为稀疏的观测结果，整个南黄海取值相同，事实上大气沉降存在一定的空间差异，以后的研究中将利用空气质量模型的模拟结果，更加准确的估算大气氮沉降对海洋初级生产力的贡献。