山西省作物生產(chǎn)藍(lán)綠水足跡核算及影響因素分析

馮變變, 劉小芳, 趙勇鋼, 李 軻

(山西師范大學(xué) 生命科學(xué)學(xué)院, 山西 臨汾 041000)

隨著人口增長(zhǎng)和社會(huì)經(jīng)濟(jì)的高速發(fā)展，我國(guó)面臨的水資源短缺壓力也愈發(fā)突出。農(nóng)業(yè)作為我國(guó)的用水大戶(hù)，據(jù)統(tǒng)計(jì)2014年用水量已占到全國(guó)總用水量63.5%[1]，合理調(diào)整農(nóng)業(yè)生產(chǎn)過(guò)程中的用水結(jié)構(gòu)、提高單位水資源的農(nóng)作物產(chǎn)量是實(shí)現(xiàn)我國(guó)水安全和糧食安全的重要途徑。水足跡理論以虛擬水理論為基礎(chǔ)，從生產(chǎn)和消費(fèi)角度揭示人類(lèi)對(duì)水資源的真實(shí)需求和占用量，為水資源評(píng)價(jià)和管理提供了新的方向[2]，近幾年已被廣泛應(yīng)用于農(nóng)業(yè)水資源管理和節(jié)水農(nóng)業(yè)研究中[3-7]。

作物生產(chǎn)水足跡(Water footprint of crop production，WFP)作為農(nóng)業(yè)用水的評(píng)價(jià)指標(biāo)不僅能夠表征農(nóng)業(yè)生產(chǎn)過(guò)程中消耗的水資源總量，而且可以反映水資源的利用類(lèi)型(綠水、藍(lán)水)[8]。藍(lán)水是指地表水或地下水，即儲(chǔ)存在淡水湖泊、河流和含水層中的水；綠水指來(lái)源于降水并通過(guò)作物蒸發(fā)蒸騰消耗掉的水[9]。目前，對(duì)作物生產(chǎn)水足跡在不同時(shí)間序列中的動(dòng)態(tài)變化以及在不同空間尺度上的分布特征已進(jìn)行了一定的研究。吳普特等[3]以省級(jí)行政區(qū)為計(jì)算單元，量化了2010年中國(guó)31個(gè)省級(jí)行政區(qū)的糧食生產(chǎn)水足跡與區(qū)域虛擬水流動(dòng)情況。操信春等[4]量化了中國(guó)灌區(qū)糧食生產(chǎn)水足跡并對(duì)其用水效率進(jìn)行了評(píng)價(jià)。史利潔等[5]量化了2000—2012年陜西省主要作物生產(chǎn)水足跡與水資源壓力指數(shù)的時(shí)空分異特征，此外還有多眾多學(xué)者就不同尺度區(qū)域作物生產(chǎn)水足跡進(jìn)行了量化[6-7]。以上研究推動(dòng)了水足跡理論的完善和進(jìn)一步發(fā)展。然而，目前的研究多集中在較大尺度區(qū)域作物生產(chǎn)水足跡的量化，對(duì)于省級(jí)以下區(qū)域尺度的研究不多，并且對(duì)于造成水足跡差異影響因素的分析還相對(duì)較少。本研究基于水足跡理論，對(duì)2005—2014年山西省11個(gè)地級(jí)市6種主要糧食作物生產(chǎn)藍(lán)水足跡和綠水足跡進(jìn)行量化，分析其時(shí)空演變特征并采用通徑分析方法對(duì)其影響因素進(jìn)行分析，以期為區(qū)域水資源的高效利用和管理以及種植結(jié)構(gòu)優(yōu)化提供一定的理論依據(jù)。

1 研究區(qū)概況與方法

1.1 研究區(qū)概況

山西省位于黃土高原東翼，地勢(shì)表現(xiàn)為東北高、西南低，地貌復(fù)雜多樣，山多川少，山地丘陵面積占全省總面積的80%以上。氣候?qū)儆跍貛Т箨懶詺夂颍媳辈町愝^大，氣溫從北向南遞增，冬季寒冷干燥，夏季炎熱多雨。全年日照時(shí)數(shù)為2 200～2 900 h，年蒸發(fā)量為1 500～2 300 mm，無(wú)霜期為120～220 d。年降雨量400～650 mm，由東南向西北遞減，干旱缺水制約著山西省農(nóng)業(yè)生產(chǎn)。

1.2 研究方法

1.2.1 作物生產(chǎn)水足跡計(jì)算 作物生產(chǎn)水足跡是指生產(chǎn)單位重量的作物所消耗的水資源，按照所消耗的水資源類(lèi)型可以分為藍(lán)水足跡和綠水足跡，本研究中作物生產(chǎn)過(guò)程中水足跡的計(jì)算按照《水足跡評(píng)價(jià)手冊(cè)》[10]所提供的方法計(jì)算，其公式如下：

式中：WFP為作物生產(chǎn)水足跡(m3/kg)；WFPg為作物生產(chǎn)綠水足跡(m3/kg)；WFPb為作物生產(chǎn)藍(lán)水足跡(m3/kg)；ETg為作物綠水蒸發(fā)蒸騰量(mm)；ETb為作物藍(lán)水蒸發(fā)蒸騰量(mm)；10為單位轉(zhuǎn)換系數(shù)；Y為作物產(chǎn)量；ETc為作物蒸發(fā)蒸騰量(mm)；ETc通過(guò)收集的氣象數(shù)據(jù)結(jié)合CROPWAT 8.0軟件求出；Pe為作物生育期的有效降雨量(mm)，可根據(jù)CROPWAT 8.0模型中內(nèi)嵌的美國(guó)農(nóng)業(yè)部土壤保護(hù)局(USDA SCS)有效降雨量公式進(jìn)行計(jì)算[11]：

(2)

式中：Pe為作物生育期的旬有效降雨量(mm)；P為作物生育期旬降雨量(mm)。

綜合作物生產(chǎn)水足跡(Water footprint of integrated crop production,IWFP)是將不同種作物生產(chǎn)水足跡對(duì)作物產(chǎn)量加權(quán)得到的綜合各作物的生產(chǎn)水足跡，能夠綜合衡量區(qū)域農(nóng)業(yè)生產(chǎn)的水資源利用狀況，計(jì)算公式如下[5]:

(3)

式中：IWFPj為區(qū)域j(為山西省的11個(gè)地級(jí)市：大同、朔州、忻州、太原、呂梁、陽(yáng)泉、晉中、長(zhǎng)治、晉城、臨汾、運(yùn)城)的綜合作物生產(chǎn)水足跡(m3/kg)；WFPi,j為地區(qū)j作物i(6種糧食作物:玉米、小麥、高粱、谷子、大豆、馬鈴薯)的生產(chǎn)水足跡(m3/kg)；Pi,j為區(qū)域j作物i的總產(chǎn)量(t)。

1.2.2 Manna-Kendall趨勢(shì)分析方法 Mann-Kendall是一種非參數(shù)統(tǒng)計(jì)檢驗(yàn)方法，用來(lái)判斷趨勢(shì)的顯著性，它無(wú)需樣本服從一定的分布，也不受少數(shù)異常值的干擾，計(jì)算公式如下[12-13]：設(shè)定時(shí)間序列Xt=(x1，x2，…，xn，n為序列長(zhǎng)度)，

(4)

s(S)=n(n-1)(2n+5)/18

(5)

式中：sgn是符號(hào)函數(shù)；當(dāng)給定顯著性水平α=5%后，若|Z|≤Z1-α/2，則接收原假設(shè)，即趨勢(shì)不顯著；若|Z|>Z1-α/2，則拒絕原假設(shè)，即趨勢(shì)顯著；在計(jì)算結(jié)果滿(mǎn)足上述二條件情況下，如果Z為正，則表明被檢測(cè)序列具有上升趨勢(shì)；若Z為負(fù)時(shí)，說(shuō)明被檢測(cè)序列趨勢(shì)是下降的；當(dāng)Z為0時(shí)，表明無(wú)顯著趨勢(shì)。

1.2.3 作物生產(chǎn)水足跡影響因素分析 各影響因素和作物生產(chǎn)水足跡關(guān)系的分析，采用SPSS 17.0軟件進(jìn)行Pearson相關(guān)分析，并在此基礎(chǔ)上采用通徑分析確定各因素對(duì)作物生產(chǎn)水足跡的影響。通徑分析能夠分析多個(gè)自變量與應(yīng)變量之間的線(xiàn)性關(guān)系，是回歸分析的拓展，可以處理較為復(fù)雜的變量關(guān)系[14]。

1.3 數(shù)據(jù)及來(lái)源

氣象數(shù)據(jù)包括2005—2014年山西省各地區(qū)逐月平均最高氣溫、月平均最低氣溫、月平均氣溫、平均風(fēng)速、相對(duì)濕度、日照時(shí)數(shù)、降雨量，數(shù)據(jù)來(lái)源于中國(guó)氣象科學(xué)數(shù)據(jù)共享網(wǎng)(http:∥cdc.cma.gov.cn/home.do)，缺失的氣象數(shù)據(jù)通過(guò)反距離權(quán)重差值方法取得。作物系數(shù)和生育期參考《北方地區(qū)主要農(nóng)作物用水灌溉定額》和聯(lián)合國(guó)糧農(nóng)組織(FAO)數(shù)據(jù)庫(kù)。11個(gè)市逐年的作物產(chǎn)量、單位面積產(chǎn)量、播種面積、化肥使用量和農(nóng)業(yè)機(jī)械總動(dòng)力等數(shù)據(jù)來(lái)源于山西省統(tǒng)計(jì)年鑒及各地區(qū)統(tǒng)計(jì)年鑒(2006—2015年)。

2 結(jié)果與分析

2.1 山西省糧食作物生產(chǎn)藍(lán)綠水足跡時(shí)間變化趨勢(shì)

2005—2014年山西省6種糧食作物生產(chǎn)消耗藍(lán)水和綠水量總體上呈現(xiàn)下降趨勢(shì)，這與作物單產(chǎn)的增加趨勢(shì)相反(圖1)。6種作物中玉米、大豆、馬鈴薯、谷子藍(lán)水足跡下降趨勢(shì)均達(dá)到顯著(M-K檢驗(yàn)值均小于-1.96)，而綠水足跡中僅有谷子下降趨勢(shì)達(dá)到顯著，且綠水足跡變異系數(shù)均小于藍(lán)水，表明作物藍(lán)水足跡下降趨勢(shì)和下降幅度均小于綠水足跡。此外，作物水足跡在作物間差異顯著，6種作物中大豆藍(lán)水足跡和綠水足跡年均值皆最大為2.74，2.11 m3/kg，玉米和小麥作為山西省最主要的糧食作物擁有較低的水足跡，生產(chǎn)單位質(zhì)量玉米所需藍(lán)水資源最小為0.42 m3/kg，小麥所需綠水資源最小為0.41 m3/kg。

山西省11個(gè)市綜合作物生產(chǎn)水足跡在研究時(shí)段內(nèi)均呈下降趨勢(shì)(圖2)。大同和呂梁的綜合作物生產(chǎn)水足跡下降幅度和波動(dòng)幅度均較大，2個(gè)區(qū)域綜合作物生產(chǎn)綠水足跡極差為0.6 m3/kg左右，其余地區(qū)極差均小于0.55 m3/kg。區(qū)域間綜合作物生產(chǎn)藍(lán)水足跡下降幅度和波動(dòng)幅度均大于綠水足跡。大同和呂梁綜合作物生產(chǎn)藍(lán)水足跡的極差為1.9，1.41 m3/kg，長(zhǎng)治和晉城綜合作物生產(chǎn)藍(lán)水足跡呈現(xiàn)相對(duì)平穩(wěn)的下降趨勢(shì)，極差均小于0.3 m3/kg，其他區(qū)域?yàn)?.5 m3/kg左右，年際波動(dòng)相對(duì)較大。區(qū)域間綜合作物整體表現(xiàn)為作物水足跡越高的地區(qū)有越大的年際波動(dòng)。

圖12005-2014年山西省糧食作物生產(chǎn)藍(lán)綠水足跡及單位面積產(chǎn)量時(shí)間演變

2.2 山西省糧食作物生產(chǎn)藍(lán)綠水足跡空間分布

山西省糧食作物生產(chǎn)藍(lán)水足跡和綠水足跡在區(qū)域間差異顯著，整體表現(xiàn)為由北向南波動(dòng)中下降的趨勢(shì)(圖3)。大同的玉米、高粱、大豆生產(chǎn)藍(lán)水足跡均為全省最高，分別為1.04，1.71，5.22 m3/kg，小麥和谷子生產(chǎn)藍(lán)水足跡為呂梁市最高，分別為1.67，2.24 m3/kg。同樣，呂梁市玉米和谷子生產(chǎn)綠水足跡也達(dá)到全省最高，分別為0.65，2.64 m3/kg。晉南地區(qū)(包括臨汾、運(yùn)城、晉城市)和晉中市作物生產(chǎn)藍(lán)水足跡和綠水足跡都比較低。此外，作物間的空間差異水平也有所不同。大豆、谷子和馬鈴薯生產(chǎn)水足跡波動(dòng)范圍較廣，谷子藍(lán)綠水足跡極差均大于1.7 m3/kg，大豆藍(lán)綠水足跡極差分別為3.68，1.97 m3/kg，空間差異顯著。馬鈴薯生產(chǎn)綠水足跡和藍(lán)水足跡極差分別為1.1，1.6 m3/kg，空間差異也較大。而小麥和玉米作為山西省最主要的糧食作物，生產(chǎn)水足跡較低，綠水足跡極差最低為0.2 m3/kg。

圖2山西省各地區(qū)綜合作物生產(chǎn)藍(lán)水和綠水足跡時(shí)間演變

圖3山西省主要糧食作物生產(chǎn)藍(lán)水和綠水足跡空間分布

2.3 作物生產(chǎn)水足跡組成

山西省農(nóng)業(yè)生產(chǎn)用水組成在區(qū)域間存在較大差異(表1)，晉南地區(qū)(包括晉城、臨汾、運(yùn)城市)糧食作物藍(lán)水消耗占比52%低于晉北地區(qū)(包括大同、朔州、忻州市)的56%，晉北地區(qū)中忻州市各作物藍(lán)水消耗比例明顯低于其他地區(qū)。作物水足跡組成在作物間差異顯著(表1)，馬鈴薯、大豆生產(chǎn)水足跡組成具有相似的空間分布，生產(chǎn)消耗均以藍(lán)水為主，占比均值為55%和56%。玉米藍(lán)綠水消耗基本持平，藍(lán)水利用率由北向南逐漸減小。谷子和高粱生產(chǎn)消耗用水除大同市藍(lán)水消耗較多外，其他地區(qū)均小于50%。小麥綠水利用率較低，水資源消耗以藍(lán)水為主。

2.4 影響因素分析

為了進(jìn)一步定性分析山西省氣候變化和農(nóng)業(yè)生產(chǎn)要素投入對(duì)于作物生產(chǎn)水足跡的影響，選取氣候因子(降雨量X1，平均風(fēng)速X2，相對(duì)濕度X3，日照時(shí)數(shù)X4，平均氣溫X5)和生產(chǎn)要素(化肥使用量X6和農(nóng)業(yè)機(jī)械總動(dòng)力X7)與綜合作物生產(chǎn)水足跡進(jìn)行相關(guān)分析。

表1 2005-2014年山西省11市糧食作物生產(chǎn)藍(lán)水足跡比例

據(jù)表2可知，氣象因子中降雨量、相對(duì)濕度、平均氣溫與綜合作物水足跡、綜合作物藍(lán)水足跡均呈負(fù)相關(guān)，且相關(guān)性達(dá)到顯著。平均風(fēng)速、日照與綜合作物生產(chǎn)水足跡呈正相關(guān)關(guān)系，且與日照達(dá)到顯著相關(guān)。而相對(duì)濕度、日照、平均氣溫與綜合作物綠水足跡也達(dá)到顯著相關(guān)。生產(chǎn)投入因子中化肥使用量、農(nóng)業(yè)機(jī)械總動(dòng)力與綜合作物總水足跡和藍(lán)水足跡呈顯著負(fù)相關(guān)。在具體作物中，氣候因子、生產(chǎn)要素投入與玉米水足跡的相關(guān)性較強(qiáng)，除風(fēng)速外其他因子均與水足跡達(dá)到顯著相關(guān)。高粱、大豆、馬鈴薯水足跡與降雨量呈顯著負(fù)相關(guān)關(guān)系。高粱水足跡和日照呈顯著正相關(guān)關(guān)系。

表2 作物生產(chǎn)水足跡與影響因素的相關(guān)分析

注：*表示顯著相關(guān)(p<0.05)，**表示極顯著相關(guān)(p<0.01)，WFP，WFPg和WFPb分別代表作物生產(chǎn)總水足跡、綠水足跡和藍(lán)水足跡。

將綜合作物生產(chǎn)水足跡與各影響因素進(jìn)行通徑分析(表3)。從各因素對(duì)于綜合作物生產(chǎn)水足跡的總體影響來(lái)看，日照對(duì)于綜合作物生產(chǎn)水足跡的影響最大，其次是平均氣溫、降水量、平均濕度、農(nóng)業(yè)機(jī)械總動(dòng)力、化肥使用量、風(fēng)速。其中風(fēng)速和日照時(shí)數(shù)對(duì)綜合作物生產(chǎn)水足跡的作用為正作用，其余因素皆為負(fù)作用，與相關(guān)分析結(jié)果相一致。降雨量和日照對(duì)綜合作物生產(chǎn)水足跡的直接影響較大。

表3 山西省作物生產(chǎn)水足跡影響因素通徑分析

注：WF為綜合作物生產(chǎn)水足跡。X1—X7分別為降雨量、平均風(fēng)速、平均濕度、日照、平均氣溫、農(nóng)用化肥總使用量、農(nóng)機(jī)總動(dòng)力。

3 討論與結(jié)論

山西省作物生產(chǎn)水足跡在研究時(shí)段內(nèi)均呈現(xiàn)下降趨勢(shì),且作物藍(lán)水足跡年際波動(dòng)和下降趨勢(shì)較綠水足跡明顯，表明灌溉水資源利用效率在逐年提高。在空間上作物生產(chǎn)水足跡呈現(xiàn)出由南向北增長(zhǎng)的趨勢(shì)，這主要與區(qū)域作物生產(chǎn)的氣候條件、用水效率、農(nóng)業(yè)生產(chǎn)技術(shù)和管理等因素的差異有關(guān)。孫世坤等[15]將中國(guó)不同區(qū)域小麥生產(chǎn)水足跡與區(qū)域降水量和作物單產(chǎn)進(jìn)行對(duì)比得出三者在空間分布上的一致性。山西省降水量南北差異顯著，從東南向西北遞減，晉南地區(qū)降雨豐富，水利設(shè)施便利，更能滿(mǎn)足作物生長(zhǎng)發(fā)育，且作物單產(chǎn)較高，其中晉南地區(qū)大豆、馬鈴薯單產(chǎn)達(dá)到晉北地區(qū)的2倍[16-17]。小麥生產(chǎn)水足跡與其他作物呈現(xiàn)不同的空間分布主要是因?yàn)闀x北地區(qū)氣候條件和土壤條件與其他地區(qū)有較大差異，不同于晉南和晉中部地區(qū)種植冬小麥，晉北部地區(qū)種植春小麥，生育期不同，降水量差異也較大，生產(chǎn)單位質(zhì)量小麥所需水分也有不同[18]。

山西省農(nóng)業(yè)生產(chǎn)用水組成在區(qū)域間和作物間差異顯著，區(qū)域間藍(lán)水比例由南向北逐漸增大，作物間小麥綠水足跡比例較低。作物生產(chǎn)綠水足跡來(lái)源于降雨，相比較藍(lán)水(灌溉水)來(lái)說(shuō)，具有較低的機(jī)會(huì)成本，且綠水在維持陸地生態(tài)系統(tǒng)和雨養(yǎng)農(nóng)業(yè)糧食安全方面具有重要意義[19]，因此分析作物水足跡組成對(duì)于明晰區(qū)域作物生產(chǎn)過(guò)程中水資源利用類(lèi)型，改善用水結(jié)構(gòu)方面具有重要意義。水足跡組成和降雨關(guān)系密切，晉南地區(qū)氣候濕潤(rùn)，降水充沛，這是晉南地區(qū)藍(lán)水足跡比例較小的主要原因。同處晉北的大同和忻州，研究時(shí)段內(nèi)大同年平均降水量為384 mm，而忻州地區(qū)年平均降水量約為433 mm，降水相對(duì)豐富是忻州地區(qū)綠水足跡比例高于大同地區(qū)的主要原因。作物對(duì)于綠水資源的利用，受當(dāng)?shù)氐慕邓?、降水時(shí)間與作物生育期吻合度影響。山西省各市降雨集中在6—8月，與玉米、高粱、谷子等作物生育期比較吻合，這是這類(lèi)作物綠水比例高于小麥的主要原因。

歸因分析結(jié)果顯示區(qū)域農(nóng)業(yè)生產(chǎn)水平和氣候因素都對(duì)作物生產(chǎn)水足跡產(chǎn)生影響，其中氣候因素中日照和降雨是造成水足跡差異的最直接因素。分析原因主要為：日照直接影響作物光合作用，從而影響作物生長(zhǎng)發(fā)育和產(chǎn)量，降雨量直接影響作物生長(zhǎng)過(guò)程中的綠水消耗[15]。日照時(shí)數(shù)通過(guò)平均氣溫對(duì)作物生產(chǎn)水足跡的間接作用較大，降水量除了直接影響外，主要通過(guò)相對(duì)濕度對(duì)作物生產(chǎn)水足跡起作用。風(fēng)速主要通過(guò)影響作物冠層的水汽擴(kuò)散速率進(jìn)而影響作物生產(chǎn)水足跡[20]。在農(nóng)業(yè)生產(chǎn)資料方面，已有研究表明化肥和農(nóng)機(jī)總動(dòng)力可大大提高農(nóng)業(yè)生產(chǎn)力，促進(jìn)作物產(chǎn)量提高[21]。

結(jié)合作物生產(chǎn)水足跡內(nèi)涵及其影響因素來(lái)看，控制和減緩作物水足跡增長(zhǎng)的途徑有：(1) 調(diào)整種植業(yè)結(jié)構(gòu)，擴(kuò)大生產(chǎn)水足跡比較優(yōu)勢(shì)區(qū)作物的播種面積，降低缺乏生產(chǎn)水足跡比較優(yōu)勢(shì)作物的播種面積，并通過(guò)虛擬水戰(zhàn)略方式進(jìn)口高耗水作物[22]；(2) 增加農(nóng)業(yè)生產(chǎn)投入，加強(qiáng)農(nóng)業(yè)機(jī)械化，提高作物單位面積產(chǎn)量；(3) 加強(qiáng)節(jié)水技術(shù)的應(yīng)用和推廣，提高降雨和灌溉水利用率。此外，化肥雖然是重要的農(nóng)業(yè)生產(chǎn)投入要素，在提高作物單產(chǎn)方面具有重要作用，但是過(guò)度使用化肥會(huì)造成土壤肥力下降和農(nóng)業(yè)面源污染[23]，同時(shí)隨著化肥和農(nóng)業(yè)機(jī)械動(dòng)力投入的增加其邊際效應(yīng)逐漸減弱，對(duì)于作物增產(chǎn)方面的作用會(huì)逐漸降低，因此，提高水資源利用效率是降低水足跡的一項(xiàng)綜合性的技術(shù)措施。

本文以地級(jí)市為計(jì)算單元研究山西省近10 a的6種糧食作物生產(chǎn)藍(lán)水和綠水足跡，尺度較為精細(xì)，時(shí)間序列較長(zhǎng)，作物種類(lèi)較全，更能綜合反映山西省農(nóng)業(yè)生產(chǎn)的用水效率，但是受限于數(shù)據(jù)的制約，本文未對(duì)灰水足跡進(jìn)行量化。此外，在對(duì)造成作物生產(chǎn)水足跡差異影響因素分析時(shí)僅考慮了氣候因子和農(nóng)業(yè)生產(chǎn)資料的影響，而未對(duì)人口、生產(chǎn)規(guī)模、農(nóng)業(yè)技術(shù)和經(jīng)濟(jì)發(fā)展等因素進(jìn)行分析，闡明這些因素對(duì)作物生產(chǎn)水足跡的影響大小及其作用機(jī)制可為降低作物生產(chǎn)水足跡提供較為全面的思路。這些都需要做進(jìn)一步的研究。

參考文獻(xiàn):

[1] 中華人民共和國(guó)水利部.中國(guó)水資源公報(bào)2014[M].北京:中國(guó)水利水電出版社,2015.

[2] Hoekstra A Y, Mekonnen M M. The water footprint of humanity [J]. Proceedings of the National Academy of Sciences of the United States of America, 2012,109(9):3232-3237.

[3] 吳普特,王玉寶,趙西寧.2010中國(guó)糧食生產(chǎn)水足跡與區(qū)域虛擬水流動(dòng)報(bào)告[M].北京:中國(guó)水利水電出版社,2012.

[4] 操信春,吳普特,王玉寶,等.中國(guó)灌區(qū)糧食生產(chǎn)水足跡及用水評(píng)價(jià)[J].自然資源學(xué)報(bào),2014,29(11):1826-1835.

[5] 史利潔,吳普特,王玉寶,等.基于作物生產(chǎn)水足跡的陜西省水資源壓力評(píng)價(jià)[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào),2015,23(5):650-658.

[6] Sun S K, Wu P T, Wang Y B, et al. Temporal variability of water footprint for maize production: The case of Beijing from 1978 to 2008 [J]. Water Resources Management, 2013,27(7):2447-2463.

[7] 蓋力強(qiáng),謝高地,李士美,等.華北平原小麥、玉米作物生產(chǎn)水足跡的研究[J].資源科學(xué),2010,32(11):2066-2071.

[8] Siebert S, D?ll P. Quantifying blue and green virtual water contents in global crop production as well as potential production losses without irrigation [J]. Journal of Hydrology, 2010,384(3/4):198-217.

[9] 戚瑞,耿涌,朱慶華.基于水足跡理論的區(qū)域水資源利用評(píng)價(jià)[J].自然資源學(xué)報(bào),2011,26(3):486-495.

[10] Hoekstra A Y, Chapagain A K, Aldaya M M, et al. The Water Footprint Assessment Manual-Setting the Global Standard[M]. London: Earthscan, 2011.

[11] D?ll P, Siebert S. Global modeling of irrigation water requirements[J]. Water Resources Research, 2002,38(4):1-8.

[12] 章誕武,叢振濤,倪廣恒.基于中國(guó)氣象資料的趨勢(shì)檢驗(yàn)方法對(duì)比分析[J].水科學(xué)進(jìn)展,2013,24(4):490-496.

[13] 張淑杰,張玉書(shū),隋東,等.東北地區(qū)參考蒸散量的變化特征及其成因分析[J].自然資源學(xué)報(bào),2010,25(10):1750-1761.

[14] 蔡甲冰,劉鈺,許迪,等.基于通徑分析原理的冬小麥缺水診斷指標(biāo)敏感性分析[J].水利學(xué)報(bào),2008,39(1):83-90.

[15] 孫世坤,王玉寶,吳普特,等.小麥生產(chǎn)水足跡區(qū)域差異及歸因分析[J].農(nóng)業(yè)工程學(xué)報(bào),2015,31(13):142-148.

[16] 趙乾坤.山西省水土保持功能分區(qū)及生態(tài)脆弱性評(píng)價(jià)[D].山東泰安:山東農(nóng)業(yè)大學(xué),2014.

[17] 朱曉磊,辛存林,盧李朋,等.山西省糧食生產(chǎn)的時(shí)空變化和驅(qū)動(dòng)因子分析[J].中國(guó)農(nóng)學(xué)通報(bào),2014,30(8):82-88.

[18] 張立生,溫輝芹,程天靈,等.山西省小麥生態(tài)區(qū)劃研究[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào),2010,18(2):410-415.

[19] Rockstr M J, Falkenmark M, Karlberg L, et al. Future water availability for global food production: the potential of green water for increasing resilience to global change [J]. Water Resources Research, 2009,45(7):142-143.

[20] 梁麗喬,李麗娟,張麗,等.松嫩平原西部生長(zhǎng)季參考作物蒸散發(fā)的敏感性分析[J].農(nóng)業(yè)工程學(xué)報(bào),2008,24(5):1-5.

[21] 馬少先.農(nóng)業(yè)機(jī)械在現(xiàn)代農(nóng)業(yè)生產(chǎn)中的地位與作用探析[J].北京農(nóng)業(yè),2012(3):172-173.

[22] Shi J, Liu J, Pinter L. Recent evolution of China′s virtual water trade: analysis of selected crops and considerations for policy [J]. Hydrology & Earth System Sciences, 2014,18(4):1349-1357.

[23] Ruttan V W. The transition to agricultural sustainability [J]. Proceedings of the National Academy of Sciences of the United States of America, 1999,96(11):5960-5967.