生态学报  2017, Vol. 37 Issue (6): 1994-2004

文章信息

韩春雪 , 刘廷玺 , 段利民 , 吕扬 , 闫雪 , 李凯旋
HAN Chunxue , LIU Tingxi , DUAN Limin , LÜ Yang , YAN Xue , LI Kaixuan .
科尔沁沙地两种植被类型土壤呼吸动态变化及其影响因子
Dynamic changes in soil respiration of two vegetation types and their influencing factors in Horqin Sandy Land
生态学报. 2017, 37(6): 1994-2004
Acta Ecologica Sinica. 2017, 37(6): 1994-2004
http://dx.doi.org/10.5846/stxb201510031998

文章历史

收稿日期: 2015-10-03
网络出版日期: 2016-08-02
科尔沁沙地两种植被类型土壤呼吸动态变化及其影响因子
韩春雪 , 刘廷玺 , 段利民 , 吕扬 , 闫雪 , 李凯旋     
内蒙古农业大学水利与土木建筑工程学院, 呼和浩特 010018
摘要: 以科尔沁沙丘-草甸相间地区为研究区,运用LI-6400土壤呼吸配套系统对沙丘-草甸过渡带人工杨树林和固定沙丘小叶锦鸡儿群落土壤呼吸及其相关因子进行观测,结果表明,(1)生长期和生长末期人工杨树林保留和去除枯枝落叶土壤呼吸的日变化表现为"多峰型"。(2)人工杨树林和小叶锦鸡儿群落保留和去除枯枝落叶土壤呼吸的季变化峰值均出现在8月份,人工杨树林的土壤呼吸明显大于小叶锦鸡儿群落的土壤呼吸。(3)人工杨树林土壤呼吸的季变化与0-10cm土壤温度相关显著(P < 0.01),保留枯枝落叶土壤呼吸小于去除枯枝落叶。(4)小叶锦鸡儿群落土壤呼吸的季变化与0-10cm土壤含水量相关显著(P < 0.01),与人工杨树林不同,保留枯枝落叶土壤呼吸大于去除枯枝落叶。(5)人工杨树林土壤呼吸与土壤温度显著呈指数关系,小叶锦鸡儿群落土壤呼吸与土壤含水量显著呈幂函数关系,并运用归一化法,通过建立土壤呼吸LnRs与土壤温度T、土壤含水量Lnθ的双因子回归模型,得出土壤温度和土壤含水量分别对保留和去除枯枝落叶土壤呼吸的贡献率。
关键词: 科尔沁沙地沙丘-草甸相间地区     枯枝落叶     土壤呼吸     土壤温度     土壤含水量    
Dynamic changes in soil respiration of two vegetation types and their influencing factors in Horqin Sandy Land
HAN Chunxue , LIU Tingxi , DUAN Limin , LÜ Yang , YAN Xue , LI Kaixuan     
Water Conservancy and Civil Engineering College, Inner Mongolia Agricultural University, Hohhot 010018, China
Abstract: The soil respiration and related correlation factors of an artificial poplar forest in the Horqin Sandy Land-meadow transitional zone and a Caragana microphylla community in a fixed dune were observed in this study using the LI-6400 soil-respiration supporting system. Results demonstrated the following five phenomena. (1) Daily changes in reserved and removed litter from the artificial poplar forest during the growth and late-growth periods have multiple peaks. (2) Seasonal peaks of reserved litter and removed litter appear in August for the artificial poplar forest and the C. microphylla community. Soil respiration of the artificial poplar forest is much higher than that of the C. microphylla community. (3) Seasonal changes of soil respiration in the artificial poplar forest are significantly correlated with the soil temperature (0-10 cm from the surface) (P < 0.01). Changes in the reserved litter are smaller than changes in the removed litter. (4) Seasonal changes of soil respiration of the C. microphylla community are significantly correlated with soil moisture content (0-10 cm from the surface) (P < 0.01). Unlike soil from the artificial poplar forest, changes in the reserved litter are larger than those in the removed litter. (5) The soil respiration of the artificial poplar forest is exponentially correlated with soil temperature, whereas soil respiration of the C. microphylla community has a significant power function relationship with the soil moisture content. By establishing a two-factor regression model of soil respiration (ln Rs) versus soil temperature (T) and soil water content (ln θ), contributions of soil temperature and soil moisture content to reserved litter and removed litter were determined using the normalization method.
Key words: Horqin sandy land     litter     soil respiration     soil temperature     soil water content    

土壤呼吸是大气CO2的重要排放源, 也是全球陆地生态系统碳循环的重要组成部分, 在区域和全球尺度上对大气CO2浓度调控起着十分关键的作用[1], 其有机碳储存量达1500Pg C, 占整个陆地生态系统碳储量的67%[2], 土壤呼吸作用的微小改变都将引起大气CO2浓度的明显改变, 从而影响全球气候[3]

干旱半干旱地区占据了地球表面的2/5以上[4], 土壤呼吸也是干旱半干旱地区土壤碳损失的主要过程之一[5], 并且因为相对较少的土壤有机碳含量, 土壤呼吸成为对气候变化最敏感的一个生态系统特征[6]。土壤呼吸是一个复杂的生物化学过程, 即受到非生物因子的影响也受到生物因子和人类活动的综合影响[7-8]。有研究表明土壤温度和土壤含水量能解释土壤呼吸大部分的变异[9], 当土壤含水量不是限制条件时, 土壤呼吸主要受土壤温度的控制, 而当受到土壤含水量胁迫时, 土壤呼吸的变化在温度的驱动下受土壤含水量控制[10]。枯枝落叶位于土壤表层, 也是影响土壤呼吸的另一关键因子, 其数量及其分解速率在很大程度上影响着土壤有机质的形成和对植物养分的供应, 或通过影响水热因子, 土壤环境以及微生物种类和数量, 直接或间接地影响土壤碳通量[11-14]。Raich和Nadehoffer[15]、田祥宇等[16]、周小刚等[17]分析得到去除枯枝落叶使土壤呼吸明显降低。但陈四清等[18]、骆土寿等[19]、王丽丽等[20]分析得到枯枝落叶对土壤呼吸起屏蔽作用, 进而使土壤呼吸降低。干旱半干旱地区关于土壤呼吸的研究较少, 有关枯枝落叶对土壤呼吸的影响研究更少, 至今仍没有相关研究量化科尔沁地区土壤呼吸, 以及相关影响因素对土壤呼吸的贡献。

本文以科尔沁沙丘-草甸相间地区沙丘-草甸过渡带人工杨树林和固定沙丘天然小叶锦鸡儿群落为研究对象, 分析保留和去除枯枝落叶土壤呼吸的变化特征, 模拟分析土壤温度、土壤含水量对保留和去除枯枝落叶土壤呼吸的影响关系。本研究估算了同一区域不同植被类型土壤碳排放量和枯枝落叶对土壤呼吸的贡献, 了解土壤温度、土壤含水量对土壤呼吸的影响以及变化机理, 为干旱半干旱地区不同植被类型土壤碳库管理提供理论依据。

1 材料和方法 1.1 研究区概况

研究区位于内蒙古自治区通辽市科尔沁左翼后旗阿古拉镇, 地理坐标122°33′00″— 122°41′00″E, 43°18′48″— 43°21′24″N, 面积55km2, 地处科尔沁沙地东南缘, 境内海拔最高232m, 最低186m。区内沙丘、农田、草甸、湖泊相间分布, 为典型的沙丘-草甸相间地区。该区年平均降水量389 mm, 且主要集中在6—9四个月;年平均蒸发量1412mm, 且主要集中在4—9月, 5月蒸发量最大;年平均相对湿度55.8%;年平均气温6.6℃, 7月月平均气温最高, 为23.8℃, 年极端最低气温为-33.9℃, 年极端最高气温36.2℃;年平均风速3—4m/s, 7—9月风速最小。区内植被群落简单, 物种多样性较差。主要植物有差巴嘎蒿 (Artemisia halodendron Turcz.et Bess.)、沙蓬 (Agriophyllum squarrosum (L.) Moq.)、羊草 (Leymus chinensis (Trin.) Tzvel.)、三芒草 (Aristida adscensionis Linn)、少花蒺藜草 (Cenchrus pauciflorus Benth)、冷蒿 (Artemisia frigida Willd. Sp. Pl)、芦苇 (Phragmites australis (Cav) Trin ex Steud.)、黄柳 (Salix gordejevii Chang et SkV.)、小叶锦鸡儿 (Caragana microphylla)、榆树 (Ulmus pumila L.)、杨树 (Populus L.) 等, 植物生长主要依赖天然降水, 草甸地部分湿生植被还依赖地下水。区内沙丘地带性土壤和非地带性土壤广泛发育, 交错分布, 砂土、砂壤土与壤砂土是主要的土壤类型。研究区地理位置见图 1

图 1 研究区地理位置、地貌类型与土壤呼吸采样点分布 Fig. 1 Geographic position, geomorphic type, and distribution of soil-respiration sampling site in the research area

本文选择沙丘-草甸过渡带低缓沙丘人工杨树与固定沙丘小叶锦鸡儿这两种植被类型区为研究对象, 对人工杨树林和小叶锦鸡儿群落选定的调查样方内植被抽样调查显示, 人工杨树林树龄介于22—28年, 平均为25年, 树高平均为11.9m, 胸径平均为11.2cm, 树的株距与行距变化在2—3m不等间距。人工杨树林内尚分布有其他低矮乔木与草本植物, 乔木主要有构树 (Broussonetia papyrifera (Linn.) L′ Hér. ex Vent.)、稠李 (Prunus padus L.) 等;草本主要有水蒿 (Artemisia selengensis Turcz. ex Bess.)、角蒿 (Incarvillea sinensis. Lamarck)、小飞蓬 (Conyza canadensis (L.) Cronq.) 等;盖度约20%左右, 枯枝落叶厚度在1.2—2.1cm, 平均为1.6cm, 现存量为672.6g/m2, 518.3g/m2为未分解, 腐殖质层平均厚度约5cm。地下水位埋深为0.996m, 地表高程为194.32m。小叶锦鸡儿群落为2—5年生, 灌丛高平均为2.2m, 根茎多变化在0.89—1.69cm之间, 丛间距多小于2.5m, 株丛间尚分布有大籽蒿 (Artemisia sieversiana Ehrhart ex Willd.)、狗尾巴草 (Setaria viridis (L.) Beauv.)、叉分蓼 (Polygonum divaricatum L.) 等植物, 盖度约13%左右, 灌丛下枯枝落叶厚度一般在0.5—0.9cm, 未分解, 现存量为104.5g/m2, 没有腐殖质层, 地下水位埋深为3.79m, 地表高程为194.59m。样方内土壤理化性质及其它相关信息见表 1

表 1 人工杨树林和小叶锦鸡儿群落的基本信息 Table1 Basic information of artificial poplar forest and Caragana microphylla community
植被类型
Vegetation
土壤类型
Soil type
土壤有机质
Soil organic
matter/
(g/kg)
全氮TN/
(g/kg)
全磷TP/
(g/kg)
全钾TK/
(g/kg)
pH 植被根系量
Vegetation root
amount/
(g/0.5×0.5m2×0.2m)
人工杨树林
Artificial poplar forest
壤砂 3.06 1.49 0.091 22.83 8.0 72.13
小叶锦鸡儿群落
Caragana microphylla
community
砂土 2.09 0.63 0.063 28.79 7.8 18.21
  土壤取样深度为0—10cm, 所给数据为平均值
1.2 实验设置

图 2, 在人工杨树林内围绕安装数据采集器 (Campbell) 的自动气象站选取20m×20m的样方, 其内布设4个土壤呼吸测点, 4个测点形成一个近似边长为10m的正方形, 每个点上并排放置两个自制的PVC塑料圈, 共8个 (内径10cm, 高7cm, 厚5mm)。每组PVC塑料圈一个里面去掉枯枝落叶, 另一个保留枯枝落叶。同样在小叶锦鸡儿群落围绕安装数据采集器 (Campbell) 的自动气象站选取20m×20m的样方, 在样方内选择边长约4m成近似等边三角形分布的3棵小叶锦鸡儿灌丛, 其内布设3个土壤呼吸测点, 每个测点在距离相近小叶锦鸡儿灌丛根部1m的位置相邻放置两个PVC塑料圈, 一个去掉枯枝落叶, 另一个保留枯枝落叶, 塑料圈均放置在小叶锦鸡儿灌丛的冠幅内, 共6个 (内径10cm, 高7cm, 厚5mm)。

图 2 两种植被 (人工杨树林 (A)、小叶锦鸡儿群落 (B)) 土壤呼吸测试样方布设示意图 Fig. 2 Layout of test sample plots of soil respiration in the artificial poplar forest (A) and Caragana microphylla community (B)
1.3 测定方法

2014年5—10月, 在人工杨树林和小叶锦鸡儿群落布设的土壤呼吸测点, 用LI-6400便携式气体分析系统 (Li-Cor Inc., NE, USA) 和LI-6400-09土壤呼吸室对土壤呼吸进行测定, 每月测1—2次。每次测量时间为8:00—18:00, 人工杨树林和小叶锦鸡儿群落的测量步长分别为2h和1h。整个测量期间均保持PVC塑料圈位置不变, 并在每次测量前24h对塑料圈内的绿色植物齐地剪除, 对去掉枯枝落叶的塑料圈内的其他杂物一并去除。同时用便携式土壤水分测定仪 (Diviner2000, Sentek, Australia) 测定每个测点10cm深度内的平均土壤含水量。利用LI-6400自带的土壤热电偶探针测量10cm深处的土壤温度。此外, 针对人工杨树林, 在生长期8月和生长末期10月还进行了24h的连续观测, 每2h测定1次。年度土壤呼吸实验结束后对各测点枯枝落叶的厚度进行测定, 试验全部结束后, 在放置土壤呼吸PVC圈的周围采取50×50cm2×20cm的土样, 计算根系量。

1.4 数据处理与分析

以同一植被类型为基本单位, 将3个 (小叶锦鸡儿群落) 或4个 (人工杨树林) 重复样点观测的所有要素进行平均, 获得1h或2h的平均值, 用于分析两种植被类型区土壤呼吸及其环境要素的日内变化特征。将小时数据进行平均, 获得日平均值, 将每个月的日平均值再次平均获得月平均值, 用于分析两种植被类型区土壤呼吸及其环境要素的月变化。

用土壤呼吸与土壤温度建立指数模型[21]

(1)

式中, Rs为平均土壤呼吸 (μmol m-2s-1), T为平均土壤温度 (℃), ab为正常数。Q10为土壤呼吸的温度敏感性指标, 是指某一温度下土壤的呼吸速率与低于该温度10℃下土壤呼吸的比值。一般由下式确定[22]

(2)

用土壤呼吸与土壤含水量建立幂函数模型:

(3)

式中, Rs为平均土壤呼吸 (μmol m-2s-1), θ为平均土壤体积含水量 (%), cd为正常数, 1>d>0。部分学者[23-24]采用二次曲线模型模拟土壤呼吸与土壤含水量的关系, 考虑到研究区小叶锦鸡儿群落受土壤含水量的胁迫影响, 不能达到土壤含水量对土壤呼吸作用的极限情况, 且通过两种模型模拟精度的比较表明, 保留和去除枯枝落叶土壤呼吸与土壤含水量的二次曲线模型R2分别为0.838、0.769, 幂函数模型R2分别为0.778、0.708, 两种模型的决定系数相差不大, 为便于统一与双因子回归模型进行对比分析, 本文选择精度略低的幂函数模型来模拟土壤呼吸与土壤含水量的关系。

由式 (1)、式 (3) 可以看出, ln RsT、Lnθ成线性关系, 为了分离土壤温度T、土壤含水量θ对土壤呼吸Rs的贡献率, 选用以下双因子回归模型进行分析:

(4)

为了得到双因子回归模型中土壤温度和土壤含水量对土壤呼吸的贡献, 对ln RsT、lnθ变量进行归一化计算, 采用下式:

(5)

式中, ln Rs*T*、ln θ*分别为土壤呼吸的对数ln Rs、土壤温度T、土壤含水量的对数lnθ的归一化变量, γ为拟合常数, αβ分别为土壤温度T、土壤含水量θ的贡献系数。Y*为归一化变量值, YYminYmax为原变量的试验值、最小试验值、最大试验值。

2 结果和分析 2.1 人工杨树林土壤呼吸以及影响因素的日变化动态

在植被不同生长期内土壤呼吸及土壤温度的日变化曲线存在一定差异, 土壤呼吸的日变化曲线表现为“多峰型”(图 3)。人工杨树林在生长期8月20号, 保留和去除枯枝落叶土壤呼吸分别在16:00、12:00出现最大值, 在24:00、8:00出现最小值;保留枯枝落叶土壤呼吸在6:00、22:00出现明显次峰值;去除枯枝落叶土壤呼吸在2:00、6:00、20:00出现明显次峰值。人工杨树林在生长末期10月2号, 保留和去除枯枝落叶土壤呼吸分别在18:00、20:00出现最大值, 均在4:00出现最小值;保留枯枝落叶土壤呼吸在12:00出现明显次峰值;去除枯枝落叶土壤呼吸在2:00、10:00出现明显次峰值。去除枯枝落叶土壤呼吸在1d内始终大于保留枯枝落叶, 并且保留和去除枯枝落叶土壤呼吸日变化的平均值, 生长期8月20号的8.097μmol m-2s-1、11.872μmol m-2s-1大于生长末期10月份的2.952μmol m-2s-1、4.429μmol m-2s-1

图 3 人工杨树林去除和保留枯枝落叶土壤呼吸与土壤温度的日动态变化 Fig. 3 Daily dynamic changes of Removed litter and Reserved litter of artificial poplar forest against soil temperature

对日变化的土壤呼吸与土壤温度进行了相关分析, 结果如表 2。生长期8月20日保留和去除枯枝落叶土壤呼吸与前移2h的土壤温度相关性最强 (P < 0.05), 生长末期10月2日保留和去除枯枝落叶土壤呼吸分别与前移2h (P < 0.05) 和无滞后 (P>0.05) 的土壤温度相关性最强。按照相关系数最大的滞后或者前移时间, 绘制出人工杨树林保留和去除枯枝落叶土壤呼吸随土壤温度的变化过程, 如图 4, 可以看出, 随着土壤温度的升高, 保留和去除枯枝落叶的土壤呼吸都是增大的, 土壤温度最多能解释土壤呼吸的52.9%。

表 2 人工杨树林日土壤呼吸与不同前移、滞后时间土壤温度的相关系数 Table2 Correlation coefficient of daily soil respiration of artificial poplar forest with soil temperature under different forward and delay times
日期 (月-日)
Date (Month-day)
处理方式
Process mode
相关系数Correlation coefficient
无滞后
On time
滞后2h
2h delay
滞后4h
4h delay
前移2h
2h forward
前移4h
4h forward
08-20 保留枯枝落叶 0.434 0.075 -0.472 0.649* 0.573
去除枯枝落叶 0.447 0.110 -0.225 0.597* 0.321
10-02 保留枯枝落叶 0.601* 0.371 0.183 0.694* 0.455
去掉枯枝落叶 0.561 0.395 0.069 0.544 0.255
  *表示相关性在0.05水平上显著

图 4 人工杨树林去除 (y1)、保留 (y2) 枯枝落叶日土壤呼吸随土壤温度的变化 Fig. 4 Daily changes of Removed litter (y1) and Reserved litter (y2) of artificial poplar forest against soil temperature
2.2 土壤呼吸以及影响因素的月变化动态

人工杨树林和小叶锦鸡儿群落生长季的土壤呼吸波动较大 (图 5)。就人工杨树林, 8月20号保留和去除枯枝落叶土壤呼吸同时达到最大值, 分别为8.10、11.87μmol m-2s-1, 在整个生长季去除枯枝落叶土壤呼吸始终大于保留枯枝落叶, 二者的相关性非常显著 (P < 0.01)。就小叶锦鸡儿群落, 8月13号保留和去除枯枝落叶土壤呼吸达到最大值, 分别为8.43、6.97μmol m-2s-1。与人工杨树林不同, 整个生长季保留枯枝落叶土壤呼吸始终大于去除枯枝落叶, 并且二者的相关性极其显著 (P < 0.01)。在整个观测期, 人工杨树林保留和去除枯枝落叶土壤呼吸平均值 (5.82、7.75μmol m-2s-1) 大于小叶锦鸡儿群落保留和去除枯枝落叶土壤呼吸平均值 (5.80、4.67μmol m-2s-1)。

人工杨树林保留和去除枯枝落叶土壤呼吸与土壤温度的相关程度更高, 小叶锦鸡儿群落保留和去除枯枝落叶土壤呼吸与土壤含水量的相关程度更高 (表 3)。对生长季人工杨树林保留和去除枯枝落叶土壤呼吸与土壤温度进行指数模型拟合, 对生长季小叶锦鸡儿群落保留和去除枯枝落叶土壤呼吸与土壤含水量进行幂函数拟合, 其土壤温度对人工杨树林保留和去除枯枝落叶土壤呼吸解释的变异系数, 土壤含水量对小叶锦鸡儿群落保留和去除枯枝落叶土壤呼吸解释的变异系数, 一并列入表 4, 人工杨树林保留枯枝落叶土壤呼吸的Q10值2.97大于去除枯枝落叶土壤呼吸的Q10值2.60(表 4)。又进一步分析土壤温度和土壤含水量对土壤呼吸的共同作用, 对各变量进行归一化处理后, 得到的土壤呼吸与土壤温度和土壤含水量的双因子回归模型结果如表 5, 其模拟效果比单因子回归模型更好, 双因子回归模型中变量前的系数即为该变量对土壤呼吸的贡献率。

图 5 去掉和保留枯枝落叶土壤呼吸与土壤温度、土壤含水量的月变化 Fig. 5 Monthly changes of Removed litter and Reserved litter against soil temperature and soil moisture content

表 3 两种植被类型去掉和保留枯枝落叶土壤呼吸与土壤温度、土壤含水量的相关性 Table3 Correlation of the Removed litter and Reserved litter of two vegetation types with the soil temperature and soil moisture content
植被类型
Vegetation
处理方式
Process mode
10cm土壤温度
10cm soil temperature
10cm土壤含水量
10cm soil moisture content
人工杨树林 去除枯枝落叶 0.903** 0.532
Artificial poplar forest 保留枯枝落叶 0.944** 0.607*
小叶锦鸡儿群落 去除枯枝落叶 0.215 0.796**
Caragana microphylla community 保留枯枝落叶 0.180 0.853**
  *表示相关性在0.05水平上显著, **表示相关性在0.01水平上显著

表 4 两种植被类型去掉和保留枯枝落叶土壤呼吸与土壤温度、土壤含水量的单因子回归模型 Table4 Single-factor regression model of the correlation of ln Rs with T and ln θ
植被类型
Vegetation
处理方式
Process
mode
模型类型 a b c d R2 P < Q10
人工杨树林
Artificial poplar forest
去除枯枝落叶 指数 1.363
(0.11)
0.095
(0.17)
0.858 0.01 2.60
(0.16)
保留枯枝落叶 指数 0.799
(0.18)
0.109
(0.21)
0.908 0.01 2.97
(0.18)
小叶锦鸡儿群落
Caragana microphylla
去除枯枝落叶 幂函数 2.726
(0.23)
0.536
(0.15)
0.708 0.01
community
保留枯枝落叶 幂函数 3.264
(0.19)
0.575
(0.20)
0.778 0.01
  数值=平均值 (标准偏差)

表 5 两种植被类型归一化处理后土壤呼吸ln Rs与土壤温度T和土壤含水量lnθ的双因子回归模型 Table5 Two-factor regression model of the correlation of ln Rs with T and ln θ after normalization
植被类型
Vegetation
处理方式
Process mode
α β γ R2 P
人工杨树林 去除枯枝落叶 0.881(0.03) 0.155(0.07) -0.0326(0.03) 0.872 < 0.01
Artificial poplar forest 保留枯枝落叶 0.846(0.06) 0.223(0.11) 0.070(0.07) 0.938 < 0.01
小叶锦鸡儿群落 去除枯枝落叶 0.181(0.05) 0.716(0.1) 0.171(0.07) 0.748 < 0.01
Caragana microphylla community 保留枯枝落叶 0.116(0.05) 0.749(0.06) 0.193(0.09) 0.812 < 0.01
  数值=平均值 (标准偏差)
3 讨论 3.1 土壤呼吸的日变化特征与土壤温度的相关性

针对人工杨树林, 对比两次土壤呼吸的日变化, 生长期8月20日保留和去除枯枝落叶土壤呼吸大于生长末期10月2日保留和去除枯枝落叶土壤呼吸, 这主要是由于8月的水热条件较好, 根系呼吸强, 微生物活动活跃[25]。两次土壤呼吸日变化最值出现的时间也有所差异, 这是由于不同季节土壤呼吸日变化可能受到季节变化显著的土壤有机碳含量和植物物候期等因素的影响[26], 从而使两次土壤呼吸日变化表现出不同的模式。两次土壤呼吸均表现为“多峰型”曲线, 这一结果与骆土寿[19, 27]日变化比较复杂的研究结果相似, 可能是由于人工杨树林为壤砂土, 有机质含量较低, 地下水位埋深较浅, 较易受到外界影响因子的影响所致。亦或是一天中土壤温度与土壤含水量在不同时期交替作为主导因子来影响呼吸速率, 同时土壤温度和土壤含水率的垂直变化也会对根系和土壤微生物产生影响, 这也可能是导致土壤呼吸与土壤温度出现不同前移或滞后现象的原因。8月20日保留和去除枯枝落叶土壤呼吸均与前移2小时10cm土壤温度的相关性最好, 而10月2日保留和去除枯枝落叶土壤呼吸与10cm土壤温度的相关性表现不一致, 并且土壤温度对土壤呼吸日变化的解释程度不高 (如图 4), 这也进一步说明土壤温度以外的其他因素对土壤呼吸的影响也是很关键的。Pierrot[28]解释土壤含水量虽然在一天中的变化不明显, 但是对微生物的活动有主要的影响, 从而间接影响着土壤呼吸。马文瑛[29]认为在整个生长季0—60cm土壤水热状况是主要的影响因素, 其次是土壤表层的环境因子 (温、湿度以及风速的影响), 最后为太阳辐射。在图 4中, 按照相关系数最大的滞后或者前移时间, 绘制了人工杨树林保留和去除枯枝落叶土壤呼吸与土壤温度的指数关系, 发现保留枯枝落叶的土壤呼吸均达到显著水平, 而去除枯枝落叶的土壤呼吸均未达到显著水平, 这是由于枯枝落叶有保温和吸持水的能力[37-38], 塑造良好的呼吸环境, 且可以减弱外界因素 (如风速、人为因素等) 对土壤呼吸的干扰。

3.2 土壤呼吸的月变化特征与土壤温度和土壤含水量的相关性

根据相关性分析, 分别用指数模型对人工杨树林保留和去除枯枝落叶的土壤呼吸和土壤温度进行了模拟分析, 用幂函数模型对小叶锦鸡儿群落保留和去除枯枝落叶的土壤呼吸与土壤含水量进行了模拟分析。温度是影响土壤呼吸的一个非常重要环境因子[21], 主要通过影响微生物活性以及植物根系生长进而影响土壤呼吸[30]。人工杨树林表层土壤由于枯枝落叶的常年积累, 有机质含量高;且地下水位埋深浅, 土壤含水量也较丰富;并且林地的根系深度、密度在0—50cm垂直范围内较集中[31], 易受温度的影响, 保留和去除枯枝落叶土壤呼吸随着温度的升高而增加, 观测日8月20日达到最大, 其后随着温度的降低而减少。土壤含水量是影响土壤呼吸的另一关键因子, 通过影响根系和微生物的生理过程以及底物和氧气的扩散进而调控土壤呼吸[32]。有研究表明只有当土壤含水量过高或过低时才会对土壤呼吸有显著影响[33, 34], 小叶锦鸡儿群落保留和去除枯枝落叶土壤呼吸与土壤含水量显著相关, 受到含水量的胁迫。固定沙丘小叶锦鸡群落土壤含水量较低, 加之地下水位埋深较深, 生长季的土壤含水量主要受降雨的影响。在植被整个生长季人工杨树林土壤呼吸大于小叶锦鸡儿群落土壤呼吸, 主要是由于人工杨树林土壤有机质含量、植被根系量等更大, 为土壤呼吸提供了更丰富的碳底物, 并且小叶锦鸡儿群落还受到土壤含水量胁迫。

运用归一化方法, 建立的土壤呼吸ln Rs与土壤温度T、土壤含水量lnθ的双因子回归模型, 土壤温度和土壤含水量共同对人工杨树林和小叶锦鸡儿群落保留和去除枯枝落叶土壤呼吸解释, 其决定系数均大于二者分别对土壤呼吸的决定系数, 表明采用土壤温度、土壤含水量的双因子模型模拟该区域的土壤呼吸有一定可行性。并且这一结果是由于各个影响因子之间是相互促进相互制约的, 在土壤含水量充足、水分不成为限制因子的条件下土壤呼吸和土壤温度呈正相关, 而当土壤含水量成为限制因子时水分含量和温度共同起作用[35]

3.3 枯枝落叶对土壤呼吸的影响

人工杨树林和小叶锦鸡儿群落保留和去除枯枝落叶处理方式下的土壤呼吸均显著相关 (P < 0.01), 枯枝落叶作为土壤生态系统独特的结构层次, 是土壤有机质输入的主要来源, 对土壤温度、湿度产生影响[7], 枯枝落叶对土壤呼吸的贡献并不完全来源于枯枝落叶自身分解释放的CO2, 并且去除枯枝落叶可以改变土壤根系和微生物所需有机碳的供应[31]。小叶锦鸡儿群落保留和去除枯枝落叶对土壤呼吸的影响结果表明, 去除枯枝落叶会显著降低土壤呼吸, 这与很多研究者的结论相同[11, 30, 32]。整个生长季枯枝落叶对土壤呼吸的贡献率为24.22%, 比陈光水等分析国内62个森林样地土壤呼吸及其相关因子数据得出枯枝落叶对土壤呼吸的贡献率为20.2%[36]相近略高, 这是由于研究区小叶锦鸡儿群落生长在固定沙丘上, 土壤有机质含量低, 土地贫瘠, 根系多分布在土壤深层, 在这种条件下枯枝落叶便成为土壤呼吸的主要来源。保留枯枝落叶土壤呼吸与10cm土壤含水量相关性更好, 并且幂函数关系更为密切, 这是由于枯枝落叶具有吸持、拦截水量, 减少地表蒸发, 改善土壤结构等作用[37]。人工杨树林保留和去除枯枝落叶土壤呼吸的研究结果与小叶锦鸡儿群落相反, 枯枝落叶对土壤呼吸的贡献率为-25.60%, 枯枝落叶对土壤呼吸的贡献为负值这一结果与陈四清对内蒙古锡林河流域大针茅草原[18]以及其他学者对山地雨林和退耕还湿地的研究结果一致[19-20]。对这一现象的解释大多集中为可能与枯枝落叶数量、土地利用方式、土壤湿润因子有关[19], 可能是土壤表层好氧性细菌因去除枯枝落叶异常活跃导致呼吸增强或枯枝落叶对土壤呼吸CO2排放的屏蔽作用[18]。孟庆涛在科尔沁沙地研究发现枯枝落叶只有当温度达到40℃时, 枯枝落叶才会向土壤及大气等环境中输入大量的碳[39], 在研究区最高温度小于40℃, 所以枯枝落叶本身对土壤呼吸的贡献很小, 而枯枝落叶对土壤温度的影响主要通过阻隔土壤与外界空气的热交换, 枯枝落叶量越多对土壤温度的影响越大[38]。在研究区人工杨树林去除枯枝落叶使土壤直接裸露被太阳直射, 土壤温度升高, 加之地表有5cm腐殖质层, 从而使土壤呼吸速率变大。去掉枯枝落叶势必会降低土壤含水量, 从而导致土壤呼吸降低[32], 但人工杨树林去掉枯枝落叶土壤呼吸与土壤含水量的关系不显著, 故受其影响不明显。王丽丽在三江平原发现屏蔽作用随着枯枝落叶输入量的增加而增大, 人工杨树林枯枝落叶 (厚1.5cm) 的输入量与其对土壤呼吸的贡献率呈负相关关系[20], 在我们的研究区人工杨树林枯枝落叶 (厚1.6cm) 对土壤呼吸的贡献率为负, 说明可能是枯枝落叶的屏蔽作用导致。Gadgil[40]的研究发现菌根真菌不从枯枝落叶中获取能量, 而是在枯枝落叶中进行休眠, 并通过利用根的分泌物生长, 抑制了土壤微生物的活性以及整个生态系统的碳循环速率[41]。随后, Olsson等也发现菌根真菌的菌丝降低了砂土中的细菌活性[42]。枯枝落叶的分解速率会受其氮浓度和磷浓度的影响[43], 因此菌根真菌还能通过改变根系凋落物的化学计量比来改变枯枝落叶的分解速率。具有外生菌根真菌侵染的植物根系比起没被侵染的根系有更多的氮[41], 这本应使之具有更高的分解速率, 但是这些氮是以难分解的形式存在的, 因此氮增加的净效应或许会降低枯枝落叶分解速率, Langley等证实了这种解释[44]。由此可见, 枯枝落叶对土壤呼吸的影响非常复杂, 不能单从宏观方面进行分析, 要深入到微生物对枯枝落叶的作用。人工杨树林保留枯枝落叶土壤呼吸与土壤温度的指数关系更为密切, 并且保留和去除枯枝落叶的Q10值分别为2.97、2.60, 比较接近中国森林土壤呼吸Q10的平均值2.65[36], 枯枝落叶的去除降低了土壤呼吸的温度敏感性, 这一结论与邓琦[45]、熊莉[32]的研究结论一致, 枯枝落叶位于土表面, 能有效降低土壤温度受气温波动的影响, 进而减缓土壤环境对外界环境尤其是全球变暖响应的敏感性[45]

4 结论

科尔沁沙丘-草甸相间地区沙丘-草甸过渡带人工杨树林土壤呼吸日变化, 生长期的8月20日和生长末期的10月2日均呈“多峰型”, 保留和去除枯枝落叶土壤呼吸均出现多次峰值。人工杨树林和固定沙丘小叶锦鸡儿群落保留和去除枯枝落叶土壤呼吸月变化波动较大, 保留和去除枯枝落叶土壤呼吸的变化趋势在5—10月表现同步, 人工杨树林去除枯枝落叶土壤呼吸大于保留枯枝落叶土壤呼吸, 而小叶锦鸡儿群落去除枯枝落叶土壤呼吸小于保留枯枝落叶土壤呼吸。土壤呼吸采用土壤温度、土壤含水量双因子回归模型取得的模拟效果要比单因子回归模型好, 双因子回归模型能更好地模拟该区域土壤呼吸对环境因子变化的响应, 人工杨树林土壤呼吸主要受土壤温度的控制, 且保留枯枝落叶土壤呼吸的Q10比去除枯枝落叶土壤呼吸的Q10大, 小叶锦鸡儿群落土壤呼吸主要受土壤含水量的控制。本试验只有一年的数据具有一定的局限性, 枯枝落叶对土壤呼吸的影响还需要在微生物方面加深研究, 在接下来的实验中需野外实验和室内实验相结合才能更清晰的把握枯枝落叶对土壤呼吸的影响机理。

致谢: 感谢内蒙古农业大学张圣微, 姚姣转老师的帮助。
参考文献
[1] 张丽华, 陈亚宁, 李卫红, 赵锐锋, 葛洪涛. 干旱区荒漠生态系统的土壤呼吸. 生态学报, 2008, 28(5): 1911–1922.
[2] Jenkinson D S, Adams D E, Wild A. Model estimates of CO2 emissions from soil in response to global warming. Nature, 1991, 351(6324): 304–306. DOI:10.1038/351304a0
[3] Schlesinger W H, Andrews J A. Soil respiration and the global carbon cycle. Biogeochemistry, 2000, 48(1): 7–20. DOI:10.1023/A:1006247623877
[4] Reynolds J F. Desertification. In: Levin, S. A. ed, Encyclopedia of Biodiversity, vol 2. San Diego: Academic Press, 2001, 61-78.
[5] Conant R T, Klopatek J M, Klopatek C C. Environmental factors controlling soil respiration in three semiarid ecosystems. Soil Science Society of America Journal, 2000, 64(1): 383–390. DOI:10.2136/sssaj2000.641383x
[6] West N E, Stark J M, Johnson D W, Abrams M M, Wight J R, Heggem D, Peck S. Effects of climatic change on the edaphic features of arid and semiarid lands of western North America. Arid Soil Research and Rehabilitation, 1994, 8(4): 307–351. DOI:10.1080/15324989409381408
[7] 张东秋, 石培礼, 张宪洲. 土壤呼吸主要影响因素的研究进展. 地球科学进展, 2005, 20(7): 778–785.
[8] 鲍芳, 周广胜. 中国草原土壤呼吸作用研究进展. 植物生态学报, 2010, 34(6): 713–726.
[9] 陈书涛, 刘巧辉, 胡正华, 刘艳, 任景全, 谢薇. 不同土地利用方式下土壤呼吸空间变异的影响因素. 环境科学, 2013, 34(3): 1017–1025.
[10] Raich J W, Potter C S. Global patterns of carbon dioxide emissions from soils. Global Biogeochemical Cycles, 1995, 9(1): 23–36. DOI:10.1029/94GB02723
[11] Sulzman E W, Brant J B, Bowden R D, Lajtha K. Contribution of aboveground litter, belowground litter, and rhizosphere respiration to total soil CO2 efflux in an old growth coniferous forest. Biogeochemistry, 2005, 73(1): 231–256. DOI:10.1007/s10533-004-7314-6
[12] Prescott C E. Do rates of litter decomposition tell us anything we really need to know?. Forest Ecology and Management, 2005, 220(1/3): 66–74.
[13] Wan S Q, Norby R J, Ledford J, Weltzin J F. Responses of soil respiration to elevated CO2, air warming, and changing soil water availability in a model old-field grassland. Global Change Biology, 2007, 13(11): 2411–2424. DOI:10.1111/gcb.2007.13.issue-11
[14] Li Y Q, Xu M, Sun O J, Cui W C. Effects of root and litter exclusion on soil CO2 efflux and microbial biomass in wet tropical forests. Soil Biology and Biochemistry, 2004, 36(12): 2111–2114. DOI:10.1016/j.soilbio.2004.06.003
[15] Raich J W, Nadelhoffer K J. Belowground carbon allocation in forest ecosystems: global trends. Ecology, 1989, 70(5): 1346–1354. DOI:10.2307/1938194
[16] 田祥宇, 涂利华, 胡庭兴, 张健, 何远洋, 肖银龙. 华西雨屏区苦竹人工林土壤呼吸各组分特征及其温度敏感性. 应用生态学报, 2012, 23(2): 293–300.
[17] 周小刚, 郭胜利, 车升国, 张芳, 邹俊亮, 张彦军, 南雅芳, 李泽, 蒲辉. 黄土高原刺槐人工林地表凋落物对土壤呼吸的贡献. 生态学报, 2012, 32(7): 2150–2157.
[18] 陈四清, 崔骁勇, 周广胜, 李凌浩. 内蒙古锡林河流域大针茅草原土壤呼吸和凋落物分解的CO2排放速率研究. 植物学报, 1999, 41(6): 645–650.
[19] 骆土寿, 陈步峰, 李意德, 林明献, 周光益, 陈德祥, 邱坚锐. 海南岛尖峰岭热带山地雨林土壤和凋落物呼吸研究. 生态学报, 2001, 21(12): 2013–2017.
[20] 王丽丽, 宋长春, 郭跃东, 刘德燕, 杨桂生. 三江平原不同土地利用方式下凋落物对土壤呼吸的贡献. 环境科学, 2009, 30(11): 3130–3135.
[21] Luo Y Q, Wan S Q, Hui D F, Wallace L L. Acclimatization of soil respiration to warming in a tall grass prairie. Nature, 2001, 413(6856): 622–625. DOI:10.1038/35098065
[22] Rey A, Pegoraro E, Tedeschi V, De Parri I, Jarvis P G, Valentini R. Annual variation in soil respiration and its components in a coppice oak forest in central Italy. Global Change Biology, 2002, 8(9): 851–866. DOI:10.1046/j.1365-2486.2002.00521.x
[23] 蒋延玲, 周广胜, 赵敏, 王旭, 曹铭昌. 长白山阔叶红松林生态系统土壤呼吸作用研究. 植物生态学报, 2005, 29(3): 411–514.
[24] 高翔, 郝卫平, 顾峰雪, 郭瑞, 夏旭, 梅旭荣, 李洁. 降雨对旱作春玉米农田土壤呼吸动态的影响. 生态学报, 2012, 32(24): 7883–7893.
[25] 李虎, 邱建军, 王立刚. 农田土壤呼吸特征及根呼吸贡献的模拟分析. 农业工程学报, 2008, 24(4): 14–20.
[26] Liu Q, Edwards N T, Post W M, Gu L, Ledford J, Lenhart S. Temperature-independent diel variation in soil respiration observed from a temperate deciduous forest. Global Change Biology, 2006, 12(11): 2136–2145. DOI:10.1111/gcb.2006.12.issue-11
[27] 方晰, 田大伦. 杉木人工林林地CO2释放量的研究. 林业科学, 1997, 33(sp. 2): 94–103.
[28] Yemadje P L, Guibert H, Blavet D, Olina J P, Chevallier T, Deleporte P, Bernoux M. 24-h variation in soil respiration after a long dry season in a Sudano-Sahelian region. Scientia Agricola, 2015, 72(5): 452–558. DOI:10.1590/0103-9016-2014-0412
[29] 马文瑛, 赵传燕, 彭守璋, 高云飞, 袁亚鹏, 李文娟. 黑河天涝池五种植被类型土壤呼吸速率动态特征及其影响因子. 生态学报, 2015, 35(17): 5654–5665.
[30] 王光军, 田大伦, 闫文德, 朱凡, 项文化, 梁小翠. 改变凋落物输入对杉木人工林土壤呼吸的短期影响. 植物生态学报, 2009, 33(4): 739–747.
[31] 陈全胜, 李凌浩, 韩兴国, 阎志丹, 王艳芬, 张焱, 袁志友, 唐芳. 温带草原11个植物群落夏秋土壤呼吸对气温变化的响应. 植物生态学报, 2003, 27(4): 441–547.
[32] 熊莉, 徐振锋, 杨万勤, 殷睿, 唐仕姗, 王滨, 徐李亚, 常晨晖. 川西亚高山粗枝云杉人工林地上凋落物对土壤呼吸的贡献. 生态学报, 2015, 35(14): 4678–5686.
[33] Davidson E A, Verchot L V, Cattânio J H, Ackerman I L, Carvalho J E M. Effects of soil water content on soil respiration in forests and cattle pastures of eastern Amazonia. Biogeochemistry, 2000, 48(1): 53–69. DOI:10.1023/A:1006204113917
[34] Maier C A, Kress L W. Soil CO2 evolution and root respiration in 11 year-old loblolly pine (Pinus taeda) plantations as affected by moisture and nutrient availability. Canadian Journal of Forest Research, 2000, 30(3): 347–359. DOI:10.1139/cjfr-30-3-347
[35] 刘绍辉, 方精云. 土壤呼吸的影响因素及全球尺度下温度的影响. 生态学报, 1997, 17(5): 469–576.
[36] 陈光水, 杨玉盛, 吕萍萍, 张亿萍, 钱小兰. 中国森林土壤呼吸模式. 生态学报, 2008, 28(4): 1748–1761.
[37] Atkin O K, Edwards E J, Loveys B R. Response of root respiration to changes in temperature and its relevance to global warming. New Phytologist, 2000, 147(1): 141–154. DOI:10.1046/j.1469-8137.2000.00683.x
[38] 潘开文, 何静, 吴宁. 森林凋落物对林地微生境的影响. 应用生态学报, 2004, 15(1): 153–158.
[39] 孟庆涛, 李玉霖, 赵学勇, 赵玉萍, 罗亚勇. 科尔沁沙地不同环境条件下植物叶凋落物CO2释放研究. 干旱区研究, 2008, 25(4): 519–524.
[40] Gadgil R L, Gadgil P D. Mycorrhiza and litter decomposition. Nature, 1971, 233(5315): 133–133.
[41] Langley J A, Dijkstra P, Drake B G, Hungate B A. Ectomycorrhizal colonization, biomass, and production in a regenerating scrub oak forest in response to elevated CO2. Ecosystems, 2003, 6(5): 424–530. DOI:10.1007/PL00021509
[42] Olsson P A, Chalot M, Bååth E, Finlay R D, Söderström B. Ectomycorrhizal mycelia reduce bacterial activity in a sandy soil. FEMS Microbiology Ecology, 1996, 21(2): 77–86. DOI:10.1111/fem.1996.21.issue-2
[43] Liu P, Sun O J, Huang J H, Li L H, Han X G. Nonadditive effects of litter mixtures on decomposition and correlation with initial litter N and P concentrations in grassland plant species of northern China. Biology and Fertility of Soils, 2007, 44(1): 211–216. DOI:10.1007/s00374-007-0195-9
[44] Langley J A, Chapman S K, Hungate B A. Ectomycorrhizal colonization slows root decomposition: the post-mortem fungal legacy. Ecology Letters, 2006, 9(8): 955–959. DOI:10.1111/ele.2006.9.issue-8
[45] 邓琦, 刘世忠, 刘菊秀, 孟泽, 张德强. 南亚热带森林凋落物对土壤呼吸的贡献及其影响因素. 地球科学进展, 2007, 22(9): 976–986.