基于地形分類的黃土丘陵區(qū)土壤有機質(zhì)插值方法

張婧 畢如田 丁皓希 許月鳳 文偉杰

摘要：黃土丘陵山區(qū)地形復(fù)雜，耕地細(xì)碎，采樣點密度不足，因此研究適宜該區(qū)域的空間插值方法具有重要意義。以山西省偏關(guān)縣為研究區(qū)，基于數(shù)字高程模型（digital elevation model，簡稱DEM）數(shù)據(jù)生成地形因子，依據(jù)地形因子將研究區(qū)劃分為不同地形單元，利用Kriging插值法對各地形單元進行分層插值，通過疊加不同地形單元的插值結(jié)果，得到全區(qū)域的土壤有機質(zhì)含量空間分布，并利用均方根誤差和決定系數(shù)評價空間插值方法的預(yù)測精度。結(jié)果表明，采用劃分地形單元分層克里金（Kriging）插值法得到的預(yù)測值與實測值的決定系數(shù)為0.375 3，明顯高于全局插值法;均方根誤差精度高于全局插值法，預(yù)測值無偏性好。因此，基于地形單元的空間插值方法可以更精確有效地獲取土壤有機質(zhì)空間分布特征，為復(fù)雜山區(qū)低密度采樣下的耕地質(zhì)量調(diào)查與評價提供技術(shù)參考。

關(guān)鍵詞：丘陵山區(qū);地形單元;土壤有機質(zhì);Kriging插值;DEM數(shù)據(jù)

中圖分類號： S159文獻標(biāo)志碼： A

文章編號：1002-1302（2021）08-0182-07

收稿日期：2020-11-25

基金項目：國土資源部公益性行業(yè)項目（編號：201411007） 。

作者簡介：張 婧（1995—），女，山西呂梁人，碩士研究生，研究方向為土地利用與信息技術(shù)。E-mail：zj972775@163.com。

通信作者：畢如田，博士，教授，博士生導(dǎo)師，主要從事資源環(huán)境信息技術(shù)、土地復(fù)墾與規(guī)劃研究。E-mail：brt@sxau.edu.cn。

有機質(zhì)是土壤的重要成分和肥力指標(biāo)，精確掌握耕地土壤有機質(zhì)的空間變異規(guī)律，有利于耕地精準(zhǔn)施肥和耕地養(yǎng)分的科學(xué)管理[1-3]。土壤在時間和空間是一個連續(xù)的變異體，在空間上具有相關(guān)度很高的異質(zhì)性，同時具有尺度依賴性[4]，而作為土壤的重要養(yǎng)分之一，有機質(zhì)也具有同樣的時空特征[5-6]。

大量研究表明，土壤有機質(zhì)的空間變異受多種因素的影響，例如地形因子、成土母質(zhì)、耕作措施、土地利用方式等，地統(tǒng)計學(xué)的研究中，常忽略這些因素對有機質(zhì)質(zhì)量比空間再分配的協(xié)同作用。付曉婷提出了基于地貌類型分類，研究土壤養(yǎng)分空間變異特征方法，但其劃分標(biāo)準(zhǔn)不是很明確，對于一個縣域來說尺度比較大，且指標(biāo)單一，這是傳統(tǒng)地形分類現(xiàn)存的問題[7]。對于地貌類型劃分，韋金麗等以數(shù)字高程模型（digital elevation model，簡稱DEM）為基礎(chǔ)數(shù)據(jù)，提取坡度，設(shè)計分級標(biāo)準(zhǔn)然后分類組合進行大致的地形因子歸類，此方法劃分結(jié)果比較粗糙，但為劃分地形單元提供了一定的參考[8]。耕地土壤有機質(zhì)空間變異規(guī)律主要受地形顯著影響，在田塊尺度上，Yanai等選擇哈薩克斯坦北部半干旱區(qū)探討影響耕地有機質(zhì)變異的因子，結(jié)果表明地形是主要因素[9];在空間多尺度上，朱洪芬等選擇太原盆地為研究區(qū)，分析不同地形部位土壤有機質(zhì)變異特征與影響因子的空間多尺度關(guān)系，并且得到了較好的預(yù)測結(jié)果[10]。趙業(yè)婷等運用描述性統(tǒng)計結(jié)合地統(tǒng)計理論，研究了耕地養(yǎng)分空間特征規(guī)律以及相互之間的影響關(guān)系[11];耿廣坡等利用克里金（Kriging）插值結(jié)合遙感的方法探討了各種地形因子與耕層土壤養(yǎng)分含量空間分布規(guī)律的關(guān)系[12]。對土壤屬性空間預(yù)測最常用的方法是Kriging插值法[13-14]，但是普通Kriging插值法比較適用于平坦、地形單一的區(qū)域，對于地形復(fù)雜地區(qū)插值效果一般。根據(jù)地理學(xué)第一定律（土壤屬性信息比較類似是因為具有相同空間位置屬性的因子特征越相近造成的，土壤屬性差異比較大是因為空間位置不同）可知，相同的地形單元里，在其地形條件相似的情況下，空間距離上呈現(xiàn)出較細(xì)微的差異，因此土壤有機質(zhì)空間變異也服從這一規(guī)律;而對于不同的地形單元，由于地形條件不同，造成空間距離很大差異，因此各種地形因子對土壤有機質(zhì)空間分布規(guī)律有很大的影響。田瑞云等基于多種地形因素，對DEM數(shù)據(jù)進行分析，分類組合劃分出典型地形單元[15]。黃魏等選取地形因子研究土壤有機質(zhì)空間變異規(guī)律特征，提出基于地形單元的空間插值方法，獲取較高精度的有機質(zhì)空間分布以及在地統(tǒng)計中地形的協(xié)同影響，為精細(xì)研究土壤屬性局部空間規(guī)律特征提供了具有實踐性的參考價值[16]。

針對黃土高原丘陵山區(qū)，由于氣候條件與耕作方式基本相同，因此，地形是引起土壤有機質(zhì)空間變異的重要原因。本研究嘗試依據(jù)地形因子特征屬性將研究區(qū)劃分為不同地形單元，單獨分析各地形單元下的土壤有機質(zhì)屬性的空間分布特征，以期為耕地土壤質(zhì)量調(diào)查、土地資源管理以及精準(zhǔn)農(nóng)業(yè)的實施提供參考依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

選擇山西省忻州市偏關(guān)縣為研究區(qū)，地處晉西北的晉蒙交界處（39°12′～39°40′N，111°22′～112°01′E），北依長城，西鄰黃河，與內(nèi)蒙古準(zhǔn)格爾旗隔河相望，屬北溫帶大陸性氣候，年均氣溫為3～8 ℃，年均降水量為425.3 mm。研究區(qū)處于黃河中游黃土丘陵區(qū)，縣域內(nèi)丘陵起伏，溝壑縱橫，地勢東高西低，地形復(fù)雜多樣，海拔高度為877～1 844 m。研究區(qū)耕地面積為37 131.03 hm2，耕地細(xì)碎，分布于全縣各種地形類型，如圖1-a所示。偏關(guān)縣成土母質(zhì)主要為黃土母質(zhì)、人工堆墊、沙質(zhì)黃土，土壤類型主要為風(fēng)沙土、黃綿土和栗褐土。收集3 007個采樣點，約每12.3 hm2 有1 個采樣點，不同地形采樣間距不均勻，如圖1-b所示。

1.2 數(shù)據(jù)來源及處理

本研究所需數(shù)據(jù)包括分辨率為12.5 m的DEM數(shù)據(jù)[來源于中國科學(xué)院計算機網(wǎng)絡(luò)信息中心地理空間數(shù)據(jù)云網(wǎng)站（www.gscloud.cn）]、土壤有機質(zhì)數(shù)據(jù)（來源于山西省偏關(guān)縣2010年測土配方項目采樣）、耕地數(shù)據(jù)（來源于2015年1 ∶10 000偏關(guān)縣土地利用變更調(diào)查數(shù)據(jù)庫以及2017年偏關(guān)縣耕地質(zhì)量等級更新成果）。以上數(shù)據(jù)均經(jīng)空間數(shù)據(jù)格式轉(zhuǎn)換、 校正并投影轉(zhuǎn)化為2000國家大地坐標(biāo)系（CGCS2000）[17]。利用ArcGIS軟件對數(shù)字高程模型進行地形因子提取;采用SPSS軟件對樣點數(shù)據(jù)進行描述性統(tǒng)計分析。

2 研究方法

2.1 地形單元劃分及設(shè)計

地形因子代表的空間特征與環(huán)境效應(yīng)相比有很大的差異，所以在土壤屬性的空間變異規(guī)律中，地形可以作為衡量整個空間特征的主要影響因素[18-19]。本研究根據(jù)研究區(qū)實際情況選取坡度（slope）、坡向（aspect）、地形位置指數(shù)（topographic position index，簡稱TPI）等3個地形指標(biāo)，利用 12.5 m 精度高程數(shù)據(jù)，在ArcGIS軟件中分析得到3個地形指標(biāo)，其中運用反距離加權(quán)修正法對地形位置指數(shù)進行修正[15]。地形因子在各類地形條件下分布特征具有一定的規(guī)律性，依據(jù)這一屬性的特征，可將偏關(guān)縣耕地劃分為不同的地形單元。本研究主要從DEM中提取出各因子，依據(jù)各因子屬性特征重分類組合，運用疊加分析劃分出地形單元。由于受自然條件的影響，研究區(qū)地形復(fù)雜，地形類型多元，依據(jù)研究區(qū)的實際情況，本著真實性、典型性、全面性的原則選取11種地形類型，運用柵格計算器按照圖2中地形因子分類組合形成目標(biāo)單元類型。其中，依據(jù)耕地坡度分級標(biāo)準(zhǔn)將坡度分為0°～6°、6°～15°、15°～25° 3級;根據(jù)光照方向?qū)⑵孪騽澐譃殛柶?、半陽坡、陰坡和半陰?個坡向;根據(jù)地形位置指數(shù)（TPI）結(jié)合坡度、坡向和土壤有機質(zhì)空間屬性，將耕地劃分為坡腳、平地、平頂峰、溝谷、陰坡、半陰坡、陽坡、半陽坡、陡坡、山脊、坡肩等共11個地形單元。根據(jù)表征出來的實際情況，提出地形位置指數(shù)值的分級標(biāo)準(zhǔn)對TPI進行分級，其中SD表示它的標(biāo)準(zhǔn)偏差[20]，劃分指標(biāo)見圖2。最后依據(jù)地形單元結(jié)果分區(qū)劃分采樣點，獲取各種地形下的土壤采樣點，為分地形單元插值作基礎(chǔ)。

2.2 全局空間自相關(guān)

全局空間自相關(guān)通常用來反映地理實體的空間關(guān)系，而耕地土壤有機質(zhì)在空間分布上變異規(guī)律也具有差異性。本研究選擇全局莫蘭指數(shù)（I）來衡量全局空間自相關(guān)，公式如下：

I=∑ni=1∑nj=1Wij（Xi-X）（Xj-X）S2∑ni=1∑nj=1Wij。（1）

式中，S2=∑ni=1（Xi-X）2n為樣本方差;Wij是要素i和j之間的空間權(quán)重;∑ni=1∑nj=1Wij為所有空間權(quán)重之和;n為采樣點數(shù)量;Xi、Xj為空間中i和j要素的屬性值;X表示土壤有機質(zhì)實測值的平均值。

將各類地形單元下的土壤有機質(zhì)作為觀測變量，通過此方法分析偏關(guān)縣不同地形單元內(nèi)的土壤有機質(zhì)空間變異特征。由于各地形單元面積大小不等，各地形單元采樣點數(shù)量不一樣，土壤有機質(zhì)蘊含內(nèi)在方差不穩(wěn)定可能不滿足空間自相關(guān)分析中隱含的數(shù)據(jù)穩(wěn)定性假設(shè)，推斷出錯誤的結(jié)果，所以進行全局自相關(guān)分析時，首先對土壤有機質(zhì)數(shù)據(jù)進行標(biāo)準(zhǔn)化[21]。

2.3 聚類分析

利用統(tǒng)計軟件SPSS采用歐氏距離和組平均法對11個地形單元進行系統(tǒng)聚類，將地形單元聚類結(jié)果在ArcGIS軟件中顯示。

2.4 基于地形單元的Kriging空間插值方法

首先采用地統(tǒng)計學(xué)中普通克里金插值法對整個研究區(qū)有機質(zhì)進行空間插值，得到預(yù)測結(jié)果;然后單獨對不同地形單元下的土壤采樣點進行分層插值，利用耕地圖斑提取與之相對范圍內(nèi)采樣點的預(yù)測結(jié)果，得到土壤有機質(zhì)空間分布圖，再將各地形單元下的空間預(yù)測圖拼接，獲得整個縣域的空間分布預(yù)測圖。最后設(shè)置驗證點，各地形單元驗證點個數(shù)不能低于采樣點個數(shù)的5%，運用驗證點位上的土壤有機質(zhì)的實測值和預(yù)測值結(jié)果驗證分地形插值和直接插值拼接獲得的預(yù)測結(jié)果精度。

南非礦山工程師Krige最早提出了克里金插值方法，然后經(jīng)由法國地理數(shù)學(xué)家Matheron進行了插值方法的優(yōu)化[22]?？死锝鸩逯凳窃谟邢迏^(qū)域內(nèi)對區(qū)域化變量進行無偏最優(yōu)估計的一種方法[23]，這主要是由變異函數(shù)和結(jié)構(gòu)分析為理論支撐。

普通克里金法由于假設(shè)條件相對較少和需求參數(shù)較簡單，所以在克里金方法中是最為普遍的，公式為

Z*（X0）=∑ni=1λiZ（Xi）。（2）

式中：Z*（X0）為估算點土壤有機質(zhì)含量的預(yù)測值;λi為參與插值的位置對預(yù)測點的權(quán)重;Xi表示土壤采樣點的實際位置;Z（Xi）表示采樣點的實測值。

2.5 預(yù)測結(jié)果準(zhǔn)確性檢驗

本研究選取均方根誤差（RMSE）、平均預(yù)測誤差（ME）作為指標(biāo)，對于研究結(jié)果合理性的驗證是通過比較不同地形單元驗證點位的預(yù)測值與實測值之間的差異，利用2個指標(biāo)分析2種插值方法的差異和好壞，展開精度評價，從而驗證本研究提出方法的可行性。其中RMSE值越小，說明所預(yù)測精度越高，ME接近0，滿足預(yù)測的無偏性最優(yōu)，反之則說明預(yù)測精度低，不滿足預(yù)測的無偏性。

3 結(jié)果與分析

3.1 地形單元結(jié)果聚類分析

首先把11個地形單元的數(shù)據(jù)標(biāo)準(zhǔn)化，然后采用歐氏距離和組平均法進行系統(tǒng)聚類，將11個地形單元分成4類。溝谷單獨為一類;坡腳、陰坡、半陰坡、平地合并為坡腳類;陡坡、陽坡、半陽坡、坡肩合并為坡肩類;山脊和平頂峰合并為山脊類，如圖3所示。根據(jù)聚類結(jié)果結(jié)合實地數(shù)據(jù)把規(guī)律具有一致性的合并為一個單元，最終將研究區(qū)耕地劃分為4個地形單元（圖4）。根據(jù)調(diào)查與實地數(shù)據(jù)相符合。

3.2 土壤有機質(zhì)含量統(tǒng)計分析

依據(jù)地形單元劃分采樣點，土壤有機質(zhì)含量描述性統(tǒng)計特征見表1。區(qū)域內(nèi)土壤有機質(zhì)含量呈正態(tài)分布特征，由平均值和標(biāo)準(zhǔn)差可知，不同地形單元內(nèi)土壤有機質(zhì)含量離散度較相近;由偏度和峰度

可知不同地形單元的有機質(zhì)含量空間分布特征規(guī)律比較明顯，符合克里金插值的要求，也為下一步對土壤有機質(zhì)含量進行半方差函數(shù)結(jié)構(gòu)研究以及相應(yīng)地形單元插值提供了依據(jù)。

3.3 土壤有機質(zhì)空間自相關(guān)分析

在ArcGIS中運用空間統(tǒng)計模塊分別計算各地形單元內(nèi)和全局土壤有機質(zhì)含量的空間自相關(guān)程度，結(jié)果見表2。

由表2可知，各地形單元下土壤有機質(zhì)含量全局Morans I 均大于0，P均小于0.001，因此可以推斷縣域內(nèi)不同地形單元下的土壤有機質(zhì)含量的空間特征并不是隨機分布的，而是空間聚集性表現(xiàn)出非常顯著的特點，但是聚集效應(yīng)不同。土壤有機質(zhì)含量在各類地形單元下有明顯的空間正相關(guān)關(guān)系，其中在坡腳的空間聚集性最高，坡肩的聚集性最低，這也可能是由于坡肩坡度較大，土壤侵蝕嚴(yán)重，造成土壤有機質(zhì)流失。由以上分析可知，不同地形單元下的土壤有機質(zhì)含量空間分布規(guī)律不同，但是各單元特點都與實際情況相符，這也說明了基于劃分地形單元研究空間插值方法的可行性。

3.4 地統(tǒng)計空間分析及Kriging空間插值

利用地統(tǒng)計GS+軟件對各地形單元采樣點土壤有機質(zhì)含量進行半方差函數(shù)及擬合參數(shù)計算[22-24]，結(jié)果見表3。

表3是根據(jù)最優(yōu)擬合模型以決定系數(shù)大、殘差小為依據(jù)，比較各模型參數(shù)，選擇土壤有機質(zhì)含量半方差函數(shù)最優(yōu)擬合指數(shù)參數(shù)的[25]。根據(jù)塊金值的相關(guān)定義，原理上采樣點之間距離等于0時，半變異函數(shù)值也是0，由于存在測量誤差以及空間變異特征的影響，土壤樣點十分相近時，半變異函數(shù)值會發(fā)生變化，就會產(chǎn)生塊金值效應(yīng)，因此最小距離尺度下變量的變異特征及誤差是由塊金值反映的。由表3可知，采樣點整體、劃分地形單元的有機質(zhì)塊金值不一樣，說明劃分地形單元空間插值與全局插值具有差異性。

C是結(jié)構(gòu)方差，表示系統(tǒng)屬性區(qū)域變量最大空間變異，是由土壤屬性本身、海拔和氣候等結(jié)構(gòu)性因素引起的變異。由于該區(qū)域氣候差異小，土壤有機質(zhì)含量的空間變異特征基本是由于地形差異等結(jié)構(gòu)因素引起的。不同地形單元下結(jié)構(gòu)方差不同，所以地形因素是影響該區(qū)域土壤有機質(zhì)空間變異規(guī)律的主要原因。C/（Co+C）即表示塊金值與基臺值的比，被稱為空間相關(guān)度，該值的大小是系統(tǒng)變量的空間相關(guān)性的大小。根據(jù)Cambardella提出的劃分依據(jù)[26]，塊金基臺比小于25%，表明空間相關(guān)性強;大于75%則變量之間相關(guān)性較弱;處于二者之中，說明具有中等強度相關(guān)性。由表3可知，溝谷空間相關(guān)度大于75%，變異差異小，可能與結(jié)構(gòu)因素密切相關(guān)性較小，其他地形單元的空間相關(guān)度處于25%～75%之間，呈中等變異。

由表3可知，與全局有機質(zhì)含量相比，坡肩塊金值、基臺值都降低，變程也降低，這說明不同地形單元土壤有機質(zhì)含量的變化相對比較均勻，這也證明對該區(qū)域土壤有機質(zhì)含量進行預(yù)測時考慮地形因子的合理性，從另一方面說明了以地形單元對偏關(guān)縣土壤有機質(zhì)含量進行Kriging插值結(jié)果的有效性和可行性。

3.5 預(yù)測精度對比分析

基于劃分地形單元驗證點Kriging插值土壤有機質(zhì)含量實測值、預(yù)測值精度與全局插值精度分析結(jié)果進行比較。根據(jù)2種方法所得到的線性回歸方程系數(shù)，散點分布圖見圖5，劃分地形單元進行Kriging分層插值得到的預(yù)測值與實測值的決定系數(shù)（0.375 3）明顯高于全局Kriging插值法的決定系數(shù)（0.132 3）。

表4是對土壤有機質(zhì)含量空間特征分地形單元和全局預(yù)測合理結(jié)果的驗證。根據(jù)均方根誤差值越小精度越高，平均預(yù)測誤差接近于0，滿足預(yù)測值的無偏性[27]，得出劃分地形單元比全局預(yù)測精度高，而且插值結(jié)果的無偏性好。

3.6 土壤有機質(zhì)含量的空間分布特征

采用Kriging插值對不同地形單元土壤有機質(zhì)含量進行空間插值，根據(jù)地形單元耕地圖斑裁剪各地形單元內(nèi)的有機質(zhì)含量空間預(yù)測結(jié)果分布，并將其拼接得到整個研究區(qū)的空間分布見圖6-A，而不劃分地形單元采用全局采樣點直接插值得到的土壤有機質(zhì)含量空間分布結(jié)果見圖6-B。

由圖6可知，劃分地形單元Kriging插值得到耕地土壤有機質(zhì)含量范圍為4.15～14.93 g/kg，全局Kriging得到的含量為4.40～11.72 g/kg?？傮w來看，研究區(qū)2種方法預(yù)測結(jié)果范圍大體一致，空間分布格局基本相同，都呈現(xiàn)出偏東部和南部區(qū)域有機質(zhì)含量較高，東西方向偏低?？梢娡寥烙袡C質(zhì)在空間上的特征規(guī)律與實際地形相符合[28]，東西方向是溝谷，受人為影響因素比較大，所以有機質(zhì)含量預(yù)測結(jié)果偏低，這也進一步驗證了地形是引起空間變異存在的必然因素。

從局部細(xì)節(jié)分布特征來看，2種預(yù)測結(jié)果之間存在著一定的差異。經(jīng)過對比分析發(fā)現(xiàn)，劃分地形單元的預(yù)測結(jié)果更能展示細(xì)節(jié)信息，而且細(xì)節(jié)信息與地形單元的分布規(guī)律一致，邊界平滑且為漸變過渡;而全局Kriging插值得到的土壤有機質(zhì)含量空間分布特征比較明晰，但是局部本身屬性細(xì)節(jié)展示不充分，這顯然與研究區(qū)的實際情況有很大的差異，結(jié)果只能反映出整個研究區(qū)的大致情況，不能精確反映研究區(qū)局部實際情況，說明平滑效應(yīng)并沒有使不同地形單元下的土壤有機質(zhì)含量差異性變小。比如個別地形單元下耕地由于個別采樣點的土壤有機質(zhì)含量偏高，造成該區(qū)域的整體偏高，全局Kriging插值并不會詳細(xì)體現(xiàn)出高的地方。全局Kriging插值是基于采樣點的空間分布特征以及樣點之間的空間位置關(guān)系來進行插值，沒有根據(jù)地形單元對采樣點進行分類，所以空間插值結(jié)果會受距離較近的個別預(yù)測樣點的影響[29]，而劃分地形單元Kriging插值因為按照地形因素與土壤采樣點之間的關(guān)系進行分區(qū)，就可以很好地消除個別樣點對插值結(jié)果的誤差影響?？傊?，劃分地形單元Kriging插值的預(yù)測結(jié)果精度較高，更能反映土壤有機質(zhì)含量的空間信息以及局部區(qū)域分布情況，而直接Kriging插值忽略了土壤有機質(zhì)含量空間分布規(guī)律平穩(wěn)過渡性以及誤差影響。

4 討論與結(jié)論

對于地形條件復(fù)雜的丘陵山區(qū)，進行高密度的土壤采樣需要大量的人力物力，提前設(shè)計好有效土壤樣點采樣的方案能夠充分獲取土壤屬性的空間信息，還能大大降低成本[30]。所以綜合更多個地形因子分類組合精確劃分地形單元，研究各個地形單元變異規(guī)律，為局部樣點更新、最佳采樣、以及測土配方精準(zhǔn)施肥方案提供地理科學(xué)依據(jù)。本研究取得的主要研究結(jié)果如下：利用DEM數(shù)據(jù)結(jié)合坡度、坡向和地形位置指數(shù)綜合來劃分地形單元，基于地形單元，利用Kriging分層插值對復(fù)雜丘陵山區(qū)土壤有機質(zhì)含量進行空間預(yù)測，可更詳盡地表征出復(fù)雜山區(qū)的特定地形條件分異特征規(guī)律，將各個地形單元預(yù)測結(jié)果圖拼接為全域的空間分布圖。設(shè)置驗證點，并將該方法與全局Kriging插值法得到的預(yù)測結(jié)果進行比較。劃分地形單元與全局預(yù)測結(jié)果范圍和規(guī)律一致，但是預(yù)測精度高，而且展示了更詳細(xì)的局部區(qū)域空間細(xì)節(jié)信息。

由于研究區(qū)成土因素和母質(zhì)的影響，土壤有機質(zhì)含量的變化較小，所以與全局插值相比溝谷和坡肩的變程影響比較大?；诘匦螁卧耐寥烙袡C質(zhì)含量空間變異規(guī)律，是一種綜合地形因素的研究方法，更是為復(fù)雜的丘陵山區(qū)精確插值提供了一種參考。

參考文獻：

[1]Massey P A，Creamer R E，Schulte R，et al. The effects of earthworms，botanical diversity and fertiliser type on the vertical distribution of soil nutrients and plant nutrient acquisition[J]. Biology and Fertility of Soils，2013，49（8）：1189-1201.

[2]韓 丹，程先富，謝金紅，等. 大別山區(qū)江子河流域土壤有機質(zhì)的空間變異及其影響因素[J]. 土壤學(xué)報，2012，49（2）：403-408.

[3]張學(xué)雷，馮婉婉，鐘國敏.豫中褐土耕地土壤性質(zhì)空間分異及質(zhì)量評價[J]. 應(yīng)用生態(tài)學(xué)報，2011，22（1）：121-128.

[4]張 澤，呂 新，呂 寧，等. 基于GIS、RS的滴灌棉田土壤養(yǎng)分精確管理分區(qū)研究[J]. 農(nóng)業(yè)機械學(xué)報，2014，45（7）：125-132.

[5]Zhu H F，Ying Z，F(xiàn)eng N，et al. Relative influence of soil chemistry and topography on soil available micronutrients by structural equation modeling[J]. Journal of Soil Science and Plant Nutrition，2016，16（4）：1038-1051.

[6]Mabit L，Bernard C. Spatial distribution and content of soil organic matter in an agricultural field in eastern Canada，as estimated from geostatistical tools[J]. Earth Surface Processes and Landforms，2010，35（3）：278-283.

[7]付曉婷.基于地貌分區(qū)的土壤養(yǎng)分空間變異性分析及采樣設(shè)計研究——以陜西省藍田縣為例[D]. 西安：西北大學(xué)，2012.

[8]韋金麗，王國波，凌子燕.基于高分辨率DEM的地形特征提取與分析[J]. 測繪與空間地理信息，2012，35（1）：33-36.

[9]Yanai J，Mishima A，F(xiàn)unakawa S，et al. Spatial variability of organic matter dynamics in the Semi-Arid croplands of northern Kazakhstan[J]. Soil Science & Plant Nutrition，2005，51（2）：261-269.

[10]朱洪芬，南 鋒，徐占軍，等. 黃土高原盆地土壤有機質(zhì)與影響因子的空間多尺度關(guān)系[J]. 生態(tài)學(xué)報，2017，37（24）：8348-8360.

[11]趙業(yè)婷，常慶瑞，李志鵬，等. 渭北臺塬區(qū)耕地土壤有機質(zhì)與全氮空間特征[J]. 農(nóng)業(yè)機械學(xué)報，2014，45（8）：140-148.

[12]耿廣坡，高 鵬，呂圣橋，等. 魯中南山區(qū)馬蹄峪小流域土壤有機質(zhì)和全氮空間分布特征[J]. 中國水土保持科學(xué)，2011，9（6）：99-105.

[13]朱會義，劉述林，賈紹鳳.自然地理要素空間插值的幾個問題[J]. 地理研究，2004，23（4）：425-432.

[14]宋 軒，李立東，寇長林，等. 黃水河小流域土壤養(yǎng)分分布及其與地形的關(guān)系[J]. 應(yīng)用生態(tài)學(xué)報，2011，22（12）：3163-3168.

[15]田瑞云，王玉寬，傅 斌，等. 基于DEM的地形單元多樣性指數(shù)及其算法[J]. 地理科學(xué)進展，2013，32（1）：121-129.

[16]黃 魏，韓宗偉，羅 云，等. 基于地形單元的土壤有機質(zhì)空間變異研究[J]. 農(nóng)業(yè)機械學(xué)報，2015，46（4）：161-167.

[17]田惠文，張欣欣，畢如田，等. 非均勻網(wǎng)格環(huán)境下基于耕地質(zhì)量與集中連片度的基本農(nóng)田劃定[J]. 農(nóng)業(yè)工程學(xué)報，2020，36（3）：254-262.

[18]李 超，李文峰.高原耕地土壤養(yǎng)分空間分布與影響因子相關(guān)性研究[J]. 土壤通報，2014，45（5）：1113-1118.

[19]Tang C Q. Forest vegetation as related to climate and soil conditions at varying altitudes on a humid subtropical mountain，Mount Emei，Sichuan，China[J]. Ecological Research，2006，21（2）：174-180.

[20]Weiss A D. Topographic position and landforms analysis[C]. San Diego：ESRI User Conference，2001.

[21]黃 婷，周冀衡，李 強，等. 不同海拔高度植煙土壤pH值分布情況及其與土壤養(yǎng)分的關(guān)系——以云南省曲靖市為例[J]. 土壤通報，2015，46（1）：105-110.

[22]Matheron G. Principles of geostatistics[J]. Economic Geology，1963，58（8）：1246-1266.

[23]楊 東，程軍奇.甘肅黃土高原近52年四季極端降水變化特征

[J]. 中國農(nóng)學(xué)通報，2013，29（11）：180-187.

[24]張素梅，王宗明，張 柏，等. 利用地形和遙感數(shù)據(jù)預(yù)測土壤養(yǎng)分空間分布[J]. 農(nóng)業(yè)工程學(xué)報，2010，26（5）：188-194.

[25]張志霞，許明祥，師晨迪，等. 黃土丘陵區(qū)不同地貌單元土壤有機碳空間變異的尺度效應(yīng)[J]. 自然資源學(xué)報，2014，29（7）：1173-1184.

[26]Cambardella C A，Moorman T B，Novak J M，et al. Field-scale variability of soil properties in central Iowa soils[J]. Soil Science Society of America Journal，1994，58（5）：1501-1511.

[27]楊之江，陳效民，景 峰，等. 基于GIS和地統(tǒng)計學(xué)的稻田土壤養(yǎng)分與重金屬空間變異[J]. 應(yīng)用生態(tài)學(xué)報，2018，29（6）：1893-1901.

[28]劉 芳，塔西甫拉提·特依拜，依力亞斯江·努爾麥麥提，等. 準(zhǔn)東煤炭產(chǎn)業(yè)區(qū)周邊土壤重金屬污染與健康風(fēng)險的空間分布特征[J]. 環(huán)境科學(xué)，2016，37（12）：4815-4829.

[29]Bhunia G S，Shit P K，Chattopadhyay R. Assessment of spatial variability of soil properties using geostatistical approach of lateritic soil （West Bengal，India）[J]. Annals of Agrarian Science，2018，16（4）：436-443.

[30]Taye M，Simane B，Selsssie Y，et al. Analysis of the spatial variability of soil texture in a tropical highland：the case of the Jema watershed，northwestern highlands of Ethiopia[J]. International Journal of Environmental Research and Public Health，2018，15（9）：1903.