表面质子储库对增强风化碱度输出的约束

《Environmental Science & Technology》：Surface Proton Reservoirs Constrain Alkalinity Export from Enhanced Weathering

　　

向土壤中添加硅酸盐和碳酸盐材料（增强风化，EW）是一种被提出的二氧化碳去除（CDR）策略，具有潜在的农学协同效益，尤其是对于土壤pH管理。然而，矿物溶解产生的碱度与表面质子储库（SPR）相互作用，在输出之前消耗碱度并使量化复杂化。由于持久性CDR需要将碱度传输至海洋储存库，土壤输出构成了EW性能的机制上可定义的上限。在此，研究人员开发了一个框架，将操作性阳离子交换测量划分为静电库和pH依赖性配位表面储库，并将该表征嵌入反应性运移模型中，以显式耦合深度依赖的SPR、阳离子交换、碱度生成和CO_2(g)交换。研究人员将该方法应用于对比硅酸盐（橄榄石）和碳酸盐（方解石）改良剂在代表性水文条件下的情况。利用基于距离的广义敏感性分析（DGSA），研究人员识别了涉及SPR大小、基线土壤pH、土壤湿度和改良剂速率的多变量控制空间，该空间控制碱度输出效率。模拟表明，改良剂衍生的碱度在早期溶解过程中大部分被中和，显著损害发生在浅层（>25 cm）土壤带下方。因此，基于溶解无机碳（DIC）或碱阳离子的代理估算碱度输出表现出系统性和深度依赖性偏差。这些发现强调了SPR是EW碱度输出的首要限制因素，并为CDR量化提出了适当的筛选和监测策略。

增强风化碱度输出的表面质子储库约束机制

**研究背景与问题**

实现全球净零排放目标需要大规模部署能够实现持久（约10³年）大气CO₂减排的二氧化碳去除（CDR）技术。增强风化（EW）作为一种在土壤中添加反应性硅酸盐或碳酸盐矿物的策略，因其理论可扩展性和农学协同效益而备受关注。当应用于农田时，粉碎的硅酸盐或碳酸盐矿物可通过养分供应和pH缓冲改善土壤肥力，同时产生碱度，经海-气交换最终增加溶解无机碳（DIC）储库，实现持久CDR。因此，碱度输出至土壤带之外是衡量EW-CDR性能的关键指标。

然而，越来越多的现场、实验和模拟研究揭示了关键挑战：可观察的矿物溶解常与可测量的碱度输出相解耦。固相分析持续证明原料溶解，但水相监测往往未能检测到阳离子、碳酸盐碱度或DIC的相应增加，有时超过90%的风化衍生阳离子滞留在土壤剖面内。这些差异引发了对当前测量、报告和核查（MRV）方法可能系统性高估浅层净大气CO₂去除量的担忧。一个关键机制是质子反应性表面位点的普遍存在：随着pH升高，矿物和有机质表面的质子反应性官能团去质子化，中和碱度并产生额外负表面电荷。在中等酸性至近中性（pH?5.5）范围内，交换性酸度最小，pH依赖性表面质子缓冲成为碱度归宿的主导控制因素。然而，操作性阳离子交换容量（CEC_op）测量通常将表面质子储库（SPR）与永久电荷静电交换库合并，掩盖了预测动态田间条件下碱度输出所需的机制划分。大多数反应性运移模型（RTMs）将阳离子交换和表面络合作为独立数值模块实施，迫切需要一种将CEC_op指标映射到化学性质不同的缓冲与交换库的方法。因此，本研究旨在开发一个将常见农学CEC_op指标划分为阳离子交换库和pH依赖性表面电荷库的两储库框架，使RTMs能够机制性表征土壤缓冲，并通过系统模拟评估碱度输出效率及其控制因素。

**研究内容与意义**

研究人员开发的框架将操作性CEC划分为永久电荷（外圈）储库（C_perm）和pH依赖性（内圈）表面电荷储库（C_coord），并在反应性运移模型CrunchFlow中实现深度依赖的SPR、阳离子交换、碱度生成和CO_2(g)交换的显式耦合。模拟采用对比性硅酸盐（橄榄石）和碳酸盐（方解石）改良剂，覆盖代表性水文条件，共进行400次模拟。研究识别了SPR作为碱度输出的首要限制因素，揭示了DIC和碱阳离子代理的系统性偏差，并为CDR量化提出了筛选与监测策略。论文发表在《Environmental Science》。

**主要技术方法**

研究的关键技术方法包括：（1）CEC-SPR两储库划分框架，将CEC_op分解为C_perm和C_coord，实现静电交换与质子缓冲的机制分离；（2）多组分反应性运移模拟，使用CrunchFlow耦合水相反应、矿物溶解-沉淀动力学、表面络合和阳离子交换与平流-弥散运移；（3）基于距离的广义敏感性分析（DGSA），识别非线性时变响应的主导控制因素；（4）碱度追踪方法，区分潜在碱度（Alk_pot）和实际碱度（Alk_act），并定义碱度损失分数（L_alk）。样本队列来源为美国农业部Kellogg土壤调查实验室（KSSL）数据库中选取的美国伊利诺伊州南部三个深度分辨农业土壤剖面（pH_1:1 ≈ 5.0–5.7）。

**主要研究结果**

**基线土壤pH及对SPR的影响**

模拟表明，稳态基线土壤pH受水饱和度（S_w）及由此决定的气体扩散系数控制，而非本地矿物动力学和表面位点密度。基线pH不仅反映长期水文生物地球化学耦合，还设定了SPR的初始物种分布和占据状态。橄榄石和方解石改良均在浅层土壤产生显著pH响应，且随深度衰减，反映了SPR在浅层剖面的优先表达。

**耦合土壤pH与碱度对改良剂的响应**

模拟结果显示强烈非线性的碱度累积路径，深度依赖性衰减凸显了SPR的影响：在25至120 cm上部剖面发生显著碱度损害，而Bt层以下损害相对较小。橄榄石风化10年累积浅层碱度输出中位数为4.5 eq m^–2（对应2.0 tCO₂e ha^–1），而Bt层以下输出衰减约一个数量级（中位数0.01 eq m^–2），且并非水相运移限制所致（保守运移时间至120 cm仅约0.4年）。深度积分SPR而非溶解供应或水文驱动因素主导有效碱度输出。

**方解石和橄榄石改良下pH演变**

两种矿物均表现出农业pH调控潜力，几乎所有模拟均产生pH ≥ 6.5的农学响应。方解石轨迹表现出更显著的峰值后pH下降，归因于快速碳酸盐溶解耗尽反应性供应后呼吸驱动酸输入导致的再酸化；橄榄石的较慢动力学则持续供给碱度并缓和后期下降。pH响应受基线pH、浅层总表面位点浓度（[>SOH]_tot）、改良剂速率和S_w共同控制，而非单一因子。

**不同改良剂的碱度效率**

碳酸盐和硅酸盐改良均表现出SPR导致的即时碱度损害，即使浅层也如此。橄榄石低损失簇在25 cm处前4年输出小于80%产生碱度；120 cm以下输出小于10%。浅层碱度损失簇按SPR容量严格分层。SPR是L_alk变异的主导控制，其次为基线pH和S_w。方解石较橄榄石L_alk下降更快，反映快速溶解后呼吸驱动酸化导致SPR再生。

**DIC和阳离子基碱度输出代理评估**

DIC输出与碱度输出对比显示系统性高估，过量分数中位数约0.5（IQR 0.23–0.72），最大偏差出现在近中性pH处，10年DIC过度计入约0.8 tCO₂ ha^–1。单一阳离子代理（Mg²⁺）因阳离子交换作用而表现弱且系统性偏差；碱阳离子总和（ΣBC）虽线性对应较好，但在低碱度输出时相对偏差达±200%。这反映电荷平衡并不确保保守的碱度输出。

**改变的土壤-大气CO_2(g)通量**

橄榄石改良通过提高pH和碱度将DIC形态从CO_2(aq)移向碳酸氢盐和碳酸盐，降低土壤pCO₂并减少扩散通量（10年减少+0.2 tCO₂ ha^–1）。相反，约四分之三方解石改良导致累积CO₂排放增加（中位数–1.4 tCO₂ ha^–1），与农用石灰构成净CO₂源的既有担忧一致。CO₂通量变化与深部碱度输出仅弱相关，显示浅层和深部信号的解耦。

**讨论与结论**

研究识别深度积分SPR、基线土壤pH和平均S_w为碱度传输的主导控制因素。所建立的RTM-DGSA框架将传统CEC_op指标转换为RTMs中的机制表征，但模拟的导出效率代表次级矿物汇存在下的上界。水文状态构成碱度输出的结构性约束，瞬时干湿循环可能诱导阈值行为。对于MRV设计，研究提出三项建议：碱度输出应在SPR衰减深度以下量化；基于输出碱度而非DIC或碱阳离子代理的量化最可辩护；基于土壤CO_2(g)通量的减排声明必须与持久碱度输出相协调。研究结论指出，虽然增强风化提供明确的农学效益，但持久CDR的可信声明需要基于地球化学质量平衡和运移约束的MRV框架，而非浅层代理等价性。