一种用于在国家层面监测碳通量的土壤传感机制

时间:2026年5月16日
来源:Soil Science Society of America Journal

编辑推荐:

摘要 土壤是地球上最大的碳储存库,根据管理方式的不同,它可以成为大气中二氧化碳(CO2)的来源或汇。土壤表面释放的CO2(FCO2)与土壤的生物和物理化学性质直接相关。我们的目标是利用可见光(400−700 nm)、近红外(700−1100 nm)、短波红外(1100−2500

广告
   X   

摘要

土壤是地球上最大的碳储存库,根据管理方式的不同,它可以成为大气中二氧化碳(CO2)的来源或汇。土壤表面释放的CO2(FCO2)与土壤的生物和物理化学性质直接相关。我们的目标是利用可见光(400−700 nm)、近红外(700−1100 nm)、短波红外(1100−2500 nm)和中红外(2500−25,000 nm,4000−400 cm−1)反射光谱技术、数字土壤制图以及机器学习方法,来估算并绘制巴西境内的净生态系统生产力(NEP)空间分布图。我们创建了以土壤物理、化学和微生物性质为协变量的FCO2和碳封存潜力预测模型,而空间分布的确定则基于裸土图像、地形、气候和土壤矿物学信息。对于FCO2,我们拟合了一个多元回归模型,其R²值为0.35;对于碳封存,我们拟合了一个空间误差模型,其R²值为0.76。空间分布的准确率在0.41到0.76之间,在外部验证中的相关性为0.56。NEP地图显示,塞拉多(Cerrado)、大西洋森林(Mata Atlântica)、卡廷加(Caatinga)和亚马逊(Amazon)生物群落的碳平衡为负,其中富含氧化铁的土壤其碳储存能力低于理论值。我们的方法可以作为农业生态系统中碳固定潜力的近似值,并有助于减缓气候变化。

**通俗语言总结**

土壤是地球上最大的碳储存库,其作为大气中CO2的来源或汇取决于管理方式。土壤表面释放的CO2(FCO2)与土壤的生物和物理化学性质密切相关。我们的目标是利用可见光、近红外-短波红外和中红外反射光谱技术、数字土壤制图以及机器学习方法,估算并绘制巴西境内的净生态系统生产力(NEP)空间分布图。我们建立了以土壤物理、化学和微生物性质为协变量的FCO2和碳封存潜力预测模型。NEP地图显示,塞拉多、大西洋森林、卡廷加和亚马逊生物群落的碳平衡为负,其中富含氧化铁的土壤其碳储存能力低于理论值。我们的方法可以作为农业生态系统中碳固定潜力的近似值,并有助于减缓气候变化。

**缩写说明**

- DSM:数字土壤制图(Digital Soil Mapping)
- EC:涡度协方差(Eddy Covariance)
- GHG:温室气体(Greenhouse Gases)
- LST:地表温度(Land Surface Temperature)
- MAOC:与矿物部分相关的有机碳(Carbon in Organic Matter Associated with the Mineral Fraction)
- MAOM:与矿物部分相关的有机质(Organic Matter Associated with the Mineral Fraction)
- MBC:微生物生物量碳(Microbial Biomass Carbon)
- MIT:Moran指数检验(Moran Index Test)
- MRA:多元线性回归(Multivariate Linear Regression)
- NEP:净生态系统生产力(Net Ecosystem Production)
- NPP:净初级生产力(Net Primary Productivity)
- PCA:主成分分析(Principal Component Analysis)
- POC:颗粒有机碳(Carbon in Particulate Organic Matter)
- POM:颗粒有机质(Particulate Organic Matter)
- RF:随机森林(Random Forest)
- RFE:递归特征消除(Recursive Feature Elimination)
- RMSE:均方根误差(Root Mean Square Error)
- RPIQ:四分位距(Interquartile Range)
- SEM:空间误差模型(Spatial Error Models)
- SOC:土壤有机碳(Soil Organic Carbon)
- SOM:土壤有机质(Soil Organic Matter)
- SRDB:全球土壤呼吸数据库(Global Soil Respiration Database)因此,本研究的目的是开发一种基于可见光近红外-短波红外(Vis NIR-SWIR)和中红外(Mid-IR)反射光谱产品、数字表面模型(DSM)以及机器学习技术的技术,通过建立和空间化土壤二氧化碳排放潜力和有机质(SOC)固存潜力模型,来估算和绘制巴西领土的净生态系统生产力分布图。

2 方法论

2.1 研究区域

研究区域覆盖了整个巴西领土(850万平方公里)。由于气候、母质和地形等形成因素的影响,该国的土壤类型多样。根据柯本气候分类(Köppen classification),巴西的气候类型包括热带气候(81.4%)、干旱气候(4.9%)和亚热带气候(13.7%)(Alvares等人,2013年)(图1)。土壤母质主要形成于新生代、中生代、新元古代和古元古代。这些土壤来源于各种类型的变质岩、火成岩和沉积岩(Gómez等人,2019年)。地形以缓坡为主,海拔范围在200至500米之间,地质活动及火山活动较少(Ross,2013年)。图1展示了研究区域和方法流程图。巴西土壤光谱库(BSSL)数据集中的土壤观测点深度为0-20厘米,数据包括可见光近红外-短波红外光谱(黄色)和中红外光谱(橙色)。含有微生物分析结果的样本(β-葡萄糖苷酶、尿素酶、酸性磷酸酶活性和微生物生物量碳)用红色标注。土壤二氧化碳排放的点测量值用紫色表示。柯本气候分类:A = 热带气候区,无干季(Af);B = 半干旱气候区,低纬度和低海拔(BSh);C = 亚热带气候区,无干季,夏季炎热(Cfa)或夏季温和(Cfb),或冬季干燥夏季炎热(Cwa)或夏季温和(Cwb),数据来源于Alvares等人(2013年)。全球土壤呼吸数据库“SRDB”由Bond-Lamberty和Thomson(2010年)以及Jian等人(2021年)开发。

2.2 土壤观测

在本研究中,我们使用了来自巴西土壤光谱库(BSSL)的包含可见光近红外-短波红外(350–2500纳米)光谱范围的遗留数据库,光谱分辨率为1纳米,其中包含69,740个观测数据,其中15,083个与0–20厘米深度相关(图1)。此外,还使用了2,866个中红外(4000–6000厘米^-1)范围内的观测数据,光谱分辨率为1厘米^-1。

2.3 酶活性、微生物生物量碳和土壤有机质(SOM)质量的光谱传递函数

2.3.1 实验室分析

为了构建β-葡萄糖苷酶、尿素酶、酸性磷酸酶和微生物生物量碳(MBC)的活性指数,在圣保罗州的皮拉西卡巴(Piracicaba)、帕拉州(Pará)的托梅-阿苏(Tomé-Açu)市、马托格罗索杜苏尔州(Mato Grosso do Sul)的班代兰特斯(Bandeirantes)市以及巴拉那州(Paraná)的贾丁奥林达(Jardim Olinda)市收集了35个样本(图1)。β-葡萄糖苷酶、酸性磷酸酶和尿素酶的酶活性测定遵循Tabatabai(1994年)和Dick等人(1996年)的方法。详细信息可见H. S. Rodríguez-Albarracín等人(2025年)的研究。MBC的量化按照Vance等人(1987年)的方法进行。此外,还对210个样本进行了物理分级,以量化颗粒有机质(POM)中的碳和与矿物组分相关的有机质(MAOM)中的碳(Cotrufo等人,2019年)。样本的选择基于MBC的预测值、β-葡萄糖苷酶、酸性磷酸酶和尿素酶的酶活性、POM中的碳(CPOM)和MAOM中的碳(CMAOM),以及真菌和放线菌的丰度。这些预测值是使用H. S. Rodríguez-Albarracín等人(2024年,2025年)提出的光谱传递方程为BSSL中的样本确定的。通过模糊聚类技术从该数据库中选择了210个样本,该技术考虑了空间位置,确保样本能够被分配到不重叠的组中,从而在不同特征空间和空间背景下进行采样。这种技术是在R语言的ClustGeo库中实现的,其中使用了自动优化的层次聚类参数(Chavent等人,2021年)。这种方法允许根据对该地区的先验知识或可用信息调整组的配置。分级方法是通过用5%的六偏磷酸钠溶液分散土壤来进行的;随后,通过53微米的筛子分离沙子和粘土/淤泥组分(Jindaluang等人,2013年)。这里,POC(颗粒有机质中的碳)指的是POM组分中的碳,而MAOC(与矿物组分相关的有机质中的碳)是指MAOM中的碳含量。

2.3.2 光谱数据采集

光谱数据使用两种传感器在可见光近红外-短波红外(Vis-NIR-SWIR)和中红外(Mid-IR)范围内采集:FieldSpec 3光谱辐射计(Analytical Spectral Devices),用于获取350至2500纳米范围内的反射数据(Vis-NIR-SWIR),其光谱分辨率在400至700纳米范围内为3纳米,在700至2500纳米范围内为10纳米。为了提高数据质量,进行了1纳米的重采样,获得了2151个光谱带(Greschuk等人,2022年)。使用该传感器进行读数时,需要将样本筛分至2毫米。我们还使用了Alpha Sample Compartment RT-DLaTGS ZnSe传感器(Bruker Optik GmbH),配备了漂移附件(Souza等人,2020年),用于获取中红外(4000–600厘米^-1)范围内的光谱信息,分辨率为1.2纳米,扫描速率为64次每秒(Terra等人,2015年)。使用该传感器进行读数时,需要将样本研磨并筛分至100目。

2.3.3 光谱传递函数

我们使用了Guyon等人(2002年)开发的递归特征消除(RFE)算法,并在R语言的“Caret”包中实现,内部模型采用随机森林(RF)机器学习方法(M. Kuhn,2021年),默认使用基尼指数(Gini index),来选择每个活性指数最重要的光谱带集合。RFE使我们能够获得每个变量最重要的Vis-NIR-SWIR光谱带,然后基于这些光谱信息构建回归模型进行量化。测试了以下模型:空间自相关模型、空间滞后模型、空间误差模型(SEM)和空间自回归-自回归模型(SARAR),后者包括响应的自回归成分和残差,从而解释残差的空间依赖性(Elhorst,2014年),详细信息见H. Rodríguez-Albarracín等人(2023年)。此外,还测试了广泛用于定义土壤传输函数的多变量线性回归(MRA)模型(Perreault等人,2022年),其中考虑了每个变量最重要的光谱带的反射值作为解释变量。通过逐步方法引入或消除预测变量来优化模型的拟合度(Dos Santos等人,2025年)。基于观测变量Y的观测值与模型估计值之间的相关系数(R)选择了最佳拟合模型(Hoge等人,2018年)。此外,还使用了Moran指数测试(MIT)来考虑残差的独立性,其中权重矩阵包括所有邻居(Liu & Chen,2021年),如果测试的p值大于0.05,则认为残差是独立的。基于这些模型,计算了遗留数据库中每个观测值的各个变量。

2.4 土壤有机质(SOM)的化学组成

与二氧化碳排放相关的SOM质量指数是根据Tavanti等人(2023年)的提议定义的,他们开发了考虑中红外光谱2853和2923厘米^-1波段的CH2指数。根据不同的中红外光谱范围和不同研究中收集的信息,定义了不同的SOM指数,详见表1。表1中显示的波段被选中,并为每个相关组绘制了地图。表1. 与土壤有机质(SOM)化学组成相关的中红外波段。

2.5 空间预测

空间预测使用了RF算法,该算法在预测和绘制各种土壤属性(如SOC)方面表现最佳(Gómez等人,2023年;Padarian等人,2020年;Zeraatpisheh等人,2020年)。由于国家级别映射涉及的数据量庞大,每个观测值计算出的结果通过Google Earth Engine Python API上传到了Google Cloud Storage(Gorelick等人,2017年)。地理空间映射遵循van den Hoogen等人(2021年)提出的流程。最终用于创建预测图的模型是基于R2(决定系数)、均方根误差(RMSE)以及校准和验证数据集之间的性能和四分位数范围(RPIQ)关系的网格搜索程序得出的最佳模型。R2系数表示解释的方差,RMSE表示准确性,RPIQ评估数据四分位数范围内的误差(Rosin等人,2023年)。作为映射的协变量定义了以下因素:使用Demattê等人(2018年)和Demattê等人(2020年)提出的地理空间土壤感测系统方法获得的合成土壤图像;地形属性(海拔、坡度、朝向、曲率、山谷、山丘和地形湿度指数)使用GEE包中的Terrain Analysis功能(Safanelli等人,2020年);以及气候(温度和年平均降水量)来自WorldClim2(Fick & Hijmans,2017年)。赤铁矿、针铁矿、高岭石和 gibbsite的相对丰度地图来自Rosin等人(2023年),并根据Dinerstein等人(2017年)报告的全球生态区地图进行了生物群落分层。对于风化指数(Ki)地图,使用了Rosas等人(2025年)提供的氧化铝(Al2O3)和二氧化硅(SiO2)地图,并应用了Embrapa Solos(1997年)提出的公式。

2.6 土壤二氧化碳通量的量化和回归分析

整合了一个全国性的土壤二氧化碳通量(FCO2)数据库(图1)。在四所大学的支持下,我们获得了马托格罗索州(Cuiabá、Diamantino、Juscimeira、Lucas do Rio Verde、Nobres、Nova Mutum、Rondonópolis、Rosário Oeste、Sinop、Sorriso、Vera、Costa Rica和Itiquira市);马托格罗索杜苏尔州(Chapadão do Sul、Aquadauana和Deodápolis市);以及伯南布哥州(Serra Talhada、Triunfo、São João和Custódia市)的特定测量数据。这些测量数据使用便携式自动化LICOR 8100和EGM-5(PP-Systems)AGA560红外二氧化碳土壤气体分析仪进行。详细的方法论描述在Della Silva等人(2025年)、Teodoro等人(2024年)、Moitinho等人(2021年)和Cicuendez等人(2020年)的研究中。此外,还使用了Bond-Lamberty和Thomson(2010年)开发的全球土壤呼吸数据库“SRDB”中的80%数据(70个观测值),该数据库已更新至Jian等人(2021年)的第五版(图1)。该数据库包含分布在巴西各地的87个观测数据,每年测量一次FCO2。SRDB数据库中的数据使用R语言的bigleaf包(Knauer,2022年)调整为µmol C m^-2 s^-1单位。鉴于FCO2的空间变化可以通过土壤的物理化学和生物特性来解释,而不是气候因素(Moitinho等人,2021年;Yang等人,2023年),我们考虑了以下解释变量来确定FCO2预测函数:MBC(微生物生物量碳)、β-葡萄糖苷酶的酶活性、酸性磷酸酶、脲酶、POC(颗粒有机碳)、MAOC(微生物可利用有机碳)、质地(沙子、粉砂和粘土)、表1中描述的与SOM(土壤有机质)化学组成相关的组别,以及Ki(土壤容重)。进行了逐步回归分析,以根据它们的相对贡献确定最重要的变量(Xiao等人,2019年)。还测试了MRA模型,同时考虑了X和Y坐标作为解释变量,如Böhner和Bechtel(2018年)所描述的。这种分析被认为是线性空间趋势函数的扩展,如果残差的独立性假设得到满足,它可以被用作预测和空间化未知值的函数。利用得到的方程,通过地图代数生成了整个巴西的FCO2地图。地图代数有助于以易于理解的方式表示和处理数据(Pullar,2003年),它基于影响单个单元属性的点特定或局部操作(Schmitz等人,2013年)。

2.7 碳封存潜力

根据Georgiou等人(2022年)的研究,最大SOC(土壤有机碳)含量取决于矿物质的量和类型,因此开发数学模型来估计与矿物质部分相关的碳是很重要的。在这方面,对于本研究,确定了潜在的C饱和度,即细小土壤颗粒(粉砂和粘土部分<20 µm)可以稳定的理论最大SOC,使用H. Rodríguez-Albarracín等人(2023年)的方法来调整这一潜力(最初是基于粘土含量确定的),同时考虑Rosin等人(2023年)报告的赤铁矿、针铁矿、高岭石和三水铝矿的相对丰度。因此,进行了空间回归分析,以潜在的C饱和度(Sat-pot)作为因变量,使用Feller和Beare(1997年)的方程(适用于热带土壤)作为解释变量。如“4.1.1”节所述,测试了空间回归模型,并且与该节中描述的程序类似,选择了满足残差独立性假设的模型(基于MIT,X. Liu & Chen,2021年),在多个模型都满足这一假设的情况下,使用的标准是观测值和估计值之间的最高相关性(R)(Hoge等人,2018年)。随后,使用地图代数计算了潜在SOC饱和度与MAOC之间的差异,以获得碳封存潜力的空间分布。

2.8 土壤二氧化碳排放(FCO2)的空间验证

Jian等人(2021年)对FCO2地图和SRDB数据库中20%的观测数据进行了相关性分析,覆盖了巴西的领土。由于Spearman相关系数在建模线性和非线性单调关系方面的优势,因此选择了它(L. Zhang & Wang,2023年)。此外,通过选择5762个点,对FCO2地图与年降水量地图、地表温度(LST)、海拔、SOC饱和潜力(C Sat-Pot)、碳封存潜力(C Satdef)、平均总生物量、三水铝矿、赤铁矿、针铁矿、高岭石、氧化铁(Fe2O3)、氧化铝(Al2O3)和二氧化硅(SiO2)进行了相关性分析。还进行了主成分分析(PCA),以分析变量之间的隐藏结构(Moitinho等人,2021年)。不同的作者指出,FCO2的时空变化受到土壤温度(J. R. D. S. Lima等人,2023年;Moitinho等人,2021年;Teodoro等人,2024年)和降水量(Zhao等人,2016年)的驱动,因为它们调节了微生物过程(W. Liu等人,2009年)。因此,我们包括了这些变量来验证FCO2地图,以解释获得的空间分布。

2.9 土壤二氧化碳排放(FCO2)地图的应用:净生态系统生产力(NEP)

NEP的确定基于Rs(根系呼吸)与植物生产之间的相关性,因为光合产物的供应和根系呼吸之间存在联系(Bond-Lamberty等人,2004年)。根据Phillips等人(2017年)和Lv等人(2022年)的研究,生态系统的碳平衡是基于NEP来评估的,NEP被定义为生态系统的光合生产力。我们使用作者提出的以下方程:
(2)
其中NPP对应于净初级生产力,单位为t C hm−2,是根据总生物量计算的,总生物量是地上和地下生物量的总和(Braswell等人,2005年;Dubey & Lal,2009年)。Rm指的是土壤中的微生物异养呼吸,由Rs确定,相当于t C hm−2的累积土壤呼吸(Bond-Lamberty等人,2004年;Phillips等人,2017年)。这种累积FCO2的计算使用了Al-Kaisi和Grote(2007年)提出的方程。常数0.865是基于Kuzyakov(2002年)的研究定义的,该研究分析了分离实际根系呼吸和Rh的四种方法。Raich和Schlesinger(1992年)指出,根系呼吸对Rs的贡献为30%。Hanson等人(2000年)报告称,在森林植被和非森林植被下,根系呼吸的贡献分别为45.8%和60.4%。Xu和Shang(2016年)指出,根系对Rs的贡献在不同生态系统之间差异很大,报告的范围从2%到90%不等。我们考虑使用Lv等人(2022年)提出的方程作为分析巴西领土上NEP分布的方法论方法;然而,需要注意的是,对于详细分析,必须根据植被情况调整方程。另一方面,NPP是从中等分辨率成像光谱辐射计(MODIS)MOD17A3HGF产品中得出的,它对应于年净光合作用的总和减去永久性木质组织中活细胞的生长和维护成本,单位为kg C m−2(Running & Zhao,2019年)。需要注意的是,MOD17算法使用的方法区分了总生产和自养呼吸的过程。后者包括活木质组织的维护呼吸以及叶子和细根的维护呼吸。从2003年1月1日到2023年12月31日的图像序列中计算了年平均值,考虑了SRDB数据库观测的时间窗口。由于这些时期的数据布局错误,没有考虑MODIS系列的2000-2002年(Running & Zhao,2019年)。这个年平均值被转换为t C hm−2。在大规模上,Rs和植物生产是相关的过程(Bond-Lamberty等人,2004年),生物和非生物因素干扰FCO2(Cicuendez等人,2020年;Yang等人,2023年)和生物量生产,这也影响了土壤微生物群落和FCO2(Giles等人,2023年)。总体而言,NEP整合了与呼吸和光合作用相关的多个过程,涉及土壤和植物(Phillips等人,2017年)。因此,我们使用RF机器学习方法来识别气候和土壤因素之间的不同关系及其对生物群落水平NEP的影响,因为这些因素的影响在生态系统和生物群落之间有很大的差异(Bond-Lamberty & Thomson,2010年;Teodoro等人,2024年)。我们决定使用这种RF方法,因为它不受预测变量之间线性依赖性的影响(Breiman,2001年)。我们使用了R包“Caret”和“train”函数,该函数允许通过重采样来评估模型调整参数对性能的影响(D. Kuhn,2019年)。RF模型允许我们对每个变量的重要性进行排名,并探索其对分析变量(NEP)的个别影响(Georgiou等人,2022年),使用部分依赖图(Friedman,2001年)。

3 结果

3.1 酶活性、MBC(微生物生物量碳)和SOM(土壤有机质)化学组成的空间预测

3.1.1 光谱传递函数

使用Vis-NIR-SWIR光谱范围对MAOC的预测进行了调整,以适应空间SEM误差模型(R = 0.87,MIT p值 = 0.24):

MAOC(mg g−1)= 0.89 + 26.95 (355) − 75.30 (356) + 61.64 (357) − 70.57 (380) + 55.23 (386) − 4.05 (566) − 483.29 (1411) + 1480.40 (1413) − 990.56 (1414) + 177.48 (1881) − 528.13 (1884) + 359.17 (1886) + 172.14 (1892) − 351.91 (1896) + 747.36 (1902) − 570.16 (1903) − 757.58 (2135) + 1036.33 (2136) − 664.13 (2138) + 396.56 (2139) − 25.33 (2159) + 393.09 (2264) − 388.24 (2265) + 12.17 (2497) − 22.84 (2498) + 9.14 (2500) + u:

W对应于权重矩阵,u与具有空间依赖性的残差向量相关,其中I是单位矩阵。POC的预测调整到了MRA模型(R = 0.64,MIT p值 = 0.052):

2. POC(mg g−1)= 0.05 + 12.99 (358) − 13.78 (359) − 16.00 (366) + 14.90 (373) + 30.42 (376) − 21.66 (377) − 56.34 (464) + 59.98 (502) − 17.19 (567) + 24.72 (616) − 25.59 (644) − 39.87 (762) + 663.35 (800) − 857.35 (806) + 304.51 (820) − 121.15 (861) + 58.00 (897) − 6.09 (1280) + 6.63 (1426) − 5.26 (1889) + 3.22 (1946) + 42.64 (2207) − 43.55 (2210) + 10.49 (2340) − 7.71 (2389) − 2.45 (2495)。对于β-葡萄糖苷酶、酸性磷酸酶和脲酶酶活性的预测,使用了H. S. Rodríguez-Albarracín等人(2025年)报告的方程,采用Vis-NIR-SWIR范围:

3. β-葡萄糖苷酶(mg PNB kg−1 soil)= 57.17 + 3140.16 (360) − 5047.09 (361) + 4507.84 (376) − 2815.09 (423) + 412.16 (535) − 129.69 (2500);R = 0.76。

4. 尿酶(µg NH4 kg−1 soil)= 12.30 + 1411.36 (410) − 1500.72 (463) + 949.22 (663) − 4750.01 (798) + 4277.31 (821) − 143.03 (1106) − 50.16 (2470);R = 0.84。

5. 酸性磷酸酶(mg PNP kg−1 soil)= 539.16 + 3340.61 (350) + 10,585.21 (356) − 13,872.10 (360) − 9709.36 (576) + 43,224.81 (667) − 536,484.72 (736) + 502,451.57 (740) − 1037.61 (2441);R = 0.68。对于MBC的预测,使用了H. S. Rodríguez-Albarracín等人(2024年)报告的方程,采用Vis-NIR-SWIR范围:

6. MBC(mg C g−1 soil)= 0.57 − 21.56 (356) + 8.74 (360) + 8.90 (371) − 2.66 (2300) + 2.73 (2467);R = 0.69。具体的Vis-NIR-SWIR(nm)波段以斜体显示,它们的反射率值被代入方程中。

3.1.2 空间预测

获得了解释变量的地图,覆盖了整个巴西领土,R2值范围从0.41到0.76,其中CH2指数的R2最低(表2)。β-葡萄糖苷酶酶活性获得了最高的R2,与与SOM化学组成相关的解释变量类似。在农场级别,CH2指数的R2值为0.62,在Tavanti等人(2023年)的研究中,SOM化学组的R2值为0.75–0.86。这些值与我们在巴西级别获得的值一致,范围从0.71到0.74,除了CH2指数,其R2为0.41。表2. FCO2解释变量映射的验证指标。解释变量
R2
RMSE
RPIQ

CH2 对称
0.719
0.007
2.668

CH2 不对称
0.713
0.007
1.833

CH3 对称 I
0.719
0.007
1.798

CH3 对称 II
0.737
0.007
2.774

CH3 不对称 I
0.741
0.007
2.703

CH3 不对称 II
0.723
0.007
2.019

CH
0.707
0.008
2.631

CH 醛
0.729
0.007
2.76

CH 脂肪族
0.708
0.007
2.456

CH 脂肪族不对称
0.708
0.007
2.547

CH2 指数
0.408
0.065
1.291

MAOC
0.490
0.404
0.857

POC
0.674
0.148
1.037

MBC
0.736
0.119
1.606

沙子
0.713
150.020
3.444

粉砂
0.504
69.284
1.338

Ki
0.460
58.860
2.130*

酸性磷酸酶
0.559
90.726
1.573

脲酶
0.692
9.268
0.986

β-葡萄糖苷酶
0.764
28.197
0.846

注:验证指标对应于Rosas等人(2025年)制作的SiO2和Al2O3地图。缩写:MAOC,与矿物质相关的有机碳;MBC,微生物生物量碳;POC,颗粒有机碳;RMSE,均方根误差。

3.2 土壤CO2通量的量化和回归分析

FCO2预测模型适配到了一个MRA模型,具有稳健的相关性R为0.68和R2为0.35(图2),其中包括β-葡萄糖苷酶、酸性磷酸酶和脲酶的酶活性;沙子含量;Ki;MAOC;对称和不对称CH2;醛类CH;以及CH2指数。最重要的变量包括MAOC(与矿物质结合的有机碳)、Ki指数、CH2指数、沙含量、对称CH2以及β-葡萄糖苷酶活性。图2(在图查看器中打开)

用于预测FCO2的多变量回归模型(MRA)中,粗体数字代表每个变量的参数,*表示p值。MAOC是指与矿物质结合的有机碳;SOM表示土壤有机质。利用得到的方程式(图2和图3),通过地图代数方法绘制了整个巴西地区的FCO2分布图,结果显示阿克里州、亚马逊西部、帕拉州北部和南部、阿马帕州东部、托坎廷斯州、皮奥伊州、马托格罗索州北部、巴拉那州南部以及南里奥格兰德州西部和北部的排放量较高。排放量最低的州位于该国东部,包括米纳斯吉拉斯州、里约热内卢州、埃斯皮里图桑托州、巴伊亚州、戈亚斯州、伯南布哥州、帕拉伊巴州以及亚马逊州中南部。图3(在图查看器中打开)

土壤CO2通量的预测模型和地图。解释变量的地图包括:沙含量(g kg−1)、β-葡萄糖苷酶(mg PNB kg−1 soil)、酸性磷酸酶(mg PNP kg−1 soil)、脲酶(µg NH4 kg−1 soil)和MAOC(mg g−1,即与矿物质结合的有机碳)。地图代数的原理是通过组合栅格地图来生成新的栅格地图(Yamamoto等人,2018年),在此过程中,涉及变量的算术运算所产生的不确定性也会被传递(Mood等人,1974年)。在本研究中,结合了图2中方程式所指定的所有地图,因此最终的FCO2地图包含了图1和图2中显示的不确定性。

3.3 碳封存潜力

C Sat-Pot数据被拟合到SEM模型中(R = 0.90,MIT p值 = 1):

对于这种类型的模型,使用R语言中的Nagelkerke函数计算伪R2值,该函数基于似然对数来指示模型的准确性(Freese & Scott Long,2006年)。最重要的变量是吉布石和高岭石的振幅,其次是赤铁矿(图4)。图4(在图查看器中打开)

碳封存潜力的建模。C Sat-pot = 潜在碳饱和度;A = 不同矿物的振幅;AKln(高岭石)、AGt(针铁矿)、AHem(赤铁矿)、AGbs(吉布石);MAOC(mg g−1)= 与矿物质结合的有机碳。粗体数字代表每个变量的参数,*表示p值。MIT表示Moran指数检验;SOC表示土壤有机质。Georgiou等人(2022年)指出,确定土壤的最大SOC封存潜力取决于矿物的数量和类型。我们的结果与作者的发现一致,如图4所示,它说明了不同矿物在预测潜在SOC饱和度方面的重要性;此外,这种矿物学差异也体现在MAOC和碳封存潜力的空间分布上(图3和图4)。

3.4 土壤二氧化碳排放量(FCO2)的空间验证

Jian等人(2021年)的SRDB数据库中20%的观测数据与FCO2地图之间的相关系数R为0.56(图5)。此外,与Moitinho等人(2021年)报告的圣保罗州巴里尼亚地区的地图也存在正相关且相关性较高,尽管他们使用了基于地质统计学的不同空间化方法。碳封存潜力(R = 0.67)、赤铁矿振幅(R = 0.22)和年降水量(R = 0.20)之间也观察到正相关;而LST(R = −0.15)和Al2O3(R = −0.29)之间则呈现负相关。PCA中的变量交互图显示了三个组:第一组将FCO2与表现出正相关的变量归为一类;第二组与土壤矿物学和潜在碳饱和度相关,这些变量与FCO2呈负相关;第三组仅包含LST这一个变量,其与其他变量的关系在不同地区有所不同,在该国西北部、东北部和南部观察到正相关(图5)。Teodoro等人(2024年)报告称,在南马托格罗索州的塞拉多生物群落中,FCO2与LST之间存在正相关,这与我们的结果一致。图5(在图查看器中打开)

空间验证:分析了Spearman相关性,包括:(A) CO2通量地图与FCO2的具体观测数据(20% Bond-Lamberty & Thomson(2010)数据库);(B) FCO2地图与年降水量、地表温度(LST)、海拔高度、潜在碳饱和度(C Sat-Pot)、碳封存潜力(C Satdef)、净初级生产力(NPP)、吉布石、赤铁矿、针铁矿、高岭石、氧化铁(Fe2O3)、氧化铝(Al2O3)和二氧化硅(SiO2)的地图;(C) 主成分分析(PCA)图,显示了变量之间的相互作用。观察到FCO2、碳封存潜力和赤铁矿含量之间存在正相关(图5),因为赤铁矿含量较高的地区具有更高的碳封存潜力,并且这些地区的酶活性也较高(图3)。在国家北部,降水量与FCO2之间也观察到正相关(R = 0.2),这些地区的生物量生产也较高。此外,Al2O3与FCO2之间呈现负相关(R = −0.29),其中氧化物含量较高的地区与Ki值较低以及高岭石和吉布石含量较高相吻合(Rosas等人,2025年),同时MAOC含量也较高(图3和图5)。

3.5 土壤二氧化碳排放量(FCO2)地图的应用:净生态系统生产力(NEP)

根据Lv等人(2022年)的研究,正的NEP值表示生态系统实现了碳封存,而负值表示二氧化碳排放。图6中,绿色表示负值,接近零的值用橙色和黄色表示,正值用红色表示。最高的排放量出现在塞拉多(北部和南部)、大西洋森林(南部)、卡廷加(北部)和亚马逊(西北部)生物群落。这些亚马逊地区的碳封存潜力较高,一些地区也呈现正平衡。然而,总体而言,正平衡地区的碳封存潜力较低,这是大多数生物群落中最重要的变量,大西洋森林除外,那里最重要的变量是Al2O3、年降水量、海拔高度和LST。Al2O3也是所有生物群落NEP的重要变量。针铁矿在潘帕地区也很重要,在这些地区,负平衡和较高碳封存潜力与这种矿物的较高含量相关,因为其具有较高的吸附能力(Pereira等人,2020年)。这些地区非常重要,因为富含氧化物的土壤如果低于其碳储存能力,则具有短期内提高土壤碳封存潜力的巨大潜力(Bramble等人,2023年)。换句话说,必须考虑土地利用方式以促进这些地区的碳封存并实现正平衡。图6(在图查看器中打开)

净生态系统生产力(NEP)。地表温度(LST)、净初级生产力(NPP)、潜在碳饱和度(C SatPot)、碳封存潜力(C Satdef)、吉布石(GB)、赤铁矿(HM)、针铁矿(GT)、高岭石(KT)、氧化铁(Fe2O3)、氧化铝(Al2O3)和二氧化硅(SiO2)。

4 讨论

4.1 酶活性、MBC以及SOM的质量和化学组成的空间预测

4.1.1 光谱传输函数

值得注意的是,方程式(1)和(2)分别对应于MAOC和POC预测模型,包含了26个光谱带。其中一些光谱带与最稳定的碳形式相关,这些碳形式会与矿物表面相互作用,例如氧化铁(Stenberg等人,2010年;Zornoza等人,2008年)。POC光谱带与微生物优先利用的易分解功能团相关,如多糖、甲基、脂肪族CH和酚类OH(Cozzolino & Morón,2003年;Fidêncio等人,2002年;Stenberg等人,2010年;Viscarra-Rossel & Behrens,2010年)。这些光谱带与H. S. Rodríguez-Albarracín等人(2024年,2025年)报告的总碳、磷酸酶和脲酶活性、POC以及MAOC的光谱带一致。SOM通常由多种化学化合物组成,包括碳水化合物、木质素和纤维素(Obeng等人,2023年),每种化合物都具有与其组成功能团相关的独特吸收特性;这解释了为什么选择这些光谱带(Stenberg等人,2010年)。然而,由于Vis-NIR-SWIR光谱包含2151个光谱带,我们的选择方法仅保留了原始特征的1.2%,这与Carnero等人(2024年)的报告相似。

4.2 土壤CO2通量的定量和回归分析

Della Silva等人(2025年)基于土壤湿度以及PlanetScope和MODIS的蓝光、绿光、红光和近红外光谱带以及LST数据,构建了预测巴西南部亚马逊地区和马托格罗索州北部FCO2的方程式,并使用偏最小二乘回归方法,其中三个模型的R2值为0.58。我们的FCO2预测模型在整个巴西地区的R2值为0.35(图2),该模型考虑了微生物酶活性、沙含量、Ki值、MAOC以及与对称和不对称CH2、CH醛和CH2指数相关的SOM化学性质。Tavanti等人(2023年)发现CH和对称CH3之间存在高但负的相关性,而与POM、对称CH2、不对称CH2和CH2指数之间存在正相关性。我们的模型证实了CH2指数以及对称和不对称CH2在预测FCO2中的重要性(图2)。此外,该模型还考虑了β-葡萄糖苷酶、酸性磷酸酶和脲酶的酶活性以及沙含量的影响。Yang等人(2023年)也强调了SOM化学性质和微生物组成在预测土壤呼吸作用中的重要性,但与我们的研究(图2)不同,他们认为纹理的贡献较低。这些与SOM化学性质之间的强相关性是因为这些化合物的分子量较低,更容易被微生物分解(Tavanti等人,2023年),因为它们是微生物的主要能量来源(Cotrufo等人,2015年)。Yang等人(2023年)还发现沙含量与FCO2之间存在强正相关,而在国家北部地区,降水量也与FCO2呈正相关,这些地区的生物量生产也较高。此外,Al2O3与FCO2之间存在负相关(R = −0.29),其中氧化物含量较高的地区与Ki值较低以及高岭石和吉布石含量较高相吻合(Rosas等人,2025年),同时MAOC含量也较高(图3和图5)。

4.2 土壤CO2通量的定量和回归分析

Della Silva等人(2025年)基于土壤湿度以及PlanetScope和MODIS的蓝光、绿光、红光和近红外光谱带以及LST数据,构建了预测巴西南部亚马逊地区和马托格罗索州北部FCO2的方程式,并使用偏最小二乘回归方法,其中三个模型的R2值为0.58。我们的FCO2预测模型在整个巴西地区的R2值为0.35(图2),该模型考虑了微生物酶活性、沙含量、Ki值、MAOC以及与对称和不对称CH2、CH醛和CH2指数相关的SOM化学性质。Tavanti等人(2023年)发现CH和对称CH3之间存在高但负的相关性,而与POM、对称CH2、不对称CH2和CH2指数之间存在正相关性。我们的模型证实了CH2指数以及对称和不对称CH2在预测FCO2中的重要性(图2)。然而,该模型还考虑了醛类CH的作用,并强调了β-葡萄糖苷酶、酸性磷酸酶和脲酶的酶活性以及沙含量的影响。Yang等人(2023年)也强调了SOM化学性质和微生物组成在预测土壤呼吸作用中的重要性,但与我们的研究不同,他们认为纹理的贡献较低。这些与SOM化学性质之间的强相关性是因为这些化合物的分子量较低,更容易被微生物分解(Tavanti等人,2023年),因为它们是微生物的主要能量来源(Cotrufo等人,2015年)。Yang等人(2023年)还发现沙含量与FCO2之间存在强正相关,而在国家北部地区,降水量也与FCO2呈正相关,这些地区的生物量生产也较高。此外,Al2O3与FCO2之间存在负相关(R = −0.29),其中氧化物含量较高的地区与Ki值较低以及高岭石和吉布石含量较高相吻合(Rosas等人,2025年),同时MAOC含量也较高(图3和图5)。

4.1 酶活性、MBC以及SOM的质量和化学组成的空间预测

4.1.1 光谱传输函数

需要注意的是,方程式(1)和(2)分别对应于MAOC和POC预测模型,包含了26个光谱带。其中一些光谱带与最稳定的碳形式相关,这些碳形式会与矿物表面相互作用,例如氧化铁(Stenberg等人,2010年;Zornoza等人,2008年)。POC光谱带与微生物优先利用的易分解功能团相关,如多糖、甲基、脂肪族CH和酚类OH(Cozzolino & Morón,2003年;Fidêncio等人,2002年;Stenberg等人,2010年;Viscarra-Rossel & Behrens,2010年)。这些光谱带与H. S. Rodríguez-Albarracín等人(2024年,2025年)报告的总碳、磷酸酶和脲酶活性、POC以及MAOC的光谱带一致。SOM通常由多种化学化合物组成,包括碳水化合物、木质素和纤维素(Obeng等人,2023年),每种化合物都具有与其组成功能团相关的独特吸收特性;这解释了为什么选择这些光谱带(Stenberg等人,2010年)。然而,由于Vis-NIR-SWIR光谱包含2151个光谱带,我们的选择方法仅保留了原始特征的1.2%,这与Carnero等人(2024年)的报告相似。

4.2 土壤CO2通量的定量和回归分析

Della Silva等人(2025年)基于土壤湿度以及PlanetScope和MODIS的蓝光、绿光、红光和近红外光谱带以及LST数据,构建了预测巴西南部亚马逊地区和马托格罗索州北部FCO2的方程式,并使用偏最小二乘回归方法,其中三个模型的R2值为0.58。我们的FCO2预测模型的R2值为0.35(图2),该模型考虑了微生物酶活性、沙含量、Ki值、MAOC以及与对称和不对称CH2、CH醛和CH2指数相关的SOM化学性质。Tavanti等人(2023年)发现CH和对称CH3之间存在高但负的相关性,而与POM、对称CH2、不对称CH2和CH2指数之间存在正相关性。我们的模型证实了CH2指数以及对称和不对称CH2在预测FCO2中的重要性(图2)。然而,该模型还考虑了醛类CH的作用,并强调了β-葡萄糖苷酶、酸性磷酸酶和脲酶的酶活性以及沙含量的影响。Yang等人(2023年)也强调了SOM化学性质和微生物组成在预测土壤呼吸作用中的重要性,但与我们的研究不同,他们认为纹理的贡献较低。这些与SOM化学性质之间的强相关性是因为这些化合物的分子量较低,更容易被微生物分解(Tavanti等人,2023年),因为它们是微生物的主要能量来源(Cotrufo等人,2015年)。Yang等人(2023年)还发现沙含量与FCO2之间存在强正相关,而在国家北部地区,降水量也与FCO2呈正相关,这些地区的生物量生产也较高。此外,Al2O3与FCO2之间存在负相关(R = −0.29),其中氧化物含量较高的地区与Ki值较低以及高岭石和吉布石含量较高相吻合(Rosas等人,2025年),同时MAOC含量也较高(图3和图5)。

4.2 土壤CO2通量的定量和回归分析

Della Silva等人(2025年)基于土壤湿度以及PlanetScope和MODIS的蓝光、绿光、红光和近红外光谱带以及LST数据,构建了预测巴西南部亚马逊地区和马托格罗索州北部FCO2的方程式,并使用偏最小二乘回归方法,其中三个模型的R2值为0.58。我们的FCO2预测模型的R2值为0.35(图2),该模型考虑了微生物酶活性、沙含量、Ki值、MAOC以及与对称和不对称CH2、CH醛和CH2指数相关的SOM化学性质。Tavanti等人(2023)发现CH和对称CH3之间存在高但负的相关性,而与POM、对称CH2、不对称CH2和CH2指数之间存在正相关性。我们的模型证实了CH2指数以及对称和不对称CH2在预测FCO2中的重要性(图2)。然而,该模型还考虑了醛类CH的作用,并强调了β-葡萄糖苷酶、酸性磷酸酶和脲酶的酶活性以及沙含量的影响。Yang等人(2023年)也强调了SOM化学性质和微生物组成在预测土壤呼吸作用中的重要性,但与我们的研究不同,他们认为纹理的贡献较低。这些与SOM化学性质之间的强相关性是因为这些化合物的分子量较低,更容易被微生物分解(Tavanti等人,2023年),因为它们是微生物的主要能量来源(Cotrufo等人,2015年)。Yang等人(2023年)还发现沙含量与FCO2之间存在强正相关,而在国家北部地区,降水量也与FCO2呈正相关,这些地区的生物量生产也较高。此外,Al2O3与FCO2之间存在负相关(R = −0.29),其中氧化物含量较高的地区与Ki值较低以及高岭石和吉布石含量较高相吻合(Rosas等人,2025年),同时MAOC含量也较高(图3和图5)。

4.3 酶活性、MBC以及SOM的质量和化学组成的空间预测

4.1.1 光谱传输函数

需要注意的是,方程式(1)和(2)分别对应于MAOC和POC预测模型,包含了26个光谱带。其中一些光谱带与最稳定的碳形式相关,这些碳形式会与矿物表面相互作用,例如氧化铁(Stenberg等人,2010年;Zornoza等人,2008年)。POC光谱带与微生物优先利用的易分解功能团相关,如多糖、甲基、脂肪族CH和酚类OH(Cozzolino & Morón,2003年;Fidêncio等人,2002年;Stenberg等人,2010年;Viscarra-Rossel & Behrens,2010年)。这些光谱带与H. S. Rodríguez-Albarracín等人(2024年,2025年)报告的总碳、磷酸酶和脲酶活性、POC以及MAOC的光谱带一致。SOM通常由多种化学化合物组成,包括碳水化合物、木质素和纤维素(Obeng等人,2023年),每种化合物都具有与其组成功能团相关的独特吸收特性;这解释了为什么选择这些光谱带(Stenberg等人,2010年)。然而,由于Vis-NIR-SWIR光谱包含2151个光谱带,我们的选择方法仅保留了原始特征的1.2%,这与Carnero等人(2024年)的报告相似。

4.2 土壤CO2通量的定量和回归分析

Della Silva等人(2025年)基于土壤湿度以及PlanetScope和MODIS的蓝光、绿光、红光和近红外光谱带以及LST数据,构建了预测巴西南部亚马逊地区和马托格罗索州北部FCO2的方程式,并使用偏最小二乘回归方法,其中三个模型的R2值为0.58。我们的FCO2预测模型的R2值为0.35(图2),该模型考虑了微生物酶活性、沙含量、Ki值、MAOC以及与对称和不对称CH2、CH醛和CH2指数相关的SOM化学性质。Tavanti等人(2023)发现CH和对称CH3之间存在高但负的相关性,而与POM、对称CH2、不对称CH2和CH2指数之间存在正相关性。我们的模型证实了CH2指数以及对称和不对称CH2在预测FCO2中的重要性(图2)。然而,该模型还考虑了醛类CH的作用,并强调了β-葡萄糖苷酶、酸性磷酸酶和脲酶的酶活性以及沙含量的影响。Yang等人(2023)也强调了SOM化学性质和微生物组成在预测土壤呼吸作用中的重要性,但与我们的研究不同,他们认为纹理的贡献较低。这些与SOM化学性质之间的强相关性是因为这些化合物的分子量较低,更容易被微生物分解(Tavanti等人,2023年),因为它们是微生物的主要能量来源(Cotrufo等人,2015年)。Yang等人(2023)还发现沙含量与FCO2之间存在强正相关,而在国家北部地区,降水量也与FCO2呈正相关,这些地区的生物量生产也较高。此外,Al2O3与FCO2之间存在负相关(R = −0.29),其中氧化物含量较高的地区与Ki值较低以及高岭石和吉布石含量较高相吻合(Rosas等人,2025年),同时MAOC含量也较高(图3和图5)。图3显示,这些Ki值低于2的区域含有更多的高岭石和三水铝矿(Rosin等人,2023年),这些区域具有更高的β-葡萄糖苷酶和酸性磷酸酶活性,更高的MAOC含量,以及更低的FCO2值。该国北部和中部地区的较高排放量是由于土壤有机质的化学性质,特别是对称和非对称CH2、CH醛的含量较低,MAOC含量为中等到低水平,酸性磷酸酶活性为中等到高水平,β-葡萄糖苷酶活性也为中等水平(图3)。根据Marchiori和Melo(1999年)的研究,β-葡萄糖苷酶活性的增加可能表明易分解碳(labile C)的增加,从而导致分解过程的加剧,这会转化为更高的FCO2值(Lv等人,2022年)。易分解碳最初会刺激微生物活动,但这种活动会导致难分解碳(MAOC)浓度的增加(Simões等人,2023年),而难分解碳对微生物分解的抵抗力更强(Six等人,2002年)。对于Lv等人(2022年)的研究,可利用磷(available P)含量和微生物生物量碳(MBC)是FCO2的关键因素。在我们的模型中,酸性磷酸酶与FCO2呈正相关,这种酶与微生物丰度有关(Banerjee等人,2012年),并且通常与与金属氧化物表面相互作用的磷酸基团有关(Omoike & Chorover,2006年;Parikh等人,2014年),并且被认为是无机磷(P)可用性的指标(Piotrowska-Dlugosz & Charzynski,2015年)。

4.3 碳封存潜力

本研究开发的用于预测土壤有机碳(SOC)饱和潜力的模型(图4)与H. Rodríguez-Albarracín等人(2023年)的结果不同,他们报告称高岭石和赤铁矿在预测圣保罗州皮拉西卡巴地区的SOC饱和潜力方面更为重要。根据Schaefer等人(2008年)的研究,大多数巴西的砖红壤含有高岭石。对于Pereira等人(2020年)的研究,在大多数热带土壤中,更常见的成土铁氧化物是针铁矿,因为三水铝矿的形成过程可能较为罕见。然而,三水铝矿是一种常见的氧化物,在不同气候条件下可以在各种土壤中存在,当Ki值<2时,其存在量可能高于高岭石,这在巴西大部分地区都是如此(图3)。根据Resende(1976年)的研究,源自富铁贫硅岩石的巴西土壤更有利于三水铝矿的形成而非高岭石。碳封存潜力(图4)显示了具有新碳保留潜力的区域(黄色和红色区域),这些区域大多与最低的MAOC值相对应,而低MAOC值通常与较高的粘土含量区域相关。根据Sokol等人(2022年)的研究,MAOC约占SOC的65%。根据Rosin等人(2023年)的研究,这些区域还含有较高的针铁矿和赤铁矿含量。总体而言,针铁矿和三水铝矿被认为是SOC吸附的关键表面(Dos Reis等人,2014年;Kaiser & Zech,2000年)。Bramble等人(2023年)强调了针铁矿的吸附能力,它每平方米矿物可以吸附多达2.5毫克的CO。因此可以说,这些矿物的吸附能力已经接近其最大值。根据Georgiou等人(2022年)的研究,当矿物接近其最大碳吸附能力时,它们吸附碳的效率会降低,因为稳定的SOC含量越高,新碳的封存潜力就越低(H. Rodríguez-Albarracín等人,2023年)。这些主要矿物含量高的区域可能接近其最大吸附能力,这就是为什么高岭石和三水铝矿含量较高的区域(Rosin等人,2023年)具有较低的碳封存能力和较高的MAOC含量的原因。该国中部和南部的一些区域具有中等到高的碳封存能力(图4),这些区域也与最高的三水铝矿、高岭石和赤铁矿含量相对应(Rosin等人,2023年)。该国北部也具有较高的碳封存潜力,尽管主要矿物的含量较低,粘土含量平均,Ki值>1.85。这些区域具有较高的地上生物量生产,因此微生物活动较强,β-葡萄糖苷酶活性为中等到高水平,脲酶活性高,酸性磷酸酶活性低(图3和图4)。东北部地区的碳封存能力较低,因为该地区的沙子含量高,Ki值高于1.85,β-葡萄糖苷酶活性低,MAOC含量低。

4.4 土壤二氧化碳排放(FCO2)的空间验证

总体而言,发现FCO2与陆地表面温度(LST)之间存在负相关(图5)。不同作者指出,在温暖的气候条件下,温度升高会导致微生物种群减少和土壤有机质(SOM)分解强度降低(Correia等人,2015年;Ferreira等人,2018年;Holanda等人,2015年;Souto等人,2005年)。然而,在热带地区的沙质土壤中,在降水量增加和温度升高的时期,微生物活动会受到刺激,加速SOM的矿化过程,从而导致更高的CO2排放(Dos Santos等人,2025年)。一般来说,这些土壤由于温度升高而更容易失去碳(Hartley等人,2021年)。然而,由于热带地区的温度变化较小,这可能对FCO2的影响较小(Maia等人,2019年),相反,降水量成为关键因素(Zhang等人,2015年),特别是在干旱和半干旱生态系统中(Bond-Lamberty & Thomson,2010年;J. R. S. Lima等人,2020年;Zhao等人,2016年)。J. R. S. Lima等人(2022年)和J. R. D. S. Lima等人(2023年)报告了Caatinga生物群落中降水量与FCO2之间的正相关关系,以及与温度的负相关关系,这与我们的结果一致(图3和图5)。根据Ferreira等人(2018年)的研究,在温暖的气候条件下,较高的温度会减少微生物种群和SOM分解强度。相反,de Campos等人(2022年)报告了FCO2与降水量之间的负相关关系,以及与LST的正相关关系。根据作者的说法,温度升高会促进光合作用和CO2捕获能力;然而,这也会增加根系和微生物的活动,从而增加土壤中的FCO2。对于Bond-Lamberty等人(2004年)的研究,根系和微生物的呼吸作用对温度变化的反应可能不同。此外,还有研究表明,微生物代谢能量的增加会导致FCO2的增加(Moitinho等人,2021年)。这与该国北部观察到的降水量与FCO2之间的正相关关系一致(图5),该地区的生物量生产较高。此外,这些地区降水量增加,可以刺激微生物呼吸作用,通过增强细胞外酶活性(Zhao等人,2016年)。另一方面,总土壤FCO2对应于根系和土壤生物释放的CO2总和(Gao等人,2018年)。

4.5 土壤二氧化碳排放(FCO2)地图的应用:净生态系统生产力(NEP)

根据Lv等人(2022年)的研究,土壤中的CO2排放受到净初级生产力(NPP)的影响,NPP在计算净生态系统生产力(NEP)时是一个关键因素,在Caatinga、大西洋森林、亚马逊和Pampa生物群落中尤其如此(图6)。Giles等人(2023年)指出,高生物量与较高的FCO2相关。相反,da Costa等人(2022年)报告称,较高的FCO2与较低的光合作用速率相关。然而,特别是在亚马逊雨林中,根系呼吸作用占FCO2的24%到35%,并且与总生物量直接相关(X. Li等人,2018年)。Mendes等人(2025年)也报告了亚马逊森林中的高呼吸速率,根据位于该生物群落的涡度协方差(EC)塔的数据,NEP值接近中性,这与Aguilar(2005年)在亚马逊大规模生物圈-大气实验(LBA)参考区域内的Rio Jaru生物保护区报告的结果相似,其NEP值接近0。图6显示,亚马逊大部分地区的NEP为中性到正值。这与Bond-Lamberty等人(2004年)的研究结果一致,他们指出当FCO2值较高时,自养呼吸作用的贡献更大。因此,在这种类型的生态系统中,高排放值可能表明生产力水平较高(Simões等人,2023年)。Phillips等人(2017年)指出,FCO2主要由相对年轻的碳组成,这些碳在一年内转化,对稳定碳的贡献很小,这与亚马逊和Caatinga生物群落的特征一致,这些生物群落的MAOC含量较低(图3和图4)。Caatinga生物群落的生物量生产较低,净平衡为负值(图6),除了南部地区,由于CO2排放量低,净平衡为正值,这与较高的MAOC含量(图3)和更多的三水铝矿等矿物的存在相关(Rosin等人,2023年),这些矿物在NEP分析中是一个重要因素,因为它们具有较高的碳保留能力(Kaiser & Zech,2000年)。Simões等人(2023年)报告称,在巴西的半干旱地区,由于植物可用水分和大孔隙度的存在,系统的输入/输出较高。在这些地区,降水量增加会导致根系和微生物活动增加,从而产生更多的呼吸作用(Phillips等人,2017年;Tavares等人,2016年),这会导致易分解碳的快速氧化(Simões等人,2023年)。Mendes等人(2025年)和da Silva等人(2024年)报告称,Caatinga西北部的NEP值为中性到正值,该地区位于帕拉伊巴州和伯南布哥州之间,基于2011-2017年和2015-2019年期间国家水和碳动态观测网络中的EC通量塔的数据。我们的结果是基于2003年至2023年的多年平均值,与作者的研究结果一致,如图6所示,NEP为正值的区域表现出较高的三水铝矿含量,FCO2较低,MAOC含量为中等到高水平,碳封存潜力较低(图3和图4),表明这些矿物接近其饱和极限,这些区域的新碳封存增加与POC有关。因此,考虑到该地区受到降雨脉冲的影响,Caatinga植被对初雨的快速响应导致雨季开始时NPP增加,进一步促进了SOM分解和微生物活动(Alberton等人,2023年)。在旱季,根系呼吸作用减少,Rh受到干旱条件的限制,导致系统呈中性(Santos e Silva等人,2024年)。根据Coutinho等人(2019年)的研究,南马托格罗索州的温室气体排放量排名第11位,图6显示该州大部分地区的排放量较高,只有北部地区的排放量减少(Pantanal生物群落),那里也观察到了正的NEP平衡。与Teodoro等人(2024年)的研究结果相反,他们指出Cerrado生物群落的FCO2较低,而Pantanal和Atlantic Forest的排放量较高,我们的结果显示Atlantic Forest的排放量增加,而Pantanal的排放量减少。这些Pantanal生物群落的区域与较高的赤铁矿和针铁矿含量相关(Rosin等人,2023年),MAOC含量高,碳封存潜力较大(图3和图4)。换句话说,这些区域的矿物尚未达到其最大吸附能力。

5 限制

重要的是要认识到,本研究提出了一种基于土壤性质和光谱代理指标确定FCO2的方法,这种方法有助于减少传统的实验室分析并监测FCO2。因此,它既适用于详细研究也适用于半详细研究。所提出的国家级地图应被视为土壤FCO2状态的近似值,因为它是通过地图代数获得的,并包含了图3和图4中所示方程式中包含的属性的累积不确定性。此外,巴西直接观测土壤CO2的数量仍然有限且空间分布不均,应当注意的是,FCO2模型的校准使用了不能完全代表该国气候、土壤类型和土地利用多样性的数据集,这解释了R2值仅为0.35的原因。此外,需要明确的是,BSSL和SRBD数据库中的观测数据分别涵盖了40多年和20多年的信息,因此本研究中呈现的结果是基于多时间段的平均值。然而,通过对更多详细尺度研究的讨论可以认识到,所提出的方法能够生成可靠且可重复的结果,这些结果能够反映像巴西这样土壤变异性极大的国家的实际情况。鉴于FCO2模型的这些局限性,NEP地图也应被视为一种结构性估计;尽管每个生物群落的固有变异性已被识别,但本研究中用于确定NEP地图的Lv等人(2022年)的方程式使用了一个通用系数来从Rs值估算Rh值。不同的研究,如Kuzyakov(2002年)、Raich和Schlesinger(1992年)、Hanson等人(2000年)以及Xu和Shang(2016年)的研究表明,根系对Rs的贡献因生态系统和土地利用方式而异;因此,对于详细分析来说,考虑与研究区域条件相关的调整是非常重要的。

**6 结论**

我们开发了一种利用反射光谱数据预测和绘制FCO2及NEP地图的方法。我们建立了光谱传输函数,能够映射影响土壤FCO2的不同解释变量,其R2值范围从0.41到0.79,并创建了一个基于物理化学和生物土壤变量的MRA方程来预测这一属性,其中关键预测因子包括风化指数Ki、土壤有机质(SOM)的化学性质和质量、沙含量以及微生物酶活性。我们制作了巴西首张FCO2地图,该地图与实地观测数据的相关性达到了0.56,其空间变异性受矿物学和MAOC含量的影响。我们还首次开发了用于估算该国潜在碳饱和度的方程,并将其调整至SEM水平(R2=0.76),这证实了铁氧化物和铝氧化物在碳封存潜力中的重要性。通过NEP分析,我们识别出了巴西土壤储存和捕获碳的潜力。我们认为有必要评估潘帕(Pampa)和塞拉多(Cerrado)生物群落的土地利用情况,尽管这些地区具有较高的碳封存潜力,但由于铁氧化物含量超过了其储存能力,仍观察到碳负平衡。相比之下,潘塔纳尔(Pantanal)生物群落中的铁氧化物含量也远未达到其最大吸附能力,但该地区碳平衡为正,且排放量低、碳捕获潜力高。高生物量生产也有助于碳平衡为正,例如亚马逊地区,那里根系呼吸作用产生的高CO2排放量被高净初级生产力(NPP)所抵消。我们的方法和结果可以作为巴西农业生态系统碳捕获能力的一般估算,如果考虑到适当的土壤管理,这一估算可以得到改进。

**作者贡献**

Heidy Soledad Rodríguez-Albarracín:概念化;形式分析;方法论;软件;可视化;初稿撰写;审稿与编辑。
José A. M. Demattê:概念化;资金获取;资源协调;监督;验证;审稿与编辑。
Miguel Palacio Pelaez Carvalho Neto:方法论;软件;审稿与编辑。
Nícolas Augusto Rosin:方法论;软件;审稿与编辑。
Carlos Eduardo Pellegrino Cerri:资金获取;资源协调;审稿与编辑。
Carlos A. da Silva Junior:审稿与编辑。
José R. de Sousa Lima:审稿与编辑。
Eduardo S. de Souza:审稿与编辑。
Mara Regina Moitinho:审稿与编辑。
Aquiles Enrique Darghan Contreras:方法论;审稿与编辑。
Paulo Eduardo Teodoro:审稿与编辑。
Rafael Felippe Ratke:审稿与编辑。
Uemeson José dos Santos:审稿与编辑。
Denizard Oresca:审稿与编辑。

**致谢**

作者感谢ESALQ-USP的土壤科学系(LSO)在本研究期间提供的支持,以及土壤科学小组(GEOCIS)(https://esalqgeocis.wixsite.com/english)的成员们的帮助。本研究的文章发表费用由巴西高等教育人员培训协调委员会(CAPES)(ROR标识符:00x0ma614)资助。

**利益冲突声明**

作者声明没有利益冲突。

生物通微信公众号
微信
新浪微博


生物通 版权所有