-
土壤有机碳是土壤中地上地下凋落物、根系沉积物、动物残体和微生物体及其分解和合成的有机物质中的碳,是评价土壤肥力和土壤质量的重要指标[1]。由于土壤有机碳库含量巨大,其微小的变化将对全球碳循环和气候变化产生显著影响[2-3]。土壤有机碳储量受林分类型、地形地貌、气候变化和土地利用方式等自然和人为因素的影响,但土壤总有机碳对这些环境因素的响应表现出一定的滞后性,难以在短期内观测其变化[4-5]。土壤活性有机碳组分是土壤中活性较高、易分解矿化的有机碳,对区域环境变化的响应更敏感[6],可以及时、准确地反映土壤有机碳对干扰和环境变化的响应[7]。碳库管理指数是陆地生态系统中土壤碳动态变化率的敏感指标,可以衡量环境因素对土壤有机碳数量和质量的影响[8-9]。因此,明确有机碳及其活性有机碳组分和碳库管理指数等有机碳库特征对预测土壤碳循环具有重要意义。
不同林分类型的根系分布、凋落物组成和微生物群落结构等存在差异,导致光照、水分和土壤理化性质等环境条件变化,进而影响土壤有机碳积累和转化速率[10-11]。Wang等[2]认为,黄土高原延河流域森林较高的凋落物输入,增加土壤氮磷含量,降低土壤容重和pH值,促进有机碳及活性组分的积累和固存。Yao等[12]对昆士兰赫伯特河流域土壤研究发现,不同植被的凋落物输入、酶活性以及微生物组成存在差异,增加潮汐森林0~30 cm土层有机碳含量。Yuan等[13]对亚热带冲击平原土壤有机碳库研究认为,相较于香樟(Cinnamomum camphora (L.) Presl)和水杉(Metasequoia glyptostroboides Hu & W. C. Cheng),落叶阔叶重阳木(Bischofia polycarpa (Levl.) Airy Shaw)人工林较低的碳氮比和较短的细根寿命加速凋落物分解,缓解土地利用方式转变初期有机碳的损失。Dawud等[14]和Ratcliffe等[15]指出,林分类型对土壤有机碳的影响在不同气候和土壤类型下存在差异。
江苏省总面积约10.72万km2,其中,平原沙土区面积约2万km2,主要分为黄河故道沙土区、沿海平原沙土区和通南高沙土区。江苏省丰县属典型黄河故道沙土区,土壤类型以潮土和风沙土为主,土壤风化程度高,结构疏松,抗侵蚀能力弱,自然肥力低,孔隙率高,保水性差[16-17]。以往对平原沙土区研究多集中于植被恢复措施的水土保持效果[18],较少关注土壤碳库特征和土壤质量提升。前人研究表明,植被恢复可以改变沙土有机碳储量[19],但造林初期不同林分类型下沙土区土壤有机碳库特征及其与环境因子的关系尚不清楚。基于此,本研究以江苏省平原沙土区丰县16年生杨树(Populus × euramevicana cl. 'Zhonglin 46')纯林、柳树(Salix × jiangsuensis cl. ‘J-172’)纯林、杨柳混交林3种林分类型及相邻的撂荒地(对照)为对象,测定土壤有机碳及活性组分、蔗糖酶、过氧化物酶、多酚氧化酶和β-葡萄糖苷酶活性,并计算土壤碳库管理指数,旨在明晰不同林分类型下土壤有机碳库特征及其影响因子,以期为平原沙土区土壤碳库的提高和生态恢复提供理论依据。
-
试验地为2006年河道整治形成的立地条件相近的防护林,宽约40 m,长约6 km。沿河道依次营造长约2 km林龄相同的杨树纯林、柳树纯林和杨柳混交林。杨树品种为“46-杨”,柳树品种为“苏柳172”,林龄16 a,初始营造密度均为500 株· hm−2。杨树纯林(PPP)、柳树纯林(WPP)、杨柳混交林(PWM)下皆有天然更新的灌木和草本,其中,灌木主要有枸杞(Lycium chinense Miller)、苎麻(Boehmeria nivea (L.) Gaudich.)等;草本植物主要有狗尾草(Setaria viridis (L.) Beauv.)、沼生水马齿(Callitriche palustris L.)等。试验地基本概况见表1。
林分类型
Stand types平均胸径
Mean DBH/cm平均树高
Mean height/m郁闭度
Canopy
density凋落物厚度
Litter
thickness/
cm细根生物量
Fine root
biomass/
(t·hm−2)细根碳储量
Fine root
carbon stock/
(t·hm−2)凋落物生物量
Litter
biomass/
(t·hm−2)凋落物碳储量
Litter
carbon stock/
(t·hm−2)杨树纯林 PPP 28.8±1.6 16.8±0.6 0.9 2.2 8.69±0.55 A 3.93±0.31 A 12.89±0.01 A 5.49±0.19 A 柳树纯林 WPP 14.9±2.3 7.2±0.4 0.4 1.3 4.66±0.63 B 2.09±0.35 B 8.85±2.61 B 3.86±0.61 B 杨柳混交林 PWM 28.67±2.67(杨树)/
18.43±1.72(柳树)17.84±0.41(杨树)/
8.31±1.01(柳树)0.7 1.5 4.12±0.74 B 1.69±0.35 B 9.99±0.76 B 4.35±0.19 AB 撂荒地 CK – – – 0.2 1.14±0.04 C 0.45±0.02 C 2.86±1.38 C 1.21±0.35 C 注:表中的数据为“平均值±标准误”。不同大写字母表示不同林分之间差异显著(P<0.05)。下同。
Notes: The data is “mean±standard error”. Capital case letters that differ indicate significant (P< 0.05) differences among stand types. The same below.Table 1. Basic information of the sample plots
采样时间为2022年3月。在每种林分类型分别设置3个间距大于300 m的面积为20 m × 20 m的标准样地,根据试验地踏勘情况,选择与3种林分类型相邻的以狗尾草等草本植物为主的立地条件相同的撂荒地为对照(CK),撂荒地同样设置3个远离林分边界的面积为20 m × 20 m的标准样地。在每个标准样地中按“S”型选取5个采样点,去除地表凋落物后挖取土壤剖面,每个剖面大小为1 m × 1.5 m × 1.2 m,分0~20、20~40、40~60、60~80、80~100 cm 5个土层分别采集土壤,将5个采样点同一土层的土壤混合成该标准样地的土壤样品。用镊子仔细去除土壤中的根系,过2 mm筛备用。将每个标准样地的土壤样品分为2份,一份作为鲜土保存在4 ℃冰箱中,用于土壤可溶性有机碳(DOC)、微生物生物量碳(MBC)和pH值的测定;另一份经自然风干过筛等处理后用于全碳(TC)、有机碳(SOC)、易氧化有机碳(EOC)、蔗糖酶、过氧化物酶、多酚氧化酶和β-葡萄糖苷酶及其它土壤理化性质的测定。使用内径9 cm、高20 cm的根钻,在5个采样点分5个土壤层次分别采集1份土芯,混合为一份装袋,带回实验室用0.2 mm细筛淘洗,筛选出直径<2 mm的细根。在每个标准样地中随机选取3个1 m × 1 m的小样方,收集各小样方地表凋落物混合为该标准样地的凋落物样。
-
试验地土壤偏碱性,土壤、凋落物和细根的全碳和全氮含量采用元素分析仪(Vario ElementⅢ,Germany)测定,SOC采用重铬酸钾-外加热法测定;DOC经K2SO4浸提,采用TOC仪(TOC-VCPH + TNM-1, Shimazu Inc., Japan)测定;MBC用氯仿熏蒸-K2SO4提取法测定;EOC采用333 mmol·L−1高锰酸钾氧化法测定。土壤团聚体平均重量直径参考Elliot[21]湿筛法测定,其它土壤理化指标参照鲍士旦[22]的分析方法,容重和含水率用环刀法测定;细根生物量和凋落物生物量用烘干称量法测定;pH值用电位法(水土比为2.5:1)测定。
土壤酶活性参考关松荫[23]的方法测定,蔗糖酶用3,5-二硝基水杨酸比色法测定;土壤过氧化物酶和多酚氧化酶以邻苯三酚为底物,用没食子素比色法测定;土壤β-葡萄糖苷酶用对硝基酚比色法测定。土壤基本理化性质见表2。
林分类型
Stand types全碳
TC content/
(g·kg−1)全氮
TN content/
(g·kg−1)碳氮比
C:N全磷
TP content/
(g·kg−1)全钾
TK content/
(g·kg−1)pH值 土壤容重
Bulk density
/(g·cm−3)土壤含水率
Water
content/%平均重量直径
Mean weight
diameter/mm杨树纯林 PPP 13.61±0.70 AB 0.60±0.04 A 25.67±0.70 B 1.45±0.03 B 18.68±0.62 B 8.27±0.10 A 1.64±0.03 A 20.75±1.42 A 0.71±0.08 A 柳树纯林 WPP 13.34±0.66 AB 0.57±0.05 A 29.97±1.04 B 1.53±0.02 AB 20.98±0.64 A 8.39±0.04 A 1.66±0.03 A 20.09±1.40 A 0.66±0.03 AB 杨柳混交林 PWM 13.97±0.74 A 0.53±0.07 B 26.60±1.39 B 1.52±0.02 AB 19.91±0.26 AB 8.21±0.13 A 1.64±0.03 A 22.19±1.05 A 0.52±0.01 BC 撂荒地 CK 12.49±0.48 B 0.39±0.04 B 46.40±3.65 A 1.56±0.03 A 20.16±0.23 AB 8.09±0.02 A 1.63±0.03 A 23.94±1.25 A 0.40±0.02 C 注:土壤理化性质数据为0~20、20~40、40~60、60~80、80~100 cm 5个土层的平均值。
Note: The data of soil physical and chemical properties were the average values in five soil layers (0-20, 20-40, 40-60, 60-80 and 80-100 cm).Table 2. Physical and chemical properties of soil in the sample plots
-
土壤有机碳储量(SOCS)计算公式如下:
SOCS= SOC × ρ × h × (1-δ)/10
式中:SOC为有机碳含量,ρ为土壤容重,h为土层厚度,δ为土壤砾石含量。
根据Blair等[8]的方法,以撂荒地为对照进行碳库管理指数(CPMI)的计算,公式如下:
CPMI=CPI × LI × 100%
式中:碳库指数(CPI)=样品SOC/对照SOC。
碳库活度指数(LI)=样品碳库活度(L)/对照碳库活度(L)。
L=EOC/(SOC-EOC)。
用Excel 2019进行数据整理和CPMI的计算,用SPSS 21对数据进行统计分析。利用Shapiro-Wilk法检验方程残差分布的正态性,SOCS、MBC、CPMI、蔗糖酶、多酚氧化酶和β-葡萄糖苷酶等指标符合正态分布,SOC、EOC、DOC和过氧化物酶等采用对数转换,利用方差分析(ANOVA)和最小显著差异法(LSD)比较各指标间差异,双因素方差分析见表3。采用Pearson法和冗余分析(RDA)进行相关性分析;用Origin 2018软件进行统计图制作。
指标
Indicator林分类型
Stand types土层
Soil layers林分类型 × 土层
Stand types × soil layersF P F P F P 有机碳含量 SOC 27.015 <0.05 35.437 <0.05 1.478 0.173 有机碳储量 SOCS 19.930 <0.05 45.566 <0.05 0.712 0.731 易氧化有机碳 EOC 18.988 <0.05 30.684 <0.05 0.853 0.598 可溶性有机碳 DOC 1.166 0.335 12.634 <0.05 0.568 0.854 微生物生物量碳 MBC 20.276 <0.05 8.193 <0.05 0.661 0.777 蔗糖酶活性 Invertase activity 4.550 <0.05 3.553 <0.05 0.809 0.639 过氧化物酶活性 Peroxidase activity 1.364 0.268 4.701 <0.05 1.493 0.167 多酚氧化酶活性 Polyphenol oxidase activity 23.941 <0.05 6.826 <0.05 1.150 0.350 β-葡萄糖苷酶活性 β-Glucosidase activity 0.470 0.705 8.490 <0.05 0.853 0.598 Table 3. Effects of stand types and soil layer on SOC, labile fractions and enzyme activities
-
由图1可知:试验地3种林分类型和对照各土层SOC和SOCS范围分别为1.03~5.88 g·kg−1和3.53~17.55 t·hm−2,在不同林分类型和土层间差异显著(P<0.05)(表3)。3种林分类型0~100 cm土层的SOC和SOCS显著高于对照,分别增加了45.2%~82.2%和47.7%~83.1%,其中,杨树纯林SOC显著高于柳树纯林和杨柳混交林(表4)。在垂直分布上,SOC随土层加深而降低,0~20 cm土层的SOC显著高于20~100 cm各土层(P<0.05)。与柳树纯林、杨柳混交林和对照相比,杨树纯林0~20、20~40、40~60、60~80和80~100 cm土层SOCS分别增加了14.6%~19.2%、24.1%~168.1%、27.9%~108.1%、20.5%~149.2%和11.3%~144.4%;0~100 cm土层杨树纯林的SOCS分别是柳树纯林、杨柳混交林和对照的1.23、1.24、1.83倍。
林分类型
Stand types有机碳含量
SOC/(g·kg−1)有机碳储量
SOCS/(t·hm−2)易氧化有机碳
EOC/(g·kg−1)可溶性有机碳
DOC/(mg·kg−1)微生物生物量碳
MBC/(mg·kg−1)杨树纯林 PPP 3.59 ± 0.20 A 57.84 ± 1.80 A 1.45 ± 0.08 A 43.63 ± 1.27 A 193.99 ± 5.78 A 柳树纯林 WPP 2.87 ± 0.33 B 46.89 ± 5.49 A 1.26 ± 0.03 A 43.45 ± 2.05 A 181.73 ± 2.08 A 杨柳混交林 PWM 2.86 ± 0.13 B 46.65 ± 3.87 A 1.39 ± 0.10 A 42.46 ± 1.25 A 176.89 ± 6.48 A 撂荒地 CK 1.97 ± 0.09 C 31.59 ± 2.23 B 0.98 ± 0.04 B 39.77 ± 2.14 A 145.44 ± 5.73 B 注:土壤有机碳及活性组分数据为0~20、20~40、40~60、60~80、80~100 cm 5个土层的平均值。土壤有机碳储量为0~100 cm土层的总和。
Notes: The data of SOC and labile fractions were the average values in five soil layers (0-20, 20-40, 40-60, 60-80 and 80-100 cm). The data of SOCS was the sum of 0-100 cm soil layer.Table 4. The content of SOC and labile fraction in 0-100 cm soil layer under various stand types
-
由图2可知:3种林分类型和对照不同土层的EOC、DOC和MBC范围分别为0.73~2.38 g·kg−1、35.96~62.53 mg·kg−1和134.94~226.84 mg·kg−1,EOC和MBC在不同林分类型和不同土层间差异显著(P<0.05),DOC仅在不同土层间差异显著(P<0.05)(表3)。在0~100 cm土层,3种林分类型的EOC和MBC均显著高于对照(P<0.05),杨树纯林、柳树纯林和杨柳混交林的EOC和MBC分别比对照平均提高了48.0%、28.6%、41.8%和33.4%、25.0%、21.6%(表4)。杨树纯林、柳树纯林、杨柳混交林0~20 cm土层的DOC显著高于对照(P<0.05),分别提高了37.9%、24.8%和21.3%。3种土壤活性有机碳组分含量在垂直分布上的规律相似,随土层深度增加呈先降低后微升趋势。
-
CMPI是基于对照土壤有机碳变化的稳定性指标,对照土壤的CMPI为100%。由图3可知:0~100 cm土层3种林分类型土壤的CPMI显著高于100%(P<0.05);0~20 cm土层CPMI为杨树纯林>柳树纯林>杨柳混交林,其中,杨树纯林显著高于杨柳混交林(P<0.05);40~100 cm土层杨柳混交林的CPMI显著高于杨树纯林和柳树纯林(P<0.05)。
-
由表5可知:与杨树纯林、柳树纯林和对照相比,杨柳混交林0~100 cm土层过氧化物酶和多酚氧化酶的活性较高,其中,杨柳混交林过氧化物酶活性显著高于柳树纯林(P<0.05),提高了26.8%,多酚氧化酶的活性显著高于杨树纯林和对照(P<0.05),分别提高了32.8%和76.2%;杨柳混交林土壤蔗糖酶活性在80~100 cm土层显著高于杨树纯林(P<0.05);各土层β-葡萄糖苷酶活性在不同林分间差异不显著。除柳树纯林多酚氧化酶活性在20~40 cm土层最高,其余各林分的4种土壤酶活性均在0~20 cm土层最高,其中,杨柳混交林0~20 cm和40~60 cm土层的过氧化物酶活性显著高于80~100 cm土层(P<0.05),0~20 cm土层的β-葡萄糖苷酶活性显著高于40~100 cm各土层(P<0.05)。双因素方差分析表明:土层对蔗糖酶、过氧化物酶、多酚氧化酶和β-葡萄糖苷酶活性均有显著影响,林分类型仅显著影响蔗糖酶和多酚氧化酶的活性,林分类型和土层的交互作用对4种酶活性影响不显著(表3)。
林分类型
Stand types土层
Soil layer/
cm蔗糖酶活性
Invertase activity/
(mg·g−1·d−1)过氧化物酶活性
Peroxidase activity/
(mg·g−1·d−1)多酚氧化酶活性
Polyphenol oxidase activity/
(mg·g−1·d−1)β-葡萄糖苷酶活性
β-Glucosidase activity/
(mg·g−1·d−1)杨树纯林
PPP0~20 17.55±0.21 Aa 11.74±0.59 Aa 35.64±2.55 ABa 10.93±1.34 Aa 20~40 16.51±0.50 Aab 7.38±1.30 Aab 32.98±3.07 Aab 3.17±1.31 Ab 40~60 16.28±0.48 Aab 4.82±0.44 Bb 19.61±3.36 Bc 5.12±0.65 Ab 60~80 15.99±0.42 Ab 4.85±1.44 Ab 18.13±3.77 Bc 3.89±1.02 Ab 80~100 15.78±0.27 Bb 6.64±1.92 Aab 22.49±3.89 BCbc 4.20±0.65 Ab 0~100 16.42±0.37 A 7.09±0.33 AB 25.77±1.32 B 5.46±0.54 A 柳树纯林
WPP0~20 16.79±0.46 Aa 7.66±0.56 Ba 35.39±3.14 ABa 8.52±1.49 Aa 20~40 16.21±0.12 Aa 4.92±0.15 Ba 35.93±3.68 Aa 4.76±1.59 Aa 40~60 16.14±0.40 Aa 5.83±0.87 Aba 29.77±3.62 ABa 4.81±0.93 Aa 60~80 16.25±0.10 Aa 6.44±1.30 Aa 34.47±1.40 Aa 6.15±0.72 Aa 80~100 16.28±0.17 ABa 5.43±0.82 Aa 30.60±1.00 ABa 5.02±1.76 Aa 0~100 16.33±0.02 A 6.05±0.41 B 33.23±1.52 A 5.85±0.22 A 杨柳混交林
PWM0~20 17.07±0.34 Aa 9.45±1.56 ABa 39.35±3.68 Aa 9.54±1.34 Aa 20~40 16.71±0.12 Aa 6.17±0.68 ABab 33.98±0.69 Aa 6.77±1.65 Aab 40~60 17.00±0.42 Aa 9.45±1.80 Aa 34.04±3.20 Aa 4.71±1.32 Ab 60~80 17.04±0.52 Aa 8.36±1.98 Aab 29.77±4.02 Aa 4.92±0.97 Ab 80~100 16.64±0.04 Aa 4.92±0.73 Ab 33.95±2.68 Aa 4.09±1.15 Ab 0~100 16.90±0.06 A 7.67±0.69 A 34.22±0.92 A 6.01±0.68 A 撂荒地
CK0~20 16.51±0.37 Aa 8.23±0.90 ABa 26.12±2.79 Ba 6.87±1.13 Aa 20~40 15.99±0.25 Aa 6.20±0.39 ABa 20.18±4.02 Ba 4.92±1.42 Aa 40~60 16.28±0.14 Aa 7.93±1.04 ABa 19.80±3.95 Ba 5.64±1.44 Aa 60~80 16.09±0.12 Aa 6.71±1.14 Aa 14.51±1.53 Ba 3.99±1.76 Aa 80~100 16.02±0.28 ABa 6.20±0.42 Aa 16.48±3.94 Ca 4.20±1.47 Aa 0~100 16.18±0.23 A 7.05±0.25 AB 19.42±0.43 C 5.12±0.24 A 注:不同大写字母表示相同土层不同林分之间差异显著(P<0.05),不同小写字母表示相同林分不同土层之间差异显著(P<0.05)。
Notes: Capital case letters that differ within same soil layers indicate significant (P< 0.05) differences among stand types, lower case letters that differ within same stand type indicate significant (P< 0.05) differences among soil layers.Table 5. Soil enzyme activities in different stand types
-
相关性分析结果(表6)表明:SOC、SOCS、EOC和CPMI与4种土壤酶均呈显著或极显著正相关,DOC与土壤多酚氧化酶和β-葡萄糖苷酶呈显著或极显著正相关,MBC与土壤蔗糖酶、多酚氧化酶和β-葡萄糖苷酶呈显著或极显著正相关。
指标
Index蔗糖酶活性
Invertase activity过氧化物酶活性
Peroxidase activity多酚氧化酶活性
Polyphenol oxidase activityβ-葡萄糖苷酶活性
β-Glucosidase activity有机碳含量 SOC 0.376** 0.326* 0.479** 0.427** 有机碳储量 SOCS 0.376** 0.379** 0.387** 0.484** 易氧化有机碳 EOC 0.418** 0.382** 0.564** 0.461** 可溶性有机碳 DOC 0.242 0.194 0.327* 0.383** 微生物生物量碳 MBC 0.284* 0.043 0.434** 0.331** 碳库管理指数 CPMI 0.432** 0.338** 0.401** 0.389** 注:**表示极显著相关(P<0.01),*表示显著相关(P<0.05)。
Note: **P<0.01, * P<0.05。Table 6. Pearson correlations between characteristics of soil organic carbon pools and soil enzyme activity
-
图4表明:在有机碳及活性组分和CPMI与环境因子的RDA中,第一轴和第二轴特征值分别为0.64和0.04,可累积解释其信息量67.31%,方差累积贡献率为99.88%。碳氮比和细根生物量对土壤有机碳及活性组分和CPMI的解释率分别为50.8%(P=0.002)和10.4%(P=0.004),表明碳氮比和细根生物量是土壤有机碳及活性组分和CPMI的显著影响因素。RDA中变量夹角的余弦值表示二者之间的相关性,因此,碳氮比与土壤有机碳及活性组分和CPMI呈负相关;细根生物量与土壤有机碳及活性组分和CPMI呈正相关。
Characteristics of Soil Organic Carbon Pools and Their Influencing Factors under Different Stand Types in the Plain Sandy Area
- Received Date: 2022-11-17
- Accepted Date: 2022-12-24
- Available Online: 2023-08-20
Abstract: