劉 勤，梅旭榮，嚴(yán)昌榮，* ，居 煇，楊建瑩

(1.中國(guó)農(nóng)業(yè)科學(xué)院農(nóng)業(yè)環(huán)境與可持續(xù)發(fā)展研究所，北京100081;2.農(nóng)業(yè)部旱作節(jié)水農(nóng)業(yè)重點(diǎn)實(shí)驗(yàn)室，北京 100081;3.農(nóng)業(yè)部農(nóng)業(yè)環(huán)境重點(diǎn)實(shí)驗(yàn)室，北京 100081)

糧食安全是當(dāng)前國(guó)內(nèi)外關(guān)注的焦點(diǎn)，尤其是近年來(lái)國(guó)內(nèi)外糧食供需矛盾越來(lái)越突出。隨著糧食主產(chǎn)區(qū)持續(xù)向缺水和生態(tài)脆弱的北方地區(qū)轉(zhuǎn)移，水土資源不相匹配、水資源匱乏對(duì)糧食生產(chǎn)的影響更加凸顯［1］。在水資源方面，我國(guó)降水呈現(xiàn)總量少、年際分布不均勻且在全球氣候變化的影響下變化愈加復(fù)雜等特性。據(jù)統(tǒng)計(jì)，1956—2005年降水總量81%分布在長(zhǎng)江流域及以南地區(qū)，而60%以上的耕地卻集中于北方地區(qū)，且以400 mm等雨量線以東地區(qū)為主。隨著我國(guó)糧食主產(chǎn)區(qū)逐漸向中部和北部轉(zhuǎn)移，單位耕地面積上水資源的占有量還將進(jìn)一步減少，可見(jiàn)，水資源不足和水土資源不相匹配成為影響我國(guó)糧食安全的主要因素。與此同時(shí)，農(nóng)業(yè)灌溉水量銳減，用水效率低，進(jìn)一步影響了我國(guó)的糧食產(chǎn)量。據(jù)統(tǒng)計(jì)，我國(guó)2/3的糧食產(chǎn)量來(lái)自于占總耕地面積1/2的灌溉面積上［2］，在農(nóng)業(yè)生產(chǎn)的規(guī)劃和實(shí)際灌溉管理中，常常通過(guò)加大灌溉力度滿足農(nóng)田用水的需求，而土壤水卻在無(wú)效消耗。據(jù)相關(guān)分析，整個(gè)華北地區(qū)降水量的55%以上轉(zhuǎn)化為土壤水資源［3］，黃淮海大部分地區(qū)的冬小麥全生育期需水量的50%—70%，夏玉米全生育期需水量的80%以上可由土壤水提供［4］，因此，研究區(qū)域尺度上降水量和作物需水量的匹配情況，摸清作物降水虧缺值變化規(guī)律，對(duì)因地制宜地采取綜合節(jié)水措施，提高農(nóng)業(yè)水資源利用效率，緩解水資源供需矛盾具有重要意義。

隨著氣候變化的加劇，農(nóng)業(yè)水資源對(duì)氣候變化尤其是全球變暖的響應(yīng)問(wèn)題，包括水循環(huán)、水量時(shí)空分布、降水極端事件與洪澇災(zāi)害等的改變逐漸受到關(guān)注［5］。作物降水虧缺是引起作物減產(chǎn)的重要原因［6］，弄清作物需水量及其空間分布，是確定作物灌溉制度以及地區(qū)灌溉用水量的基礎(chǔ)［7］。在土壤水分充分的情況下，氣象因素是影響作物需水量的主要因素［8］，同時(shí)，農(nóng)業(yè)技術(shù)措施也會(huì)對(duì)作物需水產(chǎn)生影響。多年來(lái)，對(duì)作物需水量的研究主要集中在田間或點(diǎn)的水平，主要采用經(jīng)驗(yàn)公式法、水量平衡法［9-11］和微氣象學(xué)法［12］等，并取得了較大進(jìn)展。以能量平衡原理為基礎(chǔ)的Penman公式法，只需利用常規(guī)氣象資料便可較為準(zhǔn)確地計(jì)算出參考作物的需水量，該法已成為計(jì)算參考作物需水量的一種主要方法［13-14］。我國(guó)在作物需水量的研究方面做了大量工作［15-17］，此外，有關(guān)作物系數(shù)的研究工作開(kāi)展的也比較廣泛，全國(guó)許多地方都對(duì)當(dāng)?shù)刂饕r(nóng)作物的作物系數(shù)進(jìn)行了測(cè)定，積累了比較豐富的資料。但關(guān)于區(qū)域尺度作物降水虧缺研究還相對(duì)較少，大多局限在區(qū)域全年內(nèi)水分變化或者作物全生育期內(nèi)降水盈虧方面，針對(duì)具體作物和具體生育階段降水虧缺的研究還很少見(jiàn)。楊建瑩［18］比較研究了1971—1980年和21世紀(jì)近10年兩個(gè)時(shí)期華北地區(qū)氣候資源及冬小麥生育期的變化，發(fā)現(xiàn)華北地區(qū)北部年均氣溫及≥10℃積溫增加顯著，但降水減少，暖干趨勢(shì)明顯，中部和南部年均氣溫和≥10℃積溫也呈現(xiàn)增加趨勢(shì)，但降水增多，日照下降，出現(xiàn)暖濕趨勢(shì);冬小麥生育期也發(fā)生了不同程度的變化，播期普遍推遲7—10 d，拔節(jié)期提前，其中北部地區(qū)幅度較大，約5—10 d，大部分地區(qū)冬小麥成熟期推遲5—10 d。本文旨在基于華北地區(qū)兩個(gè)不同時(shí)間階段氣象數(shù)據(jù)和冬小麥生育期資料，研究冬小麥各生育階段內(nèi)降水量與需水量匹配情況，探討氣候變化對(duì)冬小麥各生育階段降水虧缺的影響，為及早建立預(yù)警預(yù)報(bào)系統(tǒng)，合理利用農(nóng)業(yè)水資源、緩解水資源供需矛盾提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

本文采用國(guó)家“十一五”科技支撐計(jì)劃“農(nóng)田水分生產(chǎn)潛力適度開(kāi)發(fā)研究”項(xiàng)目對(duì)華北地區(qū)的劃分方法。華北地區(qū)主要指長(zhǎng)城沿線以南，淮河、秦嶺和白龍江以北，黃土高原以東，汾渭河以東地區(qū)，包括北京、天津、山西、河北、山東、河南、江蘇、安徽等省(市)的全部或部分，共計(jì)587個(gè)縣市(圖1)，位于110°12'—122°43'E，31°2'—42°42'N。據(jù)統(tǒng)計(jì)，該區(qū)耕地面積約2.44×107hm2，其中水田面積占3%，水澆地面積占54%，旱地面積占43%。年降水量在500—800 mm，季節(jié)分布不均，作物熟制以一年兩熟或兩年三熟為主，是中國(guó)小麥、玉米、棉花和花生等優(yōu)勢(shì)農(nóng)產(chǎn)品的主產(chǎn)區(qū)。

圖1 華北區(qū)主要?dú)庀笳军c(diǎn)分布圖Fig.1 Distribution of meteorological stations in Northern China

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

本文采用的氣象資料來(lái)源于國(guó)家氣象局。從華北地區(qū)及周邊篩選出107個(gè)具有1971—2009年完整觀測(cè)序列的氣象站點(diǎn)作為分析對(duì)象(華北地區(qū)內(nèi)72個(gè)氣象站點(diǎn)和華北地區(qū)周邊35個(gè)氣象站點(diǎn))。20世紀(jì)70年代冬小麥生育期數(shù)據(jù)來(lái)源于文獻(xiàn)［19］，21世紀(jì)近10年冬小麥生育期數(shù)據(jù)來(lái)源于對(duì)華北地區(qū)各縣市的冬小麥生育期調(diào)研資料。冬小麥各生育期階段的作物系數(shù)來(lái)源于文獻(xiàn)［15］。

1.3 研究方法

1.3.1 Penman-Monteith 方程

計(jì)算參考作物蒸散量的Penman-Monteith公式如下［20］

式中，ET0為參考作物蒸散量(mm/d);Rn為作物表面凈輻射量(MJ·m－2·d－1);G為土壤熱通量(MJ·m－2·d－1);γ為濕度計(jì)常數(shù)(kPa/℃);T為空氣平均溫度(℃);U2為在地面以上2m高處的風(fēng)速(m/s);es為空氣飽和水汽壓(kPa);Δ為飽和水汽壓與溫度關(guān)系曲線的斜率(kPa/℃);ea為空氣實(shí)際水壓(kPa)。

1.3.2 小麥生育期內(nèi)降水虧缺計(jì)算方法［16］

式中，Wi為生育期第i階段水分盈虧量(mm);Pei為第i階段有效降水量(mm);ETci為第i階段冬小麥需水量。

式中，ETci是第i階段冬小麥需水量(mm);kci是第i階段冬小麥的作物系數(shù)，ET0i是第i階段參考作物蒸散量(mm)。

有效降水量是根據(jù)美國(guó)農(nóng)業(yè)部土壤保持局推薦的方法，其有效性已在許多學(xué)者的研究中予以證明［21-22］。

1.3.3 敏感度分析模型

為了研究作物生育期內(nèi)降水虧缺量變化對(duì)降水量、平均溫度、日照時(shí)數(shù)、相對(duì)濕度和風(fēng)速5個(gè)基本氣象要素的變化的響應(yīng)，構(gòu)建了降水虧缺量對(duì)這5個(gè)基本要素的敏感度評(píng)價(jià)模型:

式中，Wt、Wt+1分別為基期(20世紀(jì)70年代)和末期(21世紀(jì)近10年)的降水虧缺量(mm)，xt、xt+1分別為相應(yīng)的5個(gè)基本氣象要素的值;若β＜0表明降水虧缺量對(duì)基本氣象要素變化呈反向敏感，β的絕對(duì)值越大，反向敏感性越強(qiáng);若β＞0，說(shuō)明降水虧缺量與基本氣象要素同向變化，β值越大，降水虧缺量對(duì)基本氣候要素變化的正向敏感性越高，即氣候要素的較小波動(dòng)會(huì)引起降水虧缺量的較大變化。

2 結(jié)果與分析

2.1 冬小麥全生育期降水虧缺變化

為了探討2個(gè)時(shí)間階段(20世紀(jì)70年代和21世紀(jì)近10年)冬小麥生育期內(nèi)降水虧缺變化情況，明確華北地區(qū)冬小麥生育期內(nèi)水分供需平衡變化特征，本文近10年華北地區(qū)冬小麥全生育期降水虧缺量減去70年代冬小麥全生育期降水虧缺量的差值，對(duì)其空間變化特征進(jìn)行了研究。由圖2可見(jiàn)，70年代相比，近10年來(lái)華北大部分地區(qū)冬小麥生育期內(nèi)降水虧缺愈加嚴(yán)重，研究結(jié)果與陳玉民［15］、王文峰［23］研究結(jié)論基本一致。河北省南部保定、石家莊地區(qū)，河南省鄭州以及山東省濟(jì)南地區(qū)冬小麥全生育期內(nèi)降水虧缺加重程度明顯。河南省鄭州以及山東省濟(jì)南地區(qū)70年代冬小麥全生育期內(nèi)降水虧缺量在200—300 mm，近10年來(lái)降水虧缺量增至300—350 mm，變化程度在75 mm以上。河北省唐山、秦皇島地區(qū)，山西省北部地區(qū)以及山東省膠東半島地區(qū)冬小麥全生育期內(nèi)降水虧缺程度有所緩和。70年代，這部分地區(qū)冬小麥全生育期內(nèi)降水虧缺量在300—350 mm，近10年來(lái)該地區(qū)冬小麥全生育期內(nèi)降水虧缺量下降至200—300 mm，變化幅度在50 mm以上。

圖2 冬小麥生育期內(nèi)降水虧缺變化空間分布圖Fig.2 Variation of water mismatch in whole growth stage of winter wheat

2.2 冬小麥各生育期階段降水虧缺變化

圖3 冬小麥各生育階段降水虧缺變化空間分布圖Fig.3 Variation of water mismatch of winter wheat during different periods

近10年，華北地區(qū)冬小麥播種期—返青期內(nèi)降水虧缺程度自江蘇省徐州—河南許昌一線以北都有所增加(圖3)。河北省中部及南部、山東省西北部以及山西與河南兩省的交界處冬小麥播種期—返青期內(nèi)降水虧缺變化在20 mm以上，其中以河北石家莊、衡水和山東濟(jì)南地區(qū)冬小麥播種期—返青期內(nèi)干旱程度加重，變化幅度在40 mm以上，是華北地區(qū)冬小麥播種期—返青期內(nèi)降水虧缺變化最為劇烈的地區(qū)，這部分地區(qū)需要灌溉來(lái)滿足冬小麥對(duì)水分的需求。大部分地區(qū)冬小麥返青期—拔節(jié)期內(nèi)降水虧缺程度有所緩解，只在華北中部地區(qū)有所加重(圖3)。河北省南部邯鄲、衡水地區(qū)，河南省東北部及山東省西部地區(qū)冬小麥拔節(jié)期—抽穗期內(nèi)干旱程度加劇(圖3)，增加幅度可達(dá)20 mm以上。冬小麥拔節(jié)期—抽穗期內(nèi)降水虧缺程度加重，華北北部的北京市、天津市、河北省及華北中部的山東省、河南省等省市冬小麥降水虧缺加重，降水虧缺增加量在20 mm以上，另外，70年代，江蘇省北部及安徽省北部冬小麥拔節(jié)期—抽穗期表現(xiàn)為水分盈余，而近10年來(lái)這部分地區(qū)出現(xiàn)干旱，虧缺量在20 mm以上，需要灌溉來(lái)滿足冬小麥對(duì)水分的需求。華北南部的安徽省安慶地區(qū)和江蘇省泰州地區(qū)近10年雖然也表現(xiàn)為水分盈余，但相比盈余量減少。華北中部地區(qū)冬小麥抽穗期—成熟期內(nèi)降水虧缺有所緩解，西南部地區(qū)冬小麥抽穗期—成熟期內(nèi)干旱加劇(圖3)。近10年來(lái)，北京市、天津市、河北省東部和南部、山西省、河南省北部及山東省北部冬小麥抽穗期—成熟期內(nèi)降水虧缺得到緩解，幅度可達(dá)20 mm。河南省南部及安徽省阜陽(yáng)、淮南地區(qū)冬小麥抽穗期—成熟期內(nèi)干旱加重，幅度在20 mm以上。河南省南部及安徽省南部冬小麥抽穗期—成熟期內(nèi)水分盈余面積減少，且水分盈余量減少20 mm以上。冬小麥抽穗-成熟階段是生殖生長(zhǎng)時(shí)期，是決定千粒質(zhì)量的關(guān)鍵期，在這個(gè)時(shí)段內(nèi)，河南省南部及安徽省阜陽(yáng)、淮南地區(qū)水分虧缺較重，生產(chǎn)中應(yīng)密切關(guān)注土壤墑情，做到適時(shí)灌溉，才能獲得高產(chǎn)。

2.3 冬小麥生育期內(nèi)降水虧缺量對(duì)影響因素的敏感性特征

為了查明冬小麥生育期內(nèi)降水虧缺變化原因，本文以冬小麥拔節(jié)—抽穗期為例，拔節(jié)—抽穗期是冬小麥生殖生長(zhǎng)主階段，利用構(gòu)建的敏感度分析模型，分別對(duì)華北平原內(nèi)的72個(gè)站點(diǎn)70年代和近10年冬小麥降水虧缺量對(duì)降水、平均溫度、日照時(shí)數(shù)、相對(duì)濕度以及風(fēng)速的敏感性進(jìn)行了分析計(jì)算(圖4)。

降水虧缺量對(duì)降水具有正向敏感性的氣象站點(diǎn)為53個(gè)，其中高度正向敏感(β＞0.9)的站點(diǎn)為18個(gè)，中度(0.3＜β＜0.9)13個(gè)，低度(0＜β＜0.3)22個(gè)，具有較高敏感性的站點(diǎn)主要分布在河北省的秦皇島、承德和張家口地區(qū)，河南省的鄭州、平頂山和駐馬店地區(qū)，以及安徽省的北部地區(qū)，而敏感性呈負(fù)向負(fù)向敏感的站點(diǎn)(β＜0)為19個(gè)，主要分布在山東省和主要分布在山東省以及山西省、河北省、和山西省、河南省交界地區(qū)(圖4)。降水虧缺量對(duì)平均溫度具有正向敏感性的氣象站點(diǎn)也為53個(gè)，其中高度敏感的站點(diǎn)18個(gè)，中度23個(gè)，低度12個(gè)，具有較高敏感性較高的站點(diǎn)主要分布在河北省的秦皇島、承德、張家口和衡水地區(qū)，河南省的的鄭州、平頂山和駐馬店地區(qū)，以及安徽省的北部地區(qū)，而敏感性呈負(fù)向負(fù)向敏感的站點(diǎn)有19個(gè)，主要分布在山東省以及山西省、河北省和河南省交界地區(qū)主要分布在山東省和河北省南部地區(qū)(圖4)。降水虧缺量對(duì)日照時(shí)數(shù)正向敏感性的氣象站點(diǎn)為54個(gè)，其中高度敏感的站點(diǎn)30個(gè)，中度16個(gè)，低度8個(gè)，具有較高敏感性較高的站點(diǎn)主要分布在河北省的秦皇島、承德、張家口、唐山和保定地區(qū)，安徽省的北部地區(qū)，而具有負(fù)向敏感的站點(diǎn)有18個(gè)，主要分布在河北省的邢臺(tái)和邯鄲地區(qū)和山東省和河北省的邢臺(tái)和邯鄲地區(qū)(圖4)。降水虧缺量對(duì)相對(duì)濕度具有正向敏感性呈正向的氣象站點(diǎn)為53個(gè)，其中高度敏感的站點(diǎn)17個(gè)，中度22個(gè)，低度14個(gè)，具有較高敏感性的站點(diǎn)主要分布在安徽省的北部地區(qū)和河南省中北部地區(qū)和安徽省的北部地區(qū)，而負(fù)向敏感性呈負(fù)向的站點(diǎn)有19個(gè)，主要分布在中部地區(qū)。主要分布在中部地區(qū)(圖4)。對(duì)風(fēng)速具有正向敏感性的氣象站點(diǎn)為50個(gè)，其中高度敏感的站點(diǎn)22個(gè)，中度15個(gè)，低度13個(gè)，敏感性較高的站點(diǎn)主要分布在河北省的秦皇島、承德和張家口，天津市以及安徽省的北部地區(qū)，負(fù)向敏感的站點(diǎn)有22個(gè)，主要分布在中南部地區(qū)(圖4)。

3 結(jié)論與討論

(1)70年代與近10年華北地區(qū)全生育期內(nèi)降水虧缺空間分布的總趨勢(shì)均表現(xiàn)為:由南向北降水虧缺程度逐減加重，4個(gè)不同生育階段北部地區(qū)降水虧缺嚴(yán)重，南部地區(qū)降水虧缺相對(duì)緩和的趨勢(shì)。近10年，冬小麥生育期內(nèi)降水虧缺面積較70年代略有增加。與70年代相比，近10年來(lái)華北大部分地區(qū)冬小麥生育期內(nèi)降水虧缺愈加嚴(yán)重。河北省南部保定、石家莊地區(qū)，河南省鄭州以及山東省濟(jì)南地區(qū)冬小麥全生育期內(nèi)降水虧缺加重程度明顯。

(2)近10年，華北地區(qū)冬小麥播種期—返青期內(nèi)降水虧缺程度自江蘇省徐州—河南省許昌一線往北都有所加重，返青期—拔節(jié)期內(nèi)降水虧缺程度有所緩解，只在華北中部地區(qū)有所加重，拔節(jié)期—抽穗期內(nèi)虧缺程度加重，只在華北西部山西中部太原地區(qū)水分略有盈余，抽穗期—成熟期內(nèi)降水虧缺有所緩解，西南部地區(qū)冬小麥抽穗期—成熟期內(nèi)干旱加劇，冬小麥抽穗-成熟階段生殖生長(zhǎng)時(shí)期，是決定千粒質(zhì)量的關(guān)鍵期，在這個(gè)時(shí)段內(nèi)，生產(chǎn)中應(yīng)密切關(guān)注土壤墑情，做到適時(shí)灌溉，才能獲得高產(chǎn)。

圖4 冬小麥拔節(jié)—抽穗期降水虧缺量對(duì)影響因素的敏感性格局Fig.4 Sensitivity of water deficit changes to precipitation，mean air temperature，sunshine hours，relative humidity and wind speed during jointing stage and heading stage of winter wheat

(3)冬小麥拔節(jié)期—抽穗期降水虧缺量對(duì)日照時(shí)數(shù)高度正向敏感站點(diǎn)最多，主要分布在河北省北部地區(qū)和南部安徽省，華北地區(qū)近48年太陽(yáng)輻射量減弱［24］，日照時(shí)數(shù)減少，降水虧缺程度加重，其次是風(fēng)速，降水盈虧量對(duì)溫度和降水的高度正向敏感站點(diǎn)最少。

華北區(qū)降水量南高北低，而且大部分地區(qū)呈現(xiàn)降低的趨勢(shì)［25］，而冬小麥需水量卻呈現(xiàn)北高南低的特點(diǎn)［15］。本文基于歷史氣象資料比較研究了華北區(qū)兩個(gè)時(shí)間段冬小麥全生育期和各生育期階段降水虧缺格局變化特征，將為華北區(qū)冬小麥調(diào)整合理灌溉時(shí)期和灌溉定額，達(dá)到合理灌溉提供理論指導(dǎo)。另外，在計(jì)算冬小麥生育期內(nèi)需水量時(shí)，為了選取一套完整的冬小麥作物系數(shù)值，筆者經(jīng)過(guò)思考對(duì)比，最終選用《中國(guó)主要農(nóng)作物需水量等值線圖研究》［15］一書(shū)中的冬小麥作物系數(shù)作為標(biāo)準(zhǔn)作物系數(shù)，此值基本能夠反映冬小麥不同生育階段對(duì)水分的需求特征。但是，此套數(shù)據(jù)是基于20世紀(jì)90年代初期的試驗(yàn)計(jì)算得到，至今已經(jīng)有近20年的時(shí)間，是否能夠準(zhǔn)確的表征當(dāng)前作物品種和氣候條件下作物需水特征有待進(jìn)一步研究。

［1］Yang G Y，Wang L，Wang H.Thinking of food security in China based on regional water resources and land cultivation.Transactions of the CSAE，2010，26(12):1-5.

［2］The State Council of the People's Republic of China.China's food safety plan for medium and long term(2008—2020a).2008.

［3］Liu CM.Theory and practice in hydrology water resources research.Beijing:Science Press，2004:432.

［4］Shen Z R.Su R Q.Study on countermeasures of China's agricultural water crisis.Beijing:China Agricultural Science and Technology Press，1998:235-242.

［5］Gleick P H.Climate change，hydrology and water resources，1989(03).

［6］Chen Y X，Kang SZ.Theory on deficit irrigation.Beijing:China Water Conservancy and Electricity Press，1995.

［7］Kang SZ，Cai H J.Agricultural water management.Beijing:China Agriculture Press，1996.

［8］Kang SZ，Liu X M，Xiong Y Z.Research on the model of water uptake by winter wheat root system.Acta Univ.Boreali-occidentails，1992，20(2):5-12.

［9］Feng JZ，Huang Z C.Regulation of evapotranspiration of spring wheat.Acta Agronomica Sinica，1995，21(5):544-550.

［10］Howell T A，Steiner J L，SchneiderA D，et al.Seasonal and maximum daily evaportranspiration of irrigated winter-wheat，sorghum，and cornsouthern high plains.Trans.ASAE，1997，40(3):623-634.

［11］Duan H X.Discussion on the estimating model of crop evapotranspiration.Shanxi Water Resources，2005(6):83-85

［12］Xie X Q.Experimental study for determining evapotranspiration of the field.Geographical Research，1990，9(4):94-102.

［13］Shi H B，Chen Y X，Xu Y.The application of Kriging method for ET0 optimal contour map in large area irregular sampling system.Agricultural Efficient Water Use and Water and Soil Environment Protection.Xi'an:Shaaxi Science and Technology Press，2000.

［14］Sun J S，Liu Z G，Zhang JY.Reference crop evapotranspiration calculation in windy desert area in Inner Mongolia.Irrigation and Drainage，2002，21(2):17-20.

［15］Chen Y M，Guo G S.The contour maps of major crop water requirement in China.Beijing:China Agricultural Science and Technology Press，1993.

［16］Shao X M，Yan CR.Study on the pattern of water requirements of main crops and rainfall in the Yellow River Basin.Chinese Journal of Agrometeorology，2007，28(1):40-47.

［17］Zhang SJ，Zhang Y S，Sui D，Cai F，Wu JW，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.

［18］Yang JY，Mei X R，Liu Q，Yan CR，He WQ，Liu E K，Liu S.Variations in winter wheat growth steps under climate changes in Northern China.Chinese Journal of Plant Ecology，2011，35(6):623-631.

［19］Cui D C，Liu H S，Min JR.China major agricultural climate resources altas.Beijing:China Meteorological Press，1984.

［20］Allen R G，Pereira L S，Raes D.Crop evapotranspiration guidelines for computing crop water requirements:FAOIrrigation and Drainage Paper No.56.Rome:FAO，1998.

［21］Petra Doll，Stefan Siebert.Global modeling of irrigation water requirements.Water Resources Research，2002，38(4):1-8.

［22］Li Y，Yang X G，Ye Q，Huang W H.Variation characteristics of rice water requirement in middle and lower reaches of Yangtze River during 1961—2007.Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE)，2011，27(9):175-183.

［23］Wang W F.Study on water deficit and drought risk of winter wheat in North of China.Beijing:China Agricultural University.2010.

［24］Yang JY，Liu Q，Yan C R，Mei X R.Apatial and temportal variation of solar radiation in recent 48 years in North China.Acta Ecologica Sinica，2011，31(10):2748-2756.

［25］Deng Z Y，Wang Q，Zhang Q，Qing JZ，Yang QG，Yuan ZP，Liu WJ，Xu JF.Impact of climate warming and drying on food crops in northern China and the countermeasures.Acta Ecologica Sinica，2010，30(22):6278-6288.

參考文獻(xiàn):

［1］楊貴羽，汪林，王浩.基于水土資源狀況的中國(guó)糧食安全思考.農(nóng)業(yè)工程學(xué)報(bào)，2010，26(12):1-5.

［2］中華人民共和國(guó)國(guó)務(wù)院.中國(guó)糧食安全中長(zhǎng)期規(guī)劃綱要(2008-2020年)全文.2008.

［3］劉昌明.水文水資源研究理論與實(shí)踐——?jiǎng)⒉魑倪x.北京:科學(xué)出版社，2004:432.

［4］沈振榮，蘇人瓊.中國(guó)農(nóng)業(yè)水危機(jī)對(duì)策研究.北京:中國(guó)農(nóng)業(yè)科技出版社，1998:235-242.

［6］陳亞新，康紹忠.非充分灌溉原理.中國(guó)水利水電出版社，1995.

［7］康紹忠，蔡煥杰.農(nóng)業(yè)水管理學(xué).北京:中國(guó)農(nóng)業(yè)出版社，1996.

［8］康紹忠，劉曉明，熊運(yùn)章.冬小麥根系吸水模式的研究.西北農(nóng)業(yè)大學(xué)學(xué)報(bào)，1992，20(2):5-12.

［9］馮金朝，黃子深.春小麥蒸發(fā)蒸騰的調(diào)控.作物學(xué)報(bào)，1995，21(5):544-550.

［11］段紅星.作物蒸散量計(jì)算模型探討.山西水利，2005(6):83-85.

［12］謝賢群.測(cè)定農(nóng)田蒸發(fā)的試驗(yàn)研究.地理研究，1990，9(4):94-102.

［13］史海濱，陳亞新，徐英.大區(qū)域非規(guī)則采樣系統(tǒng)ET0的最優(yōu)等值線圖Kriging法繪制應(yīng)用.農(nóng)業(yè)高效用水與水土環(huán)境保護(hù).西安:陜西科學(xué)技術(shù)出版社，2000.

［14］孫景生，劉祖貴，張寄陽(yáng).風(fēng)沙區(qū)參考作物需水量的計(jì)算.灌溉排水，2002，21(2):17-20.

［15］陳玉民，郭國(guó)雙.中國(guó)主要農(nóng)作物需水量等值線圖研究.北京:中國(guó)農(nóng)業(yè)科技出版社，1993.

［16］邵曉梅，嚴(yán)昌榮.黃河流域主要作物的降水盈虧格局分析.中國(guó)農(nóng)業(yè)氣象，2007，28(1):40-47.

［17］張淑杰，張玉書(shū)，隋東，蔡福，武晉雯，紀(jì)瑞鵬，陳鵬獅，劉慶婺.東北地區(qū)參考蒸散量的變化特征及其成因分析.自然資源學(xué)報(bào)，2010，25(10):1750-1761.

［18］楊建瑩，梅旭榮，劉勤，嚴(yán)昌榮，何文清，劉恩科，劉爽.氣候要素變化背景下華北地區(qū)冬小麥生育期變化特征研究.植物生態(tài)學(xué)報(bào)，2011，35(6):623-631.

［19］崔讀昌，劉洪順，閔謹(jǐn)如.中國(guó)主要農(nóng)作物氣候資源圖集.北京:氣象出版社，1984.

［22］李勇，楊曉光，葉清，黃晚華.1961—2007年長(zhǎng)江中下游地區(qū)水稻需水量的變化特征.農(nóng)業(yè)工程學(xué)報(bào)，2011，27(9):175-183.

［23］王文峰.北方地區(qū)冬小麥水分虧缺及干旱風(fēng)險(xiǎn)評(píng)估［D］.北京:中國(guó)農(nóng)業(yè)大學(xué).2010

［24］楊建瑩，劉勤，嚴(yán)昌榮，梅旭榮.近48年華北地區(qū)太陽(yáng)輻射量時(shí)空格局的變化特征.生態(tài)學(xué)報(bào)，2011，31(10):2748-2756

［25］鄧振鏞，王強(qiáng)，張強(qiáng)，傾繼祖，楊啟國(guó)，袁志鵬，劉文婧，徐金芳.中國(guó)北方氣候暖干化對(duì)糧食做的影響及應(yīng)對(duì)措施.生態(tài)學(xué)報(bào)，2010，30(22):6278-6288.