基于景觀及微地形特征的丘陵區(qū)土壤屬性預(yù)測

　　摘要為探討小流域尺度丘陵區(qū)的高分辨率數(shù)字土壤制圖方法，通過對景觀相分類的探索，配合應(yīng)用不同尺度的Geomorphons(GM)微地形特征數(shù)據(jù)構(gòu)成分類變量組參與高分辨率土壤pH、黏粒含量(SCC)和陽離子交換量(CEC)的預(yù)測制圖，并與傳統(tǒng)數(shù)字高程模型衍生變量和遙感變量進行組合與比較分析。此外，采用支持向量機、偏最小二乘回歸和隨機森林3種機器學(xué)習(xí)模型擇優(yōu)與殘差回歸克里金復(fù)合參與預(yù)測模型的構(gòu)建與評價。結(jié)果表明：景觀及多尺度微地形分類變量組的應(yīng)用分別提高小流域尺度丘陵地貌區(qū)pH、SCC和CEC預(yù)測精度的18.8%、8.2%和8.7%。包含植被信息的景觀相分類圖相比土地利用數(shù)據(jù)有更高的模型貢獻度;5m分辨率的GM微地形分類圖相比低分辨率的分類圖更適宜高精度的預(yù)測制圖。SCC使用隨機森林復(fù)合模型有最高的預(yù)測精度，而pH和CEC則不適宜在隨機森林模型的基礎(chǔ)上加入殘差回歸克里金模型。景觀-多尺度微地形分類變量、數(shù)字高程模型衍生變量和遙感變量三者結(jié)合的模型預(yù)測表現(xiàn)最佳，表明多元變量在起伏地形區(qū)域相比單一數(shù)據(jù)源能夠包含更多的土壤有效信息。由GM數(shù)據(jù)和地表景觀數(shù)據(jù)組成的景觀分類變量組作為主要變量能夠解釋小流域丘陵區(qū)部分土壤屬性約40%的空間變異。在同類型土壤預(yù)測制圖研究中，多分辨率GM及景觀分類數(shù)據(jù)有潛力作為環(huán)境變量參與預(yù)測模型的構(gòu)建。

　　關(guān)鍵詞景觀分類;微地形;數(shù)字土壤制圖;隨機森林;機器學(xué)習(xí)

　　土壤是一種基本的不可再生資源，在促進農(nóng)業(yè)可持續(xù)發(fā)展、維持地球生態(tài)圈穩(wěn)定和應(yīng)對氣候變化等方面都發(fā)揮著重要作用[1]。由于環(huán)境退化、糧食安全等問題日益嚴峻，農(nóng)業(yè)發(fā)展正面臨越來越大的壓力，對土壤生態(tài)系統(tǒng)的有效管理和保護需求也愈發(fā)迫切[2]。因此，如何獲取準確可靠的土壤狀況信息供給管理者做出適宜決策就顯得尤為重要。

　　近年來，隨著應(yīng)用模型的不斷發(fā)展，數(shù)字土壤制圖(digitalsoilmapping，DSM)與土壤學(xué)的有機結(jié)合越來越引起相關(guān)學(xué)者的重視[3]。DSM技術(shù)是基于土壤與其形成過程相關(guān)的環(huán)境因素之間的內(nèi)在關(guān)系構(gòu)建適宜的土壤-景觀模型，用于描繪土壤屬性的空間變異性[4]。如Scorpan模型作為Jenny方程的推廣，利用土壤信息、氣候、生物、地形、母質(zhì)以及其他因素引起的空間趨勢來預(yù)測土壤性質(zhì)，反映了土壤-景觀模型的一般性原理[5]。

　　城市景觀論文： 城市公共健康對未來住宅區(qū)景觀設(shè)計的影響

　　因此，一系列環(huán)境因子被納入DSM的建模過程中，總體可以分為兩類，一類為分類變量(又名定性變量)，如母質(zhì)類型、土壤類型和土地利用類型等;另一類為連續(xù)變量(又名定量變量)，如數(shù)字高程模型(digitalelevationmodel，DEM)衍生變量、植被指數(shù)和距離變化等[6–8]。在中國東部廣袤的丘陵山區(qū)，由于地形起伏程度較大，地形因素已經(jīng)成為土壤屬性制圖中最重要的預(yù)測因子[4]。而伴隨著高精度DEM的發(fā)展和普及，在小流域尺度上針對破碎化農(nóng)田管理進行較為準確的土壤屬性預(yù)測也基本可以實現(xiàn)[9]。因此，如何探索和應(yīng)用更適宜的多樣化地形變量，并將之與環(huán)境景觀要素相結(jié)合以構(gòu)建高精度的土壤-景觀模型將成為丘陵山區(qū)DSM研究的重要方向。其中，土壤質(zhì)地和土地利用類型作為土壤環(huán)境和人類活動的綜合反映，是小流域尺度上土壤屬性空間分布的重要控制因子。

　　之前的一些研究也表明，不同土壤質(zhì)地和土地利用方式之間的pH、有機碳、總氮等土壤屬性有顯著差異[10-11]。相對于側(cè)重人為因素影響的土地利用圖，本研究將高程、植被覆蓋、土壤質(zhì)地和土地利用方式相結(jié)合構(gòu)成景觀相分類變量，藉以更好地表達區(qū)域景觀的空間異質(zhì)性[12]。地形分類作為地上部特征的直接反映，是對區(qū)域地形信息的有力補充。相比平原區(qū)域，在以地形特征為主導(dǎo)預(yù)測變量的丘陵區(qū)DSM研究中，地形分類數(shù)據(jù)往往能提供更多的潛在土壤信息，但迄今關(guān)于小流域尺度上的DSM研究仍然較少[13]。

　　近40年來，適用于精確地形地貌分類的計算機模型在不斷地推陳出新[14]，盡管不同的模型因為算法差異會獲得不同的地形分類結(jié)果，但是大部分的算法內(nèi)核都基于相同的假設(shè)而成立，即存在兩個可比較的相鄰地形要素，且在地理信息系統(tǒng)中具有可識別可 捕捉的特征(如坡向、剖面曲率等)[15]。本研究使用了一種較為新穎的Geomorphons(GM)地形分類方法，于2013年由Jasiewicz等提出[16]。該方法主要根據(jù)目標位置的高程值在指定大小的鄰域搜索視窗內(nèi)所捕捉的相對差值來確定該處像元所對應(yīng)的地形元素類型。GM方法相對于其他地形分類模型有2個顯著特征：

　　1)能在獲得擁有準確定義的地形元素分類結(jié)果的前提下僅占用少量計算資源，更加高效;2)相比其他模型可選擇地形屬性集作為判定要素的主觀性，GM方法能有效避免分類結(jié)果的不確定性且對局部地形特征的識別擁有更高的自適應(yīng)性[17]。尺度是所有地形分類模型的關(guān)鍵參數(shù)，地形格局在不同尺度上可能存在很大差異[18]。多尺度方法有助于提高地形分類圖作為預(yù)測因子的應(yīng)用效率，并較為全面地表征復(fù)雜的地表特性。在預(yù)測模型方面，使用較為常見的3種機器學(xué)習(xí)算法：支持向量機(supportvectormachine，SVM)、偏最小二乘回歸(partialleastsquaresregression，PLSR)和隨機森林(randomforest，RF)。

　　此外，回歸克里金與最優(yōu)算法的結(jié)合也加入到模型的構(gòu)建討論中，以獲得最高的預(yù)測精度。本研究將景觀相、土地利用、土壤質(zhì)地聯(lián)合多尺度GM微地形圖作為景觀-地形分類變量組，與傳統(tǒng)DEM衍生變量組和遙感變量組形成對照，共同參與到小流域尺度丘陵區(qū)的pH、土壤黏粒含量(soilclaycontent，SCC，<0.002mm)和陽離子交換量(cationexchangecapacity，CEC)預(yù)測制圖中。

　　1研究地區(qū)與研究方法

　　1.1研究區(qū)概況

　　研究區(qū)位于江蘇省鎮(zhèn)江市黃梅鎮(zhèn)北部(32°3'20″—32°5'50″N，119°10'20″—119°12'30″E)。該區(qū)域北鄰長江，平均高度82m，北部地勢較高，最高225m，且丘陵地形突出，最大坡度69°，平均坡度12°。占地面積約5.37km2，屬北亞熱帶暖濕氣候，年均氣溫15.2℃，年均降水量1060mm，其中，70%~80%集中在6—8月。研究區(qū)丘陵地形特征主要體現(xiàn)于北部地區(qū)，而南部較平坦區(qū)域覆蓋旱地和灌溉農(nóng)業(yè)，以小麥、水稻等經(jīng)濟作物為主。

　　丘陵主體以石英砂巖巖性為主，優(yōu)勢樹種為馬尾松(Pinusmassoniana)和白櫟(Quercusfabri)。丘陵坡麓部覆蓋大片下蜀黃土，以下蜀黃土為母質(zhì)的土壤廣泛分布于長江中下游地區(qū)[19]。根據(jù)中國土壤系統(tǒng)分類檢索[20]，研究區(qū)共分布有人為土、淋溶土、雛形土和新成土4個土綱，且由于長期的稻麥輪作歷史，水耕人為土在南部農(nóng)耕區(qū)的分布較為廣泛，而淋溶土則在北部丘陵區(qū)占據(jù)主導(dǎo)地位。

　　1.2樣品采集與分析

　　為盡可能捕捉研究區(qū)內(nèi)的景觀異質(zhì)性，以景觀類型差異最大化為原則在研究區(qū)內(nèi)共采集74個土壤樣品。采樣設(shè)計以景觀分類為基礎(chǔ)，綜合與土壤形成過程密切相關(guān)的環(huán)境信息。對土壤景觀格局的分布特征進行綜合性室內(nèi)分析，最終形成的采樣點空間分布對研究區(qū)主要景觀類型如自然林地、灌木林地、旱地和水田等進行全面覆蓋。采用五點混合取樣法在1m2內(nèi)獲得每個表層土壤(0~20cm)樣品，在室內(nèi)溫度20~25℃條件下風(fēng)干14d，剔除植物殘根和石礫，經(jīng)研磨后采用2mm不銹鋼網(wǎng)篩進行篩分。土壤pH采用水浸提法(水土比2.5:1)進行測量，土壤粒徑組成采用干篩法經(jīng)過篩分和稱重獲得，土壤CEC采用乙酸鈉-火焰光度法進行測量[21]。

　　1.3景觀分類變量

　　景觀分類變量包括土壤質(zhì)地分類變量、土地利用分類變量以及景觀相分類變量。景觀相通常作為1：10000及更大比例制圖的分類等級，主要以自然或半自然區(qū)域作為分類對象，將微地貌作為基礎(chǔ)分類的優(yōu)勢因子能夠較好地反映小區(qū)域的景觀差異。對于微地貌的劃分指標因區(qū)域地質(zhì)條件、景觀特征、立地方位以及相對高程的變化等因素因地制宜加以選擇。

　　針對小流域尺度下的丘陵地形特征，將高程、坡度、土壤質(zhì)地和土地覆被類型作為分類依據(jù)：相對高度包括緩丘(20~50m)、淺丘(50~100m)、深丘(100~200m)和低山(200~500m);坡度參照包括平地(≤2°)、緩坡(2°~6°)、中緩坡(6°~15°)、中坡(15°~25°)和陡坡(≥25°);土壤質(zhì)地分為黏壤土(暗紅棕色)、粉砂壤土(棕色)、砂質(zhì)壤土(暗棕色)和壤土(紅棕色);土地覆被包括水田、旱地、茶園、桃園、馬尾松-杉木混交林群落、白櫟林群落、村莊聚落和未利用地。根據(jù)以上分類標準經(jīng)過適當(dāng)?shù)膱D斑合并和舍棄，簡化輪廓線并進行適當(dāng)?shù)念愋蜌w并，生成研究區(qū)1:25000景觀相分類圖。

　　1.4多尺度GM微地形分類變量

　　本研究使用的微地形分類變量通過GM方法獲得，打破了標準演算方法在固定尺寸的單個窗口中檢測所有地形元素的局限性。綜合可見鄰域內(nèi)的高程信息并據(jù)此判定局域幾何形態(tài)特征，最終生成10種主要的地形部位[16]。GM地形分類法有2個主要參數(shù)：檢索窗口的搜索半徑，根據(jù)以往研究，采用20像元的設(shè)置能獲得最佳分類結(jié)果[22];平坦度閾值，即視作水平區(qū)域的坡度閾值大小，由于地處丘陵，研究區(qū)平均坡度較大，因此將平坦度閾值設(shè)置為5°能獲得較好的“平地”分布，對應(yīng)于實際的耕地分布。

　　在預(yù)測模型的構(gòu)建過程中，地形分類變量對尺度的變化較敏感[23]，而預(yù)測變量與土壤屬性值之間的關(guān)系也會隨著尺度的變化而變化。尺度通常被認為是一個特定的空間分辨率，需參考輸出的制圖單元大小，受到遙感數(shù)據(jù)像元大小、空間異質(zhì)性和數(shù)據(jù)量的影響[24-25]。

　　本研究將不同尺度的GM地形分類變量相結(jié)合，探討在何種尺度下可以獲得最高的模型貢獻度和最佳的預(yù)測效果。采用通過立體像對生成的5、15、30mDEM與常用的30m分辨率ShuttleRadarTopographyMission(SRTM)DEM作為數(shù)據(jù)來源，生成多分辨率GM微地形分類圖，共同參與pH、SCC和CEC的預(yù)測制圖，其中，15、30mDEM由5mDEM數(shù)據(jù)在ArcGIS10.6中經(jīng)重采樣生成。

　　1.5遙感變量和DEM衍生變量

　　遙感影像數(shù)據(jù)源于中國資源衛(wèi)星應(yīng)用中心，包括高分二號(GF-2)多光譜衛(wèi)星影像和高分三號(GF-3)雷達衛(wèi)星影像。GF-2影像獲取時間為2019年9月，該時段研究區(qū)所處流域有較高的綠度覆蓋，適宜提取相關(guān)植被指數(shù)。GF-3數(shù)據(jù)為C頻段多極化合成孔徑雷達(SAR)標準條帶影像，采用VH、VV雙極化，幅寬130km，空間分辨率25m，成像時間為2020年1月。

　　GF-2影像的預(yù)處理在ENVI5.3軟件中完成[26]，其近紅外波段(NIR)(0.77~0.89μm)及兩種植被光譜指數(shù)被提取作為環(huán)境預(yù)測變量，包括歸一化植被指數(shù)(normalizeddifferencevegetationindex，NDVI)和花青素反射指數(shù)2(anthocyaninreflectanceindex2，ARI2)。

　　GF-3雷達衛(wèi)星影像使用PIE-SAR6.0進行預(yù)處理[27]，VH和VV數(shù)據(jù)被提取作為環(huán)境變量，分辨率為25m。本研究使用5m空間分辨率DEM進行衍生地形環(huán)境變量的提取，通過資源三號02星(ZY3-02)前后視立體像對在ENVI5.3.1軟件中利用DEM提取模塊生成，經(jīng)100個呈網(wǎng)格狀分布的實測高程值驗證后的精度誤差為1.59m，符合5mDEM的應(yīng)用精度要求[28]。所需DEM衍生變量由ArcGIS10.6和SAGAGIS軟件提取并參與后續(xù)的預(yù)測制圖。

　　2結(jié)果與分析

　　2.1土壤屬性統(tǒng)計特征及分布

　　土壤pH、SCC和CEC的偏度值分別為0.05、0.11和1.01，表現(xiàn)出不同程度的右偏趨勢。其中，pH值基本呈對稱分布，SCC值表現(xiàn)出輕微偏態(tài)分布，而CEC值表現(xiàn)出高度偏態(tài)分布，3種土壤屬性數(shù)據(jù)的標準差均低于平均值。根據(jù)Essington[36]提出的以變異系數(shù)(coefficientofvariation，CV)評估土壤屬性的均質(zhì)性和異質(zhì)性的劃分標準，即CV值低于15%為均勻分布，16%~35%為中度變異，高于36%為高度變異，則研究區(qū)內(nèi)pH值分布較為均勻，SCC和CEC則有不同程度的中度變異。

　　由于樣點個數(shù)有限，為更好地表現(xiàn)3種土壤屬性在地形起伏區(qū)域的差異，將原5m分辨率GM圖斑按地形部位相鄰原則進行聚合，具體為山峰、山脊和山肩聚合為坡頂部，凸背坡、直背坡和凹背坡聚合為坡中部，平區(qū)和山腳聚合為平坦部，溝谷和洼地聚合為凹坡部。由此可以看出，除平坦部和坡中部的pH外(P=0.034)，3種土壤屬性在4種地形部位上均表現(xiàn)出顯著差異，表明高分辨率的GM地形分類圖斑在小流域尺度的丘陵地帶能較好地表征pH、SCC和CEC之間的差異，進而應(yīng)用于土壤屬性的預(yù)測制圖。

　　2.2模型的建立與評價

　　將應(yīng)用的環(huán)境變量分為3類進行不同的組合加入到3種預(yù)測模型中，3類環(huán)境變量分別為：包含景觀變量和GM多尺度微地形變量的分類變量、DEM衍生連續(xù)變量以及從光學(xué)和SAR遙感影像中提取的連續(xù)變量，共形成5種變量組合(Ⅰ：遙感變量;Ⅱ：景觀及微地形分類變量;Ⅲ：DEM衍生變量;Ⅳ：景觀及微地形分類變量+DEM衍生變量;Ⅴ：遙感變量+景觀及微地形分類變量+DEM衍生變量)，并對其在SVM、PLSR和RF模型中的預(yù)測表現(xiàn)進行評估。RF模型對pH、SCC和CEC的預(yù)測在3種機器學(xué)習(xí)模型中都擁有最高的預(yù)測精度，在此基礎(chǔ)上對RF的預(yù)測殘差空間結(jié)構(gòu)進行回歸克里金插值(regressionKriging，RK)，獲得RF-RK預(yù)測模型精度。對預(yù)測精度的比較結(jié)果表明，3種機器學(xué)習(xí)模型在不同土壤屬性的預(yù)測中表現(xiàn)出一定的特征和差異。

　　以模組Ⅴ為例，在pH的預(yù)測結(jié)果中，模型精度由高到低依次為RF、PLSR和SVM(R2分別為0.63、0.58和0.54);在SCC的預(yù)測中，SVM的精度最低，PLSR和RF模型的預(yù)測結(jié)果相似(R2分別為0.53、0.58和0.59);對于CEC的預(yù)測，RF模型仍擁有最佳的預(yù)測結(jié)果，SVM模型的精度高于PLSR模型(R2分別為0.51、0.50和0.48)。綜合來看，相比SVM和PLSR模型，RF模型擁有總體最佳的預(yù)測表現(xiàn)，RF模型在DSM領(lǐng)域的應(yīng)用也得到過多次驗證[37]。

　　對RF模型的預(yù)測殘差進行回歸克里金插值并疊加以評估其對預(yù)測結(jié)果的影響。土壤pH由于其殘差過小導(dǎo)致預(yù)測精度并未發(fā)生明顯變化。SCC的預(yù)測精度明顯提升，CEC的預(yù)測精度低于單一RF模型，其中，SCC殘差的最佳擬合模型為高斯模型，塊金值為167.3，基臺值為334.7，塊金系數(shù)為50%，表明模型預(yù)測殘差具有一定的空間結(jié)構(gòu)性。模型的殘差可以作為預(yù)測不確定性的一種表征，當(dāng)土壤屬性呈現(xiàn)較強的空間自相關(guān)性時，殘差空間結(jié)構(gòu)與模型的結(jié)合能夠優(yōu)化模型的不確定性，使預(yù)測結(jié)果得到有效提升。

　　3討論

　　Odeh等[38]研究表明，殘差空間結(jié)構(gòu)回歸克里金應(yīng)用于不同土壤屬性的效果有較大差異，這種差異與目標屬性和地形因素之間的關(guān)聯(lián)性有關(guān)。他們的結(jié)果顯示，表土層礫石含量及亞表土層黏粒含量在應(yīng)用殘差回歸克里金模型時有較大提升，這與本研究結(jié)果一致，表明機器學(xué)習(xí)與殘差回歸克里金的復(fù)合模型適用于土壤質(zhì)地等受地形特征影響較大的土壤屬性的預(yù)測與制圖。

　　土地利用分類數(shù)據(jù)在大尺度DSM研究中已逐漸得到普及，而針對土壤調(diào)查進行劃分的景觀分類圖的應(yīng)用則鮮有涉及。小流域尺度的土壤分布與所處景觀條件有著密切聯(lián)系[39]，尺度相符且詳略得當(dāng)?shù)牡乇砀脖环诸愋畔⒂兄�于小尺度的土壤分類或相關(guān)屬性值的預(yù)測分布。GM地形分類在DSM領(lǐng)域的應(yīng)用則多見于土壤分類的研究，如Ngunjiri等[22]在肯尼亞UasinGishu高原利用GM地形分類對當(dāng)?shù)赝寥李愋瓦M行預(yù)測，總體分類精度達到了0.58。

　　由于地形分類對于尺度效應(yīng)的敏感性，不同分辨率的GM圖對于土壤景觀格局的可解釋性也有所區(qū)別，如Atkinson等[23]研究發(fā)現(xiàn)，不同尺度GM圖的地形元素對應(yīng)的土壤類型、黏粒含量和土壤深度存在顯著差異。盡管上述研究討論了土壤分布與GM地形分類特征之間的關(guān)系，卻并未將其應(yīng)用于土壤屬性的預(yù)測。

　　本研究基于土壤與所處景觀的內(nèi)在關(guān)系，在結(jié)合傳統(tǒng)DEM提取數(shù)據(jù)作為環(huán)境變量的基礎(chǔ)上，輔以高精度遙感信息數(shù)據(jù)，探討了景觀分類數(shù)據(jù)與多尺度GM微地形分類圖等類別變量在小流域尺度丘陵地貌區(qū)的應(yīng)用潛力。結(jié)果表明，詳細的地表景觀信息結(jié)合高分辨率的GM數(shù)據(jù)能夠有效提高土壤pH、SCC和CEC的預(yù)測精度。主要結(jié)論如下：

　　1)RF模型相比SVM和PLSR模型在研究區(qū)內(nèi)針對pH、SCC和CEC三種土壤屬性擁有更好的預(yù)測表現(xiàn)。對于SCC值的預(yù)測，RF模型與殘差回歸克里金插值結(jié)合形成的復(fù)合模型的預(yù)測精度最佳。2)在3種土壤屬性中，pH的制圖精度高于SCC和CEC，CEC制圖的不確定性最高。

　　3種類型環(huán)境變量共同組成的模組Ⅴ最高能解釋

　　pH值空間變異的63%、SCC值的61%和CEC值的51%。3)景觀類型數(shù)據(jù)與多尺度GM微地形分類數(shù)據(jù)組合而成的景觀-微地形分類變量組的應(yīng)用在RF模型中對土壤pH、SCC和CEC的預(yù)測精度相比單一類型的DEM衍生變量組能分別提升18.8%、8.2%和8.7%。4)在景觀類變量中，包含植被分布信息的景觀相變量相比單一土地利用變量擁有更高的模型貢獻率;在多尺度GM微地形分類變量中，5m高分辨率GM分類數(shù)據(jù)表現(xiàn)最佳。

　　參考文獻

　　[1]AmundsonR,BerheAA,HopmansJW,etal.Soilandhumansecurityinthe21stcentury.Science,2015,348:1261071

　　[2]MontanarellaL,PennockDJ,McKenzieN,etal.World’ssoilsareunderthreat.Soil,2016,2:79–82

　　[3]ScullP,FranklinJ,ChadwickOA,etal.Predictivesoilmapping:Areview.ProgressinPhysicalGeography:EarthandEnvironment,2003,27:171–197

　　[4]CameraC,ZomeniZ,NollerJS,etal.AhighresolutionmapofsoiltypesandphysicalpropertiesforCyprus:Adigitalsoilmappingoptimization.Geoderma,2017,285:35–49

　　[5]McBratneyA,MendonçaSM,MinasnyB.Ondigitalsoilmapping.Geoderma,2003,117:3–52

　　[6]Pahlavan-RadMR,AkbarimoghaddamA.SpatialvariabilityofsoiltexturefractionsandpHinafloodplain(casestudyfromeasternIran).Catena,2018,160:275–281

　　[7]DongW,WuT,LuoJ,etal.Landparcel-baseddigitalsoilmappingofsoilnutrientpropertiesinanalluvial-diluviaplainagriculturalareainChina.Geoderma,2019,340:234–248

　　作者：1魏宇宸1趙美芳2朱昌達1張秀秀1潘劍君1

轉(zhuǎn)載請注明來自發(fā)表學(xué)術(shù)論文網(wǎng)：http:///jzlw/28823.html