大數(shù)據(jù)在人工草地管理上的應(yīng)用初探

景 文，王保林，樊俊梅

(內(nèi)蒙古蒙草生態(tài)環(huán)境(集團(tuán))股份有限公司，內(nèi)蒙古呼和浩特 010030)

1 研究背景及意義

優(yōu)質(zhì)高產(chǎn)的人工草地，是草地農(nóng)業(yè)現(xiàn)代化的必要條件和重要標(biāo)志。不僅可以有效提供足量的優(yōu)質(zhì)牧草，有利于畜牧業(yè)生產(chǎn)方式的轉(zhuǎn)變和草原畜牧業(yè)的可持續(xù)發(fā)展，也在環(huán)境保護(hù)與環(huán)境產(chǎn)業(yè)中起到十分重要的作用〔1〕。要使得人工草地大力發(fā)展，可通過(guò)一系列有效的技術(shù)或管理方式實(shí)現(xiàn)整個(gè)人工牧草生產(chǎn)系統(tǒng)內(nèi)的物能良性轉(zhuǎn)化，以獲得較高的生產(chǎn)效益〔2〕。如實(shí)施灌溉，能夠?qū)Π敫珊档貐^(qū)人工優(yōu)質(zhì)飼草產(chǎn)量的提高和品質(zhì)改善起到重要的調(diào)控作用，但需加強(qiáng)灌溉管理(如灌溉方式、灌水量、灌水時(shí)間/頻次等)〔3〕，通過(guò)摸清植物的耗水規(guī)律，掌握其各個(gè)物候期土壤水分控制指標(biāo)，形成準(zhǔn)確的灌溉制度，為科學(xué)灌溉及種植管理提供技術(shù)支撐，這對(duì)人工牧草產(chǎn)業(yè)發(fā)展是極為重要的〔4〕。當(dāng)然，除了人工牧草領(lǐng)域相關(guān)技術(shù)不斷創(chuàng)新優(yōu)化外，還要依靠信息技術(shù)〔5〕。近年來(lái)，隨著現(xiàn)代農(nóng)牧業(yè)自動(dòng)化需求的不斷增長(zhǎng)和信息技術(shù)的發(fā)展，農(nóng)牧業(yè)物聯(lián)網(wǎng)在物理感知、數(shù)據(jù)傳輸、智能處理、應(yīng)用服務(wù)等領(lǐng)域取得了重要進(jìn)展〔6〕。我國(guó)作為一個(gè)農(nóng)牧業(yè)大國(guó)，應(yīng)結(jié)合實(shí)際情況，研究開發(fā)出順國(guó)情、合地情的農(nóng)牧業(yè)物聯(lián)網(wǎng)技術(shù)和設(shè)備，努力擴(kuò)展農(nóng)牧業(yè)物聯(lián)網(wǎng)技術(shù)在我國(guó)的試點(diǎn)應(yīng)用。

本研究基于物聯(lián)網(wǎng)技術(shù)，對(duì)人工草地的微氣象和土壤的溫濕度進(jìn)行監(jiān)測(cè)，分析土壤中水分的含量變化規(guī)律，以此作為判定人工草地灌溉正確時(shí)機(jī)的依據(jù)，從水分診斷與灌溉決策方面，實(shí)現(xiàn)人工草地生產(chǎn)的智能管理。

2 試驗(yàn)區(qū)概況

試驗(yàn)區(qū)位于內(nèi)蒙古自治區(qū)赤峰市阿魯科爾沁旗(以下簡(jiǎn)稱“阿旗”)，屬中溫帶半干旱大陸性季風(fēng)氣候區(qū)，四季比較分明。年均氣溫5.5℃，極端最高氣溫40.6℃，極端最低氣溫-32.7℃，全年10℃以上積溫1900-3100℃，從南向北遞減。年日照時(shí)數(shù)2767.3-3034.4小時(shí)。無(wú)霜期95-140天。年均降水量300-400mm，多集中在7-8月份。年均蒸發(fā)量1958.1-2259.0mm。年風(fēng)速3.1-2.7m/s，冬春季節(jié)多大風(fēng)，全年8級(jí)以上大風(fēng)日有20-35天。地帶性植被由北向南主要有山地草甸草原植被、沙地植被二種類型。土壤主要以黑鈣土、栗鈣土和風(fēng)沙土等為主。水資源豐富，總量3.73億m3，其中可利用水資源2.41億m3〔7〕。

3 試驗(yàn)研究方法

研究對(duì)象屬性：赤峰市阿魯科爾沁旗蒙草基地人工牧草生產(chǎn)單元。

研究設(shè)備名稱及監(jiān)測(cè)采集數(shù)據(jù)時(shí)間范圍：“天圻”(2017.5.4-9.27)、“智墑”(2017.5.4-2018.12.31)。

研究設(shè)備監(jiān)測(cè)采集數(shù)據(jù)指標(biāo)(節(jié)點(diǎn)數(shù)量，步長(zhǎng))：降雨量(1，1h)、風(fēng)速(1，1h)、風(fēng)向(1，1h)、土壤水分含量(10，1h)。

研究監(jiān)測(cè)采集數(shù)據(jù)統(tǒng)計(jì)軟件：Excel。

4 結(jié)果與分析

4.1 降雨量分析

由表1顯示，2017年6-9月間阿旗降雨主要集中在8月，8月降雨量總和占總降雨量的58.33%。

表1 2017年6-9月阿旗降雨量月份統(tǒng)計(jì)表

2017年6月1日-9月27日，119天內(nèi)，降雨出現(xiàn)10天，暴雨出現(xiàn)1天，在8月9日(54.676mm/d)，中雨出現(xiàn)1天，小雨出現(xiàn)8天。中、小雨出現(xiàn)天數(shù)占降雨出現(xiàn)天數(shù)90%。

表2 2017年6月1日-9月27日阿旗日降雨量統(tǒng)計(jì)表

4.2 風(fēng)速分析

2017年5月4日-9月27日，147天內(nèi)，六級(jí)風(fēng)(10.8-13.9m/s)出現(xiàn)2天，分別在5月5、6日(西北風(fēng))，五級(jí)風(fēng)(8-10.8m/s)出現(xiàn)8天，分別在5月4日(西南風(fēng))、5月5-7、11、28日(西北風(fēng))、7月8日(西南風(fēng))、9月5日(西北風(fēng))，四級(jí)風(fēng)(5.5-8m/s)出現(xiàn)36天，三級(jí)風(fēng)(3.4-5.5m/s)出現(xiàn)72天，三四級(jí)風(fēng)出現(xiàn)天數(shù)占比為73.47%。

4.3 風(fēng)向分析

由表3顯示，2017年5-9月間阿旗風(fēng)向數(shù)據(jù)檢測(cè)采集15089次，主導(dǎo)風(fēng)向?yàn)槲髂巷L(fēng)和西北風(fēng)，兩者發(fā)生的頻率分別為34.64%和24.06%。

表3 2017年5月4日-9月27日阿旗風(fēng)向統(tǒng)計(jì)表

4.4 土壤凍融過(guò)程分析

通過(guò)2018年的監(jiān)測(cè)數(shù)據(jù)來(lái)看，試驗(yàn)區(qū)土壤從3月13日開始進(jìn)入融化過(guò)程，從12月1日開始進(jìn)入冰凍過(guò)程，且上層土壤最先進(jìn)入融凍過(guò)程。

圖1 土壤融化過(guò)程(上)和土壤冰凍過(guò)程(下)圖

4.5 土壤十層水分分析

從圖2可看出，1-4月上旬土壤處于冰凍階段，從4月中旬至10月為種植、收割階段，從11月中旬開始土壤逐漸又進(jìn)入冰凍階段。

圖2 2017年5月20日-11月7日(上)和2018全年(下)十層水分波動(dòng)圖

4.6 土壤分層水分分析

4.6.1 土壤水分分析——10cm

10cm土層水分變化直接直觀反映出降雨或灌溉頻次和程度，每一波峰反映一次降雨或者灌溉。2017年10cm土層歷史最高含水量26.71%，出現(xiàn)在6月7日；歷史最低含水量5.67%，出現(xiàn)在7月16日。2018年10cm歷史最高含水量26.58%，發(fā)生在8月12日；歷史最低含水量4.36%，發(fā)生在8月6日。

圖3 2017年(上)和2018年(下)10cm土層水分波動(dòng)圖

4.6.2 土壤水分分析——20cm、30cm、40cm、50cm、60cm

由植物ET根系分布/分層耗水比例圖分析，在植物關(guān)鍵需水生育期內(nèi)，ET根系分布在地表至60cm左右的土層間，相對(duì)應(yīng)的耗水量總和為100%。植物ET根系分布/分層耗水量在地表至60cm土層隨土層深度增大而逐層減小。推測(cè)試驗(yàn)區(qū)最深植物ET根系為60cm。

圖4 2017年(左)和2018年(右)ET根系分布/分層耗水量比例圖

以下著重對(duì)20-60cm水分進(jìn)行分析。如圖5所示，2017年20cm土層歷史最高含水量21.84%，出現(xiàn)在9月14日；歷史最低含水量5.46%，出現(xiàn)在7月16日。2018年20cm歷史最高含水量21.27%，發(fā)生在4月25日；歷史最低含水量3.98%，發(fā)生在8月7日。

圖5 2017年(上)和2018年(下)20cm土層水分波動(dòng)圖

圖6所示，2017年30cm土層歷史最高含水量22.76%，出現(xiàn)在8月10日，歷史最低含水量6.11%，出現(xiàn)在7月17日。2018年30cm歷史最高含水量22.47%，發(fā)生在4月25日；歷史最低含水量4.85%，發(fā)生在8月6日。

圖6 2017年(上)和2018年(下)30cm土層水分波動(dòng)圖

圖7所示，2017年40cm土層歷史最高含水量20.16%，出現(xiàn)在5月22日，歷史最低含水量5.55%，出現(xiàn)在7月20日。2018年40cm歷史最高含水量19.59%，發(fā)生在4月25日；歷史最低含水量5.06%，發(fā)生在8月7日。

圖7 2017年(上)和2018年(下)40cm土層水分波動(dòng)圖

圖8所示，2017年50cm土層歷史最高含水量21.97%，出現(xiàn)在5月22日，歷史最低含水量6.82%，出現(xiàn)在7月21日。2018年50cm歷史最高含水量21.23%，發(fā)生在4月25日；歷史最低含水量6.82%，發(fā)生在8月6日。

圖8 2017年(上)和2018年(下)50cm土層水分波動(dòng)圖

圖9所示，2017年60cm土層歷史最高含水量21.06%，出現(xiàn)在5月22日，歷史最低含水量7.24%，出現(xiàn)在7月21日。2018年60cm歷史最高含水量21.19%，發(fā)生在4月25日；歷史最低含水量9.46%，發(fā)生在8月5日。

圖9 2017年(上)和2018年(下)60cm土層水分波動(dòng)圖

由以上分層水分波動(dòng)曲線可知，每一層的歷史最低含水量均發(fā)生在植物刈割前期的控水階段，且相對(duì)上層的根系先進(jìn)入水分脅迫狀態(tài)，每層根系進(jìn)入脅迫狀態(tài)的時(shí)間均集中在2017年7月16日-7月21日和2018年8月5日-8月7日中間，隨著時(shí)間推移逐層根系進(jìn)入水分脅迫狀態(tài)。

4.6.3 不同月份土層水分分析

以13#噴灌圈(種植燕麥)為例，根據(jù)2017年“智墑”監(jiān)測(cè)數(shù)據(jù)，逐月判斷灌溉是否合理。

(1)6月十層水分之和分析與分層分析

由圖10可知，最后一次灌溉為6月20日，在此之前土壤總水分增加約16mm；由于噴灌機(jī)控制的面積較小，每次灌溉基本都能到達(dá)50cm深度，較為合理；6月25日之后主要根系的活躍程度才逐漸提升，反映出之前灌溉過(guò)于頻繁，土壤水分偏高。

圖10 2017年6月土壤水分波動(dòng)圖

(2)7月十層水分之和分析與分層分析

圖11顯示出7月追加了灌溉，十層水分之和達(dá)到181mm；7月7日前無(wú)灌溉，結(jié)合6月數(shù)據(jù)可知，活躍“ET”根系已經(jīng)均在脅迫狀態(tài)，因此判斷灌溉不合理；7月7日-16日無(wú)灌溉，判斷為收割期；7月16日之后灌溉開始頻繁，出現(xiàn)了一天之中多次水分上升的情形，需結(jié)合現(xiàn)場(chǎng)分析原因。

圖11 2017年7月土壤水分波動(dòng)圖

(3)8月十層水分之和分析與分層分析

圖12顯示，8月十層水分總計(jì)減少了11mm，8月9日有效降雨為35mm，占總降雨量的64%；10cm水分曲線多次上升，且出現(xiàn)每天多次灌溉的情況，需結(jié)合現(xiàn)場(chǎng)分析原因；在10cm根系最活躍的情況下，對(duì)應(yīng)的水分值區(qū)間為11-16，而土壤大部分時(shí)間的水分值高于該區(qū)間上限；100cm曲線有三次上升，其中一次是降雨影響，另外兩次都是灌溉過(guò)于頻繁導(dǎo)致的；根據(jù)“智墑”分析，“ET”根系已經(jīng)發(fā)展至90cm，但主要部分還是集中在40cm以上。

圖12 2017年8月土壤水分波動(dòng)圖及YH圖及ET根系分布

(4)9月十層水分之和分析與分層分析

圖13顯示，9月十層水分總計(jì)增加了11mm；依然出現(xiàn)每天多次灌溉的情況，比8月還要頻繁，需結(jié)合現(xiàn)場(chǎng)分析原因；8月11-19日以及9月2日前后，作物處于較為舒服的狀態(tài)，9月2日之后，灌溉頻繁、土壤水分偏高，“ET”根系反而不活躍。

圖13 2017年9月土壤水分波動(dòng)圖

(5)10月十層水分之和分析

圖14顯示，整個(gè)10月沒(méi)有灌溉，十層水分總計(jì)減少了35mm。

圖14 2017年10月土壤水分波動(dòng)圖

4.7 活動(dòng)溫度分析

“智墑”除了監(jiān)測(cè)地下十層的土壤水分外，還可監(jiān)測(cè)地表溫度。某一段時(shí)間內(nèi)逐日平均溫度≥10℃持續(xù)期間日平均氣溫的總和，即活動(dòng)溫度總和，簡(jiǎn)稱積溫。在其他環(huán)境條件基本滿足的前提下，在一定的溫度范圍內(nèi)，溫度與生物有機(jī)體發(fā)育速度之間呈正相關(guān)。表4顯示，平均、最高和最低地表溫度隨月份的推移順序均呈現(xiàn)出先增加后逐漸降低，2017年9月26日地表均溫為10.19℃，2018年9月29日地表均溫為10.03℃，此后接近或低于10℃。降雨量和地表溫度在時(shí)間維度分布呈現(xiàn)不一致，水熱條件不是同期出現(xiàn)。

表4 地表空氣溫度月均溫統(tǒng)計(jì)表

13#圈的地表溫度在夜間-凌晨時(shí)間段頻繁出現(xiàn)低于零下的情況；10cm位置晝夜溫差也比較大，均在2℃以上；100cm位置溫度晝夜溫差極小，且溫度下降極為緩慢。

圖15 10月份地表、10cm、100cm溫度曲線

5 初步結(jié)論

根據(jù)監(jiān)測(cè)數(shù)據(jù)分析，初步總結(jié)出能代表該地生產(chǎn)單元的降雨量、風(fēng)向、風(fēng)速、植物ET根系分布/分層耗水量、土層水分、地表積溫等指標(biāo)數(shù)據(jù)的客觀實(shí)際和變化規(guī)律：

(1)降雨主要集中在8月，多出現(xiàn)三四級(jí)風(fēng)，且主導(dǎo)為西南風(fēng)和西北風(fēng)。

(2)1-4月上旬土壤處于冰凍階段，4月中旬到10月為種植收割階段，從11月中旬開始土壤逐漸又進(jìn)入冰凍階段。

(3)植物ET根系分布/分層耗水量在地表至60cm土層隨深度增大而逐層減小，其活躍程度僅能對(duì)灌溉管理起到一定程度的指導(dǎo)作用，還需結(jié)合植物生育期需求進(jìn)行生產(chǎn)管理。

(4)6-9月灌溉量基本合理，能灌溉到牧草活躍“ET”根系位置，但普遍存在灌溉周期過(guò)短、頻繁灌溉的問(wèn)題，導(dǎo)致et數(shù)據(jù)偏高，無(wú)法得出科學(xué)的作物系數(shù)，造成水肥浪費(fèi)。灌溉建議是，增大灌溉周期，即4天灌溉一次，增加灌溉量，保證每次灌溉至40cm深度。刈割后進(jìn)入10月中旬開始澆越冬水，需關(guān)注冬灌與溫度的關(guān)系。

通過(guò)對(duì)現(xiàn)存指標(biāo)數(shù)據(jù)分析，探討式分析了降雨量、土壤墑情、植物ET根系分布/分層耗水量等指標(biāo)的關(guān)聯(lián)關(guān)系。進(jìn)一步對(duì)影響植物生長(zhǎng)相關(guān)的氣象、土壤指標(biāo)數(shù)據(jù)模擬分析，增加對(duì)太陽(yáng)輻射強(qiáng)度、地上生物量等指標(biāo)監(jiān)測(cè)采集，綜合數(shù)學(xué)、統(tǒng)計(jì)學(xué)等對(duì)數(shù)據(jù)進(jìn)行跨學(xué)科分析，建立模擬地上生物量預(yù)測(cè)模型，指導(dǎo)以管控灌溉、施肥為核心的人工牧草生產(chǎn)管理，預(yù)測(cè)以水熱條件影響為主的天然草產(chǎn)量。為實(shí)現(xiàn) “智能管控”和“智慧草業(yè)”建立一定的理論、實(shí)踐基礎(chǔ)。