劉洋+姜德鋒+鄒曉霞

摘要：黃河三角洲地區(qū)是我國重要的后備耕地資源，明確土壤鹽分的空間變異特征是充分開發(fā)該地區(qū)土地資源的基礎。本試驗選取黃河三角洲地區(qū)典型土地利用類型棉田，沿黃河入海方向自內(nèi)陸至海岸布設采樣帶，研究0～100 cm土層土壤鹽分在剖面及空間的分布特征。結(jié)果發(fā)現(xiàn)，沿黃河入海方向自西向東土壤鹽分呈上升趨勢，與成土時間、距海遠近和人類活動密切相關。整個采樣區(qū)（利津縣北宋鎮(zhèn)至河口區(qū)孤島鎮(zhèn)）棉田土壤鹽分以輕鹽土為主，受人為耕作管理影響0～60 cm土壤條件趨于一致，空間變異性小于60～100 cm土層；研究區(qū)土壤鹽分剖面類型以表聚型為主，分布在研究區(qū)的近海段，隨著人類耕作熟化影響時間的加長和距海距離的增加，土壤鹽分剖面由表聚型向振蕩平均型和底聚型發(fā)展。合理的土地耕作在一定程度上可以減輕土壤鹽漬化，在后期種植管理中應參照研究區(qū)土壤鹽分等值線圖，采取恰當?shù)霓r(nóng)田管理方式，降低土壤鹽漬化，防止次生鹽漬化的發(fā)生。

關鍵詞：黃河三角洲；棉田；土壤含鹽量；時空分布

中圖分類號：S156.4文獻標識號：A文章編號：1001-4942（2017）05-0068-06

Study on Spatial and Temporal Distribution of

Soil Salinity in Cotton Field in the Yellow River Delta

Liu Yang， Jiang Defeng， Zou Xiaoxia

（School of Agriculture and Plant Protection， Qingdao Agricultural University， Qingdao 266109， China）

AbstractThe Yellow River delta is an important reserve cultivated land resources in China. To make clear the spatial variability of soil salinity is the basis of exploiting land resources in this area. In this study， the profile and spatial distribution characteristics of soil salinity in 0～100 cm soil layer in the typical land use type of cotton field were studied with sampling from inland to the coast along the Yellow River. The results found that the soil salinity along the Yellow River from inland to coastal areas was gradually increasing， which was closely related to the soil forming time， the distance from the sea and human activities. The study area was mainly light saline soil， and under the influence of artificial cultivation management， the spatial variability of 0～60 cm soil layer tended to be consistent， which was smaller than that of 60～100 cm soil layer. Surface accumulation was the main soil salinity profile type in the study area and mainly distributed in the coastal area of the study area. Along with the increase of human cultivation time and the distance from the sea， the soil salinity profile type was gradually transformed into the oscillating average type and bottom accumulation type. It showed that reasonable land management could reduce the soil salinization to some extent. In the late planting management， we should refer to the isogram of soil salinity in the study area and take appropriate farmland management ways to reduce the soil salinization and prevent the occurrence of secondary salinization.

KeywordsThe Yellow River delta； Cotton field； Soil salt content； Spatial and temporal distribution

隨著人口增長和經(jīng)濟發(fā)展，我國土地資源短缺與人民生活質(zhì)量需求逐漸提高的矛盾日益尖銳[1]，開發(fā)后備耕地資源勢在必行。作為我國乃至世界造陸速度最快的河口三角洲，黃河三角洲為我國提供了廣闊的后備耕地資源[2]。但該地區(qū)土壤鹽漬化突出，嚴重影響了農(nóng)、林、牧業(yè)的可持續(xù)發(fā)展和生態(tài)安全[3]。因而充分了解黃河三角洲地區(qū)鹽分的空間分布格局對合理開發(fā)該地區(qū)后備土地資源有重大意義。

適當?shù)霓r(nóng)藝措施和生物措施對土壤鹽漬化有一定的改良效果[4]，土壤鹽分的空間分布受植被類型的影響存在差異[5]。棉花作為耐鹽堿性較強的作物，在黃河三角洲地區(qū)有著廣泛的分布[6]。前人對棉田植被下土壤鹽分的空間變異研究主要集中在試驗田[7]或者小區(qū)域尺度[8]，跨度較小，難以體現(xiàn)自然沖擊造陸和人為耕作熟化影響下土壤鹽分的空間變異特性。而沿黃河河道自內(nèi)陸至海岸的空間序列，既體現(xiàn)了黃河自然沖積造陸的時間序列，也體現(xiàn)了人為耕作熟化的影響過程，體現(xiàn)了自然因素和人為因素的協(xié)同影響。

為了揭示黃河三角洲地區(qū)土壤鹽分在自然和人為因素協(xié)同影響下的空間特性，本試驗選取該地區(qū)典型土地利用類型棉田作為研究對象，沿黃河入海方向自內(nèi)陸至沿海布設采樣帶，揭示棉田種植方式下土壤鹽分在剖面及空間上的變異特性，為黃河三角洲地區(qū)鹽堿地合理改良、土地資源

的開發(fā)利用以及科學管理提供參考。

1材料與方法

1.1試驗地概況

采樣帶位于東營市利津縣北宋鎮(zhèn)至河口區(qū)孤島鎮(zhèn)（圖1），沿黃河河道自內(nèi)陸至入海口，介于東經(jīng)118°10′50″～119°09′02″，北緯37°26′37″～37°49′36″之間。試驗區(qū)地處北溫帶大陸性季風型氣候帶，年平均降水量在530～630 mm之間，年均蒸降比約為3.22，年均氣溫11.7～12.6℃[9，10]。地下水埋深1.5～2.4 m，該區(qū)由于沉積物的不斷積累，地面蒸發(fā)強烈，海水浸潤和黃河水側(cè)滲現(xiàn)象較嚴重，導致土壤中可溶性鹽類含量較高，加之淡水資源相對短缺，適合生長的植被類型多為抗鹽堿性植物。經(jīng)過人類長期的土地改良，由堿蓬、蘆葦、檉柳、茅草等植物群落覆蓋下的大批鹽堿地轉(zhuǎn)變成農(nóng)田，尤以棉花居多。采樣帶地塊耕作歷史包括1855年黃河改道前陸地和改道后新沖積土地。

1.2采樣設計

采樣時間為2015年10月，正值棉花成熟期以及穩(wěn)定積鹽期。選擇典型棉田，沿黃河入海方向由遠及近隨機選取11個主采樣區(qū)（圖1），每個主采樣區(qū)內(nèi)沿黃河入海方向“S”形采集5個樣本，共采集55個樣本，每個樣本以20 cm為一個層次，進行0～100 cm土層的分層取樣。

1.3測定方法

采集的土樣帶回實驗室自然風干，磨碎，過2 mm篩后備用，采用殘渣烘干—質(zhì)量法測定水溶性總鹽含量，測定方法參照魯如坤[11]的方法。

1.4數(shù)據(jù)統(tǒng)計方法

本研究采用SPSS 21.0軟件完成土壤含鹽量的統(tǒng)計特征分析，利用Office 2010完成數(shù)據(jù)的統(tǒng)計計算及繪圖，利用Surfer11軟件繪制土壤鹽分等值線圖。

2結(jié)果與分析

2.1土壤含鹽量的統(tǒng)計特征分析

研究區(qū)各土層土壤的平均含鹽量在1.61 ～2.15 g/kg范圍內(nèi)（表1），根據(jù)土壤鹽化分級標準[12]，試驗區(qū)內(nèi)0～100 cm的土層均為輕鹽土，在人類長期的耕作影響下，土壤鹽漬化有所改善，但整個采樣區(qū)的鹽分整體仍呈現(xiàn)表聚特征。

變異系數(shù)（CV）表明了樣點的離散程度，CV>100%表現(xiàn)為強變異性；10%≤CV≤100%表現(xiàn)為中等變異性；CV<10%表現(xiàn)為弱變異性[13]。分析發(fā)現(xiàn)，各土層含鹽量的變異系數(shù)在51.79%～71.14%之間，變異系數(shù)在垂直方向上有增大趨勢，但均屬于中等變異性（表1）。造成這種結(jié)果的原因可能是由于采樣田均種植棉花，人為耕作管理方式相近，導致0～60 cm土壤條件趨于一致，但隨著土層深度的增加，人為耕作影響減弱，而采樣帶跨度較長，海水浸漬側(cè)滲影響逐漸凸顯，導致60～100 cm土層鹽分的空間變異性大于0～60 cm土層。

.2土壤含鹽量的剖面特征分析

土壤鹽分剖面分布由于人類活動、地形以及氣候等因素的影響，具有一定的特點。基于土壤鹽分含量的統(tǒng)計特征分析，可以將研究區(qū)土壤剖面劃分為底聚型、表聚型和振蕩平均型三類。

表聚型鹽分剖面的主要特征是土壤表層鹽分含量較多，其它層次的鹽分含量少于該層，可認為其目前土壤脫鹽過程不太明顯，土壤鹽分的運移處于上升狀態(tài)。研究區(qū)中，這種類型的鹽分剖面占剖面總數(shù)的45.5%，在地理位置上主要分布在試驗區(qū)東部距海較近的地區(qū)（圖2a）。造成這種分布現(xiàn)象的原因，一方面是采樣區(qū)距離海洋較近，受海水的浸漬側(cè)滲較嚴重；另一方面是該區(qū)域主要為1855年黃河改道后新沖積土地，成土時間較短，人類開墾熟化程度低，加之采樣時間正值穩(wěn)定積鹽期，地下水返鹽現(xiàn)象明顯，地表蒸發(fā)強烈，低層鹽分隨水上移，呈現(xiàn)表聚型特征。

底聚型鹽分剖面的主要特征是表層含鹽量低，而下層含鹽量較高，可認為目前土壤鹽分的運移處于下行狀態(tài)。在整個研究區(qū)此類剖面占樣地總數(shù)的27.3%，在地理位置上主要分布在試驗區(qū)的西部地區(qū)（圖2b）。原因主要該區(qū)域為1855年黃河改道之前的陸地，距離海洋較遠，成土時間較長，人類開墾熟化程度高，土壤改良較明顯，呈現(xiàn)脫鹽特征。但當此類剖面中底層鹽分含量較高時，土壤次生鹽漬化的威脅性較大[14]，影響土地質(zhì)量和作物產(chǎn)量，在后期土地耕作過程中一定要注意合理管理，避免次生鹽漬化的發(fā)生。

振蕩平均型鹽分剖面的主要特征是中間層鹽分或高或低，但表層和底層的土壤含鹽量差異不大，可認為目前處于脫鹽或積鹽過程的中間狀態(tài)，土壤鹽分的運移處于動蕩狀態(tài)。在整個研究區(qū)，此類的剖面占樣地總數(shù)的27.3%，在地理位置上主要分布在試驗區(qū)的中部地區(qū)（圖2c）。原因是，該區(qū)域距海距離適中，受海水的浸漬側(cè)滲的影響介于表聚型和底聚型之間，而受人類耕作影響的程度也是介于表聚型和底聚型之間，土壤處于脫鹽和積鹽的中間狀態(tài)。此區(qū)域的土地若能合理開發(fā)改造，將有助于土體脫鹽，降低土壤鹽漬化。

綜上，采樣帶的土壤鹽分剖面類型很好地體現(xiàn)人為耕作和海水浸漬側(cè)滲的協(xié)同影響，研究區(qū)的近海段主要分布著表聚型土壤剖面，隨著人類耕作熟化影響時間的延長和距海距離的增加，土壤鹽分剖面由表聚型向振蕩平均型、進而向底聚型發(fā)展。

2.3土壤含鹽量的空間特征分析

為了得到在水平方向上的不同土層土壤含鹽量的空間分布狀況，根據(jù)Kriging空間插值方法，

利用Surfer11軟件，繪制等值線圖（圖3）。

研究區(qū)沿黃河入海方向自西向東土壤含鹽量呈遞增趨勢，主要是由于沿采樣帶自西向東，沖積成土時間逐漸縮短，人為耕作改良作用逐漸減弱，加之距海距離逐漸縮短，地下水海水浸漬側(cè)滲的程度逐漸增大，鹽分濃度整體呈增高趨勢。但是研究區(qū)表層土壤鹽分等值線變化相對較緩，表明不同區(qū)域間土壤含鹽量差異較小。

整個研究區(qū)的土壤鹽分以表聚型為主，0～20 cm土層土壤含鹽量高于其他土層，隨著土層的加深土壤含鹽量呈降低趨勢，并且鹽分的空間變異性隨著土壤深度的增加有所增加。主要原因是表層土壤受人為耕作作用的影響要大于地勢、地下水位和海水浸漬側(cè)滲的影響，采樣區(qū)均為棉田，耕作管理方式相對一致，導致土壤鹽分的空間差異性較小，隨著土層深度的加深，人為耕作的影響減弱，地下水位、海水的浸漬側(cè)滲等結(jié)構(gòu)性因素影響增強，土壤鹽分的空間變異性增加。

等值線圖在D-2區(qū)出現(xiàn)一個“魚眼”區(qū)（高鹽分含量區(qū)），經(jīng)實地查證該采樣區(qū)地勢較低洼，受地下水浸漬側(cè)滲影響較大，土壤返鹽嚴重。該采樣區(qū)的土壤鹽分剖面為振蕩平均型，但整個土層剖面的鹽分含量均高于臨近幾個采樣區(qū)，這表明雖然該區(qū)域在整個采樣帶的中部區(qū)域，成陸時間較早，受人為耕作影響時間較長，但由于地勢低洼的原因，整個土層返鹽現(xiàn)象仍較重。所以，對于此類區(qū)域，耕作利用需慎重，僅靠耕作管理很難降低土壤鹽漬化，若盲目墾殖利用難以得到理想的效果。

3討論與結(jié)論

黃河三角洲大部分面積的土壤屬于鹽漬化土壤，且鹽漬化程度相當高[15]，但本研究發(fā)現(xiàn)，整個采樣區(qū)棉田土壤鹽分以輕鹽土為主，表明鹽漬化土壤在人類長期的耕作影響下，土壤脫鹽明顯[16]。水平方向上研究區(qū)土壤鹽分呈中等變異性，但60～100 cm土層鹽分的空間變異性大于0～60 cm土層，變異程度主要與距海遠近以及人類活動等因素有關[17]。造成這種結(jié)果的原因可能是由于采樣田均種植棉花，人為耕作管理方式相近，致使0～60 cm土壤條件趨于一致，但隨著土層深度的增加，人為的耕作影響減弱，而采樣帶跨度較長，海水浸漬側(cè)滲的影響逐漸凸顯[6]。

通常將土壤鹽分剖面劃分為底聚型、表聚型和振蕩平均型，表聚型鹽分剖面和振蕩平均型鹽分剖面都屬于非穩(wěn)定狀態(tài)[18]，在采樣區(qū)中這兩種剖面類型占到72.8%。成土時間長短和土地改良狀況的差異在不同土壤鹽分剖面類型中具有一定的體現(xiàn)[19]。研究區(qū)的近海段主要分布著表聚型土壤剖面，隨著人類耕作熟化影響時間的延長和距海距離的增加，土壤鹽分剖面由表聚型向振蕩平均型和底聚型發(fā)展。李寶富等[20]研究也發(fā)現(xiàn)，未耕作土壤含鹽量高，且表層聚集現(xiàn)象明顯，隨著耕作時間的延長，土壤脫鹽明顯，且各土層土壤含鹽量逐漸趨于穩(wěn)定，這表明合理的耕作在一定程度上能夠降低土壤的鹽漬化程度。

通過土壤鹽分等值線圖可以直觀的觀察到研究區(qū)土壤鹽分在水平和垂直方向上的變化，對于及時調(diào)整耕作、灌溉等農(nóng)田管理措施以及提高作物產(chǎn)量具有十分重要的指導作用[21]。整體上由研究區(qū)的西部到東部呈鹽分遞增趨勢，呈現(xiàn)出條帶狀和斑塊狀格局，且各土層鹽分分布存在著較強的空間相似性，其原因主要是隨著采樣帶自內(nèi)陸至沿海方向，海水浸漬側(cè)滲影響增大，土壤鹽分含量增高[9]。從整個采樣帶的剖面結(jié)構(gòu)來看，隨著土壤深度的增加，土壤含鹽量呈降低趨勢，整個采樣帶土壤鹽分剖面分布以表聚型為主，說明土壤鹽漬化仍是制約當?shù)剞r(nóng)業(yè)生產(chǎn)的重要因子[10]。

綜上，研究區(qū)沿黃河入海方向自西向東土壤鹽分呈上升趨勢，與成土時間、距海遠近和人類活動密切相關。整個采樣區(qū)棉田土壤鹽分以輕鹽土為主，受人為耕作管理影響0～60 cm土壤條件趨于一致，空間變異性小于60～100 cm土層；土壤鹽分剖面類型以表聚型為主，分布在研究區(qū)的近海段，隨著人類耕作熟化影響時間的加長和距海距離的增加，土壤鹽分剖面由表聚型向振蕩平均型和底聚型發(fā)展，很好地體現(xiàn)了人為耕作和海水浸漬側(cè)滲的協(xié)同影響。鹽漬化土壤嚴重影響土地質(zhì)量和作物產(chǎn)量，合理的土地耕作在一定程度上可以減低土壤的鹽漬化，應在以后的種植管理中參照研究區(qū)土壤鹽分的等值線圖，采取恰當?shù)霓r(nóng)田耕作方式，降低土壤鹽漬化，防止次生鹽漬化的發(fā)生。

參考文獻：

[1]蔡為民. 黃河三角洲土地利用變化及其可持續(xù)性評價[D].北京：中國農(nóng)業(yè)大學，2004.

[2]Yu J， Wang Y， Li Y， et al. Soil organic carbon storage changes in coastal wetlands of the modern Yellow River Delta from 2000 to 2009[J]. Biogeosciences， 2012， 9（6）：2325-2331.

[3]繆雄誼. 黃河三角洲濱海濕地表層土壤礦物與地球化學特征研究[D].青島：中國海洋大學， 2014.

[4]丁晨曦， 李永強， 董智，等. 不同土地利用方式對黃河三角洲鹽堿地土壤理化性質(zhì)的影響[J].中國水土保持科學， 2013， 11（2）：84-89.

[5]宋玉民， 張建鋒， 邢尚軍，等. 黃河三角洲重鹽堿地植被特征與植被恢復技術[J].東北林業(yè)大學學報， 2003， 31（6）：87-89.

[6]孫運朋， 陳小兵， 張振華，等. 濱海棉田土壤鹽分時空分布特征研究[J].土壤學報， 2013， 50（5）：891-899.

[7]吳向東， 陳小兵， 呂真真，等. 黃河三角洲農(nóng)田灌溉前后土壤鹽分時空變異特征研究[J].干旱地區(qū)農(nóng)業(yè)研究，2013， 31（3）：194-199.

[8]辛承松， 董合忠， 辛梅平. 黃河三角洲鹽堿地棉花生產(chǎn)發(fā)展初探[C]//中國棉花學會2007年年會論文匯編，2007：468-470.

[9]姚榮江， 楊勁松， 姜龍，等. 基于聚類分析的土壤鹽漬剖面特征及其空間分布研究[J].土壤學報， 2008， 45（1）：56-65.

[10]翁永玲， 宮鵬. 黃河三角洲鹽漬土鹽分特征研究[J].南京大學學報（自然科學版）， 2006， 42（6）：602-610.

[11]魯如坤.土壤農(nóng)業(yè)化學分析方法[M].北京：中國農(nóng)業(yè)科技出版社，1999：85-96.

[12]吳從穩(wěn)， 陳小兵， 王松濤，等. 黃河三角洲濱海棉田土壤鹽分特征研究[J].干旱地區(qū)農(nóng)業(yè)研究， 2014， 32（5）：149-153.

[13]王遵親. 中國鹽漬土[M].西安： 陜西人民出版社， 1994：130-149.

[14]杜金龍， 靳孟貴， 歐陽正平，等. 焉耆盆地土壤鹽分剖面特征及其與土壤顆粒組成的關系[J].地球科學-中國地質(zhì)大學學報， 2008， 33（1）：131-136.

[15]范曉梅， 劉高煥， 唐志鵬，等. 黃河三角洲土壤鹽漬化影響因素分析[J].水土保持學報， 2010， 24（1）：139-144.

[16]張凌云. 黃河三角洲地區(qū)鹽堿地主要改良措施分析[J].安徽農(nóng)業(yè)科學， 2007， 35（17）：5266，5309.

[17]董芝，劉加珍，陳永金，等. 黃河三角洲典型群落土壤鹽分空間變異性分析[C]//中國自然資源學會，新疆自然資源學會.發(fā)揮資源科技優(yōu)勢 保障西部創(chuàng)新發(fā)展——中國自然資源學會2011年學術年會論文集（下冊），2011.

[18]劉慶生， 劉高煥， 范曉梅. 黃河三角洲土壤鹽分剖面類型時空分布研究[J].山東農(nóng)業(yè)科學， 2010（1）：57-62.

[19]姚榮江， 楊勁松， 姜龍. 黃河下游三角洲鹽漬區(qū)表層土壤積鹽影響因子及其強度分析[J].土壤通報， 2008， 39（5）：1115-1119.

[20]李寶富， 熊黑鋼， 張建兵，等. 不同耕種時間下土壤剖面鹽分動態(tài)變化規(guī)律及其影響因素研究[J]. 土壤學報， 2010， 47（3）：429-438.

[21]董亮， 孫澤強， 王學君，等. 黃河三角洲鹽堿耕地型中低產(chǎn)田概況及改良增產(chǎn)技術探討[J].江西農(nóng)業(yè)學報， 2014，16（2）：58-62.