陳曉旸，羅福命，曾運(yùn)東，何曉鴻，李茂軍，朱漢輝，鄧建良

（1.廣東省韶關(guān)市氣象局，廣東韶關(guān)512000；2.廣東煙草韶關(guān)市有限公司，廣東韶關(guān)512000）

根據(jù)IPCC 第六次評(píng)估報(bào)告，預(yù)計(jì)到2040年前后全球海陸表面氣溫將升高1.5~1.6 ℃，且人類活動(dòng)是導(dǎo)致極端天氣的主要驅(qū)動(dòng)力[1]。當(dāng)前，我國(guó)極端天氣頻發(fā)，給人民的生命財(cái)產(chǎn)安全造成巨大的危害。其中干旱是影響范圍最廣，持續(xù)時(shí)間最長(zhǎng)和災(zāi)情最嚴(yán)重的自然災(zāi)害之一。干旱會(huì)嚴(yán)重破壞地區(qū)水資源的供需平衡，從而制約地區(qū)的農(nóng)業(yè)生產(chǎn)和社會(huì)經(jīng)濟(jì)發(fā)展[2]。各國(guó)學(xué)者針對(duì)如何通過(guò)有效途徑進(jìn)行干旱監(jiān)測(cè)與評(píng)價(jià)都開展了大量的研究[3-5]。目前在研究中廣泛應(yīng)用的干旱指數(shù)有帕爾默干旱指數(shù)(Palmer drought severity index，PDSI)、標(biāo)準(zhǔn)化降水指數(shù)(Standardized precipitation index，SPI)、綜合氣象干旱指數(shù)(Composite index of meteorological drought，CI)和標(biāo)準(zhǔn)化降水蒸散發(fā)指數(shù)(Standardized precipitation evapotranspiration index，SPEI)等[6]。其中標(biāo)準(zhǔn)化降水蒸散指數(shù)(SPEI)是Vicente-Serrano等在SPI指數(shù)的基礎(chǔ)上綜合考慮了降水和溫度因子所提出的一個(gè)新的氣象干旱指數(shù)，是監(jiān)測(cè)干旱較優(yōu)的工具[7]，目前已成為研究干旱中使用頻率最高的干旱指數(shù)。如趙茹欣等應(yīng)用SPEI分析了黑龍江省1959-2014年干旱時(shí)空分布特征，發(fā)現(xiàn)其具有整體一致型和西北-東南對(duì)立分布兩種態(tài)勢(shì)[8]。Meza 等使用SPEI 指數(shù)研究智利的干旱特征，分析了其不同尺度下的干旱變化趨勢(shì)[9]。賴力等將SPEI 指數(shù)與植被狀態(tài)指數(shù)(VCI)結(jié)合起來(lái)，指出甘肅省的農(nóng)業(yè)干旱對(duì)3 個(gè)月尺度的氣象干旱響應(yīng)最靈敏[10]。

根據(jù)統(tǒng)計(jì)，21世紀(jì)以來(lái)，廣東省約24%的耕地經(jīng)常遭遇干旱[11]，當(dāng)前華南地區(qū)的干旱及水資源安全問(wèn)題已成為人們的關(guān)注熱點(diǎn)之一。而廣東省韶關(guān)市森林資源和耕地資源豐富，既是全國(guó)重點(diǎn)林區(qū)，也是珠三角地區(qū)重要的農(nóng)副產(chǎn)品和煙草供應(yīng)基地，因此對(duì)氣候變化背景下韶關(guān)市的干旱特征深入研究具有重要的意義。丘蔚天[12]利用日降水的持續(xù)時(shí)間作為干旱的判別標(biāo)準(zhǔn)，分析了1967-2002年韶關(guān)市的持續(xù)干旱特征，但其對(duì)干旱的定義僅考慮了降水因素。本文利用綜合降水和溫度影響的SPEI 干旱指數(shù)，采用趨勢(shì)分析、ArcGIS 空間分析和MMK 檢驗(yàn)等方法，對(duì)韶關(guān)市1971-2020年干旱時(shí)空變化特征進(jìn)行了探討，并利用FAO 推薦的作物系數(shù)法，探究了煙草生長(zhǎng)期需水量與SPEI指數(shù)之間的聯(lián)系，以期為韶關(guān)的農(nóng)業(yè)抗旱、灌溉以及煙草種植管理提供科學(xué)指導(dǎo)。

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

本次研究選取韶關(guān)市8 個(gè)國(guó)家氣象站1971-2020年逐日氣象數(shù)據(jù)，氣象要素具體包括：10 m 平均風(fēng)速、日平均氣溫、日最高，最低氣溫、日雨量、相對(duì)濕度、平均氣壓和光照時(shí)數(shù)。8 個(gè)氣象站點(diǎn)分布如圖1所示。數(shù)據(jù)經(jīng)過(guò)了較為嚴(yán)格的修訂，包括數(shù)據(jù)錯(cuò)誤訂正、遺失數(shù)據(jù)補(bǔ)漏、無(wú)效數(shù)據(jù)剔除，并經(jīng)過(guò)氣候極值、內(nèi)部一致性、空間一致性等質(zhì)量控制，確保了數(shù)據(jù)的均一性和完整性。

圖1 韶關(guān)市氣象站點(diǎn)分布圖Fig.1 Distribution map of meteorological stations in Shaoguan City

根據(jù)廣東煙草韶關(guān)市有限公司提供的資料，整理出了韶關(guān)不同產(chǎn)區(qū)的烤煙生長(zhǎng)發(fā)育期（表1）。因韶關(guān)的煙草種植主要集中在南雄和始興兩個(gè)產(chǎn)區(qū)，綜合考慮后，將3~6月定為韶關(guān)烤煙的生長(zhǎng)發(fā)育期。

表1 韶關(guān)各產(chǎn)區(qū)烤煙生長(zhǎng)發(fā)育期Tab.1 Growth period of flue cured tobacco in Shaoguan

1.1 標(biāo)準(zhǔn)化降水蒸散發(fā)指數(shù)SPEI

綜合比較了SPI、PDSI 和SPEI 三種干旱指數(shù)后，選取了融合蒸散發(fā)影響的SPEI 指數(shù)來(lái)監(jiān)測(cè)干旱的變化特征。SPEI 的具體計(jì)算步驟如下：

(1)使用彭曼—蒙蒂斯（Penman-Monteith）公式[13]對(duì)潛在蒸散發(fā)（PET）進(jìn)行計(jì)算，計(jì)算公式為：

式中：PET為潛在蒸散發(fā)量，mm/d；Rn為冠層表面凈輻射，MJ/(m2?d)；G為土壤熱通量，MJ/(m2?d)；T為平均氣溫，℃；es為飽和水汽壓，kPa；ea為實(shí)際水汽壓，kPa；Δ 為飽和水汽壓-氣溫關(guān)系曲線在T處的切線斜率，kPa/℃；γ為濕度計(jì)常數(shù)，kPa/℃；U2為近地面2 m高處的風(fēng)速，m/s。

(2)利用式(1)計(jì)算逐月降水與潛在蒸散發(fā)差值：

式中：Di為月凈降水量，mm；Pi為月降水量，mm；PETi為月潛在蒸發(fā)量，mm。

(3)基于3 參數(shù)log-logistic 概率分布，對(duì)數(shù)據(jù)序列進(jìn)行處理。其概率分布的累計(jì)函數(shù)為：

式中：α 為尺度參數(shù)；β為形狀參數(shù)；δ為Origin 參數(shù)，3 參數(shù)通過(guò)線性矩?cái)M合獲得；F(x)為累積概率分布函數(shù)。

(4)對(duì)累計(jì)概率分布函數(shù)F(x)進(jìn)行標(biāo)準(zhǔn)化正態(tài)處理，超過(guò)某個(gè)Di值的概率為P= 1-F(x)。

當(dāng)P≤0.5 時(shí)，參數(shù)標(biāo)準(zhǔn)化降水蒸散發(fā)指數(shù)SPEI計(jì)算公式為：

當(dāng)P＞0.5時(shí)，：

其他參數(shù)分別為C0= 2.515 517，C1= 0.802 853，C2=0.010 328，d1= 1.432 788，d2= 0.189 269，d3= 0.001308。

依據(jù)國(guó)家氣候中心2006年制定的《氣象干旱等級(jí)》[14]，對(duì)SPEI指數(shù)進(jìn)行干旱等級(jí)劃分，結(jié)果見表2。

表2 SPEI指數(shù)干旱等級(jí)劃分Tab.2 Drought grade classification of SPEI

本文計(jì)算干旱發(fā)生頻率的公式如下：

式中：n為干旱發(fā)生次數(shù)；N為總次數(shù)。

若相應(yīng)月份SPEI 指數(shù)小于或等于-0.5，就記為發(fā)生一次干旱。本次研究中將每年的5、8、11 和次年2月的SPEI-3 指數(shù)分別代表當(dāng)年春、夏、秋、冬的干旱狀況，將每年12月的SPEI-12指數(shù)代表當(dāng)年的年干旱狀況。

1.2 干旱站次比R

定義區(qū)域內(nèi)干旱發(fā)生的站點(diǎn)數(shù)與全部站點(diǎn)數(shù)之比為干旱站次比，計(jì)算公式為：

式中：m為發(fā)生干旱的站點(diǎn)數(shù)；M為研究區(qū)域內(nèi)的站點(diǎn)總數(shù)。當(dāng)研究區(qū)域一半以上的站點(diǎn)發(fā)生干旱時(shí)，即R≥50%時(shí)，定義為全域性干旱；當(dāng)R＜10%時(shí)認(rèn)為無(wú)明顯干旱，10%≤R＜50%則定義為區(qū)域性干旱。

1.3 Sen’s 斜率估計(jì)

Sen’s 斜率估計(jì)方法主要用于計(jì)算韶關(guān)煙草生長(zhǎng)期間降水、氣溫和SPEI 指數(shù)年際變化的斜率Qmed。對(duì)于時(shí)間序列xi=(x1,x2,…,xn)，Sen’s斜率的計(jì)算公式為：

式中：Median為取中值函數(shù)。

1.4 MMK（Modified Mann-Kendall）趨勢(shì)檢驗(yàn)方法

MMK 趨勢(shì)檢驗(yàn)方法主要用于檢驗(yàn)降水、氣溫和SPEI 趨勢(shì)變化的顯著性。原始的MK 趨勢(shì)檢驗(yàn)方法是基于假設(shè)序列獨(dú)立同分布，對(duì)于時(shí)間序列X=(x1,x2,…,xn)，定義統(tǒng)計(jì)量S為[15]：

假設(shè)各變量獨(dú)立同分布，則統(tǒng)計(jì)量S近似服從正態(tài)分布，其均值E(S)、方差Var(S)分別為[16]：

當(dāng)序列長(zhǎng)度n＞10時(shí)，計(jì)算Mann-Kendall統(tǒng)計(jì)量[17]：

通過(guò)給定的顯著性水平α，比較|Uc|與臨界值Uα的大小，即可判定序列趨勢(shì)的統(tǒng)計(jì)顯著性。

但實(shí)際中，序列的相關(guān)性會(huì)對(duì)Var(S)有明顯影響，因此本文采用Hamed 等改進(jìn)的MMK 趨勢(shì)檢驗(yàn)方法對(duì)存在相關(guān)性的序列進(jìn)行檢驗(yàn)分析，該方法主要對(duì)Var(S)進(jìn)行修正，具體方法如下：

對(duì)于序列X=(x1,x2,…,xn)，首先利用式(8)計(jì)算序列的趨勢(shì)估計(jì)量Qmed，并從序列中去除改趨勢(shì)項(xiàng)，得到與原序列對(duì)應(yīng)的平穩(wěn)序列Y=(y1,y2,…,yn)：

然后，求序列Y對(duì)應(yīng)的秩次序列Ri=1,2,..,n，并用其計(jì)算自相關(guān)系數(shù)[17]rk=1,2,..,n：

式中：Ri為yi的秩次，R是秩次的均值。然后，依據(jù)rk=1,2,..,n，計(jì)算修正系數(shù)η：

最后，改進(jìn)后的趨勢(shì)統(tǒng)計(jì)量S的方差Var*(S)為：

將Var*(S)帶入式(13)，可得到基于MMK 趨勢(shì)檢驗(yàn)方法的統(tǒng)計(jì)量，從而判斷某一顯著性水平下的序列趨勢(shì)顯著性。

1.5 煙草生長(zhǎng)期需水量計(jì)算

煙草生長(zhǎng)期需水量由生長(zhǎng)期內(nèi)逐旬需水量累積得到。逐旬需水量利用FAO 推薦的作物系數(shù)法計(jì)算，把標(biāo)準(zhǔn)條件下（長(zhǎng)勢(shì)良好，供水充足）的逐旬蒸散發(fā)量作為理論需水量，計(jì)算公式[18]：

式中：ETc為逐旬作物需水量，mm；ET0為逐旬參照作物蒸散量，mm，該變量只與氣象要素有關(guān)，反映了不同地區(qū)不同時(shí)期大氣蒸發(fā)能力對(duì)作物需水量的影響；Kc為逐旬作物系數(shù)，反映了作物蒸騰、土壤蒸發(fā)的綜合效應(yīng)，受作物類型、作物生長(zhǎng)狀況、氣候條件和土壤蒸發(fā)等多種因素影響。

因韶關(guān)與貴州地形都以山區(qū)丘陵為主，環(huán)境氣候也相似。參考袁義杰[19]等研究成果，3~6月韶關(guān)煙草生長(zhǎng)期內(nèi)共有12旬，按順序逐旬作物系數(shù)采用表3 中的參數(shù)。本文中的ET0同樣采用公式(1)計(jì)算，ET0=PET。

表3 煙草旬作物系數(shù)Tab.3 Tobacco ten-day crop coefficient

2 結(jié)果與分析

2.1 不同時(shí)間尺度SPEI指數(shù)的變化特征

本文主要基于1、3 和12 個(gè)月尺度的SPEI 指數(shù)進(jìn)行研究。韶關(guān)地區(qū)3 種不同時(shí)間尺度的SPEI 指數(shù)變化特征如圖2所示。1個(gè)月尺度（SPEI-1）數(shù)據(jù)量最大，能反映出干旱指數(shù)月際間的細(xì)微變化。根據(jù)SPEI-1，韶關(guān)市分別在2009年2月和2011年4月達(dá)到了特旱等級(jí)。3 個(gè)月尺度（SPEI-3）則可以反映干旱程度的季節(jié)性變化。由于其在計(jì)算中考慮了干旱的累積效應(yīng)，與SPEI-1 相比在某些月份反映的干旱程度會(huì)更嚴(yán)重。如SPEI-3 顯示韶關(guān)市在2019年12月為特旱，但該月的SPEI-1僅達(dá)到重旱等級(jí)。12個(gè)月尺度（SPEI-12）則主要反映干濕條件的年際變化，其波動(dòng)幅度最小。1971-2000年，韶關(guān)市共有7年發(fā)生干旱事件，無(wú)重旱，年內(nèi)發(fā)生干旱的比例為23%。而在2001-2020年，共有8年達(dá)到了干旱等級(jí)，其中有一年是重旱，年內(nèi)發(fā)生干旱的比例升至40%。因此，2000s 后，韶關(guān)市的干旱化有增強(qiáng)的趨勢(shì)，且發(fā)生重旱的可能性有所加大。

圖2 韶關(guān)市1、3和12個(gè)月尺度SPEI指數(shù)的時(shí)間變化Fig.2 Variation of SPEI at 1-month、3-month and 12-month scale in Shaoguan City

2.2 不同時(shí)間尺度干旱站次比

圖3 展示了1971-2020年韶關(guān)地區(qū)在年尺度和季節(jié)尺度上，干旱影響范圍的時(shí)間變化。從整體結(jié)果上看，四季中夏旱發(fā)生最為頻繁，且春、夏兩季較常發(fā)生全域性干旱事件。在過(guò)去的50a間，秋旱和冬旱的影響范圍都呈增加趨勢(shì)，春旱和夏旱的站次比則呈相反的下降趨勢(shì)。其中秋旱的站次比變化速率最為顯著，達(dá)到了3%/10a，說(shuō)明韶關(guān)市的秋旱越趨嚴(yán)重。另外，通過(guò)比較發(fā)現(xiàn)，除秋季外，其余三季的干旱站次比均在1990s 出現(xiàn)了最小值，而2010s 春、秋和冬季的干旱站次比為次小值。因此，需特別注意韶關(guān)市在春、冬兩季發(fā)生年際間旱澇急轉(zhuǎn)的情況。

韶關(guān)地區(qū)年尺度干旱站次比以3.22%/10a 的速率上升。2000s 以前年內(nèi)發(fā)生全域性干旱的比例為27%，而在2000s 后該比例升至45%。年干旱站次比同樣在1990s 出現(xiàn)了最小值，最大值出現(xiàn)在2000s，說(shuō)明1990s 在氣象上屬于濕潤(rùn)年。由此可見，近50a來(lái)，韶關(guān)地區(qū)年干旱影響范圍在不斷擴(kuò)大，年內(nèi)發(fā)生全域性干旱的風(fēng)險(xiǎn)在2000s后有明顯提高。

2.3 干旱頻率的時(shí)空變化特征

2.3.1 干旱頻率時(shí)間變化特征

基于SPEI-12 指數(shù)對(duì)50a 來(lái)韶關(guān)地區(qū)不同年代干旱事件的發(fā)生頻率進(jìn)行統(tǒng)計(jì)，結(jié)果如圖4所示。由圖4（a）可知，1971-2020年，韶關(guān)地區(qū)干旱事件的發(fā)生頻率達(dá)到了30%，其中主要以中度干旱為主（16%），輕度干旱次之（12%），重度干旱最少（2%)，無(wú)特旱發(fā)生。重旱發(fā)生在2000s，而各年代干旱發(fā)生頻率無(wú)顯著差別。因此，近50a韶關(guān)市的干旱事件主要以中度干旱為主，干旱發(fā)生次數(shù)無(wú)明顯增多，且在2000s后重度干旱出現(xiàn)的可能性增大。

同樣，利用SPEI-3 指數(shù)，發(fā)現(xiàn)韶關(guān)四季的干旱發(fā)生頻次基本相同，且主要以輕度干旱為主[見圖4（b）]。其中春季和冬季各發(fā)生了一次特旱事件，而春、夏和秋季都有重旱發(fā)生。所以，韶關(guān)市在春季和冬季需關(guān)注極端干旱事件的發(fā)生，其余兩季則主要關(guān)注重旱的出現(xiàn)。

圖4 1971-2020年韶關(guān)市不同年際干旱發(fā)生頻率及不同季節(jié)干旱發(fā)生頻數(shù)Fig.4 Frequency of drought in different ages and seasons in Shaoguan from 1971 to 2020

2.3.2 干旱頻率空間變化特征

韶關(guān)市季節(jié)和年干旱發(fā)生頻率的空間分布如圖5所示?？芍仃P(guān)南部的翁源和新豐是春旱和夏旱發(fā)生頻率的高值區(qū)，頻率超過(guò)了34%。秋旱和冬旱則主要發(fā)生在韶關(guān)的中西部，頻率在32%~36%之間。值得注意的是，除冬季外，韶關(guān)曲江的干旱發(fā)生頻率在其余三季都超過(guò)了34%。因此曲江是季節(jié)尺度干旱事件的易發(fā)生地。在年干旱發(fā)生頻率方面，韶關(guān)各市縣無(wú)明顯差異。從整體上看，韶關(guān)年干旱發(fā)生頻率呈現(xiàn)一個(gè)“北多南少”的分布規(guī)律。

圖5 1971-2019年韶關(guān)市季節(jié)和年干旱發(fā)生頻率的空間分布圖Fig.5 Spatial distribution of frequency of seasonal and annual drought in Shaoguan City

2.4 干旱對(duì)韶關(guān)煙草生長(zhǎng)發(fā)育期的影響

由于韶關(guān)烤煙的生長(zhǎng)發(fā)育期主要集中在3月至6月，因此將每年6月的SPEI-4 指數(shù)（以下簡(jiǎn)稱為SPEI 指數(shù)）代表當(dāng)年煙草生長(zhǎng)期內(nèi)的干旱狀況。利用Sen’s 斜率估計(jì)法統(tǒng)計(jì)了1971-2020年韶關(guān)煙草生長(zhǎng)期內(nèi)的降水，氣溫和SPEI 指數(shù)的年際變化率，如圖6所示。近50a 以來(lái)，韶關(guān)乳源和曲江的降水年際變化速率分別達(dá)到了2.20 mm/a和5.35 mm/a，降水增幅明顯；相反其余市縣的降水則呈減少趨勢(shì)。同樣SPEI 指數(shù)也只在乳源和曲江是隨年際增加。除曲江外韶關(guān)其余各地的氣溫以約0.01 ℃/a 的速率上升。這個(gè)結(jié)果表明，1971年以后，隨著降水的減少和氣溫的增加，導(dǎo)致韶關(guān)北部和南部市縣的干旱化有加重的跡象。而乳源和曲江因降水明顯增多，干旱程度有所降低。

圖6 1971-2020年韶關(guān)煙草生長(zhǎng)期內(nèi)降水、氣溫和SPEI的年際變化趨勢(shì)Fig.6 Trends in rainfall、temperature and SPEI during tobacco growing period for the period 1971-2020

MMK 檢驗(yàn)結(jié)果顯示，除曲江外韶關(guān)其他區(qū)域的氣溫變化都通過(guò)了置信度為5%的顯著性檢驗(yàn)，而曲江的降水和SPEI指數(shù)都表現(xiàn)為顯著上升。這個(gè)結(jié)果一方面印證了近半個(gè)世紀(jì)全球變暖的氣候趨勢(shì)，另一方面也表明韶關(guān)煙草生長(zhǎng)期內(nèi)的干旱化是由降水和氣溫共同作用，其中降水占主導(dǎo)因素。降水的顯著增加配合氣溫降低，會(huì)使韶關(guān)的干旱程度顯著減弱。

由于韶關(guān)的煙草種植主要集中在南雄和始興兩個(gè)產(chǎn)區(qū)，其中南雄的SPEI 指數(shù)減小趨勢(shì)更為明顯。因此，為了進(jìn)一步探究干旱指數(shù)SPEI 與煙草需水量之間的聯(lián)系，以南雄為例，采用公式(18)計(jì)算煙草生長(zhǎng)期內(nèi)的需水量。需水量和SPEI指數(shù)逐年的變化趨勢(shì)如圖7所示。從圖7 可知，南雄煙草生長(zhǎng)期需水量變化范圍為113.3~255.11 mm，多年平均值為176.11 mm，最小值發(fā)生在1993年，最大值發(fā)生在2018年，需水量呈顯著遞增趨勢(shì)，通過(guò)α= 0.05 顯著性檢驗(yàn)，遞增速率為5.98 mm/10a。運(yùn)用Pearson 相關(guān)系數(shù)法對(duì)需水量與SPEI 指數(shù)進(jìn)行相關(guān)性檢驗(yàn)，相關(guān)系數(shù)達(dá)到了-0.52，說(shuō)明兩者之間存在負(fù)相關(guān)。即當(dāng)SPEI 指數(shù)表征的干旱程度加劇時(shí)，煙草需水量越大，干旱對(duì)煙草種植的影響越嚴(yán)重。因此對(duì)于韶關(guān)南雄市，由于其干旱程度在過(guò)去50年有增強(qiáng)的趨勢(shì)，在煙草生長(zhǎng)期內(nèi)更需關(guān)注干旱災(zāi)害對(duì)當(dāng)?shù)責(zé)煵莘N植的影響。

圖7 韶關(guān)南雄市煙草生長(zhǎng)階段需水量和SPEI指數(shù)變化規(guī)律Fig.7 Changes of water requirement and SPEI during tobacco growth stages in Nanxiong

3 討論與總結(jié)

本文選用韶關(guān)市8 個(gè)國(guó)家氣象站點(diǎn)1971-2020年逐日氣象數(shù)據(jù)，利用SPEI 干旱指數(shù)對(duì)韶關(guān)市的干旱時(shí)空變化特征進(jìn)行了分析，并使用作物系數(shù)法探究了煙草生長(zhǎng)期內(nèi)韶關(guān)南雄煙草需水量與SPEI指數(shù)之間的關(guān)系。得出以下結(jié)論：

(1)韶關(guān)市SPEI指數(shù)和干旱站次比的變化趨勢(shì)表明，1971-2020年韶關(guān)市年干旱影響范圍在不斷擴(kuò)大。在進(jìn)入2000s后干旱化有增強(qiáng)的趨勢(shì)，年內(nèi)發(fā)生全域性干旱和重旱的風(fēng)險(xiǎn)明顯增大。韶關(guān)春季和夏季都易發(fā)生全域性干旱事件，另外需特別關(guān)注韶關(guān)市在春、冬兩季發(fā)生年際間旱澇急轉(zhuǎn)的情況。

(2)干旱頻率的時(shí)空變化表明，近50a韶關(guān)市的年干旱事件主要以中度干旱為主，發(fā)生次數(shù)無(wú)明顯增多。年干旱發(fā)生頻率總體上呈現(xiàn)“北多南少”的分布規(guī)律。韶關(guān)四季干旱發(fā)生次數(shù)基本相同，其中曲江是季節(jié)尺度干旱事件易發(fā)生區(qū)域。在春季和冬季韶關(guān)需關(guān)注極端干旱事件的發(fā)生。

(3)MMK 趨勢(shì)檢驗(yàn)表明，煙草生長(zhǎng)發(fā)育期間（3~6月），自1971年以來(lái)韶關(guān)除曲江外其余市縣的氣溫都為顯著增加，而曲江的降水和SPEI 指數(shù)都顯示顯著上升。因此韶關(guān)的干旱化是由降水和氣溫共同作用，其中降水占主導(dǎo)因素。降水的顯著增加配合氣溫降低，會(huì)使韶關(guān)的干旱程度顯著減弱。

(4)1971-2020年，韶關(guān)南雄煙草生長(zhǎng)期內(nèi)的需水量以5.98 mm/10a 的速率顯著遞增，與SPEI 指數(shù)呈負(fù)相關(guān)。由于過(guò)去50a 間南雄的干旱程度有增強(qiáng)的趨勢(shì)，因此更需關(guān)注南雄的干旱災(zāi)害對(duì)當(dāng)?shù)責(zé)煵莘N植的影響。

本文對(duì)韶關(guān)近50a的干旱和煙草需水量演變特征進(jìn)行了分析，但研究仍存在一些不足之處。如引起干旱的因素較多，人類活動(dòng)加劇和城市群的擴(kuò)張等要素都會(huì)對(duì)干旱產(chǎn)生一定影響。本文僅從氣象要素角度考慮了其對(duì)干旱化的貢獻(xiàn)，結(jié)果較為片面，仍有待進(jìn)一步探究。干旱對(duì)農(nóng)業(yè)生產(chǎn)影響較大，所以還需繼續(xù)研究韶關(guān)南雄未來(lái)干旱化發(fā)展趨勢(shì)對(duì)當(dāng)?shù)責(zé)煵萆L(zhǎng)的定量影響，從而為當(dāng)?shù)匚磥?lái)的煙草種植提供更細(xì)致的科學(xué)指導(dǎo)。