烟-稻轮作不同施肥土壤N2O排放对水分的响应

赵伟东，郭宝玲，郑祥洲*，汤水荣，孟磊，张玉树

（1.海南大学热带作物学院，海口 570228；2.福建省农业科学院土壤肥料研究所/福建省植物营养与肥料重点实验室，福州 350013）

近年来，随着温室气体急剧增加，极端气候等一系列重大全球性的生态环境问题愈发严重，引起人们广泛关注[1-2]。增温潜势最高的N2O 不仅能引起温室效应，还会对臭氧层造成严重破坏，甚至威胁地球生物健康[3-4]。研究表明，农田土壤是N2O 的重要排放源，每年农田的N2O 排放量占全球N2O 总排放量的43%以上[5-6]。同时，农田生态系统也具有巨大的固碳减排潜力[7]。因此，降低农田土壤N2O 排放对缓解气候变化、提高氮肥利用率等具有极其重要的意义。施肥作为农业生产中常见的农田管理措施之一，在提升土壤肥力的同时也影响着农田N2O 排放，而在农田生态系统中不同施肥种类和不同施肥量对N2O 排放的影响存在较大差异[8-11]。普遍认为施入氮肥会促进土壤N2O 的排放[8，12-13]。通过合理的施肥管理不仅能够培肥土壤，而且对缓解土壤N2O 排放也具有积极意义[11，14]。施肥是补充土壤无机氮和影响N2O排放的最直接的驱动因子。农田土壤常年高氮投入导致土壤氮过剩的问题已成为人们关注的焦点。研究发现，氮肥的过量施用显著提高nirK 和nosZ 基因拷贝数[15-16]，引起N2O 排放量增加[17]。然而，Yin 等[18]的研究却表明高氮投入降低nirS 基因拷贝数。秸秆还田具有提高土壤碳储量，增加土壤肥力[19]，降低土壤N2O 排放的作用[14，20]。但裴淑玮等[21]的研究结果表明秸秆还田显著增加了农田N2O 排放。由此可见，高氮投入和秸秆还田对土壤N2O 排放的影响均未得到共识。土壤水分是影响土壤N2O 排放的重要环境因素之一[22]。水分可以改变土壤微生物生存环境，进而影响N2O 排放。有研究报道，在低水分条件下，土壤N2O 排放主要由硝化作用主导，而在高水分条件下则是由反硝化作用主导[23-26]。徐华等[27]认为水分通过改变土壤中O2浓度、温度和底物浓度等条件，能在很大程度上改变自养型细菌的组成，进而影响土壤N2O 排放。而Qin 等[28]研究发现，增加土壤水分含量主要影响含nirK 和nosZ 基因微生物的群落结构，从而调控N2O的排放。水分含量调控土壤微生物群落结构，进而决定硝化作用和反硝化作用的相对强度及N2O 的排放。因此，不同水分条件下土壤N2O 排放的微生物机制尚未得到统一结论，还需结合土壤微生物群落结构等进一步探索。

南方地区是我国主要的水稻生产基地，高温多雨的气候条件造就了普遍的水旱轮作种植模式。烟-稻轮作是我国福建省水旱轮作的典型模式之一，为保证两季作物的产量，通常投入大量化学肥料，从而容易引起硝态氮淋溶和N2O 排放加剧等一系列生态环境问题。通常情况下，氮过量会降低土壤pH，并且增加N2O 排放，而秸秆添加能缓解土壤酸化，增加土壤有机质含量，降低N2O 排放。轮作过程中的干湿交替导致两季土壤水分存在很大差异，这不仅会改变土壤微生物群落结构，而且对N2O 排放也具有很大影响。不同水分条件下，N2O 排放对不同肥料投入的响应还需进一步验证。基于此，本试验以水-旱轮作长期定位试验土壤为研究对象，设置60%持水量（WHC）和淹水两个水分条件，用来模拟烟叶季和水稻季的田间水分状况，研究不同水分对长期定位不同施肥量土壤N2O 排放的影响，以期为亚热带地区稻田温室气体减排提供理论依据和技术支撑。

1 材料与方法

1.1 试验地概况

田间原位监测试验地建于2008 年，位于中国福建省三明市将乐县古镛镇张公村（26°44'53″N，117°26'48″E）。该地属于亚热带季风气候，年平均气温为18.9 ℃，年平均降水量为1 667~1 880 mm，年均日照约1 736 h；土壤母质为花岗岩，土壤腐殖质层15~30 cm。监测点处于山腰位置，坡度为8°。试验站点设置时耕层土壤平均pH 为5.31，有效磷20.5 mg·kg-1、有机质（SOM）25.85 g·kg-1、全氮（TN）1.39 g·kg-1，铵态氮（和硝态氮含量分别为7.12 mg·kg-1和2.20 mg·kg-1。每年2月或3月种植烟草，7月种植水稻。样地设置4 个处理，分别为不施肥（CK），推荐施肥（NPK），高氮施肥（NhPK），推荐施肥+水稻秸秆还田（NPKS），具体施肥量见表1。采用完全随机区组设计，每个施肥处理设置3 个重复地块（7 m×4 m），地块用砖框隔开。2021 年水稻季收获后采用S 型采样法多点采集表层土壤（0~15 cm），充分混匀后分为两部分：一部分土壤于室温下风干并过2 mm筛，用于土壤理化性质的测定；另一部分土样于4 ℃条件下低温储存，并在一周内进行培养试验。

表1 烟−稻轮作模式下不同施肥处理及施肥量Table 1 Different fertilization rates in fertilization treatments under tobacco-rice rotation modes

1.2 试验设计

在4 种长期不同施肥的基础上，设置60%WHC（U）和淹水（F）2 个不同的水分条件，共8 个处理，每个处理称取30 g 土壤（干基）于250 mL 三角瓶中，加入去离子水使水分含量达到50%WHC，25 ℃预培养1周，使土壤微生物充分激活。预培养结束后，将尿素溶液均匀施入到相应处理的土壤中，使得每个处理添加的氮含量达到200 mg·kg-1，同时调节含水量，使每个长期施肥处理土壤分别达到目标水分，每个处理设置3 个重复，用透气膜封口后于25 ℃条件下恒温培养。分别在培养后的第1、3、5、7、10、14、21、28、32 天取样测定N2O和无机氮浓度。

1.3 样品采集

取气时先将三角瓶连接带有软管和三通阀的橡胶塞，连接处已涂抹强力胶并晾干，保证装置气密性良好，利用配套的多孔抽真空装置，连续置换空气3次，每次不低于1 min，再次放入25 ℃条件下恒温培养6 h，然后用连接有三通阀的20 mL 针管进行取气，取气时缓慢推拉针管3 次，使瓶内气体混合均匀后再进行取气。取气完毕后，向三角瓶中加入150 mL 1 mol·L-1KCl 溶液，190 r·min-1振荡提取1 h 后过滤，得到的滤液用来测定浓度。培养结束后，均匀分取最后一天土壤样品约10.00 g，−80 ℃冷冻保存，于两周内进行荧光定量qPCR分析。

1.4 样品分析与测定

土壤pH值采用水土比为2.5∶1的电位法测定；有效磷含量用1 mol·L-1NH4F 溶液浸提、钼蓝比色法测定；SOM 采用重铬酸钾容量法-外加热法测定；TN 采用半微量开氏法测定，具体测定方法参考《土壤农业化学分析方法》。水溶性有机碳（DOC）含量采用有机碳氮分析仪（Shimadzu Corp，Kyoto，日本）测定；土壤最大持水量的测定参照文献[29]中的方法；含量采用1 mol·L-1KCl 溶液浸提后使用SKALAR 连续流动分析仪（Skalar，Breda，荷兰）测定；N2O 气体浓度采用安捷伦气相色谱仪（7890A，美国）测定。

采用MOBIO公司的土壤快速提取试剂盒（Power⁃SoilTMDNA Isolation Kits）进行DNA 提取，提取步骤参考试剂盒说明书，构建好的质粒经测序鉴定无误后用超微量紫外分光光度计（（NanoDrop2000，Thermo Fisher Scientific，美国）测定OD260的值。用1%的琼脂糖凝胶电泳检测提取DNA片段，DNA样品于-20 ℃条件下保存备用。采用土壤中氨单加氧酶（AMO）基因，亚硝酸还原酶nirK、nirS 基因和氧化亚氮还原酶nosZ基因拷贝数分别表示氨氧化古菌（AOA）、氨氧化细菌（AOB）与nirK、nirS 和nosZ 反硝化菌，PCR 扩增所用引物和反应条件如表2 所示。定量分析采用SYBR GREEN I 法，反应体系为20 μL，其中包括10 μL 2X ChamQ SYBR Color qPCR Master M，10 μmol·L-1正反向引物各0.8 μL，0.4 μL 50 X ROX Reference Dy，2 μL Template（DNA），6 μL 无菌水。并使用10-1~10-8浓度梯度的标准质粒作为模板进行荧光定量PCR 扩增（ABI 7300 型荧光定量PCR 仪，Applied Bio⁃systems，美国），最后根据扩增曲线计算基因丰度。

表2 荧光实时定量PCR扩增的引物和反应条件Table 2 Amplification primers and reaction conditions for real-time quantitative PCR

1.5 数据处理与分析

N2O排放通量的计算如公式（1）所示：

式中：F 为N2O 排放通量，μg·kg-1·h-1；ρ为标准状态下N2O-N的密度；dc/dt为单位时间内培养瓶内气体浓度增加量，10-6·h-1或10-9·h-1；V 为培养瓶中气体的有效空间体积，m3；m为培养瓶内的烘干土质量，kg；T为培养温度，℃。

气体累积排放量通过相邻2 次培养时间的平均气体排放通量与时间相乘后加权累积计算，如公式（2）所示：

式中：Ci+1为第i 次和第i+1 次采样期间的气体累积排放量，mg·kg-1；Fi为第i次采样时气体瞬时排放通量；i为采样次数；D为两次采样间隔时间，d。

数据为3 次重复的平均值，采用Excel 2019 进行数据分析，利用SPSS 27.0 软件在Duncan（SSR）方法下分析不同水肥处理各指标的差异显著性，利用Ori⁃gin Pro 2022制图。

2 结果与分析

2.1 长期不同施肥对土壤理化性质的影响

如表3 所示，与CK 相比，各长期施肥处理均显著降低土壤pH，但显著提高土壤有效磷和DOC 含量。各施肥处理下SOM 含量为NPK>NPKS>NhPK>CK（15.58~20.09 g·kg-1），相对于CK 处理，NPK 和NPKS处理均显著增加SOM 含量。土壤TN 含量为NPKS>NhPK>NPK>CK（1.05~1.26 g·kg-1），统 计 分 析 表 明NPKS 处理下TN 含量显著高于其他施肥处理。碳氮比（C/N）为NPK>NhPK>NPKS>CK（8.66~10.57），NPK显著高于其他处理。与CK 相比，长期不同施肥各处理均显著降低土壤含量，NhPK 和NPKS 均显著提高土壤含量（P<0.05）。

表3 长期不同田间施肥处理土壤理化性质Table 3 Soil properties of different treatments after long term experiment

2.2 土壤

含量 的动态变化

整个培养周期内，60%WHC 下NH+4-N 含量变化范围为13.17~196.83 mg·kg-1，各处理土壤NH+4-N 动态变化趋势基本一致（图1a），UNPK 和UNhPK 处理于培养第3 天左右达到峰值，UCK 和UNPKS 处理于培养第5 天达到峰值，5 d 后均转变为下降趋势。培养的7~28 d 时，UCK 处理NH+4-N 含量显著高于其他处理。淹水条件下NH+4-N 含量变化范围为0~138.85mg·kg-1，各处理土壤变化趋势高度趋同（图1b），除UNPKS 第3 天达到峰值外，其他处理皆于培养第5天达到峰值，5 d后转变为下降趋势。14~32 d，FCK 处理含量显著高于其他处理。氮肥添加会明显提高土壤含量，大量的由土壤硝化作用转化为。整个培养周期内，60%WHC 下含量变化范围为9.92~267.17 mg·kg-1，各处理土壤动态变化趋势大致相同（图1c），除UNPK处理于第21 天达到峰值外，其他处理含量均一直增加至培养结束。3~32 d，UCK处理含量显著低于其他处理。淹水条件下含量变化范围为4.65~162.60 mg·kg-1，淹水各处理土壤变化趋势高度相似（图1d），FCK和FNPK处理在第21天达到峰值，而FNhPK 和FNPKS 处理在第28 天达到峰值，第10天后整体变化幅度逐渐增大。10~32 d，FCK处理含量显著低于其他处理。总体而言，60%WHC 下各施肥处理和含量的变化幅度均显著大于对应的淹水各施肥处理。

图1 两种不同水分条件下各施肥处理和含量动态变化Figure 1 Dynamic changes ofandcontent per fertilization treatment under two different moisture conditions

2.3 土壤N2O排放通量和累积排放量

60%WHC 条件下，N2O 排放通量（以N 计）的变化范围为0.05~5.13 μg·kg-1·h-1（图2a）。UCK 处理在第7 天和第28 天分别出现两次峰值，且第7 天的N2O 排放通量更高；UNPK和UNhPK处理均在第28天出现峰值，而UNPKS 处理仅在第14 天显示出微弱的峰值，随后一直呈降低趋势，直至培养结束。淹水条件下，N2O 排放通量（以N 计）的变化范围为0.03~73.61 μg·kg-1·h-1；各处理N2O排放通量动态变化趋势基本一致（图2b）。除FNhPK 处理在第10 天达到峰值外，其他处理皆于培养第14 天达到峰值，随后一直处于下降状态，直至培养结束；3~21 d，FNPKS 处理略高于其他处理。总体而言，整个培养过程中，淹水条件下各施肥处理N2O 排放通量峰值均远高于60%WHC 对应的各施肥处理。

图2 两种不同水分条件下各施肥处理N2O排放动态变化Figure 2 Dynamic changes of N2O emissions per fertilization treatment under two different moisture conditions

60%WHC 条件下，至培养结束，各处理N2O 累积排放量达到0.56~1.44 mg·kg-1（图3a）。与UCK 相比，UNPK、UNhPK、UNPKS 处理的N2O 累积排放量分别降低了45%、38%、61%，均显著低于UCK 处理；UNPKS显著低于其他施肥处理。淹水条件下，至培养结束，各处理N2O 累 积排放量为14.89~20.70 mg·kg-1（图3b）。与FCK 相比，FNPK 和FNPKS 处理N2O 累积排放量分别提高了10%和18%，而FNhPK 降低了12%。总体而言，淹水各施肥处理N2O 累积排放量大于60%WHC 各施肥处理，淹水条件显著促进了各施肥处理的N2O 排放。土壤N2O 排放和无机氮相关性分析结果表明，在两种水分条件下，土壤N2O 累积排放量均与含量呈正相关（P<0.01）（图4）。

图3 两种不同水分条件下各施肥处理土壤N2O累积排放量Figure 3 Cumulative emission of N2O from per fertilization treatments under two different moisture conditions

图4 两种不同水分条件下土壤N2O累积排放量和含量的相关性Figure 4 Correlation between cumulative N2O emissions andcontent under two different moisture conditions

2.4 微生物功能基因丰度

60%WHC 培养结束后，不同施肥处理下土壤（干土）的AOA和AOB基因拷贝数分别为8.06×106~1.54×107copies·g-1和1.14×107~1.75×107copies·g-1（表4）。与UCK相比，UNPK、UNhPK、UNPKS处理的AOA基因拷贝数分别提高48%、42%、25%，AOB 基因拷贝数分别提高17%、1%、35%。相比于UCK，UNPK、UNhPK含量呈负相关，与处理显著提高AOA 基因拷贝数，UNPKS 处理显著提高AOB 基因拷贝数。淹水条件培养结束后，不同处理AOA 和AOB 基因拷贝数分别为1.45×107~3.59×107copies·g-1和1.35×107~3.00×107copies·g-1（表4）。与FCK 相比，FNPK、FNhPK、FNPKS 的AOA 基因拷贝数分别提高60%、13%、26%，AOB 基因拷贝数分别提高55%、41%、49%。相比于FCK 处理，FNPK、FNhPK、FNPKS 均显著提高了AOA 和AOB 基因拷贝数；与FNPK 相比，FNhPK、FNPKS 显著降低AOA 基因拷贝数，FNhPK 显著降低AOB 基因拷贝数。总体而言，淹水条件各施肥处理AOA 和AOB 基因拷贝数分别为60%WHC各处理的1.80~2.49倍和1.19~2.19倍。

表4 不同水分条件下各施肥处理微生物相关功能基因丰度Table 4 Abundance of microbial functional genes under different water conditions under different fertilization treatments

60%WHC 培养结束后，不同处理的nirK、nirS 和nosZ 基因拷贝数分别为8.78×106~1.23×107、1.79×107~3.42×107、5.78×107~7.90×107copies·g-1（表4）。 与UCK 相比，UNhPK、UNPKS 处理的nirK 基因拷贝数分别降低28%和10%，UNPK、UNhPK、UNPKS 的nirS 基因拷贝数分别降低19%、91%、71%；UNPK 和UNPKS处理的nosZ 基因拷贝数分别提高9% 和7%，而UNhPK 处理降低24%。各处理中nosZ 基因拷贝数均显著高于nirK和nirS基因拷贝数，分别提高0.84~0.87倍和0.52~0.74 倍。与UCK 相比，UNhPK 处理均显著降低nirK、nirS 和nosZ 基因拷贝数，说明氮肥过量施入使得反硝化相关功能基因丰度减小；UNPKS的nirS基因拷贝数显著低于UCK和UNPK处理，说明秸秆添加显著减少nirS 基因拷贝数；与此同时UCK、UNPK、UNPKS 处理的nosZ 基因拷贝数相对于nirK 和nirS 均具有较高水平。60%WHC 下，氮过量施肥均显著降低nirK、nirS 和nosZ 基因拷贝数，推荐施肥配施秸秆处理显著降低nirS基因拷贝数。淹水下培养结束后，不同处理nirK、nirS 和nosZ 基因拷贝数分别为1.01×107~1.18×107、3.41×107~5.43×107、7.06×107~8.65×107copies·g-1（表4）。与FCK 相比，FNPK、FNhPK 和FN⁃PKS 处理的nirK 基因拷贝数分别提高15%、5%和1%，FNPK、FNhPK、FNPKS 处理的nirS 基因拷贝数分别降低60%、52%、47%；FNPK 和FNhPK 的nosZ 基因拷贝数分别降低7%和7%，而FNPKS 处理提高12%。FCK 处理的nirS 基因拷贝数显著高于其他处理，说明肥料投入和秸秆配施都显著降低nirS 基因拷贝数。总体而言，淹水条件下各施肥处理nirS基因拷贝数是60%WHC的1.19~2.00倍。

3 讨论

3.1 长期施肥处理对土壤理化性质的影响

施肥是提升土壤肥力的有效措施，而土壤理化性质是评价土壤是否健康的重要指标[30]。本研究结果表明，长期不同施肥处理降低轮作土壤pH（表3），这与方凯等[31]的研究一致，长期化肥和化肥配施秸秆处理均降低土壤pH。邹湘等[32]的研究也表明，与不施肥相比，长期NPK 和NPKS 施肥处理也显著降低土壤pH。这可能是由于氮肥的施入，虽然可以暂时提高土壤pH 值，但土壤中NH+4-N 转化成亚硝态氮时，释放大量的H+，长期氮肥投入造成土壤酸化，进而土壤pH 降低。化肥配施秸秆降低土壤pH，可能由于秸秆中含有微生物代谢所需的营养物质，微生物代谢过程产生有机酸，导致土壤酸化[33]。长期NPKS 处理提高SOM 和TN 含量（表3），这与方凯等[31]和郝耀旭等[34]的研究结果一致，化肥配施秸秆处理提升SOM 和TN 含量。然而本研究中，长期NPK处理也显著增加了SOM含量（表3），这可能是因为无机氮肥的施入在增加作物产量的同时，也增加了作物凋落物、根茬残体和根分泌物，在刺激土壤微生物生长的同时，促进了土壤腐殖质的分解，故而SOM 含量增加[35]。但本研究结果显示，长期NPK处理并未显著增加土壤TN含量，这就导致长期NPKS处理C/N含量低于NPK处理。

3.2 轮作系统土壤各施肥处理N2O 排放及相关功能基因对水分的响应

土壤水分是调控土壤硝化和反硝化作用相对强度与N2O 排放的重要因素[38]。硝化作用通常偏好于好氧环境，而反硝化作用更偏好于厌氧环境[39-40]。本研究表明，淹水显著降低各施肥处理土壤的硝化作用强度，而且显著增加N2O 排放（图1、图2 和图3）。这与李平等[41]的研究结果相似，即淹水显著抑制硝化作用的进行，但显著增加N2O 的排放。本研究中，FCK、FNPK、FNhPK、FNPKS各处理和含量变化幅度均显著小于UCK、UNPK、UNhPK、UNPKS，说明淹水显著降低硝化作用的强度；但培养第32 天时FCK、FNPK、FNhPK、FNPKS 处理土壤的AOA 和AOB基因拷贝数分别是UCK、UNPK、UNhPK、UNPKS 处理土壤的1.80、2.34、2.49、1.84 倍和1.19、2.19、1.98、1.53倍，淹水各施肥土壤氨氧化微生物功能基因拷贝数远大于60%WHC 各施肥处理，这可能由于土壤中氨氧化微生物具有耐低氧条件的特性，AOA 和AOB 拷贝数较大，但氨氧化微生物活性较低[42]。有研究指出，当土壤水分大于60%WHC 时，硝化作用的强度逐渐降低，反硝化作用的强度随水分增加而增加[43]。土壤水分含量主要通过影响含nirK 和nosZ 基因微生物的群落结构，刺激部分小丰度微生物的生长，调控N2O的排放强度[28]。但本研究第32 天培养完毕时淹水处理各施肥土壤nirS 基因拷贝数显著高于60%WHC 对应的各施肥处理，这可能由于土壤C/N 及其对施肥的响应能力差异所导致。本研究中淹水各施肥土壤（nirK+nirS）/nosZ比值显著大于60%WHC对应的各施肥处理，淹水各施肥土壤透气性降低，O2向土壤中扩散被表层淹水限制，从而促进了反硝化作用，增加了土壤N2O排放[44-45]。

有研究表明，施用氮肥会增加土壤N2O排放[46-47]，秸秆添加降低土壤N2O 排放[48-49]。但本研究结果显示，UCK 处理N2O 累积排放量均大于其他施肥处理，这可能是因为长期不施肥改变了土壤质地，降低了土壤透气性，随着土壤中少量的O2被消耗，形成反硝化作用更喜好的厌氧环境，从而触发了反硝化作用的发生，从（nirK+nirS）/nosZ 比值角度也可看出，UCK 处理的（nirK+nirS）/nosZ 比值虽然低于FCK，但显著高于UNPK、UNhPK 和UNPKS 处理。淹水各施肥处理N2O排放没有显著差异，这与汤宏等[50]和朱启林等[51]的研究结果均不一致。淹水条件下N2O 的排放主要由反硝化作用主导[23]。本研究中培养第32 天时FNPK、FNhPK、FNPKS处理的（nirK+nirS）/nosZ 比值无显著差异，这表明淹水时施肥处理反硝化作用强度无显著差异，导致N2O 排放量无显著差异。李彬彬等[52]的研究表明，土壤起始C/N越低，N2O排放量越高。FCK处理N2O 排放量较高主要是CK 处理土壤C/N 较低所致。60%WHC 水分条件下，推荐施肥配施秸秆处理显著降低N2O 排放量，这与柴凯斌[14]的研究结果一致。可能原因是：（1）秸秆还田增加了土壤中SOM 含量，为微生物提供了能量，可固定更多的无机氮[53-54]；（2）秸秆还田刺激土壤中的矿质氮微生物活性，进而减少硝化和反硝化作用的底物[55-56]；（3）培养32 d 时，UNPKS处理的（nirK+nirS）/nosZ 比值显著低于其他处理，UN⁃PKS 处理减弱土壤反硝化作用，减少N2O 的排放。综上所述，淹水抑制了不同施肥的响应，且增加了土壤N2O的排放。

4 结论

（1）本研究发现，长期施肥降低土壤pH，推荐施肥配施秸秆处理增加SOM 和TN 含量；相比于推荐施肥，高氮施肥和推荐施肥配施秸秆处理降低土壤C/N。

（2）与不施肥相比，60%WHC 条件下，推荐施肥提高AOA 基因拷贝数，高氮施肥提高AOA 基因拷贝数，降低nirK、nirS 和nosZ 基因拷贝数，推荐施肥配施秸秆提高AOB 基因拷贝数，降低nirS 基因拷贝数；淹水条件下，3 种施肥处理均提高AOA 和AOB 基因拷贝数，均降低nirS的基因拷贝数。

（3）淹水条件剧烈增加N2O 排放，对温室气体减排造成不利影响，而且抑制了不同施肥的响应，尤其降低秸秆添加的减排效果。

（4）综合两种水分条件对不同定位长期施肥轮作土壤N2O 排放的影响，水稻秸秆还田对烟-稻轮作系统中烤烟季温室气体减排具有积极意义。