唐學(xué)芬, 趙文舉, 李宗禮,2, 郁 文

(1.蘭州理工大學(xué) 能源與動(dòng)力工程學(xué)院, 甘肅 蘭州 730050; 2.水利部 水利水電規(guī)劃設(shè)計(jì)總院, 北京 100120)

?

壓砂地土壤鹽分空間變異規(guī)律

唐學(xué)芬1, 趙文舉1, 李宗禮1,2, 郁 文1

(1.蘭州理工大學(xué) 能源與動(dòng)力工程學(xué)院, 甘肅 蘭州 730050; 2.水利部 水利水電規(guī)劃設(shè)計(jì)總院, 北京 100120)

摘要：[目的] 研究不同種植年限壓砂地及裸地不同土層土壤鹽分的空間變異特征,為西北干旱地區(qū)土壤鹽漬化改良和利用提供理論依據(jù)。[方法] 以甘肅省景泰縣壓砂地及裸地不同土層土壤鹽分?jǐn)?shù)據(jù)為基礎(chǔ),運(yùn)用經(jīng)典統(tǒng)計(jì)學(xué)和地統(tǒng)計(jì)學(xué)相結(jié)合的方法,分析土壤鹽分統(tǒng)計(jì)特征值,半方差模擬模型及其擬合參數(shù)和等值線分布圖。[結(jié)果] 壓砂地土壤鹽分均值明顯低于裸地,且裸地>老砂地>新砂地>中砂地,裸地與壓砂地各層土壤鹽分的均值差異較大,且土層Ⅳ(30—50 cm)>土層Ⅲ(20—30 cm)>土層Ⅱ(10—20 cm)>土層Ⅰ(0—10 cm);裸地Ⅱ,Ⅲ,Ⅳ層,中砂地各層土壤鹽分屬于弱變異性,其他各地型各土層土壤鹽分都屬于中等變異性。裸地Ⅰ,Ⅲ層,新砂地第Ⅲ層,中砂地第Ⅰ層,老砂地Ⅰ,Ⅱ,Ⅲ層土壤鹽分半方差函數(shù)最佳理論模型為指數(shù)模型,裸地第Ⅱ?qū)訛榫€性模型,其他各地型各土層均為球狀模型。裸地第Ⅱ?qū)油寥利}分具有弱空間相關(guān)性,其余各土層土壤鹽分均具有強(qiáng)相關(guān)性。試驗(yàn)區(qū)土壤鹽分北部高于南部,且各土層土壤鹽分在一定范圍內(nèi)存在空間上的變異性和相關(guān)性。[結(jié)論] 土壤表層壓砂能有效減小土壤鹽分表聚,抑制土壤次生鹽漬化,不同種植年限壓砂地不同土層土壤鹽分在一定范圍內(nèi)存在空間上的變異性和相關(guān)性。

關(guān)鍵詞：壓砂地; 土壤鹽分; 地統(tǒng)計(jì)學(xué); 空間變異性

文獻(xiàn)參數(shù)： 唐學(xué)芬, 趙文舉, 李宗禮, 等.壓砂地土壤鹽分空間變異規(guī)律[J].水土保持通報(bào)，2016,36(3)：18-23.DOI:10.13961/j.cnki.stbctb.2016.03.005

土壤表層壓砂能有效地改善土壤環(huán)境、減少水分蒸發(fā)、增加雨水入滲、蓄水保墑、增溫保溫[1]。土壤作為歷史自然體，受氣候、生物、母質(zhì)、地形、成土?xí)r間等諸多成土因素的影響，具有復(fù)雜性和高度的時(shí)空變異性，不論在大尺度還是在小尺度上，土壤的空間變異性均存在[2]。自然界土壤分布極為復(fù)雜，同一質(zhì)地的土壤在同一平面或不同深度上并不完全均質(zhì)，其它土壤特性參數(shù)在各點(diǎn)的值也不相同，這種土壤特征在空間分布上的非均一性，稱(chēng)為土壤特性的空間變異性。土壤空間變異性的研究方法由最初的經(jīng)典統(tǒng)計(jì)學(xué)方法已經(jīng)發(fā)展到了時(shí)序分析方法、地統(tǒng)計(jì)學(xué)方法、隨機(jī)模擬方法、分形和分維方法。20世紀(jì)70年代，Burgess等[3]將地統(tǒng)計(jì)學(xué)的方法引入土壤科學(xué)研究領(lǐng)域[4]，克服了經(jīng)典的Fisher統(tǒng)計(jì)理論在研究土壤性質(zhì)空間變異性規(guī)律方面的不足。地統(tǒng)計(jì)學(xué)已被證明是分析土壤特性空間分布特征及其規(guī)律的有效方法之一。

自20世紀(jì)80年代以來(lái)，國(guó)內(nèi)外已有很多利用地統(tǒng)計(jì)學(xué)理論對(duì)土壤鹽漬化程度的空間變異以及區(qū)域水、鹽空間分布動(dòng)態(tài)變化的研究。國(guó)內(nèi)學(xué)者趙成義等[5]通過(guò)研究田塊尺度下土壤水分和鹽分的空間變異性，得出土壤水分和鹽分在一定范圍內(nèi)存在空間相關(guān)性。胡克林等[6]通過(guò)對(duì)一塊麥田內(nèi)的98個(gè)觀測(cè)點(diǎn)取樣分析，測(cè)定了2個(gè)時(shí)期的土壤水分和鹽分含量，發(fā)現(xiàn)兩個(gè)時(shí)期的水分和鹽分在一定范圍內(nèi)均存在空間相關(guān)性。國(guó)外學(xué)者M(jìn) Sylla等[7]研究了非洲西部水稻農(nóng)業(yè)生態(tài)系統(tǒng)在不同尺度上的土壤鹽分空間變異性，研究結(jié)果表明地統(tǒng)計(jì)學(xué)可表明影響鹽分空間變異性的復(fù)雜因素。Mohammadi J等[8]運(yùn)用地統(tǒng)計(jì)學(xué)克里格法研究了土壤鹽分空間變異性，得出地統(tǒng)計(jì)學(xué)工具可計(jì)算土壤鹽分分布。鹽堿化是干旱區(qū)土壤的一個(gè)普遍特征，土壤鹽漬化問(wèn)題是制約干旱區(qū)農(nóng)業(yè)生產(chǎn)發(fā)展的主要障礙，也是影響綠洲生態(tài)環(huán)境穩(wěn)定的重要因素[9]。土壤鹽分的分異狀態(tài)在一定程度上反映了土壤耕作層內(nèi)的鹽漬化程度和狀態(tài)，了解其分異對(duì)于指導(dǎo)人們根據(jù)土壤鹽分空間規(guī)律和變化動(dòng)態(tài)進(jìn)行灌溉和排水及制訂防治土壤鹽漬化措施，保證土地質(zhì)量，提高農(nóng)業(yè)產(chǎn)量具有重要意義。目前，國(guó)內(nèi)外學(xué)者利用地統(tǒng)計(jì)學(xué)理論對(duì)土壤鹽漬化程度的空間變異以及水鹽空間分布開(kāi)展了大量的研究，但是針對(duì)西北干旱區(qū)壓砂地土壤特性空間變異研究較少。為此，本文采用半方差函數(shù)分析景泰地區(qū)壓砂地土壤鹽分的空間變異性，旨在為西北干旱地區(qū)土壤鹽漬化改良和利用提供一定的理論依據(jù)。

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

研究區(qū)位于甘肅省景泰縣蘭州理工大學(xué)試驗(yàn)基地附近，地處西部河西走廊東端，甘、蒙、寧3省(區(qū))交界處，黃土高原與騰格里沙漠的過(guò)渡地帶。土壤類(lèi)型主要為洪積灰棕荒漠土和灰鈣土。景泰縣地處季風(fēng)區(qū)與非季風(fēng)區(qū)過(guò)渡地帶，有降水稀少且集中，年內(nèi)分配和地區(qū)分布不均，年際變化大的特點(diǎn)，年均降水量185 mm，多集中在7—9這3個(gè)月，占全年降水量的61.4%，年均蒸發(fā)量3 038 mm，是降水量的16倍。光熱資源豐富，年日照時(shí)約為2 725 h，日照百分率62%，太陽(yáng)年平均輻射約147.8 kcal/cm2，年≥0 ℃的活動(dòng)積溫3 614.8 ℃，≥10 ℃的有效積溫3 038 ℃，無(wú)霜期141 d，年均溫度8.2 ℃，極端最高氣溫36.6 ℃，極端最低氣溫-27.3 ℃。

2材料與方法

2.1土壤采樣與分析

為了研究景泰地區(qū)壓砂地土壤鹽分空間分布特征，選取不同種植年限且為休耕土地的壓砂地，用GPS定位在樣地中心的空地作土壤剖面，確定土壤樣點(diǎn)。壓砂地砂石層厚度均為10 cm，樣地取樣面積為32 m×32 m，矩形取樣，測(cè)點(diǎn)尺寸為1 m×1 m，每1 m2內(nèi)取一點(diǎn)，測(cè)點(diǎn)中心距為4 m，共64個(gè)取樣點(diǎn)，均為單點(diǎn)采樣。分別對(duì)0—10，10—20，20—30和30—50 cm土層進(jìn)行分層采樣。土壤樣品的取樣方法為土鉆取樣法，所用工具為土鉆，樣品的采集重量為60～70 g。為提高以后根據(jù)土壤樣點(diǎn)數(shù)據(jù)分析成圖的精度，樣點(diǎn)必須比較均勻地分布在全工作區(qū)域內(nèi)?，F(xiàn)場(chǎng)調(diào)查和取樣在2013年5—8月和2014年5—8月進(jìn)行。

土壤電導(dǎo)率與土壤含鹽量密切相關(guān)，野外土壤鹽分測(cè)定時(shí)多采用測(cè)量電導(dǎo)率來(lái)表征土壤含鹽特征[10]。對(duì)所采集的土壤樣品進(jìn)行風(fēng)干，過(guò)2 mm篩備用，所有的土樣均制備5∶1水土比浸提液并測(cè)定其電導(dǎo)率。測(cè)定方法是電導(dǎo)儀法。

2.2研究方法

半方差函數(shù)，也稱(chēng)空間變異函數(shù)，是用來(lái)描述區(qū)域化變量結(jié)構(gòu)性和隨機(jī)性并存這一空間特征而提出的，是描述土壤特性空間變異結(jié)構(gòu)的一個(gè)函數(shù)。半方差函數(shù)的重要參數(shù)包括：塊金系數(shù)(Nugget)、基臺(tái)值(sill)、變程(range)，主要用來(lái)表示區(qū)域化變量在一定尺度上的空間變異和相關(guān)程度。

半方差函數(shù)γ(h)的計(jì)算式:

(1)

式中：γ(h)——半方差值； h——樣本間距； N(h)——間距為h的樣本對(duì)數(shù)； Z(xi+h),Z(xi)——xi+h,xi處的測(cè)定值。

壓砂地的劃分為以種植年限為基準(zhǔn)：新砂地

(NGM)種植年限不足10a，中砂地(MGM)種植年限為25～30a，老砂地(OGM)種植年限為45～60a的壓砂地[11]，裸地(CK)為對(duì)照組。

3結(jié)果與分析

3.1土壤鹽分的統(tǒng)計(jì)特征分布

對(duì)土壤鹽分?jǐn)?shù)據(jù)進(jìn)行統(tǒng)計(jì)特征分析是建立土壤鹽分變異模型的前提和基礎(chǔ)，因此對(duì)壓砂地與裸地各土層土壤鹽分?jǐn)?shù)據(jù)按經(jīng)典統(tǒng)計(jì)方法進(jìn)行初步的特征值計(jì)算，其統(tǒng)計(jì)特征值見(jiàn)表1。

表1　土壤鹽分統(tǒng)計(jì)特征值

由表1可得，壓砂地與裸地各土層鹽分的特征參數(shù)值均表現(xiàn)出明顯的差異性。新砂地和中砂地土壤鹽分均值小于1.0 g/kg，屬于非鹽化土；老砂地Ⅰ層土壤鹽分均值介于1.0～2.0 g/kg，屬于輕度鹽化；裸地Ⅰ層、老砂地Ⅱ，Ⅲ，Ⅳ層土壤鹽分均值在2.0～4.0 g/kg之間，屬于中度鹽化；裸地Ⅱ，Ⅲ，Ⅳ層土壤鹽分均值在5.23～7.30 g/kg之間，都大于4.0 g/kg，總體上屬于重度鹽化類(lèi)型。壓砂地土壤鹽分均值明顯低于裸地，且裸地>老砂地>新砂地>中砂地，說(shuō)明壓砂地具有明顯的抑鹽作用，許強(qiáng)等[12]研究認(rèn)為在連續(xù)種植4～5 a后壓砂地的有機(jī)質(zhì)、全氮、速氮等養(yǎng)分達(dá)到最高后不斷下降，含鹽量降低。老砂地土壤鹽分較新砂地和中砂地大，原因在于隨著種植年限的增加，壓砂地退化程度加重，較多的土?；烊肷皩又校勛饔弥饾u減弱，致使土壤水分流失嚴(yán)重，土壤表層返鹽量增大，土壤易發(fā)生鹽堿化。裸地與壓砂地各層土壤鹽分的均值差異較大，且土層Ⅳ>土層Ⅲ>土層Ⅱ>土層Ⅰ，說(shuō)明各層土壤鹽分平均含量在垂直方向上變異較強(qiáng)。裸地Ⅱ，Ⅲ，Ⅳ層，中砂地各層土壤鹽分變異系數(shù)<10%，屬于弱變異性，其他各地型各土層土壤鹽分的變異系數(shù)都在10%～100%之間，屬于中等變異性；各地型土壤鹽分偏度均為正值，各數(shù)據(jù)都呈正偏態(tài)分布；裸地Ⅱ，Ⅳ層，中砂地Ⅱ，Ⅲ層，老砂地Ⅰ，Ⅱ，Ⅲ層土壤鹽分峰度小于3呈低峰態(tài)分布外，其余數(shù)據(jù)正態(tài)分布性不好。變異函數(shù)的計(jì)算一般要求數(shù)據(jù)符合正態(tài)分布，否則可能存在比例效應(yīng)[13]。通過(guò)K-S檢驗(yàn)發(fā)現(xiàn)，除老砂地Ⅲ，Ⅳ層土壤鹽分服從正態(tài)分布外，其他均為非正態(tài)分布。因此，要對(duì)非正態(tài)分布的數(shù)據(jù)進(jìn)行對(duì)數(shù)轉(zhuǎn)換，以消除插值過(guò)程中可能存在的正比例效應(yīng)。

3.2土壤鹽分空間變異特征分析

壓砂地與裸地各層土壤鹽分表現(xiàn)出了空間等方向性。采用GS+9.0軟件分別對(duì)各地型各土層土壤鹽分進(jìn)行分析，其半方差模擬模型及其擬合參數(shù)見(jiàn)表2。以新砂地為例，地統(tǒng)計(jì)學(xué)分析的半方差函數(shù)及其模型擬合結(jié)果如圖1所示。

表2　土壤鹽分半方差函數(shù)模型參數(shù)

圖1　新砂地各土層土壤鹽分半方差

由表2可得，裸地Ⅰ，Ⅲ層土壤鹽分半方差函數(shù)最佳理論模型為指數(shù)模型，Ⅱ?qū)訛榫€性模型，Ⅳ層為球狀模型，各層變程為6.57～16.88 m；新砂地Ⅰ，Ⅱ，Ⅳ層土壤鹽分半方差函數(shù)最佳理論模型為球狀模型，Ⅲ層為指數(shù)模型，各層變程為9.83～12.35 m；中砂地第Ⅰ層土壤鹽分半方差函數(shù)最佳理論模型為指數(shù)模型，其余各層均為球狀模型，各層變程為6.44～8.76 m；老砂地Ⅰ，Ⅱ，Ⅲ層土壤鹽分半方差函數(shù)最佳理論模型為指數(shù)模型，Ⅳ層為球狀模型，各層變程為9.20～35.76 m。各模型決定系數(shù)在0.584～0.998，說(shuō)明理論模型能很好的反應(yīng)土壤鹽分含量的空間結(jié)構(gòu)特征。按照區(qū)域化變量空間相關(guān)程度的分級(jí)標(biāo)準(zhǔn)[14]，裸地第Ⅱ?qū)訅K金效應(yīng)值大于0.75，土壤鹽分具有弱空間相關(guān)性，其余各層土壤鹽分塊金效應(yīng)值均小于0.25，土壤鹽分具有強(qiáng)相關(guān)性。

3.3土壤鹽分空間分布特征分析

為準(zhǔn)確直觀地描述壓砂地與裸地各層土壤鹽分在空間上的分布，在獲得土壤鹽分的半方差函數(shù)模型后，以新砂地及老砂地為例，利用Suffer 8.0軟件進(jìn)行Kriging插值，并繪制空間分布及等值線圖(如圖2所示)。

由圖2可知，新砂地各土層土壤鹽分均呈現(xiàn)出條狀帶和斑塊狀格局。第Ⅰ層的土壤鹽分在0.13～1.15 g/kg之間變化，第Ⅱ?qū)拥耐寥利}分在0.16～2.29 g/kg之間變化，第Ⅲ層的土壤鹽分在0.13～3.22 g/kg之間變化，第Ⅳ層的土壤鹽分在0.19%～4.30%之間變化；從整體上看，試驗(yàn)區(qū)土壤鹽分北部高于南部，且各層土壤鹽分分布在空間上存在著較強(qiáng)的相關(guān)性。從局部上看，第Ⅰ層,第Ⅱ?qū)?第Ⅲ層土壤鹽分具有相似的空間分布規(guī)律，在西北部和東北部較大；第Ⅳ層土壤鹽分與Ⅰ，Ⅱ，Ⅲ層相比，分布較均勻，且各土層土壤鹽分在一定范圍內(nèi)存在空間上的變異性和相關(guān)性。從疏密程度看，各層土壤鹽分空間變異性第Ⅳ層>第Ⅲ層>第Ⅱ?qū)?第Ⅰ層。老砂地第Ⅰ層土壤鹽分介于0.52～3.46 g/kg，第Ⅱ?qū)油寥利}分介于1.03～3.94 g/kg，第Ⅲ層土壤鹽分介于1.06～5.17 g/kg，第Ⅳ層土壤鹽分介于1.57～13.27 g/kg。老砂地土壤鹽分整體上南部高于北部，與新砂地相反，這與隨著種植年限的增加土壤鹽分遷移有關(guān)。

4討論與結(jié)論

(1) 經(jīng)典統(tǒng)計(jì)分析表明，裸地與壓砂地各土層鹽分的特征參數(shù)值均表現(xiàn)出明顯的差異性。新砂地和中砂地屬于非鹽化土；老砂地Ⅰ層屬于輕度鹽化；裸地Ⅰ層、老砂地Ⅱ，Ⅲ，Ⅳ層屬于中度鹽化；裸地Ⅱ，Ⅲ，Ⅳ層屬于重度鹽化類(lèi)型。壓砂地土壤鹽分均值明顯低于裸地，說(shuō)明壓砂地具有明顯的抑鹽作用。裸地與壓砂地各層土壤鹽分的均值差異較大，各層土壤鹽分平均含量在垂直方向上變異較強(qiáng)。裸地Ⅱ，Ⅲ，Ⅳ層，中砂地各層土壤鹽分屬于弱變異性，其他各地型各土層土壤鹽分的變異系數(shù)都屬于中等變異性，且變異性均比較低。

(2) 地統(tǒng)計(jì)學(xué)分析表明，裸地Ⅰ，Ⅲ層土壤鹽分半方差函數(shù)最佳理論模型為指數(shù)模型，Ⅱ?qū)訛榫€性模型，Ⅳ層為球狀模型；新砂地Ⅰ，Ⅱ，Ⅳ層土壤鹽分半方差函數(shù)最佳理論模型為球狀模型，Ⅲ層為指數(shù)模型；中砂地第Ⅰ層土壤鹽分半方差函數(shù)最佳理論模型為指數(shù)模型，其余各層均為球狀模型；老砂地Ⅰ，Ⅱ，Ⅲ層土壤鹽分半方差函數(shù)最佳理論模型為指數(shù)模型，Ⅳ層為球狀模型。裸地第Ⅱ?qū)油寥利}分具有弱空間相關(guān)性，其余各層土壤鹽分具有強(qiáng)相關(guān)性。

(3) 對(duì)新砂地及老砂地運(yùn)用Surfer 8.0繪制土壤含水量等值線分布圖，從Kriging插值的結(jié)果可知，新砂地試驗(yàn)區(qū)土壤鹽分北部高于南部，老砂地土壤鹽分整體上南部高于北部，與新砂地相反，各層土壤鹽分分布在空間上存在著較強(qiáng)的相關(guān)性。新砂地第Ⅰ,Ⅱ,Ⅲ層土壤鹽分具有相似的空間分布規(guī)律，在西北部和東北部較大；第Ⅳ層土壤鹽分與Ⅰ，Ⅱ，Ⅲ層相比，分布較均勻，且各土層土壤鹽分在一定范圍內(nèi)存在空間上的變異性和相關(guān)性。

[參考文獻(xiàn)]

[1]王平,謝成俊,陳娟.不同種植年限砂田水鹽變化與砂田退化初探[J].水土保持通報(bào),2012,32(2):251-254.

[2]劉繼龍,馬孝義,付強(qiáng),等.不同土層土壤特性空間變異性關(guān)系的聯(lián)合多重分形研究[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2012,43(5):37-42.

[3]Burgess T M, Webster R. Optimal interpolation and isarithmic mapping of soil properties: The semivariogram and punctual Kriging[J]. Journal of Soil Science, 1980,31(2):315-341.

[4]Webster R. Quantitative spatial analysis of soil in the field[M]∥Advances in Soil Science. Springer New York, 1985:1-70.

[5]趙成義,王玉潮,李子良,等.田塊尺度下土壤水分和鹽分的空間變異性[J].干旱區(qū)研究,2003,20(4):252-256.

[6]胡克林,李保國(guó),陳德立,等.農(nóng)田土壤水分和鹽分的空間變異性及其協(xié)同克立格估值[J].水科學(xué)進(jìn)展,2001,12(4):460-466.

[7]Sylla M, Stein A, van Breemen N, et al. Spatial variability of soil salinity at different scales in the mangrove rice agro-ecosystem in West Africa[J]. Agriculture, Ecosystems and Environment, 1995,54(1):1-15.

[8]Mohammadi J. Study of the spatial variability of soil salinity in Ramhormoz Area (Khuzestan) using geostatistical theory 1: Kriging[J]. JWSS-Isfahan University of Technology, 1999,2(4):49-64.

[9]Li Jiangguo, Pu Lijie, Han Mingfang, et al. Soil salinization research in China: Advance and prospects[J]. Journal of Geographical Science, 2014,24(5):943-960.

[10]姚榮江,楊勁松,韓建均.海涂圍墾區(qū)土壤鹽分空間變異模擬的比較研究[J].土壤學(xué)報(bào),2012,49(2):275-281.

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

[12]許強(qiáng),吳宏亮,康建宏,等.旱區(qū)砂田肥力演變特征研究[J].干旱地區(qū)農(nóng)業(yè)研究,2009,27(1):37-41.

[13]楊勁松,姚榮江,劉光明,等.黃河三角洲地區(qū)突然鹽分的空間變異性及其CoKriging估值[J].干旱區(qū)研究,2006,23(3):439-445.

[14]佘冬立,邵明安,俞雙恩.黃土區(qū)農(nóng)草混合利用坡面土壤水分空間變異性[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2010,41(7):57-63.

收稿日期：2015-07-03修回日期：2015-10-09

通訊作者：趙文舉(1981—),男(漢族),甘肅省永昌縣人,博士,副教授,主要從事農(nóng)業(yè)水土工程方面的研究。E-mail:wenjuzhao@126.com。

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1000-288X(2016)03-0018-06

中圖分類(lèi)號(hào)：S156.4

Variation Pattern of Soil Salinity of Gravel-sand Mulched Field

TANG Xuefen1, ZHAO Wenju1, LI Zongli1,2, YU Wen1

(1.SchoolofEnergyandPowerEngineering,LanzhouUniversityofTechnology,Lanzhou,Gansu730050,China; 2.GeneralInstituteforWaterResourcesandHydropowerPlanningandDesign,MinistryofWaterResources,Beijing100120,China)

Abstract：[Objective] We studied the variation pattern of soil salinity at different layer of GSM(gravel-sand mulched field) and CK in order to provide basis for improvement and utilization of soil salinization in the northwest arid area in China. [Methods] Based on the statistics at different layer of GSM and CK, classical statistics and geo-statistics were used to analyze the statistical characteristic value, semi-variance simulation model and its fitting parameters, and the isocline distribution map of soil salinity in Jingtai County of Gansu Province. [Results] Mean soil salinity of GSM was lower than that of CK, and the result showed: CK>OGM(old gravel-sand mulched field)>NGM(new gravel-sand mulched field)>MGM(middle gravel-sand mulched field). There were significant differences in soil salinity among different layers of GSM and CK, and the result showed: layer Ⅳ(30—50 cm)>layer Ⅲ(20—30 cm)>layerⅡ(10—20 cm)>layerⅠ(0—10 cm). The variations of soil salinity in layerⅡ, Ⅲ, Ⅳ(CK) and MGM were weak, while variations in other layers were moderate. The best semi-variance theory model of layer Ⅱ(CK) was linear model, the optimal model for layerⅠ, Ⅲ(CK), layer Ⅲ(NGM), layer Ⅰ(MGM) and layer Ⅰ, Ⅱ, Ⅲ(OGM) was exponential model, while the best model for other layers was spherical model. The spatial auto-correlation of soil salinity at layer Ⅱ(CK) was lower, while it was high at other layers. The soil salinity in northern area was higher than that in southern area. [Conclusion] Gravel sand on soil surface can reduce soil salt accumulation and control soil secondary salinization, and there is spatial variability and correlation exists in soil salinity among different layer of GSM and CK.

Keywords:gravel-sand mulched field; soil salinity; geo-statistics; spatial variability

資助項(xiàng)目：國(guó)家自然科學(xué)基金項(xiàng)目“壓砂地土壤水分時(shí)空變異理論與模型研究”(51269008); 蘭州理工大學(xué)紅柳青年教師培養(yǎng)計(jì)劃資助項(xiàng)目(Q201310)

第一作者：唐學(xué)芬(1991—),女(漢族),甘肅省永靖縣人,碩士研究生,研究方向?yàn)檗r(nóng)業(yè)水土工程。E-mail:tang1202@163.com。