基于稳定同位素自然丰度技术的土壤氧化亚氮产生与排放过程研究进展

黄 瑾，余龙飞，李文娟，黄 平

（1. 重庆交通大学 建筑与城市规划学院，重庆 400074；2. 中国科学院 重庆绿色智能技术研究院，重庆400714；3. 中国科学院大学 重庆学院，重庆 400020；4. 暨南大学 地下水与地球科学研究院，广东 广州510632）

氧化亚氮(N2O)作为增温潜势极大的温室气体，对全球变化有着重要的影响，探究其源、汇关系意义十分重大。然而土壤N2O的源、汇具有很强的时空异质性，导致难以制定精准的N2O减排策略[1]。总的来看，土壤中N2O产生和排放过程复杂多样。目前，大多数研究认为羟胺氧化、硝化细菌反硝化和反硝化过程是N2O产生的三大主要途径[2−3]，而关于N2O还原过程的认识也还存在较大争论，其中反硝化被认为是N2O还原的主导途径[3]。为了量化上述N2O过程，有针对性地提出N2O减排策略，同位素分析技术在N2O溯源方面提供了重要支撑。同位素自然丰度法中的同位素异位体法，是一种非侵入性方法，因其不受底物同位素影响[4−5]，N2O同位素特征值δ15Nsp可被用作指示多种微生物产生N2O的作用过程[6−7]。N2O的其他同位素特征值，δ15Nbulk和δ18O也被用于指示N2O产生的微生物作用途径，但会受到N2O 前体物 [铵离子 (NH4+)，硝酸根离子 (NO3−)和水 (H2O)等 ]的同位素组成的影响[6, 8−9]。已有很多学者探索利用同位素特征值来准确分析N2O产生的微生物过程，综合研究δ15Nbulk、δ18O、δ15Nsp对揭示N2O产生机制更有重要的意义。而在N2O产生与排放过程中，不同功能微生物过程同位素分馏效应的差异构成了稳定同位素自然丰度技术分析微生物过程的基础[10]，因此，同位素分馏是稳定同位素自然丰度技术应用的理论基础，在研究中的作用不可忽视。本研究梳理了土壤N2O产生和排放过程中氮同位素分馏的效应；阐述了环境因子及微生物对土壤N2O产生和排放过程的同位素分馏效应影响；总结了稳定同位素自然丰度技术在土壤N2O源解析中的应用及进展。

1 土壤N2O产生和排放过程中的氮稳定同位素分馏

微生物参与的N2O转化过程可分为4个经典过程，分别为硝化作用 (自养/异养硝化作用和硝化细菌反硝化作用)、反硝化作用 (真菌反硝化作用和硝态氮异化还原成铵)、硝化反硝化耦合作用和共反硝化作用。借助分子技术已基本确定控制硝化作用微生物过程的酶(编码基因)，但仍需结合稳定同位素技术来分析具体的代谢过程[11]。氮同位素分馏作用是指随着生物化学反应的进行，反应底物和产物δ15N丰度发生差异的现象，这是15N稳定同位素自然丰度技术应用的理论基础[12]。本研究针对微生物驱动的硝化作用和反硝化作用的氮同位素分馏效应进行了梳理。

1.1 硝化作用

土壤硝化作用(nitrification)可分为自养硝化和异养硝化[13]。当前，自养硝化过程中的稳定同位素分馏得到了广泛的研究。自养硝化作用是指NH4+或氨气(NH3)经亚硝酸根离子(NO2−)氧化为NO3−的过程，包括羟胺氧化和硝化细菌反硝化2个部分。通过纯培养实验进行同位素检测，获得NH3或NH4+向羟胺(NH2OH)转化以及NH2OH向NO2−转化2个过程的同位素效应平均值，约(−29.6±4.9)‰。这2个过程主要涉及亚硝化单胞菌，其中欧洲亚硝化单胞菌Nitrosomonas europaea的氮同位素效应为 −38.2‰～−24.0‰；而其他亚硝化单胞菌，包括亚硝化单胞菌N. eutropha、N. marina C-113a、Nitrosospira tenuis，其氮同位素效应范围更大，为−38.2‰～−14.2‰[14−15]。仅 CASCIOTTI[16]对 NO2−向 NO3−的转化过程做了同位素效应研究，同位素效应均值为(13.0±1.5)‰。在自养硝化过程中，N2O主要通过羟胺或NO2−的化学分解产生[17]，同位素效应值也会随反应底物的不同而发生变化。NH2OH到N2O过程在纯培养实验中的同位素效应值为−26.3‰～5.7‰[18]，但后来有研究发现：将NO2−和NH2OH以1∶1混合作为反应底物，其同位素效应值变弱，为–17.8‰～0.8‰，比NH4+到N2O过程的表现得更弱，表明这些氧化过程会使同位素效应值减弱。仅以NH2OH作为底物来估算反应的同位素效应值，得出的结果为(−5.1±12.0)‰。

硝化细菌反硝化(nitrifier denitrification)是一种由氨氧化细菌(AOB)驱动的，但又有别于自养的氮转化过程[10]，该过程中，NH4+或NH3氧化产生的NO2−被同一类氨氧化细菌进一步还原为N2O和氮气(N2)[10]。该步骤涉及的酶和编码基因与硝化作用和反硝化作用中参与相应步骤的酶和基因相同，因此该过程的同位素效应值可参考硝化作用和反硝化作用的相应过程。目前，关于硝化细菌反硝化过程能否还原N2O形成N2还存在争议。有研究认为：AOB中并未发现编码N2O还原酶同源物质的基因[19]，这表明硝化细菌反硝化过程可能不是N2O的还原过程；然而，BEYER等[20]认为：尽管传统AOB中没有发现N2O还原酶，但它存在类似N2O还原酶的电子铜蛋白亚硝基(nitrosocyanin)，这可能与N2O还原成N2有关，但其具体作用机制也还有待研究证实。

1.2 反硝化作用

土壤反硝化作用(denitrification)一般是指反硝化细菌在厌氧或微氧条件下将NO3−、NO2−或者N2O作为呼吸过程的末端电子受体，并将其还原为NO2−、NO、N2O或者是N2。N2O是反硝化过程的中间产物[4]，甚至某些条件下是反硝化的最终产物[21−22]。反硝化过程的N2O还原成N2是N2O消耗的关键过程，在调节自然生态系统N2O排放中扮演着重要角色[23]。在土壤培养实验中，测得NO3−向NO2−转化过程中的同位素效应值为−52.8‰～−10.0‰，平均值约 (−31.4±11.8)‰[24]。NO2−向 NO 转化的同位素效应较弱，并且该过程以细菌和真菌参与的反硝化作用为主，同位素效应平均值约(−19.8±7.6)‰。NO向N2O转化的同位素效应值为(14.0±1.6)‰[25]。OSTROM等[26]纯培养实验的研究结果表明：N2O还原为N2过程的同位素效应平均值为(−4.1±6.8)‰，但该过程同时存在正反2种同位素效应。在不考虑反向同位素效应时，其同位素值为−9.1‰～−2.5‰，平均值为 (−6.6±2.7)‰[18, 27]。

除细菌外，某些真菌(如尖孢镰刀菌Fusarium oxysporum)也可以通过反硝化过程产生N2O，称为真菌反硝化[28−29]。但真菌一般缺乏编码N2O还原酶的基因，因而N2O是真菌反硝化的最终产物[28]。真菌反硝化作用对N2O排放总量的贡献尚未完全阐明[30]，此过程的同位素分馏效应也有待进一步研究。

硝态氮异化还原成铵态氮的过程(DNRA)是指在厌氧环境、高pH以及大量的易氧化态有机物存在下，微生物以NO3−或NO2−作为电子受体，还原成可利用性NH4+，将活性氨滞留在生态系统中的过程。NO3−还原成 NO2−后，NO2−主要还原成 NH4+并且伴随 NO2−的短暂积累和 N2O的产生[31]。其中 N2O 是该过程中NO2−还原的副产物[32−33]。实际上，自然环境中DNRA过程对N2O排放的贡献研究仍然较少。有一些研究认为：DNRA过程在N2O产生和排放中所占的比例较小，甚至可被忽略[34]。NO3−还原为 NO2−过程的同位素效应值也为(–31.4±11.8)‰[24]。但目前亚硝酸盐还原为铵盐这一关键反应的同位素效应值尚不清楚。此外，一些DNRA过程细菌，如Wolinella succinogenes和Anaeromyxo bacterdehalogenans也具有编码N2O还原酶的基因，这可能对N2O还原具有潜在贡献，但其确切过程还需进一步证实[33,35]。

2 土壤N2O产生和排放过程中与同位素分馏效应有关的主要因子

土壤产生N2O是复杂的微生物过程，研究多个影响因子协同作用对N2O排放量变化的影响，有助于了解N2O排放的内在机制[5, 36]。在N2O转化过程中，不同微生物同位素分馏的差异构成了同位素手段分析微生物过程的基础[10]。在纯培养实验中(如真菌纯培养实验)，同位素效应的可变性可能是因不同的微生物菌株或多种微生物导致，也可以由特定微生物的同位素效应来解释，例如，脱氮副球菌的值在–28.6‰～–10.0‰之间变化很大[18]。不同菌株之间的部分变异性可能来自不同的酶结构[15]，但单个微生物体测定值出现的巨大变化表明特定的实验条件(如扩散底物的可用性、中间产物的积累等)对特定的同位素效应起着重要作用[18]。

同位素分馏效应主要受到反应平衡条件以及反应速率的影响(图1)。而环境因素的改变(如土壤温度、水分等)会影响反应速率，并进一步对微生物的生理过程、群落组成和反应活性产生影响，进而影响N2O的产生和排放过程及同位素分馏效应。关于实验的环境条件(如温度、湿度、底物浓度等)与同位素效应之间的关系，有少部分学者开展了相应的研究。MARIOTTI等[37−38]研究了同位素效应对温度的依赖性，他们在土壤培养实验中观测到30 ℃下同位素分馏效应(−24.6‰)明显弱于10 ℃下的值(−31.2‰)。较高反应温度下反应速率较快，导致生成物转化效率高并使得同位素分馏效应较小[39]；而在低温下的反应进程缓慢，通过扩散传递的底物相对较高，因此表现出更强的酶同位素效应。MATHIEU等[40]还通过在底物中添加NO3−来研究同位素组成的影响。随着底物消耗量的增加，培养实验中也观察到了同位素效应的降低[41−42]。当同位素效应转移到扩散的较弱同位素效应时，会导致同位素效应减少[37]。然而，在不同湿度和温度条件下使用山地和湿地进行的土壤培养实验结果相互矛盾[39]，反应速率与同位素效应差异之间没有显著关系。其中，各处理间的反硝化速率相差1个数量级以上，但测定的同位素效应相差不到9‰。以上研究说明：同位素分馏效应(尤其是野外环境)的评估依然具有较大不确定性，除反应速率之外可能还存在其他制约因子，比如功能微生物群落构成、微生物活性等。

图1 微生物驱动的土壤N2O产生和排放过程中的同位素分馏效应Figure 1 Isotope fractionation effects in N2O production process driven by microbe

3 稳定同位素自然丰度技术在土壤N2O溯源中的应用及进展

稳定同位素自然丰度技术是指利用N2O的15N和18O的自然丰度变异来研究N2O的产生、排放过程[43]。N2O自然丰度作为N2O源和汇的一种更综合的示踪剂，已广泛用于土壤N2O源解析[44−45](图2)。然而，由于δ15Nbulk和δ18O往往依赖于底物的同位素组成，且各个过程中同位素分馏效应存在不确定性，很难通过δ15Nbulk和δ18O的变异准确解释N2O排放来源[5, 45]。

图2 稳定同位素技术在 N2O溯源中的发展Figure 2 Development of stable isotope techniques in N2O source partioning

同位素异位体法是在稳定同位素自然丰度技术上发展起来的，但不同于稳定同位素自然丰度技术，同位素异位体法(δ15NSP)[46]并不受到底物同位素丰度影响，在N2O扩散过程中相对稳定[47−48]，其仅受N2O形成的机制和途径[49]以及N2O还原过程的影响。此项优点为高精度N2O溯源提供基础，且不干扰系统、适用于田间试验，此外还提供了估算DNRA过程产生N2O的潜力[5]。N2O分子为不对称的直线型结构，其分子式为N=N=O，中间和边缘不同的氮原子被分别称为α原子和β原子[50]。N2O的各产生途径因受不同的微生物或酶调控，造成N−O键的断裂位置不同，故生成的N2O的SP值也有所差异，这被认为是运用同位素异位体法溯源N2O产生途径的理论基础[45]。同位素异位体法N2O分子中SP值可以判别两部分N2O产生过程，一部分为硝化细菌硝化作用和真菌反硝化作用产生的N2O过程(N2ON)，另一部分为硝化细菌反硝化作用和反硝化作用产生的N2O过程(N2OD)。同位素异位体法已经开始被用来区分土壤中N2O各产生过程的微生物来源[7−9]，表现出良好的应用前景。

然而，仅使用δ15NSP进行源解析也会面临一些问题：①N2O还原过程中SP值的增加，即N2O还原过程中的同位素分馏效应。②真菌反硝化作用和硝化作用的SP值重叠。③不同培养物和反应过程中SP的变异性。此外，正如MOHN等[51]所证明的，由于不同的分析技术和缺乏可用的N2O同位素标准物质，可能导致不同实验室间SP值测定结果有较大差异。虽然国际上已发布2个N2O同位素标准样品USGS51和USGS52[52]，但有关野外同位素值的测定仍有待确定。为了解决上述问题，很多学者在不断探索。问题①可以在模型中解决，因为有研究表明：在N2O还原为N2的情况下，SP的同位素效应稳定范围为−8‰～−2‰[47, 53−54]。但问题②尚未解决。有研究[25, 55−56]表明：SP在反应过程中是稳定的，但鉴于P450nor真菌一氧化氮还原酶还原氧化氮(NO)过程中SP从15‰增加到29‰，SP值的稳定性可能需要进一步证实，这也是问题③有待解答的问题，不同培养物中不同的酶和N2O之间存在的特定偶联机制。

为了更准确地测定N2O还原程度和微生物贡献，实际应用中通常采用双同位素图谱法[28, 57−58]，即使用δ18O/δ15Nbulk、δ15NSP/δ15Nbulk或δ15NSP/δ18O，即用2个同位素特征值来评估N2O不同产生途径的贡献率和N2O的还原程度。WRAGE等[17]使用了δ18O/δ15Nbulk进行硝化、硝化细菌反硝化及反硝化过程产生N2O的来源区分。但在实验环境中水的18O与土壤中O2以及NO3−存在发生氧交换的风险，导致结果分析不准确。H2O中氧原子会与NO2−和NO3−中氧原子发生氧交换，且它们对N2O中氧原子的贡献与微生物种类有关。因此，需要考虑自然条件下氧交换对δ18O值的影响。KOOL等[59]使用富集比保留值法(ERR)对反硝化过程中的氧原子交换进行了测定，尽管有待改进，但这种方法能够区分硝化、硝化细菌反硝化、硝化反硝化耦合作用等N2O生产过程，并应用于一系列土壤。TOYODA等[57]首次将δ15NSP/δ15Nbulk应用于农业土壤。之后，有研究[59−62]利用这种关系来确定N2O混合比例和N2O还原程度，发现δ18O也可以很好地示踪N2O的产生过程，这得益于细菌反硝化期间高的氧原子(O)交换量，从而使该途径的δ18O值非常稳定。基于这一发现，LEWICKA-SZCZEBAK等[63]提出了解析N2O源过程的δ15NSP/δ18O图，可用于校正反硝化作用中N2O还原的影响和定量微生物作用贡献，最近已经广为应用[64−66]。δ15NSP/δ15Nbulk和δ15NSP/δ18O已联合应用于现场研究[65]，并在计算N2O还原程度和真菌、细菌反硝化值之间显示出较好的一致性。但利用双同位素图谱法在精确量化N2O产生的路径比例和N2O还原程度方面仍存在许多限制，比如单一产生路径的同位素特征值范围广、不同产生路径的同位素特征值范围的重叠、反应底物同位素组成的变化，以及与N2O还原相关的分馏因子的可变性等。因此基于双同位素图谱法所提供的定量结果的准确性仍有待探究。

最近，LEWICKA-SZCZEBAK等[67]在双同位素图谱法的基础上，把3个同位素特征值(δ15Nbulk、δ18O、δ15NSP)结合起来计算N2O不同产生途径的贡献率和N2O的还原程度。提出了一种基于三维同位素的复杂模型，并对其结果进行了验证。这是首次在实验室和田间条件下对N2O同位素自然丰度法溯源N2O进行系统地验证研究。然而将N2O同位素特征值完全整合到模型中，需要建立一个包含同位素效应及其不确定性的综合数据库[18]。

4 总结与展望

未来关于稳定同位素技术在土壤N2O溯源中的研究需注意以下几个方面：①同位素自然丰度法在N2O溯源应用中还有很大的提升空间。ZHANG等[68]指出：需要进一步探究同位素特征值稳定性来提高对 N2O微生物过程的识别。把同位素特征值 (δ15Nbulk、δ15NSP、δ18O)与底物同位素 (如 δ15、δ15、δ18、δ18)相结合，分析反应底物在N2O生成过程中的分馏效应可以提高源解析的准确性[69]。把N2O的δ15Nbulk、δ18O、δ15NSP同位素特征值结合起来有益于识别未知源[70]，但在不同的pH、氧化还原状况、基质以及非生物N2O产生途径的条件下测得不同的δ15NSP，说明同位素效应及其变异性是制约精确度的关键，其中微生物过程及其导致的同位素分馏效应是N2O溯源的重要切入点。此外，目前对于18O富集系数ε18O的研究还相对较少，且一些研究认为：不同过程产生N2O的δ18O无明显区别[39]，这表明对δ18O的特征值研究有助于准确解析N2O的来源。并且未来还需对N2O产生过程中每一步的δ15NSP都进行测定，探究环境因素与产生N2O具体的生化过程速率之间的关系(包括共反硝化)，在排除共反硝化的影响下，对真菌反硝化的δ15NSP进行改善[56, 71]。②除了N2O转化过程中同位素分馏效应会导致误差外，N2O同位素的分析和测定也需要标定和校准[72−73]。OSTROM等[52]提出2种N2O气体标准物质USGS51和USGS52用于实验室标定和校准。未来仍需要多个实验室间进行更多平行比较并制定统一的校准方案，如对团簇同位素的精确测定有助于质谱值重叠样品的校准。③新方法的探索。稳定同位素技术与分子技术结合[74]。很多研究都提出应该把δ15N与qPCR、基因拷贝数的测定等分子技术相结合来对N2O生物过程进行区分，并且可以结合微生物生态学分析和量化N2O∶N2的产生和不同微生物N2O来源之间的分配[75−79]。另外，团簇同位素特征值的方法也开始流行。SP值不能用于区分硝化和真菌反硝化产生N2O的过程，但是可以使用团簇同位素特征值的方法，通过具体产生途径中特定的酶来判断[80]，但如何准确判断特定的酶仍需要进一步的探索。