土壤环境因子对克氏针茅草原生态系统碳通量变化的影响

1土壤环境因子对克氏针茅草原生态系统碳通量变化的影响薛红喜1，李琪2,3，黄瑜2,3，王云龙2,3，吴东丽1（1.中国气象局气象探测中心，北京100081，2.南京信息工程大学江苏省大气环境监测与污染控制高技术研究重点实验室，江苏南京210044；3.南京信息工程大学环境科学与工程学院，江苏南京210044）摘要：本文利用涡度相关法对对2008年至2010年期间克氏针茅草原生态系统的生态系统初级生产力（GEP）、生态系统呼吸（Reco）和净生态系统碳交换（NEE）的变化进行观测，探讨了三者的日尺度变化规律及土壤环境因子（土壤温度、土壤水分）对GEP、Reco、NEE的影响。结果显示：土壤环境因子是克氏针茅草原生态系统碳通量的重要驱动因子；2008至2010年各生长季内，土壤环境因子对克氏针茅草原生态系统年尺度内的GEP呈现极显著相关性；土壤环境因子也与2008至2010年土壤环境因子对克氏针茅草原生态系统Reco、NEE呈现极显著相关性；土壤温度的升高或者土壤水分的增加都能提高克氏针茅草原生态系统的初级生产力及呼吸作用及净碳通量。其中，土壤温度对这三者的作用更为明显。关键词：土壤环境因子；碳通量；涡度相关法；草原生态系统1绪论碳循环是生态系统最重要的循环过程之一。作为陆地生态系统重要的组成部分，草地生态系统是世界上分布最广的植被类型之一，它覆盖了几乎20%的陆地面积，净初级生产力约占全球陆地生物区净初级生产力的1/3[1]。草原在地球碳循环过程中有非常重要的作用。比较其它陆地生态系统，草原的年际间总生产力变化幅度更大，以致于难以确定草原生态系统是碳源或者汇。所以，研究草原生态系统碳循环的过程和作用机理对研究碳循环至关重要。土壤是影响生态系统碳通量变化最直接的因素，土壤呼吸与生物光合作用、呼吸作用积累的生物量直接导致生态系统碳通量的变化。土壤是CO2重要的源和汇，土壤呼吸是生态系统碳通量的重要组成部分和流通过程。土壤是影响生态系统碳通量变化最直接的因素。因此研究土壤温度和水分对草原生态系统碳通量变化的影响，找出两者之间的联系对研究草原生态系统碳通量的特征、研究全球碳通量及碳通量对气候变化的影响有着重要的意义。本文研究内容为土壤环境因子对克氏针茅草原生态系统碳通量变化的影响。通过整理分析2008-2010年内蒙古克氏针茅草原日尺度上的净生态系统碳交换（NEE）的变化特征和土壤环境因子的特征，对克氏针茅草原生态系统碳通量变化进行研究，找出其特点。并分别对日尺度土壤环境因子变化对比分析，找出土壤环境因子和克氏针茅草原生态系统碳通量的联系，为准确认识克氏针茅草原生态系统碳通量的环境控制机制提供科学依据。2研究区域概况和研究方法2.1研究区域概况研究区位于内蒙古自治区锡林浩特市东北24km处(44°08′31N，116°18′45E)。此处海拔高度大概1000m，下垫面开阔平坦。该区处于内蒙古典型的克氏针茅草原区，群落里以克氏针茅和羊草为优势草种，伴生有糙隐子草、冰草等禾本科草类以及菊科的冷蒿等。这片草原位于北半球中基金项目：干旱气象科学研究基金项目（编号：IAM200919）资助。作者简介：薛红喜（1978.12-），男(汉族)，山西河津人，博士，高工，主要研究方向为大气探测、农业气象。E-mail:[e-mail]通讯作者：李琪，讲师，主要研究方向农业气象、生态系统碳循环。E-mail:[e-mail]2纬度地区，接近内陆，属温带半干旱草原气候，年平均气温为2.6℃，1月-22.0℃，7月均温18.3℃。日照充足，年均日照总时数约为2600小时。春秋冬三季干燥少雨，降水主要集中在夏天6至9月，年降水量约为290毫米。土壤为典型的栗钙土。2.2观测系统本研究使用涡度相关法观测克氏针茅草原生态系统碳通量。2003年8月装配的开路式涡度相关系统由三维超声风速仪、CR5000数据采集器及开路式LI-7500CO2/H2O分析仪构成，高度为2米。使用热电偶温度传感器Model105T测定10cm土壤处的土壤温度；10cm深处的土壤水分含量则利用土壤水分时域反射仪ModelCS615进行测量。通量采集频率为0.1Hz，土壤环境要素采样频率均为2Hz，由数据采集器CR23X按30min步长计算并存储平均值。2.3数据处理方法选取了克氏针茅草原生态系统2008年至2010年3年的碳通量观测数据，进行坐标旋转、WPL校正等处理，并对由系统故障、供电、下雨或检修等因素造成的数据丢失进行了剔除。采用查表法和日平均变化法(对缺失气象数据而无法用查表法以及非线性回归法的插补，采用日平均变化法)对通量数据进行插补按30min步长计算平均值，再分别计算日统计数据[24]。利用Excel和SPSS等软件探讨土壤环境因子和克氏针茅草原生态系统碳通量变化的相关关系，从而深入理解气土壤环境因子和克氏针茅草原生态系统碳通量变化的影响方式[16,24]。2.4基本理论公式生态系统净碳通量的NEE定义式:TTNEEcvacc)/)(1()/(（I）上划线表示某时间段内的平均值，撇号表示脉动，即瞬时值与平均值的偏差。式（I）中，右边第一项为CO2湍流通量，即单位时间垂直方向上通过单位截面积的CO2的量，为垂直风速，c为CO2在此处的浓度；右边第二、第三项为WPL校正项，其中a为干空气密度，v为水汽密度。此值经过三次坐标旋转后，即为可使用的NEE值。呼吸作用是夜间生态系统植被与大气间的NEE的唯一来源，生态系统夜间的Reco可用Lloyd-Taylor方程进行插补：)11(,000TTTTErefecoecoKrefeRR（II）式（II）中，refecoR,为参考温度下的生态系统夜间呼吸，refT为参考温度，取值为283.15K，0T为温度试验常数，取值为227.13K，0E为活化能，取值为308.56K。生态系统白天的呼吸根据夜间生态系统呼吸数据建立的函数关系外延到白天后获得的。由于白天植被与大气之间的净生态系统CO2交换（NEE）是生态系统初级生产力（GEP）和生态系统呼吸（Reco）的总和，因此GEP就可定义为[25-28]：ecoRNEEGEP（III）3结果分析3.1克氏针茅草原生态系统碳通量日总值变化图1显示了克氏针茅草原生态系统2008-2010年全年NEE日总值的变化情况。从图1看，2008年生长季（4月26日-11月20日）GEP日总值变化呈“W”型，图中出现两个波谷、两个波峰。在第3140、180、240天出现波谷，第150、200天附近出现波峰。在第200天附近达最大值。2009、2010年生长季GEP日总值变化也呈现“W”型，但形状不明显。三年GEP的变化范围基本一致。-5-4-3-2-10050100150200250300350400天数GEP/g•m-2•d-1-3.5-3-2.5-2-1.5-1-0.50050100150200250300350400天数GEP/g•m-2•d-1-5-4-3-2-10050100150200250300350400天数GEP/g•m-2•d-1图12008至2010年GEP日总值变化（左上图为2008年，右上图为2009年，下同）Fig.1TheGEPchangesfrom2008to2010(lefttopis2008,rightisfor2009,thesamebelow)图2中2008至2010年Reco日总值变化均呈倒“U”型。Reco先随时间变化增大，至200天附近出现最大值，然后随着时间变化减小到0值。对比三年Reco的最大、最小值发现，2008年Reco最小值，约为0.068/g·m-2·d-1，最大值最大，约为2.133/g·m-2·d-1；而2009年Reco最小值为0.105/g·m-2·d-1,最大值为1.692/g·m-2·d-1；2010年Reco最小值为0.200/g·m-2·d-1，最大值为1.518/g·m-2·d-1。从图3可看出NEE日总值呈W型，其波峰波谷所在处同GEP基本保持一致，变化趋势与GEP的日总值变化相同，在第200天附近达最大值。从NEE值的波动范围来看，2008年波动幅度最大，2010年次之，2009年NEE值波动幅度最小。这也与三年GEP的波动幅度相符。说明克氏针茅草原生态系统固碳能力的强弱与此生态系统初级生产力的变化一致。00.511.522.5050100150200250300350400天数Reco/g•m-2•d-100.20.40.60.811.21.41.61.8050100150200250300350400天数Reco/g•m-2•d-100.20.40.60.811.21.41.6050100150200250300350400天数Reco/g•m-2•d-14图22008年至2010年Reco日总值变化Fig.2TheRecochangesfrom2008to2010-4-3-2-10120100200300400天数NEE/g•m-2•d-1-2-1.5-1-0.500.511.50100200300400天数NEE/g•m-2•d-1-4-3-2-10120100200300400天数NEE/g•m-2•d-1图32008至2010年NEE日总值变化Fig.1TheNEEchangesfrom2008to20103.2土壤环境因子日变化3.2.1土壤温度日变化从图4可知，2008至2010年土壤温度的日变化趋势基本相，呈倒“U”型。从0到第200天，土壤温度先随时间变化升高，至第200天附近达最大值，而后随着时间变化减小，呈对称图形。-30-20-100102030400100200300400天数土壤温度/℃-20-15-10-50510152025300100200300400天数土壤温度/℃-20-15-10-5051015202530350100200300400天数土壤温度/℃图42008至2010年土壤温度的日变化Fig.4Thetemperaturechangesofsoilfrom2008to201053.2.2土壤水分日总值变化由图5可知，2008至2010年土壤水分随时间变化先增后减。最小值在0.049~0.056m3·m-3间波动，其中2008、2010年土壤水分于第200天达最大值，这与碳通量变化呼应。00.050.10.150.20.250.30100200300400天数土壤水分/m3•m-300.050.10.150.20.250100200300400天数土壤水分/m3•m-300.050.10.150.20.250.30100200300400天数土壤水分/m3•m-3图52008至2010年土壤水分的日变化Fig.5Thesoilmoisturechangesfrom2008to20103.3土壤环境因子对克氏针茅草原生态系统碳通量的影响图6为2008至2010年年生长季克氏针茅草原生态系统GEP随土壤环境温度的变化图。从图中可以看出，生长季内，GEP与土壤温度呈二次方关系，显著性检验结果表明，GEP与土壤温度呈极显著的相关（p0.01，下同）关系。从拟合的趋势线上可看出，在研究区域内的常见的温度范围内，克氏针茅草原生态系统GEP绝对值随土壤温度的升高而增大，说明在生长季内，土壤温度对该生态系统的初级生产力的影响为正效应。y=0.0007x2-0.0869x-0.2868R2=0.392-5-4-3-2-101-20-10010203040土壤温度/℃GEP/g•m-2•d-1y=0.0014x2-0.057x-0.3425R2=0.2052-3.5-3-2.5-2-1.5-1-0.500.5-100102030土壤温度/℃GEP/g•m-2•d-1y=-0.0005x2-0.0381x-0.267R2=0.3319-5-4-3-2-10-10010203040土壤温度/℃GEP