科尔沁沙地沙漠化过程中土壤有机碳含量变化特征及影响因素
1.
2.
Characteristics and influencing factors of soil organic carbon in the process of desertification in Horqin Sandy Land
1.
2.
通讯作者:
收稿日期: 2022-01-23 修回日期: 2022-03-17
基金资助: |
|
Received: 2022-01-23 Revised: 2022-03-17
作者简介 About authors
闫蒙(1997—),男,内蒙古呼和浩特人,硕士研究生,主要从事草地生态学研究E-mail:
关键词:
Keywords:
本文引用格式
闫蒙, 王旭洋, 周立业, 李玉强.
Yan Meng, Wang Xuyang, Zhou Liye, Li Yuqiang.
0 引言
土地沙漠化是人地关系不协调导致的在干旱、半干旱及部分半湿润地区以风沙活动为主要标志的土地退化过程[1],分布广泛、危害严重、威胁人口众多,受到了社会和学术界的广泛关注[2-3]。土壤碳库是评价土地退化程度的重要指数,由有机碳(SOC)和无机碳(SIC)库组成,其中SOC是全球碳循环的重要组成部分,调节着陆地生态系统的平衡[4-5]。SOC主要由动植物残体、根系分泌物等有机物质的分解形成,既是植物养分循环的焦点,又可以维持和巩固土壤结构,而且会影响土壤侵蚀潜力和水分分布[6]。据不完全统计,在全球范围内大约1 550 Pg碳以有机质的形式储存于土壤中,大约是植物碳储量的3倍,是大气碳储量的2倍[7-8]。土壤有机碳库拥有如此巨大的储量,使得其微小的变化都会引起全球气候的波动。因此,SOC库在全球碳平衡中具有十分重要的作用。
目前,关于沙漠化过程中SOC的变化特征及演变规律已有不少研究。刘树林等[9]对浑善达克沙地沙漠化过程中土壤粒度与养分变化关系研究表明,沙漠化土地表层土壤中SOC含量与土壤黏粒含量和粗粉沙含量之间存在很好的相关性,尤其与土壤黏粒含量的相关性非常高;马建业等[10]通过分析毛乌素沙地逆转过程中土壤颗粒固碳效应发现,土壤颗粒均表现出显著的固碳效应,且以砂粒和粉粒为主要固碳组分;赵哈林等[11]为了了解呼伦贝尔沙质草地沙漠化过程中土壤理化特性的变化规律,调查分析了不同类型沙漠化草地土壤理化特性的差异和变化特征,发现随着沙漠化的发展SOC含量明显下降,并且与土壤黏粉粒含量、土壤硬度显著正相关。然而,前人的研究结果大多建立在点尺度上,样地之间距离较小,无法考虑到不同沙漠化阶段样点之间SOC的空间异质性,缺乏从区域尺度揭示沙漠化过程中SOC的演变规律,因此有关沙漠化过程中SOC含量区域尺度的变化特征及其影响机制有待进一步研究。
科尔沁沙地位于半干旱农牧交错区东端,水热条件相对比较优越,近年来,由于气候变化及人类不合理地过度利用,该地区的顶极植被群落类型——疏林草原植被遭到严重破坏,造成了不同类型、不同程度的土地沙漠化[1],许多学者围绕科尔沁沙地沙漠化在土壤碳、氮分布[12]、植物-土壤系统碳氮储量动态[13]等方面进行了大量的研究,然而对于科尔沁沙地在沙漠化过程中SOC含量在区域尺度变化特征及其影响因素方面的研究较少。因此本研究从区域尺度出发,以科尔沁沙地不同沙漠化阶段的土壤为研究对象,通过野外调查取样,利用空间序列代替时间演替的方法,选取不同沙漠化阶段(疏林草地、固定沙地、半固定沙地、半流动沙地、流动沙地)的土壤,分析沙漠化过程对SOC含量的影响,探究SOC变化特征及其影响机制,以期为科尔沁沙地恢复治理提供科学依据。
1 材料与方法
1.1 研究区概况
研究区位于科尔沁沙地核心区(42.3°—44.5°N、118.5°—123.5°E,海拔180—650 m),是典型的半干旱区域,总面积约5.06万km2,属于大陆性半干旱、半湿润季风气候,年平均气温5.8—6.4 ℃,年降水量343—451 mm,其中70%集中于夏季,年蒸发量1 500—2 500 mm。年平均风速3.5—4.5 m·s-1,春季平均风速4.2—5.9 m·s-1,风速≥17 m·s-1的大风日数25—40 d,其中沙尘暴天气10—15 d,主要出现在春季。土壤类型以栗钙土、黑钙土和风沙土为主,地貌以平缓的流动沙丘、半流动沙丘、半固定沙丘、固定沙丘与丘间甸子地镶嵌分布为特征。植被以沙生植物为主,主要植物种包括黄柳(Salix gordejevii)、小叶锦鸡儿(Caragana microphylla)、差巴嘎蒿(Artemisia halodendron)、猪毛菜(Salsola collina)、白草(Pennisetum centrasiaticum)、马唐((Digitaria sanguinalis)等[14]。
1.2 研究方法
1.2.1 样地选择
沙漠化程度 | 裸沙面积占比/% | 植被盖度/% |
---|---|---|
轻度沙漠化 | <5 | >60 |
中度沙漠化 | 5—25 | 30—60 |
重度沙漠化 | 25—50 | 10—30 |
极重度沙漠化 | >50 | <10 |
图1
1.2.2 土壤样品采集
2021年8—9月,在科尔沁沙地不同沙漠化阶段的每个采样点,随机设置3个5 m×5 m的小样方,分别在5 m×5 m样方内按“五点取样法”选取5个钻点,用直径2.5 cm的土钻取样,分0—10、10—20 cm两层,每个小样方内5个点相同层次的土样均匀混合组成1个混合样,每一小样方获得2个混合样品(每土层1个),每一位点3个小样方计6个混合样。同时在每个样方中用环刀(100 cm3)分层(0—10、10—20 cm)取原状土,测定不同土层土壤的容重。
1.2.3 土壤样品测定及计算
土壤在自然状态下风干,过2 mm土壤筛,并去肉眼可见根和除枯枝落叶等杂物,测定土壤基本理化性质。土壤有机碳(SOC)用重铬酸钾-硫酸外加热法;全氮(TN)用半微量凯氏定氮法;全磷(TP)采用硫酸-高氯酸消煮,钼锑抗比色法;土壤容重(SBD)采用烘干法测定,土壤样品在105 ℃烘干至恒重,测定干土的质量[17]:
式中:SBD为土壤容重(g·cm-3);M1为铝盒的重量(g);M2为铝盒+烘干土质量(g);V为环刀容积(cm3)。
土壤有机碳(SOC)密度公式如下:
式中:SOCD是土壤有机碳密度(kg·m-2);SOCi 是第i层的土壤有机碳含量(g·kg-1);SBDi 是第i层的土壤容重(g·cm-3);Di 是第i层的土层厚度(m);m是土层的数量[18]。
表2 土壤粒级
Table 2
粒径/mm | 分类 |
---|---|
2.00—0.50 | 粗砂 |
0.50—0.25 | 中砂 |
0.25—0.10 | 细砂 |
0.10—0.05 | 极细砂 |
<0.05 | 黏粉粒 |
1.2.4 环境因子数据获取
环境因子包括气象因子、地形因子及植被指数,其中气象因子包括年降水量、年积温、生长季(5—9月)平均气温和生长季降水量,由中国科学院资源环境科学数据中心(
1.3 数据分析
采用SPSS26.0进行单因素方差分析(One-way ANOVA)和最小显著差异法(LSD)分析比较不同土层SOC在不同沙漠化阶段的差异(α=0.05),采用Pearson相关系数法确定各因子间相关关系,采用主成分分析法对相关土壤因子及环境因子进行分析,采用一元线性回归分析法分析SOC与土壤各主要理化性质之间的相关性。
2 结果与分析
2.1 沙漠化过程同一土层土壤有机碳含量特征
在沙漠化过程中,除流动沙地下层(10—20 cm)土壤中SOC含量略高于上层(0—10 cm)外,其他均为上层SOC含量高于下层,且上下两层土壤SOC含量均呈显著性差异(P<0.05,图2)。同一土层土壤SOC含量均随沙漠化程度的加剧而降低,0—10 cm土层中,不同沙漠化阶段有机碳含量差异显著(P<0.05);10—20 cm土层中,疏林草地SOC含量显著高于其他阶段(P<0.05),但半固定沙地与半流动沙地土壤SOC含量无显著差异(P>0.05)。0—10 cm土壤上层SOC含量占整个剖面(0—20 cm)SOC含量的47.7%—68.6%。此外,SOC含量在疏林草地0—10 cm土层土壤达到最大值,为8.42 g·kg-1,在流动沙地0—10 cm土层土壤降至最低,为0.42 g·kg-1。
图2
图2
同一土层深度土壤有机碳含量沿沙漠化程度分布特征
不同大写字母表示同一沙漠化阶段不同土层间差异显著(P<0.05),不同小写字母表示同一土层不同沙漠化阶段间差异显著(P<0.05)
Fig.2
Distribution characteristics of soil organic carbon content along desertification degree at the same soil layer depth
2.2 不同沙漠化阶段SOC密度
在科尔沁沙地沙漠化过程中,SOC密度随沙漠化程度的加剧呈下降趋势,且不同沙漠化阶段SOC密度差异显著(P<0.05,图3),其中顶极植被群落的疏林草地SOC密度最高,为2.15 kg·m-2,而发生极重度沙漠化的流动沙地SOC密度最低,仅为0.34 kg·m-2。
图3
图3
土壤有机碳(SOC)密度随沙漠化过程变化特征
不同小写字母表示在P<0.05水平上差异显著
Fig.3
Variation characteristics of soil organic carbon (SOC) density with desertification process
2.3 土壤因子相关指标及与土壤有机碳含量相关关系
除流动沙地外,同一沙漠化阶段土壤全氮含量(STNC)和土壤全磷含量(STPC)均为上层土壤(0—10 cm)高于下层土壤(10—20 cm),其中疏林草地0—10 cm土层的全氮含量和全磷含量均最高,分别为0.91 g·kg-1和0.20 g·kg-1,并且均随沙漠化程度的加深而降低,除半流动沙地和流动沙地外,其他均呈现显著性差异(P<0.05,表3)。土壤容重随沙漠化程度的加深而增大,但同一土层所有沙漠化阶段均无显著性差异(P>0.05)。固定沙地两个土层的pH值均最高;土壤电导率随着沙漠化程度的加深而降低,疏林草地和固定沙地两个土层土壤电导率显著高于其他沙漠化阶段(P<0.05);土壤粗砂粒和中砂粒含量均为半固定沙地与半流动沙地最高,而细砂粒为流动沙地最高且显著高于其他沙漠化阶段(P<0.05),土壤极细砂粒与黏粉粒随沙漠化程度的加深而减少,并且除半固定沙地与半流动沙地外,上层土壤极细砂粒与黏粉粒均高于下层土壤,流动沙地细砂粒显著高于其他沙漠化阶段(P<0.05)。
表3 不同沙漠化阶段不同土层土壤因子相关指标
Table 3
土壤因子 | 土层深度 /cm | 沙漠化阶段 | ||||
---|---|---|---|---|---|---|
疏林草地 | 固定沙地 | 半固定沙地 | 半流动沙地 | 流动沙地 | ||
土壤全氮含量/(g·kg-1) | 0—10 | 0.91±0.11a | 0.60±0.04b | 0.33±0.14c | 0.15±0.06d | 0.07±0.01d |
10—20 | 0.37±0.05a | 0.27±0.01b | 0.17±0.07c | 0.12±0.07d | 0.06±0.01d | |
土壤全磷含量/(g·kg-1) | 0—10 | 0.20±0.01a | 0.15±0.04b | 0.11±0.03c | 0.09±0.02d | 0.06±0.02d |
10—20 | 0.15±0.01a | 0.10±0.03b | 0.08±0.02bc | 0.07±0.02c | 0.06±0.01c | |
土壤容重/(g·cm-3) | 0—10 | 1.61±0.06a | 1.63±0.06a | 1.64±0.07a | 1.63±0.06a | 1.64±0.04a |
10—20 | 1.64±0.06a | 1.66±0.04a | 1.66±0.07a | 1.67±0.06a | 1.67±0.03a | |
pH值 | 0—10 | 7.02±0.69ab | 7.21±0.69ab | 6.80±0.61ab | 6.72±0.54b | 6.90±0.66ab |
10—20 | 7.07±0.47b | 7.59±0.84a | 7.09±0.67b | 6.88±0.61b | 6.90±0.71b | |
电导率/(μS·cm-1) | 0—10 | 42.38±4.71a | 39.17±3.87a | 27.18±4.20b | 19.55±3.88c | 13.82±1.32d |
10—20 | 36.20±2.50a | 37.41±5.44a | 25.19±5.53b | 21.70±3.79bc | 17.20±5.03c | |
粗砂/% | 0—10 | 1.29±0.33b | 1.39±0.40b | 2.10±0.40a | 1.90±0.42a | 0.58±0.36c |
10—20 | 1.27±0.46ab | 1.18±0.44ab | 2.01±0.31a | 2.18±0.93a | 0.75±0.12b | |
中砂/% | 0—10 | 14.36±3.57b | 18.02±3.23a | 20.74±2.86a | 19.62±2.85a | 9.54±2.08c |
10—20 | 16.23±3.87b | 17.48±1.33b | 23.82±5.31a | 19.26±4.27b | 9.25±1.13c | |
细砂/% | 0—10 | 45.43±7.38e | 51.70±7.83d | 57.40±8.07c | 64.60±7.99b | 72.61±8.87a |
10—20 | 54.02±2.54b | 56.07±5.55b | 59.02±6.44b | 60.40±6.01b | 74.70±7.22a | |
极细砂/% | 0—10 | 16.97±3.01a | 13.47±3.01b | 9.52±2.81c | 6.62±2.18d | 7.28±2.26d |
10—20 | 14.49±3.72a | 11.27±2.51b | 9.54±1.38b | 7.64±1.36b | 7.31±1.90b | |
黏粉粒/% | 0—10 | 21.20±3.52a | 17.11±4.01b | 11.69±3.45c | 6.79±2.53d | 9.43±2.44c |
10—20 | 14.89±3.90a | 12.04±3.34b | 8.20±2.61c | 5.97±2.20c | 6.02±2.14c |
不同小写字母表示在P<0.05水平上差异显著。
不同土层SOC含量与全氮含量、全磷含量均呈极显著相关关系(P<0.01,表4),土壤电导率仅与0—10 cm土层SOC含量呈极显著负相关关系(P<0.01);不同土层SOC含量与土壤pH值无显著性相关关系(P>0.05);0—10 cm土层SOC含量与电导率呈极显著相关关系(P<0.01);土壤容重仅与10—20 cm土层SOC含量存在极显著相关关系(P<0.01);土壤细砂、极细砂及黏粉粒仅与0—10 cm上层土壤SOC含量呈显著负相关关系(P<0.05)。这说明科尔沁沙地不同沙漠化阶段0—20 cm土层SOC含量的改变依赖于全氮、全磷含量的变化,以及由于沙漠化程度的加重导致的土壤机械组成变化。
表4 不同土层土壤有机碳与土壤因子相关性
Table 4
土层深度/cm | 土壤全氮含量 | 土壤全磷含量 | 土壤容重 | pH值 | 电导率 | 粗砂 | 中砂 | 细砂 | 极细砂 | 黏粉粒 |
---|---|---|---|---|---|---|---|---|---|---|
0—10 | 0.883** | 0.784** | -0.046 | 0.02 | 0.511** | -0.020 | -0.089 | -0.389** | 0.320** | 0.586** |
10—20 | 0.503** | 0.411** | -0.097** | -0.08 | 0.051 | -0.053 | -0.102 | -0.045 | 0.102 | 0.147 |
**P<0.01。
2.4 影响土壤有机碳含量的土壤因子
本研究通过主成分分析提取影响科尔沁沙地不同沙漠化阶段土壤各层SOC的主要土壤因子(表5)。结果显示,0—10 cm土层中,影响科尔沁沙地不同沙漠化SOC含量土壤因子的第1、2、3主成分累计贡献率达到72.835%;10—20 cm土层中,第1、2、3主成分累计贡献率达到67.193%,能反映土壤因子对SOC影响效应的大部分信息。0—10 cm土层第1主成分主要受土壤全氮、全磷、电导率、细砂、极细砂以及黏粉粒的影响,累计解释变量为42.148%,第2主成分主要受粗砂和中砂的影响,累计解释变量为62.223%,第3主成分主要受pH值的影响,累计解释变量为72.835%;10—20 cm土层第1主成分与土壤全氮、全磷、极细砂有较大的相关性,累计解释变量为34.755%,第2主成分主要受中砂和细砂的影响,累计解释变量为54.515%。
表5 不同土层土壤有机碳与土壤因子主成分( F )分析
Table 5
变量 | 0—10 cm | 10—20 cm | ||||
---|---|---|---|---|---|---|
F1 | F2 | F3 | F1 | F2 | F3 | |
土壤全氮含量 | 0.853 | -0.166 | -0.299 | 0.807 | -0.254 | -0.142 |
土壤全磷含量 | 0.897 | -0.110 | -0.261 | 0.801 | -0.105 | -0.226 |
土壤容重 | -0.107 | 0.431 | 0.032 | -0.059 | 0.287 | 0.547 |
pH值 | 0.510 | 0.089 | 0.816 | 0.607 | 0.004 | 0.577 |
电导率 | 0.791 | -0.035 | 0.340 | 0.632 | -0.211 | 0.535 |
粗砂 | 0.136 | 0.814 | -0.125 | 0.150 | 0.754 | 0.113 |
中砂 | 0.013 | 0.865 | -0.052 | -0.100 | 0.840 | -0.040 |
细砂 | -0.777 | -0.533 | 0.117 | -0.628 | -0.705 | 0.237 |
极细砂 | 0.732 | -0.072 | 0.158 | 0.755 | -0.025 | 0.074 |
黏粉粒 | 0.793 | -0.271 | -0.251 | 0.642 | -0.027 | -0.450 |
特征值 | 4.215 | 2.008 | 1.061 | 3.475 | 1.976 | 1.268 |
解释变量/% | 42.148 | 20.075 | 10.612 | 34.755 | 19.761 | 12.678 |
累计解释变量/% | 42.148 | 62.223 | 72.835 | 34.755 | 54.515 | 67.193 |
2.5 不同土层土壤有机碳含量一元线性回归模型
基于单变量一元回归分析,可得出不同土层主要土壤影响因子与SOC含量的一元线性回归方程(表6)。两个垂直土层SOC含量均与全氮、全磷呈极显著正相关关系(P<0.01),能解释0—10 cm土层中SOC含量的77.8%的变异,其中电导率与SOC含量拟合R2 最大(0.995),除此之外,10—20 cm土层SOC含量与土壤容重亦呈极显著相关关系(P<0.01),其中全氮与SOC含量拟合R2最大(0.363)。
表6 不同土层土壤有机碳含量与土壤因子的一元线性回归方程
Table 6
土层深度/cm | 回归方程 | R2 | F | P |
---|---|---|---|---|
0—10 | y=8.595STNC+0.118 | 0.778 | 484.174 | <0.01 |
y=45.435STPC-0.074 | 0.612 | 218.398 | ||
y=0.011EC+0.126 | 0.995 | 48.399 | ||
y=-0.088Fsa+0.954 | 0.145 | 24.465 | ||
y=0.162Vfs+0.273 | 0.096 | 15.647 | ||
y=0.255Scl+0.130 | 0.339 | 71.625 | ||
10—20 | y=15.193STNC+0.154 | 0.363 | 79.598 | |
y=37.905STPC+0.121 | 0.239 | 44.452 | ||
y=-0.147SBD+0.692 | -0.002 | 0.659 |
y,土壤有机碳;STNC,土壤全氮含量;STPC,土壤全磷含量;EC,电导率;Fsa,细砂粒;Vfs,极细砂粒;Scl,黏粉粒;SBD,土壤容重。
2.6 影响土壤有机碳含量的环境因子
通过主成分分析提取影响科尔沁沙地不同沙漠化阶段土壤各层SOC的主要环境因子(表7)以从区域尺度分析影响科尔沁沙地SOC的环境因素。结果显示,0—10 cm、10—20 cm 土层主成分中大于1的特征值成分有4个,累计解释变量分别为78.743%、78.603%,能反映不同沙漠化土层 SOC的绝大部分信息。两个土层的第1主成分均主要受经度、高程、生长季平均气温以及NDVI均值影响,累计解释变量分别为36.112%和36.347%;第2主成分均主要受地面剖面曲率、地面粗糙度与坡度的影响,累计解释变量分别55.343%和55.721%;第3主成分均与纬度和生长季降水量的相关性最大,累计解释变量分别达69.579%和69.597%;第4主成分与年积温的相关性最大。
表7 不同土层土壤有机碳与环境因子主成分( F )分析
Table 7
变量 | 0—10 cm | 10—20 cm | ||||||
---|---|---|---|---|---|---|---|---|
F1 | F2 | F3 | F4 | F1 | F2 | F3 | F4 | |
纬度 | 0.448 | -0.048 | -0.791 | 0.314 | 0.442 | -0.053 | -0.823 | 0.188 |
经度 | 0.921 | 0.024 | 0.080 | 0.123 | 0.920 | 0.049 | 0.060 | 0.108 |
年降水量 | 0.645 | 0.236 | 0.269 | -0.168 | 0.637 | 0.257 | 0.277 | -0.106 |
年积温 | 0.480 | 0.156 | 0.367 | 0.469 | 0.496 | 0.172 | 0.309 | 0.453 |
高程 | -0.926 | 0.017 | 0.003 | -0.280 | -0.933 | -0.005 | 0.041 | -0.238 |
地面剖面曲率 | -0.073 | 0.802 | -0.102 | 0.029 | -0.091 | 0.798 | -0.128 | 0.009 |
地面粗糙度 | -0.310 | 0.883 | 0.059 | 0.206 | -0.323 | 0.872 | 0.026 | 0.255 |
坡度 | -0.398 | 0.831 | 0.029 | 0.162 | -0.413 | 0.817 | 0.007 | 0.223 |
生长季平均气温 | 0.720 | -0.053 | 0.130 | 0.427 | 0.739 | -0.037 | 0.075 | 0.380 |
生长季降水量 | 0.192 | -0.111 | 0.917 | -0.115 | 0.204 | -0.084 | 0.928 | -0.030 |
生长季EVI均值 | 0.674 | 0.391 | -0.040 | -0.482 | 0.649 | 0.414 | -0.024 | -0.554 |
生长季NDVI均值 | 0.737 | 0.367 | -0.118 | -0.454 | 0.712 | 0.389 | -0.101 | -0.520 |
特征值 | 4.695 | 2.500 | 1.851 | 1.191 | 4.729 | 2.515 | 1.804 | 1.171 |
解释变量/% | 36.112 | 19.232 | 14.236 | 9.163 | 36.374 | 19.347 | 13.876 | 9.006 |
累计解释变量/% | 36.112 | 55.343 | 69.579 | 78.743 | 36.374 | 55.721 | 69.597 | 78.603 |
3 讨论
3.1 不同沙漠化阶段土壤有机碳剖面分布特征及空间分布规律
SOC是土壤养分循环的核心,也是衡量土壤质量变化和肥力变化的重要指标[20],本研究表明,科尔沁沙地同一沙漠化阶段SOC含量垂直分布特征主要表现为除流动沙地外,土壤上层(0—10 cm)SOC含量明显高于下层(10—20 cm),呈现明显的“表聚性”特征[21],这与常帅等[22]对锡林郭勒草原土壤剖面SOC垂直分布及丁越岿等[23]就不同植被类型对毛乌素沙地SOC的影响研究结果一致。一方面可能是因为科尔沁沙地植物根系主要分布在土壤表层,同时土壤表层积累了大量的植物掉落物以及动物残体,而土壤表层SOC主要来源于植物根系分泌物、植物根系的腐殖质及枯枝落叶[21],再加上土壤表层有适合微生物生存的良好的通气性、含水量以及土壤温度等环境因素,土壤微生物活动旺盛[24-25],加速了地表掉落物及动植物残体的分解及腐殖化,使得土壤中的有机质含量增加,土壤容重降低,SOC含量不断累积;另一方面,随着土层深度的加深,土壤质地变得紧实,表层土壤的SOC向下迁移受到限制,使得下层土壤的SOC含量低于上层,但不同退化程度间SOC的降低幅度存在较大差异,由于流动沙地地表无结皮或只有物理性结皮,并且地表完全被流沙覆盖[26],因此流动沙地上下两层土壤理化性质差异较小,SOC含量无显著差异。除此之外,其他沙漠化阶段SOC含量垂直分布差异可能与不同沙漠化阶段土壤有机质及地下根系的生物量分布差异有关,土壤表层中植物根系较为富集,能较好地吸收来自土壤微生物分解地表掉落物释放的养分,同时在疏林草地中,丰富的表层根系还可以截获水分,促进植被生长状况的正向循环,进一步增加凋落物及根系向土壤的营养输入[27-28],从而对不同沙漠化阶段土壤剖面SOC含量造成影响。
3.2 不同沙漠化阶段土壤有机碳储量特征
土壤有机碳密度(SOCD)是指单位面积一定深度的土层中SOC的储量,是评价和衡量土壤中SOC储量的一个极其重要的指标[29]。本研究结果表明,在科尔沁沙地,随着沙漠化程度的增加SOCD呈现出降低的趋势,与阎欣等[30]就荒漠草原沙漠化对土壤物理和化学特性的影响的研究结果一致。科尔沁沙地沙漠化过程中,不同沙漠化阶段对应着适合本生境的典型群落类型[31],其顶级群落疏林草地植物丰富度高,植被盖度大于90%,但随着沙漠化的加剧,原生地带性植物逐渐消退,植物群落组成结构趋于简单,植被盖度下降至0—10%[32]。首先导致地表凋落物减少,有机质无法滞留,从而直接降低了SOC的输入[33],同时,土壤微生物作为土壤的重要组成部分,推动着土壤物质循环及能量流动,对SOC的形成与分解有着重要影响[34],而SOC又可以作为底物,为微生物的生长和繁殖提供碳源[35],沙漠化的加剧不利于土壤微生物的生长和繁殖,从而使SOC的积累能力减弱,导致SOCD降低。此外,土地沙漠化导致土壤的理化性质也发生变化(黏粒含量、土壤容重、土壤pH值等因素),而这些因素对SOCD的影响同样不可忽视。
3.3 不同沙漠化阶段科尔沁沙地土壤有机碳影响因素
SOC是当前全球变化和全球碳循环的研究热点,也是退化生态系统恢复重建研究的重要内容[36],在全球变化的背景下,研究科尔沁沙地不同沙漠化阶段SOC含量并揭示其影响因素,有助于增加土壤碳截存,提高土壤质量,并且对沙漠化土地的生态恢复及环境治理具有重要的意义[37]。在本研究中,虽然科尔沁沙地不同沙漠化阶段SOC含量与土壤理化性质随土层深度的增加存在一定的差异性,但不同沙漠化阶段垂直剖面SOC含量与STNC、STPC均呈极显著相关关系(P<0.01),这与相关学者的研究结果一致[38-39]。原因可能是沙漠化程度越低,地表植被覆盖度越高,局部区域风沙活动越弱,大气降尘沉积积累,同时生物小循环过程旺盛,地表枯落物在腐烂分解的过程中土壤有机质含量增加[40-41],使得STNC、STPC及SOC含量增加。
土壤容重是土壤基本物理性质,可以反映土壤的通气、透水以及松紧程度,也是评估SOC贮量的重要参数[42-43]。本研究发现,在科尔沁沙地沙漠化过程中,上层土壤(0—10 cm)SOC含量与土壤容重无显著相关关系,但下层土壤(10—20 cm)SOC含量与土壤容重呈极显著负相关关系(P<0.01),即下层土壤容重越大,SOC含量越少。这可能是由于随着土层深度的增加,土壤通气性变差,导致土壤容的增加,而土壤容重的增加会使土壤微生物因氧气供给受限而活动性减弱[23,44],从而造成土壤中有机质降低,SOC含量减少。而本研究中土壤容重仅与下层土壤(10—20 cm)显著相关,很可能是因为科尔沁沙地表层环境较适宜,微生物活性较强,土壤有机碳与土壤全氮及全磷的极显著相关关系可能掩盖了土壤容重对土壤有机碳累积的影响。
土壤机械组成代表了土壤质地的粗细,可以直接影响土壤的理化性质,是土壤最基本的物理性质[45],且土壤机械组成比例能通过其化学保护机制降低SOC的分解速率并且促进土壤有机碳的固持[46]。在本研究中,科尔沁沙地土壤砂粒含量随沙漠化程度的加剧而增加,与SOC分布变化相反,土壤极细砂和黏粉粒含量则均随沙漠化程度的加剧而减少,与SOC含量变化特征相同,这与有关学者的研究一致[47-49]。表4中,0—10 cm土层SOC含量与土壤细砂粒呈极显著性负相关关系,与土壤极细砂和黏粉粒呈极显著正相关。在不同土层土壤粒径组成中,黏粉粒含量对SOC的影响最强,可作为衡量SOC含量的重要因素,这可能是因为黏粉粒颗粒较细,表面积较大,能暴露更多的正电荷位与土壤中带负电荷的腐殖质结合[50],使得SOC更易与土壤颗粒中的粉粒结合,对土壤水分的有效性及植被生长产生影响,从而影响SOC分布[51]。
3.4 区域尺度科尔沁沙地土壤有机碳影响因素
气候和地形因子很大程度上决定着陆地植被生产力以及生物群落组成,而经纬度的不同会直接或间接地决定科尔沁沙地的降水、气温等气象因子以及土壤质地的空间异质性,土壤质地的不同也会对土壤的水、肥、 气、热的分配与循环有不同程度的作用[52],从而影响植被生长与SOC的矿化,进而对生态系统SOC的空间格局进行调控。主成分分析结果表明,经纬度、气候(生长季平均气温、年降水量、年积温)、地形(地面剖面斜率、地面粗糙程度、坡度、高程)以及植被覆盖状况均对科尔沁沙地SOC含量空间格局有不同程度的影响,其中经纬度和地形是最重要的影响因素。在科尔沁沙地区域尺度上,气温和降水均随经纬度的变化具有较为明显的空间分布特征,气温和降水决定着不同地理位置的水热条件,科尔沁沙地植被空间分布也因此不同[53]。本研究结果显示,降水对科尔沁沙地SOC的影响较气温更为显著,这与李妙宇等[54]的研究结果相似,降水作为干旱半干旱地区唯一土壤水分来源,限制着植被的生产力,较高的降水量会增加植被生产力,从而使地上植被和地下根系生物量增加,土壤微生物活性也随之提高,有机质分解速率也因此加快[55]。
除气候因子外,地形因子也显著影响着科尔沁沙地植被及SOC含量,其中高程直接或间接影响着生态系统SOC含量,随着高程的增加,气温也会随之降低,限制植被生长,从而限制SOC的分解,有利于其积累[56]。坡度代表着地面倾斜程度,是重要的地形指标。本研究结果显示,科尔沁沙地SOC含量同样受坡度影响。
4 结论
科尔沁沙地SOC含量0.42—8.42 g·kg-1,SOCD为0.34—2.15 kg·m-2,其中SOC主要分布在0—10 cm的上层土壤中;SOC含量从疏林草地(未沙漠化)、固定沙地(轻度沙漠化)、半固定沙地(中度沙漠化)、半流动沙地(重度沙漠化)到流动沙地(极重度沙漠化)逐渐降低,表现出按疏林草地—荒漠草地—流动沙地演替退化的过程,表层SOC含量流失更加明显;在科尔沁沙地沙漠化过程中,土壤机械组成和全氮、全磷含量限制着SOC的积累,而土壤黏粉粒含量更能反映土壤机械组成对SOC含量的限制作用;区域尺度中,经纬度是影响科尔沁沙地SOC含量的重要因素,其中气候因子是SOC经纬度分布格局的主控因素,未来应针对不同环境及不同地区SOC含量采取不同的区域管理方式,以提高科尔沁沙地生态系统的碳储量,也可以通过围封、灌溉、补播等适宜的措施,改变沙漠化草地的土壤结构,维持沙化草地生态系统的自我调节能力,增加其土壤碳库储量及其稳定性。
参考文献
Ecosystems and human well-bein:desertification synthesis
[R].
Revisiting no-till's impact on soil organic carbon storage in Canada
[J].
Can N2O emissions offset the benefits from soil organic carbon storage?
[J].
Evidence from chronosequence studies for a low carbon-storage potential of soils
[J].
Carbon pools and flux of global forest ecosystems
[J].
Dynamics of aeolian desertification and its driving forces in the Horqin Sandy Land,Northern China
[J].
Comparison of soil physico-chemical properties under different land-use and cover types in northeastern China's Horqin Sandy Land
[J].
A global,0.05-degree product of solar-induced chlorophyll fluorescence derived from OCO-2,MODIS,and reanalysis data
[J].
Effects of drying-rewetting frequency on soil carbon and nitrogen transformations
[J].
Estimating soil organic carbon storage and distribution in a catchment of Loess Plateau,China
[J].
Shrub encroachment in North American grasslands: shifts in growth form dominance rapidly alters control of ecosystem carbon inputs
[J].
Prediction of bulk density of soils in the Loess Plateau region of China
[J].
Predicting soil bulk density for incomplete databases
[J].
Stabilization mechanisms of soil organic matter: implications for C-saturation of soils
[J].
祁连山中段土壤有机碳剖面垂直分布特征及其影响因素
[J].
Pedogenesis,permafrost,and soil moisture as controlling factors for soil nitrogen and carbon contents across the Tibetan Plateau
[J].
How ecological restoration alters ecosystem services: an analysis of carbon sequestration in China's Loess Plateau
[J].
Storage,patterns and controls of soil organic carbon in the alpine shrubland in the Three Rivers Source Region on the Qinghai-Tibetan Plateau
[J].
Topographic influences on soil properties and aboveground biomass in lucerne-rich vegetation in a semi-arid environment
[J].
/
〈 |
|
〉 |
