徐 洲,冯 倩,王 玉,赵金磊,李常鑫,王丽梅,2,*
1 西北农林科技大学资源环境学院,杨凌 712100 2 农业部西北植物营养与农业环境重点实验室,杨凌 712100
工业革命以来,大气 CO2浓度持续升高已是不争的事实。到 21 世纪末,大气中 CO2浓度将可能达到甚至超过 700 μmol/mol。农田土壤呼吸是大气 CO2排放的主要来源之一[1-3],全球每年因土壤呼吸排放到大气中 CO2约为70—100 Pg,是化石燃料燃烧排放 CO2的12—16倍[4]。目前已有研究表明[5-6],大气 CO2浓度升高会对农田作物土壤呼吸产生一定影响。因此,土壤呼吸对大气 CO2浓度升高具有重要贡献,反之,升高了的大气 CO2浓度又势必会对农田土壤温室气体排放产生反馈作用,探明大气 CO2浓度升高条件下农田土壤呼吸速率的改变,对于明确未来气候变化条件下,农田生态系统对气候改变的响应与反馈具有重要意义。
土壤呼吸是陆地生态系统碳循环的重要环节,是陆地生态系统中碳以 CO2形式返回到大气的主要途径[7]。研究表明[8-9],植物通过光合作用将大气中的 CO2以有机碳的形式储存到地下,其中部分以枯枝落叶形式进入土壤,成为生态系统碳储量的一部分[10]。全球土壤碳储量的增加有利于减缓温室效应,而土壤呼吸的增加会加剧全球温室效应,造成一系列不可预计的后果。可见,土壤呼吸在全球碳收支中扮演着重要的角色,其微小的变化都可能给全球碳平衡带来重大影响。
农田生态系统是陆地生态系统的重要组成部分[11-13],研究农田生态系统中土壤呼吸对大气 CO2浓度增加的响应是预测和评价农田系统乃至整个陆地生态系统土壤碳周转和碳收支的重要前提与基础。玉米是全世界分布最广的农作物之一[14-16],也是我国重要的粮食作物。开展大气 CO2浓度升高及土壤温度和水分对田间玉米土壤呼吸速率的影响研究,有助于阐明农田土壤呼吸对大气 CO2浓度升高的响应和反馈。因此,本试验以玉米农田土壤为研究对象,通过改进的开顶式气室(OTC)系统自动控制大气 CO2浓度,在春玉米播前(V0)、六叶期(V6)、九叶期(V9)、吐丝期(R1)、乳熟期(R3)、蜡熟期(R5)及完熟期(R6)对土壤呼吸速率进行测定,研究玉米农田土壤呼吸对大气 CO2浓度升高的响应,以及与土壤温度、水分的互馈效应,以期为未来气候条件下农田生态系统碳循环提供基础数据,为全球气候变化背景下的农田土壤呼吸和碳固定及全球碳循环研究提供帮助。
本试验地点位于长武黄土高原农业生态试验站(35°12′N,107°40′E),该站海拔 1200 m,年均降水量为 584 mm,年均气温 9.1℃,5—9 月平均气温为 19.0℃,无霜期 171 d,地下水位 50—80 m,年蒸发量高达 1565 mm,属典型旱作农业区。地带性土壤为黑垆土,质地均匀疏松。播前耕层土壤基本理化性质如下:pH 值为8.0,容重 1.28 g/cm3,有机质 12.46 g/kg,全氮 0.98 g/kg,速效磷(Olsen-P)22.2 mg/kg,速效钾(NH4OAc-K)181.9 mg/kg,矿质氮 16.9 mg/kg。
1.2.1处理设置
本试验采用田间定位试验,共设 3 种处理:大田自然大气条件下 CO2浓度处理(CK)、设置气室但不通入 CO2使其浓度与 CK 保持一致的 OTC对照处理(OTC)及 OTC系统自动控制 CO2浓度(700 μmol/mol,OTC+CO2)处理,其中设置 OTC 对照处理是为了去除 OTC 气室可能带来的增温影响。每个处理重复 3 次,共设 9 个小区。各小区以含氮量为 46%的尿素为氮源,施纯氮 225 kg/hm2,在播前、喇叭口期(V10)和抽雄期(VT)分三次施入,施肥比为 4∶3∶3;在播前将过磷酸钙(含 P2O516%)、硫酸钾(含 K2O 50%)作为底肥一次性全部施入,施入量分别为 40、80 kg/hm2。试验种植春玉米品种为“郑单 958”,种植密度 70000株/hm2,采用双垄沟播种技术,其中宽窄行分别为 60 cm和 40 cm,低高垄分别为5 cm和 15 cm。于 2018 年 4 月 27 日播种,9月13 日收获。
玉米播种后,在每个小区中间埋设直径 20 cm、高 10 cm的用 PVC 材料制成的塑料环,将塑料环垂直插入地面深 5 cm,露出地面的部分(5 cm)使用铁丝环加地膜将其封闭,测定时取下[17]。
1.2.2土壤呼吸速率测定
春玉米从每年 4 月中旬左右开始播种,9 月中旬收获。土壤呼吸测定采用 LI-8100A 型土壤碳通量自动测量系统(美国 LI-COR 公司)田间原位测定。土壤呼吸每个生育期的测定,开始前晴天天气需持续3 d以上,从而剔除降雨因素的影响,测定时间从 9:00—11:00。该仪器在测定田间土壤呼吸的同时,仪器自身所带的温度探针会同步测定地下 5 cm的温度和环外土壤含水量。土壤呼吸速率和地下 5 cm的土温、土壤含水率可由 LI-8100A 土壤碳通量分析仪在电脑上直接导出。
数据整理和图表分析采用Excel 2016 和 Origin 8.5。采用 SPSS 20.0统计软件进行单因素方差(One-way ANOVA)检验、重复测量方差分析(Repeated-analysis ANOVA)和皮尔逊相关性分析(Univariate),采用 Duncan 法在P<0.05水平作各生育期处理间多重比较。
本试验3种处理条件下春玉米农田土壤呼吸速率自 V0 到 R6 期皆呈降-升-降趋势,其最大值出现在玉米生育中后期(图 1)。整个生育期,CK、OTC、OTC+CO23种处理条件下土壤呼吸速率平均值分别为16.59、17.36、26.64 μmol m-2s-1。V0 到 V6 期为土壤呼吸下降阶段,最低值出现在 V6 生育期,CK、OTC、OTC+CO23种处理条件下土壤呼吸速率最低值分别为 11.41、9.88、18.95 μmol m-2s-1。随后从 V6 期开始土壤呼吸速率逐渐上升并达到峰值,不同处理达到峰值时间不同,CK 在 R3 期达到峰值,为 22.57 μmol m-2s-1;OTC 在 R1 期达到峰值,为 21.04 μmol m-2s-1;OTC+CO2在 R5 期达到峰值,为 34.88 μmol m-2s-1。之后各处理土壤呼吸速率开始下降,直至 R6 期结束,CK、OTC、OTC+CO23种处理条件下土壤呼吸速率 R6 期分别为 9.79、10.96、25.07 μmol m-2s-1。
图1 三种处理条件下春玉米各生育期的土壤呼吸 Fig.1 Soil respiration at different growth stage of spring maize under three atmospheric conditionsCK:自然大气 Natural atmospheric conditions with the current CO2 concentration;OTC:开顶式气室对照 Open-top chamber system with the current atmospheric CO2 concentration;OTC+CO2:CO2浓度升高处理 OTC system with elevated CO2 concentration;V0:播前 Pre-sowing;V6:六叶期 Jointing stage;V9:九叶期 Ninth-leaf stage;R1:吐丝期 Silking stage;R3:乳熟期 Milk stage;R5:蜡熟期 Dent stage;R6:完熟期 Physiological maturity stage
相较于 CK,OTC+CO2处理土壤呼吸速率在 R3 和 R5 期分别增加 43%和 104% (P<0.05),其他生育期影响不显著。而 OTC 相较于 CK,两种处理土壤呼吸速率各生育期均无显著差异,且变化趋势基本一致,仅在播前土壤呼吸速率差异较大,CK 为 15.22 μmol m-2s-1,OTC 为 22.43 μmol m-2s-1。相较于 OTC,OTC+CO2处理土壤呼吸速率在 R3 和 R5 期分别增加 63%和 109% (P<0.05),其他生育期差异不显著。
由图2可知,CK、OTC和OTC+CO23种处理条件下土壤温度随生育期变化趋势基本一致,均呈降-升-降趋势,V6 期温度最低,R3 期温度达到峰值,之后土壤温度下降直至 R6 期。整个生育期,3种处理之间温度差异较小,在 0.21—1.94℃之间。CK、OTC和OTC+CO23种处理土壤温度变化范围分别为 16.76—32.57℃、15.32—31.25℃、18.21—31.74℃。
图2 3种处理条件下春玉米各生育期的土壤温度 Fig.2 Soil temperature at different growth stage of spring maize under three atmospheric conditions
将农田春玉米各生育期土壤温度和土壤呼吸速率用回归方法进行统计分析得到图 3,结果表明,CK、OTC和OTC+CO23种处理条件下农田春玉米土壤呼吸和土壤温度呈指数相关,皆表现出显著性,3种处理间显著性大小表现为 CK >OTC >OTC+CO2。3种处理条件下土壤温度敏感系数Q10如表1 所示,CK 和OTC 处理条件下的Q10值相差不大,分别为 1.52 和 1.54,而 OTC+CO2处理条件下的Q10值则较小,为 1.03,相较于 CK 和 OTC,OTC+CO2处理Q10值分别降低了 32%和 33%。CK 和 OTC 处理条件下的a值(0℃的呼吸速率)差异不大,分别为 5.94 和 5.92 μmol m-2s-1,而 OTC+CO2处理条件下的a值高达 13.95 μmol m-2s-1。
表1 土壤呼吸与土壤温度指数相关系数a和Q10Table 1 Correlation coefficients a and Q10 between soil respiration and soil temperature
图3 3种处理条件下土壤呼吸与土壤温度的关系Fig.3 Relationship between soil respiration and soil temperature under three atmospheric conditions
春玉米各生育期土壤含水量如图 4 所示,从 V0 期到 V6 期,OTC 和 OTC+CO2两种处理条件下土壤水分有所下降,CK 土壤水分有所上升。除去 V0 期,在其他玉米生育期,CK、OTC 和 OTC+CO23种处理条件下土壤水分变化趋势基本一致,呈先升后降趋势。春玉米土壤水分从 V6 期开始持续上升,R3 期达到峰值,之后土壤水分下降直至 R6 期。整个春玉米生育期,3种处理之间土壤水分无明显差异。CK、OTC 和 OTC+CO23种处理土壤水分变化范围分别为 19.84%—39.29%、20.46%—37.2%、19.63%—38.14%。
图4 3种处理条件下春玉米各生育期的土壤水分 Fig.4 Soil moisture at different growth stage of spring maize under three atmospheric conditions
将农田春玉米各生育期土壤水分和土壤呼吸速率用回归方法进行统计分析得到图 5,结果表明,CK、OTC 和 OTC+CO23种处理条件下农田春玉米土壤呼吸和土壤水分皆呈抛物线型相关,表现出显著性(P<0.05)。3种处理间土壤水分对土壤呼吸速率影响大小分别为 OTC >OTC+CO2>CK。
图5 3种处理条件下土壤呼吸与土壤水分的关系Fig.5 Relationship between soil respiration and soil moisture under three atmospheric conditions
农田生态系统作为陆地主要生态系统之一,研究未来气候变化对其产生的影响已有 50多年的历史[18-21],农田土壤呼吸的微小变动都会对全球碳循环产生重要影响,针对 CO2浓度升高对农田土壤呼吸及其环境因子影响的研究受到越来越多关注。
目前,国内外大部分研究表明大气 CO2浓度升高促进了土壤呼吸速率,Bachman等[22]在研究 CO2浓度升高对落叶松土壤呼吸的影响时发现,将 CO2浓度提升至 550 μmol/mol,白桦地土壤呼吸速率升高了 15%—70%。Bernhardt等[23]利用 FACE 平台将 CO2浓度提升至 700 μmol/mol研究森林土壤呼吸速率的变化,发现相比于 CK,大气 CO2浓度提升为 700 μmol/mol的试验组在前两年土壤呼吸速率显著提升了 30%—34%。但是也有研究发现 CO2浓度升高对土壤呼吸没有影响,甚至抑制了土壤呼吸。例如,马红亮等[24]在研究 CO2浓度升高对水稻土壤呼吸的影响时发现,大气 CO2浓度增加降低了土壤呼吸排放速率。同时,有研究发现[25-26],CO2浓度升高对农田土壤呼吸的影响因作物生长期不同而存在明显差异。Soe等[27]在对甜菜进行研究时发现,将大气 CO2浓度从 380 μmol/mol增加到 550 μmol/mol,甜菜作物在最旺盛生长阶段土壤呼吸速率显著增加了 34%,而在甜菜其它生长阶段则没有达到显著增加。寇太记等[28]在对冬小麦的研究中也发现,CO2浓度升高对冬小麦土壤呼吸速率的影响会随着冬小麦生长阶段不同而呈现明显的季节性差异,在冬小麦生长旺盛期,农田土壤呼吸速率达到最大。本试验研究了春玉米整个生育期 CO2浓度升高对土壤呼吸的影响后发现,将 CO2浓度提升至 700 μmol/mol,相较于对照组(OTC),试验组(OTC+CO2)土壤呼吸在玉米生殖生长为中心的阶段 R3 和 R5 期分别显著提升 63%和 109% (P<0.05)。相较于 CK,试验组(OTC+CO2)土壤呼吸在玉米生长最旺盛的 R3 和 R5 期分别显著提升 43%和 104% (P<0.05)。这与上述研究结果基本一致,表明 CO2浓度升高对土壤呼吸的影响随着作物生育期的不同而存在明显差异,这可能与 CO2浓度升高影响了植物根系生长有关。我们相关田间结果表明,R3 和 R5 期玉米根系相对发达,OTC+CO2处理玉米根系生物量(R3 和 R5 期分别为 41.38 g/m2、49.52 g/m2)显著高于 OTC 处理(R3 和 R5 期分别为 31.42 g/m2、35.46 g/m2)与 CK(R3和R5期分别为 33.07 g/m2、34.42 g/m2)。韩广轩等[29]在对玉米土壤呼吸的研究中也发现,在空间尺度上,土壤呼吸作用与根系生物量呈显著的线性关系。这表明高 CO2浓度是通过促进植物根系生长来影响土壤呼吸。
土壤温度是影响土壤呼吸速率的重要环境因素[30-32]。土壤温度主要是通过对土壤有机碳的分解和作物根系自养呼吸的影响来影响土壤呼吸,前人研究表明[33-34],在土壤水分保持在田间最大持水量 50%—80%左右时,土壤呼吸速率随着土壤温度的增加而升高,呈指数型趋势。本研究将 CK、OTC 和 OTC+CO23种处理条件下农田春玉米土壤呼吸和土壤温度进行拟合后发现,无论是 CK、OTC,还是 OTC+CO2,土壤呼吸和土壤温度皆呈现出显著的指数关系,表明土壤温度是影响土壤呼吸速率的重要环境因子。温度敏感系数Q10表示气温每增加 10℃土壤呼吸增加的倍数,OTC+CO2处理条件下的Q10值均低于 CK 和 OTC,表明大气 CO2浓度升高降低了土壤呼吸对土壤温度的敏感性。且 OTC+CO2处理条件下的a值高达 13.95 μmol m-2s-1,这可能是由于大气 CO2浓度升高和土壤温度变化两者间互相影响所致。
土壤水分主要通过调控土壤通透性及地下微生物的生理活动来影响土壤呼吸[35-36]。前人研究表明[37-39],在土壤水分充足的条件下,土壤呼吸与土壤水分呈正相关,表现为抛物线型函数、线性函数、指数函数等形式。本研究通过对 CK、OTC 和 OTC+CO23种处理条件下农田春玉米土壤呼吸和土壤水分进行拟合后表明,3种处理土壤呼吸与土壤水分皆满足抛物线型关系,表现出显著性(P<0.05)。可知土壤水分是影响土壤呼吸速率的重要环境因子。土壤水分对土壤呼吸的影响具有一定的调节范围。当土壤处于缺水状态时,微生物的活动减弱,土壤呼吸减少,这时土壤水分的增加会促进土壤呼吸。土壤含水量过高时,土壤空隙被堵塞,土壤中氧气含量减少,抑制微生物有氧呼吸从而降低土壤呼吸。因此,必然存在一个阈值,让土壤呼吸开始受到土壤水分的抑制。根据抛物线函数的特性可知,CK、OTC 和 OTC+CO23种处理土壤呼吸对土壤水分响应的阈值分别为 36.1%、33.5%、41.2%。与 CK、OTC 相比,OTC+CO2处理土壤呼吸对土壤水分响应的阈值分别提升了 14%和 23%。
基于 CK、OTC 和 OTC+CO23种处理条件下农田土壤温度、水分和呼吸数据,就土壤温度和水分与土壤呼吸的相关性进行皮尔逊分析如表2 所示。分析表明,CK 处理土壤温度和水分与土壤呼吸均呈极显著相关,OTC 处理土壤温度和水分与土壤呼吸均呈显著相关,而 OTC+CO2处理土壤温度和水分与土壤呼吸均未达显著相关水平。CK、OTC 处理条件下,土壤温度对土壤呼吸的影响程度均大于土壤水分。OTC+CO2处理土壤温度和水分与土壤呼吸均未达显著相关水平。表明大气 CO2浓度升高导致土壤温度和土壤水分对土壤呼吸解释的信息量减少,不确定因素增多。
表2 土壤温度和水分与土壤呼吸的相关性分析Table 2 Correlation of soil temperature and moisture with soil respiration
大气 CO2浓度升高、土壤水分和土壤温度都会对土壤呼吸产生一定的影响。那么,大气 CO2浓度升高是否会对土壤呼吸的水分与温度效应产生影响,即大气 CO2浓度升高与土壤温度和水分对土壤呼吸影响的耦合作用是什么呢?目前有关这方面的研究较少,寇太记等[40]在研究 CO2浓度升高对冬小麦土壤呼吸速率的影响时发现,CO2浓度升高会降低土壤呼吸对温度的敏感性。本研究也发现,与 CK、OTC 相比,OTC+CO2处理条件下Q10分别降低了 12%和 9%。此外,OTC+CO2处理条件下的a值是其他处理条件下的两倍多。分析原因可能是,CO2浓度升高促进了玉米根系的生长和土壤碳的固定,增强了同温度条件下的土壤呼吸。本研究表明,大气 CO2浓度升高确实会使土壤温度和水分对土壤呼吸产生一定的耦合作用,但该耦合作用的机理、限制因素和影响程度还需要进一步的研究加以探讨。
(1)CO2浓度升高对土壤呼吸的影响随春玉米生育期改变,在 R3 和 R5 期明显提高了土壤呼吸作用。
(2)土壤温度和水分均影响土壤呼吸,分别呈现出指数相关性和抛物线型相关性。相同条件下,土壤温度对土壤呼吸的影响大于土壤水分。
(3)大气 CO2浓度升高降低了土壤呼吸对土壤温度的敏感性,提高了土壤呼吸对土壤水分响应的阈值。