环境科学  2022, Vol. 43 Issue (1): 569-576   PDF    
秸秆还田对冬小麦-夏玉米农田土壤固碳、氧化亚氮排放和全球增温潜势的影响
万小楠1, 赵珂悦1, 吴雄伟1, 白鹤1, 杨学云1, 顾江新1,2     
1. 西北农林科技大学资源环境学院, 杨凌 712100;
2. 农业农村部西北植物营养与农业环境重点实验室, 杨凌 712100
摘要: 旱地农田温室气体净排放(以全球增温潜势表示)主要取决于土壤固碳速率和氧化亚氮(N2O)排放量.基于长期定位施肥试验,综合分析2010~2017年表层(0~20 cm)土壤有机碳含量和2014~2017年N2O排放通量的观测结果,定量评价秸秆还田对关中平原冬小麦-夏玉米农田土壤固碳速率、N2O年排放量和全球增温潜势的影响.田间试验设置传统施肥(CF)和传统施肥加玉米秸秆(CFS)处理及不施肥对照(CK).CF和CFS处理施入等量尿素,冬小麦季和夏玉米季施氮量分别为165 kg·hm-2和188 kg·hm-2;CF处理在两季作物收获时均保留本小区地上约10 cm高根茬,CFS处理在冬小麦收获时保留相同根茬,而在玉米收获时保留全部秸秆(含氮量约40 kg·hm-2).CK全年不施化肥,秸秆管理方式与CF处理一致.结果表明,CK土壤有机碳含量变化小、N2]O排放量低,其全球增温潜势变化范围为0.04~0.11 t·(hm2·a)-1.CF和CFS处理土壤有机碳含量均随施肥年限呈线性增长(p < 0.001),其固碳速率分别为0.69 t·(hm2·a)-1和0.97 t·(hm2·a)-1;两处理N2O排放量分别在1.65~5.36 kg·(hm2·a)-1和3.08~7.73 kg·(hm2·a)-1之间,CFS处理N2O年排放量较CF处理偏高43%~94%,但仅在2015~2016年有显著性差异(p < 0.05).CF和CFS处理全球增温潜势分别在-1.95~-0.28 t·(hm2·a)-1和-2.59~-0.35 t·(hm2·a)-1之间,CFS处理3 a累计全球增温潜势较CF处理偏低42%.本研究的冬小麦-夏玉米农田在传统施肥管理方式下已是温室气体的汇;尽管土壤固碳速率和N2O排放量存在消长关系,但全量玉米秸秆还田仍然更有利于温室气体减排.
关键词: 土壤固碳      氧化亚氮(N2O)      全球增温潜势(GWP)      秸秆还田      长期定位施肥      关中平原     
Effects of Stalk Incorporation on Soil Carbon Sequestration, Nitrous Oxide Emissions, and Global Warming Potential of a Winter Wheat-Summer Maize Field in Guanzhong Plain
WAN Xiao-nan1 , ZHAO Ke-yue1 , WU Xiong-wei1 , BAI He1 , YANG Xue-yun1 , GU Jiang-xin1,2     
1. College of Natural Resources and Environment, Northwest A & F and Forestry University, Yangling 712100, China;
2. Key Laboratory of Plant Nutrition and the Agri-environment in Northwest China, Ministry of Agriculture and Rural Affairs, Yangling 712100, China
Abstract: The net greenhouse gas emissions from upland soils, as indicated by global warming potential (GWP), mainly depend on the soil carbon sequestration and nitrous oxide (N2O) emissions. The annual changes in surface (0-20 cm) soil organic carbon (SOC) content from 2010 to 2017 and the N2O emissions from 2014 to 2017 were measured within a long-term fertilization experiment. The objective was to quantify the effect of stalk incorporation on the soil carbon sequestration, annual N2O emissions, and GWP of a winter wheat-summer maize field in the Guanzhong Plain. The field experiment included three treatments: conventional fertilization (CF), conventional fertilization plus maize stalks (CFS), and an unfertilized control (CK). The CF and CFS treatments received the same amount of urea per year, with nitrogen (N) input at 165 kg·hm-2 and 188 kg·hm-2 in the winter wheat season and summer maize season, respectively. The CF treatment retained the stubbles (about 10 cm above ground) when harvesting the winter wheat and summer maize crops. The CFS treatment retained the same wheat stubbles and all maize stalks (containing approximately 40 kg·hm-2 of N). The CK treatment was unfertilized throughout the year, with the stubble management the same as that in the CF treatment. The results showed that the CK treatment displayed few changes in SOC content and low N2O emissions, with GWP varying from 0.04 to 0.11 t·(hm2·a)-1. The SOC contents in the CF and CFS treatments increased linearly with the fertilization years (p < 0.001), and their SOC sequestration rates were 0.69 t·(hm2·a)-1 and 0.97 t·(hm2·a)-1, respectively. The N2O emissions from the CF and CFS treatments varied from 1.65 to 5.36 kg·(hm2·a)-1 and from 3.08 to 7.73 kg·(hm2·a)-1, respectively. The annual N2O emissions from the CFS treatment were 43%-94% higher than those from the CF treatment, whereas the difference was only significant between 2015 and 2016 (p < 0.05). The GWP of the CF and CFS treatments varied from -1.95 to -0.28 t·(hm2·a)-1 and from -2.59 to -0.35 t·(hm2·a)-1, respectively. The cumulative GWP of the CFS treatment was 42% lower than that of the CF treatment between 2014 and 2017. In summary, the studied winter wheat-summer maize field acted as a sink of greenhouse gases under the conventional fertilization regime. The stalk incorporation further favored greenhouse gas mitigation despite the trade-offs between SOC sequestration and N2O emissions.
Key words: soil carbon sequestration      nitrous oxide (N2O)      global warming potential (GWP)      stalks incorporation      long-term fertilization      Guanzhong Plain     

自工业化以来, 大气二氧化碳(CO2)、甲烷(CH4)和氧化亚氮(N2O)浓度已分别增加48%、160%和23%, 是气候变暖的主要原因[1].农业是重要的温室气体排放源, 其CO2、CH4和N2O排放量分别占人为源总量的24%、12%和60%[1].农业温室气体减排潜力巨大, 通过优化田间管理, 农田每年可减少CO2当量(CO2-eq)排放量5.5×1012 t, 约占全球温室气体年排放总量的20%[2].

秸秆还田是保护性耕作措施之一, 不仅能够有效利用农业废弃物、提高土壤肥力和防治土壤退化[3], 还能提高土壤有机碳含量[ω(SOC)], 将CO2以有机碳的形式固持在土壤圈中[4].但全面评价秸秆还田对减缓气候变暖的作用, 还必须考察CH4和N2O排放的变化.农田温室气体净排放通常以全球增温潜势(以CO2-eq计)表示, 即土壤固碳速率、CH4净交换量和N2O排放量的百年增温潜势总和[5~7].旱地土壤是CH4微弱的汇[8], 其全球增温潜势主要取决于土壤固碳速率和N2O年排放量[6, 9].

土壤释放的N2O主要由微生物硝化、反硝化过程产生, 首要影响因素为氮底物浓度和氧化还原状态[10].当土壤气体扩散较慢时, N2O可被反硝化微生物还原成氮气[10].反硝化微生物大多为异养微生物, 以土壤有机碳作为能量来源和电子供体, 因此有机碳含量也是N2O排放的限制因素之一[10].关于秸秆还田对N2O排放的影响已有较多研究, 有研究认为秸秆带入的碳底物促进了反硝化作用和N2O排放[11~13], 有研究表明秸秆还田减少N2O排放, 主要是由于秸秆碳氮比高, 分解时需要从土壤中吸附氮底物[14, 15], 也有研究显示秸秆有利于保持土壤温度和湿度, 还通过增加土壤孔隙改变氧化还原状态和气体扩散条件, 从而促进或减少N2O排放[16~18].由此可见, 秸秆可通过多种途径改变N2O的产生和消耗, 导致秸秆还田条件下N2O年排放量极高的时空异质性.针对特定农田N2O排放通量进行多年观测, 是准确评价秸秆还田减排作用的关键.

关中平原位于黄土高原南缘, 属半湿润易干旱气候区, 土壤退化风险高[19, 20].近年来, 秸秆还田在该地区逐渐受到重视.有研究显示秸秆还田提高了当地农田土壤有机碳含量和N2O年排放量[11, 19, 20], 表明土壤固碳速率和N2O年排放量存在消长关系.Gu等[9]的研究对关中平原多种施肥制度(包括秸秆还田)下农田全球增温潜势进行了评估, 但所使用的N2O年排放量为经验模型的估计值, 其结果存在较大不确定性, 仍需原位观测进行验证.本研究基于长期定位施肥试验, 综合分析2010~2017年土壤有机碳含量和2014~2017年N2O排放通量的观测结果, 定量评价秸秆还田对土壤固碳速率、N2O年排放量和全球增温潜势的影响, 以期为实现农业“碳中和”目标提供基础数据支持.

1 材料与方法 1.1 试验地概况

本研究试验田位于陕西省杨凌区(34°17′51″N, 108°00′48″E), 由国家黄土肥力与肥料效益监测基地管理.该地区位于关中平原中部, 属暖温带大陆性季风气候, 多年平均气温13℃, 年均降水量550 mm, 降水主要集中在7~10月.

长期定位施肥试验开始于1990年[11, 20], 作物为冬小麦-夏玉米一年两熟.土壤为土垫旱耕人为土, 黄土母质, 质地为黏壤土.长期试验开始时表层土壤(0~20 cm)ω(有机碳)、容重和pH分别为7.3 g ·kg-1、1.32 g ·cm-3和8.6.

1.2 试验设计

田间试验采用随机区组试验设计, 每个小区面积为14 m×14 m.选取传统施肥(CF)和传统施肥加玉米秸秆(CFS)处理及不施肥对照(CK).CF处理施化肥(尿素、过磷酸钙和硫酸钾), 冬小麦季(每年10月至次年6月)氮、磷和钾素施入量分别为165、57和68 kg ·hm-2, 播种前均匀撒施; 夏玉米季(每年6月至10月)氮、磷和钾素施入量分别为188、25和78 kg ·hm-2, 播种后3~4周穴施(约15 cm); 两季作物收获时均保留地上约10 cm高根茬, 播种时翻入土中(约20 cm).CFS处理的化肥施入量、施肥方式与CF处理相同, 在冬小麦收获时同样保留地上约10 cm高根茬, 而在玉米收获时保留全部玉米秸秆[ω(C)和ω(N)分别为42.4%和0.9%, 干重约为4.4 t ·hm-2], 破碎成约3 cm长小段平铺施入, 播种时翻入土中(约20 cm).CK全年不施化肥, 秸秆管理措施与CF处理一致.试验田全年适时灌溉2~3次, 灌溉水为地下水, 方式为漫灌.

本研究土壤有机碳含量的观测期为2010~2017年(CK部分年份缺失), 此期间冬小麦和夏玉米品种分别为小偃22和郑丹958.N2O排放通量观测期为2015年10月~2017年10月, 主要农田管理措施列于表 1.同一试验田2014年10月~2015年10月N2O排放通量数据见文献[11].

表 1 2015~2017年主要农田管理措施 Table 1 Field management practices between 2015 and 2017

1.3 样品收集与分析

每年10月休闲期间, 在每个小区内使用土钻(内径2.5 cm)随机采集表层土壤(0~20 cm)30钻, 每10钻充分混合为一个样品.仔细挑出土样中秸秆、细根等残余物, 过2 mm筛, 使用重铬酸钾氧化法测定有机碳含量[9, 20].

N2O排放通量观测采用静态箱法[9, 11].采样装置由基座(50 cm ×50 cm ×15 cm)、中段箱(50 cm ×50 cm× 50 cm)和顶箱(50 cm ×50 cm ×50 cm)组成, 均为不锈钢材质.基座顶部有水槽, 中段箱和顶箱包覆绝热反光膜.顶箱上部预留1个直径为2 cm的圆孔以平衡气压, 侧面打孔设置采样管, 采样管连接三通接口.中段箱仅在植株高度超过50 cm后使用.每个处理随机设置3个基座, 埋入土壤约15 cm.采样时向基座水槽注水并将顶箱小心放入, 使用注射器采集箱内气体(约40 mL)4次, 每次间隔时间为8~10 min.采样结束10 h内使用气相色谱仪(7890B, Agilent Technologies Inc., USA)检测N2O浓度.采样均在08:00~11:00进行.通常情况下, 每周采样1次; 施肥或灌溉后每天采样1次, 直至N2O排放通量降至背景值.

气体采集的同时测量土壤温度和湿度.使用数字温度计(JM 624, 天津今明仪器有限公司)测量土壤温度(0~10 cm), 每个处理随机测量9次.使用便携式土壤水分速测仪(ML3, Delta-T Devices, UK)测量土壤含水量体积分数(0~5 cm), 每个处理随机测量9次.

1.4 数据处理

土壤充水孔隙度(WFPS)的计算公式为[11]:

(1)

式中, VSWC和BD分别为土壤含水量体积分数(cm3 ·cm-3)和容重(g ·cm-3); 常数2.65为土壤密度(g ·cm-3).

土壤固碳速率[SOCSR, t ·(hm2 ·a)-1]的计算公式为[5, 6, 9]:

(2)

式中, d[ω(SOC)]/dt为土壤有机碳含量的年变化率[g ·(kg ·a)-1]; 常数20和10分别为土层深度(cm)和单位换算系数.

N2O排放通量[F, g ·(hm2 ·d)-1]的计算公式为[6, 11]:

(3)

式中, ρ为标准状态下N2O密度, 取值1.25 g ·L-1; H为箱体高度(m), 常数10为单位换算系数; T为密闭期间箱内温度(℃); dc/dt为箱内N2O浓度随时间的变化率[μL ·(L ·d)-1].

N2O年排放量[N2Oan, kg ·(hm2 ·a)-1]的计算公式为[6, 11]:

(4)

式中, FiFi+1分别为第i和第i+1次采样时N2O排放通量; DiDi+1分别为第i和第i+1次采样时间(d).

N2O排放系数(EF)的计算公式为[11]:

(5)

式中, N2Oan-FN和N2Oan-CK分别为施肥处理和不施肥对照N2O年排放量; FN为年施氮量[kg ·(hm2 ·a)-1].

全球增温潜势[GWP, t ·(hm2 ·a)-1]的计算公式[9]为:

(6)

式中, 常数265为N2O百年增温潜势[2], 其余常数为单位换算系数.全球增温潜势为正值或负值分别表示农田是温室气体的源或汇.

1.5 统计分析

统计分析使用R软件包(版本号3.6.3, https://www.r-project.org/)完成.分别使用线性(lm函数)和逻辑斯蒂克(nls函数)回归模型拟合施肥处理土壤有机碳含量随施肥年限的变化规律, 在默认输出的统计分析结果中, 以赤池信息准则(AIC)选择最优拟合模型, 检验土壤有机碳含量的年变化率(即线性拟合模型的斜率)在处理间的差异(p < 0.05).采用单因素方差分析检验N2O年排放量在处理间和年际间的差异(p < 0.05).

利用边界线法分析N2O排放通量对土壤温度和充水孔隙度的响应[21].首先, 将土壤温度或充水孔隙度观测值设置为自变量(X), 将N2O排放通量设置为应变量(Y), 观测期内共有444组观测值.其次, 将X按递增顺序划分为8个区间, 每个区间内有55~56组观测值, 计算第i个区间内X的平均值(Xi)和Y的99%百分位数(Yi), 并赋值给边界点(Xi, Yi), 共有8个边界点.最后, 使用线性(lm函数)或高斯(nls函数)回归模型拟合边界线(p < 0.05).边界线表示当其他条件充分满足时, N2O排放通量对土壤温度或充水孔隙度的最大响应, 在边界线以下则表明N2O排放通量受到多种因素的共同影响.

2 结果与分析 2.1 土壤固碳速率

CK土壤有机碳含量在2010~2017年期间无明显变化趋势(图 1), 变异系数仅为3%, 表明不施肥土壤固碳速率为0; CF和CFS处理土壤有机碳含量均呈逐年增加趋势, 表明两种施肥制度下土壤均有固碳作用.线性和逻辑斯蒂克回归模型均能较好地拟合两施肥处理有机碳含量随施肥年限的变化规律(p < 0.001), 其中线性模型的赤池信息准则分别为-21和-15(图 1), 略低于逻辑斯蒂克模型的赤池信息准则(分别为-20和-14).因此, 选择线性模型的拟合结果计算土壤固碳速率.CF和CFS处理土壤有机碳含量的年变化率分别为0.23 g ·(kg ·a)-1和0.34 g ·(kg ·a)-1(图 1), 均显著大于0(p < 0.01)但两者之间差异不显著(P>0.05).两施肥处理土壤固碳速率分别为0.69 t ·(hm2 ·a)-1和0.97 t ·(hm2 ·a)-1, 后者较前者偏高40%, 表明秸秆还田提高了土壤固碳速率.

实线为线性回归模型, 虚线为平均值 图 1 2010~2017年ω(SOC)的变化 Fig. 1 Dynamics of ω(SOC) in the treatments between 2010 and 2017

2.2 土壤温度、充水孔隙度和N2O排放通量

CK、CF和CFS土壤温度和充水孔隙度在观测期内均具有相似的变化特征[图 2(a)2(b)].土壤温度在0.3~32.5℃之间, 2015~2016年和2016~2017年平均值分别为16.6℃和14.9℃.充水孔隙度在10% ~97%之间, 两个试验年平均值分别为49%和53%.

土壤温度和充水孔隙度的数据点为平均值(n=9); N2O排放通量的数据点和误差线分别为平均值和标准偏差(n=3); 箭头为施肥日期 图 2 2015~2017年土壤温度、充水孔隙度和N2O排放通量的变化特征 Fig. 2 Variations in soil temperature, water-filled pore space, and N2O fluxes in the treatments between 2015 and 2017

CK的N2O排放通量低[图 2(c)], 观测期内平均值仅为0.7 g ·(hm2 ·d)-1.CF和CFS处理N2O排放通量的季节变化趋于一致, 两个试验年中共出现7个较明显的排放高峰[峰值>22 g ·(hm2 ·d)-1], 大多出现在施肥和播种后[图 2(c)].排放高峰通常持续1~3周, 其余时间内排放通量与CK保持一致.施肥处理N2O排放通量有明显的年际差异, 两个试验年中CF和CFS处理最高峰值分别出现在夏玉米季[分别为190 g ·(hm2 ·d)-1和638 g ·(hm2 ·d)-1]和冬小麦季[分别为359 g ·(hm2 ·d)-1和598 g ·(hm2 ·d)-1].

边界线分析显示N2O排放通量与土壤温度和充水孔隙度的关系分别遵循高斯曲线和线性变化规律(图 3), 相关性均达到极显著水平(p < 0.01).表明当其他条件充分满足时, N2O排放通量在土壤温度约为19℃时达到最大值, 或随着充水孔隙度的增加而提高.

图 3 N2O排放通量与土壤温度和充水孔隙度的边界线分析 Fig. 3 Boundary lines of N2O fluxes against soil temperature and water-filled pore space

2.3 N2O年排放量和排放系数

3个试验年中, CF和CFS处理N2O年排放量均显著大于CK(p < 0.001, 表 2).CFS处理N2O年排放量较CF处理偏高43% ~94%, 但仅在2015~2016试验年有显著性差异(p < 0.05).CFS处理3 a累计N2O排放量较CF处理偏高48%(p < 0.05), 表明秸秆还田促进了N2O排放.两施肥处理N2O年排放量均有显著年际间差异, 其中2016~2017试验年较其他年份偏高148% ~286%(p < 0.05).两施肥处理排放系数变化范围为0.40% ~1.90%, 平均值为1.06%.

表 2 2014~2017年N2O年排放总量1)、排放系数和全球增温潜势 Table 2 Annual N2O emissions, emission factors, and global warming potentials in the treatments between 2014 and 2017

2.4 全球增温潜势

2014~2017年, CK全球增温潜势变化范围为0.04~0.11 t ·(hm2 ·a)-1, 表明不施肥农田为温室气体微弱的源(表 2); CF和CFS处理全球增温潜势分别在-1.95~-0.28 t ·(hm2 ·a)-1和-2.59~-0.37 t ·(hm2 ·a)-1之间, 表明两处理农田均为温室气体的汇. 3个试验年中, CFS处理全球增温潜势较CF处理偏低23% ~56%; CFS处理3 a累计全球增温潜势较CF处理偏低42%, 表明秸秆还田更有利于温室气体减排.

3 讨论 3.1 旱地农田土壤固碳

土壤碳库变化是碳输入(如粪肥、秸秆、根系及其分泌物等)和输出(如土壤异养呼吸和风力、水力侵蚀等)的动态平衡过程[4].Li等[19]在泾阳(同位于关中平原)开展的试验结果显示, 全化肥(无秸秆还田)处理土壤固碳速率为负值[-0.26 t ·(hm2 ·a)-1], 表明在缺少额外碳输入的情况下由耕作导致土壤有机碳损失.本研究中CK土壤有机碳含量在过去8 a中变化小(图 1), 表明在不施肥但保留部分根茬的管理方式下, 土壤碳库已能达到收支平衡; 而CF和CFS处理土壤固碳速率为正值[0.69~0.97 t ·(hm2 ·a)-1], 表明碳输入已超出维持土壤碳库平衡所需量.其中CFS处理土壤固碳速率较泾阳[19]秸秆还田处理[0.55 t ·(hm2 ·a)-1]偏高约75%, 其原因可能有: 第一, 土壤固碳速率与秸秆还田量呈正相关关系[22].本研究CFS处理为全量玉米秸秆还田, 秸秆碳输入量高, 但泾阳试验未提供相关数据[15], 无法进行比较.第二, 秸秆还田条件下土壤固碳速率与初始有机碳含量呈负相关关系[22].本研究中土壤初始有机碳含量低, 因此固碳速率较高.

本研究CFS处理土壤固碳速率较CF处理偏高0.28 t ·(hm2 ·a)-1, 占秸秆输入碳的15.1%, 略低于全国平均值16.1%[4], 但高于东北黑钙土(12.9%)[23]、南方红土(8.1% ~10.6%)[24, 25]和华北潮土(8.5% ~9.1%)[26, 27], 这主要是由于土壤质地、气候因素、秸秆性质和施用方式等诸多因素造成土壤固碳效率的差异.

部分长期定位施肥试验结果显示, 土壤有机碳含量随施肥年限的变化规律符合逻辑斯蒂克模型[5, 9, 28], 表明土壤有机碳含量存在上限, 且同一农田在不同施肥制度下土壤有机碳含量上限也存在差异[5, 9].同时也表明土壤固碳速率可能是时间的函数, 随施肥年限的延长而逐渐降低[9, 28].本研究中两施肥处理土壤有机碳含量均随施肥年限呈线性增长(图 1), 可能是由于初始土壤有机碳含量低、施肥年限短, 尚无法评估土壤有机碳含量上限[9, 28].本研究采用线性模型拟合土壤有机碳含量变化规律仅基于统计检验结果, 其土壤固碳速率估计值仅在2010~2017年是可靠的.

3.2 N2O排放和排放系数

本研究中两施肥处理N2O排放高峰通常与施肥和播种有关(图 2), 与大多数的田间观测结果一致[11~15].施肥处理N2O年排放量在1.6~7.7 kg ·(hm2 ·a)-1之间(表 2), 基本位于我国旱地农田N2O年排放变化范围内[29].

本研究中N2O排放通量与土壤温度和充水孔隙度均无显著相关性(图 3), 主要是由于N2O排放受到多种因素(如温度、湿度、氮底物和碳底物等)的共同影响[11~15].边界线为确定田间N2O排放通量与某单一因素的相关性提供了很好的分析方法[21, 30].本研究中边界线分析结果显示N2O排放通量在土壤温度为19℃时达到最大值, 该最适土壤温度低于我国亚热带[21]及欧洲温带海洋气候[30]农田的研究结果(23~25℃), 可能与土壤微生物对气候的长期适应有关.本研究中N2O排放通量随土壤充水孔隙度的增加而提高, 表明N2O排放高峰主要由反硝化过程产生.已有研究则显示N2O排放通量的最适充水孔隙度为69% ~72%[21, 30], 当充水孔隙度高于此阈值时反硝化过程的主要产物为氮气, 因此N2O排放通量有所降低.而本研究中充水孔隙度高于80%的观测值较少, 尚无法确定N2O排放通量在高值区域的变化趋势.

本研究中秸秆还田促进了N2O排放(表 2).首先, CFS处理冬小麦季施氮量高, 在尿素与CF处理相同的基础上引入了秸秆, 提供了更多的氮底物[约40 kg ·(hm2 ·a)-1].其次, 秸秆引入了易分解有机质, 促进了反硝化微生物活动.这两点是冬小麦季CFS处理N2O排放高于CF处理的主要原因.再次, 长期秸秆还田提高了土壤有机碳含量(图 1), 有利于反硝化作用和N2O排放, 这是夏玉米季CFS处理N2O排放高的主要原因.

本研究中施肥处理N2O排放量有显著的年际间差异, 且均表现为2016~2017年最高(表 2), 这除了与该试验年土壤充水孔隙度较高(图 2)有关以外, 还可能与施肥和灌溉等管理措施的配合有关.如2015年[11]和2016年夏玉米季均在施肥一周后进行灌溉, 较高的N底物浓度和土壤湿度有利于N2O产生和排放, CF和CFS处理最高峰值分别达到190.1~280.0 g ·(hm2 ·d)-1和382.1~643.0 g ·(hm2 ·d)-1; 而2017年夏玉米季在施肥前进行灌溉, 两处理N2O排放最高峰值分别为86.1 g ·(hm2 ·d)-1和252.7 g ·(hm2 ·d)-1.尽管2017年夏玉米季两处理N2O排放量均低于2015年和2016年, 但由于2016~2017年冬小麦季N2O排放量高, 在N2O年排放量上表现为2016~2017年显著偏高(表 2).由于N2O排放存在较大的年际间差异, 准确评估秸秆还田对N2O排放的影响还需要连续多年的观测.

3个试验年中CFS处理排放系数均高于CF处理(表 2), 表明秸秆还田提高了排放系数.两处理排放系数平均值分别为0.90%和1.22%, 分别低于和高于IPCC缺省值1%[31].因此, 为准确估算区域农田N2O排放总量, 对于传统施肥农田和秸秆还田农田应采用不同的排放系数.由于排放系数具有较大的年际间差异(表 2), 准确评估排放系数也需要连续多年的观测.

3.3 秸秆还田对温室气体净排放的影响

本研究中不施肥农田为温室气体微弱的源(表 2), 主要是由于土壤不具备固碳作用(图 1)且释放N2O(表 2).在实际生产中, 关中平原几乎没有不施肥农田, 该排放源可忽略不计.在传统施肥管理方式下, 所研究的农田已是温室气体的汇, 主要是通过保留部分根茬的方式增加了碳输入.虽然全量玉米秸秆还田在提高土壤固碳速率(图 1)的同时也提高了N2O年排放量(表 2), 即土壤固碳速率和N2O年排放量存在消长关系, 但整体更有利于温室气体的减排(表 2).

有两点仍需要注意.第一, 土壤固碳速率不能长期维持较高水平.随着土壤有机碳含量接近其上限, 土壤固碳速率逐渐降低, 农田可能转变为温室气体的源[9].第二, N2O排放年际间差异大.在N2O排放量高的年份(如2016~2017年), 农田的汇功能有所减弱(表 2)甚至可能转变为源.因此, 采取适当措施提高土壤有机碳含量上限和降低N2O排放量是非常必要的.

4 结论

(1) 2010~2017年, 秸秆还田处理土壤固碳速率较传统施肥处理偏高约40%, 秸秆还田提高了土壤固碳速率.

(2) 2014~2017年, 秸秆还田处理N2O累计排放量较传统施肥处理偏高约48%, 秸秆还田促进了N2O排放.

(3) 在传统施肥管理方式下, 所研究农田已是温室气体的汇.尽管秸秆还田条件下土壤固碳速率和N2O排放量存在消长关系, 但仍然更有利于温室气体减排.

参考文献
[1] Stocker T F, Qin D, Plattner G K, et al. Climate Change 2013:the physical science basis. Contribution of working group I to the fifth assessment report of the intergovernmental panel on climate change[M]. Cambridge: Cambridge University Press, 2013.
[2] Smith P, Martino D, Cai Z C, et al. Greenhouse gas mitigation in agriculture[J]. Philosophical Transactions of the Royal Society B: Biological Sciences, 2008, 363(1492): 789-813. DOI:10.1098/rstb.2007.2184
[3] Diacono M, Montemurro F. Long-term effects of organic amendments on soil fertility. A review[J]. Agronomy for Sustainable Development, 2010, 30(2): 401-422. DOI:10.1051/agro/2009040
[4] Zhao Y C, Wang M Y, Hu S J, et al. Economics- and policy-driven organic carbon input enhancement dominates soil organic carbon accumulation in Chinese croplands[J]. Proceedings of the National Academy of Sciences of the United States of America, 2018, 115(16): 4045-4050. DOI:10.1073/pnas.1700292114
[5] Shang Q Y, Yang X X, Gao C M, et al. Net annual global warming potential and greenhouse gas intensity in Chinese double rice-cropping systems: a 3-year field measurement in long-term fertilizer experiments[J]. Global Change Biology, 2011, 17(6): 2196-2210. DOI:10.1111/j.1365-2486.2010.02374.x
[6] Cao Y S, Shan Y, Wu P C, et al. Mitigating the global warming potential of rice paddy fields by straw and straw-derived biochar amendments[J]. Geoderma, 2021, 396. DOI:10.1016/j.geoderma.2021.115081
[7] 蓝兴福, 王晓彤, 周雅心, 等. 炉渣与生物炭施加对福州平原水稻田温室气体排放的后续影响[J]. 环境科学, 2020, 41(1): 489-498.
Lan X F, Wang X T, Zhou Y X, et al. Subsequent effects of slag and biochar application on greenhouse gas emissions from paddy fields in the Fuzhou Plain[J]. Environmental Science, 2020, 41(1): 489-498.
[8] Cowan N, Maire J, Krol D, et al. Agricultural soils: a sink or source of methane across the British Isles?[J]. European Journal of Soil Science, 2021, 72(4): 1842-1862. DOI:10.1111/ejss.13075
[9] Gu J X, Yuan M X, Liu J X, et al. Trade-off between soil organic carbon sequestration and nitrous oxide emissions from winter wheat-summer maize rotations: implications of a 25-year fertilization experiment in Northwestern China[J]. Science of the Total Environment, 2017, 595: 371-379. DOI:10.1016/j.scitotenv.2017.03.280
[10] 朱永官, 王晓辉, 杨小茹, 等. 农田土壤N2O产生的关键微生物过程及减排措施[J]. 环境科学, 2014, 35(2): 792-800.
Zhu Y G, Wang X H, Yang X R, et al. Key microbial processes in nitrous oxide emissions of agricultural soil and mitigation strategies[J]. Environmental Science, 2014, 35(2): 792-800.
[11] 郝耀旭, 刘继璇, 袁梦轩, 等. 长期定位有机物料还田对关中平原冬小麦-玉米轮作土壤N2O排放的影响[J]. 环境科学, 2017, 38(6): 2586-2593.
Hao Y X, Liu J X, Yuan M X, et al. Effects of long-term organic amendments on soil N2O emissions from winter wheat-maize cropping systems in the Guanzhong Plain[J]. Environmental Science, 2017, 38(6): 2586-2593.
[12] Abalos D, Sanz-Cobena A, Garcia-Torres L, et al. Role of maize stover incorporation on nitrogen oxide emissions in a non-irrigated Mediterranean barley field[J]. Plant and Soil, 2013, 364(1-2): 357-371. DOI:10.1007/s11104-012-1367-4
[13] 许宏伟, 李娜, 冯永忠, 等. 氮肥和秸秆还田方式对麦玉轮作土壤N2O排放的影响[J]. 环境科学, 2020, 41(12): 5668-5676.
Xu H W, Li N, Feng Y Z, et al. Effects of nitrogen fertilizer and straw returning methods on N2O emissions in wheat-maize rotational soils[J]. Environmental Science, 2020, 41(12): 5668-5676.
[14] Xu C, Han X, Ru S H, et al. Crop straw incorporation interacts with N fertilizer on N2O emissions in an intensively cropped farmland[J]. Geoderma, 2019, 341: 129-137. DOI:10.1016/j.geoderma.2019.01.014
[15] Shelp M L, Beauchamp E G, Thurtell G W. Nitrous oxide emissions from soil amended with glucose, alfalfa, or corn residues[J]. Communications in Soil Science and Plant Analysis, 2000, 31(7-8): 877-892. DOI:10.1080/00103620009370484
[16] Lehtinen T, Schlatter N, Baumgarten A, et al. Effect of crop residue incorporation on soil organic carbon and greenhouse gas emissions in European agricultural soils[J]. Soil Use and Management, 2014, 30(4): 524-538. DOI:10.1111/sum.12151
[17] Frimpong K A, Baggs E M. Do combined applications of crop residues and inorganic fertilizer lower emission of N2O from soil?[J]. Soil Use and Management, 2010, 26(4): 412-424. DOI:10.1111/j.1475-2743.2010.00293.x
[18] Wu D, Wei Z J, Well R, et al. Straw amendment with nitrate-N decreased N2O/(N2O+N2) ratio but increased soil N2O emission: a case study of direct soil-born N2 measurements[J]. Soil Biology and Biochemistry, 2018, 127: 301-304. DOI:10.1016/j.soilbio.2018.10.002
[19] Li S, Hu M J, Shi J L, et al. Integrated wheat-maize straw and tillage management strategies influence economic profit and carbon footprint in the Guanzhong Plain of China[J]. Science of the Total Environment, 2021, 767. DOI:10.1016/j.scitotenv.2021.145347
[20] Yang X Y, Ren W D, Sun B H, et al. Effects of contrasting soil management regimes on total and labile soil organic carbon fractions in a loess soil in China[J]. Geoderma, 2012, 177-178: 49-56. DOI:10.1016/j.geoderma.2012.01.033
[21] Gu J X, Guo H J, Xiang H Y. The utility of a boundary line approach for simulating denitrification and nitrous oxide emissions from a Regosol under summer maize-winter wheat crop rotation in Southwest China[J]. Geoderma Regional, 2020, 20. DOI:10.1016/j.geodrs.2020.e00252
[22] Berhane M, Xu M, Liang Z Y, et al. Effects of long-term straw return on soil organic carbon storage and sequestration rate in North China upland crops: a meta-analysis[J]. Global Change Biology, 2020, 26(4): 2686-2701. DOI:10.1111/gcb.15018
[23] 高洪军, 彭畅, 张秀芝, 等. 秸秆还田量对黑土区土壤及团聚体有机碳变化特征和固碳效率的影响[J]. 中国农业科学, 2020, 53(22): 4613-4622.
Gao H J, Peng C, Zhang X Z, et al. Effects of corn straw returning amounts on carbon sequestration efficiency and organic carbon change of soil and aggregate in the black soil area[J]. Scientia Agricultura Sinica, 2020, 53(22): 4613-4622. DOI:10.3864/j.issn.0578-1752.2020.22.008
[24] 张敬业, 张文菊, 徐明岗, 等. 长期施肥下红壤有机碳及其颗粒组分对不同施肥模式的响应[J]. 植物营养与肥料学报, 2012, 18(4): 869-876.
Zhang J Y, Zhang W J, Xu M G, et al. Response of soil organic carbon and its particle-size fractions to different long-term fertilizations in red soil of China[J]. Plant Nutrition and Fertilizer Science, 2012, 18(4): 868-875.
[25] 蔡岸冬, 张文菊, 申小冉, 等. 长期施肥土壤不同粒径颗粒的固碳效率[J]. 植物营养与肥料学报, 2015, 21(6): 1431-1438.
Cai A D, Zhang W J, Shen X R, et al. Soil carbon sequestration efficiency of different particle-size fractions after long-term fertilization[J]. Journal of Plant Nutrition and Fertilizer, 2015, 21(6): 1431-1438.
[26] 高伟, 杨军, 任顺荣. 长期不同施肥模式下华北旱作潮土有机碳的平衡特征[J]. 植物营养与肥料学报, 2015, 21(6): 1465-1472.
Gao W, Yang J, Ren S R. Balance characteristics of soil organic carbon under different long-term fertilization models in the upland fluvo-aquic soil of North China[J]. Journal of Plant Nutrition and Fertilizer, 2015, 21(6): 1465-1472.
[27] 黄少辉, 杨军芳, 杨云马, 等. 长期不同施肥措施下华北潮土土壤有机碳的固存变化[J]. 华北农学报, 2019, 34(S1): 168-175.
Huang S H, Yang J F, Yang Y M, et al. Soil organic carbon sequestration under long term fertilization practices in Fluvo-Aquic soil of North China Plain[J]. Acta Agriculturae Boreali-Sinica, 2019, 34(S1): 168-175. DOI:10.7668/hbnxb.20190635
[28] Stewart C E, Paustian K, Conant R T, et al. Soil carbon saturation: concept, evidence and evaluation[J]. Biogeochemistry, 2007, 86(1): 19-31. DOI:10.1007/s10533-007-9140-0
[29] Lv F L, Yang X Y, Xu H H, et al. Effect of organic amendments on yield-scaled N2O emissions from winter wheat-summer maize cropping systems in Northwest China[J]. Environmental Science and Pollution Research, 2020, 27(25): 31933-31945. DOI:10.1007/s11356-020-09491-9
[30] Schmidt U, Thöni H, Kaupenjohann M. Using a boundary line approach to analyze N2O flux data from agricultural soils[J]. Nutrient Cycling in Agroecosystems, 2000, 57(2): 119-129. DOI:10.1023/A:1009854220769
[31] IPCC. 2006 IPCC guidelines for national greenhouse gas inventories[R]. Japan: IPCC, 2006.