孫世坤，劉文艷，劉 靜，王玉寶，陳帝伊，吳普特

（1西北農(nóng)林科技大學旱區(qū)農(nóng)業(yè)水土工程教育部重點實驗室，陜西楊凌 712100；2西北農(nóng)林科技大學旱區(qū)節(jié)水農(nóng)業(yè)研究院，陜西楊凌 712100；3河海大學水文水資源學院，南京 210098）

河套灌區(qū)春小麥生產(chǎn)水足跡影響因子敏感性及貢獻率分析

孫世坤1，2，劉文艷1，劉 靜3，王玉寶1，2，陳帝伊1，吳普特1，2

（1西北農(nóng)林科技大學旱區(qū)農(nóng)業(yè)水土工程教育部重點實驗室，陜西楊凌 712100；2西北農(nóng)林科技大學旱區(qū)節(jié)水農(nóng)業(yè)研究院，陜西楊凌 712100；3河海大學水文水資源學院，南京 210098）

【目的】農(nóng)業(yè)水資源高效利用是保障國家糧食和水資源安全的重要途徑，作物用水評價是農(nóng)業(yè)用水管理的主要研究課題之一。水足跡為農(nóng)業(yè)用水評價提供了新的指標體系，研究旨在通過對作物生產(chǎn)水足跡影響因素進行定量評價，進而進行水足跡調(diào)控，以便實現(xiàn)農(nóng)業(yè)水資源高效利用。【方法】以河套灌區(qū)為研究區(qū)域，基于水足跡概念體系，利用改進的水足跡量化方法對河套灌區(qū)春小麥生產(chǎn)水足跡進行量化并分析其在研究時段內(nèi)的演變特征；利用單因素輪換（One-At-A-Time，OAT）敏感性分析方法和貢獻率分析方法探究氣候、農(nóng)業(yè)生產(chǎn)投入因子和水資源利用效率對春小麥生產(chǎn)水足跡變化的驅(qū)動力?！窘Y(jié)果】春小麥生產(chǎn)水足跡在研究時段內(nèi)呈顯著下降趨勢，從1981年的4.71 m3·kg-1，下降到2010年的1.52 m3·kg-1；同時年際變化呈現(xiàn)出較為明顯的階段性特征，可劃分為波動性下降期（1981—1987年）、快速下降期（1988—1995年）、緩慢下降期（1996—2010年），該變化規(guī)律與灌區(qū)農(nóng)業(yè)生產(chǎn)和灌溉水平的發(fā)展特征基本一致；從水足跡的藍、綠水構成來看，河套灌區(qū)春小麥生產(chǎn)水足跡中藍水足跡比例超過90%，而綠水足跡比例不足10%，這與灌區(qū)農(nóng)業(yè)生產(chǎn)用水特征相一致。敏感性分析顯示，水足跡對各個因子的敏感性差別十分顯著，日照時數(shù)、相對濕度、降水量、灌溉水利用系數(shù)和單位面積化肥用量在±20%波動的情況下，春小麥生產(chǎn)水足跡的波動范圍分別為±30%、±24%、±2%、±63%和±4%，春小麥生產(chǎn)水足跡對灌溉水利用系數(shù)、日照時數(shù)和相對濕度的敏感性較高，對其余影響因子的敏感性較低。貢獻率分析顯示，研究時段內(nèi)相對濕度的減少和降水量的增加促使春小麥生產(chǎn)水足跡的增加。而日照時數(shù)的下降、化肥使用量以及灌區(qū)用水效率的提高促使了春小麥生產(chǎn)水足跡的降低。定量分析結(jié)果顯示，化肥和灌溉水利用系數(shù)對春小麥生產(chǎn)水足跡變化的貢獻率分別為-36.89%和-39.42%，而氣候因子的綜合貢獻率僅為 2.80%，對研究時段內(nèi)春小麥生產(chǎn)水足跡下降貢獻率最大的是灌溉水利用系數(shù)，其次為單位面積化肥用量，而相對濕度、日照時數(shù)和降水量三者的貢獻率相近，貢獻率最小的是日照時數(shù)，這主要是因為日照時數(shù)在研究時段內(nèi)的變率較小?！窘Y(jié)論】氣候、農(nóng)業(yè)生產(chǎn)資料投入和水資源利用效率是影響作物生產(chǎn)水足跡的主要因素，就河套灌區(qū)而言，農(nóng)業(yè)生產(chǎn)資料投入和水資源利用效率的提高是促使灌區(qū)春小麥生產(chǎn)水足跡下降的主要因素，而氣候因子在研究時段內(nèi)對春小麥生產(chǎn)水足跡的影響較小。

水足跡；敏感性分析；貢獻率分析；河套灌區(qū)；春小麥

0　引言

【研究意義】水資源是重要的自然資源和經(jīng)濟資源，隨著人口增長與各行業(yè)用水增加，水資源供需矛盾不斷加劇。作為人口和農(nóng)業(yè)大國，中國面臨著嚴峻的水資源和糧食安全壓力。在氣候變化的背景下，全球氣溫、降水和水資源分布格局都將發(fā)生變化，從而對區(qū)域水資源狀況和農(nóng)業(yè)生產(chǎn)帶來不確定性［1］。近年來，隨著中國經(jīng)濟與社會的快速發(fā)展，水資源短缺形勢愈來愈嚴峻。農(nóng)業(yè)是中國主要的水資源消耗行業(yè)，其用水占全國用水總量的 60%左右［2］，而灌溉用水占農(nóng)業(yè)用水的近90%，灌區(qū)水資源利用效率的提高對于提高中國農(nóng)業(yè)水資源利用效率，緩解農(nóng)業(yè)用水壓力具有重要意義。同時，灌溉農(nóng)業(yè)對保障中國糧食安全發(fā)揮著重要作用，占全國耕地面積49%的灌區(qū)生產(chǎn)出約占全國總量 75%的糧食［3］。因此，合理評價灌區(qū)作物生產(chǎn)過程中水資源利用情況可為提高灌區(qū)農(nóng)業(yè)水資源管理和調(diào)控水平提供理論依據(jù)，實現(xiàn)農(nóng)業(yè)水資源的持續(xù)、高效利用。水足跡作為一個較新的水資源評價指標體系，從水資源消耗數(shù)量、類型和效率等幾個層面評價了人類生產(chǎn)、消費活動對水資源的多重影響［4］，有別于傳統(tǒng)的用水量評價指標，它是一個體現(xiàn)消耗的水量、水源類型以及污染量和污染類型的多層面指標。同時，作物生產(chǎn)水足跡的量化是進行農(nóng)產(chǎn)品虛擬水貿(mào)易分析的基礎，開展水足跡與虛擬水的研究有助于理解淡水資源的全球?qū)傩约跋M和貿(mào)易對水資源使用的影響。對作物生產(chǎn)水足跡影響因子敏感性進行研究有助于明晰作物生產(chǎn)水足跡的主要影響因素，并為水足跡調(diào)控和農(nóng)業(yè)水資源管理提供理論依據(jù)。【前人研究進展】目前，國際上對水足跡的研究大致集中在3個方面：一是水足跡理論框架及量化方法的研究；二是水足跡的案例研究；三是基于水足跡的水資源評價和管理研究。在水足跡理論及量化方法研究方面：HOEKSTRA等［4］探討了水足跡在水資源管理中的應用前景，并初步構建了水足跡量化與可持續(xù)評價框架。MEKONNEN等［5］構建了基于柵格尺度的動態(tài)水分平衡模型，并在此基礎上對全球作物的生產(chǎn)水足跡進行了量化。在水足跡案例研究方面：龍愛華等［6］、馬靜等［7］、HOEKSTRA等［8］、CAZCARRO等［9］分別以?。▍^(qū)）、國家和全球為空間分析單元，研究了區(qū)域、國家和全球的消費水足跡以及行業(yè)水足跡（旅游業(yè)）。在基于水足跡的水資源評價方面：ZHANG等［10］指出水足跡為明確人類消費活動對水資源的影響提供了新的途徑，通過進行水足跡評價，可以幫助管理機構實現(xiàn)水資源的可持續(xù)利用。ERCIN等［11］分析了未來人口、經(jīng)濟發(fā)展及生產(chǎn)和消費結(jié)構變化情景下，全球消費和生產(chǎn)水足跡的變化情況，為政府和相關機構制定水資源管理策略提供參考［12］。上述研究對于水足跡理論和應用的發(fā)展提供了有益的借鑒，但目前對于農(nóng)業(yè)生產(chǎn)水足跡的量化方法還比較單一，針對不同尺度和評價目標的量化體系尚未建立，以往對作物生產(chǎn)水足跡的量化主要通過計算田間尺度作物蒸散量求得，沒有考慮灌溉水在輸配水及田間灌水過程中的損失，因而無法反映區(qū)域尺度農(nóng)業(yè)生產(chǎn)過程中的水資源利用狀況。同時，有關水足跡的研究大多還停留在對水足跡的量化和評價方面，對作物生產(chǎn)水足跡影響因子敏感性及貢獻率的研究還較少。ZHUO等［13］評價了降水量、作物系數(shù)、生育期等對作物生產(chǎn)水足跡的影響，開始了作物生產(chǎn)水足跡影響因素的初步探討?！颈狙芯壳腥朦c】農(nóng)業(yè)生產(chǎn)是一個復雜的自然、社會交融的生產(chǎn)過程，氣候等自然因素對作物生產(chǎn)水足跡產(chǎn)生重要的影響，但“綠色革命”以來農(nóng)業(yè)生產(chǎn)投入對作物單產(chǎn)水平的提高發(fā)揮了重要作用，而作物單產(chǎn)是影響作物生產(chǎn)水足跡的重要因素，同時灌溉水利用效率水平也將直接影響作物生產(chǎn)水足跡。因此，綜合評價氣候、農(nóng)業(yè)生產(chǎn)因子和灌溉水平對作物生產(chǎn)水足跡的影響對于全面評價作物生產(chǎn)水足跡的影響因素，進而提出水足跡調(diào)控標準和措施具有重要意義。【擬解決的關鍵問題】以河套灌區(qū)為研究區(qū)域，基于水足跡概念體系，利用改進的水足跡量化方法對河套灌區(qū)春小麥生產(chǎn)水足跡進行量化并分析其在研究時段內(nèi)的演變特征；利用單因素輪換（One-At-A-Time，OAT）敏感性分析方法和貢獻率分析方法探究氣候、農(nóng)業(yè)生產(chǎn)投入因子和灌溉水平對春小麥生產(chǎn)水足跡變化的驅(qū)動力，以期為灌區(qū)水足跡調(diào)控和水資源管理提供依據(jù)。

1　數(shù)據(jù)與方法

1.1研究區(qū)概況

河套灌區(qū)位于內(nèi)蒙古巴彥淖爾市（圖1），灌區(qū)地處干旱、半干旱荒漠草原地帶，大陸性氣候特征明顯，冬季嚴寒少雪，夏季高溫干熱，年平均氣溫6—8℃，灌區(qū)降雨量少，年平均降水量為130—215 mm，年蒸發(fā)量為2 100—2 300 mm［14］。灌區(qū)有效灌溉面積57.4萬hm2，灌區(qū)主要種植作物為春小麥、玉米和向日葵等。引黃水是灌區(qū)主要的水源，河套灌區(qū)多年平均引黃水量為50億m3左右。但由于黃河年徑流量的下降，同時隨著黃河統(tǒng)一調(diào)度，河套灌區(qū)配水量日益減少，灌區(qū)水資源供需矛盾不斷加?。?5］。

圖1　河套灌區(qū)示意圖Fig.1 Location of Hetao irrigation district

1.2數(shù)據(jù)及來源

研究所需數(shù)據(jù)包括：氣象數(shù)據(jù)、農(nóng)業(yè)統(tǒng)計資料和農(nóng)業(yè)用水資料。研究所需的氣象數(shù)據(jù)包括河套灌區(qū)1981—2010年月平均最高氣溫、月平均最低氣溫、相對濕度、風速、日照時數(shù)、降水量等參數(shù)，數(shù)據(jù)來自中國氣象科學數(shù)據(jù)共享網(wǎng)和巴彥淖爾市氣象局［16］。研究區(qū)作物產(chǎn)量、播種面積、單產(chǎn)、生產(chǎn)資料投入等農(nóng)業(yè)生產(chǎn)數(shù)據(jù)來自灌區(qū)管理年報和《內(nèi)蒙古統(tǒng)計年鑒》［17］。灌區(qū)歷年引、排水、灌溉定額、灌溉水利用系數(shù)等數(shù)據(jù)來源于河套灌區(qū)管理年報、《巴彥淖爾市水資源公報》等資料［18］。

1.3作物生產(chǎn)水足跡計算方法

作物生產(chǎn)水足跡可表述為生產(chǎn)作物單位產(chǎn)量（通常為經(jīng)濟產(chǎn)量）所消耗的廣義水資源量，它由藍水足跡（作物生產(chǎn)所消耗的地表水或地下水，即灌溉水）、綠水足跡（作物生產(chǎn)過程消耗的有效降水）和灰水足跡（用于稀釋作物生產(chǎn)過程排放的污染物，使其達到環(huán)境標準所需的水資源量）組成。其中藍水足跡與綠水足跡為作物消耗性用水，而灰水足跡為改善環(huán)境性用水。由于河套灌區(qū)所面臨的主要環(huán)境問題為土壤鹽漬化，而目前國際上關于農(nóng)業(yè)灰水足跡的污染源多界定為化肥、農(nóng)藥等造成的水體污染，因此，本研究主要關注春小麥生產(chǎn)過程中的消耗性用水（即藍水足跡與綠水足跡），暫不考慮灰水足跡［4，12］。目前國際上對作物生產(chǎn)水足跡的量化主要通過計算田間尺度作物蒸散量求得，沒有考慮灌溉水在輸配水及田間灌水過程中的損失，因而無法反映區(qū)域尺度農(nóng)業(yè)生產(chǎn)過程中的水資源利用狀況。針對該方法的不足，本研究在已有研究的基礎上，根據(jù)作物生產(chǎn)水足跡定義，提出了改進的作物生產(chǎn)藍水足跡計算方法，即將灌溉水損失量中的渠道蒸發(fā)、深層滲漏后隨地下水流出研究區(qū)域外的灌溉水量納入藍水足跡量化之中，從而反映區(qū)域作物生產(chǎn)過程中的水資源利用特征，主要計算過程如下［15］：

式中，WFgreen為春小麥綠水足跡（m3·kg-1）；WFblue為春小麥藍水足跡（m3·kg-1）；Wgreen為春小麥單位面積消耗的綠水資源量（m3·hm-2）；Wblue為春小麥單位面積消耗的藍水資源量（m3·hm-2）；10為單位換算系數(shù)；ETc為春小麥生育期蒸發(fā)蒸騰量（mm）；Pe為春小麥生育期有效降水量（mm）；IRG為作物單位面積消耗的灌溉水量（m3·hm-2）；Y為作物單位面積產(chǎn)量（kg·hm-2）。

Pe根據(jù) CROPWAT模型中內(nèi)嵌的有效降水量計算公式進行量化，該公式為美國農(nóng)業(yè)部土壤保持局（USDA SCS）提出的方法［19］：

式中，Pdec為旬降水量（mm）；Pe（dec）為旬有效降水量（mm）。根據(jù)作物生育期對Pe（dec）進行求和即可得到作物整個生育期的有效降水量，沈成［20］、AUGUSTIN［21］等人的研究表明，該方法在河套灌區(qū)具有較好的適用性。

式中，WA為灌區(qū)種植業(yè)灌溉水總消耗量（m3）；α為春小麥灌溉水消耗量占灌區(qū)灌溉水消耗總量的比例；AW為春小麥的播種面積（hm2）。

WA根據(jù)灌區(qū)水量平衡進行計算，計算過程見下式：

式中，WIn為灌區(qū)引水量（m3）；WP為灌區(qū)降水補給量（m3）；WG為地下水補給量（m3）；Wout為灌區(qū)排水量（m3）；ΔW為灌區(qū)內(nèi)部蓄水的變化量（m3）；WPe為灌區(qū)種植業(yè)消耗的有效降水量（m3）；WI為工業(yè)生產(chǎn)耗水量（m3）；WL為城鄉(xiāng)生活耗水量（m3）；WE為生態(tài)環(huán)境耗水量（m3）。

α根據(jù)下式計算［15］：

1.4水足跡敏感性分析

因子敏感性是表征一個或幾個相關影響因子發(fā)生變化時，對因變量影響程度的定量參數(shù)［22-24］。論文利用OAT（One-at-a-Time）敏感性分析法分析春小麥生產(chǎn)水足跡對其影響因素的敏感性。OAT分析是一種篩選式的試驗方法，綜合了局部和全局的敏感性分析方法，通過改變準則因子的相關數(shù)據(jù)來評估因變量對自變量的敏感性［25-27］。

1.5因子貢獻率分析

根據(jù)水足跡基本理論，作物生產(chǎn)水足跡是由作物耗水量與單位面積產(chǎn)量共同決定，作物單產(chǎn)受當?shù)責崃?、水分、輻射等農(nóng)業(yè)氣候資源，以及化肥、農(nóng)藥、機械化水平等生產(chǎn)資料因子的影響，用函數(shù)表示其關系為［15］：

式中，Ci、Si分別為影響作物單位面積產(chǎn)量的氣象因子及生產(chǎn)資料投入因子。

作物生產(chǎn)過程的水資源消耗主要與區(qū)域氣候條件及灌溉設施和技術水平有關，作物耗水量可用下式表示［15］：

式中，CWU為作物單位面積水資源消耗量（m3·hm-2）；Si為作物生產(chǎn)過程中與水資源消耗相關的生產(chǎn)變量，如灌溉水利用系數(shù)。

因此，作物生產(chǎn)水足跡可用下式表示［15］：

對等式兩邊分別取對數(shù)進行線性化可得：

對上式進行多元線性回歸得到以下回歸方程：

式中，e為回歸方程常數(shù)項，其他參數(shù)與上同。

根據(jù)彈性系數(shù)的定義，利用Eviews 8.0軟件對上式進行多元線性回歸，回歸系數(shù)即為各因子的彈性系數(shù)（一定時期內(nèi)相互聯(lián)系的兩個指標變化速度的比率），農(nóng)業(yè)生產(chǎn)資料和氣象因子變化對水足跡變化的貢獻率可按下式計算［15，28-29］：

式中，δi為第 i項生產(chǎn)資料對作物生產(chǎn)水足跡變化的貢獻率；αi為第i項生產(chǎn)資料的彈性系數(shù)；ΔSi為第i項生產(chǎn)資料投入的變化量；ΔWF為作物生產(chǎn)水足跡的變化量；γi為第i項氣象因子對作物生產(chǎn)水足跡變化的貢獻率；bi為第i項氣象因子的彈性系數(shù)；ΔCi為第i項氣象因子的變化量。

2　結(jié)果

2.1春小麥生產(chǎn)水足跡年際變化分析

春小麥生產(chǎn)過程的水資源消耗有兩種類型：藍水（灌溉水）和綠水（有效降水）。對河套灌區(qū)春小麥生產(chǎn)過程中的水資源消耗分析顯示，在研究時段內(nèi)，春小麥生育期水資源消耗主要依靠灌溉水供給，而由于降水稀少，春小麥消耗的降水較少。從年代變化來看，由于區(qū)域灌溉系統(tǒng)運行水平的提高，單位面積灌溉水消耗量呈波動下降的趨勢，而降水消耗量受降水波動性的影響，其在研究時段內(nèi)也呈現(xiàn)出較為明顯的波動趨勢（圖2）。

根據(jù)灌區(qū)氣象、農(nóng)業(yè)生產(chǎn)和灌區(qū)用水數(shù)據(jù)對河套灌區(qū)1981—2010年的春小麥生產(chǎn)水足跡進行量化。利用趨勢分析和M-K檢驗［30-31］分析了春小麥生產(chǎn)水足跡的時間演變特征。圖3中生產(chǎn)水足跡在研究時段內(nèi)呈顯著下降趨勢，從1981年的4.71 m3·kg-1，下降到2010年的1.52 m3·kg-1，年均下降率為0.11 m3·kg-1。M-K趨勢檢驗值為-6.71，說明下降趨勢達到極顯著水平（P＜0.01）。同時，春小麥生產(chǎn)水足跡年際變化呈現(xiàn)出較為明顯的階段性特征：1981—1987年間主要呈現(xiàn)出波動性下降的趨勢；1988—1995年間呈現(xiàn)出快速下降的變化趨勢；1995年之后呈現(xiàn)出較為平緩的下降趨勢。從水足跡的藍、綠水構成來看，河套灌區(qū)春小麥生產(chǎn)水足跡中藍水足跡比例超過90%，而綠水足跡比例不足10%，灌區(qū)春小麥生產(chǎn)主要依靠藍水資源。從二者在研究時段內(nèi)的變化趨勢來看，春小麥生產(chǎn)藍水足跡下降趨勢較春小麥生產(chǎn)綠水足跡更為顯著。

圖2　區(qū)域尺度春小麥單位面積耗水量年際變化Fig.2 Interannual variation of water consumption of spring wheat at regional scale

圖3　春小麥生產(chǎn)水足跡年際變化Fig.3 Interannual variation of water footprint of spring wheat

春小麥生產(chǎn)水足跡是由春小麥生產(chǎn)過程中的水資源消耗量和春小麥單位面積產(chǎn)量共同決定，通過分析研究時段二者的變化規(guī)律可以發(fā)現(xiàn)（圖 4），春小麥生產(chǎn)水足跡的快速下降階段與春小麥單產(chǎn)大幅度提高（伴隨著化肥、農(nóng)藥使用量的增加）的階段基本吻合，進入1995年之后，春小麥單產(chǎn)增速放緩，其水足跡的下降主要是在單位面積耗水量下降的帶動下實現(xiàn)的。

2.2影響因子敏感性分析

鑒于河套灌區(qū)農(nóng)業(yè)生產(chǎn)主要依靠藍水資源，春小麥藍水足跡占春小麥生產(chǎn)水足跡總量的比例在90%以上，春小麥藍水足跡與春小麥總水足跡的影響因子敏感性具有大體一致的特性。因此，本研究著重分析春小麥總水足跡的因子敏感性及貢獻率。利用 OAT方法對春小麥生產(chǎn)水足跡的影響因素進行敏感性分析，結(jié)果顯示（圖 5）：對于氣象因子而言，春小麥生產(chǎn)水足跡與日照時數(shù)和降水量呈正相關關系，與相對濕度呈現(xiàn)負相關關系，即春小麥生產(chǎn)水足跡隨著日照時數(shù)和降水量的增加而增加，隨著相對濕度的增加而減小。對于農(nóng)業(yè)生產(chǎn)因子而言，春小麥生產(chǎn)水足跡與化肥使用量和灌溉水利用效率呈負相關關系，即春小麥生產(chǎn)水足跡隨著化肥使用量的增加、灌溉水利用系數(shù)的提高而降低。

圖4　春小麥單位面積產(chǎn)量及耗水量年際變化Fig.4 Interannual variation of yield and water consumption per unit area of spring wheat

氣象因子中，水足跡對日照時數(shù)的敏感性最大，日照時數(shù)在±20%的波動的情況下，春小麥生產(chǎn)水足跡的變化范圍約為±30%；相對濕度的敏感性次之，在±20%的波動的情況下，春小麥生產(chǎn)水足跡的變化范圍約為±24%；而水足跡對降水的敏感性較低。農(nóng)業(yè)生產(chǎn)水平因子中，水足跡對灌溉水利用率的敏感性較高，灌溉水利用系數(shù)在±20%變動的幅度下，春小麥生產(chǎn)水足跡的變動幅度達到±60%以上；但是春小麥生產(chǎn)水足跡對化肥使用量的敏感性較低，單位面積化肥用量在±20%的波動的情況下，春小麥生產(chǎn)水足跡的變化范圍約為±4%左右。由此可見，水足跡對各個因子的敏感性差別十分顯著，春小麥生產(chǎn)水足跡對灌溉水利用系數(shù)、日照時數(shù)和相對濕度的敏感性較高，對其余影響因子的敏感性較低。

2.3影響因子貢獻率分析

作物生產(chǎn)水足跡影響因子在研究時段內(nèi)的變化對春小麥生產(chǎn)水足跡產(chǎn)生一定的影響，為了評價影響因子變化對春小麥生產(chǎn)水足跡變化的貢獻程度，采用貢獻率分析方法定量分析了各個影響因子對春小麥生產(chǎn)水足跡變化的貢獻程度。影響因子變化率分析結(jié)果顯示：相對濕度和日照時數(shù)在研究時段內(nèi)的變率分別為-3.04%和-0.35%，而降水量、化肥使用量和灌溉水利用系數(shù)在研究時段內(nèi)的變率分別為31.40%、406.70% 和26.22%，這說明灌區(qū)相對濕度和日照時數(shù)在研究時段內(nèi)呈現(xiàn)出下降趨勢，而降水量、化肥使用量和灌溉水利用系數(shù)在研究時段內(nèi)呈增加態(tài)勢。

上述影響因子在研究時段內(nèi)的變化將會對春小麥生產(chǎn)水足跡產(chǎn)生影響，各影響因子貢獻率分析結(jié)果表明（表1），相對濕度和降水量的貢獻率為正值，日照時數(shù)、單位面積化肥用量及灌溉水利用系數(shù)的貢獻率為負值。這說明，研究時段內(nèi)相對濕度的減少和降水量的增加促使春小麥生產(chǎn)水足跡的增加。而日照時數(shù)的下降、化肥使用量以及灌區(qū)用水效率的提高促使了春小麥生產(chǎn)水足跡的降低。定量分析結(jié)果表明，對春小麥生產(chǎn)水足跡值貢獻率最大的是灌溉水利用系數(shù)，其次為單位面積化肥用量，而相對濕度、日照時數(shù)和降水量三者的貢獻率相近，貢獻率最小的是日照時數(shù)，這主要是因為日照時數(shù)在研究時段內(nèi)的變率最小。

圖5　春小麥生產(chǎn)水足跡影響因子敏感性Fig.5 Sensitivity of influencing factors of wheat water footprint

進一步分析因子貢獻率差異的原因可以發(fā)現(xiàn)，灌溉水利用系數(shù)的彈性系數(shù)和研究時段內(nèi)的變率均較大，所以其貢獻率也最大。相對濕度和日照時數(shù)的彈性系數(shù)雖然較大，但是其研究時段內(nèi)的變率較小，所以對水足跡的貢獻也較小。此外，雖然相對濕度的敏感系數(shù)為負值，但是由于年際減少的變化趨勢，其對春小麥生產(chǎn)水足跡的貢獻率為正。

表1　春小麥生產(chǎn)水足跡影響因子貢獻率分析Table 1 Contribution rate analysis of influencing factors of wheat water footprint

3　討論

氣象因子對作物生產(chǎn)水足跡的影響主要通過其對作物單位面積產(chǎn)量和耗水量的作用實現(xiàn)。關于氣象因子對作物耗水量的影響，已有學者進行了系統(tǒng)的研究。陳素英等［32］研究指出相對濕度的變化將會影響作物蒸散和氣孔導度，濕度越大將會導致大氣蒸發(fā)能力下降，作物需水量下降。同時濕度增加將會導致作物氣孔導度減小，影響作物光合作用，降低作物生物量［32］。日照時數(shù)影響著作物生育期的太陽輻射量，從而對作物生物量的積累產(chǎn)生影響，另一方面，根據(jù)水氣擴散理論，太陽輻射將為作物蒸散提供能量，因此日照時數(shù)的改變將會影響太陽輻射從而影響作物產(chǎn)量和作物耗水量。趙璐等［27］對川中丘陵區(qū)近 60年參考作物蒸發(fā)蒸騰量變化成因進行了研究，結(jié)果顯示日照時數(shù)和風速減小對參考作物蒸發(fā)蒸騰量變化的貢獻率分別為-5.403%和-2.054%；相對濕度減小和溫度增加對參考作物蒸發(fā)蒸騰量變化的貢獻率分別為 2.085%和 0.562%，研究指出日照時數(shù)和風速減小是引起該區(qū)域參考作物蒸發(fā)蒸騰量下降的主要原因。

氣候因子對作物生產(chǎn)水足跡的影響途徑如下：隨著空氣相對濕度的增加大氣蒸發(fā)能力下降，作物耗水量也隨之下降，從而促使作物生產(chǎn)水足跡下降。日照時數(shù)決定著作物生育期的太陽輻射量，日照時數(shù)的改變將會影響太陽輻射從而影響作物產(chǎn)量和作物耗水量，在二者共同作用下對作物生產(chǎn)水足跡產(chǎn)生影響。降水作為綠水的來源，只有轉(zhuǎn)化為土壤有效水并被作物利用的那部分降水構成了作物綠水足跡，由于河套灌區(qū)地處西北內(nèi)陸干旱區(qū)，夏季高溫干旱、冬季嚴寒少雪，降水稀少，農(nóng)業(yè)生產(chǎn)主要依靠灌溉，因此該區(qū)域作物生產(chǎn)水足跡受降水的影響較小。農(nóng)業(yè)生產(chǎn)資料投入及灌溉水平對作物生產(chǎn)水足跡影響途徑如下：施肥能提高土壤肥力，是提高作物單位面積產(chǎn)量的有效措施，因此化肥使用量的增加將會促使作物生產(chǎn)水足跡的下降。而灌溉水利用系數(shù)則直接決定了藍水資源的利用效率，灌溉水利用系數(shù)的提高將會使作物生育期的水資源消耗量下降，從而促使作物生產(chǎn)水足跡的下降。

河套灌區(qū)春小麥生產(chǎn)水足跡、單位面積產(chǎn)量和耗水量的年際變化結(jié)果顯示，21世紀90年代之后，由于作物單位面積產(chǎn)量年均增加幅度放緩，因而通過提高作物單產(chǎn)水平以降低作物生產(chǎn)水足跡的空間有限，同時LIU等［33］和ZHANG等［34］的研究表明由于化肥和農(nóng)藥的“邊際效益”遞減，二者使用所帶來的單產(chǎn)水平提高幅度越來越小。這一點也從研究中水足跡對化肥的敏感性較低這一結(jié)論得到佐證，而且使用過多的化肥、農(nóng)藥將會帶來農(nóng)田生態(tài)環(huán)境的惡化。因此，河套灌區(qū)作物生產(chǎn)水足跡調(diào)控的重點和發(fā)展方向是提高農(nóng)業(yè)水資源利用率和利用效率，降低作物生產(chǎn)的水資源消耗量。春小麥生產(chǎn)水足跡影響因子敏感性分析結(jié)果表明，對于河套灌區(qū)而言，春小麥生產(chǎn)水足跡對灌溉水利用效率敏感性較高，因此通過農(nóng)田水利設施和雨水積蓄利用設施建設提高灌溉水和降水利用系數(shù)，同時合理控制化肥和農(nóng)藥用量，在維持高產(chǎn)的同時，減小化肥、農(nóng)藥對區(qū)域土壤、水環(huán)境的影響，是實現(xiàn)區(qū)域農(nóng)業(yè)水資源高效、可持續(xù)利用的主要措施。

本研究利用敏感性和貢獻率分析對水足跡影響因素的敏感性和貢獻率進行了研究，并得到了一些具有參考價值的初步結(jié)論，該結(jié)論可為灌區(qū)實施水足跡調(diào)控，實現(xiàn)水資源定量化管理提供有效的途徑和方向。但本研究還存在以下局限性：數(shù)據(jù)是水足跡評價研究的基礎，水足跡量化及評價所涉及的數(shù)據(jù)面廣、量大，對數(shù)據(jù)準確性、完整性和有效性都有較高的要求，由于數(shù)據(jù)統(tǒng)計口徑和誤差，導致基礎數(shù)據(jù)本身存在一定的不確定性，同時受限于基礎數(shù)據(jù)的時間序列較短，可能導致敏感性和貢獻率分析結(jié)果存在一定的誤差。另外，本研究利用傳統(tǒng)的水量平衡進行水足跡計算，未來可借助開發(fā)適用于灌區(qū)的分布式水文模型對水足跡進行量化，可為分析農(nóng)業(yè)水文循環(huán)過程和機理提供借鑒。

4　結(jié)論

（1）在研究時段內(nèi)，河套灌區(qū)春小麥生育期水資源消耗主要依靠灌溉水供給。從年代變化來看，由于區(qū)域灌溉系統(tǒng)運行水平的提高，單位面積灌溉水消耗量呈波動下降的趨勢，而降水消耗量受降水波動性的影響，其在研究時段內(nèi)也呈現(xiàn)出較為明顯的波動趨勢。春小麥生產(chǎn)水足跡在研究時段內(nèi)呈顯著下降趨勢，年均下降率為0.11 m3·kg-1，同時其年際變化呈現(xiàn)出較為明顯的階段性特征，可劃分為波動性下降期（1981—1987年）、快速下降期（1988—1995年）、緩慢下降期（1996—2010年），這個變化規(guī)律與灌區(qū)農(nóng)業(yè)生產(chǎn)和灌溉水平的發(fā)展特征基本一致。

（2）河套灌區(qū)春小麥生產(chǎn)水足跡與日照時數(shù)和降水量呈現(xiàn)正相關關系，與相對濕度、化肥使用量和灌溉水利用系數(shù)呈負相關關系；敏感性定量分析顯示，春小麥生產(chǎn)水足跡對灌溉水利用系數(shù)最敏感，其次是日照時數(shù)、相對濕度、化肥使用量和降水量。貢獻率分析結(jié)果顯示農(nóng)業(yè)生產(chǎn)資料（表現(xiàn)為化肥的投入）和水資源利用效率（表現(xiàn)為灌溉水利用系數(shù)）的提高是促使河套灌區(qū)春小麥生產(chǎn)水足跡在研究時段內(nèi)下降的主要因素，而氣候因素由于在研究時段內(nèi)變率較小，其對春小麥生產(chǎn)水足跡的影響較小。

（3）作物生產(chǎn)水足跡對化肥的敏感性較低表明：由于“邊際效益”遞減，化肥和農(nóng)藥的使用所帶來的單產(chǎn)水平提高幅度越來越小，而且過多的使用化肥、農(nóng)藥將會帶來農(nóng)田生態(tài)環(huán)境的惡化。根據(jù)因子敏感性和貢獻率分析結(jié)果，河套灌區(qū)作物生產(chǎn)水足跡調(diào)控的方向是提高農(nóng)業(yè)水資源利用率和利用效率，降低作物生產(chǎn)的水資源消耗量。

References

［1］ IPCC. Climate Change 2014: Impacts，Adaptation，and Vulnerability Contribution of Working Group II to the Fifth Assessment Report. Cambridge: Cambridge University Press，2014.

［2］ 中華人民共和國水利部. 2010中國水資源公報. 北京: 中國水利水電出版社，2011. Ministry of Water Resources of China. China Water Resources Bulletin 2010. Beijing: China Water & Power Press，2011. （in Chinese）

［3］ 齊學斌，黃仲冬，喬冬梅，張現(xiàn)超，李平，Andersen M N. 灌區(qū)水資源合理配置研究進展. 水科學進展，2015，26（2）: 1-9. QI X B，HUANG Z D，QIAO D M，ZHANG X C，LI P，ANDERSEN M N. Research advances on the reasonable water resources allocation in irrigation district. Advances in Water Science，2015，26（2）: 1-9. （in Chinese）

［4］ HOEKSTRA A Y，CHAPAGAIN A K，ALDAYA M M，MEKONNEN M M. The Water Footprint Assessment Manual: Setting the Global Standard. London: Earthscan，2011.

［5］ MEKONNEN M M，HOEKSTRA A Y. The green，blue and grey water footprint of crops and derived crop products. Hydrology and Earth System Sciences，2011，15（5）: 1577-1600.

［6］ 龍愛華，張志強，徐中民，蘇志勇. 甘肅省水資源足跡與消費模式分析. 水科學進展，2005，16（3）: 418-425. LONG A H，ZHANG Z Q，XU Z M，SU Z Y. Analysis of Chinese cultivated land’s spatial-temporal changes and causes in recent 30 years. Advances in Water Science，2005，16（3）: 418-425. （in Chinese）

［7］ 馬靜，汪黨獻，來海亮，王茵. 中國區(qū)域水足跡的估算. 資源科學，2005，27（5）: 96-100. MA J，WANG D X，LAI H L，WANG Y. Water footprint-an application in water resources research. Resources Science，2005，27（5）: 96-100. （in Chinese）

［8］ HOEKSTRA A Y，MEKONNEN M M. The water footprint of humanity. PNAS，2012，109（9）: 3232-3237.

［9］ CAZCARRO I，HOEKSTRA A Y，SáNCHEZ CHóLIZ J. The water footprint of tourism in Spain. Tourism Management，2014，40（6）: 90-101.

［10］ ZHANG G P，HOEKSTRA A Y，MATHEWS R E. Water footprint assessment （WFA） for better water governance and sustainabledevelopment. Water Resources and Industry，2013（1/2）: 1-6.

［11］ ERCIN A E，HOEKSTRA A Y. Water footprint scenarios for 2050: A global analysis. Environment International，2014，64（1）: 71-82.

［12］ 孫世坤，王玉寶，吳普特，趙西寧. 小麥生產(chǎn)水足跡區(qū)域差異及歸因分析. 農(nóng)業(yè)工程學報，2015，31（13）: 142-148. SUN S K，WANG Y B，WU P T，ZHAO X N. Spatial variability and attribution analysis of water footprint of wheat in China. Transactions of the Chinese Society of Agricultural Engineering，2015，31（13）: 142-148. （in Chinese）

［13］ ZHUO L，MEKONNEN M M，HOEKSTRA A Y. Sensitivity and uncertainty in crop water footprint accounting: A case study for the Yellow River basin. Hydrology and Earth System Sciences，2014，18（6）: 2219-2234.

［14］ 白崗栓，張蕊，耿桂俊，任志宏，張沛琪，史吉剛. 河套灌區(qū)農(nóng)業(yè)節(jié)水技術集成研究. 水土保持通報，2010，31（1）: 149-154. BAI G S，ZHANG R，GENG G J，REN Z H，ZHANG P Q，SHI Z G. Integrating agricultural water-saving technologies in Hetao irrigation district. Bulletin of Soil and Water Conservation，2010，31（1）: 149-154. （in Chinese）

［15］ SUN S K，WU P T，WANG Y B，ZHAO X N，LIU J，ZHANG X H. The impacts of interannual climate variability and agricultural inputs on water footprint of crop production in an irrigation district of China. Science of the Total Environment，2013，444（8）: 498-507.

［16］ 國家氣象局信息中心. 中國氣象科學數(shù)據(jù)共享服務網(wǎng)［DB/OL］. ［2014-06-02］. http://cdc.cma.gov.cn/home.do，2010. National Meteorological Information Center. China Meteorological Data Sharing Service System［DB/OL］. ［2014-06-02］. http://cdc.cma. gov.cn/home.do，2010. （in Chinese）

［17］ 內(nèi)蒙古自治區(qū)統(tǒng)計局. 內(nèi)蒙古統(tǒng)計年鑒 1981-2010. 北京: 中國統(tǒng)計出版社，1982-2011. Inner Mongolia Statistical Bureau. Inner Mongolia Statistical Yearbook 1981-2010. Beijing: China Statistics Press，1982-2011. （in Chinese）

［18］ 巴彥淖爾市水務局. 巴彥淖爾市水資源公報［DB/OL］. ［2014-03-02］. http://www.htgq.gov.cn/news/show.asp？id=4379. Water Authority in Bayannur. Water Resources Bulletin in Bayannur ［DB/OL］. ［2014-03-02］. http://www.htgq.gov.cn/news/show.asp？id =4379. （in Chinese）

［19］ DOLL P，SIEBERT S. Global modeling of irrigation water requirements. Water Resources Research，2002，38（4）: 1037-1048.

［20］ 沈成. 基于 CROPWAT模型河套灌區(qū)套種模式非充分灌溉制度模擬與區(qū)域決策［D］. 內(nèi)蒙古: 內(nèi)蒙古農(nóng)業(yè)大學，2012. SHEN C. Insufficient irrigation system simulation and regional strategy under interplanting mode based on CROPWAT model in Hetao irrigation area［D］. Inner Mongolia: Inner Mongolia Agricultural University，2012. （in Chinese）

［21］ AUGUSTIN L K. 基于GIS和CROPWAT的內(nèi)蒙古河套灌區(qū)作物的最優(yōu)灌溉制度研究［D］.北京: 中國農(nóng)業(yè)大學，2015. AUGUSTIN LK. Optimal irrigation scheduling for crops: Based on GIS and CROPWAT application in Hetao district of Inner Mongolia Autonomous Region，China［D］. Beijing: China Agricultural University，2015. （in Chinese）

［22］ 黃清華，張萬昌. SWAT模型參數(shù)敏感性分析及應用. 干旱區(qū)地理，2010，33（1）: 8-15. HUANG Q H，ZHANG W C. Application and parameters sensitivity and analysis of SWAT model. Arid Land Geography，2010，33（1）: 8-15. （in Chinese）

［23］ 田雨，雷曉輝，蔣云鐘，李薇. 水文模型參數(shù)敏感性分析方法研究評述. 水文，2010，30（4）: 9-12. TIAN Y，LEI X H，JIANG Y Z，LI W. Comment on parameter sensitivity analysis of hydrological model. Journal of China Hydrology，2010，30（4）: 9-12. （in Chinese）

［24］ HAMBY D M. A review of techniques for parameter sensitivity analysis of environmental models. Environmental Monitoring and Assessment，1994，32（2）: 135-154.

［25］ SUN X Y，NEWHAM L. Three complementary methods for sensitivity analysis of a water quality model. Environmental Modelling & Software，2012，37（1）: 19-29.

［26］ 曹雯，申雙和，段春鋒. 西北地區(qū)生長季參考作物蒸散變化成因的定量分析. 地理學報，2011，66（3）: 407-415. CAO W，SHEN S H，DUAN C F. Quantification of the causes for reference crop evapotranspiration changes in growing season in Northwest China. Acta Geographica Sinica，2011，66（3）: 407-415. （in Chinese）

［27］ 趙璐，梁川，崔寧博，劉利文. 川中丘陵區(qū)參考作物蒸發(fā)蒸騰量近60年變化成因研究. 水利學報，2013，44（2）: 183-190. ZHAO L，LIANG C，CUI N B，LIU L W. Attribution analyses of ET0change in hilly area of central Sichuan in recent 60 years. Journal of Hydraulic Engineering，2013，44（2）: 183-190. （in Chinese）

［28］ 劉忠，黃峰，李保國. 2003-2011年中國糧食增產(chǎn)的貢獻因素分析.農(nóng)業(yè)工程學報，2013，29（23）:1-8. LIU Z，HUANG F，LI B G. Investigating contribution factors to China’s grain output increase in period of 2003 to 2011. Transactions of the Chinese Society of Agricultural Engineering，2013，29（23）: 1-8 （in Chinese）

［29］ 王丹. 氣候變化對中國糧食安全的影響與對策研究［D］. 武漢: 華中農(nóng)業(yè)大學，2009. WANG D. Impact of climate change on Chinese grain security and countermeasures［D］. Wuhan: Huazhong Agricultural University，2009. （in Chinese）

［30］ HAMED K H. Trend detection in hydrologic data: The Mann-Kendall trend test under the scaling hypothesis. Journal of Hydrology，2008，349（3/4）: 350-363.

［31］ 張淑杰，張玉書，隋東，蔡福，武晉雯，紀瑞鵬，陳鵬獅，劉慶婺.東北地區(qū)參考蒸散量的變化特征及其成因分析. 自然資源學報，2010，25（10）:1750-1761. ZHANG S J，ZHANG Y S，SUI D，CAI F，WU J W，JI R P，CHEN P S，LIU Q W. Changes in reference evapotranspiration and its causes in Northeast China. Journal of Natural Resources，2010，25（10）: 1750-1761. （in Chinese）

［32］ 陳素英，張喜英，邵立威，孫宏勇. 農(nóng)業(yè)技術和氣候變化對農(nóng)作物產(chǎn)量和蒸散量的影響. 中國生態(tài)農(nóng)業(yè)學報，2011，19（5）: 1039-1047. CHEN S Y，ZHANG X Y，SHAO L W，SUN H Y. Effects of climate change and agricultural technology improvement on evapotranspiration and crop yield. Chinese Journal of Eco-Agriculture，2011，19（5）: 1039-1047. （in Chinese）

［33］ LIU X J，ZHANG Y，HAN W X，TANG A H，SHEN J L，CUI Z L，VITOUSEK P，ERISMAN J W，GOULDDING K，CHRISTIE P，F(xiàn)ANGMEIER A，ZHANG F S. Enhanced nitrogen deposition over China. Nature，2013，494（7306）: 459-462.

［34］ ZHANG F S，CHEN X P，VITOUSEK P. An experiment for the world. Nature，2013，497（7447）: 33-35.

（責任編輯 楊鑫浩）

Sensitivity and Contribution Rate Analysis of the Influencing Factors of Spring Wheat Water Footprint in Hetao Irrigation District

SUN Shi-kun1，2，LIU Wen-yan1，LIU jing3，WANG Yu-bao1，2，CHEN Di-yi1，WU Pu-te1，2
（1Key Laboratory of Agricultural Soil and Water Engineering in Arid and Semiarid Areas of Ministry of Education，Northwest A&F University，Yangling 712100，Shaanxi；2Institute of Water Saving Agriculture in Arid regions of China，Northwest A&F University，Yangling 712100，Shaanxi；3College of Hydrology and Water Resources，Hohai University，Nanjing 210098）

【Objective】The efficient utilization of agricultural water resources is a key measure to guarantee national food and water security. Crop water use efficiency evaluation is one of the main research fields of agricultural water management. Water footprint provides a new index for agricultural water use evaluation，and the quantitative evaluation on the influencing factors of crop water footprint will be helpful to the implementation of water footprint control and improvement of the agricultural water use efficiency. 【Method】Based on the concept of water footprint，the water footprint of spring wheat in Hetao irrigation district was quantified by using an improved calculation method and the temporal variation was analyzed. Sensitivity and contribution rate analysis were used to quantify the relationship between crop water footprint and its influencing factors. 【Result】The results show that water footprint of wheat declined significantly during the study period. It decreased from 4.71 m3·kg-1in 1981 to 1.52 m3·kg-1in 2010. The variation of water footprint of wheat displayed an obvious stage characteristic. It can be divided into three stages: fluctuate declining stage （1981-1987），rapid declining stage （1988-1995） and slow declining stage （1996-2010）. And this variation characteristic was consistent with the variation of agricultural production and irrigation level in Hetao irrigation district. The blue water footprint accounted for the larger proportion （more than 90%），while for the share of green water footprint it was relatively small. Therefore，the production of wheat in the Hetao irrigation district mainly depended on blue water （irrigation water）. Sensitivity analysis shows that the difference of sensitivity between the influencing factors was significant. The variations of water footprint of spring wheat was ±30%，±24%，±2%，±63% and ±4% when sunshine hours，relative humidity，precipitation，irrigation water use coefficient and fertilizer rate per unit area varied at ±20%. Irrigation water use coefficient is the most sensitivity factor of wheat water footprint，following by sunshine hours，relative humidity，fertilizer usage and precipitation. Contribution analysis results show that the decline of relative humidity and the increase of precipitation led to the increase of wheat water footprint. On the contrary，the decrease of sunshine hours combined with the increase of fertilizer usage and irrigation water use efficiency led to the decrease of wheat water footprint. The contribution rates of fertilizer and irrigation water use coefficient were -36.89% and -39.42%，respectively，while the contribution rate of the total climatic factors was 2.80%. The increase of utilization coefficient of irrigation water had the largest contribution rate to the decrease of wheat water footprint during the study period，followed by fertilizer usage. The relative humidity，sunshine hours and precipitation had similar contribution rate to the variation of wheat water footprint. 【Conclusion】 The main kinds of influencing factors of crop water footprint are climate，agricultural production inputs and water use efficiency. As for Hetao irrigation district，the improvement of agricultural production and water use efficiency are the major driving forces that cause the variation of wheat water footprint in Hetao irrigation district，while the climate factors have little effect on wheat water footprint. The results of this study could provide reference for the water footprint control.

water footprint； sensitivity analysis； contribution rate analysis； Hetao irrigation district； spring wheat

2015-12-14；接受日期：2016-05-27

國家自然科學基金（51409218）、中央高?；究蒲袠I(yè)務費專項資金（2014YB050）、陜西省自然科學基礎研究計劃（2016JQ5092）

聯(lián)系方式：孫世坤，E-mail：ssksci@163.com。通信作者王玉寶，E-mail：wyb0406@sina.com