張可心 張謀草 劉翔 張洪芬 張?zhí)旆?/p>

摘要：基于1971—2018年隴東地區(qū)的逐月氣象資料，使用標(biāo)準(zhǔn)化降水指數(shù)（SPI）作為干旱指標(biāo)，采用線性趨勢(shì)和Mann-Kendall突變檢驗(yàn)，多尺度定量分析隴東地區(qū)近48年的氣象干旱特征，并結(jié)合冬小麥產(chǎn)量資料，研究SPI表征的氣象干旱對(duì)冬小麥產(chǎn)量的影響。結(jié)果顯示，1971—2018年隴東地區(qū)發(fā)生干旱的概率約為30%，發(fā)生中等以上強(qiáng)度干旱的概率為15.8%;年SPI分析結(jié)果顯示，隴東地區(qū)具有逐漸濕潤(rùn)化的趨勢(shì)，20世紀(jì)90年代（1990s）最干旱，2000年后干旱顯著緩解;季節(jié)SPI分析結(jié)果表明，春季干旱化趨勢(shì)明顯，且干旱化從1990年開始，其余季節(jié)均具有濕潤(rùn)化趨勢(shì)，其中冬季的濕潤(rùn)化速度最快;春季、秋季干旱對(duì)冬小麥產(chǎn)量的影響最大，尤其是干旱強(qiáng)度對(duì)產(chǎn)量的影響最為顯著，通過了0.001水平的顯著檢驗(yàn)。SPI所表征的氣象干旱指標(biāo)與冬小麥的相關(guān)性高于降水量與其相關(guān)性，表明SPI對(duì)于農(nóng)業(yè)干旱分析具有一定的指導(dǎo)意義。

關(guān)鍵詞：隴東黃土高原;標(biāo)準(zhǔn)化降水指數(shù)（SPI）;干旱特征;冬小麥產(chǎn)量;影響

中圖分類號(hào)： S162.5 ?文獻(xiàn)標(biāo)志碼： A ?文章編號(hào)：1002-1302（2020）23-0233-08

干旱是對(duì)人類社會(huì)影響最為嚴(yán)重的氣象災(zāi)害之一，造成干旱的直接原因有偶發(fā)性或周期性的降水量減少、蒸發(fā)量增大等[1]。在全球氣候變暖的背景下，發(fā)生大范圍干旱、高溫的時(shí)間、頻率和強(qiáng)度明顯增加，干旱化趨勢(shì)也隨之加劇，不但造成水資源短缺、荒漠化加劇、沙塵暴頻發(fā)等諸多不利的生態(tài)影響，同時(shí)也對(duì)農(nóng)業(yè)生產(chǎn)造成巨大損失[2]。西北地區(qū)處于干旱半干旱地區(qū)，對(duì)于全球氣候變化的響應(yīng)非常敏感，加強(qiáng)對(duì)干旱演變特征的研究有利于認(rèn)識(shí)全球氣候變化背景下的區(qū)域性響應(yīng)，對(duì)生態(tài)文明建設(shè)及防災(zāi)減災(zāi)工作具有現(xiàn)實(shí)意義[3-4]。

目前，國(guó)內(nèi)外研究者關(guān)于干旱指標(biāo)已做了大量研究，各個(gè)指標(biāo)本身各有利弊，研究區(qū)域也具有其自身的自然特征及社會(huì)特征。因此，根據(jù)研究目的，揚(yáng)長(zhǎng)避短地選擇合適的干旱指標(biāo)是研究區(qū)域性干旱的基礎(chǔ)[5]?；诮邓康母珊抵笜?biāo)有標(biāo)準(zhǔn)化降水指標(biāo)（SPI）、Z指數(shù)[6]，基于多因子的干旱指標(biāo)有地表濕潤(rùn)指數(shù)[7]、Palmer指數(shù)等[8]。帕默爾干旱指標(biāo)（palmer drought severity index，簡(jiǎn)稱PDSI）和標(biāo)準(zhǔn)化降水指標(biāo)（standardized precipitation index，簡(jiǎn)稱SPI）是全球干旱研究中運(yùn)用最廣泛的2個(gè)指標(biāo)。帕默爾干旱指標(biāo)PDSI計(jì)算過程復(fù)雜，參數(shù)區(qū)域性較強(qiáng)，并且具有固定的時(shí)間尺度，而標(biāo)準(zhǔn)化降水指標(biāo)只須以降水作為輸入量，計(jì)算簡(jiǎn)便易行且具有多種時(shí)間尺度，使同一個(gè)指標(biāo)可以反映不同時(shí)間尺度的干濕變化特征[9]。

基于SPI指數(shù)的干旱研究一般側(cè)重于對(duì)干旱時(shí)空分布特征的研究[10-12]，與糧食作物產(chǎn)量變化相結(jié)合的干旱特征檢測(cè)與影響研究較少。劉維等研究了SPI與東北春玉米產(chǎn)量之間的關(guān)系發(fā)現(xiàn)，玉米成災(zāi)率隨SPI6的變大而變小，當(dāng)SPI6處于干旱或者過濕情況下將處于低產(chǎn)量水平[13];朱玲玲等以SPI作為干旱指標(biāo)對(duì)河南冬小麥生育期變化特征進(jìn)行分析，發(fā)現(xiàn)河南省返青-抽穗-成熟期和全生育期干旱對(duì)冬小麥氣候產(chǎn)量影響較大[14]。山西省關(guān)于玉米產(chǎn)量與SPI的研究發(fā)現(xiàn)，對(duì)其產(chǎn)量影響較大的是夏季干旱，中等以上強(qiáng)度的夏季干旱對(duì)玉米產(chǎn)量的影響更為顯著[15]。甘肅省慶陽市地處隴東黃土高原，屬于黃土高原丘陵溝壑區(qū)的生態(tài)脆弱區(qū)，以雨養(yǎng)農(nóng)業(yè)為主，由降水量減少導(dǎo)致的干旱對(duì)隴東地區(qū)農(nóng)業(yè)生產(chǎn)和生態(tài)環(huán)境的影響顯著[16]。冬小麥?zhǔn)请]東地區(qū)最主要的糧食作物之一，在慶陽市種植面積最大，由于生長(zhǎng)周期長(zhǎng)，不同時(shí)段的干旱變化對(duì)冬小麥產(chǎn)量都有一定的影響，而對(duì)冬小麥的影響同樣也可以反映對(duì)其他作物的影響。本研究以SPI作為干旱評(píng)價(jià)指標(biāo)，分析1971—2018年隴東黃土高原不同時(shí)間尺度的氣象干旱特征及對(duì)冬小麥產(chǎn)量的影響，以期為氣象干旱檢測(cè)、干旱災(zāi)害預(yù)防及冬小麥安全生產(chǎn)提供參考依據(jù)。

1 資料與方法

1.1 研究區(qū)域概況

隴東黃土高原屬于黃河中游的內(nèi)陸地區(qū)，東倚子午嶺，西接六盤山，所在行政區(qū)域主要為甘肅省慶陽市（106°20′～108°45′E、35°15′～37°10′N），屬于暖溫帶半干旱半濕潤(rùn)氣候，光熱條件適宜，年降水量為400～700 mm，是黃土高原區(qū)降水最多的地區(qū)。慶陽市馬嶺以北的環(huán)縣、華池及慶城的中北部為黃土丘陵溝壑區(qū)，溝壑縱橫，水土流失嚴(yán)重，馬嶺以南塬面平坦，其中董志塬面積為910.7 km2，是我國(guó)最大、較完整的一片黃土塬地。隴東黃土高原塬地面積廣，約占土地總面積的1/3，坡度小，盛產(chǎn)小麥、玉米、油料作物及特色小雜糧，有“隴上糧倉”的美譽(yù)。

1.2 站點(diǎn)選擇與資料

本研究采用慶陽市8個(gè)國(guó)家氣象站1971—2018年的逐月降水量數(shù)據(jù)與慶陽市統(tǒng)計(jì)局1978—2018年冬小麥年產(chǎn)量數(shù)據(jù)，通過MATLAB計(jì)算得出SPI值及四季SPI的Mann-Kendall（M-K）突變檢驗(yàn)，SPI時(shí)間變化趨勢(shì)用氣候傾向率進(jìn)行計(jì)算，用MATLAB、Excel等軟件繪制圖表。在統(tǒng)計(jì)季節(jié)降水量時(shí)，3—5月為春季，6—8月為夏季，9—11月為秋季，12月至次年2月為冬季。

1.3 標(biāo)準(zhǔn)化降水指數(shù)的計(jì)算

標(biāo)準(zhǔn)化降水指數(shù)的計(jì)算原理是McKee等在1993年提出的1個(gè)相對(duì)簡(jiǎn)便易行的干旱指數(shù)，該指數(shù)基于自然降水的Γ分布，基于正態(tài)分布標(biāo)準(zhǔn)化得出，計(jì)算過程詳見文獻(xiàn)[17-18]。SPI的計(jì)算有多種尺度（1、3、6、12、24月等），對(duì)于時(shí)間尺度短的SPI，由于受短時(shí)間降水的影響，導(dǎo)致SPI上下波動(dòng)頻繁，值域也較大，從而可反映短期旱澇變化，而時(shí)間較長(zhǎng)的SPI對(duì)短期降水的響應(yīng)減慢，旱澇變化比較穩(wěn)定，可反映長(zhǎng)期旱澇的變化趨勢(shì)[19]。本研究根據(jù)CB/T 20481—2006《氣象干旱等級(jí)標(biāo)準(zhǔn)》[20]，結(jié)合隴東地區(qū)的實(shí)際情況，劃分出干旱等級(jí)（表1）。

1.4 干旱強(qiáng)度

某段時(shí)間內(nèi)的干旱強(qiáng)度可以用SPI表示，在發(fā)生干旱（SPI≤-0.5）的情況下，SPI絕對(duì)值越大，表示干旱越嚴(yán)重。某區(qū)域的多年干旱程度可由下式表示：

式中：|SPIi|為發(fā)生干旱時(shí)的SPI絕對(duì)值，具體的干旱分級(jí)見表1;m為發(fā)生干旱的時(shí)間（月數(shù)）。

2 結(jié)果與分析

2.1 不同時(shí)間尺度的SPI特征

不同時(shí)間尺度的SPI可以反映不同類型的旱澇情況，通過多時(shí)間尺度的SPI綜合分析可以進(jìn)行綜合旱澇情況分析[21]。為了對(duì)隴東地區(qū)的干旱情況進(jìn)行綜合評(píng)價(jià)，本研究基于慶陽市8個(gè)氣象站點(diǎn)48年的逐月降水量資料，著重分析SPI的1、3、6、12月尺度的干旱變化情況，并對(duì)各個(gè)時(shí)間尺度的干旱等級(jí)進(jìn)行劃分。分析圖1-a的SPI1指數(shù)可知，出現(xiàn)干旱的概率為31.9%，其中中等以上等級(jí)干旱分布如下：極旱10個(gè)月（1.7%），重旱15個(gè)月（2.6%），中旱67個(gè)月（11.6%）。分析圖1-b的SPI3指數(shù)可知，出現(xiàn)干旱的概率為29.1%，其中極旱13個(gè)月（2.3%），重旱29個(gè)月（5.0%），中旱47個(gè)月（8.2%）。分析圖1-c的SPI6指數(shù)可知，干旱發(fā)生概率為32.1%，其中極旱9個(gè)月（1.6%），重旱29個(gè)月（5.0%），中旱51個(gè)月（8.9%）。分析圖1-d的SPI12指數(shù)可知?干旱發(fā)生概率為32.9%?7個(gè)月為極旱（1.2%），33個(gè)月為重旱（5.8%），52個(gè)月為中旱（9.2%）。由以上分析可以看出，隴東干旱頻發(fā)，平均30%以上的月份發(fā)生不同程度的干旱，用不同時(shí)間尺度表征的中等以上干旱狀況較為一致，中等以上程度的干旱月份占比約為15.8%，SPI指數(shù)起伏較大，表明氣候干濕變化波動(dòng)幅度大，旱澇時(shí)間分布不均勻。

甘肅省氣象災(zāi)害統(tǒng)計(jì)結(jié)果表明，1986—1987年隴東出現(xiàn)了數(shù)十年一遇的持續(xù)性干旱，旱情嚴(yán)重，50%以上的人畜飲水困難。SPI1分析結(jié)果顯示，1986年7—10月和1987年8—9月出現(xiàn)了輕度以上干旱;SPI3分析結(jié)果顯示，1986、1987年8—10月均出現(xiàn)中旱;SPI6分析結(jié)果顯示，1986年8月至次年2月出現(xiàn)了重旱;而SPI12分析結(jié)果顯示，1986年9月至1987年底出現(xiàn)了連續(xù)性重旱。1995年隴東地區(qū)出現(xiàn)春夏連旱的情況，農(nóng)作物大量減產(chǎn)，冬小麥單位面積產(chǎn)量（簡(jiǎn)稱單產(chǎn)）降至882 kg/hm2，SPI1分析結(jié)果顯示，2—7月出現(xiàn)了中等以上的干旱，5月達(dá)到特旱;SPI3分析結(jié)果顯示，3—7月出現(xiàn)了特旱。1997年隴東地區(qū)出現(xiàn)了嚴(yán)重的夏秋連旱，SPI3分析結(jié)果顯示，1997年6—10月出現(xiàn)了中等以上干旱，其中重旱以上的時(shí)間有3個(gè)月。不同時(shí)間尺度的SPI分析結(jié)果表明，1995、1997年的干旱與1986—1987年的干旱具有相同特征，即隨著時(shí)間尺度的增大，干旱開始時(shí)間和結(jié)束時(shí)間有相應(yīng)推遲，并且干旱等級(jí)有所變化，反映了前期降水積累的影響因素。對(duì)于類似1986—1987年這種持續(xù)時(shí)間長(zhǎng)的干旱，SPI12表征的時(shí)間尺度與干旱程度較準(zhǔn)確，而對(duì)于1995、1997年這種季節(jié)性干旱，與實(shí)況對(duì)應(yīng)最好的是SPI3。不同時(shí)間尺度的SPI表征的干旱年際變化與隴東歷史上實(shí)際發(fā)生的干旱年對(duì)應(yīng)的年份一致，可見用SPI分析隴東地區(qū)的干旱變化是可行的[22]。

2.2 旱澇的年際變化特征

2.2.1 SPI的年際變化特征 1971—2018年隴東地區(qū)的年SPI總體變化趨勢(shì)為波動(dòng)升高（圖2），傾向率為0.061/10年，有8年出現(xiàn)中等程度及以上的干旱，其中發(fā)生中旱的時(shí)間為5年（1972、1979、1982、1987、2008年），發(fā)生重旱的時(shí)間為3年（1986、1995、1997年）。有8年出現(xiàn)中等程度以上的偏濕，其中中濕4年（1983、1988、2017、2018年），重濕4年（1975、1990、2003、2013年）。從干旱程度年代際的分布來看，1990s（20世紀(jì)90年代）干旱程度最為嚴(yán)重，出現(xiàn)了2年重旱，其次是1980s（20世紀(jì)80年代），出現(xiàn)過1年重旱與2年中旱，1970s（20世紀(jì)70年代）出現(xiàn) 2年中旱;到了21世紀(jì)，中等程度及以上的干旱僅出現(xiàn)過1年。根據(jù)以上分析可知，隴東地區(qū)1990s最干旱，其次是1980s，1970s相對(duì)較濕潤(rùn)，2000s（21世紀(jì)00年代）以后干旱程度大幅度減弱。中度及以上干濕年份占比為33%，最干、最濕的年份分別為1997年（-1.882，重旱）、2003年（1.944，重濕）。本研究結(jié)果與張可心等得出的“慶陽市20世紀(jì)80年代濕潤(rùn)，90年代最干旱，2000后恢復(fù)濕潤(rùn)”的結(jié)論[23]基本一致。由以上分析可知，隴東地區(qū)氣候干濕程度的發(fā)展趨勢(shì)具有波動(dòng)中逐漸變濕潤(rùn)的趨勢(shì)，這與黃小燕等的研究結(jié)論[24]一致。

2.2.2 SPI的季節(jié)變化特征 由圖3可知，除春季以外，其余季節(jié)的SPI均呈現(xiàn)升高的趨勢(shì)。春季SPI呈緩慢下降的趨勢(shì)，變化傾向率為0.047/10年，表明春季有干旱化趨勢(shì)，干旱發(fā)生頻率為27%，中等及以上干旱季與濕潤(rùn)季分別為8、5年，在1992—2011年間出現(xiàn)高頻率的干旱，其中最嚴(yán)重的是1995年，干旱等級(jí)為極旱（SPI=-3.01），最濕潤(rùn)的是1998年（SPI=2.11）。夏季SPI的變化傾向率為0.05/10年，干旱發(fā)生頻率為27%，中等及以上干旱季、濕潤(rùn)季均為8年?1997、2018年的SPI值分別為-1.92（重旱）、1.99（重濕），這2年分別為夏季最干旱、最濕潤(rùn)的年份。秋季SPI的變化傾向率為0.057/10年，干旱發(fā)生頻率為31%，中等以上干旱季、濕潤(rùn)季分別為5、10年，秋季最干旱、最濕潤(rùn)年份分別為1986、1975年，SPI值分別為-1.97（重旱）、2.20（重濕）。冬季SPI的變化傾向率為0.077/10年，干旱發(fā)生頻率為21%，中等以上干旱季、濕潤(rùn)季分別7、8年，1992、1999年出現(xiàn)過2次極旱，1992、2008年分別為最干旱、最濕潤(rùn)的年份，SPI分別為 -3.92（極旱）、1.50（中濕），2000年之后僅出現(xiàn)過1次輕旱，可見21世紀(jì)以后隴東地區(qū)冬季的干旱次數(shù)大幅度減少。以上分析結(jié)果表明，春季具有干旱化趨勢(shì)，其余季節(jié)均具有濕潤(rùn)化趨勢(shì)，冬季的濕潤(rùn)化趨勢(shì)最為顯著;干旱發(fā)生頻率最高的是秋季，其次是春季、夏季，最低的是冬季。秋季的干濕變化跨度大，正值冬小麥播種期，需要防范干旱的不利影響;冬季整體偏旱，但是由于冬季的農(nóng)業(yè)活動(dòng)少，加上濕潤(rùn)化趨勢(shì)明顯，因此對(duì)隴東地區(qū)農(nóng)業(yè)的影響也相對(duì)較小。

[8]王勁松，黃玉霞，馮建英，等. 徑流量Z指數(shù)與Palmer指數(shù)對(duì)河西干旱的監(jiān)測(cè)[J]. 應(yīng)用氣象學(xué)報(bào)，2009，20（4）：471-477.

[9]Mishra A K，Singh V P. A review of drought concepts[J]. Journal of Hydrology，2010，391（1/2）：202-216.

[10]龔艷冰，張繼國(guó)，劉高峰，等. 基于SPI指數(shù)與R/S分析的曲靖市干旱特征研究[J]. 干旱地區(qū)農(nóng)業(yè)研究，2015，33（3）：273-277.

[11]那音太. 基于SPI指數(shù)的近50a內(nèi)蒙古地區(qū)干旱特征分析[J]. 干旱區(qū)資源與環(huán)境，2015，29（5）：161-166.

[12]王江濤. 基于華北地區(qū)氣象指數(shù)SPI干旱時(shí)空分析[J]. 水土保持研究，2019，26（4）：203-207.

[13]劉 維，李祎君，何 亮，等. 基于SPI判定的東北春玉米生長(zhǎng)季干旱對(duì)產(chǎn)量的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2018，34（22）：121-127.

[14]朱玲玲，張竟竟，李治國(guó)，等. 基于SPI的河南省冬小麥生育期干旱時(shí)空變化特征分析[J]. 灌溉排水學(xué)報(bào)，2018，37（5）：51-58.

[15]郭 偉，李 瑩，杜莉麗. 基于SPI的山西省1972—2012年春夏干旱特征及對(duì)玉米產(chǎn)量的影響分析[J]. 干旱地區(qū)農(nóng)業(yè)研究，2018，36（1）：230-236.

[16]張耀宗，張 勃，劉艷艷，等. 近50年隴東黃土高原干旱特征及未來變化趨勢(shì)分析[J]. 干旱地區(qū)農(nóng)業(yè)研究，2017，35（2）：263-270.

[17]Edwards D C，McKee T B.Characteristics of 20th century drought in the United States at multiple time scales[D]. USA：Colorado State University，1997：18-21.

[18]袁文平，周廣勝. 干旱指標(biāo)的理論分析與研究展望[J]. 地球科學(xué)進(jìn)展2004（6）：982-991.

[19]車少靜，李春強(qiáng)，申雙和，等. 基于SPI的近41年（1965—2005）河北省旱澇時(shí)空特征分析[J]. 中國(guó)農(nóng)業(yè)氣象，2010，31（1）：137-143.

[20]中國(guó)氣象局. 氣象干旱等級(jí)：GB/T 20481—2006[S]. 北京：中國(guó)標(biāo)準(zhǔn)出版社，2006.

[21]張利利，周俊菊，張恒瑋，等. 基于SPI的石羊河流域氣候干濕變化及干旱事件的時(shí)空格局特征研究[J]. 生態(tài)學(xué)報(bào)，2017，37（3）：996-1007.

[22]董安祥. 中國(guó)氣象災(zāi)害大典：甘肅卷[M]. 北京：氣象出版社，2005.

[23]張可心，張謀草，杜 軍. 1981—2016年慶陽市地表濕潤(rùn)狀況時(shí)空變化特征分析[J]. 中國(guó)農(nóng)學(xué)通報(bào)，2019，35（13）：107-106.

[24]黃小燕，張明軍，賈文雄，等. 中國(guó)西北地區(qū)地表干濕變化及影響因素[J]. 水科學(xué)進(jìn)展，2011，22（2）：151-159.李躍飛，李 彬，陶加樂，等. 宿豫保護(hù)地土壤硝酸鹽污染與土壤酸化、鹽漬化狀況調(diào)查與分析[J]. 江蘇農(nóng)業(yè)科學(xué)，2020，48（23）：241-245.