不同采樣幅度和間距下壓砂地棗樹土壤水分的空間變異性研究

趙文舉，崔珍，馬孝義，馬宏，于平

（1. 蘭州理工大學(xué)能源與動力工程學(xué)院,甘肅 蘭州 730050;

2. 西北農(nóng)林科技大學(xué)旱區(qū)農(nóng)業(yè)水土工程教育部重點實驗室,陜西 楊凌 712100）

不同采樣幅度和間距下壓砂地棗樹土壤水分的空間變異性研究

趙文舉1，崔珍1，馬孝義2，馬宏1，于平1

（1. 蘭州理工大學(xué)能源與動力工程學(xué)院,甘肅 蘭州 730050;

2. 西北農(nóng)林科技大學(xué)旱區(qū)農(nóng)業(yè)水土工程教育部重點實驗室,陜西 楊凌 712100）

趙文舉, 崔珍, 馬孝義, 馬宏, 于平. 不同采樣幅度和間距下壓砂地棗樹土壤水分的空間變異性研究[J]. 農(nóng)業(yè)現(xiàn)代化研究, 2016, 37(6): 1181-1189.

Zhao W J, Cui Z, Ma X Y, Ma H, Yu P. Effects of sampling extent and spacing on spatial variability of soil water content in a field with gravel-mulched jujube tree[J]. Research of Agricultural Modernization, 2016, 37(6): 1181-1189.

為了探明西北地區(qū)特有的壓砂地棗樹土壤水分的空間變異及其尺度效應(yīng)，本文基于野外試驗，選取32 m×32 m區(qū)域，并在此基礎(chǔ)上改變采樣幅度和采樣間距，基于經(jīng)典統(tǒng)計學(xué)和地統(tǒng)計學(xué)理論，研究了不同采樣幅度和間距條件下0-50 cm土層土壤含水量的空間分布特征及其空間變異性。結(jié)果表明，對于所有5種采樣幅度（32 m×32 m、28 m×28 m、24 m×24 m、20 m×20 m和16 m×16 m），隨土層深度的增加，土壤含水量呈減小的趨勢，而變異系數(shù)呈增大的趨勢，空間變異強度基本表現(xiàn)為弱變異至中等偏弱變異；當(dāng)采樣幅度增大時，土壤含水量的變異系數(shù)Cv、塊金值C0及變程A均不斷增大。對于4 m、8 m和12 m這3種采樣間距，當(dāng)采樣間距增大時，土壤含水量的塊金值C0不斷增大，變程A不斷減小，而變異系數(shù)Cv不受影響。在不同尺度內(nèi)，土壤含水量均存在強烈的空間自相關(guān)性。各土層土壤含水量在不同采樣間距下的空間分布形態(tài)相似，多處出現(xiàn)明顯的“隆起”與“凹陷”，受地形影響顯著，并隨著采樣間距的增大，逐漸平坦化，8 m為較合理的采樣間距。

壓砂地；土壤含水量；尺度效應(yīng)；空間分布

土壤水分是土壤-植物-大氣連續(xù)體（soil-plantatmosphere continuum，SPAC系統(tǒng)）中大氣降水、地表水、土壤水和地下水相互轉(zhuǎn)化的一個重要環(huán)節(jié)[1]。土壤水分的空間格局及空間變異是空間合理配置和人工植被恢復(fù)的科學(xué)基礎(chǔ)。隨季節(jié)和研究尺度不同，土壤水分表現(xiàn)出較強的空間變異性[2-5]。土壤水分的空間異質(zhì)性主要受到土壤質(zhì)地、植被、地形地貌、降雨及人類活動等非線性因子的影響[6]。有研究表明，不同尺度條件下對土壤水分空間變異規(guī)律的研究可為提高土壤水分利用率、建立水分管理信息系統(tǒng)以及發(fā)展和驗證遙感反演提供基礎(chǔ)的地面觀測數(shù)據(jù)支持，是實現(xiàn)田間“精準灌溉”的基礎(chǔ)[7-8]。Bl?schl[9]認為，采樣尺度的變化將使土壤水分的表觀變異與真實變異出現(xiàn)偏差。胡偉等[10]通過改變采樣幅度和間距的“再采樣”方法，研究了尺度大小對土壤水分空間變異的影響，得出在一定范圍內(nèi)隨采樣幅度的增大，土壤水分變異的特征參數(shù)都呈不同程度地增大；隨采樣間距的增大，相關(guān)指數(shù)沒有變化，而相關(guān)距離卻減小。魏新光等[11]研究指出土壤含水量時空變異性顯著，垂直方向上隨著土層深度的增加，變異系數(shù)逐漸降低。邢旭光等[12]認為對于40 m×40 m、8 m×8 m和2 m×2 m 3種尺度，0-60 cm深度各土層土壤含水量在水平方向上的變異強度表現(xiàn)為弱變異至中等（偏弱）變異，且隨尺度減小和土層深度增加而減小，且所有取樣點0-60 cm深度內(nèi)土壤含水量在垂直方向上的變異強度表現(xiàn)為弱變異至中等（偏弱）變異。趙文舉等[13]研究得出不同種植年限壓砂地的變異系數(shù)隨著采樣幅度的增大而增大。史文娟等[14]研究了3種微尺度（0.25 m、1 m、4 m）條件下棉田土壤水分的空間變異性，指出土壤水分變異性為中等變異，半方差函數(shù)大部分可采用高斯模型模擬。同時也有研究表明[15-17]，微尺度條件下的土壤含水量半方差函數(shù)模型和空間變異特性與傳統(tǒng)的大中尺度之間存在較大差異，因此，深入了解尺度變化對土壤水分空間變異的影響，定量分析土壤水分空間變異特征及其尺度效應(yīng)，針對特定地區(qū)進行的微尺度土壤特性的空間變異性研究亟需進一步開展。

壓砂地是西北旱區(qū)勞動人民利用砂石作為覆蓋材料，依賴天然降水創(chuàng)造出的一種以砂石覆蓋和免耕為核心的保護性耕作模式，距今已有300多年的發(fā)展歷史[18]。土壤表層覆蓋砂石具有明顯的增滲、減蒸、保溫、抗蝕、減小地表徑流等作用[19-20]，但會使得土壤水分的入滲、蒸發(fā)機理比均質(zhì)土特殊復(fù)雜。同時鑒于土壤水分分布狀況對植被生長有至關(guān)重要的影響[21-22]，對土壤水分空間變異性的研究具有重要意義。為此，基于野外試驗，擬對西北地區(qū)特有的壓砂地棗樹土壤水分在水平方向和垂直方向上的空間變異特性及分布特征進行微尺度研究，以期為壓砂地田間水分管理提供理論依據(jù)。

1 研究材料與方法

1.1 試驗區(qū)概況

研究區(qū)在景泰縣蘭州理工大學(xué)試驗研究基地附近。地處河西走廊東端，北依寧夏，黃土高原與騰格里沙漠的過渡地帶。土壤類型主要為洪積灰棕荒漠土和灰鈣土。景泰縣地處季風(fēng)區(qū)與非季風(fēng)區(qū)過渡地帶，縣境內(nèi)氣候呈現(xiàn)出明顯的大陸性氣候特征，年均降水量185 mm，多集中在7-9三個月，占全年降水量的61.4%，年均蒸發(fā)量3 038 mm，是降水量的16倍。無霜期141 d，年均溫度8.2 ℃。

1.2 試驗方法

用GPS定位在樣地中心的空地作土壤剖面，確定土壤樣點。在東西方向與南北方向，采用規(guī)格網(wǎng)格布點法，以網(wǎng)格最小尺寸為32 m×32 m的矩形網(wǎng)格取樣作為每個試驗塊，網(wǎng)格的尺度取4 m，測點中心距為4 m，共64個取樣點，采樣時將表層砂石小心鏟除，取樣深度為0-10 cm，10-20 ㎝，20-30 cm和30-50 cm。土樣的取樣方法為土鉆取樣法，樣品的采集重量為60-70 g，采用烘干法測定土壤含水量。采樣點均勻地分布在研究區(qū)域內(nèi)，取樣點分布如圖1所示?，F(xiàn)場調(diào)查和取樣在2015年10月14日進行。

圖1 土壤含水量取樣點分布圖Fig. 1 Sampling point distribution of soil water content

1.2.1 改變采樣幅度 對所有的測定數(shù)據(jù)進行分析，并分別以28 m×28 m、24 m×24 m、20 m×20 m、16 m×16 m的“滑動窗”從測定區(qū)域的西北角向東南角滑動，計算每個滑動窗內(nèi)的變異參數(shù)。并利用Sufer8.0對已有土壤含水量數(shù)據(jù)通過克里格插值分別繪出其空間分布圖。

1.2.2 改變采樣間距 首先對所有的數(shù)據(jù)進行空間變異分析，然后在原始采樣點的基礎(chǔ)上分別在東西和南北向每隔一個和兩個點對原始采樣點進行隔行抽取，來實現(xiàn)采樣密度的改變，即增大采樣間距，并使在每種間距下，所有的原測數(shù)據(jù)都取到，并將相同采樣點數(shù)的樣本變異參數(shù)的平均值作為該采樣密度內(nèi)的變異結(jié)果。需要指出的是，在地統(tǒng)計學(xué)分析時，如果間距太大可能造成數(shù)據(jù)太少，從而使數(shù)據(jù)失去可靠性，因此只用采樣間距為4、8、12 m這3種情況分析采樣間距對土壤含水量空間變異的影響。

1.3 數(shù)據(jù)處理與分析

數(shù)據(jù)的空間變異研究均采用經(jīng)典統(tǒng)計學(xué)和地統(tǒng)計學(xué)相結(jié)合的方法。采用SPSS20.0進行經(jīng)典統(tǒng)計計算，采用GS+9.0對土壤含水率進行半方差分析和理論模型的擬合，采用克里格法對各尺度0-50 cm土層土壤含水量進行插值，最后借助Sufer8.0和Origin9.0完成圖形繪制。

1.3.1 變異函數(shù)的理論模型擬合 球狀模型：

式中：r（h）為半方差值；h為樣本間距；C0為塊金常數(shù)；C0+C為基臺值；C為拱高；a為變程。

指數(shù)模型：

注：3 a為變程。

高斯模型：

半方差函數(shù)模型中塊金值（C0）反映隨機因素引起的空間變異程度，越接近于0表明空間相關(guān)性越強；基臺值（C0+C）反映區(qū)域化變量在研究范圍內(nèi)總的空間變異程度；變程（A）反映區(qū)域化變量空間自相關(guān)變異的尺度范圍，在變程之內(nèi)具有空間自相關(guān)性，反之，則是獨立的；基臺比C0/(C0+C)為空間相關(guān)度，反映可度量的空間隨機因素引起的變異所占的比例，當(dāng)C0/(C0+C)＜0.25，表示變量具有強烈的空間自相關(guān)性，0.25≤C0/(C0+C)≤0.75，表示變量具有中等的空間自相關(guān)性，C0/(C0+C)＞0.75，表示變量空間自相關(guān)性很弱[24]。

2 結(jié)果與分析

2.1 采樣幅度對空間變異結(jié)果的影響

2.1.1 不同采樣幅度尺度條件下土壤含水量統(tǒng)計特征值分析 變異系數(shù)Cv與采樣幅度的關(guān)系如圖2所示。由圖2可知，對于16 m×16 m、20 m×20 m、24 m×24 m、28 m×28 m和32 m×32 m 5種采樣幅度0-50 cm土層內(nèi)土壤水分在水平方向上的變異系數(shù)均值分別為13.0%、13.3%、14.0%、14.8%和16.5%，呈現(xiàn)出隨采樣幅度增大而增大的趨勢。這主要是由于隨著研究區(qū)域的增大不斷引入新的變異因素的緣故，即某些在較大尺度內(nèi)對土壤水分分布起影響作用的因素在較小有限尺度內(nèi)，可能具有相對一致性；當(dāng)尺度增大時，它對水分分布的影響才得到體現(xiàn)，從而使得其變異系數(shù)增大，變異增大。對變異系數(shù)Cv與采樣幅度的關(guān)系進行擬合，得到式（4）：

式中：y為變異系數(shù)；x為采樣幅度。由擬合結(jié)果可知，在研究區(qū)域內(nèi)，所選尺度均較小，故變異系數(shù)差別不大；隨著采樣幅度x的增大，變異系數(shù)y值也隨之增大。

圖2 采樣幅度對土壤水分變異系數(shù)的影響Fig. 2 The effect of scale for soil moisture variation coefficient

通過對16 m×16 m到32 m×32 m采樣幅度的壓砂地棗樹區(qū)水平方向的土壤含水量和0-50 cm土層垂直方向的土壤含水量統(tǒng)計分析（表1），對于不同采樣幅度，0-50 cm土層含水率均表現(xiàn)出隨深度增加而減小的趨勢；對于全部取樣點，各土層土壤含水量變異系數(shù)在水平方向上呈現(xiàn)隨深度增加而增加的趨勢，介于8.1%-22.5%，除0-10 cm土層空間變異強度在24 m×24 m、20 m×20 m及16 m×16 m尺度時表現(xiàn)為弱變異外，其余均表現(xiàn)為中等偏弱變異。

2.1.2 不同采樣幅度條件下土壤含水量空間變異特性分析 以地統(tǒng)計學(xué)方法為基礎(chǔ)，對不同采樣幅度尺度條件下壓砂地棗樹區(qū)土壤水分進行半方差分析，繪制采樣幅度為32 m×32 m時的土壤水分半方差函數(shù)圖（圖3）。由圖3可以看出，采樣幅度為32 m×32 m時的土壤水分半方差函數(shù)曲線變化較平穩(wěn)，模型擬合的相關(guān)系數(shù)較高，表明在整個研究區(qū)域內(nèi)土壤水分的相關(guān)性較高。各采樣幅度尺度的半方差函數(shù)模型及相關(guān)參數(shù)見表2。

從表2可以看出，在32 m×32 m尺度條件下，各土層（除10-20 cm外）均采用指數(shù)模型模擬；在28 m×28 m尺度條件下，各土層的最佳理論模型均為指數(shù)模型；在24 m×24 m尺度條件下，各土層可采用指數(shù)、高斯和球狀模型模擬；在20 m×20 m尺度條件下，各土層可采用球狀和指數(shù)模型模擬；在16 m×16 m尺度條件下，各土層可采用高斯和球狀模型進行模擬。采樣幅度從32 m×32 m減小到16 m×16 m時，C0均值分別為0.394、0.356、0.333、0.187和0.083，表現(xiàn)為隨采樣幅度的減小，塊金值不斷減小的趨勢，這是因為隨著采樣幅度的減小，短距離內(nèi)的變異和測量誤差不斷減小的緣故。變程（A）均值分別為22.00 m、18.54 m、15.62 m、14.80 m和8.74 m，表明空間自相關(guān)距離隨采樣幅度減小而減小，且最小變程為8.12 m，大于采樣間距4 m，表明采樣間距合理?；_值(C0+C)呈現(xiàn)出隨土層深度增加而增大的趨勢，表明在研究區(qū)域內(nèi)，50 cm土層土壤含水量的空間變異程度也與土層深度密切相關(guān)。在各采樣幅度尺度內(nèi)，基臺比C0/(C0+C)大部分小于0.25，其余介于0.25-0.50之間，表明在各研究幅度尺度內(nèi)，0-50 cm土層土壤含水量整體表現(xiàn)出較強的空間自相關(guān)性。

表1 不同采樣幅度0-50 cm深度各土層水平方向上土壤水分統(tǒng)計特征值Table 1 Statistics parameters of horizontal soil water content in 0-50 cm layers under different sampling extents

圖3 32 m×32 m時的土壤水分半方差函數(shù)圖Fig. 3 Semi-variance function diagram of soil water content of the scale of 32 m×32 m

表2 不同采樣幅度0-50 cm深度各土層土壤水分半方差函數(shù)模型及參數(shù)擬合Table 2 Semi-variance models and fitted parameters of soil water content in 0-50 cm layers under different sampling extents

2.2 采樣間距對空間變異結(jié)果的影響

2.2.1 不同采樣間距條件下土壤含水量統(tǒng)計特征值分析 當(dāng)采樣間距從4 m變化至12 m時，對壓砂地棗樹區(qū)水平方向的土壤含水量和0-50 cm土層垂直方向的土壤含水量進行統(tǒng)計分析（表3）。當(dāng)采樣密度發(fā)生變化時，各土層土壤含水量和變異系數(shù)均圍繞一個固定值波動。其變異系數(shù)的均值分別為11.4%、14.2%、19.4%和23.2%。由此可以看出，研究區(qū)域內(nèi)土壤含水量變異系數(shù)受土層深度的影響，而采樣間距對變異系數(shù)基本不產(chǎn)生影響，這也說明采樣幅度一定時，采樣密度的改變并不能改變影響土壤含水量變異的因素。由此可見，在一定的研究區(qū)域內(nèi)，適當(dāng)增大采樣間距仍可以得到土壤水分的實際變異系數(shù)。

2.2.2 不同采樣間距條件下土壤含水量空間變異特性分析 在不同采樣間距下，0-10 cm土層采用指數(shù)模型模擬；10-20 cm土層采用球狀模擬；20-30 cm土層和30-50 cm土層采用指數(shù)模型模擬（表4）。

隨采樣間距從4 m增大到12 m時，塊金值（C0）基本呈現(xiàn)增大的趨勢，這可能是由于短距離內(nèi)微地域或植被間差異雖表現(xiàn)出一定的結(jié)構(gòu)特征，但卻被較大距離上影響土壤水分變異過程的局部微地貌、土壤質(zhì)地等因素所掩蓋?；_值無明顯變化規(guī)律，變程（A）隨采樣間距的增大呈明顯減小的趨勢，說明適當(dāng)減小采樣間距能使土壤含水量表現(xiàn)出更好的空間連續(xù)性，采樣間距為4 m和8 m時，變程均大于相應(yīng)的采樣間距，表明所取采樣間距合理，而當(dāng)采樣間距為12 m時，變程均小于采樣間距，采樣間距不再合理，說明本研究中8 m的采樣間距是較為合理的，可以在取相對較少采樣點的情況下仍可獲得較為合理的結(jié)果。不同采樣間距時，各土層基臺比C0/（C0+C）大部分小于0.25，其余介于0.25-0.50之間，表明當(dāng)采樣間距不同時，0-50 cm土層內(nèi)土壤含水量整體表現(xiàn)出較強的空間自相關(guān)性。

表3 不同采樣間距0-50 cm深度各土層水平方向上土壤水分統(tǒng)計特征值Table 3 Statistics parameters of horizontal soil water content in 0-50 cm layers of different sampling spacings

表4 采樣間距對0-50 cm土層土壤水分半方差參數(shù)的影響Table 4 Effects of sampling spacing on semi-variogram parameter of soil water content in 0-50 cm layers

2.2.3 不同采樣間距條件下土壤含水量空間分布特征 為便于直觀了解研究區(qū)域內(nèi)土壤含水量的空間分布特征，分別繪制了各采樣間距尺度條件下0-50 cm深度內(nèi)各土層土壤含水量在水平方向和豎直方向的空間分布圖，如圖4所示。

由圖4可知，各尺度土壤水分均表現(xiàn)出“凹凸不平”的分布情況，這可能與取樣區(qū)域內(nèi)的地形地貌有關(guān)，這與邢旭光等[12]的研究結(jié)果相似。在不同采樣間距條件下，0-20 cm土層土壤含水量均較20-50 cm土層波動幅度大，原因在于淺層土壤易受蒸發(fā)和降水作用影響，由于表層（0-10 cm）土壤極易受人為因素影響，導(dǎo)致10 cm深度土壤含水量分布出現(xiàn)多處明顯“隆起”與“凹陷”，而30-50 cm深度土壤含水量分布則相對“平坦”。由圖4進一步可知，各土層土壤含水量在不同采樣間距下的空間分布總體趨勢相同，但隨著采樣間距的增大，其分布形態(tài)趨于“平坦”。采樣間距為8 m時，其圖形雖趨于平坦化，但與采樣間距為4 m的空間分布圖差別較小，可以較好的表征土壤含水率的空間分布特征，而當(dāng)采樣間距為12 m時，相較于前兩個采樣間距，其分布特征過于平坦化，不再能表征土壤含水率的實際空間變異特征。所以本研究中采樣間距可取為8 m，這可以在獲得可靠的實驗結(jié)果的基礎(chǔ)上大大減小采樣的工作量。

3 討論

本文通過改變采樣幅度和采樣間距探討了不同采樣尺度對壓砂地棗樹區(qū)0-50 cm土層土壤水分空間變異性的影響。研究結(jié)果表明，隨著土層深度的增加，土壤含水量呈減小的趨勢，而變異系數(shù)呈增大的趨勢，這與張繼光等[25]、姚雪玲等[26]的研究結(jié)果相同；而邢旭光等[12]在對獼猴桃果園土壤含水量進行研究時得出，在0-60 cm土層內(nèi)隨土層深度增加，土壤含水量先增大后減小，變異系數(shù)逐漸減小的結(jié)論；魏新光等[11]對丘陵區(qū)棗林土壤水分動態(tài)進行了研究，得出隨土層深度增加，變異系數(shù)逐漸降低的結(jié)論；趙文舉等[13]在對不同種植年限壓砂地0-50 cm土層進行研究時得出隨著采樣幅度尺度的增大，土壤含水量和變異系數(shù)均增大的結(jié)論。這些不同的變化規(guī)律可能與各研究區(qū)域地形地貌、管理措施、土壤質(zhì)地以及植被類型和種植結(jié)構(gòu)等有關(guān)，這也決定了不同研究區(qū)域土壤水分和變異系數(shù)的差異性。

圖4 不同采樣間距0-50 cm土層土壤含水量空間分布圖Fig. 4 Spatial distribution of soil water content in 0-50 cm layers under different sampling spacings

3.1 采樣幅度對空間變異的影響

郭德亮等[15]在對不同采樣幅度黑河中游綠洲區(qū)不同土地利用類型表層土壤水分空間變異進行研究時，林地在濕潤條件下土壤水分變異的特征參數(shù)顯示出與本研究相似的規(guī)律，即塊金值（C0）和變程（A）均隨采樣幅度的增大而增大，而農(nóng)地的變程（A）則無固定的變化規(guī)律；胡偉等[10]對黃土高原退耕坡地土壤水分的空間變異進行尺度性研究時得出在一定范圍內(nèi)隨采樣幅度的增大，土壤水分變異的特征參數(shù)都呈不同程度地增大的結(jié)論，與本研究中變異系數(shù)的變化規(guī)律相同；而邢旭光等[12]在對不同采樣幅度下獼猴桃果園土壤含水量進行研究時得出，變程（A）隨采樣幅度的增大而增大的結(jié)論，這與本研究結(jié)果一致，而塊金值（C0）較穩(wěn)定且變化程度減小的結(jié)論卻與本研究中塊金值隨采樣幅度增大而增大的結(jié)論有所差異。變程（A）隨采樣幅度增大而增大的現(xiàn)象可能與土壤水分分布格局的等級巢式結(jié)構(gòu)有關(guān)[10]，即隨采樣幅度的增大，大尺度上結(jié)構(gòu)因素引起的變異相關(guān)性會掩蓋小尺度上結(jié)構(gòu)因素和隨機因素引起的變異相關(guān)性，使土壤水分在較大尺度內(nèi)變化平穩(wěn)。這些相似的和不同的變化規(guī)律可能與各研究區(qū)域氣候和環(huán)境條件不同有關(guān)。通過比較不同采樣幅度下土壤水分的變異結(jié)果，我們發(fā)現(xiàn)在采樣和模擬時需要根據(jù)實際需求調(diào)整采樣幅度的大小。較大的采樣幅度可以改善土壤的分類系統(tǒng)，提高土壤的調(diào)查和測繪質(zhì)量，中、小尺度可以為合理制定作物布局、有效管理土地資源及提高土壤利用率進而實現(xiàn)精確灌溉奠定基礎(chǔ)。因此，在今后研究中需要根據(jù)研究區(qū)實際情況和研究目的確定采樣幅度，不斷總結(jié)不同區(qū)域在不同采樣幅度下土壤水分的空間變異規(guī)律，為微尺度采樣幅度下土壤水分的合理監(jiān)測提供理論依據(jù)。

3.2 采樣間距對空間變異的影響

有關(guān)采樣間距對土壤含水量空間變異特性的研究也較多，而研究結(jié)果也有一定的差異性。Western和Bl?schl[27]的研究指出變程（A）隨采樣間距的增大而增大；胡偉等[10]對黃土高原退耕坡地土壤水分的空間變異進行尺度性研究時得出隨采樣間距的增大，相關(guān)距離卻減小的結(jié)論，這與本研究中得出的結(jié)論一致，這可能與土地類型和模型的擬合有關(guān)。因此，有關(guān)采樣間距對壓砂地棗樹土壤含水量空間變異特征的影響有待進一步研究。同時，Western和Bl?schl[27]指出，在空間變異研究中，非常小的采樣體積、非常大的采樣幅度和采樣密度是最理想的采樣尺度，只有樣本數(shù)據(jù)足夠多才可以代表研究區(qū)域的真實變異特征。但實際應(yīng)用中受土壤水分測定時人力和物力的影響，采樣幅度和采樣密度不可能無限增大。

因此，本文中對壓砂地棗樹區(qū)土壤水分空間變異的尺度性研究，在節(jié)省資源的前提下可為確定最優(yōu)的采樣幅度和密度提供有效的參考。但是，本文只對特定時間里土壤水分的空間變異性規(guī)律進行了探討，如要全面了解景泰地區(qū)壓砂地棗樹土壤水分的空間變異規(guī)律，還需要對土壤水分在不同尺度上的空間分布隨時間的變化規(guī)律做更進一步的研究。

4 結(jié)論

1）在不同采樣幅度條件下，土壤含水量隨土層深度增加而減?。蛔儺愊禂?shù)（Cv）隨土層深度和采樣幅度的增加均呈現(xiàn)出增大的趨勢，表現(xiàn)為弱變異或中等偏弱變異；塊金值（C0）基本表現(xiàn)為隨采樣幅度的減小而不斷減小的趨勢；變程（A）隨采樣幅度減小而減小；基臺值（C0+C）呈現(xiàn)出隨土層深度增加而增大的趨勢，且土壤含水量存在強烈的空間自相關(guān)性。

2）在不同采樣間距條件下，變異系數(shù)變化不大，這說明采樣幅度一定時，采樣密度的改變并不能改變影響土壤含水量變異的因素。隨采樣間距的增大，塊金值（C0）基本呈現(xiàn)增大的趨勢；變程（A）隨采樣間距的增大呈明顯減小的趨勢，且8 m為較合理的采樣間距，可在減少采樣點的情況下仍可獲得較為合理的結(jié)果；土壤含水量整體表現(xiàn)出較強的空間自相關(guān)性。

3）在不同采樣間距條件下，各土層土壤含水量的空間分布均表現(xiàn)出“凹凸不平”的形態(tài)，0-20 cm土層土壤含水量較20-50 cm土層波動幅度大，多處出現(xiàn)明顯“隆起”與“凹陷”。各土層土壤含水量在不同采樣間距下的空間分布總體趨勢相同，但隨著采樣間距的增大，其分布形態(tài)趨于“平坦”。8 m的采樣間距，仍可以較好的表征土壤含水率的空間分布特征，是較為合理的采樣間距。

[1] 高曉東, 吳普特, 張寶慶, 等. 黃土丘陵區(qū)小流域土壤有效水空間變異及其季節(jié)性特征[J]. 土壤學(xué)報, 2015, 52(1): 57-67.

Gao X D, Wu P T, Zhang B Q, et al. Spatial variability of available soil moisture and its seasonality in a small watershed in the hilly region of the loess plateau[J]. Acta Pedologica Sinica, 2015, 52(1): 57-67.

[2] 宗路平, 角媛梅, 李石華, 等. 哈尼梯田景觀水源區(qū)土壤水分時空變異性[J]. 生態(tài)學(xué)雜志, 2015, 34(6): 1650-1659.

Zong L P, Jiao Y M, Li S H, et al. Spatial and temporal variability of soil moisture in water source region of Hani terrace landscape[J]. Chinese Journal of Ecology, 2015, 34(6): 1650-1659.

[3] Brocca L, Tullo T, Melone F, et al. Catchment scale soil moisture spatial-temporal variability[J]. Journal of Hydrology, 2012, 422/423: 63-75.

[4] Gao X D, Wu P T, Zhao X N, et al. Soil moisture variability along transects over a well-developed gully in the Loess Plateau, China[J]. Fuel & Energy Abstract, 2011, 87(3): 357-367.

[5] 王艷莉, 劉立超, 高艷紅, 等. 人工固沙植被區(qū)土壤水分動態(tài)及空間分布[J]. 中國沙漠, 2015, 35(4): 942-950.

Wang Y L, Liu L C, Gao Y H, et al. Dynamic and spatial distribution of soil moisture in an artificially re-vegetated desert area[J]. Journal of Desert Research, 2015, 35(4): 942-950.

[6] 李小雁. 干旱地區(qū)土壤-植被-水文耦合、響應(yīng)與適應(yīng)機制[J].中國科學(xué):地球科學(xué), 2011, 41(12): 1721-1730.

Li X Y. Mechanism of coupling, response and adaptation betweensoil, vegetation and hydrology in arid and semiarid regions[J]. Science China: Earth Science, 2011, 41(12): 1721-1730.

[7] Navas A, Maehin J. Spatial distribution of heavy metals and arsenic in soils of Aragón (Northeast Spain): Controlling factors and environmental implications[J]. Applied Geochemistry, 2002, 17(8): 961-973.

[8] Wang J K, Zhao Y C, Zhang X D, et al. Studies on the spatial variability of heavy metal contents in black soils in central areas of Hailun[J]. Chinese Journal of Soil Science, 2003, 34(5): 398-403.

[9] Bl?schl G. Scale and Scaling in Hydrology: A Framework for Thinking and Analysis[M]. Chichester: Wiley, 1999.

[10] 胡偉, 邵明安, 王全九. 黃土高原退耕坡地土壤水分空間變異的尺度性研究[J]. 農(nóng)業(yè)工程學(xué)報, 2005, 21(8): 11-16.

Hu W, Shao M A, Wang Q J. Scale-dependency of spatial variability of soil moisture on a degraded slope-land on the Loess Plateau[J]. Transactions of the Chinese Society of Agricultural Engineering, 2005, 21(8): 11-16.

[11] 魏新光, 聶真義, 劉守陽, 等. 黃土丘陵區(qū)棗林土壤水分動態(tài)及其對蒸騰的影響[J]. 農(nóng)業(yè)機械學(xué)報, 2015, 46(6): 130-140.

Wei X G, Nie Z Y, Liu S Y, et al. Soil moisture characteristics and its influence on jujube tree transpiration in Loess hilly region[J]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(6): 130-140.

[12] 邢旭光, 趙文剛, 柳燁, 等. 獼猴桃果園不同采樣密度下土壤含水量空間變異性研究[J]. 農(nóng)業(yè)機械學(xué)報, 2015, 46(8): 138-145.

Xing X G, Zhao W G, Liu Y, et al. Spatial variability of soil moisture in Kiwi field under different sampling density conditions[J]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(8): 138-145.

[13] 趙文舉, 李娜, 李宗禮, 等. 不同種植年限壓砂地土壤水分空間變異規(guī)律研究[J]. 農(nóng)業(yè)現(xiàn)代化研究, 2015, 36(6): 1067-1073.

Zhao W J, Li N, Li Z L, et al. Study on spatial variability of soil water content in gravel-mulched fields with different planting years[J]. Research of Agricultural Modernization, 2015, 36(6): 1067-1073.

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

Shi W J, Ma Y, Xu F, et al. Spatial variability of soil moisture and salt content in cotton field on microscales under mulch drip irrigation[J]. Advances in Water Science, 2014, 25(4): 585-593.

[15] 郭德亮, 樊軍, 米美霞. 黑河中游綠洲區(qū)不同土地利用類型表層土壤水分空間變異的尺度效應(yīng)[J]. 應(yīng)用生態(tài)學(xué)報, 2013, 24(5): 1199-1208.

Guo D L, Fan J, Mi M X. Scale-dependency of spatial variability of surface soil moisture under different land use types in Heihe oasis, China[J]. Chinese Journal of Applied Ecology, 2013, 24(5): 1199-1208.

[16] 郭麗俊, 李毅, 李敏, 等. 塿土土壤水力特性空間變異的多重分形分析[J]. 農(nóng)業(yè)機械學(xué)報, 2011, 42(9): 50-58.

Guo L J, Li Y, Li M, et a1.Multifractal study on spatial variability of soil hydraulic properties of Lou soil[J]. Transactions of the Chinese Society for Agricultural Machinery, 2011, 42(9): 50-58.

[17] 楊磊, 衛(wèi)偉, 陳利頂, 等. 黃土丘陵溝壑區(qū)深層土壤水分空間變異及其影響因子[J]. 生態(tài)與農(nóng)村環(huán)境學(xué)報, 2012, 28(4): 355-362.

Yang L, Wei W, Chen L D, et a1. Spatial variation of deep soil water in the loess hilly and gully regions and its influencing factors[J]. Journal of Ecology and Rural Environment, 2012, 28(4): 355-362.

[18] 王占軍, 蔣齊, 何建龍, 等. 寧夏環(huán)香山地區(qū)壓砂地土壤水分特征曲線及入滲速率的特征分析[J]. 土壤通報, 2013, 44(6): 1364-1368.

Wang Z J, Jiang Q, He J L, et al. Analysis of soil water characteristic curve and infiltration rates of gravel-mulched land around Xiangshan Mountain area in Ningxia[J]. Chinese Journal of Soil Science, 2013, 44 (6): 1364-1368.

[19] Li X Y. Gravel-sand mulch for soil and water conservation in the of Northwest China[J]. Catena, 2003, 52(2): 105-127.

[20] Qiu Y, Wang Y J, Xie Z K. Long-term gravel-sand mulch affects soil physicochemical properties, microbial biomass and enzyme activities in the semi-arid Loess Plateau of North-western China[J]. Acta Agriculturae Scandinavica, Section B-Soil & Plant Science, 2014, 64(4): 294-303.

[21] Xie Z K, Wang Y J, Cheng G D, et al. Particle-size effects on soil temperature, evaporation, water use efficiency and watermelon yield in fields mulched with gravel and sand in semi-arid Loess Plateau of Northwest China[J]. Agricultural Water Management, 2010, 97(6): 917-923.

[22] 馬波, 田軍倉. 壓砂地西瓜水肥耦合模型及優(yōu)化組合方案[J].干旱地區(qū)農(nóng)業(yè)研究, 2010, 28(4): 24-29.

Ma B, Tian J C. Model of coupling water with fertilizer in gravelmulched watermelon field and its optimum combination scheme[J]. Agricultural Research in the Arid Areas, 2010, 28(4): 24-29, 35.

[23] 雷志棟, 楊詩秀, 謝森傳. 土壤水動力學(xué)[M]. 北京: 清華大學(xué)出版社, 1988.

Lei Z D, Yang S X, Xie S Z. Soil Water Dynamics[M]. Beijing: Tsinghua University Press, 1988.

[24] 王云強, 邵明安, 劉志鵬. 黃土高原區(qū)域尺度土壤水分空間變異性[J]. 水科學(xué)進展, 2012, 23(3): 310-316.

Wang Y Q, Shao M A, Liu Z P. Spatial variability of soil moisture at a regional scale in the Loess Plateau[J]. Advances in Water Science, 2012, 23(3): 310-316.

[25] 張繼光, 陳洪松, 蘇以榮, 等. 喀斯特峰叢洼地坡面土壤水分空間變異研究[J]. 農(nóng)業(yè)工程學(xué)報, 2006, 22(8): 54-58.

Zhang J G, Chen H S, Su Y R, et al. Spatial variability of soil moisture on hillslope in cluster-peak depression areas in Karst region[J]. Transactions of the Chinese Society of Agricultural Engineering, 2006, 22(8): 54-58.

[26] 姚雪玲, 傅伯杰, 呂一河. 黃土丘陵溝壑區(qū)坡面尺度土壤水分空間變異及影響因子[J]. 生態(tài)學(xué)報, 2012, 32(16): 4961-4968.

Yao X L, Fu B J, Lü Y H. Spatial patterns of soil moisture at transect scale in the Loess Plateau of China[J]. Acta Ecologica Sinica, 2012, 32(16): 4961-4968.

[27] Western A W, Bl?schl G. On the spatial scaling of soil moisture[J]. Journal of Hydrology, 1999, 217(3/4): 203-224.

（責(zé)任編輯：王育花）

Effects of sampling extent and spacing on spatial variability of soil water content in a field with gravel-mulched jujube tree

ZHAO Wen-ju1, CUI Zhen1, MA Xiao-yi2, MA Hong1, YU Ping1
（1.College of Energy and Power Engineering, Lanzhou University of Technology, Lanzhou, Gansu 730050, China; 2. Key Laboratory for Agricultural Soil and Water Engineering in Arid Area of Ministry of Education, Northwest A&F University, Yangling, Shaanxi 712100, China）

To explore the effects of sampling extent and spacing on spatial variability of soil moisture in a gravelsand mulched jujube orchard, spatial analyses including classical statistics and geo-statistics were conducted with soil water content (SWC) of 0-50 cm in the 32 m×32 m field by changing the sampling extent and spacing. Irrespective of sampling extents (i.e. 32 m×32 m, 28 m×28 m, 24 m×24 m, 20 m×20 m and 16 m×16 m), the SWC decreased with the increase of soil depth while the variation coefficient increased with the increasing soil depth. The spatial variability of SWC was weak and moderate for all scales; The variation coefficient (Cv), nugget (C0), and variation range (A) of SWC all increased with the increase of sampling extent. When three sampling spacings (i.e. 4, 8 and 12 m) were concerned, the nugget (C0) increased with the increase of sampling spacing; The variation range (A) decreased with the increase of sampling spacing; the Cvwas not affected by the sampling spacing. A strong spatial autocorrelation characteristic of SWC was found for all scales. The distribution patterns of SWC in the same soil layer at different sampling spacings were similar and tended to be "flat" with the increase of the sampling spacing. The reasonable sampling spacing was found to be 8 m in our case.

gravel-sand mulched field; soil water content; scale effect; spatial distribution

ZHAO Wen-ju, E-mail: wenjuzhao@126.com.

10.13872/j.1000-0275.2016.0114

S152.7

A

1000-0275（2016）06-1181-09

國家自然科學(xué)基金項目（51269008，51279167）；蘭州理工大學(xué)紅柳青年教師培養(yǎng)計劃項目(Q201310)。

趙文舉（1981-），男，甘肅永昌人，博士，副教授，主要從事農(nóng)業(yè)水土工程方面的研究，E-mail：wenjuzhao@126.com。

2016-04-15，接受日期：2016-09-07

Foundation item: National Natural Science Foundation of China (51269008, 51279167); HongLiu Young Teachers Training Program of Lanzhou University of Technology (Q201310).

Received 15 April, 2016;Accepted 7 September, 2016