張瑜瑜,栗現(xiàn)文
(1.新疆水利水電勘測(cè)設(shè)計(jì)研究院,新疆 烏魯木齊 830000;2.西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院,陜西 楊凌 712100)
水資源短缺和土壤鹽漬化是干旱區(qū)農(nóng)業(yè)可持續(xù)發(fā)展面臨的兩個(gè)主要問題。我國(guó)西北干旱區(qū)氣候干燥,蒸發(fā)強(qiáng)烈,降水稀少,以新疆南部庫爾勒地區(qū)為例,其多年平均降水量?jī)H58.6 mm[1],蒸降比極大,對(duì)鹽漬化風(fēng)險(xiǎn)尤為敏感。膜下滴灌可防止灌水深層滲漏,而這也使鹽分具備了在土壤淺部積聚的條件。栗現(xiàn)文等[2,3]的研究表明,微咸水膜下滴灌生育期結(jié)束后會(huì)存在大量影響下一期棉花出苗的鹽斑,為可持續(xù)利用微咸水,采取有效措施調(diào)控土壤水鹽狀況成為微咸水安全利用的基礎(chǔ)[4]。田間氣象特征、土壤水分運(yùn)動(dòng)及重分布、地下水位水質(zhì)等對(duì)土壤剖面鹽分分布具有重要影響,這些控鹽因子漫灌前后的動(dòng)態(tài)演化過程直接決定著田間土壤鹽分的空間分布特征。春季漫灌會(huì)造成地下水位大幅抬升,當(dāng)?shù)叵碌V化度較高且水位較淺時(shí),毛細(xì)水上升作用,是土壤次生鹽漬化的重要誘因[5]。因此,對(duì)漫灌前后土壤控鹽因子動(dòng)態(tài)演化特性進(jìn)行系統(tǒng)分析,是田間土壤鹽分調(diào)控的重要基礎(chǔ),具有重要理論與實(shí)際意義。
研究區(qū)設(shè)有新疆塔里木河流域巴音郭勒管理局水利科研所,為國(guó)家水利部重點(diǎn)灌溉試驗(yàn)站,耕地面積8000 余畝,位于典型干旱區(qū)新疆巴音郭楞蒙古自治州首府庫爾勒市境內(nèi),在該市經(jīng)濟(jì)技術(shù)開發(fā)區(qū)西尼爾鎮(zhèn)以南3 km 處,距市區(qū)23 km。地理坐標(biāo)東經(jīng)86°10′,北緯41°35′,為天山南麓塔里木盆地邊緣孔雀河沖積平原帶,地勢(shì)平緩,海拔900 m 左右;屬暖溫帶大陸性荒漠氣候,干旱少雨,蒸發(fā)強(qiáng)烈,晝夜溫差大,多年平均降水量58.6 mm,年最大蒸發(fā)量2788.2 mm(E20 蒸發(fā)皿),年平均風(fēng)速2.4 m/s,最大風(fēng)速22 m/s,年平均氣溫11.5℃,最低氣溫-30.9℃,最高氣溫42.2℃,無霜期191 d。灌溉用微咸地下水溶解性總固體(TDS)2.4 g/L 左右,漫灌洗鹽用河渠淡水TDS 為0.6 g/L 左右。試驗(yàn)田以壤土(美國(guó)農(nóng)業(yè)部土壤質(zhì)地分類)為主[15],土質(zhì)較均一,砂粒、粉粒和粘粒含量分別為46.81%、45.96%和7.23%,干密度為1.63 g/cm3。地下水位埋深在5.8 m左右。
1.2.1 漫灌洗鹽設(shè)置
試驗(yàn)區(qū)自2008 年起,部分田塊已連續(xù)8 年采用微咸水進(jìn)行膜下滴灌,為控制田間土壤鹽分均采用在非生育期漫灌方式淋洗鹽分。本研究在2016 年4 月29 日至4 月30 日持續(xù)開展淡水漫灌洗鹽實(shí)驗(yàn)。為增加可控性及避免對(duì)大田實(shí)際生產(chǎn)造成影響,試驗(yàn)采用田間打埂方式圈定10 m×10 m 小區(qū),埂高25 cm左右。參照大田漫灌定額,確定試驗(yàn)小區(qū)漫灌量為200 m3/畝,即300 mm。
1.2.2 土壤樣品采集與測(cè)定
用土鉆采取不同深度分層土樣,每10 cm 采集一組。其中,取樣深度為100 cm(設(shè)定作物對(duì)水鹽最大潛在利用深度位100 cm)和60 cm(生育期滴灌主要影響范圍)的隨機(jī)采樣點(diǎn)各18 個(gè),單次采集288 組,6 個(gè)監(jiān)測(cè)日共采集1728 組土樣(漫灌前:4 月26 日;漫灌后:5 月6 日、5 月10 日、5 月15 日、5 月20 日及6月14 日。將采集土壤樣品烘干后測(cè)定含水率。
1.2.3 氣象監(jiān)測(cè)
為明確灌水洗鹽至播種前氣象條件,利用Davis Vantage Pro2 自動(dòng)氣象站采集氣象數(shù)據(jù),數(shù)據(jù)采集間隔為1 小時(shí),主要監(jiān)測(cè)項(xiàng)目有:降水、空氣溫濕度、風(fēng)向、風(fēng)速、太陽輻射等。
1.2.4 地下水動(dòng)態(tài)監(jiān)測(cè)
布設(shè)地下水位監(jiān)測(cè)井,采用Levelogger 采集地下水位信息。
選取2016 年4 月至2016 年6 月采集的相關(guān)氣象監(jiān)測(cè)數(shù)據(jù),開展播前漫灌前后氣象因子動(dòng)態(tài)變化特征分析。潛在蒸發(fā)量ET0(mm/h)根據(jù)國(guó)際糧農(nóng)組織FAO 推薦的Penman-Monteith公式計(jì)算[6]:
式中,ET0為參考作物潛在蒸散發(fā)量,mm/hour;Rn為植被層凈輻射,MJ/m2;G 為土壤地?zé)嵬浚∕J/m2)/hour;Thr為小時(shí)平均氣溫,℃;△為氣溫在Thr時(shí)的飽和水汽壓梯度,kPa/℃;γ 為干濕球常數(shù),kPa/℃;e0(Thr) 為氣溫為Thr時(shí)的飽和水汽壓,kPa;ea為小時(shí)平均實(shí)際水汽壓,kPa;u2為小時(shí)平均風(fēng)速,m/s。各變量的詳細(xì)確定方法見文獻(xiàn)[16]。基于監(jiān)測(cè)田間氣象站監(jiān)測(cè)頻率計(jì)算ET0,日內(nèi)求和得出當(dāng)日潛在蒸發(fā)量。
ET0與氣溫和太陽輻射的動(dòng)態(tài)變化趨勢(shì)基本一致,5 月份開始逐漸增大,4、5、6 月日內(nèi)Ep均值分別為3.0 mm、3.8 mm 及5.8 mm,4 月~6 月累計(jì)為386.2 mm。監(jiān)測(cè)期降水量較小,4 月~6 月累計(jì)降水量為31.7 mm,次最大降水量13.1 mm。潛在蒸發(fā)強(qiáng)度遠(yuǎn)大于降水影響,蒸降比為12.2。因此,在氣象因子作用下,灌水結(jié)束后田間土壤蒸發(fā)強(qiáng)烈,是促進(jìn)鹽分隨蒸發(fā)向表層遷移的重要條件。
圖1 2016 年漫灌前后氣象因子動(dòng)態(tài)變化
對(duì)監(jiān)測(cè)小區(qū)0 cm~100 cm 取樣點(diǎn)各層位體積含水量進(jìn)行統(tǒng)計(jì)分析,不同監(jiān)測(cè)時(shí)間剖面土壤水分動(dòng)態(tài)變化如圖2 所示,其中40 cm~60 cm 土壤含水量數(shù)據(jù)離散性相對(duì)較大。在試驗(yàn)區(qū)采用圍框法測(cè)定的田間持水率為27.1%[7]。為度量剖面土壤水分的數(shù)量變化,采用土壤水儲(chǔ)存量作為指標(biāo)進(jìn)行計(jì)算,其在數(shù)值上為土壤體積含水量分布函數(shù)θ(z,t)在作物潛在可利用深度上的積分[8]:
式中,W(t)為t 時(shí)刻的土壤水儲(chǔ)存量,mm;z 為垂向坐標(biāo),地面為0 點(diǎn);d 為作物潛在可利用深度(此處為絕對(duì)值),mm。本研究實(shí)測(cè)各深度土壤體積含水量為對(duì)應(yīng)土層的平均含水量,其分布如圖2 所示。采用分段計(jì)算求和后,可得不同監(jiān)測(cè)時(shí)間土壤剖面儲(chǔ)存量變化(表2)。
為比較不同深度范圍各層位在不同時(shí)刻土壤水儲(chǔ)存量變化,根據(jù)式(3),本研究定義δ 為土壤水儲(chǔ)存變化強(qiáng)度這一指標(biāo)(物理意義等同于這兩個(gè)時(shí)刻土壤含水率變化的均值,無量綱):
式中符號(hào)意義同前,t1,t2為選定比較的兩個(gè)時(shí)刻,t2>t1。
漫灌前土壤含水量較低,漫灌后第6 日含水率在各層位均有大幅增加。但在土面蒸發(fā)和深層滲漏作用下,隨時(shí)間變化,除100 cm 深度外,剖面含水量呈不斷降低趨勢(shì)。灌后100 cm 深度在灌后20 日基本維持在0.37 左右,僅在灌后20 日~45 日略有降低,而田間實(shí)測(cè)飽和水量約為0.4,表明其受土面蒸發(fā)作用影響較弱,應(yīng)接近土壤剖面不受蒸散發(fā)影響的最深零通量面位置。
從隨時(shí)間變化水分分布動(dòng)態(tài)及剖面土壤水儲(chǔ)存變化量看,灌后6 日~10 日,100 cm 以上各土層含水量均有大幅降低。值得注意的是,在漫灌后45 日(6 月15 日)水分剖面分布與漫灌前(4 月23 日),在40 cm~90 cm 有較好的一致性,即在持續(xù)蒸發(fā)作用下達(dá)到剖面水分的近似穩(wěn)定狀態(tài);而0 cm~30 cm 土壤水分小于灌前,與氣象條件有一定關(guān)系,漫灌前的冬春季節(jié)有部分降水,而氣溫與太陽輻射均相對(duì)較低,潛在蒸發(fā)作用較弱,凍土狀態(tài)有利于水分保持;而漫灌后蒸散發(fā)作用逐漸增強(qiáng),表層土壤水分的散失也就更為強(qiáng)烈。
(數(shù)據(jù)點(diǎn)為同一層位均值,誤差棒為標(biāo)準(zhǔn)差)
表2 2016 年漫灌前后不同監(jiān)測(cè)時(shí)間剖面土壤水儲(chǔ)存量(mm)
利用三參數(shù)水位計(jì)采集的地下水位、水質(zhì)及溫度數(shù)據(jù),可獲取春季漫灌時(shí)段各地下水相關(guān)因子的動(dòng)態(tài)變化(圖3)。綜合分析來看,地下水對(duì)漫灌是有近似即時(shí)響應(yīng)的,而地下水位介于5.83 m~6.2 m,在監(jiān)測(cè)時(shí)段內(nèi)并未表現(xiàn)出明顯的滯后響應(yīng),這應(yīng)歸因于田間非飽和帶是存在優(yōu)先流對(duì)地下水補(bǔ)給的,但田間地下水動(dòng)態(tài)還是受控于區(qū)域地下水流場(chǎng)變化。
圖3 漫灌時(shí)段地下水位、水質(zhì)及溫度變化
本文對(duì)漫灌前后監(jiān)測(cè)時(shí)段(計(jì)45 日)的控鹽因子:氣象要素、剖面土壤含水率與地下水位、水質(zhì)動(dòng)態(tài)進(jìn)行系統(tǒng)分析,主要得出以下結(jié)論:
(1)氣象要素動(dòng)態(tài)顯示,試驗(yàn)時(shí)段(2016 年4 月~6 月)氣溫與太陽輻射均波動(dòng)性較強(qiáng),太陽輻射從5 月中下旬開始有大幅增加,氣溫和太陽輻射的大幅增加對(duì)出苗及苗期生長(zhǎng)較為有利。潛在蒸發(fā)量ET0與氣溫及太陽輻射的動(dòng)態(tài)變化趨勢(shì)基本一致,遠(yuǎn)大于降水影響,蒸降比達(dá)12.2,氣象因素表層鹽分運(yùn)移的重要驅(qū)動(dòng)因子。
(2)40 cm 處土壤容重較大,持水性相對(duì)較好。
(3)地下水位處于蒸發(fā)最大影響深度以下,在播前非生育期不存在潛水蒸發(fā)返鹽對(duì)表層土壤鹽分的作用。