呂田田, 陳艷梅, 宋 婷, 丁疆輝, 張 揚(yáng), 劉銘華

(1.河北師范大學(xué) 資源與環(huán)境科學(xué)學(xué)院, 河北省環(huán)境變化遙感識(shí)別技術(shù)創(chuàng)新中心, 河北省環(huán)境演變與生態(tài)建設(shè)實(shí)驗(yàn)室,石家莊 050024; 2.中國環(huán)境科學(xué)研究院, 國家環(huán)境保護(hù)區(qū)域生態(tài)過程與功能評(píng)估重點(diǎn)實(shí)驗(yàn)室, 北京 100012)

2017年4月1日，中共中央、國務(wù)院決定在河北省雄縣、容城、安新三縣及周邊部分區(qū)域設(shè)立雄安新區(qū)，用以轉(zhuǎn)移北京的非首都職能[1-2]，雄安新區(qū)所在區(qū)域生態(tài)環(huán)境保護(hù)與修復(fù)工作引起各界重視。雄安新區(qū)生態(tài)支撐區(qū)為大清河流域，是全國水土流失嚴(yán)重區(qū)域之一[3]，研究該區(qū)域水土保持功能時(shí)空格局對(duì)于有針對(duì)性進(jìn)行區(qū)域生態(tài)保護(hù)、修復(fù)以及生態(tài)建設(shè)工作具有重要現(xiàn)實(shí)意義。

水土保持功能是生態(tài)系統(tǒng)通過其結(jié)構(gòu)與過程相互作用所具有的減少土壤水蝕作用的能力，即單位面積潛在土壤侵蝕量與實(shí)際土壤侵蝕量的差值，受到自然地理要素及人類生態(tài)建設(shè)活動(dòng)、農(nóng)業(yè)生產(chǎn)活動(dòng)、采礦業(yè)擾動(dòng)等多種因子影響，是生態(tài)系統(tǒng)提供的重要功能之一[4-5]。學(xué)者們一直進(jìn)行著水土保持算法的研究，并取得了豐碩的成果，其中，使用最成熟的就是通用水土流失模型(USLE)[6-7]。1965年在美國農(nóng)業(yè)部正式開始使用該模型[8]；1978年W.H.Wischmeier和D.Smith針對(duì)原有模型在實(shí)際應(yīng)用中存在的問題進(jìn)行了修正，形成RUSLE模型[9]；2009年Terranova等[10]采用RUSLE模型定量分析了意大利南部Calabria地區(qū)的水土流失情況，為當(dāng)?shù)卣耙?guī)劃部門治理水土流失提供了科學(xué)的參考和依據(jù)；Park等[11]采用RUSLE模型從流域視角對(duì)韓國1985—2005年的土壤侵蝕狀況進(jìn)行評(píng)估，發(fā)現(xiàn)其年均土壤侵蝕模數(shù)逐漸增加。Khalid chadli等[12]通過優(yōu)化RUSLE模型，對(duì)摩洛哥塞布流域的土壤侵蝕量進(jìn)行評(píng)估，其結(jié)果在該區(qū)域土地規(guī)劃方面發(fā)揮了重要作用。國外學(xué)者經(jīng)過多年的研究與改進(jìn)，該模型已經(jīng)成為評(píng)估精度較高、實(shí)用性較強(qiáng)的水土保持功能定量研究模型之一。在國內(nèi)，學(xué)者們也不斷對(duì)該模型進(jìn)行修正和改進(jìn)，王萬中、畢小剛、寧婷等[13-15]對(duì)降水侵蝕因子進(jìn)行了修正和改進(jìn)；高江波等[16]利用該模型評(píng)價(jià)了京津冀農(nóng)業(yè)生態(tài)系統(tǒng)的水土保持功能及價(jià)值，王嬌等[17]利用此模型并添加溝壑密度因子對(duì)河北太行山區(qū)水土流失進(jìn)行評(píng)價(jià)，揭示了不同地理背景下的水土流失敏感區(qū)的空間分布特征。以上研究不斷完善RUSLE模型中各因子取值與算法，使該方法在國內(nèi)應(yīng)用越來越廣泛。

該研究在ArcGIS 10.4軟件中采用RUSLE模型，對(duì)雄安新區(qū)生態(tài)支撐區(qū)2000—2018年水土保持功能時(shí)空格局進(jìn)行評(píng)價(jià)與分析，研究結(jié)果為確定雄安新區(qū)生態(tài)支撐區(qū)生態(tài)保護(hù)、建設(shè)與修復(fù)重點(diǎn)提供科學(xué)依據(jù)。

1 研究區(qū)概況

本研究確定的雄安新區(qū)生態(tài)支撐區(qū)為新區(qū)所在大清河流域全域，見圖1，對(duì)雄安新區(qū)生態(tài)環(huán)境具有重要影響。該流域地處113°39′—117°47′E，38°01′—40°03′N，面積約4.74萬km2。全區(qū)包括太行山區(qū)及大清河沖積平原，地形以山地、丘陵和平原為主，地勢西部高，東部低。氣候?qū)儆谂瘻貛О霛駶櫞箨懠撅L(fēng)性氣候，年降水量在400～650 mm，降水集中于夏季。區(qū)域地帶性植被為溫帶落葉闊葉林，太行山區(qū)植被覆蓋度高，平原地區(qū)以農(nóng)田、綠地等人工植被為主。土壤主要為褐土、潮土、粗骨土、栗褐土和棕壤等。2018年底，研究區(qū)內(nèi)常住人口約2 944.05萬人。

圖1 雄安新區(qū)生態(tài)支撐區(qū)位置與范圍

2 數(shù)據(jù)來源與方法

2.1 數(shù)據(jù)來源

本研究使用的數(shù)據(jù)包括氣象數(shù)據(jù)、土壤類型數(shù)據(jù)、遙感數(shù)據(jù)、DEM數(shù)據(jù)以及土地利用類型數(shù)據(jù)等。其中，氣象數(shù)據(jù)來源于中國氣象數(shù)據(jù)網(wǎng)(http:∥data.cma.cn/)的中國地面氣候標(biāo)準(zhǔn)值數(shù)據(jù)集，主要收集2000—2018年北京市、天津市、河北省和山西省氣象站點(diǎn)數(shù)據(jù)，包括年降水量、月降水量等，并利用ArcGIS 10.4進(jìn)行反距離權(quán)重插值分析。土壤類型數(shù)據(jù)，來源于世界土壤數(shù)據(jù)庫(HWDS)，該數(shù)據(jù)庫由聯(lián)合國糧農(nóng)組織(FAO)和維也納國際應(yīng)用系統(tǒng)研究所(ⅡASA)構(gòu)建。遙感NDVI數(shù)據(jù)來源于MODIS網(wǎng)站，通過ENVI 4.7軟件進(jìn)行校正和格式轉(zhuǎn)化，在ArcGIS 10.4中進(jìn)行合成、裁切和重采樣，并通過柵格統(tǒng)計(jì)得到。DEM高程數(shù)據(jù)為SRTM地形數(shù)據(jù)，空間分辨率為90 m，來源于中國科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心國際科學(xué)數(shù)據(jù)鏡像網(wǎng)站。土地利用類型數(shù)據(jù)為ArcInfo coverage格式，比例尺為1∶10萬，來源于中國科學(xué)院南京地理與湖泊研究所湖泊—流域數(shù)據(jù)集成與模擬中心。

2.2 研究方法

2.2.1 評(píng)估方法 修正通用水土流失模型(RUSLE)考慮了降水、土壤的可蝕性、地表植被覆蓋情況、地形坡度坡長和水土保持措施五大因子[18-19]。其表達(dá)式為：

AP=R×K×LS×C×P

(1)

Ar=R×K×LS

(2)

A=AP-Ar

(3)

式中:Ap為潛在土壤侵蝕量;Ar為現(xiàn)實(shí)土壤侵蝕量;A為水土保持量;R為降雨侵蝕力因子;K為土壤可侵蝕因子;LS為地形因子;C為植被覆蓋因子;P為水土保持措施因子。

降雨侵蝕力因子R，參考王萬中等人關(guān)于R值的計(jì)算[13]，并結(jié)合京津冀晉的實(shí)際情況對(duì)參數(shù)進(jìn)行調(diào)整，得到各氣象站點(diǎn)的R值，并通過ArcGIS 10.4掩膜提取雄安新區(qū)生態(tài)支撐區(qū)的R值。該R值能夠反映降水量與降水集中程度[20]。土壤可侵蝕因子K，采用威斯奇邁爾等人EPIC模型的K值估算方法[21-22]。地形因子LS，利用黃炎和[23]建立的LS計(jì)算方程，并結(jié)合京津冀晉數(shù)字高程模型(DEM)，通過ArcGIS 10.4掩膜提取得到。植被覆蓋因子C，利用馬超飛等[24]人關(guān)于植被覆蓋因子的計(jì)算公式，并結(jié)合京津冀晉植被覆蓋度提取得到。水土保持措施因子P，根據(jù)研究區(qū)土地利用實(shí)際情況，參考張憲奎等人[25]的研究進(jìn)行賦值：林地賦值1，草地賦值1，農(nóng)田賦值0.4，水面和建設(shè)用地賦值0，未利用地賦值1。

2.2.2 時(shí)空格局分析方法 本研究使用時(shí)間序列數(shù)據(jù)變化趨勢分析方法[26]，來研究雄安新區(qū)生態(tài)支撐區(qū)每個(gè)像元水土保持量隨時(shí)間的變化趨勢。采用相關(guān)系數(shù)的統(tǒng)計(jì)檢驗(yàn)方法對(duì)每個(gè)像元水土保持變化趨勢的相關(guān)系數(shù)進(jìn)行顯著性檢驗(yàn)，得到相關(guān)系數(shù)R，對(duì)比相關(guān)系數(shù)顯著性檢驗(yàn)表，得到與皮爾森相關(guān)系數(shù)0.01，0.05相對(duì)應(yīng)的R值[27]，并利用ArcGIS 10.4進(jìn)行分級(jí)和重分類處理，得到研究區(qū)水土保持功能顯著性檢驗(yàn)結(jié)果。

3 結(jié)果與分析

3.1 水土保持功能空間格局分析

水土保持量多年均值大于240 t/(hm2·a)的區(qū)域分布在研究區(qū)西部和北部山區(qū)，呈零星狀分布，具體包括房山區(qū)、阜平縣、淶水縣西部，靈壽縣西北部，易縣西北部和南部；水土保持量多年均值在120～240 t/(hm2·a)的分布區(qū)域呈東北—西南走向、面狀分布在最高值區(qū)的周圍；60～120 t/(hm2·a)的區(qū)域呈面狀分布在研究區(qū)太行山及山前丘陵，具體包括房山區(qū)、阜平縣、靈壽縣東部，曲陽縣、行唐縣、靈丘縣全域；低于60 t/(hm2·a)的區(qū)域分布在研究區(qū)東部平原，見圖2。

圖2 2000-2018年雄安新區(qū)生態(tài)支撐區(qū)水土保持功能均值空間分布

3.2 水土保持功能時(shí)空格局變化

3.2.1 水土保持量年均值變化 2000—2018年，雄安新區(qū)生態(tài)支撐區(qū)水土保持量年均值呈顯著增加趨勢，見圖3。研究區(qū)2000—2018年水土保持量各年均值變化較大，介于39.78～135.62 t/(hm2·a)，總體趨勢顯著增強(qiáng)。其中，2018年水土保持量均值最大，約為135.62 t/(hm2·a)；2001年水土保持量均值最小，約為39.78 t/(hm2·a)，這與當(dāng)年降水量大小有直接關(guān)系。

圖3 水土保持量各年均值變化

3.2.2 像元尺度水土保持功能時(shí)空格局動(dòng)態(tài)變化 像元尺度上，2000—2018年水土保持功能呈增強(qiáng)趨勢的面積比例較大，見表1。19 a間研究區(qū)水土保持功能呈增強(qiáng)趨勢的區(qū)域占研究區(qū)的63.89%，其中通過顯著性檢驗(yàn)的占研究區(qū)28.75%，未能通過顯著性檢驗(yàn)的占研究區(qū)35.14%；呈退化趨勢的區(qū)域占研究區(qū)36.11%，其中通過顯著性檢驗(yàn)的占研究區(qū)8.2%，未能通過顯著性檢驗(yàn)的占研究區(qū)27.91%。

表1 2000-2018年水土保持功能變化顯著性檢驗(yàn)

研究區(qū)水土保持功能通過0.01，0.05顯著性檢驗(yàn)的增強(qiáng)區(qū)面積分別為8 656.58 km2，4 971.00 km2，主要分布在研究區(qū)太行山北部，涉及房山區(qū)、淶水縣、易縣、淶源縣、靈丘縣、阜平縣、順平縣、滿城區(qū)等區(qū)縣。其中，通過0.01顯著性檢驗(yàn)的增強(qiáng)區(qū)占轄區(qū)面積較大的是房山區(qū)、淶水縣、淶源縣、易縣等；通過0.05顯著性檢驗(yàn)的增強(qiáng)區(qū)占轄區(qū)面積比例較大的是阜平縣、易縣等，詳見表2。

表2 重點(diǎn)生態(tài)功能區(qū)水土保持功能增強(qiáng)情況

統(tǒng)計(jì)分析發(fā)現(xiàn)，房山區(qū)顯著性檢驗(yàn)增強(qiáng)區(qū)面積為1 490.21 km2，占房山區(qū)75.45%；顯著性檢驗(yàn)退化區(qū)面積為25.52 km2，占房山區(qū)1.29%；未能通過顯著性檢驗(yàn)的面積為459.25 km2，占房山區(qū)23.26%。房山區(qū)是研究區(qū)內(nèi)顯著性檢驗(yàn)增強(qiáng)區(qū)占比最大的區(qū)縣，是研究區(qū)內(nèi)水土保持功能增強(qiáng)典型區(qū)。

研究區(qū)水土保持功能通過0.01，0.05顯著性檢驗(yàn)的退化區(qū)域面積分別為2 148.59 km2，1 738.15 km2，主要分布在研究區(qū)內(nèi)太行山區(qū)西南部及東部平原，因平原區(qū)水土保持功能不突出，不作為重點(diǎn)分析區(qū)域，主要分析太行山區(qū)水土保持功能變化。在太行山區(qū)，顯著性檢驗(yàn)退化區(qū)占轄區(qū)面積比例較大的是行唐縣、曲陽縣、靈丘縣、唐縣、靈壽縣、阜平縣、淶源縣等，其中通過0.01顯著性檢驗(yàn)的退化區(qū)占轄區(qū)面積比例較大的是行唐縣、曲陽縣、靈丘縣等，通過0.05顯著性檢驗(yàn)的退化區(qū)占轄區(qū)面積比例較大的是曲陽縣、行唐縣等(表3)。

表3 重點(diǎn)生態(tài)功能區(qū)水土保持功能退化情況

統(tǒng)計(jì)分析發(fā)現(xiàn)，行唐縣顯著性檢驗(yàn)退化區(qū)面積為99.69 km2，占行唐縣10.35%；顯著性檢驗(yàn)增強(qiáng)區(qū)面積為159.67 km2，占行唐縣16.57%；未能通過顯著性檢驗(yàn)的面積為704.20 km2，占行唐縣73.08%。曲陽縣顯著性檢驗(yàn)退化區(qū)面積為110.48 km2，占曲陽縣10.31%；顯著性檢驗(yàn)增強(qiáng)區(qū)面積為219.80 km2，占曲陽縣20.51%；未能通過顯著性檢驗(yàn)的面積為633.28 km2，占曲陽縣65.72%。按照顯著性檢驗(yàn)退化區(qū)面積占比，行唐、曲陽是太行山區(qū)諸縣中的第一、二名，是研究區(qū)內(nèi)水土保持功能退化典型區(qū)。

3.3 水土保持功能時(shí)空格局變化人為驅(qū)動(dòng)因素分析

水土保持功能是生態(tài)系統(tǒng)重要功能之一，其功能強(qiáng)弱受到自然和人類活動(dòng)等諸多因素共同影響。降水、土壤等自然因素對(duì)水土保持功能影響巨大，但人類在短時(shí)間內(nèi)難以改變，本研究選擇人為驅(qū)動(dòng)因素進(jìn)行分析。根據(jù)雄安新區(qū)生態(tài)支撐區(qū)水土保持功能時(shí)空格局研究結(jié)果及實(shí)際調(diào)研，發(fā)現(xiàn)人為驅(qū)動(dòng)因素主要有生態(tài)建設(shè)活動(dòng)，農(nóng)業(yè)生產(chǎn)與建設(shè)用地?cái)U(kuò)張，采礦業(yè)發(fā)展等。

3.3.1 生態(tài)建設(shè)活動(dòng)提高植被覆蓋度，增強(qiáng)水土保持功能 利用土地利用轉(zhuǎn)移矩陣法，分析2000—2018年水土保持功能典型增強(qiáng)區(qū)(房山區(qū))土地利用類型構(gòu)成與變化特征，詳見表4。19 a間由耕地向林地、草地等生態(tài)用地轉(zhuǎn)化面積為65.89 km2，說明該區(qū)域退耕還林還草等措施增加生態(tài)空間的占比，增強(qiáng)了區(qū)域水土保持功能。另外，該區(qū)對(duì)已有的生態(tài)空間進(jìn)行了生態(tài)修復(fù)，據(jù)北京市水土保持公報(bào)[28]顯示，2000—2018年房山區(qū)累計(jì)完成水土保持生態(tài)修復(fù)面積達(dá)320.65 km2，生態(tài)建設(shè)與生態(tài)修復(fù)利于地表植被自然恢復(fù)，植被覆蓋度明顯提高，增強(qiáng)區(qū)域生態(tài)空間的水土保持功能。

表4 2000-2018年水土保持功能典型增強(qiáng)區(qū)(房山區(qū))土地利用轉(zhuǎn)移矩陣

3.3.2 農(nóng)業(yè)生產(chǎn)活動(dòng)和建設(shè)用地?cái)U(kuò)張，降低水土保持功能 利用土地利用轉(zhuǎn)移矩陣法，分析2000—2018年水土保持功能典型退化區(qū)(行唐和曲陽兩縣)土地利用類型構(gòu)成與變化特征，詳見表5。由于兩縣耕地占土地總面積比例較大，耕地水土保持功能退化面積也較多，為168.20 km2，占兩縣水土保持功能退化區(qū)面積的45.05%。耕地水土保持功能退化主要原因是該典型區(qū)位于低山丘陵和山前平原區(qū)，長期以來該區(qū)域農(nóng)業(yè)生產(chǎn)活動(dòng)不規(guī)范，甚至出現(xiàn)很多撂荒和人為破壞現(xiàn)象，水土保持功能有所降低。另外，水土保持功能退化與建設(shè)用地面積增加有關(guān)。截至2018年，建設(shè)用地面積為15.56 km2，自2000年來各類用地向建設(shè)用地轉(zhuǎn)化面積為4.96 km2，分析發(fā)現(xiàn)顯著性退化區(qū)多分布于建設(shè)用地以及向建設(shè)用地轉(zhuǎn)化的區(qū)域周圍，建設(shè)活動(dòng)使地表受到人類擾動(dòng)，導(dǎo)致其水土保持能力下降。

表5 2000-2018年水土保持功能典型退化區(qū)(行唐和曲陽兩縣)土地利用轉(zhuǎn)移矩陣

3.3.3 采礦業(yè)擾動(dòng)生態(tài)系統(tǒng)穩(wěn)定性，降低區(qū)域水土保持能力 采礦用地上的生產(chǎn)建設(shè)活動(dòng)，包括開挖土石方、采掘礦石、修建工程設(shè)施、堆置礦渣等，若缺乏植被保護(hù)，極易造成水土保持功能退化。2000年研究區(qū)采礦用地面積為56.10 km2，2018年為84.63 km2，19 a間采礦用地面積增加51.86%。采礦活動(dòng)擾動(dòng)生態(tài)系統(tǒng)穩(wěn)定性，降低區(qū)域生態(tài)空間的水土保持功能。研究發(fā)現(xiàn)，有采礦用地的縣級(jí)行政區(qū)內(nèi)，水土保持功能退化區(qū)占轄區(qū)面積比例與新增采礦用地面積呈顯著正相關(guān)關(guān)系(R=0.463，p<0.05)，新增采礦用地面積越多，退化區(qū)占轄區(qū)面積比例越大(圖4)，如曲陽縣采礦用地增加面積最多，面積為4.26 km2，該縣水土保持功能退化面積占轄區(qū)面積比例也較大，為10.31%。

圖4 水土保持功能退化區(qū)占轄區(qū)面積比例與轄區(qū)新增采礦用地面積相關(guān)關(guān)系

4 討 論

本研究對(duì)雄安新區(qū)生態(tài)支撐區(qū)2000—2018年水土保持功能進(jìn)行評(píng)測，發(fā)現(xiàn)該區(qū)水土保持功能整體呈增強(qiáng)趨勢。其他學(xué)者多是針對(duì)京津冀開展研究，如李素曉[29]和潘梅等[30]也有較為一致的研究結(jié)果，本研究主要針對(duì)大清河流域開展研究，研究不僅發(fā)現(xiàn)區(qū)域水土保持功能整體增強(qiáng)，且水土保持功能時(shí)空變化異質(zhì)性明顯，研究區(qū)太行山區(qū)北部水土保持功能明顯增強(qiáng)，研究區(qū)西南部太行山區(qū)部分區(qū)域和東部平原區(qū)呈現(xiàn)退化趨勢。同時(shí)，有學(xué)者針對(duì)太行山水土保持功能退化的人為驅(qū)動(dòng)因素進(jìn)行研究，提出人口增長、森林遭到破壞、過度放牧以及生產(chǎn)建設(shè)活動(dòng)不合理[31]是水土保持功能退化的原因，但大部分研究成果主要對(duì)水土保持功能退化因素進(jìn)行定性分析。本研究通過調(diào)研并結(jié)合實(shí)際，采用轉(zhuǎn)移矩陣法、相關(guān)分析法，針對(duì)水土保持功能變化典型區(qū)的人為驅(qū)動(dòng)因素進(jìn)行定量分析，證明生態(tài)建設(shè)活動(dòng)是區(qū)域水土保持功能增強(qiáng)的主要原因，而農(nóng)業(yè)生產(chǎn)活動(dòng)和建設(shè)用地?cái)U(kuò)張及采礦業(yè)擾動(dòng)會(huì)導(dǎo)致研究區(qū)內(nèi)水土保持功能退化。研究成果能夠直接服務(wù)于雄安新區(qū)生態(tài)安全格局構(gòu)建，有利于提升新區(qū)人居環(huán)境適宜性。

基于水土保持功能時(shí)空格局變化和人為驅(qū)動(dòng)因素分析結(jié)果，本研究提出以下三點(diǎn)建議：一是將顯著性檢驗(yàn)增強(qiáng)區(qū)作為重點(diǎn)生態(tài)保護(hù)區(qū)，堅(jiān)持保護(hù)優(yōu)先，自然恢復(fù)為主，人工措施為輔的方針進(jìn)行保護(hù)，對(duì)現(xiàn)有森林、草地、未利用土地等進(jìn)行維護(hù)、管理，適當(dāng)控制人類活動(dòng)，維持生態(tài)系統(tǒng)穩(wěn)定性。二是將顯著性檢驗(yàn)退化區(qū)作為重點(diǎn)生態(tài)修復(fù)區(qū)，以人工方式為主對(duì)該區(qū)域進(jìn)行重點(diǎn)修復(fù)，加強(qiáng)農(nóng)田水土保育，規(guī)范農(nóng)業(yè)生產(chǎn)活動(dòng)，對(duì)坡耕地進(jìn)行退耕還林還草，增加林草地面積；嚴(yán)格控制建設(shè)用地?zé)o序擴(kuò)張，禁止生產(chǎn)生活用地?cái)D占生態(tài)用地。三是嚴(yán)格管控生態(tài)支撐區(qū)的采礦用地，從規(guī)劃到閉礦全過程進(jìn)行礦區(qū)生態(tài)建設(shè)和生態(tài)修復(fù)，提高礦區(qū)植被覆蓋度，提升該區(qū)域水土保持功能。

5 結(jié) 論

(1)研究區(qū)多年平均水土保持量呈西高東低分布格局。西部太行山區(qū)多年均值為191.10 t/(hm2·a)，中部丘陵為119.25 t/(hm2·a)，東部平原為30.25 t/(hm2·a)。

(2)2000—2018年，研究區(qū)水土保持功能整體有所增強(qiáng)。研究區(qū)水土保持年均值呈顯著增加趨勢，介于39.78～135.62 t/(hm2·a)；像元尺度上，水土保持功能呈增強(qiáng)趨勢的區(qū)域占研究區(qū)的63.89%，其中通過顯著性檢驗(yàn)的有28.75%，面積為13 627.58 km2，主要分布在研究區(qū)太行山北部；呈退化趨勢的區(qū)域占研究區(qū)的36.11%，其中通過顯著性檢驗(yàn)的有8.20%，面積為3 886.74 km2，主要分布在研究區(qū)太行山西南部和東部平原。

(3)研究區(qū)水土保持功能時(shí)空格局變化人為驅(qū)動(dòng)因素主要包括：一是生態(tài)建設(shè)活動(dòng)提高植被覆蓋度，增強(qiáng)水土保持功能；二是農(nóng)業(yè)生產(chǎn)活動(dòng)不合理和建設(shè)用地?zé)o序擴(kuò)張，導(dǎo)致水土保持功能退化；三是采礦業(yè)擾動(dòng)生態(tài)系統(tǒng)穩(wěn)定性，降低區(qū)域水土保持能力。

基于上述分析，提出應(yīng)將水土保持功能顯著增強(qiáng)區(qū)作為重點(diǎn)生態(tài)保護(hù)區(qū)，水土保持功能顯著退化區(qū)作為重點(diǎn)生態(tài)修復(fù)區(qū)，嚴(yán)格管控生態(tài)支撐區(qū)的采礦用地，從規(guī)劃到閉礦進(jìn)行全過程嚴(yán)格管理。