基于田间高分辨率地球化学监测与数值模拟的农业场地氮转化同位素效应定量研究

《ACS ES&T Water》：Field-Based Quantification of Isotope Effects during Nitrogen Transformations Using High-Resolution Geochemical Monitoring and Numerical Modeling at an Agricultural Site

　　

农业土壤中的氮（N）转化强烈影响硝酸盐（NO3–）淋滤、温室气体排放和水质。然而，由于微生物反应与水文传输在变化的氧化还原电位、温度和pH条件下相互作用，量化这些过程的田间条件仍然具有挑战性。稳定同位素是氮循环的强大示踪剂，但控制这些过程的同位素效应通常源自实验室研究或从聚合的田间观察中推断，两者均未完全捕捉自然地下环境的复杂性。本研究提出了一种基于田间的框架，整合高分辨率地球化学监测和基于物理的数值模拟，以约束自然条件 下氮转化的同位素效应。结果表明，反硝化同位素效应系统地受氧化还原电位、温度和pH影响。在强还原条件下，反硝化快速进行，效应较小；而在更氧化的条件下，效应较大。本研究中得出的同位素效应通常低于典型的实验室估计值，反映了地下环境异质性的影响。相比之下，与硝化相关的同位素效应无法得到稳健约束，这凸显了对额外同位素和过程水平测量的需求。总体而言，该方法通过提供更能反映自然变异性的基于田间的同位素效应估计，有助于改进对非均质地下环境中硝酸盐衰减的评估。

论文解读

研究背景与问题提出

农业景观中的氮（N）循环是导致水质退化的关键途径，主要源于施肥土壤中的硝酸盐（NO₃^–）淋滤进入地下水和地表水。过量的氮输入超过植物吸收能力，引发富营养化、一氧化二氮（N₂O）排放及饮用水资源污染。理解控制硝酸盐产生、传输和转化的过程对于缓解污染至关重要。稳定同位素技术，特别是分析硝酸盐中的δ¹⁵N和δ¹⁸O，被广泛用于研究氮生物地球化学循环，例如区分硝酸盐来源和约束硝化与反硝化等转化过程。这些同位素信号由同位素分馏塑造，导致剩余底物中重同位素逐渐富集，通常用同位素效应（ε）表示，即轻、重同位素体反应速率的相对差异。

尽管近期的模型进展提升了概念理解，但同位素效应通常仍源自实验室实验，其中微生物动力学在受控条件下被检测。然而，此类实验设置可能无法完全代表田间环境中发生的复杂相互作用，因为水文动力变异性、微生物群落动态和土壤异质性会影响反应速率和同位素效应。因此，从实验室条件推断的同位素效应可能与田间观察到的不同，反映了动力学和环境影响的共同作用。此外，基于田间的富集因子通常是从空间或时间聚合的观察中估算的，这可能限制其过程水平的解释。这凸显了需要更严格确定的同位素效应来告知和评估氮循环的概念与数值模型。环境条件还可以通过改变微生物活性和反应路径间接影响同位素分馏。温度、pH和氧化还原电位调节反硝化动力学，从而调节ε的大小。由于反硝化通常需要还原条件，氧化还原波动可以通过改变反硝化微生物群落的活性和路径来改变反应速率和同位素分馏的表达。因为这些控制因素与开放系统中的水文传输和混合共同发生，从受控设置推断的ε值可能无法直接转化为田间条件，这促使研究人员寻求在耦合的传输-反应框架内约束同位素效应的方法。

本研究旨在通过直接从原位田间测量中推导同位素效应来解决这一差距，而不是仅仅依赖受控的基于实验室的估计或聚合观察。通过将17个月的高分辨率硝酸盐浓度和同位素组成深度剖面与HYDRUS中的数值模拟相结合，研究人员在自然条件下的包气带和浅层地下水中估计了同位素分馏。研究人员的方法在基于物理的传输-反应框架内推导了内在的同位素效应。从概念上讲，该模型通过明确包含平流和弥散，将先前用于氮转化的同位素公式扩展到开放系统条件，从而将反应驱动的同位素效应与异质农业环境中的溶质传输动力学联系起来。

主要技术方法

研究人员在德国黑森州里德农业区的一个场地进行了为期17个月（2022年12月至2024年4月）的监测，每月采集水样和土壤样。水样分析了主要离子、水稳定同位素（δ²H-H₂O和δ¹⁸O-H₂O）以及硝酸盐稳定同位素（δ¹⁵N-NO₃^–和δ¹⁸O-NO₃^–）。土壤样测定了氮含量和δ¹⁵N_bulk组成。同时，研究人员开发了一个数值模型，使用HYDRUS软件模拟4米深剖面的水流、热传输和根系吸水。模型校准采用PEST与CMA-ES优化器。随后应用粒子追踪计算轨迹，模拟平流传输，并通过引入随机游走分量来处理弥散效应。在质量平衡方面，研究人员使用同位素质量平衡校正土壤氮库的同位素组成，扣除溶解氮的贡献。同位素效应计算则基于闭系统模型，通过分析单个粒子轨迹来分离反应动力学与混合过程，并根据氧化还原电位（Eh）确定主导过程是硝化还是反硝化，进而分别计算相应的同位素效应。

研究结果

同位素和水文地球化学组成

土壤氮库沿剖面呈现明显的分层。在0.85米至2.75米深度之间，δ¹⁵N-NO₃^–和溶解氮浓度保持相对稳定，表明净硝酸盐变化有限，可能反映了较低的就地氮周转或硝酸盐供应/传输与同时发生的生产和消耗之间的平衡。此区间氧化还原电位（Eh）值表明主要为氧化条件，不利于净反硝化。pH接近中性至微碱性。随着深度增加，沉积物变得更致密湿润，δ¹⁵N-NO₃^–和δ¹⁸O-NO₃^–开始急剧增加，而溶解氮浓度下降，表明发生反硝化作用。这个反应区与氧化还原条件的显著转变相吻合，Eh从上部的氧化值下降到明显更还原的值。pH在此区间略有系统性增加，这与质子消耗的异养反硝化一致。在更大深度（>3.55米），硝酸盐的同位素组成表明来自不同池的硝酸盐发生了混合，这些池受到不同程度的反硝化影响。很可能部分硝酸盐源自侧向流入。这种混合进一步得到了该区域硝酸盐浓度的轻微增加以及Eh并未随深度变得更还原的支持。整体而言，同位素和浓度数据的垂直分布表明，反硝化作用在中等深度最为显著，而最深区域（>3.7米）主要受混合过程影响，而非持续的同位素富集。

基于模型的同位素效应推导

硝化作用

研究人员探索了模型参数空间的广泛部分，以调查整个合理参数行为范围。总共评估了约3.63亿个参数组合。其中约64%满足模型的接受标准（δ¹⁸O值的残差小于0.5‰）并被保留用于进一步分析，表明模型具有相当大的灵活性，许多不同的参数集产生合理的结果（非唯一性）。对接受的参数集的分析揭示了模型校准内的明显约束模式。五个参数中的三个（¹⁸ε_AMO-O₂、¹⁸ε_AMO-H₂O和¹⁸f_eq）在其各自的先验范围内表现出广泛的、近乎均匀的后验分布，表明它们仅受模型结构和应用选择标准的弱约束。参数α_AMO-O₂和α_AMO-H₂O描述了氨氧化过程中的同位素效应，在测试范围内对模型输出的影响似乎有限。它们的弱约束可能反映了可用于解析硝化途径此步骤的同位素信息有限，也可能指向（结构性）不可识别性，因为这两个项可以相互补偿，因此在同位素模型应用中通常被视为复合的氨氧化同位素效应。参数¹⁸f_eq也显示出均匀的后验分布，并且似乎与这些同位素效应相互作用，可能补偿其影响，并加强了这组参数内的等效性。这种行为可能源于缺乏亚硝酸盐同位素（δ¹⁸O-NO₂^–）测量，否则将为此阶段硝化过程中的同位素分馏和交换过程提供约束。相比之下，其余两个参数（¹⁸ε_NXR-NO₂^–、¹⁸ε_NXR-H₂O）始终显示明显的偏斜后验分布，接受值始终偏向其各自先验范围的较低值，表明更强的结构敏感性或更严格的观测约束。参数¹⁸ε_NXR-NO₂^–收敛于-6.6‰左右，而¹⁸ε_NXR-H₂O集中在12.5‰附近。这两个值与Nitrobacter sp. Nb 355的实验结果非常匹配，后者报告的¹⁸ε_NXR-NO₂^–和¹⁸ε_NXR-H₂O值分别为-6.5‰和12.3‰。类似的幅度也在海洋亚硝酸盐氧化细菌Nitrococcus mobiliz和厌氧氨氧化（ANAMMOX） biomass孵化实验中观察到，尽管这些系统在微生物组成和环境背景上有所不同。模型与实验同位素效应之间的相似性表明，该模型可以约束亚硝酸盐氧化过程中的氧同位素分馏。然而，推断的参数范围应解释为在监测活动捕获的田间条件范围内的有效值。环境变量如pH可以调节表达的同位素效应和¹⁸f_eq。由于pH在硝化模型方程系统中未表示为显式驱动因素，此类效应无法直接解析，而是会被纳入推断的有效参数范围中。因此，环境条件的异质性可能会扩大可接受参数组合的集合，并降低参数的可识别性。

反硝化作用

由于3.55米以下的同位素组成可能受到与受反硝化影响较小的富硝酸盐水的混合影响，反硝化作用的同位素效应仅针对此深度以上采集的样本进行计算，从而确保推导的富集因子主要反映反硝化驱动的过程。15ε_NAR和¹⁸ε_NAR均呈现右偏分布，大多数值低于20‰，少数高达100‰的高异常值。¹⁵ε_NAR的平均值（13.2‰）和中位数（7.8‰），¹⁸ε_NAR的平均值（11.3‰）和中位数（5.8‰）表明存在这种偏斜。总体而言，结果表明ε值较低，与先前的研究一致，显示自然条件下表达的同位素效应通常小于受控实验室研究中测量的值。相反，¹⁸ε_NAR/¹⁵ε_NAR的比率显示出更对称的分布，中心接近1，平均值为0.95，中位数为0.89，表明硝酸盐还原过程中氧和氮同位素分馏之间存在广泛一致的1:1关系。四分位距（0.63-1.12）非常适合通常预期的比率范围（约0.5-1.0），反映了硝酸盐还原过程中N和O同位素分馏的特征耦合。它也符合先前研究报告中的更广泛数据集（0.07-2.19，中位数=0.90），强调了研究人员的结果与发表值的总体