區域鹽漬化土壤表層水鹽變化特征及空間變異分析

趙文娟， 黃凌

(1.寧夏大學土木與水利工程學院，寧夏銀川 750021； 2.寧夏節水灌溉與水資源調控工程技術研究中心，寧夏銀川 750021)

?

區域鹽漬化土壤表層水鹽變化特征及空間變異分析

趙文娟1，2， 黃凌1，2

(1.寧夏大學土木與水利工程學院，寧夏銀川 750021； 2.寧夏節水灌溉與水資源調控工程技術研究中心，寧夏銀川 750021)

摘要表層土壤水鹽的變化特征和空間變異性是研究改良土壤鹽漬化的基礎。將試驗與統計學結合，分析寧夏銀北試驗區0~25 cm土層的水鹽變化、空間變異性。結果表明，春夏季表層土壤水分含量、鹽分含量總體隨時間逐漸增大；水分變異系數處于中等變異，鹽分變異系數處于弱變異區間內。通過統計計算，認為試驗測定數目是達到精度要求的。利用Matlab軟件擬合半方差函數求參，顯示出土壤水分的空間自相關程度要明顯高于土壤鹽分。

關鍵詞區域化變量；空間變異；鹽漬土；土壤表層

Characteristics and Spatial Variability of Surface Soil Salinity and Moisture in Saline Area

ZHAO Wen-juan1,2, HUANG Ling1,2(1.School of Civil and Hydraulic Engineering, Ningxia University, Yinchuan, Ningxia 750021; 2.Ningxia Research Center of Technology on Water Saving Irrigation and Water Resources Regulation, Yinchuan, Ningxia 750021)

Abstract The characteristics and spatial distribution of surface saline soil is important to study in Yinbei area of Ningxia Province.This paper uses the methods of geographic statistics to analyze the variation of soil water(0-25 cm) and salt, spatial variability in experimental plots.The results show that surface soil moisture content and salt content was gradually increased over time in summer.Second, the coefficient of variation of soil moisture is higher than the soil salt’s.According to statistical calculation, the number of test measurement is to achieve the required precision.Last, the spatial autocorrelation of soil moisture level is significantly higher than that of soil salinity through adopting Matlab software parameters by fitting semivariogram.

Key words Regionalized variables; Spatial variability; Saline soil; Topsoil

寧夏銀北地區位于寧夏回族自治區銀川平原的北部和鄂爾多斯臺地西緣的高階地上，屬于沖湖積平原潛水-承壓水區。由于地形低洼、地下水位高以及人為的農業生產方式，該地區形成大面積的鹽漬化土地。鹽分通常積聚在25 cm以上的土層中。深翻土層、大水溝灌及優化施肥等傳統降鹽措施無法完全解決現階段的土壤鹽漬化問題，所以掌握土壤水鹽的分布特征和時空變異可為制定合理的土壤鹽漬化調控模式提供參考[1]。筆者就研究區耕作層土壤在春夏季的水分和鹽分含量變化進行比較，分析不同時期表層土壤水鹽的演變過程，為后續的土壤水鹽運移數值模擬研究提供可靠的試驗數據。

1材料與方法

1.1試驗區域概述土壤水鹽監測區域位于寧夏省銀川市賀蘭縣南梁臺子鐵西村。賀蘭縣屬于溫帶干旱氣候區，夏季受東南暖濕氣流的影響，雨水較多，蒸發強烈[2]。監測時段為2012年5~10月，平均氣溫為18.77 ℃，月平均降水量為41.2 mm，水面蒸發量為124.00 mm，月平均地下水埋深為1.32 m。使用HJ16-TFC-1BⅢ型土肥測試儀，對試驗區土壤的化學性質進行測定，發現全鹽量平均值為1.42 g/kg，pH為8.56。

試驗區為一條形區，面積約1 000 m2, 以五行十列布置形式將50根TDR(Time Delay Reflectometry，簡稱TDR)探測管埋設土層內。將時間延遲反射儀與TDR管配合使用，分時段監測該試驗區的50個測點0~25處土壤含水率和電阻值。TDR探測管是用特殊塑料制成，需使用支架、土鉆、固定鋼管、防振錘對其進行施工安裝。探管安裝圖見圖1。

1.2監測方法將TDR探頭放置于地表10 cm深的探測管內，可直接獲得土壤含水率和土壤電導率。通過多點測到的土壤含水率變異系數均小于10%，認為該測定值是可靠的。通過TDR設備可以測定出土壤含水率和一個與電導率相關的參數LEVEL值，但不能直接獲得土壤的含鹽量值。所以，需要建立參數LEVEL值、土壤含水率與土壤電導率三者的函數關系，通過顯著性經驗最終確定上述3個值的線性函數關系式[3]。

Y=1.326 4+2.125 4X1-0.006 9X2

(1)

式中，Y為土壤含鹽量；X1為TDR所測土壤體積含水率；X2為TDR所測 LEVEL值。

2結果分析

2.1土壤表層水鹽變化特征表層土壤平均含水率在監測期以7月最低，以9月最高。月平均含水率最大差值為8.6%，隨著時間的遷移呈由高到底而后逐漸增加的變化特征。各月平均全鹽量最高值出現在8月份，最低值則出現在5月份。結合監測期內不同月份的平均蒸發量、降雨量，與土壤全鹽量、含水率進行對比分析。由圖2可知，7月份降雨量達到最高，蒸發量達到最低，此時土壤水分含量也達到最低，并且土壤內部開始積鹽；隨著蒸發量的逐漸增加，降雨減少，土壤內部的鹽分含量在8月達到最大，隨后蒸發量逐漸減少，降雨又有所增加，土壤鹽分含量隨之降低，但橫斷面上全鹽量未出現明顯的回落。

2.2不同時間土壤水鹽統計特征值由表1可知，5個采樣時間的土壤含水率最小值與最大值相差范圍在48%~71%。變化量逐月增加，到8月28日土壤表層含水率最大值為38%，最小值為19.9%，且整個橫斷面的平均土壤含水率達到最高值。土壤樣本的變異系數是樣本的標準值與平均數的百分比。將變異系數分為不同的區間范圍，用于反映樣本的變異程度，即0~15%為弱變異，16%~35%為中等變異，大于36%為強變異[4]。從變異系數來看，5個時期的土壤水分變異系數中只有7月26日監測的數值為弱變異強度，其他時期土壤含水率的變異系數在0.2~0.4之間，屬于中等變異強度；秋末土壤水分的變異系數大于36%，變異性達到最強，說明數據受土壤含水率影響較大。從取樣數目看，當計算精度為5%時，土壤表層含水率的合理取樣數目在9~18個。

表1　20 cm深土層不同采樣時間土壤水分的統計特征值

由表2可知，土壤全鹽量變幅比土壤含水率變幅大，范圍在20%~70%，土壤全鹽量在7月26日的監測值變幅最小，標準差也是5個時期的最小值，且土壤表層全鹽量變異系數為0.084，屬于弱變異強度。變異系數高值也出現在秋末，但屬于中度變異強度的范圍。土壤含水率的變異系數在不同時期均高于土壤全鹽量的變異系數，說明土壤鹽分的變異性弱于土壤含水率。土壤表層全鹽量的取樣數目大于含水率的合理取樣數目，約為48個，可見測點數目滿足取樣精度。

表2　20 cm深土層不同采樣時間土壤鹽分的統計特征值

2.3半方差函數的結構半方差函數是描述土壤特性空間變異結構的一個函數[5]，其基本方程為：

(2)

式中，h為土壤樣本之間的距離；γ(h)為半方差,用于描述h和γ(h)之間函數關系；Z(xi)為土壤樣本在xi處的測定值;Z(xi+h)為在距離為xi+h處的測定值；n為以h為間距的所有監測點的成對數目(若取樣點有m個, 則n=m-1)。

半方差函數可以通過變差圖來表示。該圖中存在3個重要的參數，分別是基臺值C0(sill)、范圍值α(Range)和塊金值C(Nugget)。上述參數決定了半方差函數的結構。基臺值是指監測點間隔距離不同時半方差函數的穩定值；塊金值是監測間隔為0時的半變異函數值；范圍值是塊金值達到基臺值的橫坐標差值。C0/(C0+C)為樣本特性參數的空間差異性強度，表示可度量空間隨機因素引起的變異所占的比例。當空間差異性強度C0/(C0+C)值小于0.25，屬于強空間相關性；當C0/(C0+C)處于0.25～0.75之間，屬于中等空間相關性；當C0/(C0+C)大于0.75，屬于弱空間相關性[6]。半方差函數的非負性確定了現階段只可采用理論模型進行擬合求參。常用模型有線性模型(Linear)、球狀模型(Spherical)、指數模型(Exponential)、高斯模型(Gaussian)等[5]。球狀模型對擬合淺層土壤含水率和全鹽量是最為常用的[7-10]，因此該研究通過調用Matlab軟件中非線性曲線擬合函數Isqcurvefit對球狀模型進行擬合求參(表3)。

從空間相關性來看，6月和7月土壤全鹽量的C0/(C0+C)值在0.25～0.75之間，屬于中等強度的空間自相關性；5月、8月及9月的C0/(C0+C)值均小于0.25，則具有較強的空間自相關性；9月土壤含水量的C0/(C0+C)值為0.33，其余各月C0/(C0+C)的值均小于0.25。由此可知，在夏季，表層土壤含水率的空間自相關程度要明顯高于土壤全鹽量。隨時間的推移，土樣含水率自相關性逐漸變弱，而土壤全鹽量的空間自相關性呈現出由強到弱再到強的變化態勢。從相關距離上看，不同時期的含水率的自相關范圍均處在26m；土壤全鹽量處于8～25m，監測期土壤含水率的自相關距離大于土壤全鹽量的，說明表層含水率受到外界的影響因素多于表層全鹽量。

表3　20 cm深土層不同采樣時間土壤含水率和全鹽量半方差函數擬合參數

3結論

以寧夏銀北地區的鹽漬土水鹽分布為研究對象，著重分析夏季淺層土壤含水率、全含量的變化特征及空間變異性。研究表明，該地區土壤表層含水率在8月、9月較高，淺層土壤全鹽量隨著時間的遷移而逐漸增加，其值隨著降雨量的增加又會降低；不同時期土壤含水率的變異系數均高于土壤全鹽量，其合理取樣次數低于土壤全鹽量的取樣次數；采用球狀模型擬合半方差函數所得參數表明，表層鹽漬土含水率的空間變異程度高于土壤全鹽量的，鹽分的變異強度以中等變異為主，水分的變異強度以弱變異為主；土壤含水率的自相關距離對受外界因素的影響較大。

參考文獻

[1] 邢旭光,趙文剛, 馬孝義，等.覆膜滴灌條件下棉花根層土壤鹽分時間穩定性研究[J].農業機械學報，2015,46(7):146-154.

[2] 杜歷,周華.淺層地下水利用與鹽堿地治理[M].銀川:寧夏人民出版社,2003.

[3] 馬文敏,康金虎.TDR在土壤鹽分測試中的試驗研究[J].中國農村水利水電，2006(6):10-13.

[4] 管孝艷,王少麗,高占義，等.鹽漬化灌區土壤鹽分的時空變異特征及其與地下水埋深的關系[J].生態學報,2012,32(4):1202-1210．

[5] 徐建華.現代地理學中的數學方法[M].北京:高等教育出版社,2002.

[6] CAMBARDELLA C A，MOORMAN T B，NOVAK J M, et al.Field-scale variability of soil properties in central low a soils[J].Soil Sci Soc Am J,1994, 58: 1501-1511.

[7] 史文娟,馬媛,徐飛，等.不同微尺度膜下滴灌棉田土壤水鹽空間變異特性[J].水科學進展,2014,25(4):287-596.

[8] 于淑會,白璐,劉慧濤，等.河北濱海平原土壤鹽分空間異質性與淋排均質化作用研究[J].土壤,2014,46(6): 1135-1141.

[9] 周在明,張光輝, 王金哲.環渤海微咸水區土壤鹽分及鹽漬化程度的空間格局[J].農業工程學報,2010,26(10):15-20.

收稿日期2015-10-19

作者簡介趙文娟(1980- )，女，寧夏銀川人，講師，博士，從事旱區水資源與調控方面的研究。

基金項目國家自然科學基金資助項目(91230111)；寧夏高等學校科學技術研究項目(NGY2014005)。

中圖分類號S 152.7

文獻標識碼A

文章編號0517-6611(2015)33-204-03