王 瑩，張曉月，張 琪，焦 敏，黃 巖，馬 林，王若男

(1.遼寧省氣象科學(xué)研究所，沈陽(yáng) 110166；2.遼寧省氣象裝備保障中心，沈陽(yáng) 110166)

遼寧省地處歐亞大陸東岸，中緯度地區(qū)。由于受長(zhǎng)白山脈的遼東半島和山東半島山地相隔，東南濕潤(rùn)季風(fēng)不易深入，遼寧東部和中部地區(qū)呈溫帶大陸季風(fēng)氣候。遼寧西部與蒙古高原相接，受西風(fēng)環(huán)流和西伯利亞大陸性氣候影響，西部地區(qū)呈大陸性氣候。遼寧省是東北地區(qū)重要的糧食主產(chǎn)區(qū)，受氣候和復(fù)雜地貌影響，干旱災(zāi)害頻繁發(fā)生制約著糧食產(chǎn)業(yè)的發(fā)展[1,2]。有關(guān)研究表明[3-5]，從空間角度來(lái)看遼西地區(qū)是遼寧省干旱災(zāi)害的頻發(fā)、重發(fā)區(qū)；從時(shí)間角度上，遼寧省春旱發(fā)生頻率遠(yuǎn)遠(yuǎn)大于夏旱和秋旱，春旱對(duì)農(nóng)業(yè)的影響也最為嚴(yán)重。因此，掌握遼西地區(qū)春季農(nóng)業(yè)干旱的發(fā)生規(guī)律，做好春旱預(yù)測(cè)服務(wù)工作，對(duì)于預(yù)防和減輕春旱危害以及合理安排春耕生產(chǎn)具有重要意義。

進(jìn)入21世紀(jì)，由于頻繁的干旱災(zāi)害給農(nóng)業(yè)生產(chǎn)帶來(lái)嚴(yán)重影響，開展農(nóng)業(yè)干旱監(jiān)測(cè)、預(yù)報(bào)及預(yù)警研究成為農(nóng)業(yè)氣象服務(wù)與研究的熱點(diǎn)[6]。學(xué)者一般采用降水量指標(biāo)、土壤含水量指標(biāo)、作物旱情指標(biāo)及綜合類指標(biāo)等來(lái)做評(píng)價(jià)和研究[7-9]，考慮到農(nóng)作物生長(zhǎng)發(fā)育需要的水分絕大部分來(lái)自于土壤，因此以土壤濕度作為評(píng)價(jià)農(nóng)業(yè)氣象干旱的指標(biāo)，可以直觀的表現(xiàn)出農(nóng)田干旱程度。李輯等[10]通過(guò)研究土壤相對(duì)濕度與氣象因子、環(huán)流指數(shù)等的關(guān)系，建立了月尺度的土壤相對(duì)濕度動(dòng)態(tài)預(yù)測(cè)模型，預(yù)報(bào)農(nóng)業(yè)干旱。李潤(rùn)春等[11]分析了土壤相對(duì)濕度變化的影響因子，運(yùn)用統(tǒng)計(jì)學(xué)方法得到不同土層逐旬土壤相對(duì)濕度預(yù)報(bào)模型，預(yù)報(bào)效果較好。還有一些學(xué)者[12-13]以土壤水分平衡原理為出發(fā)點(diǎn)，計(jì)算土壤水分變化的各影響參數(shù)，預(yù)報(bào)農(nóng)田土壤水分含量。然而，季節(jié)尺度上的土壤相對(duì)濕度的預(yù)測(cè)方法卻鮮見報(bào)道，相應(yīng)地關(guān)于季節(jié)尺度上的定量化農(nóng)業(yè)干旱預(yù)測(cè)也少有研究。

本文通過(guò)分析遼西地區(qū)土壤相對(duì)濕度變化特征，揭示遼西地區(qū)春季農(nóng)業(yè)干旱發(fā)生發(fā)展規(guī)律。選取遼西地區(qū)具有完整觀測(cè)資料的錦州作為代表站，利用相關(guān)性分析各影響要素與土壤相對(duì)濕度的關(guān)系，從中尋找影響土壤相對(duì)濕度的主要因子，建立錦州春季(3-5月)土壤相對(duì)濕度預(yù)測(cè)模型，實(shí)現(xiàn)定量化的春季農(nóng)業(yè)干旱預(yù)測(cè)，為該地區(qū)有效應(yīng)對(duì)春旱發(fā)生、合理安排灌溉提供科學(xué)依據(jù)。

1 資料與方法

1.1 資 料

選取遼西地區(qū)具有完整資料的錦州、朝陽(yáng)、阜新、建平為代表站點(diǎn)，所用資料包括土壤相對(duì)濕度資料、地溫資料和氣象資料。

(1)土壤相對(duì)濕度。土壤相對(duì)濕度采用人工土壤水分觀測(cè)站1980-2015年每旬逢3、8觀測(cè)數(shù)據(jù)。

(2)地溫。地溫采用遼寧省氣象局提供的1980-2015年80 cm地溫的歷史資料。

(3)氣象資料。氣象資料為1980-2015 年逐日平均氣溫和日降水量的觀測(cè)數(shù)據(jù)。

1.2 方 法

1.2.1 標(biāo)準(zhǔn)差標(biāo)準(zhǔn)化[14]

由于土壤相對(duì)濕度與氣象因子、地溫等數(shù)據(jù)單位相差較大，因此利用以下公式對(duì)本文涉及的數(shù)據(jù)資料，為相關(guān)性分析及回歸模型的建立做數(shù)據(jù)準(zhǔn)備。

序列Xi(i=1,2,…,n)的標(biāo)準(zhǔn)差為：

(1)

距平變量除以標(biāo)準(zhǔn)差后得到標(biāo)準(zhǔn)化變量序列，即為標(biāo)準(zhǔn)差標(biāo)準(zhǔn)化。

(2)

1.2.2 多元回歸方程[14]

設(shè)因變量y隨m個(gè)子自變量x1，x2，…，xm的變化而變化，且有線性關(guān)系：

y=β0+β1x1+β2x2+…+βmxm+ε

(3)

有n組觀測(cè)資料 (x1t，x2t，…,xmt;yt)它們滿足回歸方程，即滿足下面n個(gè)等式：

yt=β0+β1x1t+β2x2t+…+βmxmt+εt

(4)

t=1,2，…，n，εt是隨機(jī)誤差，它仍需滿足一元回歸中所述的三個(gè)假定，即“誤差εt相互獨(dú)立地服從正態(tài)分布N(0σ2)”。

本文采用此方法，建立春季(3-5月)土壤相對(duì)濕度預(yù)測(cè)模型，實(shí)現(xiàn)定量化的春季農(nóng)業(yè)干旱預(yù)測(cè)。

1.2.3 農(nóng)業(yè)干旱等級(jí)劃分

根據(jù)遼寧省多年研究和業(yè)務(wù)實(shí)際應(yīng)用，以土壤相對(duì)濕度劃分農(nóng)業(yè)干旱等級(jí)(表1)，可以直觀的表現(xiàn)旱澇狀況。

表1 土壤相對(duì)濕度劃分農(nóng)業(yè)干旱等級(jí)Tab.1 Different grades of agricultural drought disaster by soil relative humidity

2 結(jié)果與分析

2.1 遼西地區(qū)春季土壤相對(duì)濕度變化特征

分析遼西地區(qū)春季土壤相對(duì)濕度變化的時(shí)間特征(如圖1)，可以看出：遼西地區(qū)土壤相對(duì)濕度年際變化存在很大的波動(dòng)性。總體上，土壤相對(duì)濕度呈波動(dòng)下降趨勢(shì)，下降幅度較小，傾向率為-0.23%/a。結(jié)合表1發(fā)現(xiàn)，2000年以后伴隨著土壤相對(duì)濕度減少，農(nóng)業(yè)干旱年份頻現(xiàn)。土壤濕度的時(shí)間變化規(guī)律揭示了遼西地區(qū)春季農(nóng)業(yè)干旱越來(lái)越重的現(xiàn)象，這主要與遼西地區(qū)春季氣候變化特征相關(guān)。孫鳳華等[15]研究發(fā)現(xiàn)，遼西地區(qū)總體氣候變化趨勢(shì)為氣溫升高、降水減少。暖干的氣候有利于干旱趨勢(shì)的發(fā)展，加之春季風(fēng)大，土壤水分蒸發(fā)加快，從而出現(xiàn)了遼西地區(qū)春季土壤相對(duì)濕度減少的現(xiàn)象，農(nóng)業(yè)干旱增加。

圖1 遼西地區(qū)春季土壤相對(duì)濕度年際變化曲線Fig.1 Inter-annual variation of soil relative humidity in spring in western Liaoning Province

2.2 土壤相對(duì)濕度預(yù)測(cè)方法研究

本文選取遼西地區(qū)具有完整觀測(cè)資料的錦州作為代表站點(diǎn)，建立春季(3-5月)土壤相對(duì)濕度預(yù)測(cè)模型，結(jié)合農(nóng)業(yè)干旱等級(jí)劃分(表1)，實(shí)現(xiàn)春季農(nóng)業(yè)干旱定量化預(yù)測(cè)。

2.2.1 預(yù)測(cè)因子分析

有研究表明[16]，土壤相對(duì)濕度的變化不僅與降水有關(guān)，還與同期各氣象要素、前期的降水多少、土壤底墑狀況、灌溉條件、地形、土質(zhì)等諸多因素有關(guān)。本文綜合前人研究成果，通過(guò)分析遼西地區(qū)春旱發(fā)生特點(diǎn)及可能影響因子，計(jì)算了錦州地區(qū)春季土壤相對(duì)濕度與前一年秋季降水、封凍前底墑、前一年冬季80 cm地溫及同期降水、氣溫等要素之間的相關(guān)性，尋找相關(guān)顯著的因子，進(jìn)一步分析高相關(guān)因子如何影響土壤濕度，從而篩選出預(yù)測(cè)因子。

計(jì)算錦州地區(qū)春季土壤相對(duì)濕度與影響因子的相關(guān)性(表2)，從結(jié)果中可以看出，錦州地區(qū)春季土壤相對(duì)濕度與前一年秋季降水呈顯著正相關(guān)，通過(guò)α=0.01(rc=0.449)的顯著性檢驗(yàn)。這說(shuō)明，前秋降水偏多時(shí)，來(lái)年春季土壤水分含量多，土壤濕潤(rùn)減少春旱發(fā)生，有利于進(jìn)行春耕春播；反之，前秋降水偏少時(shí)，第二年春季土壤偏干，容易發(fā)生春旱。前秋降水是北方特有的土壤相對(duì)濕度影響因子，是指上一年9-11月份的降水量，在北方將其稱為封凍雨。秋季降水儲(chǔ)存在土壤里，隨著冬季到來(lái)，氣溫下降到冰點(diǎn)以下，水分被封凍在土壤里，當(dāng)?shù)诙甏杭镜絹?lái)，伴隨著氣溫升高，土壤開始解凍，便會(huì)形成明顯的返漿現(xiàn)象，土壤變得濕潤(rùn)，從而減輕春旱發(fā)生的程度。

底墑是指前一年土壤封凍前最后一次測(cè)墑的平均土壤相對(duì)濕度。本文分析了土壤相對(duì)濕度與底墑的相關(guān)性(表2)，可以看出二者之間存在明顯的正相關(guān)關(guān)系，通過(guò)了α=0.01(rc=0.449)的顯著性檢驗(yàn)。在所有前期影響因子中，底墑與春季土壤相對(duì)濕度的相關(guān)系數(shù)明顯高于其他因子，這說(shuō)明了底墑所體現(xiàn)的土壤相對(duì)濕度是最終封凍在土壤里影響第二年春季土壤干旱與否最重要的因子，春季土壤濕度與底墑的關(guān)系更為直接。底墑大，則代表封凍在土壤里的水分多，第二年春季土壤解凍，土壤相對(duì)濕潤(rùn)濕度大，不易發(fā)生春旱，反之亦然。

表2 春季土壤相對(duì)濕度與影響因子相關(guān)統(tǒng)計(jì)表Tab.2 Correlation coefficients between soil relative humidity and influencing factors

注:**表示通過(guò)α=0.01的顯著性檢驗(yàn)，*表示通過(guò)α=0.05的顯著性檢驗(yàn)。

地溫是衡量地表土壤熱能的物理量，在時(shí)間演變過(guò)程中，地溫的變化會(huì)對(duì)大氣中天氣氣候變化產(chǎn)生反饋?zhàn)饔?。深層地溫的時(shí)空演變具有大尺度特征，與大氣環(huán)流之間有很好的相互對(duì)應(yīng)關(guān)系。因此有學(xué)者[17]對(duì)深層地溫與長(zhǎng)期降水趨勢(shì)的相關(guān)性進(jìn)行研究，并將深層地溫作為降水長(zhǎng)期預(yù)測(cè)的一個(gè)因子，編制預(yù)報(bào)方程，得到了很好的預(yù)報(bào)效果。據(jù)統(tǒng)計(jì)，80 cm深處地溫振動(dòng)持續(xù)期大致為3個(gè)月左右。因此從時(shí)間尺度上的對(duì)應(yīng)來(lái)說(shuō)，前一年冬季80 cm地溫可作為第二年春季干旱的預(yù)報(bào)因子。通過(guò)計(jì)算前冬80 cm地溫與春季土壤相對(duì)濕度的相關(guān)性可知(表2)，二者之間存在較好的正相關(guān)性，通過(guò)α=0.01(rc=0.449)的顯著性檢驗(yàn)。這表明，前冬80 cm地溫越高，則第二年春季土壤相對(duì)濕度越大，春旱發(fā)生的可能性越??；反之前冬80 cm地溫越小，則發(fā)生春旱的可能性越大。

春季土壤相對(duì)濕度與同期平均氣溫相關(guān)性不顯著。一方面因?yàn)楸狈降貐^(qū)春季回暖慢，錦州地區(qū)土壤化通時(shí)間一般在3月中旬至4月中旬，對(duì)于還沒(méi)有完全解凍的土壤，氣溫并不能造成明顯的蒸發(fā)失墑現(xiàn)象；另一原因是，錦州地區(qū)春季氣溫相對(duì)較低且年際變化較小，不會(huì)對(duì)土壤造成明顯的失墑?dòng)绊懀圆荒茏鳛橥寥老鄬?duì)濕度變化的影響因子。

春季土壤相對(duì)濕度與同期降雨量具有顯著正相關(guān)，即春季降水偏多時(shí)，同期土壤相對(duì)濕度偏大，會(huì)減少春旱的發(fā)生；反之，春季降水偏少，同期土壤偏干，易發(fā)生春旱，不利于進(jìn)行春播春耕等農(nóng)事活動(dòng)。值得注意的是，同期降水與土壤濕度的正相關(guān)性小于底墑與土壤濕度的正相關(guān)性，即同期降水對(duì)土壤濕度的影響仍未超過(guò)底墑對(duì)土壤濕度的影響，這也說(shuō)明了同期降水并不能對(duì)土壤相對(duì)濕度起決定性作用，農(nóng)業(yè)干旱還是要以土壤相對(duì)濕度作為定量評(píng)估指標(biāo)才更具有實(shí)際意義，對(duì)春播農(nóng)業(yè)生產(chǎn)更具有指導(dǎo)價(jià)值，這也是本文選用土壤相對(duì)濕度作為春旱預(yù)測(cè)指標(biāo)的意義所在。

通過(guò)上述分析發(fā)現(xiàn)，春季土壤濕度的主要影響因子是底墑、前秋降水、前冬80 cm地溫和同期降水。因此，選用以上具有預(yù)測(cè)意義的4個(gè)因子，可為春旱預(yù)測(cè)提供有力的參考。

2.2.2 預(yù)測(cè)模型建立

根據(jù)以上分析，選取相關(guān)性較高的底墑、前秋降水、前冬80 cm地溫和春季降水作為預(yù)測(cè)因子(表3)。

表3 預(yù)報(bào)春季土壤相對(duì)濕度關(guān)鍵因子Tab.3 The key factors for forecasting of soil relative humidiy in spring

選取1981-2012年32各樣本，利用多元線性回歸分析得到錦州地區(qū)春季土壤相對(duì)濕度的預(yù)測(cè)方程：

Y=-3.36×10-16+0.201X1+0.203X2+

0.222X3+0.446X4

方程通過(guò)α=0.01顯著性檢驗(yàn)，復(fù)相關(guān)系數(shù)達(dá)到0.78。預(yù)測(cè)模型中，春季降水可用預(yù)報(bào)值，前冬80 cm地溫、前秋降雨量和底墑可用實(shí)測(cè)值。這種方法當(dāng)土壤相對(duì)濕度觀測(cè)資料比較系統(tǒng)完整，且資料年代較長(zhǎng)時(shí)，是比較實(shí)用的方法。模型中的因子物理意義明確，所以預(yù)報(bào)結(jié)果比較穩(wěn)定。

根據(jù)上述回歸模型得到春季土壤相對(duì)濕度預(yù)測(cè)值，找到相應(yīng)干旱等級(jí)(表1)，便可預(yù)測(cè)春旱的發(fā)生。

2.2.3 預(yù)測(cè)結(jié)果檢驗(yàn)與應(yīng)用分析

為驗(yàn)證錦州地區(qū)春季土壤相對(duì)濕度回歸模型的準(zhǔn)確性，選取未參加計(jì)算的錦州地區(qū)2013-2015年春季土壤相對(duì)濕度數(shù)據(jù)進(jìn)行驗(yàn)證。采用統(tǒng)計(jì)相對(duì)誤差的方法進(jìn)行預(yù)測(cè)結(jié)果檢驗(yàn)，令相對(duì)誤差為：

R=|Y1-Y|/Y×100%

式中:Y1為模擬的土壤相對(duì)濕度；Y為實(shí)測(cè)土壤相對(duì)濕度。根據(jù)公式得到的春季土壤相對(duì)濕度相對(duì)誤差結(jié)果見表4。結(jié)果可以看出，春季土壤相對(duì)濕度預(yù)測(cè)結(jié)果較好，準(zhǔn)確率達(dá)到93.3%以上，模擬精度較高。

應(yīng)用春季土壤相對(duì)濕度預(yù)測(cè)模型對(duì)錦州地區(qū)1981-2015年春季土壤相對(duì)濕度進(jìn)行擬合，等到的數(shù)據(jù)與實(shí)測(cè)值進(jìn)行對(duì)比(圖2)，可以看出模型的擬合效果較好，擬合值與實(shí)測(cè)值相差不大。根據(jù)農(nóng)業(yè)干旱等級(jí)劃分(表1)，分析模型在預(yù)測(cè)春旱發(fā)生上的應(yīng)用效果。從整體上看，錦州地區(qū)常年春季土壤相對(duì)濕度實(shí)測(cè)值多集中在60%～90%，土壤表現(xiàn)為適宜；在這一區(qū)域內(nèi)，回歸模型得到的擬合值與實(shí)測(cè)值差別不大，同樣土壤濕度為適宜。

表4 模擬春季土壤相對(duì)濕度相對(duì)誤差Tab.4 Relative error of simulated soil relative humidiy in spring

圖2 春季土壤相對(duì)濕度模擬值與觀測(cè)值比較Fig.2 The comparison of simulation and observation of soil relative humidiy in spring in Jinzhou

應(yīng)用模型模擬春旱發(fā)生，漏報(bào)1年(1983年)、正確5年(1981、1982、1984、2003、2014年)、錯(cuò)報(bào)1年(2001年)。其中， 1983年模擬結(jié)果為土壤適宜，而實(shí)測(cè)結(jié)果為土壤輕旱，分析預(yù)報(bào)影響因子發(fā)現(xiàn)：底墑?shì)p旱，同期降水偏多；2001年模擬結(jié)果為土壤輕旱，而實(shí)測(cè)結(jié)果為土壤適宜，分析預(yù)報(bào)影響因子發(fā)現(xiàn)：底墑適宜，同期降水偏少。由此可見，相比同期降水量而言，底墑是最具有預(yù)測(cè)意義的因子，與春季土壤濕度相關(guān)尤為顯著，這同樣驗(yàn)證了前文分析中提到的底墑與春季土壤相對(duì)濕度相關(guān)性最為顯著的結(jié)論。底墑是春旱預(yù)測(cè)的主要持續(xù)性影響因子，對(duì)于春旱的發(fā)生具有指示意義。綜上所述，應(yīng)用該模型對(duì)春旱進(jìn)行預(yù)報(bào)時(shí)，要綜合考慮土壤相對(duì)濕度預(yù)測(cè)值和底墑的情況，即當(dāng)春季土壤相對(duì)濕度模擬結(jié)果在50%～60%左右，則要參考底墑情況，若底墑為旱，則考慮關(guān)注春旱的發(fā)生情況。

圖2中可以看出，對(duì)于土壤相對(duì)濕度實(shí)測(cè)值大于90%的年份(1994、2005年)，預(yù)測(cè)模型的擬合值均在85%以上。鑒于底墑對(duì)春季土壤相對(duì)濕度影響的顯著性和持續(xù)性，分析這兩年的底墑情況發(fā)現(xiàn)：1994年底墑為90%，土壤飽和；2005年底墑為88%，土壤相對(duì)濕度接近飽和。綜合分析可知，當(dāng)模擬值接近飽和時(shí)，若底墑值接近飽和或已達(dá)到飽和，則錦州地區(qū)的春季土壤飽和的可能性較大。

3 結(jié) 語(yǔ)

(1)分析遼西地區(qū)春季土壤相對(duì)濕度年際變化特征發(fā)現(xiàn)，近35年土壤相對(duì)濕度呈遞減趨勢(shì)，這揭示了遼西地區(qū)春季農(nóng)業(yè)干旱越來(lái)越重的現(xiàn)象，特別是2000年以后，春季農(nóng)業(yè)干旱年份頻現(xiàn)。

(2)遼西地區(qū)春季土壤濕度主要與底墑、前秋降水、前冬80 cm地溫和同期降水呈顯著正相關(guān)，且與底墑相關(guān)性最為顯著。

(3)選取相關(guān)性顯著的各因子作為預(yù)報(bào)因子，基于多元線性回歸建立錦州春季土壤相對(duì)濕度預(yù)報(bào)模型，復(fù)相關(guān)系數(shù)達(dá)到0.78，擬合準(zhǔn)確率在93.3%以上，預(yù)報(bào)模型模擬準(zhǔn)確度較高。根據(jù)回歸模型得到的春季土壤相對(duì)濕度預(yù)測(cè)值，結(jié)合農(nóng)業(yè)干旱等級(jí)劃分，實(shí)現(xiàn)對(duì)春旱發(fā)生的預(yù)測(cè)。

(4)鑒于底墑與土壤相對(duì)濕度的相關(guān)性最為顯著，在應(yīng)用預(yù)測(cè)模型對(duì)春季干旱進(jìn)行預(yù)報(bào)時(shí)，要綜合考慮土壤相對(duì)濕度預(yù)測(cè)值和底墑的情況，即當(dāng)春季土壤相對(duì)濕度模擬結(jié)果在50%～60%左右，則要參考底墑情況，若底墑為旱，則可能發(fā)生春旱。

(5)春季土壤墑情預(yù)測(cè)模型，作為一種定量化春旱預(yù)測(cè)方法，可以推廣到遼寧省其他站點(diǎn)和地區(qū)進(jìn)行春季干旱預(yù)測(cè)，為有效應(yīng)對(duì)春旱發(fā)生、合理安排灌溉提供科學(xué)依據(jù)。應(yīng)該注意的是，影響因子的選擇及模型的實(shí)際應(yīng)用要做本地化分析及處理。