張冬明，張 文，鄭道君，吉清妹，吳宇佳，潘孝忠，符傳良，王鐘友，謝良商（海南省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)環(huán)境與土壤研究所/農(nóng)業(yè)部海南耕地保育科學(xué)觀測(cè)試驗(yàn)站，海口 571100）

?

海水倒灌農(nóng)田土壤鹽分空間變異特征①

張冬明，張 文，鄭道君，吉清妹，吳宇佳，潘孝忠，符傳良，王鐘友，謝良商*
（海南省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)環(huán)境與土壤研究所/農(nóng)業(yè)部海南耕地保育科學(xué)觀測(cè)試驗(yàn)站，?？?571100）

摘 要：以海南省?？谑腥r(nóng)場(chǎng)為研究區(qū)域，采集并化驗(yàn)海水倒灌后土壤樣品 191份，利用地統(tǒng)計(jì)學(xué)方法和空間插值法對(duì)區(qū)域土壤鹽分空間變異特征進(jìn)行了研究。結(jié)果表明：①4個(gè)土層中土壤鹽分的均值最大的為0 ～ 5 cm土層，土壤鹽分平均含量達(dá)到9.63 g/kg，最小的是20 ～ 40 cm土層，鹽分平均含量?jī)H為4.36 g/kg；②4個(gè)土層中第一層（0 ～ 5 cm）土壤鹽分含量屬弱相關(guān)性，第四層（20 ～ 40 cm）土層屬于劇烈空間自相關(guān)性；③各土層土壤鹽分分布格局差異明顯，含鹽量 ＜ 1 g/kg的各土層均沒(méi)有分布；含鹽量1 ～ 3和3 ～ 5 g/kg兩個(gè)級(jí)別隨著土層深度的增加面積逐漸增大，＞ 9 g/kg的級(jí)別則面積逐漸減少。不同深度土壤層土壤鹽分含量差異明顯，空間變異特征各不相同，總體呈現(xiàn)土壤鹽分表聚特征。

關(guān)鍵詞：土壤鹽分；地統(tǒng)計(jì)學(xué)；空間變異；三江農(nóng)場(chǎng)；海水倒灌

海南島屬熱帶島嶼季風(fēng)性氣候，沿海岸線有不少面積的農(nóng)田分布，受東北和西南季風(fēng)影響，熱帶風(fēng)暴和臺(tái)風(fēng)頻繁，海水倒灌耕地現(xiàn)象時(shí)有發(fā)生；尤其是受較強(qiáng)臺(tái)風(fēng)影響的時(shí)候，海水倒灌現(xiàn)象明顯，受害農(nóng)田面積大，土壤恢復(fù)耕性周期長(zhǎng)。土壤鹽分過(guò)高對(duì)作物的危害大致有3點(diǎn)：①造成作物的生理干旱；②鹽分離子的毒害；③破壞作物的正常代謝。土壤鹽分的空間變異一定程度上反映了耕層土壤鹽漬化狀況，土壤資源的科學(xué)合理利用是建立在正確掌握區(qū)域土壤屬性空間變異特征基礎(chǔ)之上的［1］。研究土壤鹽分空間變異特征是土壤改良利用和鹽漬化防治的基礎(chǔ)，因此，研究和摸清海水倒灌農(nóng)田土壤表層土壤鹽分含量及分布特征對(duì)于指導(dǎo)人們根據(jù)土壤鹽分動(dòng)態(tài)變化規(guī)律和空間變異特性進(jìn)行改良和防治，提升耕地質(zhì)量，提高農(nóng)業(yè)效益和農(nóng)業(yè)的可持續(xù)發(fā)展乃至為海南的綠色崛起具有重要意義［2-4］。

近年來(lái)，隨著3S技術(shù)的迅猛發(fā)展，許多國(guó)內(nèi)外學(xué)者運(yùn)用地統(tǒng)計(jì)學(xué)對(duì)土壤鹽分的空間變異性進(jìn)行了大量研究，取得了一定的成果。Jordan和 Navarro［5］對(duì)土壤鹽分在干旱和半干旱地區(qū)地質(zhì)和環(huán)境因素影響下的空間變異進(jìn)行了研究；李亮亮等［6］對(duì)地統(tǒng)計(jì)學(xué)基本原理及在土壤物理特性、土壤養(yǎng)分、土壤鹽分、土壤重金屬等土壤空間變異研究中的應(yīng)用進(jìn)行了論述；姚榮江等［7］針對(duì)黃河三角洲地區(qū)存在的土壤鹽漬化問(wèn)題，運(yùn)用經(jīng)典統(tǒng)計(jì)學(xué)和地統(tǒng)計(jì)學(xué)相結(jié)合的方法研究了不同深度土層鹽分含量的空間變異特征，并繪制了各土層鹽分的隨機(jī)性和結(jié)構(gòu)性的半方差圖以及空間分布圖；王勇輝等［8］分析了艾比湖主要補(bǔ)給河流下游河岸帶土壤鹽分特征；張同娟等［9］運(yùn)用電磁感應(yīng)儀EM38，結(jié)合GIS技術(shù)和地統(tǒng)計(jì)方法對(duì)長(zhǎng)江河口地區(qū)典型地段土壤含鹽量的空間變異性進(jìn)行研究；孫運(yùn)朋等［10］對(duì)濱海棉田土壤鹽分時(shí)空分布特征進(jìn)行了研究；冉啟洋等［11］對(duì)塔里木河上游綠洲土壤表層鹽分特征進(jìn)行了研究，認(rèn)為在研究區(qū)內(nèi)鹽分指標(biāo)影響最大的是Cl-和Na+。本研究通過(guò)野外調(diào)查采樣、室內(nèi)分析化驗(yàn)、綜合運(yùn)用地統(tǒng)計(jì)軟件，對(duì)海水倒灌三江農(nóng)場(chǎng)土壤鹽分空間特征進(jìn)行分析，以揭示鹽漬化狀況和空間分布規(guī)律，為當(dāng)?shù)佧}漬化土壤的綜合提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

?？谑腥r(nóng)場(chǎng)，位于海南島北部地區(qū)，土地總面積5 680 hm2，其中，陸域面積4 654 hm2，淺海灘涂地1 026 hm2。屬濱海平原地形，地勢(shì)東南高，西北低，西至北為大洋田，地形平坦；東南部地形為玄武巖和淺海沉積物形成的臺(tái)地。主要土壤為水稻土、磚紅壤、濱海鹽土和潮汐土四大類。氣候?qū)贌釒Ъ撅L(fēng)氣候，年平均氣溫 23.4℃，年降水量1 560 mm，光照充足適合農(nóng)作物生長(zhǎng)。目前，三江農(nóng)場(chǎng)水稻田以稻菜輪作模式為主，即春夏季種植水稻、秋冬季種植瓜菜，而旱坡地則以發(fā)展花卉和熱帶水果為主。據(jù)2014年6月采集的農(nóng)田土壤分析結(jié)果可知，農(nóng)場(chǎng)農(nóng)田土壤肥力中等，主要養(yǎng)分含量為：水溶性鹽0.92 g/kg、pH 5.10、有機(jī)質(zhì)1.97 g/kg、堿解氮 85.63 mg/kg、有效磷 13.52 mg/kg、速效鉀55.14 mg/kg。2014年7月18日，受超強(qiáng)臺(tái)風(fēng)“威馬遜”和天文大潮的共同影響，導(dǎo)致農(nóng)場(chǎng)三分之一的農(nóng)田倒灌海水，加上水利設(shè)施的落后，倒灌后的海水不能及時(shí)排掉，造成多數(shù)田洋土壤鹽分含量超標(biāo)，農(nóng)作物根系腐爛變黑，水稻嚴(yán)重減產(chǎn)或絕收，直接影響到當(dāng)?shù)剞r(nóng)民的吃飯問(wèn)題。

1.2 土壤樣品采集與分析

本次采樣是在2014年第九號(hào)超強(qiáng)臺(tái)風(fēng)“威馬遜”過(guò)后的 8月 4—8日，以 2014年 2月 4日拍攝的MODIS衛(wèi)星影像數(shù)據(jù)為參考，結(jié)合三江農(nóng)場(chǎng)土地利用現(xiàn)狀圖，應(yīng)用GPS定位技術(shù)布設(shè)采樣點(diǎn)（圖1），觀測(cè)點(diǎn)的布置綜合考慮了當(dāng)?shù)赝恋乩梅绞?、土壤類型和農(nóng)田水利情況等因素，并在此基礎(chǔ)上盡可能規(guī)則分布，以便于進(jìn)行地統(tǒng)計(jì)分析?？偣膊荚O(shè)了191個(gè)采樣點(diǎn)，相鄰兩個(gè)采樣點(diǎn)之間距離不超過(guò)1.3 km，每個(gè)樣點(diǎn)用鋤頭挖出一個(gè)微型土壤剖面，分層采集土壤樣品，取樣深度為0 ～ 5、5 ～ 10、10 ～ 20和20 ～ 40 cm，每層取3個(gè)樣品，混合均勻后用四分法取大約500 g的土樣帶回實(shí)驗(yàn)室。土壤樣品室內(nèi)分析方法，采用文獻(xiàn)［12］的方法進(jìn)行分析化驗(yàn)，化驗(yàn)前先對(duì)土壤樣品進(jìn)行自然風(fēng)干、磨碎、過(guò)篩處理，制備1:5的土水質(zhì)量比浸提液，測(cè)定土壤全鹽含量。

1.3 數(shù)據(jù)處理

采用SPSS16.0軟件對(duì)研究區(qū)全鹽的數(shù)據(jù)進(jìn)行經(jīng)典統(tǒng)計(jì)分析，運(yùn)用GS+7.0軟件進(jìn)行半方差函數(shù)計(jì)算，而全鹽空間插值以及分布圖繪制均采用 ArcGIS 9.3進(jìn)行。

2 結(jié)果與分析

2.1 不同土層土壤鹽分的統(tǒng)計(jì)特征

由表1中偏度和峰度系數(shù)可以判斷，土壤10 ～ 20、20 ～ 40 cm兩個(gè)土層的土壤鹽分含量均不符合正態(tài)分布類型，在變異函數(shù)結(jié)構(gòu)分析中一般需要將數(shù)據(jù)轉(zhuǎn)化為符合正態(tài)分布或近似正態(tài)分布的數(shù)據(jù)類型，避免產(chǎn)生比例效益［13］，本研究中10 ～ 20、20 ～ 40 cm兩個(gè)土層的土壤鹽分含量經(jīng)對(duì)數(shù)轉(zhuǎn)換后符合正態(tài)分布；而0 ～ 5和5 ～ 10 cm兩個(gè)土壤鹽分含量符合正態(tài)分布類型，無(wú)需進(jìn)行數(shù)據(jù)轉(zhuǎn)換處理。4個(gè)土層中土壤鹽分的均值差異較大，最大的為0 ～ 5 cm土層，土壤鹽分平均含量達(dá)到9.63 g/kg，最小的是20 ～ 40 cm土層，鹽分平均含量?jī)H為4.36 g/kg，鹽分含量大小順序依次為：0 ～ 5 cm＞5 ～ 10 cm＞10 ～ 20 cm＞20 ～ 40 cm，即隨土層深度的增加土壤鹽分含量逐漸降低。從平均值來(lái)看，受本次海水倒灌的農(nóng)田表層土壤鹽漬化較為嚴(yán)重，從單個(gè)土壤采樣點(diǎn)的鹽分含量來(lái)看，仍有近12.5% 的采樣點(diǎn)土壤鹽分超過(guò)了20 g/kg。從變異系數(shù)來(lái)看，4個(gè)土層均超過(guò)了 100%，屬于強(qiáng)度變異性［14］，說(shuō)明受海水倒灌影響，土壤鹽分的差異性較大，近海地區(qū)海拔相對(duì)較低，海水浸泡時(shí)間長(zhǎng)，滲入到土壤表層的土壤鹽分離子相對(duì)多；反之亦然。

圖1 研究區(qū)土壤采樣點(diǎn)分布圖

2.2 半變異函數(shù)分析

半方差函數(shù)模型的選取步驟是：首先根據(jù)公式計(jì)算出Y（h）-h的散點(diǎn)圖；然后用不同類型的模型進(jìn)行擬合，以選擇離差平方最小為原則選擇最佳模型，最后用交叉驗(yàn)證法來(lái)修正所選模型的參數(shù)［15］。本研究得到的半方差函數(shù)模型參數(shù)見(jiàn)表2。

表1 土壤鹽分的描述性統(tǒng)計(jì)特征

表2 土壤鹽分變異函數(shù)模型參數(shù)

一般，半方差函數(shù)值都是隨著樣點(diǎn)間距的增加而增大，并在一定的間距 （稱為變程，arrange） 升大到一個(gè)基本穩(wěn)定的常數(shù) （稱為基臺(tái)，sill）；C0/C0+C為塊金比（塊金效應(yīng)），表示空間異質(zhì)性程度，即空間自相關(guān)程度。如果比值＜25%，說(shuō)明系統(tǒng)具有劇烈空間自相關(guān)性，空間變異主要受結(jié)構(gòu)性因素的影響；如果比值在25% ～ 75% 之間，表明系統(tǒng)具有中等空間自相關(guān)性，空間變異是結(jié)構(gòu)性因素和隨機(jī)性因素共同作用的結(jié)果；＞75% 說(shuō)明系統(tǒng)空間自相關(guān)性很弱，空間變異主要是隨機(jī)因素引起。由表2可知，5 ～ 10 cm 和10 ～ 20 cm兩個(gè)土層土壤鹽分含量屬于中等空間自相關(guān)性，即空間變異是由結(jié)構(gòu)性因素（成土母質(zhì)、地形、地下水位和土壤類型等）和隨機(jī)性因素（耕作模式、施肥方法等）共同作用的結(jié)果；0 ～ 5 cm土壤鹽分含量屬弱相關(guān)性，主要受隨機(jī)因素引起的（如天氣和農(nóng)事活動(dòng)等），即可把臺(tái)風(fēng)天氣的海水倒灌理解成一次人為不合理措施；而20 ～ 40 cm深土層屬于劇烈空間自相關(guān)性，由于未直接與外界接觸，其土壤鹽分含量主要受地下水位和地表水滲漏的共同影響。從變程來(lái)看，20 ～ 40 cm土層土壤鹽分趨于一個(gè)穩(wěn)定數(shù)值的樣點(diǎn)間距最短，僅139 m；此外，5 ～ 10 cm和10 ～ 20 cm土壤鹽分含量的半方差函數(shù)模型均為球形模型。

2.3 空間插值分析

為了更直觀地反映整個(gè)農(nóng)場(chǎng)土壤鹽分的空間分布格局，依據(jù)2.2節(jié)所得到的半方差函數(shù)模型，采用Kriging最優(yōu)內(nèi)插法，編繪了研究區(qū)不同深度土壤鹽分含量的空間分布圖。

圖 2顯示了各層土壤不同含鹽量范圍的面積分布情況?？傮w而言，隨著土壤深度的增加土壤含鹽量逐漸下降，第一層（0 ～ 5 cm）土壤含鹽量相對(duì)較高，含量＞9 g/kg的土壤有3 938.16 hm2，占農(nóng)場(chǎng)總面積的36.39%；第二層（5 ～ 10 cm）土壤含鹽量較第一層有所降低，尤其是＞9 g/kg的區(qū)域有較大幅度的減少，只剩西部沿海地區(qū)有連片分布，此層土壤鹽分含量集中在 5 ～ 7 g/kg，面積為 3 737.61 hm2，占總面積的65.86%；從第三層（10 ～ 20 cm）土壤鹽分的插值圖可以看出，在西部沿海呈現(xiàn)一個(gè)三角形狀的較高含鹽量（7 ～ 9 g/kg）區(qū)域，面積為1 710.27 hm2，占總面積的30.13%，此區(qū)域的南部地區(qū)含鹽量低于北部地區(qū)；第四層土壤鹽分含量表現(xiàn)出與第三層土壤類似的變化規(guī)律，即以本農(nóng)場(chǎng)西南沿海地區(qū)為中心，由西南向東北階梯狀下降的圓弧變化規(guī)律，只有極少面積的土壤含鹽量＞9 g/kg，面積僅為24.38 hm2，僅占總面積的0.43%，絕大部分土壤含鹽量＜5 g/kg。各土層的不同含鹽量面積分布情況見(jiàn)表3。

同樣，從表3還可以看出，土壤含鹽量＜1 g/kg的各土層均沒(méi)有分布。含鹽量1 ～ 3和3 ～ 5 g/kg兩個(gè)級(jí)別隨著土層深度的增加，面積逐漸增大。0 ～ 5 和5 ～ 10 cm兩個(gè)土層中含鹽5 ～ 7 g/kg的級(jí)別面積增加最大，第二層較第一層增加了3 204.90 hm2；面積減少最多的是土壤含鹽量 ＞ 9 g/kg的級(jí)別，第二層較第一層減少了3 872.09 hm2。

3 結(jié)論與討論

采集的?？谑腥r(nóng)場(chǎng)土壤剖面 4個(gè)土層中土壤鹽分的均值差異較大，最大的為0 ～ 5 cm土層，土壤鹽分平均含量達(dá)到9.63 g/kg，最小的是20 ～ 40 cm土層，鹽分平均含量?jī)H為4.36 g/kg，鹽分含量大小順序依次為：0 ～ 5 cm＞5 ～ 10 cm＞10 ～ 20 cm＞20 ～ 40 cm，表層土壤第一、二層次（0 ～ 5，5 ～ 10 cm）的鹽分含量服從正態(tài)分布；第三、四層土壤（10 ～ 20，20 ～ 40 cm）的鹽分含量不服從正態(tài)分布。4個(gè)層次土壤鹽分含量的變異系數(shù)均超過(guò)了100%，屬于強(qiáng)變異程度。

圖2 土壤鹽分含量空間插值分布圖

表3 不同深度土層各鹽分級(jí)別面積統(tǒng)計(jì)表

從空間自相關(guān)性來(lái)看，第一層（0 ～ 5 cm）土壤鹽分含量屬弱相關(guān)性，第四層（20 ～ 40 cm）深土層屬于劇烈空間自相關(guān)性；第四層（20 ～ 40 cm）變程最短，僅139 m；5 ～ 10 cm和10 ～ 20 cm土壤鹽分含量的半方差函數(shù)模型均為球形模型。

4個(gè)土層鹽分含量分布格局差異明顯，含鹽量＜1 g/kg的各土層均沒(méi)有分布。含鹽量1 ～ 3和3 ～ 5 g/kg兩個(gè)級(jí)別隨著土層深度的增加，面積逐漸增大。0 ～ 5和5 ～10 cm兩個(gè)土層中含鹽5 ～ 7 g/kg的級(jí)別面積增加最大，第二層較第一層增加了3 204.90 hm2；面積減少最多的是土壤含鹽量 ＞ 9 g/kg的級(jí)別，第二層較第一層減少了3 872.09 hm2。

一般認(rèn)為，土壤鹽分的空間分布特征是自然條件和人為因素綜合作用的結(jié)果，但主導(dǎo)因素各地區(qū)有所不同，受土地利用方式和地形地貌的影響較大。本研究中，海水倒灌是最主要因素，由于海拔較低又是沿海地區(qū)，地下水也是不可忽略的因素；同時(shí)，研究區(qū)內(nèi)的水利設(shè)施欠完善，也是導(dǎo)致本次鹽害不能及時(shí)排除或降低的重要因素之一，7—10月是海南的雨季，這4個(gè)月集中了全年60% 的降雨量，按常理經(jīng)過(guò)這幾個(gè)月的降雨洗鹽作用，土壤鹽分應(yīng)該有較大幅度下降，但是研究人員幾次到實(shí)地調(diào)研，發(fā)現(xiàn)多數(shù)田洋還是積水，但土壤鹽分含量還有待第二次采樣調(diào)查，在此未能進(jìn)行研究討論。

參考文獻(xiàn)：

［1］ 連綱， 郭旭東， 傅伯杰， 等.黃土高原小流域土壤養(yǎng)分空間變異特征及預(yù)測(cè)［J］.生態(tài)學(xué)報(bào)， 2008， 28（3）: 946-954

［2］ 姚榮江， 楊勁松， 姜龍， 等.基于聚類分析的土壤鹽漬剖面特征及其空間分布研究［J］.土壤學(xué)報(bào)， 2008， 45（1）: 56-65

［3］ 管孝艷， 王少麗， 高占義， 等.鹽漬化灌區(qū)土壤鹽分的時(shí)空變異特征及其與地下水埋深的關(guān)系［J］.生態(tài)學(xué)報(bào)，2012， 32（4）: 1 002-1 010

［4］ 田長(zhǎng)彥， 周宏飛， 劉國(guó)慶.21世紀(jì)新疆土壤鹽漬化調(diào)控與農(nóng)業(yè)持續(xù)發(fā)展研究建議［J］.干旱區(qū)地理， 2000， 23（2）: 178-181

［5］ Jordan M M， Navarro P J.spatial dynamics of soil salinity under arid and semiarid conditions［J］.Environmental Geology， 2004， 45: 448-456

［6］ 李亮亮， 依艷麗， 凌國(guó)鑫， 等.地統(tǒng)計(jì)學(xué)在土壤空間變異研究中的應(yīng)用［J］.土壤通報(bào)， 2005， 36（2）: 265-268

［7］ 姚榮江， 楊勁松， 劉廣明， 等.黃河三角洲地區(qū)典型地塊土壤鹽分空間變異特征研究［J］.農(nóng)業(yè)工程學(xué)報(bào)， 2006，22（6）: 61-66

［8］ 王勇輝， 馬蓓， 海米提，等.艾比湖主要補(bǔ)給河流下游河岸帶土壤鹽分特征［J］.干旱區(qū)研究， 2013， 30（2）: 196-201

［9］ 張同娟， 楊勁松， 劉廣明.長(zhǎng)江河口地區(qū)典型地段土壤鹽分特征研究［J］.灌溉排水學(xué)報(bào)， 2010， 29（3）: 32-35

［10］ 孫運(yùn)朋， 陳小兵， 張振華， 等.濱海棉田土壤鹽分時(shí)空分布特征研究［J］.土壤學(xué)報(bào)， 2013， 50（5）: 891-899

［11］ 冉啟洋， 貢璐， 韓麗， 等.塔里木河上游綠洲土壤表層鹽分特征［J］.中國(guó)沙漠， 2013， 33（4）: 1 098-1 103

［12］ 鮑士旦.土壤農(nóng)化分析［M］.北京: 中國(guó)農(nóng)業(yè)出版社，2005

［13］ Yao R J， Yang J S， Quntitative analysis of spatial distribution patern of soil salt accumulation in plough layer and shallow groundwater in the Yellow River Delta［J］.Transactions of the CSAE.， 2007， 23（8）: 45-51

［14］ 胡克林， 李保國(guó)， 林啟美， 等.農(nóng)田土壤養(yǎng)分的空間變異性特征［J］.農(nóng)業(yè)工程學(xué)報(bào)， 1999， 15（3）: 33-38

［15］ 顏安， 蔣平安， 盛建東， 等.瑪納斯河流域表層土壤鹽分空間變異特征研究［J］.土壤學(xué)報(bào)， 2014， 52（2）: 410-414

中圖分類號(hào)：S156

DOI:10.13758/j.cnki.tr.2016.03.030

基金項(xiàng)目：①海南省耕地改良專項(xiàng)（HNGDhs201501）和海南省科研院所技術(shù)開(kāi)發(fā)專項(xiàng)（KYYS-2015-12，KYYS-2016-02）資助。

* 通訊作者（lshxie@163.com）

作者簡(jiǎn)介：張冬明（1982—），男，江西井岡山人，碩士，副研究員，主要從事土壤改良與土壤信息系統(tǒng)研究。E-mail: purvery@163.com

Spatial Variability of Surface Soil Salinity of Farmland by Seawater

ZHANG Dongming， ZHANG Wen， ZHENG Daojun， JI Qingmei， WU Yujia， PAN Xiaozhong ，

FU Chuanliang，WANG Zhongyou， XIE Liangshang*
（Agricultural Environment and Soil Research Institute of Hainan Academy of Agricultural Sciences /Scientific Observing and Experimental Station of Arable Land Conservation （Hainan）， Ministry of Agriculture， Haikou 571100， China）

Abstract:Characteristics of the soil salinity of the SanJiang farm in Haikou city was studied by using geostatistics and spatial interpolation methods on the bases of collecting and measuring soil samples with seawater intrusion.The results showed that the maximum mean soil salinity was 9.63 g/kg happened in 0-5 cm layer， and the minimum one was 4.36 g/kg， happened in 20-40 cm layer.The 0-5 cm layer was a weak correlation in soil salt content while the 20-40 cm layer was intense spatial correlation.The distribution pattern of soil salinity had obvious difference between different layers， and no layer with less than 1 g/kg of salt content.With the increase of soil depth， the areas of 1-3 g/kg and 3-5 g/kg salt contents increased gradually， while the area of greater than 9 g/kg salt content reduced.The differences of soil salt contents in different soil depths were obvious， and spatial variability were also different， generally there is a surface accumulation of salt.

Key words:Soil salinity； Geostatistics； Spatial variation； Sanjiang farm； Seawater intrusion