基于GIS的寧夏土壤侵蝕敏感性與景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)

鄭雪慧, 楊 志, 任正龑, 魏小燕, 程金花

(1.北京林業(yè)大學(xué) 水土保持學(xué)院, 北京 100083; 2.寧夏水土保持監(jiān)測(cè)總站, 銀川 750021)

土壤侵蝕導(dǎo)致區(qū)域生態(tài)環(huán)境惡化與生態(tài)平衡失調(diào)，嚴(yán)重阻礙社會(huì)經(jīng)濟(jì)的可持續(xù)發(fā)展，一直是全球普遍關(guān)注的環(huán)境問(wèn)題[1-2]。近年來(lái)，隨著“3S”技術(shù)的發(fā)展，使突破徑流小區(qū)[3]、坡面[4]尺度開展更大尺度的土壤侵蝕評(píng)價(jià)研究成為熱點(diǎn)，基于USLE，RUSLE[5]、CSLE等[6]土壤侵蝕經(jīng)驗(yàn)?zāi)Ｐ偷耐寥狼治g敏感性評(píng)價(jià)應(yīng)用廣泛。

土壤侵蝕敏感性評(píng)價(jià)是以土壤侵蝕方程為基礎(chǔ)來(lái)識(shí)別自然條件下容易形成土壤侵蝕的區(qū)域[7]。國(guó)外學(xué)者Terranova[8]，Nekhay[9]等借助GIS和RS技術(shù)和通用土壤侵蝕方程分別在地中海區(qū)域和西班牙山區(qū)開展土壤侵蝕定量評(píng)價(jià)研究；國(guó)內(nèi)學(xué)者龔雪梅[10]、楊青林[11]、張恩偉[7]等也分別在岷江上游、川西北地區(qū)、滇池流域開展土壤侵蝕敏感性的研究。另一方面，隨著景觀生態(tài)學(xué)[12]發(fā)展，景觀格局指數(shù)逐漸應(yīng)用到土壤侵蝕的研究中。景觀生態(tài)風(fēng)險(xiǎn)[13]是根據(jù)土地利用的空間格局探究人類社會(huì)因素所造成的生態(tài)環(huán)境風(fēng)險(xiǎn)。高艷麗等[1]在漢江流域開展景觀格局變化對(duì)土壤侵蝕影響的研究；王計(jì)平等[14]在黃土丘陵區(qū)開展景觀格局對(duì)土壤侵蝕過(guò)程影響的研究，發(fā)現(xiàn)流域景觀格局和土壤侵蝕過(guò)程特征存在明顯空間分異和相對(duì)差異；劉希朝等[15]基于土地利用數(shù)據(jù)構(gòu)建黃河流域生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)模型，探究景觀格局變化下的生態(tài)風(fēng)險(xiǎn)變化趨勢(shì)。土壤侵蝕是過(guò)程，土地利用空間分布為格局，格局決定過(guò)程，同時(shí)過(guò)程也會(huì)促進(jìn)格局的改變[16]。目前，對(duì)于土壤侵蝕敏感性和景觀生態(tài)風(fēng)險(xiǎn)的研究大多是獨(dú)立展開的，無(wú)法識(shí)別自然和社會(huì)綜合因素下的土壤侵蝕風(fēng)險(xiǎn)。因此，將兩者結(jié)合構(gòu)建土壤侵蝕和景觀綜合生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)體系，開展區(qū)域土壤侵蝕及其景觀生態(tài)風(fēng)險(xiǎn)研究，不僅能夠綜合反映自然與人為相互作用下對(duì)生態(tài)環(huán)境產(chǎn)生的不利影響，而且能夠?yàn)閰^(qū)域精準(zhǔn)化識(shí)別土壤侵蝕的敏感區(qū)和進(jìn)行高質(zhì)量的綜合治理提供理論支撐。

寧夏全境位于黃河流域，是我國(guó)西北的綠色屏障，承擔(dān)著維護(hù)西北乃至全國(guó)生態(tài)安全的重要使命。據(jù)《2019寧夏回族自治區(qū)水土保持公報(bào)》顯示，目前寧夏仍有土壤侵蝕面積1.59×104km2，占國(guó)土面積的23.93%。為了探究寧夏全區(qū)潛在土壤侵蝕區(qū)域的空間差異，識(shí)別土壤侵蝕敏感性區(qū)域，在寧夏開展土壤侵蝕敏感性和景觀生態(tài)風(fēng)險(xiǎn)的綜合研究，實(shí)現(xiàn)對(duì)寧夏全區(qū)土壤侵蝕情況的基礎(chǔ)摸排，為推進(jìn)流域精準(zhǔn)治理提供科學(xué)支撐。同時(shí)，為寧夏水土保持規(guī)劃、水土保持措施布設(shè)提供參考依據(jù)，對(duì)寧夏實(shí)現(xiàn)黃河流域生態(tài)保護(hù)與高質(zhì)量發(fā)展先行區(qū)的目標(biāo)有重要意義。

1 材料與方法

1.1 研究區(qū)概況

寧夏回族自治區(qū)總面積6.64×104km2，地處黃河上游及沙漠與黃土高原的交接地帶，東、西、北三面被毛烏素、騰格里、烏蘭布和三大沙漠包圍，是我國(guó)西北重要的生態(tài)屏障區(qū)。地理位置介于35°14′—39°23′N，104°17′—107°39′E，地勢(shì)從西南向東北逐漸傾斜，平均海拔1 000 m以上。地貌復(fù)雜，山地迭起，盆地錯(cuò)落，大體可分為：黃土高原，鄂爾多斯臺(tái)地，洪積沖積平原和六盤山、羅山、賀蘭山南北中三段山地。屬溫帶大陸性干旱半干旱氣候，年均氣溫為6～10℃，年均降水量南多北少，大都集中在夏季，年降雨量在300 mm左右，蒸發(fā)量卻在1 000 mm以上。除黃河干流過(guò)境外，還有清水河、苦水河、葫蘆河、涇河、祖厲河等主要河流。土壤類型主要為灰漠土、灌淤土、綿土、灰褐土、栗鈣土等。全區(qū)草原綜合植被蓋度為53.5%，森林覆蓋率為15.2%。土地利用以草地、耕地為主，兩者合計(jì)面積占比為69.7%[17]。

1.2 數(shù)據(jù)來(lái)源

氣象數(shù)據(jù)源于中國(guó)氣象數(shù)據(jù)網(wǎng)(http:∥data.cma.cn)下載的中國(guó)地面累年值日值數(shù)據(jù)集；土壤質(zhì)地?cái)?shù)據(jù)來(lái)源于中國(guó)科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心(http:∥www.resdc.cn)；DEM數(shù)據(jù)(Digital elevation model，數(shù)字高程數(shù)據(jù))為地理空間數(shù)據(jù)云(http:∥www.gscloud.cn)提供的30 m分辨率原始高程數(shù)據(jù)；寧夏行政區(qū)劃、土地利用、植被覆蓋等基礎(chǔ)數(shù)據(jù)由寧夏水土保持監(jiān)測(cè)總站提供。其中，土地利用數(shù)據(jù)解譯所需的遙感影像時(shí)間為2019年，空間分辨率30 m；植被覆蓋度所需的遙感影像數(shù)據(jù)為2019年MODIS歸一化植被指數(shù)(NDVI)產(chǎn)品，時(shí)間分辨率為每半月1期，共24期，空間分辨率為30 m。

1.3 研究方法

1.3.1 土壤侵蝕敏感性評(píng)價(jià) 基于RUSLE通用土壤流失方程中降雨侵蝕力、土壤可蝕性、地形起伏度和植被覆蓋度因子構(gòu)建土壤侵蝕敏感性評(píng)價(jià)方程[18](公式1)，經(jīng)柵格計(jì)算獲得評(píng)價(jià)結(jié)果。

(1)

式中:SSi為土壤侵蝕敏感性指數(shù)；Ri為降雨侵蝕力；Ki為土壤質(zhì)地因子；LSi為地形起伏度；Ci為植被覆蓋度。各評(píng)價(jià)因子的獲取方法以及具體分級(jí)標(biāo)準(zhǔn)[19]見表1—2。

表1 土壤侵蝕敏感性評(píng)價(jià)因子獲取方法

表2 土壤侵蝕敏感性評(píng)價(jià)因子分級(jí)標(biāo)準(zhǔn)

1.3.2 景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià) 景觀格局的變化是決定土壤侵蝕程度的關(guān)鍵因素[15]，診斷寧夏景觀格局及生態(tài)風(fēng)險(xiǎn)對(duì)促進(jìn)黃河流域高質(zhì)量發(fā)展具有重要意義。景觀生態(tài)風(fēng)險(xiǎn)由外部的擾動(dòng)強(qiáng)度和內(nèi)部的脆弱性來(lái)衡量，景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)包括景觀格局指數(shù)選取、生態(tài)風(fēng)險(xiǎn)小區(qū)劃分、生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)模型構(gòu)建3個(gè)過(guò)程。首先，在參照相關(guān)景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)研究[23-24]和綜合考慮寧夏實(shí)際情況及各景觀格局指數(shù)的生態(tài)學(xué)意義的基礎(chǔ)上，選取景觀破碎度(Ci)、景觀分離度(Ni)、景觀分維數(shù)(FDi)、景觀干擾度(Si)、景觀脆弱度(Fi)、景觀損失度(Ri)共6種景觀格局指數(shù)。各景觀指數(shù)的數(shù)值可由預(yù)處理后的土地利用數(shù)據(jù)經(jīng)Fragstats 4.2軟件計(jì)算獲得，具體計(jì)算公式和參數(shù)含義見表3。其次，在充分考慮寧夏景觀空間異質(zhì)性、斑塊大小和面積的基礎(chǔ)上，利用ArcGIS 10.2軟件數(shù)據(jù)管理模塊下的Create Fishnet功能，將研究區(qū)劃分為10 km×10 km的網(wǎng)格單元，共劃分風(fēng)險(xiǎn)小區(qū)617個(gè)。最終，根據(jù)景觀干擾度指數(shù)和景觀脆弱度指數(shù)構(gòu)建評(píng)價(jià)模型，計(jì)算得到各小區(qū)景觀生態(tài)風(fēng)險(xiǎn)指數(shù)ERIk(公式2)，并采用ArcGIS 10.2的Geostatistical Analyst模塊下Kriging插值對(duì)小區(qū)樣點(diǎn)的風(fēng)險(xiǎn)值進(jìn)行球狀擬合，獲得研究區(qū)的景觀生態(tài)風(fēng)險(xiǎn)指數(shù)。

表3 景觀格局指數(shù)計(jì)算方法

(2)

式中:ERIk為評(píng)價(jià)單元k的生態(tài)風(fēng)險(xiǎn)指數(shù);n為景觀類型的數(shù)量;Aki為景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)單元k中第i類景觀的面積(km2);Ak為景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)單元k的面積(km2);Ri為景觀損失度指數(shù)。

1.3.3 土壤侵蝕景觀綜合生態(tài)風(fēng)險(xiǎn)評(píng)價(jià) 基于土壤侵蝕敏感性和景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果構(gòu)建區(qū)域綜合生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)模型(SERI)，即土壤侵蝕敏感性指數(shù)與景觀損失度指數(shù)的關(guān)系特征函數(shù)式(公式3)[18]。運(yùn)用ArcGIS 10.2中柵格計(jì)算器計(jì)算，并采用自然分界法劃分為5個(gè)風(fēng)險(xiǎn)等級(jí)，獲得寧夏綜合生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果。

SERI=SSi·ERIk

(3)

式中:SERI為土壤侵蝕的景觀生態(tài)風(fēng)險(xiǎn)值；SSi和ERIk同上。

2 結(jié)果與分析

2.1 寧夏風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果的空間分布

總的來(lái)說(shuō)，寧夏土壤侵蝕敏感性等級(jí)以無(wú)敏感和輕度敏感為主，二者共計(jì)面積占比在60%以上；中度及以上敏感區(qū)隨敏感性等級(jí)升高面積占比減小。其中，無(wú)敏感區(qū)面積19 555 km2，面積占比為29.45%；輕度敏感區(qū)面積20 817.88 km2，面積占比為31.35%；中度敏感區(qū)面積13 371.70 km2，面積占比為20.14%；高度敏感區(qū)面積8 163.35 km2，面積占比為12.29%；極敏感區(qū)面積4 492.07 km2，面積占比為6.77%。此外，從空間格局上看(圖1A)，高度、極敏感區(qū)主要分布在寧夏南部山區(qū)；中度敏感區(qū)分布相對(duì)分散，主要分布在鹽池、海原、沙坡頭以及大武口區(qū)等地；無(wú)敏感和輕度敏感區(qū)則主要分布在中部干旱帶和北部的引黃灌區(qū)。

由圖1B可知，寧夏景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)空間差異顯著，較低風(fēng)險(xiǎn)區(qū)面積占比最高，達(dá)28.10%，主要分布在中部的鹽池、利通等地；中等風(fēng)險(xiǎn)區(qū)面積占比次之，為25.42%，主要分布在北部平羅周邊、中部中寧海原等地；較高風(fēng)險(xiǎn)區(qū)面積占比第三，為21.22%，分布在原州、同心以及金鳳區(qū)周圍；高風(fēng)險(xiǎn)和低風(fēng)險(xiǎn)面積占比較小且接近，分別為12.64%，12.62%，其中低風(fēng)險(xiǎn)區(qū)在大武口、靈武、沙坡頭西北部等地，高風(fēng)險(xiǎn)區(qū)則主要集中在彭陽(yáng)、西吉。

整體看來(lái)(圖1C)，土壤侵蝕與景觀生態(tài)風(fēng)險(xiǎn)綜合評(píng)價(jià)分區(qū)面積由高到低依次為：較低風(fēng)險(xiǎn)區(qū)、低風(fēng)險(xiǎn)區(qū)、中等風(fēng)險(xiǎn)區(qū)、較高風(fēng)險(xiǎn)區(qū)、高風(fēng)險(xiǎn)區(qū)，依次為22 280.94，23 217.2，9 979.05，7 681.28，3 241.53 km2。其中，高風(fēng)險(xiǎn)區(qū)面積占比不足5%，集中分布在西吉、彭陽(yáng)、隆德部分地區(qū)；較高風(fēng)險(xiǎn)區(qū)面積占比為11.57%，分布在南部山區(qū)及鹽池最南部；中等風(fēng)險(xiǎn)區(qū)面積占比15.03%，主要分布在海原、同心等地；低、較低風(fēng)險(xiǎn)分布范圍廣泛，共計(jì)占比為68.52%。

圖1 土壤侵蝕敏感、景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)、綜合生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)空間分布

2.2 寧夏風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果的空間差異

從縣域尺度上看，土壤侵蝕敏感性等級(jí)在各縣(區(qū))分布差異較大(圖2)。無(wú)敏感區(qū)所占面積的比重大于60%的地區(qū)有金鳳區(qū)、興慶區(qū)、利通區(qū)、紅寺堡區(qū)，面積占比依次為69.41%，62.78%，62.61%，61.30%；輕度敏感區(qū)面積占比最大的前5名依次為：鹽池縣、惠農(nóng)區(qū)、靈武市、青銅峽市、永寧縣，面積占比均在40%以上；中度敏感區(qū)面積占比前5名為大武口區(qū)、彭陽(yáng)縣、西吉縣、海原縣、原州區(qū)；高度敏感區(qū)面積占比前5名為原州區(qū)、西吉縣、海原縣、彭陽(yáng)縣、涇源縣；存在極敏感區(qū)的縣(區(qū))有隆德縣、涇源縣、彭陽(yáng)縣、原州區(qū)、西吉縣、海原縣、同心縣，其面積占比依次為70.93%，48.21%，26.69%，21.57%，19.67%，6.34%，1.28%。

圖2 寧夏各縣(區(qū))土壤侵蝕敏感性等級(jí)分布

從圖3各縣(區(qū))景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果看，高風(fēng)險(xiǎn)區(qū)僅存在于西吉縣、彭陽(yáng)縣、原州區(qū)、隆德縣、鹽池縣、海原縣地區(qū)，其中，西吉縣面積占比為88.39%，彭陽(yáng)縣面積占比為79.94%；較高風(fēng)險(xiǎn)區(qū)面積占比前5名為金鳳區(qū)、原州區(qū)、興慶區(qū)、賀蘭縣、永寧縣，面積比重均達(dá)到50%，其中，金鳳區(qū)較高風(fēng)險(xiǎn)區(qū)面積比重為96.89%；中等風(fēng)險(xiǎn)區(qū)面積占比前5名為平羅縣、西夏區(qū)、中寧縣、永寧縣、海原縣；低、較低風(fēng)險(xiǎn)區(qū)面積共計(jì)所占比重大于70%的地區(qū)有大武口區(qū)、涇源縣、靈武市、利通縣、青銅峽市、沙坡頭區(qū)，占比依次為92.48%，80.13%，78.24%，77.81%，75.32%，70.75%。

圖3 寧夏各縣(區(qū))景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)等級(jí)分布

從圖4綜合評(píng)價(jià)結(jié)果來(lái)看，高風(fēng)險(xiǎn)區(qū)僅存在于西吉縣、隆德縣、彭陽(yáng)縣、原州區(qū)、海原縣和鹽池縣，面積占比依次為36.19%，27.10%，25.98%，12.31%，1.82%，0.72%。較高風(fēng)險(xiǎn)區(qū)面積占比大于40%的地區(qū)有原州區(qū)、隆德縣、西吉縣。高、較高風(fēng)險(xiǎn)區(qū)共計(jì)面積占比達(dá)60%以上的有西吉縣、隆德縣、彭陽(yáng)縣、原州區(qū)，占比依次為76.68%，74.78%，63.89%，60.65%。中等風(fēng)險(xiǎn)區(qū)面積占比前3名為涇源縣、海原縣、彭陽(yáng)縣，面積占比分別為45.35%，42.38%，35.48%。低、較低風(fēng)險(xiǎn)區(qū)面積占比達(dá)90%的地區(qū)有大武口區(qū)、利通區(qū)、青銅峽市、靈武市、惠農(nóng)區(qū)、平羅縣、西夏區(qū)、沙坡頭區(qū)、中寧縣、紅寺堡寺區(qū)、永寧縣、興慶區(qū)、金鳳區(qū)。其中，利通區(qū)和大武口區(qū)低、較低風(fēng)險(xiǎn)區(qū)面積共計(jì)占比為100%，無(wú)更高風(fēng)險(xiǎn)等級(jí)；西吉縣、彭陽(yáng)縣、隆德縣、原州區(qū)則無(wú)低風(fēng)險(xiǎn)區(qū)，土壤侵蝕與景觀綜合生態(tài)風(fēng)險(xiǎn)較大。

圖4 寧夏各縣(區(qū))土壤侵蝕與景觀綜合生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)等級(jí)分布

2.3 土壤侵蝕景觀綜合生態(tài)風(fēng)險(xiǎn)與土地利用的關(guān)系

在耕地、林地、草地、水域、建設(shè)用地、未利用土地6個(gè)一級(jí)地類中，寧夏土地利用類型面積最大的為草地，面積為26 834.99 km2，面積占比為40.41%；面積第二大的是耕地，面積為19 447.55 km2，占比為29.29%；其次為林地，面積為11 843.59 km2，占比為17.84%；建設(shè)用地、水域、未利用土地面積共計(jì)占比為12.46%，面積依次為5 081.64，1 926.43，1 265.78 km2。

將2019年土地利用空間分布圖和綜合生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果進(jìn)行空間疊置分析，獲得研究區(qū)不同土地利用類型的綜合生態(tài)風(fēng)險(xiǎn)程度信息(圖5)。由圖5可知，不同土地利用類型土壤侵蝕與景觀的綜合生態(tài)風(fēng)險(xiǎn)面積占比差異較大。耕地、草地、建設(shè)用地中各風(fēng)險(xiǎn)區(qū)隨等級(jí)升高面積占比減小，并且主要以低、較低風(fēng)險(xiǎn)為主，面積占比在50%以上。林地和水域這兩種土地利用類型中除較低風(fēng)險(xiǎn)區(qū)面積占比最大外也符合上述規(guī)律，但是未利用土地與之相反，其面積占比最大的風(fēng)險(xiǎn)區(qū)為高風(fēng)險(xiǎn)區(qū)，占比達(dá)63.16%。

圖5 寧夏不同土地利用類型的綜合生態(tài)風(fēng)險(xiǎn)分布

3 討 論

由人為干擾引起的土地利用變化往往造成嚴(yán)重的土壤侵蝕，分析不同土地利用類型與土壤侵蝕的關(guān)系意義重大[10]。對(duì)寧夏而言，未利用土地地區(qū)土壤侵蝕綜合評(píng)價(jià)的風(fēng)險(xiǎn)等級(jí)更高，這與何君等[25]認(rèn)為未利用地的土壤侵蝕敏感性最強(qiáng)結(jié)論一致。而鄒琴英等[12]的研究則認(rèn)為土壤侵蝕主要發(fā)生在耕地、有林地、草地，這與不同研究區(qū)內(nèi)各類土地利用類型的面積占比有關(guān)，耕地、林地、草地在研究區(qū)內(nèi)面積占比大，發(fā)生土壤侵蝕的絕對(duì)面積大，但是對(duì)同一土地利用來(lái)說(shuō)，耕地、林地、草地中未發(fā)生土壤侵蝕的面積更大。

前人的研究往往通過(guò)比較單因子與綜合因子敏感性評(píng)價(jià)結(jié)果，來(lái)確定土壤侵蝕的主要的影響因子，但土壤侵蝕是個(gè)復(fù)雜的過(guò)程，不同地區(qū)差異顯著。例如，張恩偉等[7]在云南滇池流域和莫建飛等[26]在廣西巖溶區(qū)的研究中都認(rèn)為降雨侵蝕力因子是影響土壤侵蝕敏感性的最關(guān)鍵因子；而冉濤等[27]在北方生態(tài)脆弱區(qū)的研究中則認(rèn)為土壤侵蝕力因子的敏感性被充分識(shí)別，降雨侵蝕力因子受站點(diǎn)數(shù)量和空間分布的影響作用被削弱。根據(jù)研究區(qū)自身特點(diǎn)建立不同的評(píng)價(jià)模型，成為今后研究趨勢(shì)，莫建飛等[26]在研究中增加石漠化等級(jí)變化為影響因子構(gòu)建了巖溶區(qū)水土流失敏感性評(píng)價(jià)模型。楊青林等[11]對(duì)川西北地區(qū)的土壤侵蝕敏感性評(píng)價(jià)研究時(shí)增加了溝壑密度為評(píng)價(jià)因子。結(jié)合寧夏干旱多風(fēng)的特點(diǎn)增加“風(fēng)”這一影響因子進(jìn)行評(píng)價(jià)發(fā)現(xiàn)，增加風(fēng)力因子后，無(wú)敏感、輕度、中度等級(jí)區(qū)面積分別減少1 661.87，2 115.83，388.72 km2，高、極高敏感區(qū)面積增加2 508.77，1 657.65 km2。目前的研究水力侵蝕與風(fēng)力侵蝕是分開計(jì)算的，沒有考慮水力侵蝕、風(fēng)力侵蝕同時(shí)發(fā)生時(shí)的交互作用。寧夏存在大面積的水風(fēng)蝕交錯(cuò)區(qū)，在今后研究中應(yīng)對(duì)比總結(jié)風(fēng)水交錯(cuò)的交互關(guān)系，同時(shí)增加研究的時(shí)間尺度，從而獲得更加精準(zhǔn)的結(jié)果，為水土流失的綜合防治提供理論支撐。

4 結(jié) 論

(1) 從數(shù)值上看，在寧夏土壤侵蝕與景觀生態(tài)風(fēng)險(xiǎn)綜合評(píng)價(jià)中，高風(fēng)險(xiǎn)區(qū)面積占比不足5%，較高風(fēng)險(xiǎn)區(qū)面積占比為11.57%，中等風(fēng)險(xiǎn)區(qū)面積占比15.03%，低、較低風(fēng)險(xiǎn)分布范圍廣泛，共計(jì)占比為68.52%。

(2) 從空間分布上看，高風(fēng)險(xiǎn)區(qū)集中分布在西吉、彭陽(yáng)、隆德部分地區(qū)；較高風(fēng)險(xiǎn)區(qū)分布在南部山區(qū)及鹽池最南部；中等風(fēng)險(xiǎn)區(qū)主要分布在海原、同心等地，較低等級(jí)以下風(fēng)險(xiǎn)區(qū)在各縣(區(qū))均有分布。

(3) 在不同的土地利用類型中，各風(fēng)險(xiǎn)等級(jí)占比不同。耕地、草地、建設(shè)用地中低風(fēng)險(xiǎn)區(qū)面積占比最大，林地、水域中較低風(fēng)險(xiǎn)區(qū)占比最大，未利用土地中高風(fēng)險(xiǎn)區(qū)占比最大。