-
森林作为重要的生态系统,其碳储量是陆地碳的46%,但是近几十年来,森林对全球变暖的正反馈作用,导致森林土壤中碳的大量流失[1]。土壤碳库的动态变化主要体现在土壤活性碳库中。土壤活性有机碳受植被和土壤微生物影响强烈,并与大气二氧化碳存在较频繁的相互作用,因此在全球气候变暖的背景下,研究环境变量和土壤有机碳储量对未来气候变化的影响具有重要意义[1-2]。溶解性有机碳(dissolved organic carbon, DOC)、微生物生物量碳(microbial biomass carbon, MBC)和易氧化有机碳(readily oxidization carbon, ROC)作为土壤活性碳库的3个重要指标,由于其周转率高且快,所以相较于总有机碳(total organic carbon, TOC),更能反映森林环境随海拔梯度的变化[3]。DOC不仅能够从森林底部向矿质土壤进行水文输送,而且是矿质土壤中有机碳的重要来源,是全球碳循环中的重要碳通量[4]。MBC是土壤有机质的重要组成部分,占土壤有机碳的2%~5%,在养分循环和整个有机质动态中起着重要作用[5]。ROC指能被333 mmol·L-1高锰酸钾氧化的碳,在土壤碳库中其含量通常大于MBC,能代表驱动反硝化和甲烷生成过程所需的碳源[6]。由于碳和氮循环之间的相互作用,森林生态系统氮素的储存与碳库之间紧密耦合,在全球气候变化中起重要作用[7]。溶解性有机氮(dissolved organic nitrogen, DON)是微生物和植物的一种重要氮源,而微生物和植物通过产生和吸收各种温室气体从而对当地气候变化起作用[8]。微生物生物量氮(microbial biomass nitrogen, MBN)不仅在土壤氮矿化等过程中起重要作用,而且对气候变化具有敏感的生物指示作用[9]。土壤酶在许多涉及土壤有机质降解、转化和矿化的生化过程中起重要作用[10]。已有研究表明:土壤酶活性与土壤养分动态密切相关,尤其是DON[11]。土壤活性有机碳不仅依赖于TOC,也与土壤脲酶(urease, UE)、过氧化氢酶(catalase, CAT)等酶活性密切相关[12]。土壤酶如β-葡萄糖苷酶(β-glucosidase, βG)与土壤DOC和DON储量密切相关,DOC和DON储量对土壤碳通量和循环至关重要[8]。黄山具有明显的植被和土壤垂直地带性[13],关于该区域的森林土壤碳氮库及相关酶活性的研究报道较少。本研究以黄山4种不同海拔典型植被带的土壤为研究对象,分析不同植被带土壤碳氮及其活性组分与相关土壤酶活性关系,以期为评价黄山森林生态系统的土壤碳储量和生产力以及保护区森林土壤管理等提供依据。
HTML
-
黄山位于安徽省南部(30°01′~30°18′N,118°01′~118°17′E),景区面积为160.6 km2,地势高差达1 400 m,地处亚热带季风气候带。黄山植物群落分布自下而上为次生林、常绿阔叶林、常绿落叶阔叶混交林、落叶阔叶林、山地矮林与灌丛及山地草甸;土壤自下而上分别为黄红壤、黄壤、暗黄棕壤、酸性棕壤,山顶局部分布有山地草甸土及山地草甸沼泽土[13]。
-
2018年7月,沿黄山云谷寺-光明顶步道,以海拔、植被类型及实际情况为指标,按约300 m的间隔、远离步道处选择4个海拔植被带样地,样地基本情况见表 1。每个海拔植被带设置3个10 m × 10 m的样方,每个样方内按S型采集10个表层10 cm深的土样,用冰盒运回实验室。一份土样过2 mm筛后保存于4 ℃冰箱中用于溶解性有机碳、氮以及微生物生物量的测定,另一份风干过筛后用于其他指标的测定。
代码 植被类型 海拔/m 土壤类型 主要树种 EBF 常绿阔叶林 870 黄壤 细叶香桂 Cinnamomum subavenium,青冈 Cyclobalanopsis glauca,小叶青冈 Cyclobalanopsis myrsinifolia,薄叶山矾 Symplocos anomala,甜槠 Castanopsis eyrei EDFⅠ 常绿落叶阔叶混交林Ⅰ 1 100 棕壤 青冈,小叶青冈,窄基红褐柃 Eurya rubiginosa,马银花 Rhododendron ovatum,浙江新木姜子 Neolitsea aurata,白木乌桕 Sapium japonicum,野桐 Mallotus japonicus,板栗 Castanea mollissima EDFⅡ 常绿落叶阔叶混交林Ⅱ 1 440 棕壤 黄山杜鹃 Rhododendron maculiferum,映山红 Rhododendron simsii,青冈,马银花,薄叶山矾,多脉榆 Ulmus castaneifolia,红豆杉 Taxus chinensis AM 高山草甸 1 760 草甸土 龙须草 Eulaliopsis binata,地榆 Sanguisorba officinalis,野古草 Arundinella anomala Table 1. Basic situations of the sampling sites
-
总有机碳(TOC)用重铬酸钾氧化-外加热法测定;全氮(total nitrogen, TN)用开氏法测定;pH值用电位法测定,水土质量比为2.5:1.0[14]。微生物量碳(MBC)和微生物量氮(MBN)用氯仿熏蒸浸提法测定[15],熏蒸土和未熏蒸土均用0.5 mol·L-1的硫酸钾浸提30 min(土液质量比为1:4),提取液中的有机碳质量分数用有机碳分析仪测定,全氮质量分数用开氏消化法测定。wMBC=EC/0.45;wMBN=EN/0.45。其中:EC和EN分别为熏蒸和未熏蒸土壤有机碳、全氮质量分数之差,转换系数为0.45。溶解性有机碳(DOC)和溶解性有机氮(DON)的测定[16]:用去离子水振荡浸提(水土质量比为5:1)30 min,高速(15 000 r·min-1)离心10 min,用0.45 μm滤膜抽滤,用有机碳分析仪测定滤液碳浓度得到DOC,用过硫酸钾氧化-紫外分光光度法测定滤液TN,靛酚蓝比色法测定NH4+-N,紫外分光光度法测NO3--N,则wDON=wTN-(wNH4+-N+wNO3--N)。易氧化有机碳(ROC)用高锰酸钾氧化法测定,过0.25 mm的风干土加333 mol·L-1高锰酸钾溶液25 mL,震荡1 h后离心,取上清液稀释250倍,在分光光度计565 nm波长下比色[16]。
-
脲酶(UE)用靛酚蓝比色法测定,酶活性以24 h土中产生NH3-N表示(mg·g-1·d-1);β-葡萄糖苷酶(βG)用硝基酚比色法测定,酶活性以土中产生对硝基苯酚表示(mg·g-1·d-1);过氧化氢酶(CAT)用高锰酸钾滴定法测定,酶活性以土中消耗的高锰酸钾表示(mL·g-1·h-1)[17]。
-
调查样地林下植被种类、株数或多度、盖度、高度等,并计算其多样性和重要值。Shannon-Wiener多样性指数(H):H=-∑PilnPi;Simpson优势度指数(D):D=∑Pi2。其中:Pi为样方内物种的重要值,Pi=(相对盖度+相对高度)/2。
-
用SPSS 20软件统计分析,单因素方差分析(one-way ANOVA)和Duncan法比较不同海拔植被带间的差异。Origin 2016软件制图。用R 3.3.4软件Performance Analytics包做Pearson相关矩阵图;Vegan包做酶活性和活性有机碳、氮组分冗余分析(redundancy analysis,RDA),对酶的数据进行Hellinger转化,999次蒙特卡罗置换检验(Monte-Carlo permutation test)考察环境因子的显著性。
1.1. 研究区概况
1.2. 研究方法
1.2.1. 样品采集
1.2.2. 土壤活性组分及基本性质
1.2.3. 土壤酶活性
1.2.4. 林下植被多样性
1.2.5. 数据分析
-
不同植被带土壤碳及其活性组分呈现随海拔变化的特征(表 2)。各海拔植被带的土壤TOC和ROC质量分数差异均不显著。高山草甸的土壤MBC平均质量分数最高,EDFⅡ的DOC平均质量分数最高,且各海拔植被带的MBC和DOC质量分数存在极显著差异(P<0.01)。说明不同树种覆盖的植被带,导致土壤活性碳及微生物生物量的差异。
植被带 TOC/(g·kg-1) MBC/(mg·kg-1) DOC/(mg·kg-1) ROC/(g·kg-1) TN/(g·kg-1) MBN/(mg·kg-1) DON/(mg·kg-1) EBF 58.43 ± 17.16 a 664.36 ± 232.03 b 45.66 ± 19.72 c 21.51 ± 11.30 a 3.01 ± 0.70 b 82.28 ± 71.41 b 39.35 ± 18.71 b EDFⅠ 47.99 ± 7.64 a 614.31 ± 81.71 b 50.67 ± 5.4 9 bc 13.22 ± 1.88 a 2.85 ± 0.43 b 70.66 ± 46.52 b 54.41 ± 11.40 ab EDFⅡ 49.89 ± 9.78 a 627.00 ± 91.62 b 102.29 ± 17.58 a 16.27 ± 4.69 a 2.85 ± 0.60 b 65.70 ± 28.53 b 65.55 ± 15.77 ab AM 56.67 ± 4.13 a 1 331.57 ± 103.41 a 74.36 ± 10.41 b 17.77 ± 3.05 a 4.487 ± 0.20 a 260.99 ± 80.22 a 74.01 ± 3.01 a F 0.67 18.34 9.64 0.88 7.06 7.36 3.65 P ns ** ** ns ** ** ns 说明:不同字母表示差异显著(P<0.05)。**表示差异极显著(P<0.01),ns表示差异不显著 Table 2. Soil carbon and nitrogen and their labile fractions in different vegetation zones along an altitudinal gradient on Mount Huangshan
由表 2可知:高山草甸土TN,MBN及DON质量分数均最高,分别为4.487,260.99,74.01 mg·kg-1。常绿落叶阔叶混交林下土壤TN和MBN质量分数较低,分别为2.85和65.70 mg·kg-1;常绿阔叶林下的土壤DON质量分数最低,为39.35 mg·kg-1。说明高海拔土壤氮素有效性较高。不同海拔植被带土壤TN和MBN质量分数存在极显著差异(P<0.01),而各海拔植被带间的土壤DON质量分数无显著差异。
-
MBC/TOC称为土壤微生物熵[18]。从表 3可见:高山草甸土的微生物熵(2.35%)明显高于其他植被带,而常绿阔叶林的土壤微生物熵(1.12%)最低;常绿落叶阔叶混交林Ⅱ的DOC平均分配比例(20.65%)显著高于其他植被带,常绿阔叶林的DOC平均分配比例最低(7.95%),不同海拔植被带的微生物熵以及DOC/TOC存在显著差异(P<0.01)。总体上常绿落叶阔叶林的ROC分配比例最大(35.11%),常绿落叶阔叶混交林Ⅰ则最小(27.66%)。不同植被带土壤活性有机氮分配比例呈现不同的海拔变化特征。高山草甸土的MBN/TN(5.79%)明显高于其他植被带,但是高山草甸的DON/TN(1.66%)较低。相较于其他植被类型,常绿落叶阔叶混交林MBN的分配比例较低,DON的分配比例较高。说明不同植被类型影响土壤微生物生物量和活性,从而在一定程度上影响土壤DON分配比例。
植被带 MBC/TOC/% DOC/TOC/% ROC/TOC/% MBN/TN/% DON/TN/% EBF 1.12 ± 0.08 b 0.08 ± 0.03 c 35.11 ± 10.66 a 2.47 ± 2.07 b 1.41 ± 0.86 a EDFⅠ 1.29 ± 0.11 b 0.11 ± 0.01 bc 27.66 ± 2.72 a 2.38 ± 1.26 b 1.96 ± 0.60 a EDFⅡ 1.26 ± 0.06 b 0.21 ± 0.02 a 32.26 ± 2.98 a 2.42 ± 1.22 b 2.32 ± 0.41 a AM 2.35 ± 0.10 a 0.13 ± 0.01 b 31.28 ± 3.88 a 5.79 ± 1.53 a 1.66 ± 0.14 a F 119.12 25.59 0.78 3.51 1.44 P ** ** ns ns ns 说明:不同字母表示差异显著(P<0.05)。**表示差异极显著(P<0.01),ns表示差异不显著 Table 3. Percentages of labile organic carbon and nitrogen in different vegetation zones along an altitudinal gradient on Mount Huangshan
-
如图 1所示:脲酶活性在常绿阔叶林下达最高值,为1.81 mg·g-1·d-1,在常绿落叶阔叶混交林Ⅱ下达最低值,为0.73 mg·g-1·d-1,且这2种植被类型的脲酶活性差异显著(P<0.05)。β-葡萄糖苷酶和过氧化氢酶均在高山草甸表现为最高活性,分别是2.05 mg·g-1·d-1和19.12 mL·g-1·h-1;β-葡萄糖苷酶在常绿落叶阔叶混交林Ⅱ下的活性最低,为0.94 m·g-1·d-1,过氧化氢酶在常绿阔叶林下的活性最低,为11.30 mL·g-1·h-1。不同海拔植被带的β-葡萄糖苷酶(F=23.85, P<0.01)和过氧化氢酶(F=23.24, P<0.01)均呈极显著差异,脲酶呈显著差异(F=4.47, P<0.05)。
-
如表 4所示:不同海拔植被带的土壤碳氮及其活性组分与土壤酶及其他土壤因子之间存在一定的相关性。就土壤碳组分而言,TOC与TN(R=0.65,P<0.05)、ROC(R=0.90,P<0.001)存在显著正相关关系;MBC与TN(R=0.96)、MBN(R=0.92)、β-葡萄糖苷酶活性(R=0.86)呈极显著正相关关系(P<0.001),与过氧化氢酶(R=0.66,P<0.05)显著正相关;DOC与pH值显著负相关(R=-0.65,P<0.05)。就土壤氮组分而言,TN与MBN(R=0.89,P<0.001)、β-葡萄糖苷酶活性(R=0.76,P<0.01)极显著正相关;MBN与β-葡萄糖苷酶(R=0.82,P<0.01)、过氧化氢酶(R=0.66,P<0.05)活性均显著正相关;DON与过氧化氢酶活性(R=0.58,P<0.05)显著正相关,与脲酶活性(R=-0.97,P<0.01)呈极显著负相关关系。此外,本研究发现所选样地的林下植被多样性与各因子之间均无相关关系。
变量 TOC MBC DOC ROC TN MBN DON pH UE βG CAT H D TOC 1 0.51 0.14 0.90** 0.65* 0.50 -0.25 -0.38 0.47 0.16 -0.12 -0.04 0.02 MBC 1 0.20 0.34 0.96** 0.92** 0.39 -0.04 -0.08 0.86** 0.66* 0.10 0.01 DOC 1 0.06 0.21 0.26 0.27 -0.65* -0.49 0.01 0.34 0.14 -0.23 ROC 1 0.45 2.31 -0.18 -0.44 0.30 0.05 -0.31 -0.05 -0.04 TN 1 0.89** 0.37 -0.08 0.02 0.76** 0.57 0.14 -0.03 MBN 1 0.19 0.00 0.04 0.82** 0.66* -0.09 0.21 DON 1 -0.04 -0.79** 0.48 0.58* 0.53 -0.51 pH 1 0.14 0.26 0.16 -0.13 0.24 UE 1 -0.23 -0.40 -0.23 0.30 βG 1 0.81** 0.09 0.08 CAT 1 0.35 -0.17 H 1 -0.96** D 1 说明:**表示在0.01水平(双侧)上显著相关,*表示在0.05水平(双侧)上显著相关 Table 4. Pearson correlation analysis between environmental variables
以土壤酶活性作为响应变量,活性有机碳、氮组分作为解释变量,并结合不同海拔植被带样点进行RDA分析。由图 2可知:土壤酶活性与活性有机碳、氮组分关系显著,可以解释86.96%的变异,其中第1轴解释了75.11%的变异,第2轴解释了11.85%的变异。从各酶活性与环境因子夹角可知,土壤脲酶与DON有很强的负相关性,过氧化氢酶与DOC以及β-葡萄糖苷酶与MBC,MBN均有很强的正相关性,这基本上与相关性检验结果一致。蒙特卡罗检验结果显示:仅DON对于土壤酶的分布的影响达到显著水平(r2=0.695 7,P=0.006),土壤活性有机碳氮组分对土壤酶分布的决定系数从大到小依次为DON,MBC,DOC,MBN和ROC。