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

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

遼寧省地處歐亞大陸東岸，中緯度地區(qū)。由于受長白山脈的遼東半島和山東半島山地相隔，東南濕潤季風(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]，從空間角度來看遼西地區(qū)是遼寧省干旱災(zāi)害的頻發(fā)、重發(fā)區(qū)；從時間角度上，遼寧省春旱發(fā)生頻率遠(yuǎn)遠(yuǎn)大于夏旱和秋旱，春旱對農(nóng)業(yè)的影響也最為嚴(yán)重。因此，掌握遼西地區(qū)春季農(nóng)業(yè)干旱的發(fā)生規(guī)律，做好春旱預(yù)測服務(wù)工作，對于預(yù)防和減輕春旱危害以及合理安排春耕生產(chǎn)具有重要意義。

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

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

1 資料與方法

1.1 資 料

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

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

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

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

1.2 方 法

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

由于土壤相對濕度與氣象因子、地溫等數(shù)據(jù)單位相差較大，因此利用以下公式對本文涉及的數(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個子自變量x1，x2，…，xm的變化而變化，且有線性關(guān)系：

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

(3)

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

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

(4)

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

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

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

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

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

2 結(jié)果與分析

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

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

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

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

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

2.2.1 預(yù)測因子分析

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

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

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

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

注:**表示通過α=0.01的顯著性檢驗，*表示通過α=0.05的顯著性檢驗。

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

春季土壤相對濕度與同期平均氣溫相關(guān)性不顯著。一方面因為北方地區(qū)春季回暖慢，錦州地區(qū)土壤化通時間一般在3月中旬至4月中旬，對于還沒有完全解凍的土壤，氣溫并不能造成明顯的蒸發(fā)失墑現(xiàn)象；另一原因是，錦州地區(qū)春季氣溫相對較低且年際變化較小，不會對土壤造成明顯的失墑影響，所以不能作為土壤相對濕度變化的影響因子。

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

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

2.2.2 預(yù)測模型建立

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

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

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

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

0.222X3+0.446X4

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

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

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

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

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

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

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

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

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

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

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

3 結(jié) 語

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

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

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

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

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