毛萌++任理++韓琳琳

摘要：干旱是黑龍港地區(qū)最嚴(yán)重的農(nóng)業(yè)氣象災(zāi)害，評估冬小麥各生育階段的干旱風(fēng)險可為合理制定灌溉制度、采取相應(yīng)的防旱抗旱措施提供依據(jù)?；隍炞C后的APSIM模型所模擬的1961年-2005年31個點的土壤含水量動態(tài)，計算得到了冬小麥各生育期的土壤相對濕度干旱指數(shù)，分析了農(nóng)民方案和優(yōu)化方案的該地區(qū)的農(nóng)業(yè)干旱時空分布特征。結(jié)果表明，干旱頻率發(fā)生較高的時期農(nóng)民方案中是播種越冬期和返青拔節(jié)期，分別為4889%和2760%；基于“吳橋模式”的優(yōu)化方案是越冬返青期、返青拔節(jié)期和灌漿收獲期，干旱發(fā)生的頻率為5161%、3978%和6222%。易發(fā)生干旱的地方多位于邢臺、衡水南部和邯鄲北部。在冬小麥各生育階段，干旱指數(shù)年際變化顯示出兩種方案的變化趨勢基本一致。盡管總體說來優(yōu)化方案發(fā)生干旱的風(fēng)險要高于農(nóng)民方案的，但基于“吳橋模式”的優(yōu)化方案充分發(fā)揮了根系帶的土體水庫功能，減少了灌溉水，能夠高效利用土壤水，如果全區(qū)31個點在1961年-2005年間，均按照優(yōu)化方案灌溉，相比農(nóng)民方案總共可節(jié)水61 785 mm，表明“吳橋模式”可在黑龍港地區(qū)推廣，它在一定程度上可規(guī)避干旱風(fēng)險，同時能達(dá)到節(jié)水、穩(wěn)產(chǎn)的目標(biāo)。

關(guān)鍵詞：土壤相對濕度；農(nóng)業(yè)干旱；黑龍港地區(qū)；冬小麥；吳橋模式

中圖分類號：TV122；S152.7 文獻(xiàn)標(biāo)志碼：A 文章編號：16721683（2016）06001809

Risk assessment of agricultural drought in Heilonggang District

MAO Meng，REN Li，HAN Linlin

（College of Resources and Environmental Sciences，China Agricultural University；

Key Laboratory of Arable Land Conservation（North China），Ministry of Agriculture，Beijing 100193，China）

Abstract：Drought is the most severe agricultural meteorological disaster in Heilonggang District，assessment of drought risk during different growing stages of winter wheat may provide basis to determine rational irrigation schedule and to prevent and control drought.Based on the simulated soil water dynamics of 31 profiles for 1961-2005 using the calibrated APSIM model，relative soil moisture drought indexes were calculated for different growth stages of winter wheat，and the spatial and temporal characteristics of agricultural drought of two scenarios，F(xiàn)arm Scenario（FS） and Optimized Scenario（OS） were analyzed.The results showed that，from the point of view of different growth stages of winter wheat，highfrequency drought occurred in sowingoverwintering stage and turning greenjointing stage for FS，with frequencies of 4889% and 276% respectively，and in overwinteringturning green stage，turning greenjointing stage and grain fillingharvesting stage for OS，with drought frequencies of 5161%，3978% and 6222% respectively.High drought risk frequently occurred in central region and Xingtai，south part of Hengshui and north part of Handan.The time series of relative soil moisture drought index in the different growth stages of winter wheat showed the same change trend in the time series of relative soil moisture drought index of two scenarios.This study suggested that，in general，the relative soil moisture drought index was higher in FS than OS，but the OS，which based on ′Wuqiao Mode′，brought the soil root zone reservoir into play enough，reduced the irrigation water and could use soil water with high efficiency.The total water amount of 61 785 mm could be saved in OS compared to that in FS，if irrigation schedule of the OS was adopted in 31 profiles through the district in the period of 1961-2005.It showed that ′Wuqiao Mode′ can be extended to application in Heilonggang District，and this mode may avoid drought risk in some degree，and can reach the goal of saving water and maintaining high yield.

Key words：relative soil moisture；agricultural drought；Heilonggang District；winter wheat；Wuqiao mode

干旱是一個全球性問題，人類面臨的第一個生態(tài)問題就是水分不足[1]。農(nóng)業(yè)是受干旱發(fā)生影響首當(dāng)其沖的部門，因為它依賴于作物生長不同生育期內(nèi)的水資源和土壤水的保持[2]。據(jù)民政部1949年-2003年的統(tǒng)計，中國平均每年受旱耕地面積約2 2316萬hm2，約占各種氣象災(zāi)害影響耕地面積的60%，因旱災(zāi)每年損失糧食100億kg[3]。2009年春天的黃淮海平原大旱，以及2009年秋冬至2010年春天的西南特大重度干旱給我國農(nóng)業(yè)生產(chǎn)與糧食安全構(gòu)成了嚴(yán)重威脅。華北地區(qū)近20多年來干旱不斷加劇的形勢十分嚴(yán)峻，從1970年代后期開始至今，華北的干旱不斷加劇，20世紀(jì)90年代后期以來華北地區(qū)更是連年出現(xiàn)大旱[3]?；糁螄萚4]也指出華北地區(qū)最嚴(yán)重的農(nóng)業(yè)氣象災(zāi)害是干旱，受旱和成災(zāi)面積分別占全國受旱和成災(zāi)面積的47%和51%。雖然不同學(xué)科對干旱有著不同的理解和定義，但農(nóng)業(yè)干旱總是以一定時間和空間范圍內(nèi)降水缺少，并造成農(nóng)業(yè)減產(chǎn)為特征的[5]。準(zhǔn)確地評估與監(jiān)測干旱的發(fā)生發(fā)展，是進(jìn)一步采取科學(xué)的對策以減少其不良影響的關(guān)鍵，目前國內(nèi)外研究農(nóng)業(yè)干旱的方法多是制定出各種各樣的干旱評估指標(biāo)。關(guān)于干旱指數(shù)的計算和干旱評估方法的研究，在華北平原或黃淮海平原范圍內(nèi)有過大量的報道[623]，但這些研究成果中，所采用的干旱指數(shù)多為氣象干旱指標(biāo)，如以Z指數(shù)[67]、降水距平和減產(chǎn)率[812]、標(biāo)準(zhǔn)化降水指數(shù)[1316]、Palmer干旱指數(shù)[1718]、干旱綜合指數(shù)值[1920]等居多；其他的，劉建棟等[21]建立了具有明確生物學(xué)機理的華北農(nóng)業(yè)干旱預(yù)測數(shù)值模式，袁國富和吳連海[22]使用減產(chǎn)的程度衡量作物干旱程度，康西言等[23]采用作物水分生產(chǎn)函數(shù)，嚴(yán)登華等[24]結(jié)合累積相對濕潤度指數(shù)和模糊集對評價法研究了干旱問題。Wu等[25]指出精確的土壤含水量模擬對于評價和預(yù)報農(nóng)業(yè)干旱是必要的，也是重要的，土壤含水量要比降雨重要。在采用土壤含水量作為干旱指標(biāo)方面，房穩(wěn)靜等[26]、郭晶等[27]、尹正杰等[28]、陳曉楠等[29]、Hunt等[30]、Wu等[25]在不同的地區(qū)做過一些研究。在這些以土壤含水量作為農(nóng)業(yè)干旱指標(biāo)的研究中，尚未見到對黑龍港地區(qū)內(nèi)農(nóng)業(yè)干旱的時空分析。

黑龍港地區(qū)氣候干旱少雨，水資源嚴(yán)重匱乏，是華北平原的干旱中心，目前其聞名于世的是由于深層地下水超采而形成的深層地下水大漏斗，并且這些漏斗正在逐年擴大加深，節(jié)水的任務(wù)艱巨而迫切[31]。針對水資源不足和作物高產(chǎn)之間的矛盾，中國農(nóng)業(yè)大學(xué)的研究人員們在河北省吳橋縣進(jìn)行了多年試驗研究，形成了“冬小麥節(jié)水省肥簡化高產(chǎn)四統(tǒng)一栽培技術(shù)”，在澆好底墑水和生育期間澆1～2水（每次灌溉75 mm）的條件下，產(chǎn)量達(dá)到6 000～7 500 kg/hm2，成功地實現(xiàn)了節(jié)水、高產(chǎn)、高資源效益三者的統(tǒng)一。河北省將這套節(jié)水技術(shù)體系定為“吳橋模式”，并大面積推廣應(yīng)用[3233]。本研究嘗試用土壤相對濕度干旱指數(shù)作為農(nóng)業(yè)干旱指標(biāo)，分析農(nóng)民方案和基于“吳橋模式”的優(yōu)化方案下黑龍港地區(qū)冬小麥-夏玉米輪作體系下冬小麥生育階段的干旱發(fā)生風(fēng)險及其時空變化規(guī)律，探討長期使用“吳橋模式”在該地區(qū)應(yīng)用的農(nóng)業(yè)干旱風(fēng)險，以期為準(zhǔn)確評估冬小麥各生育階段的干旱風(fēng)險和合理制定灌溉制度、有效的采取相應(yīng)的防旱抗旱預(yù)案提供科學(xué)依據(jù)。

1 材料與方法

1.1 資料來源

根據(jù)收集的來自氣象、國土、水利和農(nóng)業(yè)等部門的監(jiān)測與統(tǒng)計資料，應(yīng)用經(jīng)過驗證的農(nóng)業(yè)生產(chǎn)系統(tǒng)模型（APSIM）的模擬結(jié)果[34]來進(jìn)行分析研究，有關(guān)站點上模型率定的過程，可參考文獻(xiàn)[35]。為了驗證在區(qū)域上模型參數(shù)的可靠性，我們將在黑龍港地區(qū)所采集的31個2 m深的土壤剖面（簡記為31個點）上模擬的產(chǎn)量進(jìn)行空間插值，然后與縣域尺度的統(tǒng)計年鑒產(chǎn)量進(jìn)行了對比。在模擬結(jié)果中選出了兩種情景，一種是農(nóng)民情景，是指管理措施參照農(nóng)民的施肥量和灌溉量及其他日常管理，冬小麥的灌溉時期概化為越冬、拔節(jié)、抽穗和灌漿期。施肥管理概化為：冬小麥?zhǔn)┎デ胺屎桶喂?jié)肥，施氮肥的基追比為1∶2。該方案的1985年-2005年冬小麥產(chǎn)量模擬值與縣域上年鑒產(chǎn)量[36]的平均相對誤差為624%，證明我們對農(nóng)民方案的概化是可靠的。優(yōu)化情景是基于“吳橋模式”設(shè)置的9種節(jié)水省氮方案所選出的保證糧食產(chǎn)量、水分利用效率和氮肥偏生產(chǎn)力都較優(yōu)的方案基礎(chǔ)上進(jìn)行晚播增密設(shè)置后的方案，即在現(xiàn)有冬小麥播種日期的基礎(chǔ)上晚播10天，對冬小麥的播種密度增至750萬株/hm2。在灌溉上，冬小麥一般在播前、拔節(jié)、開花期灌溉，每次灌溉量為75 mm；施肥時間與農(nóng)民方案一致，但施肥量方面進(jìn)行了優(yōu)化。具體的，可參考文獻(xiàn)[37]中的節(jié)水省氮情景的設(shè)定。在模擬中我們考慮了冬小麥-夏玉米輪作，但在此僅分析展示冬小麥生育期的農(nóng)業(yè)干旱情況。

2 結(jié)果與分析

2.1 冬小麥生育期內(nèi)水分變化特征

圖1顯示了1961年-2005年冬小麥生育期內(nèi)的作物蒸散量（ETa）、降水量（P）和土壤相對濕度的變化動態(tài)。整個冬小麥生育期降雨量較少，45年來平均降水量為151 mm，而平均的實際蒸散量農(nóng)民方案為269 mm，優(yōu)化方案為347 mm，降水量不能滿足作物需水量。農(nóng)民方案的蒸散量模擬的偏小，可能跟早期灌溉量較少、小麥產(chǎn)量模擬的較低有關(guān)。由農(nóng)民方案和優(yōu)化方案兩種情景下0～50 cm土層的平均土壤相對濕度可以看出，兩種土壤相對濕度基本與降雨量成正相關(guān)，通常降雨較多的年份，土層的土壤相對濕度也較高；優(yōu)化情景中的土壤相對濕度明顯小于農(nóng)民情景，但兩者的變化趨勢一致。冬小麥生育期內(nèi)兩種方案的多年平均土壤相對濕度分別為83%和71%。按照表1中的標(biāo)準(zhǔn)，無論哪種情景的平均土壤相對濕度均屬于無旱等級，但實際上在冬小麥的有些生育階段有干旱的發(fā)生，所以還需要分生育期來討論。

2.2 冬小麥各生育階段干旱發(fā)生頻率變化

特征

利用公式（1）和（2）計算了在黑龍港地區(qū)31個點1961年-2005年冬小麥各生育期干旱的發(fā)生頻率，統(tǒng)計結(jié)果見表2。從表2可以看出，對于冬小麥整個生育期來說，農(nóng)民方案中發(fā)生干旱頻率較高的時期是播種越冬期，約有49%，并且特旱發(fā)生的頻率達(dá)到了10%，這是因為農(nóng)民方案中沒有灌溉播前水而是在越冬時期灌溉，從而在這一時期發(fā)生了干旱。返青到拔節(jié)時期由于耗水量大，所以有大約28%的站點發(fā)生了不同程度的干旱；其他生育期干旱發(fā)生的頻率基本在4%以下。優(yōu)化方案由于澆足了底墑，所以播種越冬期內(nèi)發(fā)生干旱的頻率較農(nóng)民方案有所下降，大約為17%，越冬返青期、返青拔節(jié)期及灌漿收獲期是干旱發(fā)生頻率較高的時期，分別為52%、40%和62%。雖然在灌漿收獲期發(fā)生了干旱，但不會影響冬小麥的產(chǎn)量（限于篇幅，在此未列出）?？傮w來說，由于優(yōu)化方案的總灌溉量遠(yuǎn)小于農(nóng)民方案的，優(yōu)化方案的干旱發(fā)生頻率要高于農(nóng)民方案的；盡管如此，優(yōu)化方案在產(chǎn)量和水分利用率方面的表現(xiàn)要優(yōu)于農(nóng)民方案，所以是可取的。

2.3 基于土壤相對濕度干旱指數(shù)的冬小麥

各生育階段干旱空間分布特征

根據(jù)公式（1）計算了31個點的1961年-2005年每年冬小麥各生育期的土壤相對濕度，并按表1的分級標(biāo)準(zhǔn)進(jìn)行了干旱等級劃分，然后對這31個點的多年平均土壤相對濕度值進(jìn)行kriging插值得到其空間分布并在區(qū)域上進(jìn)行分級顯示（圖2）。比較圖2（a）和圖2（b）可見，農(nóng)民方案在播種到越冬期有許多地方發(fā)生了干旱，是因為農(nóng)民方案沒有灌溉播前水，而我們的優(yōu)化方案，在這一時期僅有2個點發(fā)生了輕旱。在越冬返青期（圖2（c）和2（d）），由于農(nóng)民方案有越冬水，所以這個時期的土壤相對濕度要高于優(yōu)化方案的。雖然優(yōu)化方案里邢臺市周圍有7個點發(fā)生了中旱，但據(jù)吳東麗等[10]的研究，在苗期和越冬期，冬小麥處于營養(yǎng)生長階段，水分虧缺并不直接影響產(chǎn)量的形成，加之氣溫較低，土壤蒸發(fā)能力較弱，小麥耗水較少，故對水分虧缺并不十分敏感。在返青拔節(jié)期（圖2（e）和2（f）），農(nóng)民方案中有6個點發(fā)生了輕旱，但在區(qū)域上內(nèi)僅邢臺和邯鄲的小部分地區(qū)發(fā)生了輕旱；而優(yōu)化方案中有6個點發(fā)生了中旱，位于研究區(qū)內(nèi)的邢臺市和邯鄲市的縣域的土壤相對濕度都在50%～60%之間，屬于輕旱類型。在拔節(jié)期，冬小麥由營養(yǎng)生長逐漸轉(zhuǎn)入生殖生長階段，葉面積逐漸增大，加之土壤蒸發(fā)能力日益增強，小麥耗水量也日益加大，作物產(chǎn)量對水分虧缺越來越敏感[10]。兩個方案都有拔節(jié)水，故而拔節(jié)孕穗這一時期的土壤相對濕度要高于上一生育期，優(yōu)化方案雖然有一個點發(fā)生了輕旱，但整個區(qū)域的土壤相對濕度也達(dá)到75%以上（圖2（g）和2（h））。由于拔節(jié)孕穗時期是小麥生長的關(guān)鍵期，對穗數(shù)和穗粒數(shù)影響最大[32]，缺水會造成減產(chǎn)，所以要保證這一段時期的土壤墑情。但農(nóng)民方案的灌溉量明顯過大，使得整個區(qū)域的土壤相對濕度基本都高于田持的90%。在孕穗抽穗期（圖2（i）和2（j）），這一時段水分供應(yīng)較充分，只有優(yōu)化方案中有1個點發(fā)生了輕旱，兩種方案在這個研究區(qū)上的土壤相對濕度都大于70%，但土壤相對濕度小于上一個生育期的。在抽穗灌漿期（圖2（k）和2（l）），由于農(nóng)民方案這一時期有灌溉，所以土壤貯水量繼續(xù)增高，優(yōu)化方案由于上一時期的灌溉水分繼續(xù)向下運動，這一時期0～50 cm土體的土壤貯水量較上一時期有所下降，但整個研究區(qū)的土壤相對濕度也都大于70%，這說明農(nóng)民方案中的抽穗水可以省略。吳東麗等[10]指出，在抽穗開花期，冬小麥葉面積最大，作物產(chǎn)量對水分虧缺最為敏感；而農(nóng)民方案和優(yōu)化方案都保證了這一時期的充分供水。在灌漿收獲期（圖2（m）和2（n）），農(nóng)民方案由于有灌溉，0～50 cm土層中的土壤相對濕度大于90%，沒有干旱現(xiàn)象的發(fā)生；而優(yōu)化方案的土壤含水量較上一時期有明顯的下降。農(nóng)民方案的土壤相對濕度仍然處于超過田持90%的狀態(tài)，優(yōu)化方案中有多于1/2的點出現(xiàn)了干旱，其中有8個點還發(fā)生了中旱，整個研究區(qū)大部分地區(qū)的土壤相對濕度處于50%～60%，屬于輕旱。然而在灌漿成熟期，冬小麥籽粒趨于飽滿，水分虧缺對產(chǎn)量形成的影響漸小[10]，所以我們的優(yōu)化方案在這個時期雖然出現(xiàn)了干旱，但對產(chǎn)量的影響不大。

2.4 冬小麥各生育階段土壤相對濕度干旱

指數(shù)的年際變化特征

將整個研究區(qū)31個點每年的冬小麥生長的7個階段的土壤相對濕度取平均后，分析冬小麥各生育時期的土壤相對濕度干旱指數(shù)在1961年-2005年間的年際變化。由圖3可以看出，總的來說，農(nóng)民方案的土壤相對濕度的年際變化較優(yōu)化方案的年際變化要平緩，這可能是由于農(nóng)民方案供水較充足，在拔節(jié)期以后土壤相對濕度一般在80%～100%間波動，農(nóng)民方案中，各生育期年際間變異較大的為播種越冬期、返青拔節(jié)期和孕穗抽穗期；農(nóng)民方案由于灌溉量大，反映不出枯水年年份，而優(yōu)化方案中土壤相對濕度較低的年份一般對應(yīng)的是枯水年和特枯水年。從5年移動平均線可以看出，兩種方案在各生育期內(nèi)的平均土壤相對濕度變化趨勢基本一致。在播種越冬期，優(yōu)化方案的土壤相對濕度高于農(nóng)民方案的，兩者都是在1967年前呈下降趨勢，1967年-1978年間呈緩慢增加趨勢，1978年-1983年呈急速下降趨勢，1983年以后波動平緩。在越冬返青期，農(nóng)民方案由于有越冬水，土壤相對濕度的年際間變化不明顯，而優(yōu)化方案的土壤相對濕度在1967年前呈下降趨勢，1967年-1978年有緩慢增加趨勢，1978年-1985年下降幅度較大，1985年-2005年間呈波動狀態(tài)。返青拔節(jié)期，兩種方案的土壤相對濕度都是在1990年和1997年最高，而在1999年和2002年最低。在拔節(jié)期以后，兩種方案都是在1961年-1964年呈上升趨勢、1964年-1968年呈明顯的下降趨勢，1968年以后成波動狀態(tài)，但總體上有略微下降趨勢。綜合這幾個時期土壤相對濕度的變化情況，在拔節(jié)期之前，1978年之前為年際間變化相對平緩的時期，1978年-1985年有一明顯的下降趨勢，1985年以后變化也相對平緩，但波動幅度比1961年-1978年間的幅度大。拔節(jié)期以后，1961年-2005年間兩種方案均呈緩慢下降趨勢，1998年-2003年間的下降幅度增大，但2003年-2005年稍有所回升。以干旱較嚴(yán)重的位于邢臺市廣宗縣的站點為代表，在1961年-2005年間，優(yōu)化方案相比農(nóng)民方案可節(jié)水2 010 mm，這些水如果按照優(yōu)化灌溉方案，大約可提供該站點9季冬小麥的生長。如果全區(qū)31個點在1961年-2005年間，均按照優(yōu)化方案灌溉，相比農(nóng)民方案總共可節(jié)水61 785 mm。

3 結(jié)論

基于模型模擬的1961年-2005年黑龍港地區(qū)31個點的土壤含水量動態(tài)，計算得到了冬小麥各生育期的土壤相對濕度干旱指數(shù)，分析了兩種方案下該地區(qū)農(nóng)業(yè)干旱的時空分布特征。結(jié)果表明，（1）對于冬小麥整個生長發(fā)育階段來說，農(nóng)民方案中發(fā)生干旱頻率較高的時期是播種越冬期和返青拔節(jié)期，優(yōu)化方案中越冬返青期、返青拔節(jié)期和灌漿收獲期是干旱發(fā)生頻率較高的時期；干旱均以輕旱為主，中旱次之。雖然優(yōu)化方案發(fā)生干旱的風(fēng)險高于農(nóng)民方案的，但并不影響優(yōu)化方案的產(chǎn)量和水分利用率，還是值得推廣的。（2）從兩種方案的土壤相對濕度干旱指數(shù)的空間分布圖可知，對于冬小麥各生育期，黑龍港地區(qū)的西部地區(qū)的土壤相對濕度干旱指數(shù)要小于東部，說明西部發(fā)生干旱的風(fēng)險要大于東部，易發(fā)生干旱的地方多位于邢臺、衡水南部和邯鄲北部。（3）從冬小麥不同生育階段土壤相對濕度干旱指數(shù)的年際變化特征看，在冬小麥各生育階段，優(yōu)化方案和農(nóng)民方案的變化趨勢基本一致，優(yōu)化方案的土壤相對濕度干旱指數(shù)雖然小于農(nóng)民方案的，但優(yōu)化方案45年間在全區(qū)31個點可節(jié)水61 785 mm。

當(dāng)然，此研究也有一些局限性：由于收集到的評價標(biāo)準(zhǔn)、資料中未涉及到其他層次的土壤相對濕度對應(yīng)的農(nóng)業(yè)干旱的評判標(biāo)準(zhǔn)，所以只選擇了較淺的50 cm土層用于農(nóng)業(yè)干旱的評價。受收集的資料限制，在模擬時簡單的將整個黑龍港地區(qū)分為南北兩個分區(qū)，以保定和滄州為界，北部的滄州、保定、廊坊市的縣（市）采用滄州吳橋的農(nóng)民方案、作物品種，南部的衡水、邢臺和邯鄲的縣（市）采用邯鄲曲周的農(nóng)民方案、作物品種；并且在45年的模擬時段里，由于缺乏比較詳細(xì)的品種變化資料，沒有考慮實際生產(chǎn)中的品種更替；在模擬時，不能滿足每一個模擬站點所在縣域有相對應(yīng)的氣象站點，基本上是每個地級市一個氣象站，這些都會對模擬結(jié)果有影響?？傮w說來優(yōu)化方案發(fā)生干旱的風(fēng)險要高于農(nóng)民方案的，但基于“吳橋模式”的優(yōu)化方案充分發(fā)揮了根系帶的土體水庫功能，減少了灌溉水，能夠高效利用土壤水，并創(chuàng)造前期和后期上層土壤適度水分虧缺的環(huán)境，從而可減少氮素?fù)p失。優(yōu)化方案在黑龍港地區(qū)的應(yīng)用表明“吳橋模式”可在黑龍港地區(qū)推廣，它在一定程度上可規(guī)避干旱風(fēng)險，同時能達(dá)到節(jié)水、穩(wěn)產(chǎn)的目標(biāo)。參考文獻(xiàn)（References）：

[1] 山侖.科學(xué)應(yīng)對農(nóng)業(yè)干旱[J].干旱地區(qū)農(nóng)業(yè)研究，2011，19（2）：15.（SHAN Lun.To cope rationally with agricultural drought[J].Agricultural Research in the Arid Areas，2011，19（2）：15.（in Chinese））

[2] Narasimhan B，Srinivasan R.Development and evaluation of Soil Moisture Deficit Index（SMDI）and Evapotranspiration Deficit Index（ETDI）for agricultural drought monitoring[J].Agricultural and Forest Meteorology，2005，133：6988.

[3] 秦大河.氣候變化與干旱[J].科技導(dǎo)報，2009，27（11）：3.（QIN Dahe.Climate change and drought[J].Science & Technology Review，2009，27（11）：3.（in Chinese））

[4] 霍治國，李世奎，白月明，等.冬小麥中輕度水分脅迫的增產(chǎn)節(jié)水效應(yīng)研究[J].自然資源學(xué)報，2003，18（1）：5866.（HUO Zhiguo，LI Shikui，BAI Yueming，et al.The experimental research on slight and moderate water stress effects on increasing yield and saving water mechanism of winter wheat[J].Journal of Natural Resources，2003，18（1）：5866.（in Chinese））

[5] 朱自璽，劉榮花，方文松，等.華北地區(qū)冬小麥干旱評估指標(biāo)研究[J].自然災(zāi)害學(xué)報，2003，12（1）：145150.（ZHU Zixi，LIU Ronghua，F(xiàn)ANG Wensong，et al.Evaluation indices of drought of winter wheat in North China[J].Journal of Natural Disasters，2003，12（1）：145150.（in Chinese））

[6] 邵曉梅，劉勁松，許月卿.河北省旱澇指標(biāo)的確定及其時空分布特征研究[J].自然災(zāi)害學(xué)報，2001，10（4）：133136.（SHAO Xiaomei，LIU Jinsong，XU Yueqing.Determination of droughtflood index and its temporal and spatial distribution in Hebei Province[J].Journal of Natural Disasters，2001，10（4）：133136.（in Chinese））

[7] 阮新，劉學(xué)鋒，李元華.河北省近40年干旱變化特征分析.干旱地區(qū)資源與環(huán)境，2008，22（1）：5053.（RUAN Xin，LIU Xuefeng，LI Yuanhua.Drought climate changes in Hebei province during 19612000[J].Journal of Arid Land Resources and Environment，2008，22（1）：5053.（in Chinese））

[8] 王石立，婁秀榮.華北地區(qū)冬小麥干旱風(fēng)險評估的初步研究[J].自然災(zāi)害學(xué)報，1997，6（3）：6368.（WANG Shili，LOU Xiurong.Preliminary study on the drought risk assessment of winter wheat in North China[J].Journal of Natural Disasters，1997，6（3）：6368.（in Chinese））

[9] 劉榮花，朱自璽，方文松，等.華北平原冬小麥干旱區(qū)劃初探[J].自然災(zāi)害學(xué)報，2003，12（1）：140144.（LIU Ronghua，ZHU Zixi，F(xiàn)ANG Wensong，et al.Preliminary study on drought zonation of winter wheat in North China Plain[J].Journal of Natural Disasters，2003，12（1）：140144.（in Chinese））

[10] 吳東麗，王春乙，張雪芬，等.華北冬小麥作物氣候干旱指數(shù)研究[J].科技導(dǎo)報，2009，27（7）：3236.（WU Dongli，WANG Chunyi，ZHANG Xuefen，et al.Study on the crop climate drought index of winter wheat in North China[J].Science & Technology Review，2009，27（7）：3236.（in Chinese））

[11] 李樹巖，劉榮花，馬志紅.基于降水距平的黃淮平原夏玉米干旱評估指標(biāo)研究[J].干旱地區(qū)農(nóng)業(yè)研究，2012，30（3）：253256.（LI Shuyan，LIU Ronghua，MA Zhihong.Research on drought assessment index of summer maize base on precipitation anomaly[J].Agricultural Research in the Arid Areas，2012，30（3）：253256.（in Chinese））

[12] 康西言，李春強，代立芹.河北省冬小麥生產(chǎn)干旱風(fēng)險分析[J].干旱地區(qū)農(nóng)業(yè)研究，2012，30（6）：232237.（KANG Xiyan，LI Chunqiang，DAI Liqin.Risk assessment of drought for winter wheat production in Hebei Province[J].Agricultural Research in the Arid Areas，2012，30（6）：232237.（in Chinese））

[13] 張強.華北地區(qū)干旱指數(shù)的確定及其應(yīng)用[J].災(zāi)害學(xué)，1998，13（4）：3438.（ZHANG Qiang.Research on determination of drought index in North China and its application[J].Journal of Catastrophology，1998，13（4）：3438.（in Chinese））

[14] 車少靜，李春強，申雙和.基于SPI的近41年（19652005）河北省旱澇時空特征分析[J].中國農(nóng)業(yè)氣象，2010，31（1）：137143.（CHE Shaojing，LI Chunqiang，SHEN Shuanghe.Analysis of droughtflood spatialtemporal characteristics based on standard precipitation index（SPI）in Hebei Province during 19652005[J].Chinese Journal of Agrometeorology，2010，31（1）：137143.（in Chinese））

[15] 閆峰，王艷姣，吳波.近50年河北省干旱時空分布特征[J].地理研究，2010，29（3）：423430.（YAN Feng，WANG Yanjiao，WU Bo.Spatial and temporal distributions of drought in Hebei Province over the past 50 years[J].Geographical Research，2010，29（3）：423430.（in Chinese））

[16] 趙林，武建軍，呂愛鋒，等.黃淮海平原及其附近地區(qū)干旱時空動態(tài)格局分析—基于標(biāo)準(zhǔn)化降雨指數(shù)[J].資源科學(xué)，2011，33（3）：468476.（ZHAO Lin，WU Jianjun，LV Aifeng，et al.Spatial and temporal analysis of drought over the HuangHuaiHai Plain and its surroundings based on the standardized precipitation index[J].Resources Science，2011，33（3）：468476.（in Chinese））

[17] 衛(wèi)捷，陶詩言，張慶云.Palmer干旱指數(shù)在華北干旱分析中的應(yīng)用[J].地理學(xué)報，2003，58（增刊）：9199.（WEI Jie，TAO Shiyan，ZHANG Qingyun.Analysis of drought in northern China based on the Palmer severity drought index[J].Acta Geographica Sinica 2003，58（sup.）：9199.（in Chinese））

[18] 盧洪健，莫興國，胡實.華北平原19602009 年氣象干旱的時空變化特征[J].自然災(zāi)害學(xué)報，2012，21（6）：7282.（LU Hongjian，MO Xingguo，HU Shi.Spatiotemporal variation characteristics of meteorological droughts in North China Plain during 19602009[J].Journal of Natural Disasters，2012，21（6）：7282.（in Chinese））

[19] 張婧，梁樹柏，許曉光，等.基于CI指數(shù)的河北省近50 年干旱時空分布特征[J].資源科學(xué)，2012，34（6）：10891094.（ZHANG Jing，LIANG Shubo，XU Xiaoguang，et al.Temporal and spatial distribution characteristics of droughts for recent 50 years in Hebei province based on meteorological drought composite index[J].Resources Science，2012，34（6）：10891094.（in Chinese））

[20] 張金龍，劉學(xué)鋒，于長文.河北省干旱分布特征和變化規(guī)律分析[J].干旱區(qū)研究，2012，29（1）：4146.（ZHANG Jinlong，LIU Xuefeng，YU Changwen.Analysis on distribution and variation of drought in Hebei province[J].Arid Zone Research，2012，29（1）：4146.（in Chinese））

[21] 劉建棟，王馥棠，于強，等.華北地區(qū)農(nóng)業(yè)干旱預(yù)測模型及其應(yīng)用研究[J].應(yīng)用氣象學(xué)報，2003，14（5）：593604.（LIU Jiandong，WANG Futang，YU Qiang，et al.A drought prediction model in North China Plain and its application[J].Journal of Applied Meteorological Science，2003，14（5）：593604.（in Chinese））

[22] 袁國富，吳連海.京津冀平原地區(qū)一年兩熟條件下農(nóng)業(yè)干旱時空分析[J].中國農(nóng)業(yè)氣象，2000，21（2）：59，30.（YUAN Guofu，WU Lianhai.Spatial and temporal analysis on agricultural drought under double cropping system in BeijingTianjinHebei Plain area[J].Chinese Journal of Agrometeorology，2000，21（2）：59，30.（in Chinese））

[23] 康西言，李春強，馬輝杰，等.基于作物水分生產(chǎn)函數(shù)的冬小麥干旱評估模型[J].中國農(nóng)學(xué)通報，2011，27（8）：274279.（KANG Xiyan，LI Chunqiang，MA Huijie，et al.Study on drought evaluation model of winter wheat based on crop water production function[J].Chinese Agricultural Science Bulletin，2011，27（8）：274279.（in Chinese））

[24] 嚴(yán)登華，袁喆，楊志勇，等.1961年以來海河流域干旱時空變化特征分析[J].水科學(xué)進(jìn)展，2013，24（1）：3441.（YAN Denghua，YUAN Zhe，YANG Zhiyong，et al.Spatial and temporal changes in drought since 1961 in Haihe River Basin[J].Advances in Water Science，2013，24（1）：3441.（in Chinese））

[25] Wu Z Y，Lu G H，Wen L，et al.Reconstructing and analyzing China′s fiftynine year（19512009）drought history using hydrological model simulation[J].Hydrol Earth Syst Sci Discuss，2011，8：18611893.

[26] 房穩(wěn)靜，張雪芬，趙龍，等.河南省冬小麥干旱災(zāi)害的空間分布研究[J].氣象與環(huán)境科學(xué)，2007，30（2）：6163.（FANG Wenjing ，ZHANG Xuefen，ZHAO Long，et al.Study on the winter wheat drought spatial distribution in Henan province[J].Meteorological and Environmental Sciences，2007，30（2）：6163.（in Chinese））

[27] 郭晶，景元書，王春林，等.基于土壤水分平衡模型的廣東干旱時空分布特征[J].中國農(nóng)業(yè)氣象，2008，29（3）：353357.（GUO Jing，JING Yuanshu，WANG Chunlin，et al.Spatial and temporal variations of drought in Guangdong province based on soil moisture balance model[J].Chinese Journal of Agrometeorology，2008，29（3）：353357.（in Chinese））

[28] 尹正杰，黃薇，陳進(jìn).基于土壤墑情模擬的農(nóng)業(yè)干旱動態(tài)評估[J].灌溉排水學(xué)報，2009，28（3）：58.（YIN Zhengjie，HUANG Wei，CHEN Jin.Dynamic evaluation of agricultural drought based on soil moisture simulation[J].Journal of Irrigation and Drainage，2009，28（3）：58.（in Chinese））

[29] 陳曉楠，段春青，劉昌明，等.基于兩層土壤計算模式的農(nóng)業(yè)干旱風(fēng)險評估模型[J].農(nóng)業(yè)工程學(xué)報，2009，25（9）：555.（CHEN Xiaonan，DUAN Chunqing，LIU Changming，et al.Model of risk assessment for agricultural drought based on twolayer soil computing model[J].Transactions of the CSAE，2009，25（9）：555.（in Chinese））

[30] Hunt E D，Hubbard K G，Wilhite D A，et al.The development and evaluation of a soil moisture index[J].Int J Climatol，2009，29：747759.

[31] 王慧軍.河北省糧食綜合生產(chǎn)能力研究[M].石家莊：河北科學(xué)技術(shù)出版社，2010.（WANG Huijun.Research of grain comprehensive productivity in Hebei Province[M].Shijiazhuang：Hebei Science & Technology Press，2010.（in Chinese））

[32] 王志敏，蘭林旺，王樹安.詮釋“吳橋模式”——冬小麥節(jié)水省肥簡化高產(chǎn)“四統(tǒng)一”栽培技術(shù)[J].中國農(nóng)資，2004（z1），5557.（WANG Zhimin，LAN Linwnag，WANG Shuan.The annotation of “Wuqiao Mode"cultivation technique system of winter wheat with unify of watersaving，fertilizersaving，highyielding and simple cultivation.China Agricultural Means of Production，2004（z1），5557.（in Chinese））

[33] 王樹安，蘭林旺，周殿璽，等.華北地區(qū)冬小麥節(jié)水高產(chǎn)高效栽培模式——吳橋模式[EB/OL].http：//www.jsgg.com.cn/Index/Display.asp？NewsID=8194，20140909.（WANG Shuan，LAN Linwang，ZHOU Dianxi，et al.High yield and high efficient cultivation techniques of winter wheat in northern China“Wuqiao Mode"[EB/OL].http：//www.jsgg.com.cn/Index/Display.asp？NewsID=8194，20140909.（in Chinese））

[34] 韓琳琳.應(yīng)用農(nóng)業(yè)生產(chǎn)系統(tǒng)模型分析海河流域平原區(qū)農(nóng)田冬小麥-夏玉米輪作體系的水氮利用效率[D].北京，中國農(nóng)業(yè)大學(xué)，2010.（HAN Linlin.Using agricultural production systems simulator to analyse water and nitrogen use efficiency of the winter wheatsummer maize rotation system in Hairiver Basin Plain[D].Beijing，China Agricultural，2010.（in Chinese））

[35] 韓琳琳，任理.應(yīng)用農(nóng)業(yè)生產(chǎn)系統(tǒng)模型分析冬小麥夏玉米輪作體系的農(nóng)田水氮利用效率I.模型參數(shù)的靈敏度分析與標(biāo)定[J].水利學(xué)報，2011，42（9）：10441050.（HAN Linlin，REN Li.Analysis of water and nitrogen use efficiency of winter wheatsummer maize rotation using agricultural production systems simulator I.Parameter sensitivity analysis and calibration[J].Journal of Hydraulic Engineering，2011，42（9）：10441050.（in Chinese））

[36] 河北省人民政府辦公廳，河北省統(tǒng)計局.河北農(nóng)村統(tǒng)計年鑒19952006[M].北京：中國統(tǒng)計出版社，19952006.（The General Office of the Hebei Provincial Government，Hebei Provincial Bureau of Statistics.Hebei Rural Statistical Yearbook 19952006[M].Beijing：China Statistics Press，19952006（in Chinese））

[37] 韓琳琳，任理.應(yīng)用農(nóng)業(yè)生產(chǎn)系統(tǒng)模型分析冬小麥夏玉米輪作體系的農(nóng)田水氮利用效率II.模型的模擬驗證與情景分析[J].水利學(xué)報，2011，42（10）：11851191.（HAN Linlin，REN Li.Analysis of water and nitrogen use efficiency of winter wheatsummer maize rotation using agricultural production systems simulator II.Model validation and scenario analysis[J].Journal of Hydraulic Engineering，2011，42（10）：11851191.（in Chinese））

[38] 中華人民共和國國家質(zhì)量監(jiān)督檢驗檢疫總局，中國國家標(biāo)準(zhǔn)化管理委員會.農(nóng)業(yè)干旱等級[S].2015.（General Administration of Quality Supervision，Inspection and Quarantine of the People′s Republic of China，Standardization Administration of the People′s Republic of China.Grade of agricultural drought[S].2015（in Chinese））