巢湖地區(qū)土壤相對濕度氣象預(yù)測模型研究

侍永樂 章超 靳青春

摘要?[目的]描述氣象影響因子與土壤相對濕度的關(guān)系，預(yù)測巢湖地區(qū)逐旬10～100 cm土壤相對濕度，定量評估土壤旱澇情況。[方法]通過分析巢湖地區(qū)各層土壤相對濕度逐月分布特征，討論各氣象影響因子與土壤相對濕度的相關(guān)關(guān)系，并分別基于多元線性回歸和BP神經(jīng)網(wǎng)絡(luò)方法建立預(yù)測模型。[結(jié)果]氣象因子和土壤相對濕度存在相關(guān)關(guān)系，相關(guān)系數(shù)從大到小分別為蒸發(fā)量、氣溫、日照時(shí)數(shù)、相對濕度和降水量。隨著土層的深入，兩者的相關(guān)系數(shù)和顯著程度呈下降趨勢。[結(jié)論]2種方法均能較好地預(yù)測出巢湖地區(qū)各層土壤相對濕度逐旬動(dòng)態(tài)變化，可為了解巢湖地區(qū)秈稻主生長周期內(nèi)的土壤相對濕度變化提供依據(jù)，為該地區(qū)作物旱澇評估提供參考。

關(guān)鍵詞?巢湖;土壤相對濕度;多元線性回歸;BP神經(jīng)網(wǎng)絡(luò)

中圖分類號?P426.1?文獻(xiàn)標(biāo)識碼?A?文章編號?0517-6611（2021）01-0218-05

doi：10.3969/j.issn.0517-6611.2021.01.058

Abstract?[Objective]In order to describe the relationship between meteorological factors and soil relative humidity， predict the soil relative humidity of 10-100 cm of ten days in Chaohu Region and evaluate the soil drought and flood quantitatively. [Method]The relationship between meteorological factors and soil relative humidity was discussed by analyzing the monthly distribution characteristics of soil relative humidity in Chaohu Region， and the prediction models were established based on multiple linear regression and BP neural network. [Result]There is a correlation between meteorological factors and soil relative humidity. The correlation coefficients from large to small are evaporation， temperature， sunshine duration， relative humidity and precipitation. With the depth of soil layer， the correlation coefficient and significance of the two have a downward trend. [Conclusion]The two methods can predict the dynamic change of soil relative humidity in every layer of Chaohu Region， and provide the basis for the change of soil relative humidity in the main growth cycle of Indica Rice in Chaohu Region， and provide a reference for crop drought and flood assessment in this area.

Key words?Chaohu;Relative soil moisture;Multiple linear regression;BP neural network

土壤濕度即土壤水分，是氣候系統(tǒng)關(guān)鍵變量之一，作為陸面過程中一個(gè)重要的物理量，是影響作物生長、發(fā)育及產(chǎn)量的重要因素[1]，它通過改變地表反照率、熱容量、地表蒸發(fā)、植被生長狀況等途徑影響著氣候的變化[2]。土壤相對濕度是土壤水分狀況的綜合體現(xiàn)[3]，大氣降水、灌溉、地下水等有助于土壤相對濕度增加，而土壤蒸散、植物吸收、蒸騰、土壤水分滲漏和徑流則加快了土壤相對濕度降低[4]。關(guān)于土壤相對濕度的預(yù)測模型研究，大多學(xué)者采用降水量、土壤含水量、作物旱情等指標(biāo)來研究[5-6]，運(yùn)用模擬、替代、實(shí)測資料從空間分布、季節(jié)變化、垂直分布等不同層面、不同尺度上分析了土壤濕度變化特征及其影響因素。李輯等[7]通過研究土壤相對濕度與氣象因子、環(huán)流指數(shù)等的關(guān)系，建立月尺度的土壤相對濕度動(dòng)態(tài)預(yù)測模型用于預(yù)報(bào)農(nóng)業(yè)干旱;左志燕等[2，8-9]分析了我國土壤濕度的垂直分布、時(shí)空變化、季節(jié)循環(huán)特征、年際變化及其與氣候變率的關(guān)系;李潤春等[10]運(yùn)用統(tǒng)計(jì)學(xué)方法得到不同土層的逐旬土壤相對濕度預(yù)報(bào)模型，預(yù)測效果較好。

巢湖流域位于安徽省中部，東南瀕臨長江，北部依江淮分水嶺，屬長江下游左岸水系。流域均屬副熱帶季風(fēng)氣候區(qū)，氣候溫和濕潤[11]，降水分布不均。21世紀(jì)以來，頻繁的旱澇災(zāi)害給農(nóng)業(yè)生產(chǎn)帶來嚴(yán)重的影響，開展農(nóng)業(yè)旱澇監(jiān)測、預(yù)報(bào)和預(yù)警工作成為農(nóng)業(yè)氣象服務(wù)的熱點(diǎn)[12-13]。巢湖流域的土壤質(zhì)地以壤土和黏壤土為主，種植的作物多為常規(guī)秈稻，在水稻田水、肥、氣、熱等四大因素中，水為矛盾的主要方面。因此，了解巢湖流域土壤相對濕度發(fā)生規(guī)律，對預(yù)防和減輕旱澇危害及合理安排農(nóng)業(yè)生產(chǎn)具有重要意義。

國內(nèi)已有的土壤濕度統(tǒng)計(jì)學(xué)方法預(yù)報(bào)研究對多層土壤濕度預(yù)測方法比對不夠全面。筆者以巢湖流域?yàn)檠芯繀^(qū)域，選取巢湖作為代表站，分析巢湖流域土壤相對濕度分布特征，并利用相關(guān)性分析各個(gè)氣象影響因子與土壤相對濕度的關(guān)系，找出影響土壤相對濕度的主要影響因子，再基于上一旬土壤相對濕度及氣象因子數(shù)據(jù)，采用逐步線性回歸、BP神經(jīng)網(wǎng)絡(luò)2種方法預(yù)測模擬本旬10～100 cm土壤相對濕度，旨在為該地區(qū)有效應(yīng)對旱澇的發(fā)生、合理安排灌溉提供科學(xué)依據(jù)。

1?資料與方法

1.1?資料來源

研究所用的氣象資料以及土壤水分資料來源于巢湖市國家級基本氣象站及其土壤水分站的近4年（2014年11月7日—2018年10月31日）觀測數(shù)據(jù)，氣象站觀測變量包括氣溫、降水量、日照時(shí)數(shù)、相對濕度等，土壤水分站觀測變量包括體積含水率、相對濕度、重量含水率、有效水分貯存量，土壤觀測縱深為10～100 cm，共8層，埋深分別為10、20、30、40、50、60、80、100 cm。數(shù)據(jù)來源于中國氣象局。

1.2?數(shù)據(jù)處理

為方便分析討論，將近4年氣象和土壤逐小時(shí)觀測資料按需換算為旬均值和月均值，計(jì)算土壤相對濕度、氣溫、相對濕度和日照時(shí)數(shù)的平均值，計(jì)算蒸發(fā)量和降水量的累計(jì)值。因?yàn)榭紤]到不同埋深的土壤水分傳感器觀測到的土壤相對濕度具有較強(qiáng)的垂直連續(xù)性，所以對逐層土壤相對濕度數(shù)據(jù)進(jìn)行垂直一致性檢驗(yàn)，以保證土壤水分?jǐn)?shù)據(jù)準(zhǔn)確可用。對于10～100 cm的土壤觀測縱深，通過設(shè)置差異閾值，比較相鄰兩層土壤相對濕度觀測數(shù)值差值，來判斷并剔除超過閾值范圍的異常數(shù)據(jù)，規(guī)則如下：

①-80%<10～50 cm各層土壤相對濕度與上層差值≤90%;

②-65% <50～100 cm各層土壤相對濕度與上層差值≤70%。

2?結(jié)果與分析

2.1?土壤相對濕度時(shí)間變化特征

由圖1可知，巢湖地區(qū)土壤相對濕度的時(shí)間分布存在較大波動(dòng)，總體上數(shù)值變化趨勢呈“單谷”分布。1—5月各層土壤相對濕度變化較平緩，6月起隨著日照時(shí)數(shù)和氣溫逐漸升高，導(dǎo)致各層土壤相對濕度呈下降趨勢，于9月降至全年最低值，10月又開始回升?？傮w上來說，各層土壤相對濕度與其觀測的土壤深度成正比，淺層土壤相對濕度最低，土壤深度越深，土壤濕度逐漸升高;表層的土壤相對濕度受氣象條件影響明顯，其中10 cm土壤相對濕度受氣象條件影響尤為顯著。

2.2?土壤相對濕度與氣象因子的相關(guān)系數(shù)

研究表明，土壤相對濕度與降水、氣溫、蒸發(fā)、日照、灌溉條件、土壤墑情等因素有關(guān)[14]。結(jié)合前人研究成果，針對氣象因子的影響，分析巢湖土壤相對濕度與氣溫、相對濕度、降水、蒸發(fā)和日照的相關(guān)關(guān)系，選取影響土壤相對濕度的主要?dú)庀笠蜃樱瑸榻㈩A(yù)測模型做準(zhǔn)備。

由表1可知，各氣象因子中，日照時(shí)數(shù)、氣溫、蒸發(fā)量與土壤相對濕度的相關(guān)性更好，降水量和相對濕度次之;從土壤觀測縱深方面看，10 cm土壤相對濕度與氣象因子的相關(guān)性最高，與氣溫、相對濕度、降水量、蒸發(fā)量和日照時(shí)數(shù)的相關(guān)系數(shù)分別為-0.33、0.25、0.27、-0.43和-0.37。隨著觀測埋深的增加，逐層土壤相對濕度和氣象因子的相關(guān)關(guān)系總體呈波動(dòng)減弱趨勢。另一方面，土壤相對濕度主要與氣溫、蒸發(fā)量和日照時(shí)數(shù)呈負(fù)相關(guān)關(guān)系，表明氣溫升高、日照時(shí)數(shù)變長以及地表水分蒸發(fā)量的提高，將一定程度地降低土壤含水量，造成土壤相對濕度下降;而空氣相對濕度與地表降水量的增加，有助于提升土壤含水量，與土壤相對濕度呈正相關(guān)。對于小于20 cm深度的土壤相對濕度，與各氣象因子的相關(guān)關(guān)系均通過了顯著水平檢驗(yàn)（α=0.05）;對于大于20 cm的土壤相對濕度，與各氣象因子的相關(guān)關(guān)系顯著性水平存在減弱趨勢，其中30 cm土壤深度僅降水量通過α=0.05檢驗(yàn)，表明表層土壤相對濕度與氣象因子相關(guān)關(guān)系更顯著。

通過比較各氣象因子與土壤相對濕度的相關(guān)關(guān)系，選取氣溫、相對濕度、降水量、蒸發(fā)量和日照時(shí)數(shù)作為影響因子，同時(shí)因土壤水分變化存在連續(xù)性和滯后性，故選擇上一時(shí)次對應(yīng)土壤層的相對濕度作為自相關(guān)因子，綜合上述因子來建立關(guān)于土壤相對濕度的統(tǒng)計(jì)預(yù)測模型。

2.3?預(yù)測模型的建立

2.3.1?多元線性回歸法。

為建立巢湖地區(qū)不同土壤深度的相對濕度預(yù)測模型，結(jié)合氣象因子以及土壤相對濕度的相關(guān)關(guān)系，選取上一旬的氣溫、相對濕度、降水量、蒸發(fā)量、日照時(shí)數(shù)以及對應(yīng)層土壤相對濕度作為預(yù)測因子，采用基于最小二乘法的多元線性回歸方法，對當(dāng)旬逐層土壤相對濕度進(jìn)行模擬預(yù)測，選取近3年氣象和土壤觀測資料旬值（2014年11月—2017年8月）作為輸入數(shù)據(jù)，近1年旬值（2017年9月—2018年10月）作為預(yù)測結(jié)果驗(yàn)證，輸入輸出的樣本數(shù)分別為99和40，具體多元回歸系數(shù)見表2。

2.3.2?人工神經(jīng)網(wǎng)絡(luò)法。

因?yàn)闅庀笠蜃雍屯寥老鄬穸戎g的影響關(guān)系相當(dāng)復(fù)雜且存在非線性特征，所以采用基于人工神經(jīng)網(wǎng)絡(luò)的方法建立巢湖地區(qū)不同土壤深度的相對濕度預(yù)測模型。為比較2種預(yù)測方法，采用控制變量法，將其結(jié)果與多元線性回歸方法進(jìn)行對比，控制輸入和輸出變量、樣本數(shù)以及驗(yàn)證數(shù)據(jù)與多元線性回歸方法相同。模型建立過程中，根據(jù)BP神經(jīng)網(wǎng)絡(luò)需要，對輸入數(shù)據(jù)進(jìn)行歸一化處理，建立前饋網(wǎng)絡(luò)對象，設(shè)置的參數(shù)包括訓(xùn)練次數(shù)、學(xué)習(xí)率以及訓(xùn)練目標(biāo)，數(shù)值分別為100、0.1以及0.000 04，進(jìn)一步對神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練，計(jì)算仿真并通過反歸一化輸出基于BP神經(jīng)網(wǎng)絡(luò)方法對逐層土壤相對濕度的預(yù)測結(jié)果。

2.4?結(jié)果檢驗(yàn)

為驗(yàn)證對各層土壤相對濕度的預(yù)測準(zhǔn)確性，選擇2017年9月—2018年10月的氣象和土壤水分觀測資料作為結(jié)果檢驗(yàn)，基于絕對誤差百分率來評估模型預(yù)測結(jié)果，分別從2種預(yù)測方法的總體結(jié)果、各層土壤相對濕度的逐月預(yù)測分布以及預(yù)測結(jié)果較觀測值的離散情況，綜合比較多元線性回歸和BP神經(jīng)網(wǎng)絡(luò)方法在巢湖地區(qū)土壤相對濕度預(yù)測中的結(jié)果表現(xiàn)。

由表3可知，多元線性回歸和BP神經(jīng)網(wǎng)絡(luò)對0～100 cm平均土壤相對濕度預(yù)測值的絕對誤差百分率分別為4.31%和4.87%，10 cm土層相對濕度受氣象等其他要素的影響導(dǎo)致變化的趨勢較復(fù)雜，預(yù)測誤差最大，分別為17.61%和25.05%，其余土層的誤差均小于10 cm土層。除40 cm土層多元線性回歸法外，總體的絕對誤差在5%以內(nèi)。比較2種預(yù)測方法，除了40、50、100 cm土層的土壤相對濕度預(yù)測絕對誤差是BP神經(jīng)網(wǎng)絡(luò)法小于多元線性回歸法，其余各層是多元線性回歸的絕對誤差小。表明在建模預(yù)測巢湖地區(qū)各層土壤相對濕度的情況下，2種方法結(jié)果接近，均能較好地預(yù)測出土壤相對濕度的動(dòng)態(tài)變化?？傮w而言，多元線性回歸在大多數(shù)土層的預(yù)測結(jié)果絕對誤差百分率較低，整體上多元回歸方法優(yōu)于BP神經(jīng)網(wǎng)絡(luò)方法。

由表4可知，從季節(jié)來看，夏季溫度升高，降雨天氣過程頻發(fā)，土壤含水量波動(dòng)劇烈，導(dǎo)致各層土壤相對濕度的最大預(yù)測絕對誤差百分率較大;秋季受臺風(fēng)、連陰雨及秋旱等影響，各層土壤相對濕度的最大預(yù)測絕對誤差百分率僅次于夏季。夏季多元線性回歸法和BP神經(jīng)網(wǎng)絡(luò)法預(yù)測結(jié)果均值分別為11.96%和14.72%;秋季次之，分別為9.11%和9.94%;春季和冬季的預(yù)測誤差較小，春季分別為3.29%和3.76%，冬季分別為4.16%和4.04%。由于巢湖地區(qū)秈稻的主要生長發(fā)育階段在6—8月，收獲時(shí)節(jié)在10月前后（秋季），因此上述預(yù)測模型能夠較好地預(yù)測巢湖地區(qū)土壤相對濕度的時(shí)間變化規(guī)律，為水稻等作物的旱澇評估提供參考。

采用線性回歸和BP神經(jīng)網(wǎng)絡(luò)法分別預(yù)測各層土壤相對濕度月分布，結(jié)果見圖2。由圖2可知，2種方法均能較好地預(yù)測出各深度土壤濕度的逐月變化趨勢，模擬效果良好，兩者誤差標(biāo)準(zhǔn)差分別為5.90%和6.96%;設(shè)置絕對誤差百分率閾值范圍分別為0～5%、5%～20%、20%～50%和大于50%，其中多元線性回歸方法對應(yīng)誤差閾值范圍分別為57.81%、32.19%、8.75%和1.25%，BP神經(jīng)網(wǎng)絡(luò)對應(yīng)為59.06%、29.69%、9.38%和1.19%。對比可知，BP神經(jīng)網(wǎng)絡(luò)計(jì)算出的誤差百分率低于20%和高于50%的結(jié)果比率值均優(yōu)于多元線性回歸結(jié)果，表明BP神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果的誤差離散比率更小，可為定量評價(jià)土壤墑情預(yù)測方法提供數(shù)據(jù)支撐。

由圖3可知，10 cm深度的2種方法預(yù)測值散度大，其他各層的散度對應(yīng)良好;20和30 cm的土壤相對濕度變化范圍不大，其余各層變化范圍較大;從60 cm土層的2種方法預(yù)測值來看，在這一層用BP神經(jīng)網(wǎng)絡(luò)法對應(yīng)不好，跳變點(diǎn)較用多元線性回歸法多。

3?結(jié)論

（1）巢湖地區(qū)土壤相對濕度的逐月分布存在較大的波動(dòng)，總體上數(shù)值變化趨勢呈“單谷”分布，1—5月各層土壤相對濕度變化較平緩，自6月開始各層土壤相對濕度呈下降趨勢，于9月降至全年最低值。

（2）相關(guān)系數(shù)能夠較好地反映各層土壤相對濕度和氣象因子之間的關(guān)系，其中日照、氣溫和蒸發(fā)與土壤相對濕度的相關(guān)性更好，降水和相對濕度次之。表層土壤相對濕度受氣象因子影響最為明顯，10 cm土壤相對濕度與氣象因子的相關(guān)性最高，對于大于20 cm的土壤相對濕度，隨著土層的深入，土壤相對濕度和氣象因子的相關(guān)系數(shù)和顯著程度呈下降趨勢。

（3）關(guān)于建立巢湖地區(qū)各層土壤相對濕度的逐旬預(yù)測模型，一方面，基于氣象因子和土壤相對濕度的多元線性回歸預(yù)測方法，可定量闡述氣象因子與土壤相對濕度的權(quán)重值。另一方面，氣象因子和各層土壤相對濕度之間的影響相當(dāng)復(fù)雜，存在非線性關(guān)聯(lián)，BP神經(jīng)網(wǎng)絡(luò)方法可建立基于人工神經(jīng)網(wǎng)絡(luò)的非線性關(guān)聯(lián)模型，預(yù)測土壤相對濕度。

（4）關(guān)于多元線性回歸和BP神經(jīng)網(wǎng)絡(luò)預(yù)測巢湖地區(qū)逐旬各層土壤相對濕度的預(yù)測結(jié)果檢驗(yàn)，總體上看，2種方法均能較好地預(yù)測出各層土壤相對濕度的動(dòng)態(tài)變化趨勢，其中多元線性回歸平均絕對誤差百分率較小。從各土壤層結(jié)果來看，10 cm土壤相對濕度受氣象等其他要素的影響導(dǎo)致變化的趨勢較復(fù)雜，模擬誤差最大，除40 cm土層多元線性回歸法外，總體的絕對誤差在5%以內(nèi)。從季節(jié)來看，夏季的預(yù)測誤差最大，秋季、冬季次之，春季預(yù)測誤差最小;從兩者預(yù)測穩(wěn)定性來看，BP神經(jīng)網(wǎng)絡(luò)預(yù)測誤差離散比率更小。上述模擬預(yù)測結(jié)果可對應(yīng)秈稻的主要生長周期，能夠較好地預(yù)測巢湖地區(qū)秈稻主生長周期內(nèi)的土壤相對濕度變化，為土壤定量旱澇評估提供一定的參考。

參考文獻(xiàn)

[1]孫倩倩.東北地區(qū)土壤濕度的時(shí)空分布特征及其預(yù)報(bào)方法的研究[D].北京：中國氣象科學(xué)研究院，2013.

[2]左志燕，張人禾.中國東部春季土壤濕度的時(shí)空變化特征[J].中國科學(xué)（D輯：地球科學(xué)），2008，38（11）：1428-1437.

[3]CHAHINE M T.The hydrological cycle and its influence on climate[J].Nature，1992，359（6394）：373-380.

[4]吳娜，石培基，朱國鋒，等.長江流域土壤相對濕度變化及其影響因素[J].長沙流域資源與環(huán)境，2017，26（7）：1001-1010.

[5]丘寶劍，盧其堯.農(nóng)業(yè)氣候條件及其指標(biāo)[M].北京：測繪出版社，1990.

[6]王密俠，馬成軍，蔡煥杰.農(nóng)業(yè)干旱指標(biāo)研究與進(jìn)展[J].干旱地區(qū)農(nóng)業(yè)研究，1998，16（3）：119-124.

[7]李輯，李雨鴻，胡春麗，等.遼西干旱區(qū)春播期土壤相對濕度動(dòng)態(tài)預(yù)測模型研究[J].土壤通報(bào)，2014，45（4）：830-834.

[8]張秀芝，吳迅英，何金海.中國土壤濕度的垂直變化特征[J].氣象學(xué)報(bào)，2004，62（1）：51-61.

[9]馬柱國，魏和林，符淙斌.土壤濕度與氣候變化關(guān)系的研究進(jìn)展與展望[J].地球科學(xué)進(jìn)展，1999，14（3）：299-305.

[10]李潤春，張秀芝，高俊壽，等.山西省代表站不同土層逐旬土壤相對濕度預(yù)報(bào)模型[J].中國農(nóng)業(yè)氣象，2009，30（3）：354-359.

[11]陳實(shí)，高超，黃銀蘭.不同季節(jié)劃分尺度下巢湖流域氣候變化趨勢分析[J].長江流域資源與環(huán)境，2013，22（5）：582-587.

[12]王素萍，張存杰，宋連春，等.多尺度氣象干旱與土壤相對濕度的關(guān)系研究[J].冰川凍土，2013，35（4）：865-873.

[13]袁媛，王心源，李祥，等.巢湖流域旱澇時(shí)空特性分析[J].災(zāi)害學(xué)，2007，22（2）：97-100.