劉松，王效琴,2，胡繼平，李強，崔利利，段雪琴，郭亮,2

?

施肥與灌溉對甘肅省苜蓿碳足跡的影響

劉松1，王效琴1,2，胡繼平1，李強1，崔利利1，段雪琴1，郭亮1,2

（1西北農林科技大學資源環(huán)境學院，陜西楊凌 712100；2農業(yè)部西北植物營養(yǎng)與農業(yè)環(huán)境重點實驗室，陜西楊凌 712100）

【目的】評估甘肅省苜蓿不同種植模式的碳足跡，探明溫室氣體的主要排放環(huán)節(jié)。分析施肥與灌溉對甘肅省苜蓿碳足跡的影響，評估可能的減排潛力，為甘肅省苜蓿生產的可持續(xù)發(fā)展提供理論依據?！痉椒ā繎蒙芷谠u價理論和IPCC（2006）田間溫室氣體計算方法，建立苜蓿碳足跡評估方法。通過調研甘肅省苜蓿主產區(qū)隴東、隴中和河西地區(qū)10個縣區(qū)的苜蓿種植農戶和農場，收集苜蓿的產量、化肥、灌溉等生產數(shù)據。根據施肥和灌溉水平及灌溉水源將甘肅省苜蓿種植模式分為4種，使用建立的苜蓿碳足跡評估方法和甘肅省苜蓿生產投入/產出數(shù)據，分析4種種植模式的碳足跡構成特點、施氮水平與灌溉對苜蓿干物質產量和碳足跡的影響特征。使用情景分析方法揭示氮肥減施、改進氮肥生產工藝、無機有機肥混施、滴灌和噴灌技術降低苜蓿種植過程溫室氣體排放的潛力?！窘Y果】甘肅省苜蓿4種種植模式的碳足跡（以CO2eq計）從小到大依次為0.02（不施肥不灌溉模式，NFNI）、0.19（施肥不灌溉模式，SFNI）、0.22（施肥+河水灌溉模式，SFRI）、0.64（施肥+井灌模式，SFWI）kg CO2eq·kg-1苜蓿干物質（DM）。SFNI模式除與SFRI模式之間碳足跡差異不顯著外，與其他種植模式之間的差異均顯著。不同種植模式之間的碳足跡構成環(huán)節(jié)和各環(huán)節(jié)對碳足跡的貢獻率存在差異。NFNI模式的碳足跡主要由苜蓿殘茬和農機使用兩部分的排放組成；SFNI模式和SFRI模式碳足跡的主要產生環(huán)節(jié)是化肥生產和氮肥田間施用的排放，其次是農機使用；SFWI模式碳足跡的最大來源是灌溉耗電，其次為化肥生產和氮肥田間施用的排放。通過氮肥減施、無機有機肥混施、降低氮肥生產過程中排放可使甘肅省SFNI、SFRI和SFWI模式的苜蓿碳足跡分別降低10.0%—18.0%、-3.0%—8.0%和1.8%—5.8%。如果不考慮節(jié)水管材生產的額外溫室氣體排放，節(jié)水灌溉（噴灌和滴灌）可減少SFWI模式苜蓿碳足跡的12.7%—38.5%。【結論】甘肅省4種苜蓿種植模式的產量和碳足跡均存在差異，高投入高產出的SFWI模式的苜蓿產量最高，但碳足跡也顯著高于其他模式。除NFNI模式外，其他3種模式均存在過量施肥現(xiàn)象。減施氮肥和降低氮肥生產過程中的溫室氣體排放均可降低甘肅省苜蓿生產的碳足跡；無機有機肥混施可以降低SFNI和SFRI模式的碳足跡，但同等施氮水平下，無機有機肥混施短期內會降低產量。因此，在保證一定產量同時降低溫室氣體排放量的目標下，有機肥和化肥混施的最佳比例及實際減排潛力仍需通過田間的長期試驗進一步驗證；節(jié)水灌溉是SFWI模式的主要減排措施，但節(jié)水灌溉所能帶來的綜合減排潛力仍需針對具體區(qū)域的田間節(jié)水試驗和額外耗材帶來的溫室氣體排放進一步評估。

苜蓿；碳足跡；施肥；灌溉；甘肅省

0 引言

【研究意義】苜蓿被譽為牧草之王，是優(yōu)質的奶牛粗飼料[1]。隨著近幾年消費者對牛奶品質以及消費需求的提高，中國奶牛飼料正在由秸稈+精飼料型向牧草+青貯+精飼料型轉變[2]，對苜蓿的需求量日益增加，推廣高產優(yōu)質苜蓿種植已成為中國振興奶業(yè)的一個重要舉措[3]。實現(xiàn)苜蓿的高產和優(yōu)質多數(shù)情況需要肥料和灌溉，但是肥料生產、施用以及灌溉過程會消耗大量能源，增加溫室效應。在全球控制溫室效應的大背景下，施肥和灌溉是否會以及如何影響苜蓿生產的溫室氣體排放也是目前所必須考慮的。甘肅省是中國苜蓿傳統(tǒng)種植區(qū)域，多年苜蓿留床面積約占全國種植面積的1/3，穩(wěn)居全國之首[4]，且因其境內地域和降雨差異，苜蓿種植存在不同的施肥和灌溉模式。因此以甘肅省為研究區(qū)域，評估施肥和灌溉對苜蓿碳足跡的影響，對于尋求苜蓿種植的最佳模式，促進其可持續(xù)發(fā)展具有重要意義。【前人研究進展】“碳足跡”是生產單位目標產品的過程中，直接或間接排放的溫室氣體的總和[5]，通常用生產單位產品的二氧化碳當量（CO2eq）表示。農業(yè)系統(tǒng)的碳足跡研究主要用于定量評估不同農產品對全球暖化的貢獻、溫室氣體的主要排放環(huán)節(jié)及減少其排放量的有效措施。國內外對農作物碳足跡的研究主要集中在主糧作物和油料作物。史磊剛等[6]評價了華北平原小麥-玉米的碳足跡，并對不同農作措施下冬小麥-夏玉米的碳足跡及優(yōu)化潛力進行了評價；Gan等研究了大麥、油菜等碳足跡，指出了氮肥施用量與作物碳足跡呈正相關關系[7-8]，劉松等[9]同樣發(fā)現(xiàn)氮肥施用量與小麥和玉米的碳足跡呈顯著正相關關系。然而關于苜蓿的碳足跡國內外只有少數(shù)文獻報道[10]，劉松等[9]研究指出關中平原苜蓿的碳足跡為0.38 kg CO2eq·kg-1苜蓿干草，主要貢獻環(huán)節(jié)是化肥生產、氮肥田間施用的排放和灌溉，但由于樣本較少而沒有分析化肥及灌溉與碳足跡的關系；Adom等[11]研究表明美國生產的苜蓿干草碳足跡僅為 0.17 kg CO2eq，碳足跡主要來源為作物殘茬、磷肥施用、石灰、柴油和電，苜蓿的碳足跡存在較大區(qū)域差異，該研究比較了苜蓿與其他作物之間的碳足跡差異，但沒有分析主要因素與碳足跡之間的變化趨勢?！颈狙芯壳腥朦c】已有的研究表明，氮肥施用量與小麥、玉米等其他農作物的碳足跡之間存在線性正相關關系，而其對苜蓿碳足跡的影響尚不清楚；苜蓿碳足跡存在區(qū)域差異，該差異可能是由于種植模式、施肥和灌溉量的差異引起。甘肅省作為中國苜蓿的主產區(qū)且省內存在不同的施肥和灌溉模式，而其碳足跡特征及減排潛力尚不清楚?！緮M解決的關鍵問題】采用生命周期方法和IPPC[12]提供的溫室氣體計算方法，定量評估并明確甘肅省苜蓿不同施肥和灌溉模式的碳足跡大小及其主要影響環(huán)節(jié)；分析能夠降低苜蓿碳足跡的施肥和灌溉模式，評估其減排潛力，為甘肅省苜蓿生產的可持續(xù)發(fā)展提供理論依據。

1 材料與方法

1.1 研究區(qū)域與數(shù)據來源

于2016年6—9月調研了甘肅省苜蓿主產區(qū)隴東、隴中和河西地區(qū)10個縣區(qū)14個鄉(xiāng)鎮(zhèn)苜蓿的產量、化肥、農機使用及灌溉等情況，每個鄉(xiāng)鎮(zhèn)調研5—6個農戶或農場，共獲得73個有效樣本數(shù)據，共覆蓋3 063 hm2苜蓿種植面積，約占甘肅省苜蓿種植面積的4%。調研點分布見圖1，分別為慶陽市鎮(zhèn)原縣、白銀市會寧縣、定西市安定區(qū)、金昌市永昌縣、張掖市（甘州區(qū)、民樂縣、臨澤縣）、酒泉市（玉門市、肅州區(qū)、金塔縣），這些市（縣區(qū)）是甘肅省苜蓿主產區(qū)，能夠代表甘肅省的苜蓿生產情況。

調研數(shù)據采用SPSS 18.0、Excel 2013軟件處理分析。

圖1 甘肅省苜蓿調研地點分布

1.2 評估方法

本研究依據生命周期評價的一般程序進行，包括目標定義和范圍、清單分析、碳足跡評估（影響評估）和結果分析[15]。

1.2.1 目標定義和范圍 評估甘肅省苜蓿的碳足跡，即生產1 kg苜蓿干物質（DM）從原材料到苜蓿的整個生產過程中溫室氣體的排放。功能單位為1 kg苜蓿干物質。評估范圍根據苜蓿的生產活動確定，從化肥等的生產到苜蓿的收獲，包括（1）化肥生產過程的二氧化碳排放；（2）氮肥施用引起的氧化亞氮直接和間接排放，以及尿素的田間二氧化碳排放；（3）灌溉能耗帶來的二氧化碳排放；（4）播種和收獲等農機使用引起的二氧化碳排放；（5）苜蓿殘茬引起的氧化亞氮直接和間接排放；（6）農藥和農膜生產引起的二氧化碳排放。

1.2.2 清單分析 清單分析是對定義的產品生命周期階段資源、能源和環(huán)境污染物的排放進行數(shù)據量化分析。本文苜蓿生產過程的農資消耗和產量來自于實地調研數(shù)據，代表了2011—2015年甘肅省苜蓿的生產情況。甘肅省隴東、隴中與河西地區(qū)的氣候與水文條件存在較大差異，加之苜蓿種植的布局、規(guī)模與種植習慣的影響，使得苜蓿種植模式差異也較大。根據施肥和灌溉情況甘肅省苜蓿種植主要存在4種模式：（1）不施肥、不灌溉的粗放模式（non-fertilization and non-irrigation，NFNI）；（2）施肥、不灌溉模式（spreading fertilizer and non-irrigation，SFNI）；（3）施肥、河水灌溉模式（spreading fertilizer and river irrigation，SFRI）；（4）施肥、地下水灌溉的高投入高產出模式（spreading fertilizer and well irrigation，SFWI）。前兩種模式主要存在于隴東和隴中地區(qū)，后兩種模式主要存在于河西地區(qū)，4種種植模式的投入產出見表1。SFRI模式中的部分農戶屬于以河水灌溉為主，少量井水灌溉補充。

表1 甘肅省苜蓿不同種植模式的年投入和產出

同列數(shù)據后不同字母表示差異顯著（＜0.05）。尿素指化肥中尿素用量，化肥N中包含了所有化肥來源的氮

Different lowercases in the same column mean significant difference (＜0.05). Urea refers to the amount of urea in the fertilizer, and the fertilizer N contains all the sources of nitrogen fertilizer

使用的參數(shù)依據“就近原則”，首選國內與本研究區(qū)域一致或相近地區(qū)的參數(shù)，其次是國內其他地區(qū)的參數(shù)，再次選擇國外參數(shù)。苜蓿地上殘余物的干物質含量及殘余物含氮量的計算依據IPCC（2006）國家溫室氣體清單指南提供的方法。氮肥（N）、磷肥（P2O5）、鉀肥（K2O）生產的溫室氣體排放系數(shù)（kg·kg-1）分別為4.77[16]、2.02、0.49[17]；苜蓿施用氮肥田間N2O直接排放系數(shù)（kg N2O-N·kg-1N）和氨揮發(fā)系數(shù)（kg NH3-N·kg-1N）分別為0.0105[18]和0.17[19]；硝態(tài)氮淋溶系數(shù)參考趙榮芳等[20]的研究，選取值為0.25；大氣氮沉降的N2O排放系數(shù)（kg·kg-1）、氮滲漏或徑流損失引起的N2O排放系數(shù)（kg·kg-1）、尿素田間CO2排放系數(shù)（kg·kg-1）和苜蓿殘余物含氮量（kg·kg-1）均來自IPCC（2006），分別為0.01、0.0075、0.2和0.027；柴油、電力、農藥、地膜生產或使用的CO2排放系數(shù)（kg·L-1、kg·kwh-1、kg·kg-1、kg·kg-1）分別為2.76[21]、0.95[22]、18.04[23]、6.41[24]。

1.2.3 碳足跡評估 苜蓿生產涉及的溫室氣體為CO2和N2O，為便于匯總排放量，統(tǒng)一的計量單位二氧化碳當量（CO2eq）被用來計量不同種類溫室氣體的排放量。苜蓿碳足跡計算公式如下：

=E/Y

數(shù)據經核對校驗無誤后，應用SPSS 19.0軟件進行統(tǒng)計分析。計數(shù)資料用百分率(%)表示，組間比較采用χ2檢驗，檢驗水準 α=0.05。

式中，為苜蓿碳足跡（kg CO2eq·kg-1DM）；E為生產苜蓿生命周期內的以CO2為當量計算的總的溫室氣體排放量（kg）；Y為苜蓿干物質產量（kg）。

(F+C×C)×(ef+F×ef)×44/28×265

式中，F為生產苜蓿施用第種肥料的量（kg）；ef為生產第種肥料的CO2排放系數(shù)（kg CO2·kg-1）；為苜蓿農機柴油消耗量（L）；ef為柴油燃燒的CO2排放系數(shù)（kg CO2·L-1）；E為生產第種作物的耗電量（kWh）；ef為生產電的CO2排放系數(shù)（kg CO2·kWh-1）；為苜蓿尿素施用量（kg）；ef為尿素排放因子（kg·kg-1）；44/12為CO2-C換化為CO2的系數(shù)；F為苜蓿施氮量（kg）；C為苜蓿殘余物量（kg）；C為苜蓿殘余物中含氮量（kg N·kg-1）；ef氮肥N2O直接排放系數(shù)（kg N2O-N·kg N-1）；F為苜蓿氨揮發(fā)系數(shù)，ef為大氣氮沉降排放N2O的系數(shù)；44/28為N2O-N轉化為N2O的系數(shù)；265[25]為N2O增溫潛趨勢值。

2 結果

2.1 甘肅省苜蓿不同種植模式的碳足跡

甘肅省苜蓿不同種植模式生產1 kg苜蓿（以干物質計）的溫室氣體排放量（以kg CO2eq計）由小到大依次為NFNI 模式0.021、SFNI模式0.190、SFRI模式0.219、SFWI模式0.639。SFNI模式除過與SFRI模式之間碳足跡差異不顯著外，與其他兩種植模式之間的差異均顯著（表2）。

2.2 甘肅省苜蓿碳足跡的主要構成環(huán)節(jié)

甘肅省苜蓿碳足跡構成如圖2所示，不同種植模式之間的構成環(huán)節(jié)和各環(huán)節(jié)對碳足跡的貢獻率存在差異。NFNI模式的碳足跡主要由苜蓿殘茬和農機使用兩部分的排放組成，其對碳足跡的貢獻率分別為56%和44%。SFNI模式和SFRI模式碳足跡的主要產生環(huán)節(jié)是化肥生產和氮肥田間排放，對苜蓿碳足跡貢獻率別占38%—41%和28%—42%；農機使用次之，貢獻率為13%—16%；其他環(huán)節(jié)的貢獻率均小于10%。SFWI模式碳足跡的最大來源是灌溉耗電，對苜蓿碳足跡貢獻約占67%；其次為化肥生產約占15%，氮肥田間排放占12%，其他環(huán)節(jié)的貢獻率均低于6%。雖然不同種植模式的苜蓿碳足跡不同，但對于施肥型的種植模式，化肥生產、氮肥田間排放、灌溉用電和農機使用是溫室氣體主要的排放環(huán)節(jié)。

表2 甘肅省苜蓿的碳足跡

化肥生產和氮肥田間排放都可歸類為與化肥有關的排放，化肥對苜蓿碳足跡的貢獻率為69%—81%（NFNI、SFWI模式除外），其中由氮肥生產和田間施用的排放占化肥總排放的80%—91%。這主要由于生產單位氮肥的溫室氣體排放量遠高于生產磷肥和鉀肥，約是磷肥的2.4倍，鉀肥的10倍；另外，氮肥施到田間后會排放N2O，而N2O的溫室效應是CO2的265倍，同時如果施用的氮肥是尿素，尿素除了會排放N2O，還會直接釋放CO2，這些都導致氮肥成為影響苜蓿碳足跡的重要因素。

苜蓿生長季需水量大，而河西地區(qū)降雨少蒸發(fā)量大，因此苜蓿生長期需要灌溉大量的水，加上使用機井灌溉需要消耗大量的電，這導致灌溉耗電成為河西井灌區(qū)苜蓿碳足跡的首要來源。根據本文調研數(shù)據，農機已經廣泛應用于甘肅省苜蓿種植。苜蓿地翻耕、播種、施肥、噴灑農藥、收割、摟草、打捆、翻曬等環(huán)節(jié)都需要使用農機，同時甘肅苜蓿一般收割2—4茬，每一茬都需要收割、摟草等，因而苜蓿農機（一般使用柴油）使用柴油消耗較高，使其成為甘肅省苜蓿碳足跡不可忽視的因素。

2.3 氮肥施用量對苜蓿干物質產量和碳足跡的影響

調研獲得的農戶實際氮肥施用量和苜蓿干物質產量之間的變化趨勢如圖3，當?shù)适┯昧浚?0 kg·hm-2時，隨著氮肥施用量的增加，苜蓿干物質產量有增長的趨勢；但當?shù)适┯昧浚?0 kg·hm-2時，干物質產量不再隨著施氮量的增加而增長，即調研區(qū)域苜蓿的最優(yōu)施氮量是70 kg·hm-2。該最優(yōu)施氮水平顯著低于SFNI、SFRI和SFWI模式的平均施氮水平。

甘肅省苜蓿施氮量對碳足跡的影響見圖4，施氮量與苜蓿碳足跡呈正相關線型關系（=0.72）。氮肥對碳足跡的影響主要有2個環(huán)節(jié)，一個是氮肥生產，另一個是氮肥田間施用。減少氮肥對苜蓿碳足跡的影響有2個途徑：合理減施氮肥和降低氮肥生產過程的溫室氣體排放。因此，本研究設置了2個氮肥減排情景：一是調研數(shù)據模擬的最優(yōu)施氮水平（70 kg N·hm-2）；二是氮肥生產工藝改進。氮肥減施方案假設甘肅省3種施肥模式（SFNI、SFRI和SFWI）的平均施氮水平通過合理的田間管理，降到70 kg·hm-2，而產量能保持原來的水平。氮肥生產工藝改進假設中國生產氮肥的排放系數(shù)由4.77 kg CO2·kg-1N降至美國水平3.87 kg CO2·kg-1N[11]。2種減排情景下3種苜蓿種植模式的碳足跡如圖5，平均施氮水平降為70 kg·hm-2時，SFNI、SFRI和SFWI模式的平均碳足跡分別降低10%、18%和11%，改進氮肥生產工藝可使3種模式的碳足跡分別降低5.8%、4.3%和1.8%。

圖3 甘肅省苜蓿氮肥施用量對干物質產量的影響

圖4 甘肅省苜蓿氮肥施用量對碳足跡的影響

FS為調研數(shù)據碳足跡；PINP為改進氮肥生產工藝的碳足跡；N70為最優(yōu)施氮水平的碳足跡

2.4 施用有機肥對苜蓿碳足跡影響

在本文調研中，甘肅苜蓿種植基本上施用化肥，很少施用有機肥。本文根據楊園園等[26]田間試驗研究的單施氮肥和牛糞氮肥混施（60%牛糞N，40%化肥N）的田間N2O排放量和對產量的影響數(shù)據，并考慮牛糞作為有機肥所分擔的奶牛養(yǎng)殖過程溫室氣體排放量[27]，評估牛糞氮肥混施對施肥的3種模式SFNI、SFRI和SFWI碳足跡的影響。根據楊園園等兩年的田間試驗數(shù)據，在同樣的施氮水平下，牛糞氮肥混施（60%牛糞N，40%化肥N）的田間N2O排放系數(shù)比單施氮肥高12%，產量低9.4%[26]，不過牛糞作為有機肥所分擔的奶牛養(yǎng)殖過程溫室氣體排放量2.15 kg CO2eq·kg-1N[27]卻低于生產N肥的溫室氣體排放量4.77。牛糞氮肥混施（60%牛糞N，40%化肥N）使SFNI和SFRI模式的碳足跡分別降低8%和6%，但對于SFWI模式由于其碳足跡的主控因素是灌溉能耗和產量，使用牛糞減少的N肥生產過程的排放對總溫室氣體排放量的影響較小，但卻較大地影響了產量，導致碳足跡比單施氮肥增加了3%。

2.5 灌溉對苜蓿干物質產量和碳足跡的影響

甘肅省年降雨量少而苜蓿生長季節(jié)需水量大，灌溉是苜蓿產量的控制因素，灌溉量與苜蓿干物質產量之間存在正相關線型關系（=0.75，圖6），隨著灌溉量的增加，干物質產量有增長趨勢。苜蓿灌溉對碳足跡的影響主要由灌溉是否消耗能源及消耗多少決定。灌溉耗電量與碳足跡有顯著的線型正相關關系（=0.93，圖7）。甘肅省苜蓿生產的4種模式中只有SFRI和SFWI灌溉，SFRI主要是河水自流灌溉，灌溉耗電量小，其碳足跡受灌溉影響小，SFWI模式為井水灌溉，耗電量大，因此受灌溉量影響大。

本文調研的SFWI模式均為大水漫灌，大水漫灌的水資源有效利用率只有30%—40%[28]，因而苜蓿漫灌浪費了大量的水資源。郭學良等[29]研究表明，與漫灌相比，噴灌產量增加了21.6%，每公頃需灌溉9 100 m3，與SFWI模式漫灌相比節(jié)水2 150 m3，節(jié)約用電1 075 kwh。王東等[30]的研究指出地埋滴灌種植苜蓿，畝產較常規(guī)灌溉高40%以上，每公頃需灌溉水4 500 m3，與SFWI模式漫灌相比可節(jié)水6 750 m3，節(jié)約用電3 375 kwh。若SFWI模式采用噴灌和滴灌，假設產量保持原來的水平，若不考慮額外的管材生產帶來的溫室氣體排放，則河西井灌區(qū)苜蓿碳足跡將降至0.558和0.383 kg CO2eq·kg-1DM，與漫灌相比分別降低了12.7%和38.5%。

圖6 甘肅省灌溉水量對苜蓿干物質產量的影響

圖7 灌溉耗電量對甘肅省苜蓿碳足跡的影響

3 討論

本研究獲得的甘肅省苜蓿生產的SFNI模式和SFRI模式的碳足跡（0.190和0.219 kg CO2eq·kg-1DM）高于美國的研究（0.17 kg CO2eq·kg-1DM）[8]，主要原因為美國苜蓿種植很少施氮肥，其碳足跡主要來源于磷肥施用，而磷肥對碳足跡的影響遠低于氮肥。甘肅省的施肥河水灌溉（SFRI）模式與劉松等[9]研究的關中平原苜蓿種植模式一致，但碳足跡（0.219 kg CO2eq·kg-1DM）低于關中平原（0.382 CO2eq· kg-1DM），主要原因是甘肅省SFRI模式的施氮水平（113 kg N·hm-2）低于陜西省關中平原（179 kg N·hm-2），而氮肥施用量是施肥河水灌溉模式碳足跡的主要因素，對碳足跡的貢獻率為80%—91%。已有文獻未對不施肥不灌溉（NFNI）和施肥井灌（SFWI）模式的苜蓿碳足跡做過研究，無法對甘肅省的這兩類苜蓿種植模式進行同類研究比較。

甘肅存在的4種苜蓿種植模式中，NFNI模式較為粗放，既不施肥也不灌溉，是4種模式中產量最低的，但由于生產過程沒有引入大量的溫室氣體排放源，使得溫室氣體排放量非常少，導致其碳足跡最低。SFNI模式也不灌溉，但施肥，其產量顯著高于前者，而由于施肥帶來的溫室氣體排放量的增量高于產量的增量，導致SFNI模式碳足跡高于NFNI模式，且差異顯著。SFRI模式與SFWI模式均屬于既施肥又灌溉的資源高投入模式，雖然產量都顯著高于前兩者，但高投入帶來的產量增量尚不足以抵消溫室氣體排放量的增量，因此碳足跡高于前兩者，尤其是SFWI模式屬于井水灌溉，抽取地下水灌溉消耗大量電能導致溫室氣體排放量顯著升高，其碳足跡不僅顯著高于前兩者也顯著高于河水灌溉的SFRI模式。

氮肥施用是甘肅省苜蓿生產SFNI、SFRI和SFWI模式碳足跡的主要來源，占前兩類模式碳足跡的80%—91%，通過對實際生產數(shù)據模擬分析，甘肅省多數(shù)苜蓿種植農戶存在過量施肥的現(xiàn)象。研究表明，通過提高氮肥利用效率可以減少氮肥施用量，西方國家氮肥利用率為40%—60%[31]，而中國一般只有30%—35%[32]；另外，苜蓿是固氮作物，適當減少氮肥使用量也有利于促進與充分利用苜蓿的生物固氮作用，這些都表明甘肅省通過減少氮肥施用量從而減少碳足跡仍有較大空間。中國氮肥生產是煤基型的，生產氮肥的溫室氣體排放量高于美國等其他發(fā)達國家[11]，氮肥生產工藝改進往往需要較長時間，因此降低氮肥生產溫室氣體排放系數(shù)來降低甘肅省苜蓿碳足跡需要在較長時間內實現(xiàn)。

甘肅苜蓿種植基本上施用化肥，很少施用有機肥。有研究表明，有機無機肥配合施用能夠促進苜蓿固氮能力，減少氮素損失，提高苜蓿氮肥利用率[33-35]。合理地將畜禽糞便引入到苜蓿生產中，可降低苜蓿生產中的化學氮肥用量，同時也可以實現(xiàn)農牧之間的養(yǎng)分循環(huán)[36]。本文根據楊園園等[26]的試驗數(shù)據的評估結果表明，有機無機肥配合施用可以降低SFNI和SFRI模式的碳足跡，但施用有機肥對土壤質量和作物營養(yǎng)是一種長效機制，在同樣的施氮水平下，用有機肥替代化肥，短期內會降低產量，因此，在保證一定產量同時降低溫室氣體排放量的目標下，如何確定有機肥和化肥混施的最佳比例及實際減排潛力仍需通過田間的長期試驗進一步驗證。

灌溉耗電是井灌區(qū)SFWI模式碳足跡的首要影響因素，井水灌溉還是河水灌溉主要由當?shù)氐乃Y源分布和地形決定，井灌區(qū)改為河水灌溉并不現(xiàn)實，井灌區(qū)應主要依靠節(jié)水灌溉減少碳足跡。本研究結果表明，苜蓿生產采用噴灌和滴灌均可降低碳足跡，尤其是滴灌可顯著降低灌溉耗水量和碳足跡。張愛寧等[37]研究表明，適宜苜蓿生長的最佳灌溉方式為滴灌，對苜蓿的傷害較小，測得葉片中葉綠素含量最高，脯氨酸含量和丙二醛含量最低?？梢钥闯龅喂嗖粌H大量節(jié)約了水資源，還增加苜蓿產量，減少井水灌溉能耗帶來的溫室氣體排放。但實現(xiàn)滴灌需要鋪設管道，管道材料在生產和鋪設過程中也會產生溫室氣體排放，滴灌所能帶來的綜合減排潛力仍需針對具體區(qū)域的節(jié)水試驗和滴灌額外耗材帶來的溫室氣體排放進一步評估。

4 結論

甘肅省苜蓿主要種植區(qū)域的氣候和水文條件存在較大的差異，形成了種植模式和碳足跡的差異。高投入高產出種植模式雖然增加了苜蓿產量，但高投入帶來的產量的增量尚不足以平衡其造成的溫室氣體排放的增量，SFWI模式雖然獲得了4種模式中最高的產量，但其碳足跡也顯著高于其他3種模式。

通過氮肥減施和降低氮肥生產過程中的溫室氣體排放均可降低甘肅省苜蓿生產的碳足跡。在同等施氮水平下，牛糞氮肥混施（60%牛糞N，40%化肥N）雖然在短期內降低了苜蓿產量，卻可以降低SFNI和SFRI模式的碳足跡。

灌溉是井水灌溉區(qū)SFWI模式碳足跡的最大影響因素。如果不考慮管道材料在生產和鋪設過程的溫室氣體排放，通過節(jié)水灌溉（噴灌和滴灌）可減少SFWI模式苜蓿的碳足跡。

[1] 王顯國, 韓建國, 祝美俊. 美國的苜蓿種子產業(yè). 世界農業(yè), 2004(6): 42-44.

Wang X G, Han J G, Zhu M J. Alfalfa seed industry of the United States., 2004(6): 42-44. (in Chinese)

[2] 中國奶業(yè)年鑒編委會. 中國奶業(yè)年鑒2013. 北京: 中國農業(yè)出版社, 2013.

Editorial Committee of China dairy yearbook.. Beijing: China Agriculture Press, 2013. (in Chinese)

[3] 聶迎利, 馮艷秋. 加快實施振興奶業(yè)苜蓿發(fā)展行動為奶業(yè)轉型升級和現(xiàn)代奶業(yè)建設提供有力支撐——全國振興奶業(yè)苜蓿發(fā)展行動現(xiàn)場會在遼寧錦州召開. 中國乳業(yè), 2013(138): 4-7.

Nie Y L, Feng Y Q. Accelerate the implementation of the development of the dairy industry to revitalize the dairy industry to upgrade and promote the transformation and upgrading of modern dairy industry to provide a strong support——National dairy alfalfa development initiative to revitalize the field has been held in Jinzhou, Liaoning., 2013(138): 4-7. (in Chinese)

[4] 朱新強, 王曉力, 王春梅, 張茜. 甘肅省苜蓿種植現(xiàn)狀及成本收益分析. 中國草食動物科學, 2014, 34(6): 63-67.

Zhu X Q, Wang X L, Wang C M, Zhang Q. Analysis on the status quo and cost benefit of alfalfa in Gansu Province., 2014, 34(6): 63-67. (in Chinese)

[5] 王占彪, 王猛, 陳阜. 華北平原作物生產碳足跡分析. 中國農業(yè)科學, 2015, 48(1): 83-92.

Wang Z B, Wang M, Chen F. Carbon footprint analysis of crop production in North China Plain., 2015, 48(1): 83-92. (in Chinese)

[6] 史磊剛, 陳阜, 孔凡磊, 范士超. 華北平原冬小麥-夏玉米種植模式碳足跡研究. 中國人口·資源與環(huán)境, 2011, 21(9): 93-98.

Shi L G, Chen F, Kong F L, Fan S C. The carbon footprint of winter wheat-summer maize cropping pattern on North China Plain., 2011, 21(9): 93-98. (in Chinese)

[7] Gan Y T, Liang C, May W, Malhi S S, Niu J Y, Wang X Y. Carbon footprint of spring barley in relation to preceding oilseeds and N fertilization., 2012, 17(5): 635-645.

[8] Gan Y T, Liang C, Huang G B, Malhi S S, Brandt S A, Katepa-Mupondwa F. Carbon footprint of canola and mustard is a function of the rate of N fertilizer., 2012, 17(1): 58-68.

[9] 劉松, 王效琴, 崔利利, 段雪琴, 趙加磊. 關中平原飼料作物生產的碳足跡及影響因素研究. 環(huán)境科學學報, 2017, 37(3): 1201-1208.

Liu S, Wang X Q, Cui L L, Duan X Q, Zhao J L. Carbon footprint and its impact factors of feed crops in Guanzhong Plain., 2017, 37(3): 1201-1208. (in Chinese)

[10] Baek C Y, Lee K M, Park K H. Quantification and control of the greenhouse gas emissions from a dairy cow system., 2014, 70: 50-60.

[11] Adom F, Maes A, Workman C, Clayton-Nierderman Z, Thoma G, Shonnard D. Regional carbon footprint analysis of dairy feeds for milk production in the USA., 2012, 17(5): 520-534.

[12] IPCC 2006. 2006 IPCC guidelines for national greenhouse gas inventories[EB/OL]. [2017-04-26]. http://www.ipcc-nggip.iges.or.jp/ public/2006gl/vol4.html.

[13] 韓建民. 甘肅苜蓿草業(yè)產業(yè)化發(fā)展研究. 中國農村經濟, 2000(3): 38-42.

Han J M. Study on the development of alfalfa industry in Gansu., 2000(3): 38-42. (in Chinese)

[14] 徐廣, 刁治民, 肖前青, 王堅. 苜蓿深層次開發(fā)現(xiàn)狀及前景. 青海草業(yè), 2008, 17(2): 18-20.

Xu G, Diao Z M, Xiao Q Q, Wang J. deep-rooted alfalfa development statue and prospects., 2008, 17(2): 18-20. (in Chinese)

[15] ISO 14040. Environmental management-life cycle assessment- principles and framework[EB/OL]. [2017-04-26]. https://www.iso.org/ standard/37456.html.

[16] 王效琴, 梁東麗, 王旭東, 彭莎, 鄭金正. 運用生命周期評價方法評估奶牛養(yǎng)殖系統(tǒng)溫室氣體排放量. 農業(yè)工程學報, 2012, 28(13): 179-184.

Wang X Q, Liang D L, Wang X D, Peng S, Zheng J Z. Estimation of greenhouse gas emissions from dairy farming systems based on LCA., 2012, 28(13): 179-184. (in Chinese)

[17] Ecoinvent Centre, 2010. Ecoinveny Data v2.2. Swiss Cebtre for Life Cycle Inventories, Dübendorf, Switzerland[DB/OL]. [2017-04-26]. http://www.ecoinvent.ch.

[18] 張強, 巨曉棠, 張福鎖. 應用修正的IPCC2006方法對中國農田N2O排放量重新估算. 中國生態(tài)農業(yè)學報, 2010, 18(1): 7-13.

Zhang Q, Ju X T, Zhang F S. Re-estimation of direct nitrous oxide emission from agricultural soils of China via revised IPCC2006 guideline method., 2010, 18(1): 7-13. (in Chinese)

[19] Larsson L, Ferm M, Kasimir-Klemedtsson A, Klemedtsson L. Ammonia and nitrous oxide emissions from grass and alfalfa mulches., 1998, 51(1): 41-46.

[20] 趙榮芳, 陳新平, 張福鎖. 華北地區(qū)冬小麥-夏玉米輪作體系的氮素循環(huán)與平衡. 土壤學報, 2009, 46(4): 684-697.

Zhao R F, Chen X P, Zhang F S. Nitrogen cycling and balance in winter wheat-summer maize rotation system in Northern China plain., 2009, 46(4): 684-697. (in Chinese)

[21] 謝鴻宇. 生態(tài)足跡評價模型的改進與應用. 北京: 化學工業(yè)出版社, 2008: 52-53.

Xie H Y.. Beijing: Chemical Industry Press, 2008: 52-53. (in Chinese)

[22] 國家發(fā)展和改革委員會應對氣候變化司. 2015中國區(qū)域電網基準線排放因子[EB/OL]. (2016-09-11)[ 2017-04-26]. http://cdm.ccchina. gov. cn/archiver/cdmcn/UpFile/Files/Default/20160606120244478242.pdf.

National Development and Reform Commission secretary for responding to climate change. China regional power grid baseline emission factor in 2015[EB/OL]. (2016-09-11)[2017-04-26]. http:// cdm.ccchina.gov.cn/archiver/cdmcn/UpFile/Files/Default/20160606120244478242.pdf. (in Chinese)

[23] WEST T O, MARLAND G. A synthesis of carbon sequestration, carbon emissions, and net carbon flux in agriculture: comparing tillage practices in the United States., 2002, 91(3): 217-232.

[24] 陳紅, 郝維昌, 石鳳, 王天民. 幾種典型高分子材料的生命周期評價. 環(huán)境科學學報, 2004, 24(3): 545-549.

Chen H, Hao W C, Shi F, Wang T M. Life cycle assessment of several typical macromolecular materials., 2004, 24(3): 545-549. (in Chinese)

[25] IPCC Working Group I . Climate change 2013: the physical science basis[EB/OL]. [2017-04-26]. http://www.climatechange2013.org/.

[26] 楊園園, 高志嶺, 王雪君. 有機、無機氮肥施用對苜蓿產量、土壤硝態(tài)氮和溫室氣體排放的影響. 應用生態(tài)學報, 2016, 27(3): 822-828.

Yang Y Y, Gao Z L, Wang X J. Impact of organic and inorganic fertilization on alfalfa yield, soil nitrate and greenhouse gas emissions., 2016, 27(3): 822-828. (in Chinese)

[27] Wang X, Kristensen T, Mogensen L, Knudsen M T, Wang X. Greenhouse gas emissions and land use from confinement dairy farms in the Guanzhong plain of China—using a life cycle assessment approach., 2016, 113: 577-586.

[28] 姚素梅, 康躍虎, 劉海軍, 馮金朝, 王君. 噴灌與地面灌溉條件下冬小麥光合作用的日變化研究. 農業(yè)工程學報, 2005, 21(11): 16-19.

Yao S M, Kang Y H, Liu H J, Feng J C, Wang J. Diurnal changes of photosynthesis of winter wheat under sprinkler irrigation and surface irrigation conditions., 2005, 21(11): 16-19. (in Chinese)

[29] 郭學良, 李衛(wèi)軍. 不同灌溉方式對紫花苜蓿產量及灌溉水利用效率的影響. 草地學報, 2014, 22(5): 1086-1090.

Guo X L, Li W J. Effects of different irrigation methods on alfalfa yield and irrigation water use efficiency., 2014, 22(5): 1086-1090. (in Chinese)

[30] 王東, 曹新成, 李富先, 王小莉. 地下滴灌苜蓿栽培的試驗效果分析. 新疆農墾經濟, 2004(5): 65-66.

Wang D, Cao X C, Li F X, Wang X L. Experimental result analyses of underground drip irrigation to alfalfa growth., 2004(5): 65-66. (in Chinese)

[31] Ladha J K, Pathak H, Krupnik T J, SIX J, VAN KESSEL G. Efficiency of fertilizer nitrogen in cereal production: retrospects and prospects., 2005, 87: 85-156.

[32] 巨曉棠, 張福鎖. 中國北方土壤硝態(tài)氮的累積及其對環(huán)境的影響. 生態(tài)環(huán)境, 2003, 12(1): 24-28.

Ju X T, Zhang F S. Nitrate accumulation and its implication to environment in north China., 2003, 12(1): 24-28. (in Chinese)

[33] 張小莉, 孟琳, 王秋君, 羅佳, 黃啟為, 徐陽春, 楊興明, 沈其榮. 不同有機無機復混肥對水稻產量和氮素利用率的影響. 應用生態(tài)學報, 2009, 20(3): 624-630.

Zhang X L, Meng L, Wang Q J, Luo J, Huang Q W, Xu Y C, Yang X M, Shen Q R. Effects oforganic-inorganic mixed fertilizers on riceyield and nitrogen use efficiency., 2009, 20(3): 624-630. (in Chinese)

[34] 彭佩欽, 仇少君, 侯紅波, 李恩堯, 丘亞群.15N交叉標記有機與無機肥料氮的轉化與殘留. 生態(tài)學報, 2011, 31(3): 858-865.

Peng P Q, Qiu S J, Hou H B, Li E Y, Qiu Y Q. Nitrogen transformation and its residue in pot experimentsamended with organic and inorganic15N cross labeled fertilizers., 2011, 31(3): 858-865. (in Chinese)

[35] 朱菜紅, 董彩霞, 沈其榮, 徐陽春. 配施有機肥提高化肥氮利用效率的微生物作用機制研究. 植物營養(yǎng)與肥料學報, 2010, 16(2): 282-288.

Zhu C H, Dong C X, Shen Q R, Xu Y C. Microbial mechanism on enhancement of inorganicfertilizer-N useefficiency for combined use of inorganic and organic fertilizers., 2010, 16(2): 282-288. (in Chinese)

[36] 韓思訓, 王森, 高志嶺, 王殿武. 不同施肥條件下苜蓿產量、氮素累積量及肥料氮素利用率研究. 華北農學報, 2014, 29(6): 220-225.

Han S X, Wang S, Gao Z L, Wang D W. Impact of N application on the yield of alfalfa, N accumulation and nitrogen use efficiency., 2014, 29(6): 220-225. (in Chinese)

[37] 張愛寧, 王玉祥, 隋曉青, 張博. 不同灌溉方式對苜蓿形態(tài)及生理的影響. 中國農學通報, 2014, 30(18): 161-165.

Zhang A N, Wang Y X, Sui X Q, Zhang B. Effects of different irrigation methods on morphology and physiology of alfalfa., 2014, 30(18): 161-165. (in Chinese)

（責任編輯 岳梅）

Effects of Fertilization and Irrigation on the Carbon Footprint of Alfalfa in Gansu Province

LIU Song1, WANG XiaoQin1,2, HU JiPing1, LI Qiang1,CUI LiLi1, DUAN XueQin1, GUO Liang1,2

(1College of Resources and Environment, NorthwestA&F University, Yangling 712100, Shaanxi;2Key Laboratory of Plant Nutrition and Agri-environment in Northwest China, Ministry of Agriculture, Yangling 712100, Shaanxi)

【Objective】The objective of the study is to evaluate the carbon footprint (CF) of alfalfa production in different cropping patterns in Gansu Province, clarify the main process emitting greenhouse gases (GHG), analyze the effects of fertilization and irrigation on the CF, and investigate the potential measures for decreasing CF of alfalfa production in this province. 【Method】Life cycle assessment methodology combined with IPCC (2006) guide for calculating greenhouse gases emission in the field were used to quantify the CF of alfalfa production. Data were collected from farm survey in major regions of alfalfa cultivation in Gansu Province which covered 10 counties in the districts of Longdong, Longzhong and Hexi. Cropping pattern of alfalfa production in Gansu province was classified into four types according to the level of nitrogen application and irrigation, and the irrigation water source. The CF composition characteristics of four cropping patterns and the effect of nitrogen application rate and irrigation on productivity and CF of alfalfa were analyzed based on CF calculation model and input and output data of alfalfa production. Scenario analysis was used to estimate the GHG abatement potential through strategies such as decreasing chemical fertilizer application, improving fertilizer production technique, combining chemical fertilizer and organic fertilizer, and utilizing sprinkle and dripping irrigation. 【Result】CF of alfalfa production was 0.02, 0.19, 0.22 and 0.64 kg CO2eq·kg-1DM under four cropping patterns including NFNI (non-fertilization and non-irrigation), SFNI (spreading fertilization and non-irrigation), SFRI (spreading fertilization and river irrigation) and SFWI (spreading fertilization and well irrigation), respectively. The CF of alfalfa production in SFNI cropping pattern was significantly different from other cropping patterns except for SFRI. The main contributors of CF and their proportion to the total GHG emissions varied among different cropping patterns. For NFNI pattern, CF was mainly composed of GHG emissions from crop residues and the use of agricultural machinery. For SFNI and SFRI patterns, the major contributors of CF were chemical fertilizer production and nitrogen fertilizer application, followed by the use of agricultural machinery. The greatest GHG emissions source for SFWI pattern was electricity for irrigation, followed bychemical fertilizer production and nitrogen fertilizer application. Reducing the amount of applied nitrogen, combining chemical fertilizer with manure, and decreasing the GHG emission from nitrogen fertilizer production would decrease CF by 10.0%-18.0%, -3.0%-8.0% and 1.8%-5.8% in SFNI, SFRI and SFWI patterns, respectively. The water-saving irrigation (sprinkler irrigation and drip irrigation) would decrease the CF by 12.7%-38.5%.【Conclusion】 The alfalfa production and CF of four cropping patterns in Gansu province are different. The SFWI pattern got the highest production with high input, while its CF was significantly higher than that of other patterns. With the exception of NFNI pattern, excessive fertilizer was used in other three patterns. Reducing the amount of applied N fertilizer and the GHG emission from N fertilizer production will decrease CF of alfalfa production in Gansu province. The combination of chemical fertilizer and manure will decrease CF in SFNI and SFRI patterns, but also the yield of alfalfa hay. So, the optimal ratio of chemical fertilizer to manure should be investigated further. The water-saving irrigation will be the major mitigation measure in SFWI pattern, while the integrated abatement potential of GHG emission should be evaluated including emissions from the production of pipes used in sprinkler and drip irrigation.

alfalfa; carbon footprint; fertilization; irrigation; Gansu Province

2017-04-26：2017-07-01

國家自然科學基金（41201588）、“十二五”國家科技支撐計劃（2012BAD14B11）、中央高?；究蒲袠I(yè)務費專項基金（Z109021538）

劉松，E-mail：ls2014050964@nwsuaf.edu.cn。

王效琴，Tel：029-87080282；E-mail：xqwang@nwsuaf.edu.cn