工业革命以来，人类活动向大气释放活性氮（Nr）的数量和强度急剧增加，这改变了氮循环的自然过程，进而诱发了一系列生态环境效应。据估算，从1860年到2005年，全球人类活动排放的Nr由15 Tg N·a^-1增加到187 Tg N·a^-1，这些Nr有一半以上通过大气沉降进入地表，在数量上相当于全球陆地生物固氮量（90~130 Tg N·a^-1）^[1-4]。过量的氮沉降是全球变化的重大问题，会导致生物多样性降低、水体富营养化和森林氮饱和以及土壤酸化等诸多不利影响^[5-7]。

氨气/铵盐（NH₃/NH₄⁺，NH_x）是Nr的重要组成部分，而农业活动是区域大气中NH_x的重要来源，这些Nr经过一系列化学反应会造成大气颗粒物污染，危害人体健康。最终，大气中的NH_x会通过降水、气溶胶和气体干沉降三种途径返回地表，即氮沉降过程^[8-10]。中国和美国的观测研究表明，NH_x是大气氮干湿沉降的主要形态^[8]。我国以往氮沉降研究主要集中于湿沉降观测，干沉降研究多以模式估算为主。干沉降观测研究较少且尚未考虑NH₃在大气/植被界面的交换过程，对氮沉降总量（湿沉降、颗粒物和气体干沉降总和）存在系统性高估的可能^[8]。因此，研究大气/植被界面NH_x的动态变化，有利于明确植被NH₃排放和吸收在农田氮循环中的作用，对于深入研究农田大气氮沉降通量的测算方法及变化规律具有重要的科学意义。

植被叶片是大气NH₃的源或汇取决于环境和叶片气孔间气态NH₃的浓度差。当气孔内NH₃浓度与大气环境相同，叶片与大气间无净交换，此时的NH₃浓度称为气孔NH₃补偿点^[11]。补偿点受植物种类、生长阶段和土壤含氮量等因素影响^[12-14]。国际上对植物与大气间NH₃交换的方向、强度、发生过程和生理机制等方面开展了系统研究，但植被大气界面NH₃交换对NH_x沉降通量的影响还鲜见报道^[15]。

华北平原人口密集且农业集约化程度高，经过近几十年城市化的迅速发展，已经成为中国乃至全球大气Nr排放和沉降的热点地区^[8]。本研究选取了区域内典型玉米农田生态系统，探索研究大气/植被界面NH_x的动态变化，将植被冠层大气假设为一个闭合系统（图 1），在没有明显传输源影响的情况下，系统的氮输入和输出主要包括：①植被与大气交换的输入和输出过程；②气态和颗粒态氮干沉降的输入；③气态和颗粒态氮湿沉降的输入；④土壤向大气排放的输出。本文针对以上过程开展综合观测实验，同步采集并测量了植被冠层大气中NH₃和NH₄⁺（气溶胶和降水）浓度、玉米叶片和土壤中NH₄⁺含量，分别根据玉米叶片的氨气补偿点、植被冠层的气态及颗粒态氮浓度、分段降水样品中的氮浓度及土壤剖面中氮随土壤深度的变化来评估植被冠层大气的氮输入和输出，探讨植被大气界面NH₃交换对NH_x沉降通量的影响，为深入了解农田大气-植被-土壤氮循环过程提供科学参考。

1 材料与方法 1.1 采样点介绍

采样点选择在河北香河南郊农田区，属暖温带大陆性季风气候，多年平均气温11.5 ℃，多年平均降雨量近600 mm，降雨集中在夏季，无霜期180 d左右。开展实验的农田区域在6月20日种植玉米，9月30日左右收获。本研究于2016年8月和9月开展了两次综合观测实验，分别处于玉米的抽穗期和花粒期。观测前一年（即2015年8月到2016年8月），农田共施肥两次：一次为农业堆肥，施肥日期为2015年11月15日，含氮量约18%，施肥量约375 kg·hm^-2；另一次为复合肥，施肥日期为2016年6月20日，含氮量33%，施肥量约450 kg·hm^-2。两次施肥输入到农田土壤的氮总量约430 kg N·hm^-2·a^-1。

1.2 实验观测和样品采集

植被冠层的大气Nr浓度（NH₃、气态硝酸HNO₃、NH₄⁺和颗粒态硝酸盐NO₃^-）采用Delta系统测定（Denuder for long-term atmospheric sampling，英国生态水文中心）。采样持续时间为24 h（上午9:00至次日上午9:00）。Delta采样系统抽取空气的流速0.3~0.4 L·min^-1，采样泵与流量表相连以计算采样体积（图 1）。空气进入采样链后，气态和颗粒态Nr分别被扩散管和滤膜吸附/收集。收集HNO₃的扩散管长15 cm，涂抹1%KOH+1%甘油甲醇溶液作为吸附剂。吸附NH₃的扩散管长10 cm，吸附剂为5%柠檬酸甲醇溶液。上层滤膜收集NO₃^-，吸附剂为5%KOH+10%甘油甲醇溶液；下层滤膜涂有13%柠檬酸甲醇溶液用以收NH₄^+[16]。8月10—13日和9月15—17日观测期间共获得5条Delta采样链用以测量大气Nr浓度，每条采样链的测定可以同步得到HNO₃、NH₃、NO₃^-和NH₄⁺四种成分的浓度。HNO₃、NH₃、NO₃^-和NH₄⁺浓度的准确度分析见参考文献[17]。

在观测植被冠层大气Nr浓度期间，同步采集玉米植株叶片，清晨傍晚各一次，用以考察植株NH₃挥发的昼夜变化特征。选取距Delta采样系统0.5、1.0 m和1.5 m的3株玉米，每株分别采集上（距地面1.5~ 1.7 m）、中（0.9~1.1 m）、下（0.5~0.7 m）3个高度的叶片。8月10—13日（上午10:00和下午8:00）和9月15—17日（上午6:00和下午6:00）观测期各采集5次植株叶片，共获得90个样品。

针对观测期间8月12日的降水过程，用直径36 cm的PVC容器间隔10~30 min分段收集降水样品，用量筒量取样品体积后保存于50 mL的PET瓶。此次降水事件前期每10 min收集一次，后期每30 min一次。共采集6个分段降水样品，样品冷藏保存，采集3 d天内即进行化学分析。

在距离Delta采样点3 m的农田取深度为50 cm土壤剖面，以5 cm为间隔划分土壤层，每层取一个混合均匀的土壤样品。另在农田中随机采集了3个表层土壤样品，共采集土壤样品13个。

1.3 化学分析

Delta采样链的NH₃扩散管和NH₄⁺滤膜用去离子水浸提，HNO₃扩散管和NO₃^-滤膜用H₂O₂浸提。浸提液用离子色谱仪测定NH₄⁺和NO₃^-浓度，计算得到植被冠层大气中不同形态Nr浓度。

玉米叶片样品先用去离子水清洗表面附着物，晾干后称取1.5~2.0 g剪碎研磨，用30 mL去离子水浸提可溶性成分。浸提液用2 μm滤膜过滤后取一部分测定pH，另一部分用离子色谱仪测定NH₄⁺浓度。

分段降水样品取一部分测定pH和电导率（EC），另取一部分用离子色谱仪测定NH₄⁺和NO₃^-浓度。

土壤样品自然晾干后，称取约0.05 g溶解于40 mL去离子水中，过滤后用离子色谱测定NH₄⁺和NO₃^-浓度。离子色谱仪的测定方法和检出限等见参考文献[18-20]。

1.4 玉米农田的NH₃补偿点估算

大气/植被界面NH₃双向交换过程取决于植株叶片中可排放的NH₃浓度：当大气NH₃浓度高于叶片NH₃补偿点时，植株叶片吸收NH₃；相反，当大气中NH₃浓度低于补偿点时，植株叶片向大气中释放NH₃。NH₃补偿点的计算公式为^[21]：

式中：[NH_x]和pH分别为叶片质外体中NH_x浓度和质外体pH值，本研究用整个叶片中的对应物质浓度和pH值进行NH₃补偿点估算。R为气体常数8.207 5× 10^-2 atm^-1·mol^-1·K^-1，T为K氏温度，γ_NH3和γ_NH4⁺分别为NH₃和NH₄⁺的活度系数，这里取为1，亨利常数H和NH₄⁺的电离常数K_NH4⁺均与温度有关，分别依据以下两个公式计算：

本研究利用整个叶片组织的[NH_x]浓度和pH值替代质外体估算了NH₃补偿点，该估算方法具有一定的不确定性。为了考察方法的可靠性，我们后续进行了叶片组织和质外体方法的同步测算。结果表明，叶片组织的[NH_x]浓度低于质外体，会造成补偿点的低估；而pH值高于质外体，会造成补偿点的高估；整体上利用叶片组织估算的补偿点略高于质外体方法（7.2%），两者的变化趋势一致。此外，采用气温代替叶温在一定程度上会低估植物叶片的氨气补偿点。

2 结果与讨论

大多数植物具有“双峰型”的NH₃挥发趋势，即在开花期和成熟期各出现一次NH₃挥发高峰^[22-23]。考虑到植株与大气交换活动较为强烈的时间段，本研究选择在玉米抽穗期（8月10—13日）和花粒期（9月15— 17日）分别开展综合观测实验。从国家气象局获取了2016年8—9月香河气象参数的时间变化序列。8月平均温度、相对湿度和降水量分别为25.7 ℃、77.6%和78.7 mm，这些气象参数在8月10—13日研究时段分别为28.5 ℃、82.6%和47.0 mm；9月平均温度、相对湿度和降水量分别为20.7 ℃、77.6%和15.7 mm，在9月15—17日研究时段分别为22.8 ℃、82.6%和1.0 mm。可以看出，两个研究时段代表了区域内夏秋季典型的气象条件。

2.1 植株叶片NH₄⁺浓度变化及其与大气NH₃交换

叶片NH₄⁺含量从玉米植株顶端到底部呈逐渐下降趋势，8月10—13日采集的3个高度叶片NH₄⁺含量从上到下依次为0.25、0.21 mg·g^-1和0.17 mg·g^-1，9月15—17日采集3个高度的叶片NH₄⁺含量从上到下依次为1.20、0.99 mg·g^-1和0.74 mg·g^-1。这种差异反映了同一植株不同部位的NH₃挥发量可能不同^[24]。清晨和傍晚叶片NH₄⁺含量分别受到夜间和白天NH₃交换过程的直接影响，反映了植株NH₃交换的昼夜变化特征。8月10—13日清晨采集的植株叶片NH₄⁺浓度整体上高于傍晚，反映出植株NH₃挥发具有昼高夜低的特征，可能受控于叶片NH₄⁺同化和气孔导度的变化^[25]。白天温度较高，植物光合/蒸腾作用活跃，植株气态NH₃分压升高^[26]，叶片NH₃挥发明显高于夜间，造成傍晚采集叶片的NH₄⁺浓度较低，而清晨采集叶片的NH₄⁺浓度较高。此外，植株叶片NH₄⁺浓度也在一定程度上受到采样时刻温度的影响。叶片NH₄⁺浓度与温度的相关性分析表明，随着温度升高，叶片NH₄⁺浓度整体呈增高趋势，两个研究时段的决定系数R²分别为0.67和0.79（P < 0.05）。

为进一步考察植株叶片与大气中NH₃的交换过程，测定了植被冠层的NH₃浓度和植株叶片中NH₄⁺浓度，并计算了植株叶片NH₃补偿点（图 2）。可以看出，8月10—13日植株NH₃补偿点为1.9±1.4 μg·m^-3，远低于大气NH₃浓度31.4±4.8 μg·m^-3，整体呈现叶片从大气吸收NH₃的趋势，NH_x沉降通量测算若不考虑植被交换过程将被显著低估；9月10—17日植株氨气补偿点52.0±30.7 μg·m^-3，高于大气NH₃浓度29.3±11.7 μg·m^-3，整体呈现叶片向大气排放NH₃的趋势，NH_x沉降通量测算若不考虑植被交换过程将被显著高估。NH₃补偿点很大程度上受pH驱动，花粒期叶片pH值平均为6.4，高于抽穗期的平均值5.3，导致植株的花粒期NH₃补偿点高于抽穗期。此外，作物生长后期植物叶片以氮素转运为主，尤其是到了衰老期，蛋白质水解产生大量的NH₄⁺，为叶片NH₃挥发提供了直接来源，也导致较高的叶片NH₃补偿点。

2.2 植被冠层大气NH₃和NH₄⁺气溶胶浓度及氮的干沉降量估算

8月10—13日和9月15—17日两次观测期间植被冠层大气NH₃浓度的变化范围分别为27.1~36.6μg·m^-3和21.1~37.6 μg·m^-3（图 3），这与北京周边农田地区夏季NH₃浓度比较接近（14.9~29.8 μg·m^-3和27~ 41 μg·m^-3）^[27-28]，而远高于欧美一些农田地区（低于10 μg·m^-3）^[29-31]。8月10—13日植被冠层大气NH₄⁺的变化范围为2.1~8.0 μg·m^-3，较大的波动范围主要与湿清除造成NH₄⁺浓度的降低有关。8月12日13:00— 15:00出现一次短时降水过程，植被冠层的NH₃浓度没有明显下降（30.1 μg·m^-3），可能是植被或土壤排放补充了被降水清除的NH₃。降水结束后颗粒物NH₄⁺浓度却由8.0 μg·m^-3显著降低到2.1 μg·m^-3，表明降水对颗粒物NH₄⁺的湿清除效果大于气体NH₃。

利用农田Nr年平均干沉降速率（Vd）^[9]与本研究期间测量的气态和颗粒态Nr浓度，估算本研究期间（8月10—13日，3 d）NH₃、HNO₃、NH₄⁺和NO₃^-的干沉降量分别是33.5、2.6、2.5 mg N·m^-2和0.9 mg N·m^-2；干沉降中还原性氮与氧化性氮的比值是10:1。

2.3 大气NH_x的湿清除过程及氮沉降总量估算

观测期间分段收集的降水中NH₄⁺浓度整体呈现下降趋势（图 4）。降水开始后的前10 min，NH₄⁺浓度高达2.6 mg N·L^-1，这是因为大气中含有大量气体NH₃和颗粒物NH₄⁺。随着降水过程的持续，降水开始后的10~40 min内（13:25—13:55）采集到的分段样品NH₄⁺浓度在1.5 mg N·L^-1上下波动。降水结束前的1 h内（13:55—14:55），分段降水样品NH₄⁺浓度稳定在0.5 mg N·L^-1。降水NO₃^-浓度的时间变化规律与NH₄⁺基本相似。

此次降水过程NH₄⁺和NO₃^-的湿沉降量分别是22.6 mg N·m^-2和11.9 mg N·m^-2，还原性氮与氧化性氮的比值是2:1，与华北地区典型农田氮湿沉降的形态构成相似^[8]。由此，结合上一节的干沉降通量，估算8月10—13日研究期间干、湿沉降无机氮总量为39.6+ 34.5=74 mg N·m^-2（0.74 kg N·hm^-2），这与山东禹城和河北栾城观测到的无机氮沉降的通量范围一致（分别是70.7 mg N·m^-2和75.5 mg N·m^-2）^[8]。

2.4 土壤剖面无机氮分布和NH₃排放量估算

观测期间，土壤剖面NH₄⁺主要集中在表层0~10 cm，NH₄⁺含量随着土壤深度增加呈现出幂指数下降趋势（图 5）。表层土壤NH₄⁺含量平均值为0.03±0.01 mg·g^-1，略高于华北的河北栾城、南皮和山东惠民等地（无机氮为0.006~0.029 mg·g^-1，0~20 cm）^[32]。土壤NH₄⁺含量在0~25 cm深度内由0.41 mg·g^-1剧烈下降至0.03 mg·g^-1。土壤深度25~50 cm范围内，NH₄⁺含量变化幅度较小，维持在0.03 mg·g^-1左右。土壤剖面NO₃^-含量的变化与NH₄⁺略有不同，NO₃^-在0~15 cm深度内略有下降趋势，15~50 cm各层中NO₃^-含量稳定在0.07~ 0.10 mg·g^-1。

假设深层土壤中所含有的物质含量代表母质背景，表层土壤的物质含量代表母质背景上叠加的人为源输入，依据表层和深层物质含量的差异可估算人为源输入量^[33]。由此估算人为源输入到土壤中NH₄⁺-N和NO₃^--N总和约420 kg N·hm^-2·a^-1，与施肥量（430 kg N·hm^-2·a^-1）相当。大气氮沉降是农田土壤和作物可利用氮的重要来源，需要在农田氮管理中充分考虑氮沉降输入量的动态变化。

3 结论

本研究以华北典型农田玉米生态系统为研究对象，综合研究了植被冠层大气、玉米叶片、土壤剖面和分段降水样品中NH_x浓度的动态变化，讨论大气-植被-土壤界面间NH_x的交换通量，得到以下主要结论：

（1）植被冠层的NH₃和NH₄⁺浓度变化范围分别是27.1~36.6 μg·m^-3和2.1~8.0 μg·m^-3，降水对NH₄⁺的清除作用显著高于NH₃。

（2）植株叶片中的NH_x浓度从顶端到底部呈逐渐下降趋势，且日变化特征一定程度上受温度影响，这些变化可能受制于植株NH₃挥发的潜势。

（3）抽穗期植株叶片NH₃补偿点低于植被冠层NH₃浓度，而花粒期NH₃补偿点整体高于植被冠层NH₃浓度，大气与植株叶片NH_x交换过程在不同生长阶段呈现不同方向。

（4）基于观测结果估算大气Nr通过降水、气体和颗粒物输入到农田生态系统的沉降量分别为34.5、36.1 mg N·m^-2和3.4 mg N·m^-2，约合74 mg N·m^-2。