ACTAECOLOGICASINICA
生态学报
Vol.40,No.2Jan.,2020
DOI:10.5846/stxb201808031659
张彦军,郁耀闯,牛俊杰,龚兰兰.秦岭太白山北坡土壤有机碳储量的海拔梯度格局.生态学报,2020,40(2):629⁃639.EcologicaSinica,2020,40(2):629⁃639.
ZhangYJ,YuYC,NiuJJ,GongLL.TheelevationalpatternsofsoilorganiccarbonstorageonthenorthernslopeofTaibaiMountainofQinling.Acta
秦岭太白山北坡土壤有机碳储量的海拔梯度格局
张彦军1,2,郁耀闯1,2,∗,牛俊杰2,3,龚兰兰1,2
1宝鸡文理学院地理与环境学院,宝鸡 721013
2宝鸡文理学院陕西省灾害监测与机理模拟重点实验室,宝鸡 7210133太原师范学院汾河流域科学发展研究中心,晋中 030619
摘要:土壤碳库是陆地生态系统的重要碳库,山区大约占全球25%的陆地表面积,因此研究山区土壤有机碳储量的变化特征及其影响因素对丰富陆地生态系统土壤碳循环理论具有重要的意义。在秦岭太白山北坡上,海拔高度每隔50m设置一个采样点,研究土壤有机碳储量的海拔梯度变化特征及其影响因素。海拔梯度对太白山北坡的SOCD影响显著,且不同土层厚度的SOCD均呈现出随着海拔梯度的增加而增加的趋势,增加幅度高达10%—88%。不同气候带上的SOCD差异显著,呈现出亚寒带(3.63kg/m2)大于寒温带(3.40kg/m2)大于温带(3.39kg/m2)大于暖温带(3.30kg/m2)的趋势。SOCD也因植被带发生显著差异,最小值出现在低山扰动带仅为2.15kg/m2,而最大值出现在高山草甸带高达3.63kg/m2,平均2.94kg/m2。SOCD随着cm减少到中山区的33cm和高山区的10cm,相应的总SOCD从低山区的35.45kg/m2减少到中山区的28.84kg/m2和高山区的因此,在秦岭太白山北坡上,除了海拔梯度、气候带、植被带、土层厚度等因素以外,土层深度对土壤有机碳储量的影响不容忽视,这对准确预测该区域的土壤碳储量具有重要的理论和实际意义。关键词:太白山;土壤有机碳密度;土层深度;土层厚度;海拔梯度
土层厚度的增加呈现出减少的趋势,减少幅度为14%—70%。在太白山北坡上,随着海拔梯度的增加,土层深度从低山区的7812.29kg/m2,但是对应单位土层深度上的SOCD却从低山区的0.51kg/m2增加到中山区的0.87kg/m2和高山区的1.23kg/m2。
TheelevationalpatternsofsoilorganiccarbonstorageonthenorthernslopeofTaibaiMountainofQinling
1GeographyandEnvironmentalEngineeringDepartment,BaojiUniversityofArtsandSciences,Baoji721013,China
ZHANGYanjun1,2,YUYaochuang1,2,∗,NIUJunjie2,3,GONGLanlan1,2
2ShaanxiKeyLaboratoryofDisastersMonitoringandMechanismSimulation,BaojiUniversityofArtsandSciences,Baoji721013,China3ResearchCenterforScientificDevelopmentinFenheRiverValley,TaiyuanNormalUniversity,Jinzhong030619,China
Abstract:Soilcarbonpoolisanimportantcarbonpoolofterrestrialecosystems.Themountainousareasaccountforabout25%oftheworld′ssurfacearea.Studyonvariationsofthesoilorganiccarbonstorageinmountainousareasandtheirinfluencingfactorsisofgreatsignificanceforenrichingthesoilcarboncycletheoryofterrestrialecosystems.OnthenorthernslopeoftheTaibaiMountainintheQinling,thesamplingsitesweresetatevery50metersabovesealeveltoexplorecharacteristicsofthesoilorganiccarbonstoragealonganaltitudegradient.Theresultsshowedthataltitudegradienthadasignificantimpactonthesoilorganiccarbondensity(SOCD).TheSOCDofdifferentsoilthicknessesincreasedby10%—88%withtheincreasingaltitudegradient.TheclimatezonealsopresentedasignificantimpactontheSOCD.Theorderwas
基金项目:国家自然科学基金项目(41801069);陕西省科技厅项目(2018JQ3023);宝鸡文理学院博士科研启用费项目(ZK2017041);陕西省灾害监测与机理模拟重点实验室项目(2010JS074)收稿日期:2018⁃08⁃03; 网络出版日期:2019⁃11⁃04
∗通讯作者Correspondingauthor.E⁃mail:zhangyanjun@bjwlxy.cn;yuyaochuang@126.com
http://www.ecologica.cn
630 生 态 学 报 40卷
subfrigidzone(3.63kg/m2)>coldtemperatezone(3.40kg/m2)>temperatezone(3.39kg/m2)>warmtemperatezone(3.30kg/m2).TheSOCDshowedsignificantlyvariationalongdifferentvegetationzones,withanaveragevalueof2.94kg/m2.Theminimumvalueof2.15kg/m2occurredinthelowmountaindisturbancezone,whilethemaximumvalueof3.6kg/area,whereasitwas33cmand10cminthemiddleandhighmountainousarea,respectively.TotalSOCDwas35.45kg/m2appearedinthealpinemeadowzone.TheSOCDshowedadecreasingtrendwiththeincreaseofsoilthickness,witha
reductionrangeof14%—70%.OnthenorthernslopeofTaibaiMountain,thesoildepthwas78cminthelowmountainousm2,28.84kg/m2,and12.29kg/m2inthelow,middle,andhighmountainousarea,respectively.Asopposed,theSOCDrespectively.Theresultsindicatedthat,inadditiontotheelevationgradient,climaticzone,vegetationzone,andsoilstorageinthisarea.
thickness,theinfluenceofsoildepthonsoilorganiccarbonstoragecannotbeignoredforaccuratepredictionofsoilcarbonKeyWords:TaibaiMountain;soilorganiccarbondensity;soildepth;soilthicknesses;altitudegradient
atthedepthoftheunitsoilinthelow,middle,andhighmountainousareawas0.51kg/m2,0.87kg/m2,and1.23kg/m2,
半的人类提供生活所需的水分[1⁃2],其较大的海拔梯度为研究环境变化影响生物地球化学循环过程提供了优越的自然条件[3]。森林生态系统是陆地生态系统的主体,是陆地上最大的碳储存库和碳吸收汇[4],森林生物量碳贮量达283GtC,约占全球植被碳储量的77%,森林土壤的碳贮量约占全球土壤碳储量的39%[5]。因此,研究山区森林土壤有机碳储量的变化特征及其影响因素对应对未来气候变化具有重要的理论和实际意义。
减少[8⁃9]的趋势,土壤有机碳储量的这种海拔梯度格局可能与植被类型[10]、气候带[11]、土壤理化性质[12]、土壤微环境[12]、土层厚度[10]、土壤微生物群落结构和活性[13]等因素密切相关。除了上述几种因素以外,土层深度也是影响土壤有机碳储量的一个重要因素[14]。例如,在全球尺度上,0—100cm土层内的SOC储量大约2344×1015PgC[15],因此土层深度对土壤有机碳储量的影响不容忽视,这在土层深度受限的山地生态系统中显得尤为重要[6]。在山地生态系统中,土层深度随着海拔梯度的增加呈现出或者减少[16⁃17]或者增加[6,18]的趋势,但是在该区域中目前很少有研究在较大海拔梯度下探讨土层深度对土壤有机碳储量的影响[6],这对全面理解山地生态系统的土壤有机碳储量的海拔梯度格局具有作用。因此,在山地生态系统中,在已有的研究基础上进一步深入研究土壤有机碳储量的海拔梯度变化特征及其驱动因素对准确估算该区域的土壤碳储量具有重要的理论和实际意义。
理、地质分界线[19]。太白山作为秦岭的主峰,是我国东部的最高峰(海拔高度3767m),由于山体高大,由山麓到山顶有着明显的气候垂直带谱和复杂完整的植被垂直带谱[20⁃21]。在秦岭太白山上,针对海拔梯度的研究主要集中在植被的垂直分布和群落结构调查[21⁃22]、叶片光合特性和分解特征研究[23⁃24]、土壤肥力状况和成土因素调查[19,25]、土壤有机碳储量和土壤呼吸速率特征等方面的研究[26⁃27],然而对于土壤有机碳储量的海拔梯度变化特征及其影响因素的研究相对较少。基于此,本研究在秦岭太白山的北坡上,从山麓到山顶海拔高度每隔50m设置一个采样点,研究土壤有机碳储量的海拔梯度变化特征及其影响因素,以期为准确估算该区域的土壤有机碳储量提供重要的数据支撑和理论依据。1 材料与方法1.1 研究区域概况
太白山位于陕西省西南部眉县、太白县和周至县三县的交界,最高峰拔仙台海拔高度3767.2m,垂直落差
http://www.ecologica.cn
秦岭导致我国南北地区在气候、河流、植被、土壤、农业等方面存在显著的差异,是我国南北最重要的地为1502×1015PgC,而0—200cm土层内却含有1993×1015PgC,但是在0—300cm土层内的SOC储量却高达
已有的研究结果表明,在山地生态系统中随着海拔梯度的增加土壤有机碳储量呈现出或者增加[6⁃7]或者
山区大约占全球25%的陆地表面积,其至少拥有全球1/3的陆地生态系统植被物种多样性,并为全球一
2期 张彦军 等:秦岭太白山北坡土壤有机碳储量的海拔梯度格局 631
3000多米,地跨北纬33°49′01″—34°08′11″,东经107°41′23″—107°51′40″。太白山位于秦岭山脉中段,是秦岭的主峰,也是我国地区中东部第一高峰,其东西长约61km,南北宽约39km,山体近东西向展布。太白山自然地理位置特殊,居暖温带之南缘,亚热带之北界,冬季受蒙古冷气团的控制,夏季受太平洋副高压带的影响。该地区年平均气温1.8—2.1℃,降雨量800—900mm,属性季风气候区。太白山地处中国西北部温带至暖温带的过渡区,四季分明,是典型的内陆季风气候区。太白山海拔梯度大,海拔高度每上升100m,气温下降约0.21—0.58℃,具有显著的气候垂直地带性差异,由高至低依次为亚寒带、寒温带、温带、暖温带,具体详见表1。
表1 不同气候带特征
Table1 Characteristicsofdifferentclimaticzones
海拔高度/mAltitudegradient800—13001300—26002600—33503350—3767
气候带Climaticzone暖温带温带 寒温带亚寒带
年均降水/mmAverageannualprecipitation650—800900—1000800—900750—800
年均气温/℃Averageannualtemperature8.7—12.71.7—8.7-2.1—1.8-4.4—2.1
特征
Characteristic
温和湿润,干湿季分明,蒸发量小于等于降雨量
温凉湿润,冬长夏短,晚霜频繁,植物生长期易遭雪压、雪折和霜冻,蒸发量小于降雨量寒冷潮湿,风大雾多,日夜温差大,蒸发量小于降雨量
寒冷半湿润,风大雾雪多,冬长而无夏,蒸发量小于等于降雨量
太白山是由东太白(拔仙台)、西太白山(鳌山)及其间主脊(跑马梁)和两侧南北延伸的多条支脉构成的组合体,属于褶皱断块高山。受地质构造的控制,太白山北坡极为陡峻,多深切峡谷或嶂谷,由下到上可分为低山区、中山区、高山区三种地貌类型。其中,营头稍南至刘家崖属于黄土覆盖的石质低山区,海拔高度在800—1500m之间,起伏的地形,具有黄土地貌和石质山地地貌的综合特征,相对高差不大,黄土覆盖;刘家崖至放羊寺属于石质中山区,海拔高度在1500—3000m之间,以峰岭地貌发育为特征;放羊寺至太白山顶峰是以第四纪冰川地貌为特征的高山区,海拔高度在3000m以上。
太白山植被类型丰富,植被垂直带谱完整,垂直地带性明显,根据植物群落的外貌、种类组成、群落结构及环境特征,由下向上依次可以分为以下6个植被带[19]:
1)低山扰动带
太白山山麓海拔高度约为600—1000m为低山带,人为扰动较大,农田主要以小麦(Triticumaestivum)
和玉米(Zeamays)为主,果园以猕猴桃(Actinidiachinensis)为主,有散生的栓皮栎(Quercusvariabilis)林和侧柏(Platycladusorientalis)林,陡坡有栎林、混交林和侧柏林等,其中栓皮栎林是山地系统较稳定的植物群落,其群落的均匀度、多样性指数和物种丰富度都比较低。在丘陵农田区常见的散生树木有杨树(Populusspp.)、刺槐(Robiniapseudoaca⁃cia)、柳树(Salixspp.)、椿树(AilanthusaltissimaSwingleh)等。
2)锐齿栎林带
本带海拔高度约为1000—2000m,栎林带在太白山北坡自上而下主要为辽东栎(Quercuswutaishanica)
林、锐齿槲栎(Quercusalienavar.acuteserrata)林和栓皮栎(Quercusvariabilisvar.variabilis)林。辽东栎是一种落叶乔木,喜光,多分布在土层较厚和稍湿润的向阳山坡上,纯林较多,锐齿槲栎和栓皮栎纯林较少。
3)混交林带
本带海拔高度约为2000—2500m,是太白山北坡树种最复杂、植被类型最多的地区,主要植物群落有华华山松林主要分布在海拔高度1400m以上的山地,有纯林,也有与辽东栎林组成的混交林。华山松为常绿乔木,喜光、喜温和凉爽湿润的气候,能适应多种土壤,但对土壤水分要求严格,多分布于阴坡、山脊和山
http://www.ecologica.cn
山松(Pinusarmandii)林、栎林和落叶阔叶混交林。
632 生 态 学 报 40卷
图1 太白山北坡气候和植被垂直带谱
Fig.1 VerticalspectrumofclimateandvegetationalonganaltitudegradientonthenorthernslopeofTaibaiMountainofQinling
顶上。
落叶阔叶混交林分布于本带山谷地区,组成种类复杂,层次较多,且层次之间区分不明显。混交林下有不同高度的乔木和灌木;同时在林缘,尤其是沟边,藤本植物非常丰富。
4)桦木林带
太白山北坡海拔高度约在2300—2600m之间,生长着红桦(Betulaalbo⁃sinensis)林及其变种牛皮桦
(Betulaalbo-sinensisvar.septentrionalis)林。这两种桦树各自组成纯林,上段与巴山冷杉(Abiesfargesii)交错分布,下段与华山松交错分布。
红桦林多分布在气候温凉湿润、土壤偏酸性、地形起伏剧烈、植物生长立地条件较差的地区。在土壤非常瘠薄的地方,红桦常常被适应性更强的牛皮桦取代,后者是原生裸地上的先锋树种,能形成原生林。影响牛皮桦林分布的主导因子是坡向和海拔高度。
5)高山针叶林带
高山针叶林带主要分布于海拔高度2500—3300m之间,上段为太白红杉(Larixchinensis)林,下段主要
是巴山冷杉(Abiesfargesii)林。林下灌木有松花竹、蒲乐杜鹃、金背枇杷。与巴山冷杉相伴生的乔木有桦树(Betulaspp.)及槭树(Acersaccharum)等。太白红杉主要分布在3000—3300m的崩塌陡坡,林下植物较少,以苔草、藓类和地衣为主。
6)高山灌丛草甸
太白山北坡海拔高度在3300m以上即无乔木分布,此界限以上为20—30cm高的矮型或匍匐性高山灌
木。头花杜鹃(Rhododendroncephalanthum)灌丛是太白山高山灌丛的主要类型,与其相间分布的是高山柳20%。近山顶多为危崖峭壁,岩石分化崩塌甚烈,在岩石表面缀生许多壳状地衣,在稍有碎土之地,则有枝状地衣(Cladoniaspp.)。1.2 测定项目及方法1.2.1 土壤样品的采集
m×20m的样地作为采样点,为了保证样地的代表性,减少局地地形条件对土层深度和土壤有机碳积累等因
http://www.ecologica.cn
于2016年5月至2016年9月,在太白山北坡上沿着海拔梯度的变化,海拔高度每升高50m设置一个20(Salixcupularis)灌丛和高山绣线菊(Spiraeaalpina)灌丛。灌丛群落中常有禾叶蒿草密生,其覆盖度可达
2期 张彦军 等:秦岭太白山北坡土壤有机碳储量的海拔梯度格局 633
素的影响[6,15],选取的样地尽可能保证坡位、坡向、坡度等地形条件基本接近一致。在每块样地内按S型或对10—20cm、20—40cm、40—60cm、60—80cm和80—100cm(具体采样深度以实际采样为主,不够100cm的采集到母岩为止)6个层次分层采集各层土壤样品,装袋、标记编号。在采样的过程中利用GPS定位采样点的海拔高度、地理位置等信息。同时在每个样地内,挖一个深度约为100cm的土壤剖面(具体采样深度以实际采样为主,不够100cm的采集到母岩为止),除去土壤表层的植物体,利用体积约为100cm3的环刀,按照0—10cm、10—20cm、20—40cm、40—60cm、60—80cm和80—100cm分层采集土壤,立即带回实验室进行称重,并计算土壤容重(环刀法)。由于是土石山区,土壤中肯定含有不少的石粒,本研究中土壤容重的计算包括石粒,如果在采样过程中遇到母岩,则该层次的土壤样品和土壤容重均不采集。1.2.2 样品分析
新鲜样品混合均匀后,在室内进行自然风干,风干样品过0.15mm筛后,测定SOC(H2SO4⁃K2Cr2O7外加
角线法随机选取3—5个能代表整个样地的样点进行土壤样品的采集。用土钻(内径3cm)按0—10cm、
热法)含量(g/kg),土壤容重的测定利用环刀法。1.2.3 土壤有机碳密度的估算
土壤有机碳密度(SOCD)利用下列方程进行计算:
SOCD=
SOCiBDiDi∑i
m
式中,SOCD是土壤有机碳密度(kg/m2),SOCi是第i层的土壤有机碳含量(g/kg),BDi是第i层的土壤容重1.3 统计分析
(gcm-3),Di是第i层的土层厚度(m),m是土层的数量。
(1)
利用Excel软件进行数据的预处理,利用Sigmaplot软件进行有关基础图件的制作,并用SAS软件的UNIVARIATEplotnormal程序包对SOCD数据进行正态分布检验。利用SAS软件对不同海拔梯度、气候带和植被带下的SOCD进行了方差分析(GLM),当F检验显著时再进行3种海拔梯度(低山区、中山区和高山区)、4种气候带(亚寒带、寒温带、温带、暖温带)以及6种植被带(低山扰动带、锐齿栎林带、混交林带、桦木林带、高山针叶林带和高山草甸带)均值间的Duncan检验。此外,利用SAS软件的PROCreg程序包对分析获得的SOCD(y)和土壤理化性质(x)进行回归分析(线性和非线性)。2 结果分析
2.1 SOCD变异的影响因子分析
正态分布(W=0.8404,N=652,P<0.0001)。进一步的数据分析表明,气候带、山区类型、海拔梯度、土层厚度、气候带和土层厚度的交互作用、气候带与山区类型和土层厚度的交互作用以及气候带与海拔梯度和土层厚度的交互作用均对太白山北坡的SOCD影响显著(P<0.0001,表2)。2.2 海拔梯度和土层厚度对SOCD的影响
海拔梯度对太白山北坡的SOCD影响显著,且不同土层厚度的SOCD均呈现出随着海拔梯度的增加而增加的趋势(图2)。在0—10cm土层,SOCD从低山区的3.23kg/m2增加到中山区的3.71kg/m2和高山区的3.55kg/m2,分别增加了15%和10%;在10—20cm土层,SOCD从低山区的2.03kg/m2增加到中山区的2.3588%(图3)。
kg/m2,增幅达到16%;在20—40cm土层,SOCD在低山区和中山区分别为3.08kg/m2和5.78kg/m2,增加了
除了海拔梯度以外,土层厚度也是影响SOCD的一个不可忽视的因素。在太白山北坡上,SOCD基本呈太白山北坡的SOCD的变化范围为0.24—14.99kg/m2,均值为2.23kg/m2,变异系数高达76%,数据接近
现出随着土层厚度的增加而减少的趋势(图2)。SOCD在0—10cm、10—20cm、20—40cm、40—60cm、60—80cm和80—100cm土层依次为3.41kg/m2、2.11kg/m2、3.64kg/m2、1.83kg/m2、1.26kg/m2和1.08kg/m2,
http://www.ecologica.cn
634 生 态 学 报 40卷
减少幅度分别为68%、49%、70%、41%和14%。
表2 太白山北坡SOCD变异的影响因子方差分析
Table2 VarianceanalysisofinfluencefactorsofSOCDvariationonthenorthernslopeofTaibaiMountain
项目ItemsQSHTQ×T
df3525713222
SS194.6190930822.6665129306.915869099.86859687.624281741.1473662
MS64.873031020.573683115.820509961.383173814.26694243.81214082.5832506
F359.51114.01340.1779.0621.1314.3287.67
Pr>F<0.0001<0.0001<0.0001<0.0001<0.0001<0.0001<0.0001
Q×S×T
Q代表气候带,S代表山区分类,H代表海拔梯度,T土层厚度,df代表自由度,SS代表平方和,MS代表均方
Q×H×T340.9890784
图2 海拔梯度对SOCD的影响
Fig.2 EffectofaltitudegradientonSOCDonthenorthernslopeofTaibaiMountain
2.3 气候带对SOCD的影响
在太白山北坡上,气候带对SOCD的影响显著,且基本呈现出亚寒带大于寒温带大于温带大于暖温带的趋势(图4)。在0—10cm土层,SOCD在不同气候带间的差异不显著(P>0.05),但也呈现出亚寒带(3.63kg/m2)>寒温带(3.42kg/m2)>温带(3.39kg/m2)>暖温带(3.30kg/m2)的趋势;在10—20cm土层,寒温带(2.24kg/m2)和温带(2.13kg/m2)的SOCD分别是暖温带(1.86kg/m2)SOCD的1.20和1.15倍;在20—40cm土层,寒温带(6.12kg/m2)和温带(3.19kg/m2)的SOCD依次较暖温带(2.65kg/m2)的SOCD提高了层,SOCD在暖温带和温带之间的差异不显著。
131%和20%;在40—60cm土层,SOCD在暖温带和温带依次为2.12kg/m2和1.73kg/m2;在60—100cm土
http://www.ecologica.cn
2期 张彦军 等:秦岭太白山北坡土壤有机碳储量的海拔梯度格局 635
图3 山区对SOCD的影响
Fig.3 EffectofmountaintypeonSOCDonthenorthernslopeofTaibaiMountain
图4 气候带对SOCD的影响
Fig.4 EffectofclimaticzoneonSOCDonthenorthernslopeofTaibaiMountain
2.4 植被带对SOCD的影响
带(3.38kg/m2)、高山针叶林带(3.30kg/m2)和高山草甸带(3.63kg/m2)之间的差异不明显(P>0.05),但是与混交林带(2.83kg/m2)和桦木林带(4.13kg/m2)的差异明显(P<0.05);在10—20cm土层,SOCD在桦木林带最高高达2.46kg/m2,在低山混交林带最低仅为1.86kg/m2,而在锐齿栎林带(2.00kg/m2)、混交林带(2.15kg/m2)和高山针叶林带(2.23kg/m2)之间的差异不明显(P>0.05);在20—40cm土层,较低山扰动带的SOCD而言(2.65kg/m2),锐齿栎林带(3.06kg/m2)、混交林带(3.28kg/m2)、桦木林带(4.89kg/m2)和高山针叶林带(5.97kg/m2)的SOCD分别提高了15%、24%、85%和125%;在40—60cm土层,低山扰动带(2.12kg/m2)、混交林带(2.35kg/m2)和桦木林带(2.33kg/m2)的SOCD分别是锐齿栎林带(1.42kg/m2)SOCD依次较锐齿栎林带(1.18kg/m2)的SOCD提高了10%和35%;在80—100cm土层,混交林带(1.52kg/m2)的SOCD分别是低山扰动带(0.99kg/m2)和锐齿栎林带(1.04kg/m2)SOCD的1.54和1.46倍。SOCD的1.49、1.65和1.64倍;在60—80cm土层,低山扰动带(1.30kg/m2)和混交林带(1.59kg/m2)的
植被带对太白山北坡的SOCD影响显著,在0—10cm土层,SOCD在低山扰动带(3.28kg/m2)、锐齿栎林
图5 植被带对SOCD的影响
Fig.5 EffectofvegetationzoneonSOCDonthenorthernslopeofTaibaiMountain
http://www.ecologica.cn
636 生 态 学 报 40卷
2.5 土层深度对SOCD的影响
在太白山北坡上,随着海拔梯度的增加,土层深度发生了显著的改变(P<0.05),平均从低山区的78cm,减少到中山区的33cm,然而在高山区仅有10cm,减少幅度分别达87%和70%(图6)。随着土层深度的改变,总SOCD随着海拔梯度的增加呈现出减少的趋势(P<0.05),且在低山区(35.45kg/m2)和中山区(28.84kg/m2)的总SOCD分别较高山区(12.29kg/m2)的总SOCD减少了65%和57%(图6)。但是,相应的单位土层深度上的SOCD却呈现出增加的趋势,从低山区的0.51kg/m2,增加到中山区的0.87kg/m2,进而增加到高山区的1.23kg/m2,增加幅度依次为141%和41%(图6)。
图6 海拔梯度对土层深度和SOCD的影响
Fig.6 EffectofaltitudegradientonsoildepthsandSOCDonthenorthernslopeofTaibaiMountain
土层深度、总SOCD和单位土层深度上的SOCD在不同植被带下也具有类似的变化规律,土层深度在低山扰动带下平均90cm、在锐齿栎林带下平均为80cm、在混交林带下平均为66cm、在桦木林带下平均为3858%。与此相对应的,总SOCD在不同植被带下呈现出减少的趋势,且总SOCD在低山扰动带下平均为40.75kg/m2、在锐齿栎林带下平均为34.41kg/m2、在混交林带下平均为33.25kg/m2、在桦木林带下平均为30.69cm、在高山针叶林带下平均为24cm、而在高山草甸下仅有10cm,降低幅度分别高达89%、88%、85%、74%和
kg/m2、在高山针叶林带下平均为20.93kg/m2,而在高山草甸下平均仅为14.50kg/m2。但是,单位土层深度上的SOCD在不同植被带下却呈现出增加的趋势,且SOCD在低山扰动带下、锐齿栎林带下、混交林带下、桦木林带下、高山针叶林带下和高山草甸下平均依次为0.48kg/m2、0.51kg/m2、0.55kg/m2、0.86kg/m2、0.91kg/m2和1.45kg/m2,增加幅度分别为6%、15%、79%、90%和202%。
http://www.ecologica.cn
2期 张彦军 等:秦岭太白山北坡土壤有机碳储量的海拔梯度格局 637
3 讨论
R2=0.51,P<0.0001),这与已有的研究结果基本相类似[6⁃7]。随着海拔梯度的改变(800—3800m),土壤微气太白山北坡上的SOCD在不同气候带下呈现出亚寒带>寒温带>温带>暖温带的趋势,这是因为不同气候带下的土壤微环境差异明显,例如年平均土壤温度在四个气候带下依次为0.22℃、2.78℃、8.21℃、13.72℃,而较低的土壤温度可能通过降低土壤微生物活动而有利于土壤有机质的积累[28⁃29]。此外,太白山北坡上的SOCD因植被带差异明显,而植被主要通过凋落物的输入(地上和地下)和凋落物的生化性质(木质素含量和碳氮比等因素)两个方面影响土壤有机碳的积累[30⁃32],且进一步研究发现,SOCD与地上部凋落物现存量呈现出抛物线关系(y=-0.00005x2+0.09x-1.69,R2=0.89,P<0.05)。气候带和植被带的改变必然会导致土壤的理化性质也会发生相应的改变(粘粒含量、土壤容重、土壤温度等因素),而这些因素对SOCD的影响不容忽P<0.001)、土壤pH(y=22.34x-105.53,R2=0.41,P<0.05)、土壤容重(y=20.17x+12.22,R2=0.89,P<0.01)、土壤温度(y=1.87x+17.36,R2=0.90,P<0.0001)均呈线性关系,这与已有的研究结果相类似[33⁃34]。视[33]。通过对数据的进一步挖掘发现,太白山北坡上的SOCD与土壤粘粒含量(y=48.28x+2.22,R2=0.70,
本研究结果表明,在太白山北坡上,随着海拔梯度的增加SOCD呈现出增加的趋势(y=0.001x+0.52,
候、植被带、土壤理化性质也会发生相应的改变,而这些因素均会对SOCD产生影响[3,10,12]。研究结果表明,
除了海拔梯度、气候带、植被带、土壤理化性质以外,土层厚度和土层深度也是影响SOCD的一个重要因
素。在本研究中,SOCD基本呈现出随着土层厚度的增加而减少的趋势(图2),且降低幅度高达14%—70%,这与已有的研究结果相类似[30,35]。SOCD在不同土层厚度下的分布特征可能与地下根系凋落物[30]、土壤微例如,在砖窑沟小流域内,有研究发现地下根系生物量随着土层厚度的增加呈现出减少的趋势,平均有超过64%的地下根系生物量分布在0—40cm土层,而仅有36%的地下根系生物量分布在40—100cm土层[30]。
在太白山北坡上,随着海拔梯度的增加土层深度呈现出减少的趋势,平均从低山区的78cm,减少到中山生物群落结构[36]、土壤水分含量[17]、土壤温度[17]以及土壤养分含量[17]在土壤剖面中的分布特征密切相关。
区的33cm,然而在高山区仅有10cm,减少幅度分别达87%和70%(图6),这与已有的研究结果相类似[16⁃17,37]。例如,在鼎湖山自然保护区内,土层深度的分布范围大概为18—105cm,集中分布于45—80cm,且土层深度最低分布在海拔高度最高的地点,仅有20—60cm,土层深度最高的分布在山脚(海拔高度最低),高达45—100cm,基本呈现出随着海拔梯度的增加而减少的趋势[37]。与此相对应,总SOCD随着海拔梯度的增加呈现出减少的趋势,平均总SOCD在中山区和高山区分别减少了65%和57%,但是单位土层深度上的SOCD却呈现出增加的趋势,平均单位土层深度上的SOCD在高山区和中山区分别增加了141%和41%(图6),因此在土层深度受限的土石山区土层深度对土壤有机碳储量的影响不容忽视[6]。例如,在全球尺度0—100cm土层内(1502×1015PgC)的1.2倍和1.6倍[15],这一结论也间接证明了土层深度在估算土壤有机碳储量的重要性[14]。
土层深度受局地地形条件(坡度、坡向、地形部位)影响显著[38],例如张文太课题组以的草地生态系上,土壤有机碳储量在0—300cm土层内(2344×1015PgC)分别是在0—200cm土层(1993×1015PgC)和
统为研究对象,利用GIS技术和相关分析的方法,探讨了海拔高度、坡度、坡向、平面曲率和剖面曲率等5个地形特征参数与土层深度之间的相关性,研究结果表明在塔里木盆地的北部,海拔高度、坡度和坡向与土层深度相关关系均为极显著水平(P<0.01)[39]。此外,土层深度的这种异质性可以通过土壤水分、土壤养分以及地下空间可利用性等方面对地上植物群落的物种组成和分布产生直接影响[40],从而影响土壤有机碳储量的空系凋落物量、土壤微生物群落结构和土壤养分含量等因素的数量差异有关:1)地下根系凋落物量[35],例如在现出对数减少的趋势(P<0.05),且在灌木植被、乔木植被、灌草植被以及草本植被下,表层(0—20cm)的细
http://www.ecologica.cn
间分布特征[30]。不同海拔梯度下,土层深度与土壤有机碳储量之间的关系可能与不同土层深度内的地下根黄土高原丘陵沟壑区的灌木、乔木、灌草以及草本植被下研究发现,地下根系凋落物量随着土层深度的增加呈
638 生 态 学 报 40卷
根凋落物量占剖面(0—100cm)根系凋落物量的比例依次为65%、63%、62%和81%[41];2)土壤微生物的群落结构[42],例如在黄土高原的人工刺槐林内,土壤微生物数量(放线菌、真菌和细菌)均随着土层深度的增加而呈现出增加的趋势,且表层(0—20cm)放线菌的数量分别是中层(20—50cm)和底层(50—90cm)放线菌数量的3倍和24倍,表层真菌的数量分别是中层和底层真菌数量的2倍和7倍,表层细菌的数量分别是中层和底层细菌数量的18倍和185倍[36];3)土壤养分含量,例如在莽山国家森林公园的常绿阔叶林生态系统中,底层(60—100cm)的水解氮、有效磷和有效钾占整个剖面的(0—100cm)的比例依次50%、45%和91%[43]。除了上述的地下根系凋落物量、土壤微生物群落结构和土壤养分含量以外,不同土层深度内的土壤但具体原因有待我们作进一步的试验进行论证。4 结论
太白山北坡上SOCD的变化范围为0.24—14.99kg/m2,均值为2.23kg/m2,变异系数高达76%,且SOCD温度[16]、土壤水分含量[17]以及土壤氧气浓度[36]的差异也可能导致不同土层深度内的土壤有机碳储量不同,
呈现出随着海拔梯度的增加而增加的趋势。太白山北坡上SOCD的海拔梯度格局与气候带、植被带、土壤理化性质、土层厚度等因素有关。此外,在太白山北坡上土层深度和总SOCD随着海拔梯度的增加呈现出减少的趋势,但是单位土层深度上的SOCD却呈现出增加的趋势。上述研究结果对准确估算土石山区的土壤碳储量具有重要的理论和实际意义。
致谢:感谢王阳妮、张瑜、封旭升对采样和试验给予的帮助,感谢张蓓蓓教授对写作的帮助。
参考文献(References):
[1] BarthlottW,LauerW,PlackeA.Globaldistributionofspeciesdiversityinvascularplants:towardsaworldmapofphytodiversity.Erdkunde,[2] ViviroliD,WeingartnerR,MesserliB.Assessingthehydrologicalsignificanceoftheworld′smountains.MountainResearchandDevelopment,[3] KörnerC.Theuseof′altitude′inecologicalresearch.TrendsinEcology&Evolution,2007,22(11):569⁃574.
2003,23(1):32⁃40.1996,50(4):317⁃327.
[4] 刘世荣,王晖,栾军伟.中国森林土壤碳储量与土壤碳过程研究进展.生态学报,2011,31(19):5437⁃5448.
Press,2000.
[5] IPCC.LandUse,LandUseChange,andForestry.SpecialReport,Inter⁃GovernmentalPanelonClimateChange.Cambridge:CambridgeUniversity[6] ZimmermannM,LeifeldJ,ConenF,BirdMI,MeirP.Cancompositionandphysicalprotectionofsoilorganicmatterexplainsoilrespiration[7] 朱凌宇,潘剑君,张威.祁连山不同海拔土壤有机碳库及分解特征研究.环境科学,2013,34(2):668⁃675.
decompositioninMediterraneanhighmountainsoils.Geoderma,2015,237⁃238:1⁃8.PlantNutritionandSoilScience,2010,173(1):30⁃38.
temperaturesensitivity?Biogeochemistry,2012,107(1/3):423⁃436.
[8] Gutiérrez⁃GirónA,Díaz⁃PinésE,RubioA,GavilánRG.Bothaltitudeandvegetationaffecttemperaturesensitivityofsoilorganicmatter[9] DjukicI,ZehetnerF,TatzberM,GerzabekMH.Soilorganic‐matterstocksandcharacteristicsalonganAlpineelevationgradient.Journalof[10] 程浩,张厚喜,黄智军,徐自坤,杨强,刘爱琴.武夷山不同海拔高度土壤有机碳含量变化特征.森林与环境学报,2018,38(2):[11] TianQX,HeHB,ChengWX,BaiZ,WangY,ZhangXD.Factorscontrollingsoilorganiccarbonstabilityalongatemperateforestaltitudinal[12] 李丹维,王紫泉,田海霞,和文祥,耿增超.太白山不同海拔土壤碳、氮、磷含量及生态化学计量特征.土壤学报,2017,54(1):160⁃170.[13] 厉桂香,马克明.土壤微生物多样性海拔格局研究进展.生态学报,2018,38(5):1521⁃1529.
Catena,2017,150:291⁃301.
[14] PremrovA,CumminsT,ByrneKA.Assessingfixeddepthcarbonstocksinsoilswithvaryinghorizondepthsandthicknesses,sampledbyhorizon.[15] JobbágyEG,JacksonRB.Theverticaldistributionofsoilorganiccarbonanditsrelationtoclimateandvegetation.EcologicalApplications,2000,
10(2):423⁃436.
gradient.ScientificReports,2016,6(1):18783.135⁃141.
http://www.ecologica.cn
2期 张彦军 等:秦岭太白山北坡土壤有机碳储量的海拔梯度格局 639
[16] 黄群山,叶源忠.浅析武夷山山地土壤理化性质的垂直分异.太原师范学院学报:自然科学版,2008,7(3):128⁃131.[17] 周运超,罗美.喀斯特小流域土壤厚度的影响因素.山地农业生物学报,2017,36(3):1⁃5.
[18] 耿增超,姜林,李珊珊,佘雕,侯磊.祁连山中段土壤有机碳和氮素的剖面分布.应用生态学报,2011,22(3):665⁃672.[20] 韩福琦,李兆元.太白山的气候.陕西气象,1980,(6):3⁃7.
[19] 雷梅,陈同斌,冯立,常庆瑞,阎湘.太白山北坡成土因素及不同土壤垂直带谱的比较.地理研究,2001,20(5):583⁃592.[21] 岳明.秦岭植物垂直带谱完整复杂.森林与人类,2015,(2):76⁃81.
[22] 任学敏.太白山主要植物群落数量分类及其物种组成和丰富度的环境解释[D].杨凌:西北农林科技大学,2012.
[23] 曹晓娟,刘建军,杨梅.太白山5种杜鹃属植物叶片光合特性及解剖结构的生态适应性研究.西北植物学报,2009,29(12):2483⁃2491.
学版,2005,33(1):64⁃68.
[24] 尚廉斌,王得祥,樊璐,刘建军,曹支敏,陈增宏.秦岭太白山北坡五种林分凋落叶分解状况初步研究.西北农林科技大学学报:自然科[25] 佘雕,朱盼盼,程苗苗,李浩宏,耿增超.太白山不同林分土壤肥力状况比较研究.西北林学院学报,2011,26(4):7⁃12.
[26] 郭建明,胡理乐,林伟,郑博福,李俊生.秦岭太白红杉林土壤有机碳密度研究.环境科学研究,2010,23(12):1464⁃1469.
[27] 屈冉,李俊生,罗遵兰,吴晓莆,赵彩云,汤博.秦岭太白山土壤CO2的释放特征及温度的影响.华北农学报,2010,25(2):188⁃193.
Biochemistry,2018,124:218⁃226.
[28] OldfieldEE,CrowtherTW,BradfordMA.Substrateidentityandamountoverwhelmtemperatureeffectsonsoilcarbonformation.SoilBiologyand[29] WangQK,LiuSE,TianP.Carbonqualityandsoilmicrobialpropertycontrolthelatitudinalpatternintemperaturesensitivityofsoilmicrobial[30] 张彦军,郭胜利.复杂地形条件下根系对土壤有机碳的贡献.环境科学,2019,40(2):961⁃969.
hydrolysisapproach.Geoderma,2002,107(1/2):109⁃141.
respirationacrossChineseforestecosystems.GlobalChangeBiology,2018,24(7):2841⁃2849.
[31] RoviraP,VallejoVR.Labileandrecalcitrantpoolsofcarbonandnitrogeninorganicmatterdecomposingatdifferentdepthsinsoil:anacid[32] HuYL,ZengDH,MaXQ,ChangSX.Rootratherthanleaflitterinputdrivessoilcarbonsequestrationafterafforestationonamarginal[33] 薛丽佳,高人,杨玉盛,尹云锋,马红亮,刘燕萍.武夷山土壤有机碳和黑碳的分配规律研究.林业科学研究,2011,24(3):399⁃403.[34] 姜泽宇.长白山森林土壤有机碳空间分布及其矿化特征研究[D].延吉:延边大学,2016.
[35] WangZQ,GuoSL,SunQQ,LiNN,JiangJS,WangR,ZhangYJ,LiuQF,WuDF,LiRJ,DuLL,ZhaoM.Soilorganiccarbon[36] 张社奇,王国栋,田鹏,郭满才.黄土高原刺槐林地土壤微生物的分布特征.水土保持学报,2004,18(6):128⁃131.
sequestrationpotentialofartificialandnaturalvegetationinthehillyregionsofLoessPlateau.EcologicalEngineering,2015,82:547⁃554.cropland.ForestEcologyandManagement,2016,362:38⁃45.
[37] 方运霆,莫江明,BrownS,周国逸,张倩媚,李德军.鼎湖山自然保护区土壤有机碳贮量和分配特征.生态学报,2004,24(1):135⁃142.[38] 李晓莎,许中旗,赵娱,李校.冀北山地干旱阳坡土壤厚度及土壤养分的空间异质性.河北农业大学学报,2018,41(1):24⁃30,48⁃48.[39] 柴强,盛建东,武红旗,刘耘华,张文太.草地土壤厚度及与地形特征的相关性.中国草地学报,2015,37(5):18⁃23.
(10):2023⁃2032.
[40] 李周,赵雅洁,宋海燕,张静,陶建平,刘锦春.喀斯特土层厚度异质性对草地群落结构和优势种生长的影响.草业科学,2017,34[41] 王彦丽.不同植被恢复措施下剖面根系与SOC的分布特征.水土保持研究,2013,20(6):19⁃23,43⁃43.
andBiochemistry,2018,121:154⁃164.
[42] HsiaoCJ,SassenrathGF,ZeglinLH,HettiarachchiGM,RiceCW.Verticalchangesofsoilmicrobialpropertiesinclaypansoils.SoilBiology[43] 谢寄托.莽山常绿阔叶林土壤养分分布规律研究[D].长沙:中南林业科技大学,2014.
http://www.ecologica.cn
因篇幅问题不能全部显示,请点此查看更多更全内容
Copyright © 2019- nryq.cn 版权所有 赣ICP备2024042798号-6
违法及侵权请联系:TEL:199 1889 7713 E-MAIL:2724546146@qq.com
本站由北京市万商天勤律师事务所王兴未律师提供法律服务