渝西北土壤有機質空間變異及影響因素分析

鮑麗然，周 皎，李 瑜，賈中民

(重慶市地質礦產勘查開發局川東南地質大隊，重慶 400038)

渝西北土壤有機質空間變異及影響因素分析

鮑麗然，周 皎，李 瑜，賈中民

(重慶市地質礦產勘查開發局川東南地質大隊，重慶 400038)

【目的】深入認識重慶西北部丘陵區土壤有機質空間變異特征及影響因素，為農業生產、土地利用等提供科學依據。【方法】采用地統計學和GIS相結合的方法，進行有機質空間變異和基于土地利用類型修正的空間插值研究。【結果】研究區土壤有機質平均含量為16.26 g/kg，屬中等偏下水平；變異系數為33.8 %，存在中等程度的空間變異性；土壤有機質空間分布呈現中間高，兩端低的趨勢，合川區整體偏低。【結論】土壤有機質具有中等空間自相關性，空間變異性受結構因素(地形坡度、土壤類型)和隨機因素(土地利用類型、耕作措施等)影響，土地利用類型是主控因素。

地統計學；土壤有機質；空間變異；影響因素

【研究意義】有機質是土壤的重要組成部分，包括各種動植物的殘體、微生物體及其分解和合成的各種有機質。土壤有機質在耕地養分管理、農業發展等方面具有重要意義，是土壤學、農學和地球化學研究的重點問題之一[1-2]。查明土壤有機質的空間變異特征及其影響因素，有助于制定合理的作物種植規劃及農田配方施肥，促進耕地的高效利用，保障農業可持續發展[3]。【前人研究進展】土壤有機質空間變異的研究始于20世紀70年代[4]，隨著GIS技術的發展，80年代后多位學者相繼運用地統計學結合GIS方法，探討了農田土壤有機質空間變異特征[5-7]，如Moorman等[6]對包括土壤有機質在內的多種養分進行了空間變異研究，Mishra等[7]曾針對紅壤地區的有機質和土壤酸堿度展開空間變異特征分析。近年，國內學者對土壤有機質空間變異性在農田尺度、區域尺度上展開了大量研究，研究區涉及平原區、干旱荒漠區、丘陵區、典型喀斯特區、高原流域、山區小流域、城市近郊等[8-13]，主要探討了自然因素(地形地貌、氣候和土壤類型等)和人為因素(土地利用和耕作方式等)對有機質空間變異的影響。【本研究切入點】2003年以來，重慶市多目標區域地球化學調查獲得了大量高精度土壤有機質數據，為研究其空間變異特征提供了基礎。渝西北地區地勢平緩、土壤肥沃，耕地廣泛分布，農業是該區最重要的經濟支柱，因此土壤有機質的分布特征至關重要。【擬解決的關鍵問題】基于GIS操作平臺，研究了重慶西部丘陵區土壤有機質空間變異特征及影響因素，為全面、準確掌握該區土壤有機質空間變異規律，更好的進行農業生產活動、合理規劃土地利用等提供科學依據。

1 材料與方法

1.1 研究區概況

研究區位于重慶市西北部，地理坐標為東經105°17′ ～ 106°39′，北緯29°15′ ～ 30°26′，總面積約7130 km2(圖1)。地勢起伏變化不大，海拔一般為200 ～ 500 m，廣大平緩淺丘遍布全區。研究區出露地層主要為侏羅系沙溪廟組、遂寧組和蓬萊鎮組，其次為三疊系須家河組，第四系零星分布于山間河谷地區。礦產有煤、天青石、石灰巖、鹽和天然氣等。土地利用類型以耕地為主，占土地總面積的60.1 %。土壤類型主要為廣大丘陵區的紫色土和水稻土，其次為黃壤、石灰土。該區屬亞熱帶季風氣候，氣候溫和，雨量充沛，年均降雨量1022 mm。區內分布嘉陵江水系和沱江水系兩大長江一級支流，嘉陵江水系包括渠江、瓊江和涪江等，沱江水系包括瀨溪河及其支流。糧食作物主要為水稻、玉米、油菜、大豆、薯類，建有大型果蔬基地，是重慶市重要的糧油果蔬生產基地。

1.2 樣品采集與測試

按照1件/km2的密度采集0 ～ 20 cm表層土壤，樣點布設以代表性為主要原則，同時兼顧均勻性和合理性，最大限度控制調查面積。采樣時除去地表雜物，樣品加工處理后重量不低于800 g。將大格內(4 km2)樣品組合為一個分析樣，全區共采集表層土壤樣品7135件，組合分析樣1783件。

樣品分析測試按照DD2005-01《多目標區域地球化學調查規范》(1∶250000)[14]的相關要求，由安徽省地質實驗研究所完成。實驗室人員采用重鉻酸鉀容量法分析檢測，檢出限0.05 %，報出率100 %，合格率100 %。

1.3 數據處理

運用SPSS19.0軟件對土壤有機質進行經典描述性統計分析；依據峰度、偏度值和非參數檢驗的單樣本K-S值判別統計數據是否符合正態分布；利用單因素方差分析F檢驗和回歸分析，分別判斷有機質含量差異顯著性和不同因素對土壤有機質的影響程度。

借助地統計學軟件GS+9.0進行半變異函數計算和模型擬合，結合專業知識分析土壤有機質的空間變異性及其影響因素。

圖1 地理位置Fig.1 Geographical position of the study area

指標Index樣本Samples含量范圍Contentrange中位數Median均值Mean標準差Std.dev.變異系數c.v.偏度Skewness峰度KurtosisK-S檢驗K-Stest土壤有機質(g/kg)17836.36～63.7115.5216.265.5033.833.2528.810.0對數轉換0.3941.540.165

基于Arcgis10.0的Geostatistic analysis模塊，采用Kriging法對土壤有機質殘差進行空間插值，運用柵格計算將土地利用圖層和殘差預測圖層相加，得到土壤有機質含量空間分布圖。

2 結果與分析

2.1 土壤有機質經典統計分析

從一般的經典統計分析結果(表1)看，土壤有機質含量介于6.36 ～ 63.71 g/kg，平均值為16.26 g/kg，根據《土地質量地球化學評價規范》屬較缺乏水平[15]；變異系數33.8 %，揭示了有機質含量的離散程度，屬中等變異強度[16]。峰度值為28.8，較陡峭，并且K-S檢驗中P值為0，明顯不符合正態分布。經對數轉換后，偏度和峰度值分別為0.394和1.54，再經K-S檢驗表明符合正態分布(Pgt;0.05)，滿足地統計學分析的基本要求。

2.2 土壤有機質的空間變異結構特征

運用GS+9.0軟件對土壤有機質含量進行各向同性變異函數擬合(表2)，按照決定系數(R2)最大、殘差(RSS)最小的原則選擇最優半方差函數擬合模型，最佳理論模型為指數模型。

空間變異主要包括隨機性和結構性兩部分變異。塊金值C0表示區域化變量內部隨機性的可能程度，由測量誤差和最小取樣間距內土壤性質的人為隨機因子(施肥、耕作措施、土地利用方式等)引起。結構方差C是由自然因素(地質背景、地形、土壤類型等)引起土壤性狀的結構變異[17]。基臺值C0+C通常表示系統總變異，包括隨機變異和結構變異。塊金系數，即塊金值和基臺值之比C0/C0+C，表示空間異質性程度，比值高表明隨機因素引起的空間變異程度較大；相反，則結構性因素引起的空間變異程度較大[18]。C0/C0+Clt;25 %，說明系統具有強烈的空間相關性；25 %～ 75 %，說明系統具有中等空間相關性；C0/C0+Cgt;75 %，說明系統空間相關性弱[19]。變程反映在一定觀測尺度下，變量的空間自相關距離，只有變程以內的空間變量才具有空間自相關性。

研究區土壤有機質的塊金系數C0/C0+C為46.2 %，表現出中等空間相關性，說明結構性因素和人為隨機因素土壤有機質的空間變異影響作用相當。變程為184 km，說明土壤有機質在此空間范圍內存在空間自相關性。

2.3 土壤有機質空間變異影響因素分析

2.3.1 土地利用類型 不同土地利用類型下，人類對土地的使用、保護和改造活動不同，土壤有機質含量隨之產生了較大差別[20]。研究區土地利用主要為旱地和水田，其次為林地。林地土壤有機質平均含量最高(表3)，為23.82 g/kg，水田為17.15 g/kg，旱地只有14.72 g/kg，相當于水田的85 %；旱地的變異系數為0.35，是水田的1.24倍。研究區水田更易于土壤有機質的累積，且含量較穩定，與許泉[21]等發現的不同土地利用類型下農田土壤有機碳密度區域差異規律相似。方差分析(表4)顯示，組間均方大于組內均方，表明組間差異大于組內差異，F檢驗結果(Plt;0.05)說明不同土地利用類型土壤有機質含量具有顯著性差異。

表2 土壤有機質空間變異的理論模型和相應參數

表3 不同土地利用類型土壤有機質含量

表4 不同土地利用類型土壤有機質含量方差分析

注：F檢驗Plt;0.05。

2.3.2 地形(海拔、坡度) 海拔、坡度等地形因子也會影響土壤有機質含量。高海拔區，人類活動較少，易于積累土壤有機質。坡度較大地區，強烈的水土流失會攜帶走大量的土壤有機質。

將研究區1073個樣點有機質含量和海拔高度進行相關分析，得出相關系數為0.170，兩者關系極弱，認為不相關[22]。進一步統計了不同海拔高度上各等級有機質含量百分比。由圖2可知，不同海拔高程帶上土壤有機質均以10 ～ 20 g/kg為主，占所在高程帶的80 %左右，其次為20 ～ 30 g/kg。隨著海拔高度的變化，土壤有機質沒有表現出較大的差異性。本區屬低緩丘陵區，不存在高山或洼地，地勢起伏、氣候差異較小，土地利用、人類耕作活動少受海拔影響。因此，海拔高度對土壤有機質的空間變異影響不明顯。

從地形的坡度因素考慮，根據各采樣點坡度信息，得到不同坡度范圍土壤有機質均值(表5)，土壤有機質的空間分布和坡度密切相關。坡度小于5°時，有機質含量最高，為19.77 g/kg；大于35°時，有機質含量最低，為12.23 g/kg；隨著坡度增大，土壤有機質含量呈逐漸降低趨勢。這與李婷等[23]對沱江流域中游土壤有機質含量與坡度關系研究結果相符，不同坡度土地利用類型和土地管理水平不同，低緩坡度土地利用主要為水田，比較高陡的坡度土地利用為旱地和林地。另一方面研究區屬紫色丘陵區，坡度大的區域水土流失較為嚴重，大量土壤有機質被攜帶流出。

圖2 不同海拔土壤有機質分布特征Fig.2 Distributions of soil organic matter in different altitudes

2.3.3 土壤類型 土壤類型也是有機質含量的重要影響因素之一。研究區土壤類型以紫色土和水稻土為主，分別占總面積的70 %和21.6 %，黃壤和石灰土分布較少。紫色土主要分布于北部潼南、合川、大足地區，水稻土于銅梁、榮昌地區分布較集中。統計分析(表6)表明，土壤有機質平均含量表現為石灰土gt;黃壤gt;水稻土gt;紫色土，最高值石灰土是最低值紫色土的1.79倍，水稻土和黃壤差異較小。這與成都平原不同類型表層土壤中有機質分布規律相同[24]，均為石灰土、黃壤中較高，紫色土最低。由表7可知，不同土壤類型間，組間均方大于組內均方，組間差異較大，F檢驗結果(Plt;0.05)說明不同土壤類型中有機質含量具有顯著性差異。

表5 不同坡度土壤有機質含量

表6 不同土類有機質含量組成

表7 不同土類土壤有機質含量方差分析

注：F檢驗Plt;0.05。

表8 不同影響因子與土壤有機質相關性

2.4 土壤有機質的空間分布格局

土壤有機質空間分布格局可借助空間插值模型實現從離散的點狀信息向連續面狀信息的轉換，用有限的樣點預測整個研究區域的土壤屬性空間分布特征。在空間尺度較小、地勢平坦且土地利用類型單一的區域，通常用普通克里格插值進行預測[25]，而研究區空間尺度較大，地形條件復雜、土地利用類型多樣，應利用土壤有機質相關因子與克里格插值相結合的方法來展示土壤有機質的空間分布。目前，結合相關因子的空間插值的研究尚不多見。劉靜等[26]在研究梁平地區土壤有機質空間分布時，結合地形因素改進了空間插值方法。顧成軍、文雯等[27-29]探討了基于土地利用類型修正的殘差法空間插值問題。

本文首先通過回歸分析定量判斷了土地利用、土壤類型、地形坡度等因素對土壤有機質空間變異的影響程度(表8)。研究區土壤有機質含量與不同影響因素的相關系數大小依次為土地利用類型gt;土壤類型gt;地形(坡度)，土地利用類型是該區土壤有機質分布最為重要的影響因素。因此，為提高土壤有機質空間分布預測精度，本文選擇結合土地利用類型的克里格插值法展開研究。首先將樣本土壤有機質含量去除土地利用類型均值，得到殘差，并將其作為一個新的區域變量行普通克里格插值。然后，將插值過后的殘差值圖層按土地利用類型的不同，分別加上其地類均值，即得出土壤有機質密度最終預測值，由此生成該區土壤有機質含量空間分布圖(圖3)。

研究區土壤有機質含量插值結果范圍為1.96

圖3 研究區土壤有機質空間分布Fig.3 Spatial distribution of soil organic matter in the study area

～ 71.21 g/kg，空間分布上整體呈中間高，南北兩端低的趨勢。土壤有機質高值區主要分布于大足區東部、銅梁區西部、潼南區南部和榮昌區最南端等地區。大足區境內高值區占總面的60 %左右。這些高值區土地利用類型以水田為主，并有林地分布。合川區土壤有機質整體偏低，榮昌區北部含量最低。合川區和潼南區北部分布嘉陵江及其支流渠江、涪江，水網交錯，水土流失相對嚴重，這可能是土壤有機質低的主要因素。

3 討 論

研究區土壤有機質的空間變異主要受土地利用、土壤類型的影響，不同土地利用或土壤類型的有機質含量具有顯著差異。

從土地利用類型看，林地區植被遍布，源源不斷的樹木殘枝落葉進入土壤，有機質來源充足，因此有機質含量最豐富。耕地中水田土壤有機質高于旱地，主要是由于水田區海拔較低，地形平緩和土地平整，這些都屬于有機質聚集的有利的條件，并且水田一般為滯水狀態，有機質的分解速率較低，容易積累[23]。另一方面，相比旱地而言，水田作物秸稈殘存量多，大部分被微生物轉化為腐殖質，構成了土壤有機質的重要來源。

從土壤類型看，該區石灰土、紫色土、水稻土和成土母巖不同，物質來源不同，加之土壤自身的理化性質不同，導致了不同土類中有機質含量的差異。石灰土的成土母巖為三疊系碳酸鹽巖(灰巖和白云巖)，這類巖石中生物化石豐富，攜帶了較多的有機質成分[30]；同時石灰土本身的物質成分也有利于有機質的貯存[31]。紫色土、水稻土和黃壤由侏羅系或三疊系的紫色砂巖、泥巖及粉砂巖風化而成，該類巖石形成于陸相環境，有機質成分較少。紫色土呈中性或微堿性，易風化、墾殖率高、水土流失嚴重，不利于有機質的保存[32]。

4 結 論

(1)重慶西北部丘陵區土壤有機質平均含量為16.26 g/kg，屬中等偏下水平。有機質的空間變異系數為33.8 %，屬中等變異程度。空間結構采用指數型模型擬合較好，有機質的塊金系數C0/C0+C為46.2 %，表明土壤有機質具有中等空間相關格局，其空間變異性受結構性因素和隨機因素影響程度相當。

(2)研究區土壤有機質分布受土地利用類型、地形坡度、土壤類型等因素影響。不同土地利用類型下，林地土壤有機質含量最高，水田次之。不同坡度下，低緩坡度中有機質含量較高；海拔高度對有機質含量影響不明顯。不同土壤類型，石灰土有機質含量最高，黃壤和水稻土次之，紫色土含量最低。土地利用類型是有機質空間分布最重要的影響因素。

(3)基于土地利用類型修正的殘差法空間插值預測顯示，研究區土壤有機質含量空間分布呈中間高，南北兩端低的趨勢。合川區水土流失較嚴重，土壤有機質含量整體偏低。

[1]黃昌勇. 土壤學[M]. 北京: 中國農業出版社, 2000: 32.

[2]韓 丹, 程先富, 謝金紅, 等. 大別山區江子河流域土壤有機質的空間變異及其影響因素[J]. 土壤學報, 2012, 49(2): 403-407.

[3]苑小勇, 黃元仿, 高如泰, 等. 北京市平谷區農用地土壤有機質空間變異特征[J]. 農業工程學報, 2008, 24(2): 70-75.

[4]WARRICK A W. Analytical solutions to the one-dimensional linarized moisture flow equation for arbitrary input[J]. Soil Science, 1975, 120: 79-84.

[5]CANN M D, HUMMEL J W, BROUER B H. Spatial analysis of soil fertility for site specific crop management[J]. Soil Sci Soc Am J, 1994, 58: 1240-1258

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

[7]MISHRA T K, BANERJEE S K. Spatial variability of soil pH and organic matter under Shorea robusta in lateritic region[J]. Indian Journal Forest, 1995, 18(2): 144-152.

[8]趙明松, 張甘霖, 王德彩, 等. 徐淮黃泛平原土壤有機質空間變異特征及主控因素分析[J]. 土壤學報, 2013, 50(1):1-9.

[9]楊 東, 劉 強. 河西地區土壤全氮及有機質空間變異特征分析-以張掖市甘州區為例[J]. 干旱地區農業研究, 2010, 28(4): 183-186.

[10]楊 葳. 紫色丘陵區土壤養分空間變異及肥力評價研究——以重慶市銅梁縣為例[D]. 重慶: 西南大學, 2012.

[11]王建鋒, 謝世友. 喀斯特石漠化地區土壤養分空問異質性研究構想[J]. 湖北農業科學, 2010, 49(7): 1605-1608.

[12]邱 揚, 傅伯杰, 王 軍, 等. 黃土高原流域土壤養分的時空變異及其影響因子[J].自然科學進展, 2004, 14: 294-299.

[13]宋 莎, 李廷軒, 王永東, 等. 縣域農田土壤有機質空間變異及其影響因素分析[J]. 土壤, 2011, 43(1): 44-49.

[14]中國地質調查局. DD2005-01《多目標區域地球化學調查規范》(1∶250000)[S]. 2005: 10.

[15]中華人民共和國國土資源部. DZ/T 0295-2016,土地質量地球化學評價規范[S]. 2016: 41.

[16]陳延良, 雷國平. 松嫩平原南部土壤養分空間變異規律分析及其等級評價[J]. 揚州大學學報(農業與生命科學版), 2010, 31(2): 43-47.

[17]蔣勇軍, 袁道先, 謝世友, 等. 典型巖溶流域土壤有機質空間變異——以云南小江流域為例[J]. 生態學報, 2007, 27(5):2041-2045.

[18]劉愛利, 王培法, 丁 園. 地統計學概論[M]. 北京: 科學出版社, 2012: 60-67.

[19]孔 毅, 劉國順, 臘貴曉, 等. 小尺度下基于GIS和地統計學的煙田土壤養分空間變異[J]. 江蘇農業科學, 2012, 40(4): 317-320.

[20]SCHUIP C J E, VELDKAMP A. Long-term landscape-land use interactions as explaining factor for soil orgnic matter variability in Dutch agricultural landspace[J]. Geodema, 2008, 146: 457-465.

[21]許 泉, 芮雯奕, 何 航, 等. 不同土地利用方式下中國農田土壤有機碳密度特征及區域差異[J]. 中國農業科學, 2006, 39(12): 2505-2510.

[22]孫逸敏. 利用 SPSS軟件分析變量間的相關性[J]. 新疆教育學院學報, 2007, 23(2):120-122.

[23]李 婷, 張世熔, 劉 潯. 沱江流域中游土壤有機質的空間變異特點及其影響因素[J]. 土壤學報, 2011,48(4): 864-867.

[24]奚小環, 楊忠芳, 廖啟林, 等. 中國典型地區土壤碳儲量研究[J]. 第四紀研究, 2010, 30(3):573-583.

[25]申廣榮, 葛曉燁, 黃秀梅. 上海崇明島表層土壤有機碳密度的空間分布特征及碳儲量估算[J]. 上海交通大學學報(農業科學版), 2011, 29(6): 61-66.

[26]劉 靜, 蔡國學, 劉洪斌. 西南丘陵地區土壤有機質含量的空間插值法研究[J]. 西南大學學報(自然科學版), 2008, 30(3): 107-111.

[27]文 雯, 周寶同, 汪亞峰, 等. 基于輔助環境變量的土壤有機碳空間插值——以黃土丘陵區小流域為例[J]. 生態學報, 2013, 33(19): 6390-6395.

[28]顧成軍, 史學正, 于東升. 結合土地利用與克里格插值的區域土壤有機碳空間表征方法[J]. 水土保持研究, 2014, 21(2):40-42.

[29]范勝龍, 黃 炎, 林金石. 表征土壤有機碳區域分布的優化空間插值模型研究[J]. 水土保持研究, 2011, 18(6): 1-5.

[30]杜遠生，童金南. 古生物地史學概論[M]. 北京: 中國地質大學出版社, 1998: 150-194.

[31]曹建華, 袁道先, 潘根興. 巖溶生態系統中的土壤[J]. 地球科學進展, 2003, 18(1): 37-44.

[32]倪九派, 袁道先, 謝德體, 等. 重慶巖溶區土壤有機碳庫的估算及其空間分布特征[J]. 生態學報, 2009, 29(11):6292-6300.

(責任編輯 李 潔)

SpatialVariabilityofSoilOrganicMatterandItsInfluenceFactorsofHillyAreainNorthwestChongqing

BAO Li-ran , ZHOU Jiao, LI Yu , JIA Zhong-min

(Southeast Sichuan Geological Group, Chongqing Bureau of Geology and Minerals Exploration, Chongqing 400038, China)

【Objective】The spatial variability of soil organic matter and its influence factors of hilly area in Northwest Chongqing ws discussed for providing scientific basis on agricultural production and land uses. 【Method】The geostatistics and GIS methods were used to study, spatial variability of soil organic matter and spatial interpolation based on land uses amendment.【Result】The soil organic matter average content was 16.26 g/kg, in the lower middle level and the variation coefficient was 33.8 %, in the middle level. The spatial distribution of soil organic matter was higher in middle, lower on both ends, and the whole Hechuan district was lower. 【Conclusion】Geostatistical analysis suggested that soil organic matter had medium spatial correlation; The spatial variability were resulted from structure factors(topography slope and soil types) and random factors(land uses, farming measures, ect.), and land uses was the dominated factor.

Geostatistics;Soil organic matter; Spatial variability;Influence factors

P596

A

1001-4829(2017)11-2541-07

10.16213/j.cnki.scjas.2017.11.026

2016-08-26

全國土壤現狀調查及污染防治 ([2008]GZTR01-11)；重慶市國土房管科技計劃項目(K J-2015043)

鮑麗然(1983-)，女，河北石家莊人，工程師，碩士，從事環境地球化學研究，Email:305946963@qq.com。