中国西北荒漠区藻结皮覆盖下土壤有机碳垂直分布特征及其影响因素
1.
2.
3.
Vertical distribution characteristics and influencing factors of soil organic carbon under cyanobacterial crusts in the Northwest desert region of China
1.
2.
3.
通讯作者:
收稿日期: 2025-03-24 修回日期: 2025-04-24
基金资助: |
|
Received: 2025-03-24 Revised: 2025-04-24
作者简介 About authors
杜蕃(1998—),女,山东泰安人,硕士研究生,从事生物土壤结皮碳循环研究E-mail:
关键词:
Keywords:
本文引用格式
杜蕃, 杨军刚, 郭星, 陆永兴, 陶冶, 尹本丰, 荣晓莹, 李永刚, 张元明, 周晓兵.
Du Fan, Yang Jungang, Guo Xing, Lu Yongxing, Tao Ye, Yin Benfeng, Rong Xiaoying, Li Yonggang, Zhang Yuanming, Zhou Xiaobing.
0 引言
土壤有机碳(Soil Organic Carbon,SOC)是陆地碳库的核心组成部分,其储量分别约为大气和植被碳库储量的2倍和3倍[1],在调控陆地环境稳定和全球气候方面至关重要。尽管旱区SOC含量较低,但仍占全球土壤碳储量的27.9%[2]。其中,植物光合作用、动植物残骸以及土壤自养微生物的光合固碳,共同构成了旱地SOC的潜在来源[3-4]。生物土壤结皮作为旱地生态系统中最具代表性的生物成分,是主要由藻类、地衣和藓类等光合自养群落与异养细菌、真菌和古菌等生物类群与土壤颗粒紧密结合形成的土壤复合体[5]。生物土壤结皮约占全球旱区面积的30%[6],约占中国旱区面积的13.9%,主要分布在中国西北荒漠区,其中藻结皮约占中国旱地面积的10.7%[7]。藻结皮具有极强的环境适应性和生存能力,能够在极端干旱、贫瘠的环境中迅速定居和繁殖,为地衣和藓类结皮的形成提供条件。同时由于藻结皮具备一定的碳固定能力[8],使其成为干旱半干旱地区土壤的重要碳源,对于维持该区域生物多样性和生态系统功能具有不可替代的作用。为此,研究中国西北荒漠区藻结皮对SOC的影响,对深入理解土壤有机碳固存和碳循环具有重要的意义。
SOC的固存与土壤矿物间相互作用密切相关,通过形成物理和化学屏障调节SOC周转[9]。因此,通常将SOC划分为颗粒态有机碳(Particulate Organic Carbon,POC)和矿物结合态有机碳(Mineral-Associated Organic Carbon,MAOC),分别研究SOC的形成机制和稳定性[10]。其中,POC主要源自植物残体的分解过程,易被微生物利用,平均停留时间较短[11],代表不稳定易利用的土壤碳;MAOC的形成涉及微生物残体、分解产物和可溶性植物残留物吸附到土壤矿物表面[12-14],可指示相对稳定且周转周期较长的碳,因此POC/MAOC值一定程度上反映了SOC的稳定性[15]。由于POC和MAOC在形成和稳定性方面存在差异,其分布受不同因素驱动。研究表明气候是POC和MAOC储存的重要影响因素[16],气候变暖和降水量增加促使POC含量增加,而MAOC则呈下降或稳定趋势[17]。此外,土壤理化性质对土壤微生物的活跃程度以及有机质分解速率具有显著影响,同时土壤中氮磷元素可通过影响植物生长和微生物代谢调控SOC组分的输入量与分解速率[18]。因此,解析SOC的稳定性,需要综合考虑气候和土壤特征等影响因子。
国内外关于生物土壤结皮SOC周转方面已取得诸多进展。研究表明生物土壤结皮对荒漠SOC的积极影响强于草原[19];此外,在黄土高原北部小流域的研究表明,生物土壤结皮显著提高了SOC含量,并增加了表层土壤碳储量[20]。目前相关研究大多局限于特定的点位,对于区域生物土壤结皮覆盖对SOC组分含量及其驱动因子的研究较少。中国西北荒漠区具有特殊的光-温-湿的气候梯度分布格局,这些因子如何影响藻结皮覆盖下SOC含量及其垂直分布特征仍缺乏深入的理解。本研究调查了西北荒漠区具有生物土壤结皮分布的12个样点,并结合气候数据和土壤理化特征,分析了藻结皮覆盖下0~20 cm不同土层中SOC组分含量的分布特征及其影响因素,旨在揭示藻结皮在该区域SOC动态变化中的重要作用,为区域碳汇评估与管理提供科学依据。
1 材料与方法
1.1 研究区概况
中国西北荒漠区地处中纬度地带的欧亚大陆腹地,以山地和荒漠为主体,自西北向东南经过的主要沙漠为古尔班通古特沙漠、巴丹吉林沙漠、腾格里沙漠以及毛乌素沙地等。该地区降水稀少,日照强烈,气温变化强烈,其植被覆盖主要是以盐生旱生和超旱生以及耐盐碱的草本、灌木及小乔木为建群种构成的稀疏植被,同时生物土壤结皮广泛分布,藻结皮是西北荒漠区主要的结皮类型[21]。
1.2 样品采集和处理
2021年7月在西北荒漠区开展了大规模的采样工作,选择古尔班通古特沙漠、腾格里沙漠、库布其沙漠和毛乌素沙地中具有生物土壤结皮分布的12个研究样点(图1,表1)。样品采集遵循随机采样原则,首先在每一个样点周围随机选择地势平坦且藻结皮发育良好的区域划定5个30 m×30 m的大样方,并确保各样方之间保持至少50 m的间距以减少相互影响。在每个大样方中随机选择裸沙及藻结皮并根据对角线划定3个1 m×1 m的小样方,在每个小样方采集藻结皮层样品,同时对藻结皮下0~2、2~5、5~10、10~20 cm 4个土层采集土壤样品,将每个小样方中采集到的土壤样品混匀为一个样品。将采集的样品带回实验室,过2 mm筛去除植物根系和石块,随后将每份土壤样品分为2份,1份样品保存-20 ℃冰箱,另1份样品风干后保存。
图1
图1
中国西北荒漠区采样示意图
Fig.1
Sampling schematic diagram of the northwest desert region of China
表1 中国西北荒漠区采样点位详细信息
Table 1
点位 | 主要沙区 | 纬度 N/(°) | 经度 E/(°) | 年均气温 MAT/(℃) | 年均降水量 MAP/mm | 干旱指数 AI | 植被类型 |
---|---|---|---|---|---|---|---|
G1 | 古尔班通古特沙漠 | 45.30109 | 87.15127 | 7.56 | 162 | 0.1022 | 温带矮半乔木/灌木荒漠 |
G2 | 古尔班通古特沙漠 | 45.24653 | 87.67413 | 6.86 | 164 | 0.1069 | |
G3 | 古尔班通古特沙漠 | 45.13156 | 88.30210 | 6.11 | 185 | 0.1245 | |
G4 | 古尔班通古特沙漠 | 44.94260 | 88.65117 | 6.68 | 167 | 0.1085 | |
G5 | 库布其沙漠 | 40.38708 | 108.63548 | 7.06 | 236 | 0.1479 | 温带半灌木荒漠 |
G6 | 库布其沙漠 | 40.04151 | 108.67425 | 6.60 | 267 | 0.1724 | |
G7 | 腾格里沙漠 | 38.88572 | 102.72562 | 7.84 | 118 | 0.0704 | 温带半灌木、矮半灌木荒漠 |
G8 | 腾格里沙漠 | 37.96888 | 103.37158 | 8.25 | 159 | 0.0998 | |
G9 | 腾格里沙漠 | 38.26532 | 104.05925 | 8.68 | 146 | 0.0882 | |
G10 | 腾格里沙漠 | 38.52553 | 104.97848 | 8.97 | 164 | 0.0956 | |
G11 | 腾格里沙漠 | 38.78888 | 105.46609 | 8.74 | 168 | 0.0973 | |
G12 | 毛乌素沙地 | 38.86390 | 107.91095 | 7.37 | 292 | 0.1908 | 温带丛生禾草典型草原 |
1.3 土壤理化性质测定和气候因子获取
土壤pH和电导率(EC)采用电位法通过pH计和电导率仪测定(土壤与水混合比例分别为1∶2.5和1∶5);土壤全氮(TN)含量采用高氯酸-硫酸消化浸提液后,通过连续流动分析仪(SEA1 Auto Analyzer 3,德国)测定;土壤全磷(TP)含量通过高氯酸-硫酸消化法得到浸提液后,采用钼锑抗比色法进行测定;土壤速效磷(AP)含量采用碳酸氢钠浸提-钼锑抗比色法进行测定;土壤硝态氮(NO
年均气温(MAT)数据采用MODIS数据中的MOD11A全球气温格网数据(
1.4 SOC及其组分含量测定
SOC的测定采用盐酸-干烧法测定。具体方法为:称取0.2 g土壤风干样品于燃烧舟中,加入0.5 mol·L-1的盐酸溶液去除无机碳后置于烘箱中烘干,在此过程中加入蒸馏水去除盐酸残留,最后用碳氮分析仪(Analytik Jena Multi N/C 3100 TOC analyzer, 德国)测定。
POC和MAOC含量采用六偏磷酸钠分散法分离测定[22]。取10 g风干土壤样品放入100 mL塑料瓶中,加入60 mL 0.5%六偏磷酸钠溶液,置于摇床(90 r·min-1)振摇18 h分散样品,以破坏聚集体。将土壤悬液过53 µm筛,反复用蒸馏水漂洗,直至水流清澈。筛下的细质部分(˂53 µm)为矿物结合态有机质(MAOM),留在筛上的粗质部分(>53 µm)为颗粒态有机质(POM),将筛上下组分在60 ℃烘箱中烘48 h至恒重,称重后收集,将烘干土壤研磨,使用碳氮分析仪测定其碳含量。
土壤中POC和MAOC含量计算公式如下:
式中:M1和M2分别为上层(POM)和下层(MAOM)干重(g);POC测和MAOC测分别为测定的POC和MAOC含量,单位均为g·kg-1。
1.5 统计分析
每个样点取其均值作为一个重复,共有12个重复进行后续数据分析。采用双因素方差分析(Two-way ANOVA)探究不同覆盖类型、土层深度及其交互作用对SOC组分含量和土壤理化性质的影响。采用单因素方差分析(ANOVA)确定藻结皮和裸沙SOC组分含量和土壤理化性质在不同土层之间的差异;通过T检验进行藻结皮和裸沙SOC组分含量之间差异显著性检验,使用Origin 2021(Origin Lab, Hampton MA, USA)作图。为探究POC和MAOC含量的影响因素,利用R 4.2.3(R Development Core Team, 2023)软件的“LinKET”包对POC和MAOC含量与气候和土壤理化性质进行Person相关性分析,解析SOC组分与影响因子及各因子之间的关系,通过随机森林分析确定影响SOC组分含量的主要预测因子。
2 结果与分析
2.1 裸沙和藻结皮覆盖下土壤的基本性质及影响因素
通过对裸沙与藻结皮覆盖下土壤理化性质的分析,结果表明藻结皮的覆盖使得pH、EC和SWC的值高于裸沙,同时土壤养分TN和TP含量均高于裸沙,但裸沙与藻结皮之间差异并不显著(P>0.05)。随着土层深度的增加,藻结皮和裸沙TN、NO
表2 西北荒漠区裸沙和藻结皮覆盖下土壤理化性质
Table 2
类型 | 土层 /cm | 酸碱度 pH | 电导率 EC /(μS·cm-1) | 全氮 TN /(mg·kg-1) | 硝态氮 NO /(mg·kg-1) | 铵态氮 NH /(mg·kg-1) | 全磷 TP /(mg·kg-1) | 速效磷 AP /(mg·kg-1) | 土壤含水率SWC /% |
---|---|---|---|---|---|---|---|---|---|
裸沙 | 0~2 | 8.15±0.09a | 105.91±14.38a | 98.59±19.56a | 12.21±6.93a | 8.78±1.31a | 259.89±53.25a | 3.14±0.45a | 0.12±0.04b |
2~5 | 8.22±0.10a | 120.67±34.85a | 94.66±17.73a | 9.82±4.51a | 7.38±1.28a | 243.56±46.20a | 3.20±0.36a | 0.20±0.07b | |
5~10 | 8.27±0.09a | 90.91±8.53a | 95.14±18.01a | 11.28±7.06a | 6.80±1.21a | 237.19±46.82a | 2.84±0.26a | 0.31±0.09b | |
10~20 | 8.32±0.09a | 113.14±28.72a | 79.85±11.63a | 11.34±7.60a | 7.53±1.11a | 233.01±42.80a | 2.72±0.25a | 0.84±0.28a | |
藻结皮 | BSC | 7.72±0.09c | 194.91±37.15a | 202.71±30.34a | 10.90±2.79a | 11.43±1.50a | 365.15±43.00a | 3.37±0.40a | 0.18±0.03a |
0~2 | 8.14±0.08b | 161.83±50.30a | 147.62±23.26a | 12.11±4.41a | 6.56±1.41b | 282.09±39.65a | 3.18±0.36a | 0.19±0.02a | |
2~5 | 8.31±0.07ab | 165.78±65.19a | 134.15±23.42a | 12.66±5.73a | 6.15±1.25b | 269.44±41.53a | 3.20±0.33a | 0.34±0.11a | |
5~10 | 8.44±0.06a | 148.29±54.10a | 125.24±23.03a | 7.84±2.95a | 5.87±1.05b | 245.19±35.44a | 3.13±0.25a | 0.28±0.04a | |
10~20 | 8.51±0.06a | 136.24±36.39a | 124.15±24.78a | 6.87±2.35a | 6.98±1.25b | 263.09±45.28a | 2.88±0.27a | 0.35±0.08a |
双因素方差分析结果表明,覆盖类型对SOC、POC、MAOC、TN均具有显著影响,而对pH、EC、NO
表3 覆盖类型和土层深度对SOC组分以及土壤理化性质影响的双因素方差分析( F 值)
Table 3
指标 | 有机碳 SOC | 颗粒态有机碳 POC | 矿物结合态有机碳 MAOC | 酸碱度 pH | 电导率 EC | 全氮 TN | 硝态氮 NO | 铵态氮 NH | 全磷 TP | 速效磷 AP | 含水量 SWC |
---|---|---|---|---|---|---|---|---|---|---|---|
覆盖类型 | 16.41*** | 7.94** | 26.02*** | 3.63 | 2.52 | 6.93* | 0.12 | 1.88 | 0.48 | 0.27 | 1.00 |
土层深度 | 6.56*** | 2.25 | 11.99*** | 15.13*** | 0.31 | 2.28 | 0.12 | 3.46* | 1.16 | 0.59 | 4.33** |
覆盖类型×土层深度 | 0.38 | 0.17 | 0.58 | 0.66 | 0.08 | 0.07 | 0.20 | 0.16 | 0.02 | 0.08 | 3.21* |
2.2 裸沙和藻结皮覆盖下SOC组分的分布特征
图2
图2
裸沙和藻结皮覆盖下不同土层SOC、POC和MAOC含量分布特征
Fig.2
Distribution characteristics of SOC, POC, and MAOC content in different soil depths under bare sand and cyanobacterial crust
图3
图3
裸沙和藻结皮覆盖下POC和MAOC含量占比与稳定性
Fig.3
Proportions and stability of POC and MAOC under bare sand and cyanobacterial crust
2.3 环境因子与土壤POC和MAOC含量的相关性分析
裸沙与藻结皮覆盖下不同土层POC和MAOC含量与环境因子之间的相关性存在差异(图4)。在裸沙中,土壤养分(TN、NO
图4
图4
裸沙和藻结皮POC和MAOC与环境因子的相关关系
Fig.4
Correlation between environmental factors and POC and MAOC in bare sand and cyanobacterial crust
随机森林分析表明,MAT和MAP是预测POC含量的重要变量,其中MAT对藻结皮POC含量预测的贡献度为45.50%;MAP是预测MAOC含量的首要变量,影响裸沙MAOC含量的其次变量为TN,而影响藻结皮MAOC含量的其次变量则是NO
图5
图5
随机森林预测影响POC和MAOC含量的关键因子
Fig.5
Random forest model predicts the key factors affecting POC and MAOC content
3 讨论
3.1 中国西北荒漠区藻结皮覆盖对土壤碳含量的影响
随着土层深度的增加,裸沙和藻结皮的SOC、POC和MAOC含量均呈下降趋势,且同土层SOC组分含量表现为藻结皮高于裸沙。这与已有研究结果一致,在黄土高原,生物土壤结皮POC和MAOC含量随土层增加而降低[23]。对毛乌素沙地藻结皮研究发现,其覆盖下SOC含量为裸沙的2.1倍[24]。藻结皮发育在干旱区土壤发育和稳定中具有重要的生态功能,一方面,藻结皮作为生物土壤结皮初级演替阶段,其固碳作用主要通过原核生物和光合生物固碳途径完成[25],为有机质的累积提供了基础。另一方面,干旱区土壤普遍缺乏氮素,藻结皮中的固氮蓝藻在生物固氮中发挥重要作用,可为其他植物及土壤微生物供给充足的氮素资源[26],不仅促进了植物生长和土壤微生物活性,还增强了土壤中有机碳输入量,逐步增加生物土壤结皮的稳定性。
随着土层深度的增加,藻结皮层和0~5 cm土层SOC及其组分含量显著高于下层土壤,表明藻结皮层存在养分聚集效应。藻结皮层具有丰富的藻类,可通过光合作用合成有机物,加之藻结皮的覆盖,可以在一定程度上减少土壤表层养分流失,进而使得表层SOC组分含量高于下层[27]。同时,表层SOC增加又为土壤微生物生长提供适宜的条件,导致土壤微生物数量随土层深度增加呈现下降趋势。微生物是SOC形成与周转的主要参与者和调控者,在其增殖、生长直至死亡的生命周期中,会产生大量的细胞残留物及其分解过程中释放的小分子代谢产物,部分微生物残体与土壤矿物质发生紧密结合,促进MAOC的积累[10]。这一系列过程共同作用使得藻结皮在固定与储存SOC方面展现出显著的能力,成为干旱区中重要的碳汇。
裸沙与藻结皮中MAOC占SOC含量的比例呈现出随深度增加而下降的趋势。具体而言,从表层开始,裸沙MAOC含量占SOC的比例由16.85%降至12.69%,而藻结皮中这一比例则由42.70%降至28.95%。这些观测值显著低于先前研究中报道的比例,即SOC的60%以上被储存在MAOC库中[10]。西北荒漠区MAOC整体含量较低,一方面,荒漠区降水稀少使得土壤SWC处于较低水平,进而微生物活动受到限制,导致其对植物残体的分解速度相对较低,土壤有机质周转速度减慢,从而有利于POC含量的积累,导致MAOC占比降低。另一方面,荒漠区土壤中粗砂颗粒含量占据优势,其细砂颗粒含量相对较少,有机物质难与矿物表面充分接触和结合,不利于MAOC的形成[28]。由于POC易于被微生物分解的特性以及对环境变化的高度敏感性,通过POC/MAOC的值可以反映出SOC稳定性的相对大小[15]。本研究中,藻结皮SOC稳定性高于裸沙,主要原因是藻结皮能够增加土壤水分保持能力,为微生物生长和活动提供更适宜的条件,更有利于MAOC含量积累,进而提高了SOC的稳定性[29]。
3.2 藻结皮中POC和MAOC含量的影响因素
在区域尺度下,SOC的变化主要受气候、土壤和植被条件的影响[16]。本研究也表明,裸沙和藻结皮SOC组分变化的重要预测变量为MAP和MAT。多数研究认为,生物土壤结皮只有在外界温度和水分条件达到一定程度才能正常生长,其生长发育所需要的水分来自近地表土壤水分,引起微生物的快速响应[24,30],进而影响其SOC含量。在干旱区,降水对SOC的影响主要通过土壤微生物活性而非植被[31]。温度升高可能会增加微生物生长和周转,提高微生物代谢活性,从而会加速POC分解和损失,使得SOC含量减少。另外,变暖会引起土壤真菌相对丰度增加,促进MAOC分解,从而导致MAOC储量减少[15]。因此,二者与MAT呈负相关关系。小尺度上,一般来说,适宜的SWC有利于SOC积累,但是干旱区SWC比较低,SWC增加使得藻结皮增加植物碳输入,从而减轻微生物抑制,使得POC和MAOC含量增加[32]。因此,SWC与POC和MAOC含量呈正相关。
土壤理化特征变化会影响土壤微生物的活性以及植物凋落物的输入,这些变化会对SOC含量产生影响[33]。本研究结果表明,土壤养分(TN、NH
4 结论
在西北荒漠区,藻结皮覆盖不仅改变了土壤理化性质,还显著提升了SOC及其组分含量和SOC的稳定性。随着土层深度的增加,裸沙和藻结皮SOC、POC以及MAOC含量均呈现出逐渐降低的趋势。裸沙和藻结皮POC和MAOC含量与气候变量(MAT和MAP)和土壤理化性质之间存在显著相关性关系,同时MAP是预测POC和MAOC含量的重要变量。藻结皮的存在显著提高了西北荒漠区土壤有机碳的稳定性,有利于荒漠区碳固存。
参考文献
Soil carbon sequestration impacts on global climate change and food security
[J].
Environmental drivers of soil carbon and nitrogen accumulation in global drylands
[J].
Carbon sequestration in dryland ecosystems
[J].
The contentious nature of soil organic matter
[J].
The potential roles of biological soil crusts in dryland hydrologic cycles
[J].
Biocrust as one of multiple stable states in global drylands
[J].
Mapping biocrust distribution in China’s drylands under changing climate
[J].
Effects of biological soil crusts on soil labile organic carbon of patchy alpine meadows in the Source Zone of the Yellow River,West China
[J].
Dynamic interactions at the mineral-organic matter interface
[J].
Conceptualizing soil organic matter into particulate and mineral‐associated forms to address global change in the 21st century
[J].
SOM fractionation methods:relevance to functional pools and to stabilization mechanisms
[J].
Quantitative assessment of microbial necromass contribution to soil organic matter
[J].
Dual,differential isotope labeling shows the preferential movement of labile plant constituents into mineral-bonded soil organic matter
[J].
Direct evidence for microbial-derived soil organic matter formation and its ecophysiological controls
[J].
Soil organic carbon fractions in China: spatial distribution,drivers,and future changes
[J].
Distinct,direct and climate-mediated environmental controls on global particulate and mineral-associated organic carbon storage
[J].
Soil organic carbon response to global environmental change depends on its distribution between mineral-associated and particulate organic matter:a meta-analysis
[J].
Soil microbial carbon and nitrogen limitation constraints soil organic carbon stability in arid and semi-arid grasslands
[J].
Soil nitrogen and climate drive the positive effect of biological soil crusts on soil organic carbon sequestration in drylands:a meta-analysis
[J].
Contributions of three types of biocrusts to soil carbon stock and annual efflux in a small watershed of Northern Chinese Loess Plateau
[J].
Particulate Soil organic-matter changes across a grassland cultivation sequence
[J].
Biocrusts enhance soil organic carbon stability and regulate the fate of new-input carbon in semiarid desert ecosystems
[J].
Global turnover of soil mineral-associated and particulate organic carbon
[J].
Spatiotemporal dynamics of soil organic carbon in desert region and its response to climate change:a case study of Heihe river
[J].
Impacts of altered precipitation regimes on soil communities and biogeochemistry in arid and semi‐arid ecosystems
[J].
Nitrogen input alleviates the priming effects of biochar addition on soil organic carbon decomposition
[J].
Latitudinal patterns and drivers of plant lignin and microbial necromass accumulation in forest soils:disentangling microbial and abiotic controls
[J].
The amounts and ratio of nitrogen and phosphorus addition drive the rate of litter decomposition in a subtropical forest
[J].
Soil organic matter priming:the pH effects
[J].
/
〈 |
|
〉 |
