(太原理工大學(xué)水利科學(xué)與工程學(xué)院，山西 太原 030024)

土壤水分是土壤最重要組成部分之一，是作物生長吸收水分的主要來源[1]。土壤水分具有一定的時(shí)空變異性[2]。土壤水分空間變異性受土壤水分控制因素的影響，隨研究尺度的不同而不同[3-4]。土壤水分的變異性與土壤含水率有顯著關(guān)系，土壤含水率高則呈現(xiàn)出強(qiáng)變異性，土壤含水率低則呈現(xiàn)弱變異性[5]。國內(nèi)外學(xué)者利用經(jīng)典統(tǒng)計(jì)學(xué)、地統(tǒng)計(jì)學(xué)等方法對土壤水分的時(shí)變性進(jìn)行了較多研究[6]。秦海琴等[7]應(yīng)用地統(tǒng)計(jì)學(xué)方法對白芨灘淺層土壤水分的空間變異特征進(jìn)行了研究，結(jié)果表明10cm、30cm的土層含水率屬中等變異；馮向星等[8]的研究表明，100cm以上土層，草地的土壤水分變異性最大，林地的較?。恍闲窆獾萚9]的研究表明0～60cm深土壤含水率在水平和垂直方向上均表現(xiàn)為弱變異和中等變異；李建波等[10]分析了小尺度溫室土壤含水率的空間變異性，發(fā)現(xiàn)土壤含水率的空間變異性受土壤結(jié)構(gòu)、耕地深度、植物生長以及土壤微生物活動(dòng)的影響較大。綜上，土壤含水率的空間變異性研究多限于人工土鉆取土以及中子儀定位測試的土壤含水率，目前基于數(shù)據(jù)云平臺+物聯(lián)網(wǎng)環(huán)境墑情網(wǎng)絡(luò)遠(yuǎn)程智能監(jiān)控技術(shù)條件下的土壤含水率空間變異性的研究還較少。本研究以物聯(lián)網(wǎng)環(huán)境墑情監(jiān)測系統(tǒng)在線監(jiān)測的土壤含水率為依據(jù)，對霍泉灌區(qū)小尺度區(qū)域上的土壤含水率的空間變異性進(jìn)行研究，為有限供水情況下精準(zhǔn)灌溉提供理論依據(jù)。

1 材料與方法

1.1 試區(qū)概況

試驗(yàn)區(qū)位于山西省洪洞縣廣勝寺鎮(zhèn)，區(qū)域?qū)儆谂瘻貛О敫珊荡箨懶约撅L(fēng)氣候。多年平均氣溫12.4℃，年日照時(shí)數(shù)1982.3h，無霜期為243天，年最大凍土深度為31cm，冬春寒冷干燥，夏秋暖熱多雨。多年平均降水量439.3mm，雨季多集中在5—8月，多年平均蒸發(fā)量1420.5mm，地下水埋深4～6m。灌區(qū)土壤質(zhì)地為輕壤土，土壤容重1.46g/cm3，土壤飽和含水率30.8%，田間持水率24.6%，凋萎含水率6.0%，有機(jī)質(zhì)含量1.23，全氮0.11%，有效磷19.14mg/kg，速效鉀141.82mg/kg。

2019年11月至2020年11月試驗(yàn)區(qū)降水量蒸發(fā)量及氣溫分布見圖1，測試期間平均氣溫12.79℃，最高氣溫28.40℃，最低氣溫-4.80℃?？偨邓?48.90mm，日平均降水量9.10mm，日平均蒸發(fā)量2.90mm。降水量集中在5—8月，其中，以8月的降水量最高。

圖1 霍泉灌區(qū)降水量、蒸發(fā)量及日平均氣溫分布

1.2 試驗(yàn)設(shè)計(jì)

試驗(yàn)采用兩套物聯(lián)網(wǎng)環(huán)境墑情監(jiān)測系統(tǒng)，分別安置于洪洞縣廣勝寺鎮(zhèn)東安村和三條溝試驗(yàn)區(qū)，用于在線監(jiān)測大田土壤墑情。每套系統(tǒng)有一個(gè)基站，基站控制2～3個(gè)節(jié)點(diǎn)(Ⅰ、Ⅱ代表試驗(yàn)站第一個(gè)、第二個(gè)節(jié)點(diǎn)的土壤含水率，Ⅲ、Ⅳ、Ⅴ代表三條溝村第一個(gè)、第二個(gè)、第三個(gè)節(jié)點(diǎn)的土壤含水率)。每個(gè)節(jié)點(diǎn)均有3個(gè)傳感器，其中有兩個(gè)是土壤水分傳感器，分別埋設(shè)于地下30cm、60cm深度處；有一個(gè)溫度傳感器，埋設(shè)于地下30cm深度處。用戶可通過手機(jī)或電腦上的科百云數(shù)據(jù)平臺實(shí)時(shí)查看監(jiān)測點(diǎn)土壤墑情狀況。灌區(qū)灌溉面積6733.33hm2，試區(qū)面積約80.00hm2。兩臺基站相距約有3km左右，系統(tǒng)布置見圖2。

圖2 試驗(yàn)區(qū)基站及節(jié)點(diǎn)布置示意

1.3 主要測試指標(biāo)

1.3.1 平均含水率

傳感器監(jiān)測的是30cm、60cm處逐時(shí)土壤含水率，以每天逐時(shí)的平均含水率作為逐日土壤含水率，利用式(1)計(jì)算出0～60cm的日平均含水率。時(shí)間節(jié)點(diǎn)為2019年11月12日至2020年11月13日。

θt=a1θ30+b1θ60-c1

(1)

式中：θt為0～60cm的日平均含水率；θ30為30cm處的日平均含水率；θ60為60cm處的日平均含水率。根據(jù)李泳霖等[11]的研究結(jié)果：a1為0.687，b1為0.328，c1為-0.179。

1.3.2 土壤水分空間變異系數(shù)

標(biāo)準(zhǔn)差和變異系數(shù)反映土壤含水率的空間變異特性，變異系數(shù)值越小表示土壤含水率變化越小。按式(2)、式(3)計(jì)算5個(gè)測點(diǎn)30cm、60cm處土壤含水率變化、0～60cm平均土壤含水率的標(biāo)準(zhǔn)差和變異系數(shù)值。

計(jì)算公式分別為[8]

(2)

(3)

1.4 數(shù)據(jù)分析

利用Spss17.0對土壤水分的極大值、極小值、標(biāo)準(zhǔn)差、方差、峰度、偏度S-K檢驗(yàn)進(jìn)行描述性統(tǒng)計(jì)分析，Excel軟件進(jìn)行數(shù)據(jù)計(jì)算并作圖。

2 結(jié)果分析與討論

2.1 監(jiān)測點(diǎn)土壤含水率隨時(shí)間的變化過程分析

土壤水分傳感器測得的5個(gè)測點(diǎn)的30cm、60cm處土壤含水率變化和0～60cm的日平均土壤含水率隨時(shí)間的變化見圖3。

圖3 監(jiān)測點(diǎn)含水率隨時(shí)間的變化

從圖3可以看出，各土層含水率與降水量的變化趨勢基本一致。2020年8月總降水量高達(dá)288.90mm，該觀測期間的土壤含水率相對較高。降水和灌溉是土壤含水率隨時(shí)間變化的兩個(gè)最主要的控制因素。受降水或灌溉的影響，5個(gè)測點(diǎn)土壤含水率的變化比較大。測點(diǎn)Ⅲ的上層土壤含水率大于下層土壤含水率，土壤含水率變化差異較大；其他測點(diǎn)則是下層土壤的含水率較大。遇強(qiáng)降水或者灌水后，土壤含水率整體有顯著增大趨勢，且表層含水率大于下層含水率；長時(shí)間無降水或灌溉，土壤含水率整體有持續(xù)降低的趨勢，且表層土壤水分隨著蒸發(fā)以及土壤水分再分布逐漸低于下層土壤含水率。各測點(diǎn)0～60cm土壤含水率變化范圍分別為16.62%～30.24%、19.48%～30.77%、17.67%～35.42%、20.95%～34.06%、22.02%～32.84%，測點(diǎn)Ⅰ、Ⅱ、Ⅲ的土壤含水率變化幅度較大，最小含水率與最大含水率相差在11.29%～17.75%左右。測點(diǎn)Ⅳ、Ⅴ的土壤含水率變化較為平緩，土壤最低含水率均在20.00%以上，大小差異在10.82%～13.11%之間。

2.2 測點(diǎn)土壤含水率描述性統(tǒng)計(jì)及秩相關(guān)分析

2.2.1 土壤含水率描述性統(tǒng)計(jì)分析

通過Spss17.0統(tǒng)計(jì)分析，給出霍泉灌區(qū)5個(gè)監(jiān)測點(diǎn)不同土層土壤含水率的極大值、極小值、平均值、方差、標(biāo)準(zhǔn)差、偏度系數(shù)、峰度系數(shù)、K-S檢驗(yàn)以及變異系數(shù)的變化范圍，見表1。

表1 土壤含水率特征值

從表1中可以看到，5個(gè)測點(diǎn)不同土層土壤含水率極大值的變化范圍分別為22.54%～35.32%、26.60%～35.30%、22.60%～35.42%，土壤含水率極小值的變化范圍分別為16.34%～28.70%、14.20%～29.05%、16.62%～29.26%。偏度和峰度系數(shù)表明，5個(gè)測點(diǎn)土壤含水率曲線有的呈左偏，有的呈右偏，有的陡峭、有的平緩。對5個(gè)測點(diǎn)的數(shù)據(jù)進(jìn)行K-S檢驗(yàn)，其P值均大于0.05，表明不同土層土壤含水率符合正態(tài)分布。變異系數(shù)表明，測點(diǎn)30cm處土壤含水率變異系數(shù)的變化范圍為0.01～0.26；60cm處土壤含水率變異系數(shù)的變化范圍為0.02～0.29；0～60cm平均土壤含水率變異系數(shù)的變化范圍為0.01～0.27。各層土壤含水率變異系數(shù)的最小值均小于0.1、最大值均小于1，均呈弱變異和中等變異。

2.2.2 測點(diǎn)間土壤含水率的秩相關(guān)分析

利用Spearman秩相關(guān)系數(shù)對同測點(diǎn)土壤含水率相關(guān)性進(jìn)行分析。表2分別表示30cm、60cm、0～60cm土壤含水率秩相關(guān)系數(shù)。5個(gè)測點(diǎn)的土壤含水率的秩相關(guān)系數(shù)在0.01水平上總體上呈顯著相關(guān)。0～60cm土壤含水率秩相關(guān)系數(shù)變化范圍為0.21～0.83，平均值為0.39。30cm土層含水率秩相關(guān)系數(shù)，變化于0.03～0.76，均值為0.40；60cm土層含水率秩相關(guān)系數(shù)，變化于0.02～0.62，均值為0.34。

表2 各測點(diǎn)間土壤含水率的秩相關(guān)系數(shù)

續(xù)表

2.3 土壤含水率特征值隨時(shí)間的變化過程分析

5個(gè)測點(diǎn)土壤含水率的最大值、最小值和平均值以及方差、標(biāo)準(zhǔn)偏差、變異系數(shù)隨時(shí)間的變化過程見圖4。

從圖4可以看出各層土壤含水率及其變異系數(shù)均呈波動(dòng)性變化。土壤含水率越高，其變異系數(shù)越小。降雨強(qiáng)度的大小在一定程度上影響土壤含水率的變化。30cm、60cm、0～60cm處的逐日土壤含水率之間差異較大，土壤含水率的變化幅度分別為18.98%、21.10%、18.80%，60cm處土壤含水率的變幅最大。受灌溉或降水以及作物生長的影響，土壤含水率變幅較大，土壤含水率的變異系數(shù)存在較大差異，30cm與0～60cm平均土壤含水率的變異系數(shù)變化趨勢基本一致，土壤含水率變異系數(shù)曲線波動(dòng)較大，變異系數(shù)值有大有小。60cm處土壤含水率的變異系數(shù)變化相對比較穩(wěn)定，只有20個(gè)觀測日期的變異系數(shù)值是小于0.1的，而其余觀測日期的變異系數(shù)值均在0.10～0.29之間，說明60cm深度土層的土壤水分變異程度基本為中等變異，其變異系數(shù)普遍高于30cm和0～60cm土壤含水率的變異系數(shù)，各測點(diǎn)之間土壤含水率變幅較大，土壤含水率最低為14.20%，最高為35.30%。

3 結(jié) 論

山西省霍泉灌區(qū)試驗(yàn)區(qū)土壤含水率及其變異系數(shù)隨時(shí)間呈波動(dòng)性變化。60cm土層的土壤含水率總體大于30cm土層的土壤含水率。兩個(gè)土層間的土壤含水率具有一定差異。土壤水分變異系數(shù)隨土層深度增加而增加。表層土壤的含水率及變異系數(shù)波動(dòng)較大，下層土壤的含水率及變異系數(shù)之間差異性較小，變化相對比較穩(wěn)定。灌溉時(shí)期的土壤含水率及其變異系數(shù)變幅較大。降水和灌溉對土壤含水率的影響較大。灌溉或降水后的土壤含水率較大，其變異系數(shù)較小。30cm和60cm的土壤含水率的變異系數(shù)均在0.01～1.00之間，屬弱變異和中等變異。土壤含水率的秩相關(guān)性分析表明5個(gè)測點(diǎn)相互間具有顯著相關(guān)性。綜上所述，研究土壤含水率隨時(shí)間的變化和空間變異性以及各監(jiān)測點(diǎn)間不同土層含水率的相關(guān)性，可對山西省霍泉灌區(qū)利用云平臺和物聯(lián)網(wǎng)技術(shù)實(shí)施智能化、數(shù)字化的精準(zhǔn)灌溉提供理論基礎(chǔ)，為鄉(xiāng)村振興數(shù)字化轉(zhuǎn)型貢獻(xiàn)力量。

圖4 土壤含水率特征值隨時(shí)間的變化