管理重塑全球资源梯度下的地下生产力分配

《Ecosystems》：Management Reshapes Belowground Productivity Allocation Across Global Resource Gradients

　　

地下碳分配调控土壤碳输入与生态系统韧性，然而关于管理如何沿资源梯度改变分配方式的全球证据仍然有限。研究人员利用全球陆地净初级生产档案，检验人类管理是否改变总净初级生产向地下分配的占比，以及该效应是否与气候和土壤因子对总净初级生产的控制相当。研究人员将456个陆地站点（涵盖森林、草地、泥炭地、灌丛、苔原和耕地）的重复野外测量聚合为站点均值，为293个站点计算地下生产占比，并拟合纳入存档方法学不确定度的加权最小二乘模型。非管理草地比管理草地分配了显著更多的生产到地下（中位数差异=0.182，95%置信区间0.097–0.238，p<0.001），而森林未检测到管理差异（中位数差异=0.003，95%置信区间?0.057至0.059，p=0.967）。在最终分配模型中，非管理×开放生态系统交互项为正（估计值=0.154，95%置信区间0.054–0.254，p=0.003），交叉验证R2从气候基线的0.160升至0.243。对于总净初级生产，气候与资源协变量将交叉验证R2从0.223提升至0.269，年均温、干旱指数（AI）和大气氮沉降呈正向效应，而在排除耕地与处理实验的保守子集中管理项减弱。结果表明，管理比总生产力更强地重组碳分配，尤其在开放生态系统中，为生态系统科学与模型评估提供了基于野外的基准。

研究背景与问题提出

陆地生态系统通过植物生产决定生物量积累与土壤有机质输入，介导着显著的碳–气候反馈。全球陆地碳汇在1960年代至2010年代近乎翻倍，其中地下净初级生产（Belowground Net Primary Production，BNPP）约占陆地碳固定的46%，是土壤有机碳与养分循环的核心通路。然而地下生产占比（f_B=BNPP/NPP_tot）的野外观测稀缺、方法异质，且管理（放牧、收获、施肥、物种组成改变）如何沿气候与资源梯度重组碳分配仍缺乏多生物群系（biome）综合证据。既有大尺度研究多聚焦总生产或地上绿度，未将管理效应与站点方法学不确定度联合传播。研究人员基于Rodal等人2025年发布的全球陆地NPP数据库，旨在检验“管理对NPP分配的重组作用是否强于对总NPP量级的作用”这一核心假设。该文发表于《Ecosystems》。

主要技术方法概括

研究人员采用Rodal等人汇编的456个陆地站点档案（森林202、草地149、泥炭地、灌丛、苔原、耕地合计105），将各站点重复测量按算术均值折叠为站点均值；地下占比由BNPP_s/NPP_tot,s计算，总NPP取自然对数。协变量来自WorldClim 2、CRU TS、Ackerman氮沉降栅格及Zomer全球干旱指数（AI=P/PET₀，FAO-56 Penman–Monteith），土壤变量按HWSD v2.0单元份额加权聚合。质量权重w_sk由观测时长L_sk与地上/下方法学折减因子（RF^AG、RF^BG）均值平方反比构成；预测变量标准化后，分配模型拟合加权最小二乘（式10：f_B~AI^*+非管理U+开放生态系统O+U×O），总NPP模型拟合式11（含年均温T^*、AI^*、氮沉降N^*、有效含水量AWC^*、pH^*、U与biome哑变量），采用HC3异方差稳健协方差，以10折交叉验证（种子42）计算R²_CV、RMSE、MAE，并对草地/森林分别做Mann–Whitney与10000次bootstrap中位数差异检验，另设去耕地、去处理实验等敏感性子集。

研究结果

研究区域与数据可用性

456站中455站有管理标签（152管理、303非管理），293站可算f_B，344站有总NPP；最终分配模型n=289，总NPP模型n=325。耕地全为管理态，仅8个耕地站同时具备地下与总生产。

地下分配模型结果

草地非管理与管理站点f_B bootstrap中位数差异0.182（95%CI 0.097–0.238，p<0.001，n=134），森林差异0.003（95%CI ?0.057–0.059，p=0.967，n=102）。嵌套模型R²_CV由气候基线0.160升至加管理0.210、最终交互模型0.243；最终模型中U×O估计0.154（95%CI 0.054–0.254，p=0.003），标准化AI系数?0.034（p=0.037），主效应U与O不显著。

分配验证与敏感性

去耕地后U×O=0.125（p=0.021），去处理实验后0.221（p<0.001），去两者后0.161（p<0.001）；草地观测均值差0.151，森林0.011。交互项方向在所有子集一致为正。

总生产力模型结果

总NPP模型调整R² 0.591，R²_CV 0.269；年均温系数0.277（p<0.001）、AI 0.207（p<0.001）、氮沉降0.145（p=0.026）、pH ?0.171（p=0.047），AWC不显著；以管理耕地为参照，草地/灌丛/苔原拟合对数NPP更低，非管理项?0.208（p=0.043）。

总NPP验证与敏感性

保守子集（去耕地+处理实验）中温度0.335、AI 0.165仍显著，氮沉降降至0.089（p=0.162），非管理项?0.166（p=0.073），调整R² 0.494，表明管理项依赖人为干预最强站点。

讨论部分总结

研究人员指出三条主线：管理信号集中于分配而非总量；该效应具生态系统结构依赖性——开放生态系统（尤其草地）驱动对比，森林因木质结构需求与冠层缓冲未现可检测差异；气候与氮沉降对总NPP的解释力稳定强于二值管理。这与Gherardi & Sala干旱主控下f_B下降、Li等人草地上下部资源限制分化、Chomel等人集约草地削弱干旱后光合产物向根–菌根转移的实验证据一致。机制上，开放系统根系直接竞争表层水养、管理偏好快速地上再生；森林 woody 结构需求掩盖细微分配变化。局限性包括二值管理不能区分放牧强度/施肥/弃荒时长、站点均值压缩时间变异、地下测量残余方法异质。意义在于为陆地生物圈模型（TBM）评估提供“分配基准”：若模型仅校总通量而忽略f_B对管理的响应，将误估土壤碳输入与地下韧性。

结论部分翻译

本研究利用全新全球野外数据库检验管理是否改变陆地生态系统在上下部 compartment 间分配NPP。三点结论：其一，管理信号在NPP地下占比中最强而非总量，非管理草地比管理草地显著更多分配至地下；其二，关键分配结果依赖生态系统结构——非管理×开放生态系统交互项在保守子集中仍为正，森林观测f_B无管理差异；其三，总NPP与温度、长期水分可用性与氮沉降的关联比与二值管理的关联更一致。研究主要贡献非高精度预测模型，而是基于野外证明人类管理可在全球尺度重组碳分配路径，即便对总生产力效应较弱或依情境而定。该区分对评估至关重要，因地下分配支配根系生产、土壤碳输入及模型对陆地碳循环的表征；未来工作应解析何种管理强度与干扰史在特定biome中最强偏移分配。