碳足迹是指某项活动或某个产品在其整个生命周期内直接或间接产生的二氧化碳(CO2)排放总量, 以二氧化碳当量表示(CO2-eq)[1-2]。它是一种用来评估某生产过程或活动在整个生命周期内的温室气体排放的工具。农业中的碳足迹能够系统地评价农作物在生产过程农资投入(柴油、电力、化肥、农药、农膜)引起的间接碳排放和农田直接碳排放总量[3-4], 对实现农业的清洁生产具有重要理论指导意义。

众所周知稻田生态系统是大气中CH4和N2O的主要排放源。我国作为主要水稻种植区, 每年稻田CH4和N2O总排放量的大约分别为7.7—8.0 Tg CH4和88.0—98.1 Gg N2O [5-6]。在未来十年, 随着人口的增加, 对粮食作物水稻的产量需求必然增加, 因此, 将消耗更多的能源、化学肥料、杀虫剂和塑料薄膜[7], 这将直接和间接地促进农田温室气体排放。如何实现稻田生态系统的粮食作物生产低碳耗, 成为科研工作者所努力的方向。近年来, 水资源短缺, 灌溉能源成本高, 稻米价格波动, 以及对玉米等牲畜饲料需求的增加, 促使农民从水稻连作种植转向春玉米-晚稻水旱轮作种植[8-9]。在亚洲地区玉米-水稻轮作系统种植面积超过了3 Ma hm2[9]。联合国粮农组织把春玉米-晚稻模式作为一种对农业集约化生产具有可持续性和战略意义的模式[10]。随着春玉米在我国南方双季稻区大面积推广, 春玉米-晚稻模式必将成为稻田种植业战略性调整的重心。相比于小麦—晚稻、早稻-晚稻和双季玉米, 春玉米-晚稻具有增产潜力大, 自然资源利用率高和经济效益高等优势[9, 11-12]。在全球气候变暖的背景下, 有必要探讨春玉米-晚稻模式对农田碳足迹的影响。目前众多学者对稻田系统碳足迹的研究主要集中于水稻连作系统及轮作模式的研究集中在冬小麦-夏玉米。Xue等[13]评价了我国南方地区双季稻的碳足迹, 并对不同农业管理措施下双季稻的碳足迹及优化进行了评估。Zhang等[14]研究了耕作方式对华北平原冬小麦-夏玉米碳足迹的影响。但是, 目前对春玉米-晚稻轮作模式的碳足迹研究较少。研究表明我国水稻生产过程的温室气体排放和碳足迹高于玉米和小麦[15]。基于此, 本研究提出以下假设, 相对于早稻-晚稻连作模式, 春玉米-晚稻轮作模式不仅可以降低农田温室气体排放而且可以降低产量尺度的碳足迹。因此, 本研究利用生命周期法, 以早稻-晚稻连作模式作为对照, 研究春玉米-晚稻轮作模式对稻田系统碳足迹的影响, 以期获得春玉米-晚稻轮作的农田生态系统的碳排放清单, 为稻田生产系统的低碳清洁化生产以及减缓气候变化提供有效的科学依据。

本试验2017年在浙江省杭州市农业科学研究院试验田(120.16°E, 30.13°N)实施。试验田土壤类型属于粉粘壤土, 土壤全氮含量为2.32 g/kg, pH值为5.52, 总有机碳含量为22.2 g/kg。试验开始前试验田主要种植模式是水稻连作(>20年)。

本试验方案设置了两个作物种植模式处理：早稻-晚稻, 春玉米-晚稻, 每个处理重复3次, 每一个试验小区的大小为8 m× 8 m, 完全随机区组排列。春玉米品种为美玉(加甜糯)7号, 早稻和晚稻的品种分别为甬籼115及秀水134。

早稻采用直播方式, 从2017年4月11日播种, 2017年7月20收获; 第二季晚稻采用移栽方式, 移栽时间为2017年7月22日, 2017年11月19日收割。早、晚稻磷肥(过磷酸钙, 41 kg P2O5/hm2)和钾肥(氯化钾, 47 kg K2O/hm2)全部基施, 氮肥(尿素, 225 kg N/hm2)分为3个阶段早晚稻季将氮肥按3:4:3(基肥, 一次追肥和二次追肥)的比例施到试验组中, 分别是分蘖肥(占总量的30%)、孕穗肥(占总量的40%)以及抽穗肥(占总量的30%)。为防止小区之间水肥串流, 小区之间采用PVC管和薄膜隔离。灌溉方式采用连续灌溉方式。早稻和晚稻收获时每个小区各选取3个1 m2测产。

春玉米采用厢作, 厢宽和沟宽分别为1.3 m和0.3 m, 每个小区6厢, 每厢播种两行春玉米。春玉米采用移栽方式, 三叶期(2017年4月13日)开始移栽, 每行移栽23株。玉米季磷肥(过磷酸钙, 75 kg P2O5/hm2)和钾肥(氯化钾, 47 kg K2O/hm2)全部基施, 氮肥(尿素, 225 kg N/hm2)分为两次施肥, 分别是苗期肥(4月23日施, 占总量40%), 喇叭口肥(6月3日追肥, 占总量60%); 2017年7月9日作为鲜食玉米收获。每个小区分别选取3和4厢中的两个连续20株玉米果穗, 烘干脱粒称重, 并计算含水量。晚稻种植同1.2.1之晚稻。

两种模式的农资投入数据通过实地调研和本研究投入数据的平均值获得, 共调研当地农户37户。调研内容包括水稻或春玉米生产过程中农资投入情况：肥料和农药的投入、农业机械油耗(播种、耕作和粮食运输)和灌溉的电力消耗。玉米、早稻和晚稻的农资投入如表 1。

根据PAS2050[16]的定义, 农业生产碳足迹计算为基于生命周期评价的产品体系中的所有直接和间接温室气体排放量总和, 以CO2当量(CO2-eq)表示。因此, 在本研究中, 玉米和水稻生产的碳足迹是通过农业投入物原材料的初始生产到农场大门的最终产品(如肥料、农药、机械、电力灌溉及收获的稻米和玉米)的间接温室气体排放及稻田直接温室气体排放来确定的。温室气体排放量包括农资投入, 稻田非CO2即CH4和N2O排放量。由于种子在生产过程中产生的碳足迹相对较小, 往往不足农作物生命周期中整体碳足迹的千分之一, 故在此不将其纳入计算范围; 因此, 本研究将水稻或玉米生产的温室气体排放量包括如下：(1)化肥(N, P和K)和农药的生产, 加工和运输及其施用引起的碳排放; (2)耕种和收获机械作业的柴油燃料消耗的碳排放; (3)用于灌溉的电力消耗的碳排放; (4)水稻或春玉米整个生育期稻田CH4和N2O排放。

农资投入(如农药、化肥、电力、柴油等消耗)的温室气体排放量使用下列公式估算。

式中, CFinput是指相关农资投入品所导致的温室气体排放总和(kg CO2-eq/hm2); i是指某一种农资投入品; Ai是指第i个单个农业投入的强度或数量(农药化肥, kg/hm2; 电力, kwh/hm2; 柴油, L/hm2); δi是指制造或应用投入时第i个单个农资投入的温室气体排放系数。排放系数通过数据开放性软件eBalance v3.0(IKE Environmental Technology CO., Ltd., China)获得, 包括Chinese Life Cycle Database 0.7(CLCD 0.7)和Ecoinvent 2.2, 各农资生产投入的温室气体排放因子见表 1。

使用密闭箱气体收集法进行温室气体的采集, 密闭箱的制作材料为PVC板, 分为底座和主箱体两部分, 主箱体外部铺一层铝箔纸进行反光, 以防止夏秋季节采样时主箱体长时间暴露在太阳光中而造成了采样箱内部的气体温度发生较大波动, 从而影响采样结果的有效性。静态箱的尺寸大小(长×宽×高)为50 cm × 50 cm × 90 cm, 以确保水稻完全成熟以后还能将其全部装在静态箱中进行温室气体采集; 底座的尺寸大小(长×宽×高)为50 cm × 50 cm × 30 cm, 底座安装时要埋入地下10 cm深, 且在水稻生产的过程中一直保留, 内部种植9株水稻或4株玉米。采样时把静态箱体盖在每一个试验小区的固定底座上, 同时采用水封的方法来确保箱子内部处于完全密闭状态, 不与外界的进行气体交换。为防止测定区域和非测定区域水肥串流, 用塑料薄膜隔开。从播种或移栽后第5天开始持续到收获前一周每周采集测定一次温室气体(CH4和N2O)通量[17], 采样时间在9:00—11:00之间, 在40 min里每间隔10 min采集一次气体样本, 使用容量为60 mL的注射器从三通口管道中抽取气体, 并立即将其注射到100 mL的铝膜采气袋中, 完成收集样本后需在一周内进行样本的测量。

温室气体(CH4和N2O)采用气相色谱仪(Agilent 7890A, 中国上海)检测。然后根据采集得到的4个气体样本(0、10、20、30 min)的浓度进行线性回归模型测算得到气体样本的浓度增长速率, 在经过温室气体释放速率公式[18]计算得到温室气体释放速率。

式中, F是CH4或N2O的排放速率(mg m-2 min-1), dC/dt是40 min内温室气体浓度的变化速率, m是气体的分子量(g/mol), P是大气压强(Pa), H是静态箱体的高度(m), T是静态箱内温度(℃), R是气体常数(为8.134 g mol-2 k-1), 273.15是绝对温度(K)。

稻田中累积的温室气体(CH4和N2O)碳足迹由公式(3)计算得到。

式中, EC是春玉米或晚稻积累CH4和N2O的CO2排放当量(kg CO2-eq/hm2); α代表潜在影响全球变暖的系数：CH4(34)和N2O(298)[19]。Fi是CH4或N2O第i次取样时间段排放速率(mg/m2/h); Di是第i次和第i-1次采样之间的时间间隔(d); 1440是分钟(min)转换为天(d)的系数; 1/100是mg/m2转换为kg/hm2的系数。

采用单位作物产量碳足迹(kg CO2-eq/kg)评价法, 即稻田系统生产每千克玉米或水稻所导致的碳排放总和, 计算公式如下。

式中, CFa指单季作物(春玉米或水稻)在生长期的碳足迹(kg CO2-eq/kg籽粒产量); EN2O和ECH4是指将CH4和N2O累积排放量转换成的CO2-eq(kg CO2-eq/hm2); Y指当年玉米或水稻的籽粒产量(kg/hm2); CFb是指玉米-晚稻或早晚-晚稻模式的碳足迹(kg CO2-eq/kg籽粒产量), ∑Ei指春玉米-晚稻或早晚-晚稻模式所有温室气体排放量和(kg CO2-eq/hm2); ∑Yi春玉米-晚稻或早稻-晚稻模式所有籽粒产量和(kg/hm2)。

采用Excel 2013软件对数据进行处理, 应用OriginPro 9.4软件对数据进行绘图。

与早稻-晚稻模式碳排放(17848 kg CO2-eq/hm2)相比, 春玉米-晚稻轮作模式降低碳排放6724 kg CO2-eq/hm2; 春玉米比早稻少排放6229 kg CO2-eq/hm2; 与早稻-晚稻模式中晚稻碳排放相比, 春玉米-晚稻轮作模式降低了晚稻碳排放497 kg CO2-eq/hm2(表 2)。

早稻-晚稻连作模式的碳足迹为1.32 kg CO2-eq/kg籽粒产量, 比春玉米-晚稻轮作模式(0.86 kg CO2-eq/kg)高0.56 kg CO2-eq/kg籽粒产量。春玉米的碳足迹(0.50 CO2-eq/kg籽粒产量)较早稻的碳足迹(1.27 CO2-eq/kg籽粒产量)降低0.77 CO2-eq/kg籽粒产量。对于两种模式晚稻碳足迹而言, 春玉米-晚稻模式降低了晚稻碳足迹0.24 CO2-eq/kg籽粒产量(表 3), 同时该模式提高了晚稻的产量(14.8%)。

早稻-晚稻种植模式的碳足迹主要来源于CH4、氮肥使用及灌溉电力消耗的碳排放。在早稻-晚稻种植系统中, CH4的碳排放为9776 kg CO2-eq/hm2(54.8%), 氮肥的碳排放为2871 kg CO2-eq/hm2(16.1%), 灌溉电力消耗的碳排放为2849 kg CO2-eq/hm2(16.0%)。早稻和晚稻的碳足迹组成情况与早稻-晚稻模式碳足迹相似(图 1)。

春玉米-晚稻轮作模式的碳足迹主要来源CH4、氮肥和灌溉电力消耗的碳排放。在该模式中, CH4的碳排放为4442 kg CO2-eq/hm2(39.9%), 氮肥的碳排放为2871 kg CO2-eq/hm2(25.81%), 灌溉电力消耗的碳排放为1508 kg CO2-eq/hm2(13.6%)。晚稻的碳足迹组成情况与春玉米-晚稻模式的碳足迹相似。但是, 对于春玉米而言, 其碳足迹主要来源为氮肥, N2O和CH4。春玉米的氮肥碳排放为1435 CO2-eq/hm2(50.1%), N2O的碳排放为579 kg CO2-eq/hm2(20.2%), CH4的碳排放为378 CO2-eq/hm2(13.2%)(图 1)。

本研究早稻和晚稻的碳足迹为1.14—1.38 kg CO2-eq/kg籽粒产量(表 3), 这与前人研究结果相似[15, 20-21]。Xu等[17]研究表明我国长江下游地区水稻的碳足迹为1.34 kg CO2-eq/kg籽粒产量; Cheng等[11]利用国家统计数据发现我国水稻的平均碳足迹为1.36 kg CO2-eq/kg籽粒产量。但是, 本研究春玉米的碳足迹低于Cheng等[15]的研究结果, 可能主要归因于数据的来源, 系统边界的定义, 农业投入的排放因子以及碳足迹计算方法的差异。春玉米-晚稻水旱轮作模式较早稻-晚稻连作模式降低了0.56 kg CO2-eq/kg, 这也印证本研究的假设。主要归因于较早稻而言春玉米的温室气体排放较低, 其中包括CH4低排放和灌溉电力的零消耗(在春玉米种植期间充足的雨水使得春玉米无需灌溉)。Cheng等[15]结果也表明玉米生产的碳足迹低于水稻, 其认为推广玉米种植将有利于我国农业的低碳发展, 若玉米替代水稻将每年减少100 Mt CO2-eq排放。此外, 春玉米-晚稻模式提高了晚稻产量, 比早稻-晚稻连作模式中晚稻的产量高14.8%。可能的原因是种植春玉米期间土壤好氧的条件加速了土壤碳氮循环, 导致微生物对土壤有机质分解的速度加快, 这有助于促进土壤养分的释放, 提高了后季稻养分利用率[22]。钾素在土壤中活性和流动性较大, 是影响晚稻产量的一个重要因素[23]。尽管春玉米季施钾量高于早稻季, 但是两者收获后土壤的全钾和速效钾含量没有显著性差异(P>0.05, 数据没有在文中展示), 因为钾素主要是当季利用为主[23]。然而增产的原因是否是干湿交替引起的稻田土壤物理化学性质或者生物学特性的交替变化引起的, 还需进一步研究。本研究结果表明早稻-晚稻模式高碳排放, 比春玉米-晚稻模式高60.5%, 同时早稻-晚稻模式产量潜力(13473 kg/hm2)比春玉米-晚稻模式(12965 kg/hm2)高4.00%。主要归因于本研究的春玉米产量较低(5694 kg/hm2), 原因有2个方面：1)试验地地处南方地区, 玉米生长期降雨量大, 土壤中的氮素容易流失, 导致土壤中的氮素损失, 从而造成产量低; 2)本研究的玉米品种为甜糯型, 在没有完全成熟时作为鲜食玉米收获, 相对于饲料型玉米[12], 春玉米果穗偏小, 产量偏低。不过通过更换春玉米品种为饲料型和加强春玉米季排水都可以解决春玉米-晚稻产量潜力相对于低的问题。另外, 研究显示相对于双季稻, 春玉米-晚稻可以提高光温资源利用率和产量潜力, 还可以提高经济效益[11-12]。因此, 综上得知相对于水稻连作模式春玉米-晚稻轮作模式具有低碳排放、高资源利用率和高生产效益的优势。

与早稻-晚稻中的晚稻相比, 春玉米-晚稻模式降低了晚稻CH4的排放, 可能是因为春玉米的种植加速了土壤有机碳的矿化[24], 导致了晚稻季土壤氨氮含量较高, 土壤铵态氮含量过高会加速CH4的氧化[25], 进而降低了稻田CH4的排放。此外, 对于两种种植模式中的水稻而言, CH4是水稻生长期间碳排放的主要成分, 对总温室气体的贡献率为49.2%—57.0%(图 1), 这与Xue等[21]研究类似。本研究结果表明减少稻田的CH4排放对于减少水稻生产的碳足迹至关重要。研究表明干湿交替灌溉可以有效减少稻田CH4排放, 同时也可以减少灌溉次数和灌溉用水量[26], 从而减少灌溉的电力能源消耗和CH4排放。另外, 与传统耕作相比, 免耕是可以显著降低稻田CH4的排放[27-28], 同时也可减少机械柴油的能源消耗。

春玉米-晚稻轮作模式的春玉米碳足迹清单中, 氮肥引起的碳足迹所占的比重最大, 超过了50%(图 1)。Cheng等[29]通过国家统计数据报道了氮肥施用产生的碳足迹占我国农作物生产的碳足迹53.1%—57.3%。基于田间调研数据, Yan等[30]发现我国氮肥引起的温室气体排放占作物生产的碳足迹44.2%—79.6%。主要归因于氮肥生产和运输过程中需要大量消耗的能源, 另外, 为保证作物的生产力往往施用大量的氮肥, 从而造成氮肥引起温室气体排放相对量高。众多研究[15, 18, 31]表明减少氮肥的施用量和提高氮肥的利用率是降低作物生产的碳足迹关键所在。然而, 在我国目前家庭农场式农业经营模式背景下, 通过降低氮肥的施用和提高氮肥利用率来实现农田系统的低碳足迹是一项缓慢艰巨的任务[32]。目前我国正处于工业供给侧改革去产能的时期, 而通过政府驱动政策改善氮肥的生产工艺, 降低氮肥引起的温室气体排放是一项行之有效的措施。Zhang等[32]通过模拟发现如果氮肥生产采用新的工艺技术将每年减少氮肥引起的温室气体排放的20.0%—63.7%, 相当于102—357 Tg CO2-eq。如果我国氮肥的碳排放系数(6.38 kg CO2-eq/kg N)降到美国的水平(3.87 kg CO2-eq/kg N), 两种种植模式的碳排放均将降低1130 kg CO2-eq/hm2, 碳足迹降低8.38%(早稻-晚稻模式)和8.71%(春玉米-晚稻模式)。

综上, 与水稻连作相比, 春玉米-晚稻轮作模式可降低农田系统的碳足迹, 同时可提高晚稻的产量, 是一种具有推广潜力的种植模式。另外, 减少氮肥生产过程的能源消耗和在晚稻种植过程中采用间歇灌溉是降低春玉米-晚稻模式的碳足迹的关键策略。

本研究碳足迹评价参照PAS2050标准[16], 该标准并不包括由土壤有机碳变化引起的碳排放。然而, 农田系统既是温室气体排放的源也是库。Lu等[33]报道了我国农田土壤碳固持量为101 Tg/a, 相当于化石燃料碳排放的13.3%。考虑到土壤碳固存可抵消部分温室气体排放并减少碳足迹[34], 土壤碳变化是否应该包括在粮食产品碳足迹的计算中需要进一步研究。研究表明相对于早晚-晚稻连作, 春玉米-晚稻模式可以加速土壤碳矿化, 增加碳排放[24]。如果考虑到土壤碳变化, 本研究春玉米-晚稻模式的碳排放可能会升高, 这或许会增加本研究的不确定性。尽管如此, 本研究为稻田系统温室气体减排提供了有效的措施, 为南方稻田地区推广春玉米-晚稻模式提供了理论依据。

(1) 早稻-晚稻连作系统单位面积碳排放为17848 kg CO2-eq/hm2, 单位产量的碳足迹为1.32 kg CO2-eq/kg; 春玉米-晚稻轮作系统单位面积碳排放为11123 kg CO2-eq/hm2, 单位产量的碳足迹为0.86 kg CO2-eq/kg。尽管春玉米-晚稻轮作模式的产量潜力比早稻-晚稻连作模式低3.84%, 但是早稻-晚稻连作模式碳排放却比春玉米-晚稻轮作模式高60.5%。

(2) 早稻-晚稻连作系统粮食生产的碳足迹主要来自于稻田CH4的排放(54.8%); 春玉米-晚稻轮作系统粮食生产的碳足迹主要来自于CH4的碳排放(39.9%)。但是, 春玉米生产的碳足迹主要来自氮肥的生产和施用(50.1%)。

(3) 春玉米-晚稻轮作模式是降低连作稻田高碳排放的一种农田管理方式。