張世文,卜中原,沈 強(qiáng),楊邵文,胡青青,周 妍, 羅明

1 安徽理工大學(xué)地球與環(huán)境學(xué)院,淮南 232001 2 自然資源部土地整治重點(diǎn)實(shí)驗(yàn)室,北京 100035 3 安徽理工大學(xué)測(cè)繪學(xué)院, 淮南 232001 4 自然資源部國(guó)土整治中心,北京 100035

生態(tài)系統(tǒng)的穩(wěn)定不僅是社會(huì)經(jīng)濟(jì)可持續(xù)發(fā)展的前提,也是人類賴以生存的基礎(chǔ)。長(zhǎng)期以來(lái)自然資源的過(guò)度開(kāi)發(fā),對(duì)我國(guó)生態(tài)系統(tǒng)的平衡狀態(tài)造成了巨大破壞,部分關(guān)系國(guó)家發(fā)展前景的核心地區(qū)的生態(tài)系統(tǒng)已經(jīng)出現(xiàn)了不同程度的退化[1-7]?；洷蹦蠋X是國(guó)家生態(tài)文明建設(shè)和生態(tài)環(huán)境保護(hù)的重要規(guī)劃和試點(diǎn)區(qū),生態(tài)退化將影響“粵港澳大灣區(qū)”和“一帶一路”等國(guó)家區(qū)域戰(zhàn)略目標(biāo)的實(shí)現(xiàn)。山水林田湖草是一個(gè)生命共同體,具有相互聯(lián)系、相互依存、相互制約的特點(diǎn),單要素治理往往顧此失彼,不僅達(dá)不到很好效果,還會(huì)導(dǎo)致生態(tài)系統(tǒng)穩(wěn)定性的破壞。因此,生態(tài)系統(tǒng)的保護(hù)和修復(fù)需要綜合治理[8],需要將過(guò)去單一的要素保護(hù)轉(zhuǎn)變?yōu)橐远嘁貥?gòu)成的生命系統(tǒng)共治共管、統(tǒng)一保護(hù)和修復(fù)。

圖1 研究區(qū)位置圖Fig.1 Location of study area

國(guó)家高度重視對(duì)“山水林田湖草”的綜合治理和整體修復(fù)。2016年10月,由財(cái)政部、國(guó)土資源部、環(huán)境保護(hù)部聯(lián)合發(fā)表的《關(guān)于推進(jìn)山水林田湖生態(tài)保護(hù)修復(fù)工作的通知》中,明確要求全國(guó)各地開(kāi)展山水林田湖生態(tài)修復(fù),并在2017年8月,中央全面深化改革領(lǐng)導(dǎo)小組第三十七次會(huì)議上又將“草”納入,自此構(gòu)成“山水林田湖草”的生命共同體[9-11]。遵循物質(zhì)與能量的源匯理念,按流域尺度,揭示其內(nèi)部的物質(zhì)運(yùn)移與轉(zhuǎn)換以及對(duì)整個(gè)流域的影響,分析其修復(fù)阻力和優(yōu)先級(jí),可為修復(fù)項(xiàng)目在時(shí)空尺度上落地提供科學(xué)依據(jù),從而達(dá)到科學(xué)治理的目的。目前該領(lǐng)域的研究尚存一定的不足。研究方向上局限于風(fēng)險(xiǎn)源的歸類識(shí)別,較少考慮由于風(fēng)險(xiǎn)源所表現(xiàn)出的內(nèi)外空間分布差異,包括對(duì)于植被覆蓋、地形地貌、水質(zhì)環(huán)境等多種因素間的綜合作用[12-16]。研究手段上直接將小尺度區(qū)域的方法應(yīng)用到大尺度區(qū)域,很難達(dá)到滿意的效果,降低了研究方法的實(shí)踐價(jià)值。

綜合考慮當(dāng)前研究進(jìn)展和不足,以粵北南嶺大寶山礦流域?yàn)檠芯繀^(qū),基于遙感和GIS軟件,選取代表山水林田綜合質(zhì)量和治理可行性的地質(zhì)災(zāi)害、植被覆蓋狀況等8種指標(biāo)作為阻力面分級(jí)指標(biāo)。利用最小累積阻力模型(minimal cumulative resistance,MCR)構(gòu)建阻力面,并根據(jù)阻力大小進(jìn)行修復(fù)優(yōu)先級(jí)分區(qū),分區(qū)提出修復(fù)重點(diǎn)和對(duì)策。研究結(jié)果可以為阻力識(shí)別、項(xiàng)目立項(xiàng)優(yōu)選、山水林田湖草生態(tài)共治提供技術(shù)支持和理論依據(jù)。

1 研究區(qū)概況與數(shù)據(jù)獲取

1.1 研究區(qū)概況

大寶山礦流域?qū)儆趶V東粵北南嶺山區(qū)山水林田湖草生態(tài)保護(hù)修復(fù)試點(diǎn)中礦山及土壤生態(tài)修復(fù)的重點(diǎn)區(qū)域。大寶山礦位于研究區(qū)西北部,是以鐵、銅、硫、鉬為主的大型多金屬礦山,歷史上私挖亂采活動(dòng)頻繁,且老舊礦山較多,部分礦產(chǎn)資源開(kāi)發(fā)利用至今仍然保持著粗放開(kāi)發(fā)方式,技術(shù)設(shè)施落后,地表裸露嚴(yán)重,嚴(yán)重影響了礦區(qū)及其周邊的生態(tài)環(huán)境。研究區(qū)地形總體呈現(xiàn)北高南低,以丘陵為主,流經(jīng)礦區(qū)的水域流向大致由北向南,沿地勢(shì)從東南部流出(圖1)。

1.2 阻力評(píng)價(jià)指標(biāo)體系建立與數(shù)據(jù)獲取

基于山水林田湖草全要素系統(tǒng)理論,分析研究區(qū)生態(tài)問(wèn)題和礦業(yè)開(kāi)采情況,選取研究區(qū)水質(zhì)、地表植被、重金屬污染等主要的生態(tài)阻力因素,建立了涵蓋生態(tài)修復(fù)5個(gè)要素類型、5個(gè)阻力類型和8個(gè)指標(biāo)的生態(tài)阻力評(píng)價(jià)體系(表1)。

基于ArcGIS和ENVI軟件完成數(shù)據(jù)統(tǒng)計(jì)與分析。數(shù)據(jù)柵格大小全部為30 m,實(shí)地采樣調(diào)查于2017、2018年完成,測(cè)試方法按《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838—2002)和《土壤環(huán)境質(zhì)量 農(nóng)用地土壤污染風(fēng)險(xiǎn)管控標(biāo)準(zhǔn)》(GB 15618—2018)執(zhí)行。綜合遙感影像、礦區(qū)復(fù)墾、實(shí)地調(diào)查采樣等材料,獲取研究區(qū)用地類型分布、土地?fù)p毀程度、工程難易程度、地表水質(zhì)、土壤重金屬含量。綜合考慮空間預(yù)測(cè)精度和不確定性,基于經(jīng)驗(yàn)貝葉斯方法獲取土壤重金屬綜合污染指數(shù)和地表水綜合污染指數(shù)的空間分布格局。

2 研究方法

2.1 土壤可蝕性

土壤可蝕性是指土壤對(duì)侵蝕的敏感性,是反映地質(zhì)災(zāi)害程度的重要指標(biāo),是土壤對(duì)侵蝕抵抗力的倒數(shù),一般用K表示[17-19]。采用Williams等人建立的EPIC(erosion productivity impact calculator)公式進(jìn)行計(jì)算[20-21]。

(1)

式中,san、sil、cla和c分別代表土壤中砂粒、粉粒、黏粒、有機(jī)碳的含量(%),snl= 1-san/100。

表1 研究區(qū)生態(tài)阻力評(píng)價(jià)指標(biāo)

2.2 植被覆蓋度

為消除部分輻射誤差,利用多光譜遙感影像提取歸一化植被指數(shù)(normalized difference vegetation index, NDVI)[22-24]。

(2)

式中,NIR為紅外光譜反射值,R為紅光反射值。

利用歸一化植被指數(shù),采用像元二分模型估算植被覆蓋度(vegetation fraction coverage, VFC)[25-27]。

(3)

(4)

(5)

式中,NDVIsoil、NDVIveg分別表示裸地或無(wú)植被覆蓋和完全植被覆蓋的NDVI值,VFCmax、VFCmin分別表示在一定置信范圍內(nèi)的最大值與最小值。

2.3 土壤重金屬綜合污染指數(shù)

參考《土壤環(huán)境質(zhì)量 農(nóng)用地土壤污染風(fēng)險(xiǎn)管控標(biāo)準(zhǔn)》(GB 15618—2018)和《土壤環(huán)境監(jiān)測(cè)技術(shù)規(guī)范》(HJ/T 166—2006)[28-29],采用內(nèi)梅羅法計(jì)算土壤重金屬綜合污染指數(shù)。

(6)

式中,PImean和PImax分別是平均單項(xiàng)污染指數(shù)和最大單項(xiàng)污染指數(shù)。

2.4 地表水綜合污染指數(shù)

按照《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838—2002)[30],進(jìn)行各單項(xiàng)組分評(píng)價(jià),對(duì)各項(xiàng)指標(biāo)分別按地表水單項(xiàng)組分評(píng)分值賦值Fi,并帶入公式計(jì)算綜合評(píng)價(jià)分值F。單項(xiàng)組分評(píng)分值賦值分為五級(jí)(Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ),對(duì)應(yīng)的Fi值分別為0、1、3、6和10。

(7)

(8)

2.5 生態(tài)阻力模型的建立

在一定流域內(nèi)的生態(tài)環(huán)境問(wèn)題中,對(duì)區(qū)域生態(tài)環(huán)境安全造成威脅的景觀稱之為“風(fēng)險(xiǎn)源”。針對(duì)風(fēng)險(xiǎn)源,研究區(qū)的各指標(biāo)因素之間存在著不同程度的連接,多個(gè)連接共同產(chǎn)生影響,進(jìn)而產(chǎn)生生態(tài)風(fēng)險(xiǎn),而這些由連接產(chǎn)生的生態(tài)風(fēng)險(xiǎn)是需要通過(guò)克服阻力實(shí)現(xiàn)的。本文以流域內(nèi)的大寶山礦區(qū)為風(fēng)險(xiǎn)源。

最小累積阻力模型(MCR)是耗費(fèi)距離模型的衍生應(yīng)用,最初用來(lái)反映物種從源到目的地運(yùn)動(dòng)過(guò)程中所需耗費(fèi)的最小代價(jià)[31-32]。后來(lái)被廣泛應(yīng)用于生態(tài)領(lǐng)域,如物種保護(hù)、景觀格局分析等方面。該模型考慮源、空間距離和阻力基面3方面因素,表達(dá)式為：

(9)

式中,MCR為最小累積阻力；f是一未知負(fù)函數(shù),為最小累積阻力與生態(tài)適宜性的負(fù)相關(guān)關(guān)系；min為某景觀單元對(duì)不同源取累積阻力最小值；Dij為從源j到景觀單元i的空間距離；Ri為景觀單元i對(duì)運(yùn)動(dòng)過(guò)程的阻力系數(shù)。f函數(shù)是未知的,但Dij·Ri的累積值,被認(rèn)為是從源到空間某一點(diǎn)路徑的相對(duì)易達(dá)性的衡量,其中從所有源到某點(diǎn)阻力的最小值被用來(lái)衡量該點(diǎn)的易達(dá)性。采用Jenks自然斷點(diǎn)法,劃分阻力級(jí)別(表2)。

3 結(jié)果與分析

3.1 不同指標(biāo)阻力分級(jí)

基于ArcGIS和ENVI軟件,按照表2,構(gòu)建各指標(biāo)的阻力系數(shù)分級(jí)圖,并計(jì)算研究區(qū)綜合阻力系數(shù)值(圖2)。

由圖2可以看出,除地表水綜合污染指數(shù)外,其他評(píng)價(jià)指標(biāo)阻力系數(shù)空間分布格局具有一定的相似性,均表現(xiàn)出西北部高于東南部。研究區(qū)土壤可蝕性值變化范圍不大,位于0.25—0.32之間。土壤可蝕性阻力系數(shù)三級(jí)以下占總面積的77.32%,五級(jí)以上分布較少,僅占總面積的6.28%。研究區(qū)近75%的區(qū)域植被覆蓋度達(dá)到了70%以上,研究區(qū)西北部和東南部植被覆蓋度普遍低于50%。植被覆蓋度以一級(jí)為主,占總面積比為64.97%,三級(jí)和五級(jí)分別集中在研究區(qū)東南部和西北部,其占比率依次為7.92%和7.56%。礦區(qū)由于長(zhǎng)期開(kāi)采,造成礦區(qū)土壤重金屬嚴(yán)重超標(biāo),礦區(qū)土壤重金屬綜合污染指數(shù)阻力系數(shù)為五級(jí),占總面積的6.91%,周邊受其影響阻力系數(shù)為四級(jí),其余區(qū)域阻力系數(shù)均在三級(jí)以下,合占總面積的81.20%；研究區(qū)地勢(shì)北高南低,礦區(qū)排放的污染物經(jīng)雨水沖刷滲入土壤并流入河流,大部分地區(qū)土壤pH明顯低于韶關(guān)平均值(5.8),研究區(qū)西部土壤pH和南部地表水綜合污染指數(shù)較低,阻力系數(shù)均以四級(jí)和五級(jí)為主,其中五級(jí)分別占總面積的7.39%和32.52%。土地?fù)p毀程度四級(jí)主要集中在研究區(qū)西北部礦區(qū),占研究區(qū)總面積的9.61%,礦區(qū)周邊受到開(kāi)采和排污影響,損毀系數(shù)分別為二級(jí)和三級(jí),均占總面積的12.49%；從坡度阻力和工程難易程度來(lái)看,其一級(jí)分布區(qū)域較為相似,主要集中在研究區(qū)東南部和東北部,坡度阻力五級(jí)和工程難易程度三級(jí)分布較少,依次占總面積的9.49%和3.75%。研究區(qū)綜合阻力系數(shù)處于9—32之間,西北部高,東南部低。

表2 阻力分級(jí)

圖2 阻力系數(shù)分級(jí)圖Fig.2 Map of classification of resistance coefficient

3.2 不同評(píng)價(jià)指標(biāo)阻力面特征

利用最小累積阻力模型計(jì)算出單項(xiàng)指標(biāo)的生態(tài)阻力值。其中生態(tài)阻力值越高,區(qū)域生態(tài)環(huán)境的穩(wěn)定性越高,安全水平越高,生態(tài)修復(fù)阻力越小(圖3)。

圖3 生態(tài)阻力面Fig.3 Ecological resistance surface

由圖3可知,不同評(píng)價(jià)指標(biāo)阻力值均呈現(xiàn)出東南部大,西北部小的空間格局。研究區(qū)土壤可蝕性阻力值位于0—33589,礦區(qū)及周邊地區(qū)較低。土壤可蝕性反映了土壤受侵蝕的速度和能力,土壤可蝕性越低,生態(tài)阻力值越大。研究區(qū)地表水綜合污染指數(shù)阻力值位于0—40361,表現(xiàn)為上游阻力值低,下游阻力值高。地表水綜合污染指數(shù)阻力主要取決于高程和距離礦山的遠(yuǎn)近,水分匯聚到地勢(shì)較低的地區(qū),不容易產(chǎn)生分流,風(fēng)險(xiǎn)源的擴(kuò)散能力降低,阻力值增大。研究區(qū)植被覆蓋度阻力值位于0—13960,呈現(xiàn)礦區(qū)內(nèi)及其下游河流沿岸相對(duì)較低,其他區(qū)域較高的空間格局。植被覆蓋度反映了生態(tài)修改和質(zhì)量工作的質(zhì)量情況和環(huán)境的惡化程度,植被覆蓋度越高,生態(tài)阻力值越大,生態(tài)修復(fù)阻力越小。大寶山礦有近70年的開(kāi)創(chuàng)歷史,且位于研究區(qū)上游,大量污染物滲透到土壤、水體、大氣中并遷移到研究區(qū)中下游的其他地區(qū),造成了一系列重金屬污染和土壤酸化等問(wèn)題[33-34]。土壤重金屬綜合污染指數(shù)阻力值位于0—19049,pH阻力值位于0—23523,重金屬污染指數(shù)越低,阻力值越大,研究區(qū)pH明顯低于區(qū)域自然背景值[35],pH越高,阻力越大。研究區(qū)土地?fù)p毀程度阻力值位于0—14054,與礦業(yè)開(kāi)采區(qū)域的空間分布格局一致,常年采礦形成大面積的尾礦庫(kù)和排土場(chǎng),破壞了流域范圍內(nèi)的地質(zhì)環(huán)境,土地?fù)p毀程度越小,生態(tài)阻力越大,修復(fù)阻力越小。研究區(qū)工程難易程度阻力值位于0—13878,呈現(xiàn)礦區(qū)內(nèi)及上游河流沿岸相對(duì)較低,其他區(qū)域較高的空間格局,工程難易程度越低,生態(tài)阻力越大。研究區(qū)坡度阻力值位于0—15792,研究區(qū)主要以丘陵為主,地表起伏大,西北部坡度阻力值低,東南部坡度阻力值高。坡度越低,生態(tài)阻力值越大。

4 優(yōu)先級(jí)分區(qū)分析

圖4 修復(fù)優(yōu)先級(jí)分區(qū)Fig.4 Repair priority partition

根據(jù)景觀生態(tài)學(xué)中的源-匯理論,綜合分析由于礦業(yè)生產(chǎn)活動(dòng)造成生態(tài)環(huán)境問(wèn)題,疊加各單項(xiàng)生態(tài)阻力值,計(jì)算綜合阻力值。依據(jù)研究區(qū)不同區(qū)域格局優(yōu)化所針對(duì)的生態(tài)過(guò)程決定生態(tài)系統(tǒng)功能發(fā)揮的作用不同,并利用Jenks自然斷點(diǎn)法將研究區(qū)進(jìn)行修復(fù)實(shí)施的優(yōu)先級(jí)分區(qū)(圖4)。

整個(gè)流域共分成4個(gè)區(qū)。由圖4可知,Ⅰ—Ⅳ區(qū)分別占整個(gè)研究區(qū)面積的27.36%、33.63%、25.05%、13.96%。Ⅰ區(qū)位于在研究區(qū)的西北部礦區(qū),為整個(gè)流域的風(fēng)險(xiǎn)源,由于采礦活動(dòng)導(dǎo)致植被覆蓋度下降,地表破壞,土地?fù)p毀和污染嚴(yán)重。植被覆蓋度、重金屬含量、土地?fù)p毀是其主導(dǎo)阻力因素。Ⅱ—Ⅳ區(qū)位于流域中下游,依次遠(yuǎn)離大寶山礦區(qū)。礦山生產(chǎn),特別是早期的私挖亂采活動(dòng)導(dǎo)致大量污水排放和礦渣侵滲,污染了下游的水質(zhì)和土壤,造成了區(qū)域重金屬超標(biāo)、土壤酸化、土地破壞等問(wèn)題,致使研究區(qū)生態(tài)系統(tǒng)服務(wù)功能遭到破壞、生態(tài)壞境退化。Ⅱ—Ⅳ區(qū)的主要影響因素是重金屬含量和水質(zhì)。

阻力值越低,生態(tài)修復(fù)越具有緊迫性和優(yōu)先性。研究區(qū)修復(fù)優(yōu)先性依次為Ⅰ區(qū)>Ⅱ區(qū)>Ⅲ區(qū)>Ⅳ區(qū)。這也與廣東省山水林田湖草生態(tài)保護(hù)修復(fù)試點(diǎn)項(xiàng)目強(qiáng)調(diào)“源頭控制,末端治理”的理念相一致。廣東省山水林田湖草生態(tài)保護(hù)修復(fù)試點(diǎn)框架下的首選落地項(xiàng)目應(yīng)為Ⅰ區(qū)。

5 分區(qū)治理重點(diǎn)與對(duì)策

基于阻力面分析和修復(fù)優(yōu)先級(jí)分區(qū),分區(qū)提出治理重點(diǎn)與對(duì)策。

(1)Ⅰ區(qū)為礦區(qū)集中區(qū),也是流域的源頭,植被覆蓋度低、重金屬含量高、土地?fù)p毀嚴(yán)重是其主導(dǎo)的修復(fù)障礙,治水、降低重金屬含量和恢復(fù)植被是生態(tài)修復(fù)的核心內(nèi)容。①治水。采用清污分流措施(設(shè)置截洪溝、排洪隧洞、排洪豎井等)、攔水渠技術(shù),從根本上治理酸性與重金屬超標(biāo)礦水外排的問(wèn)題。②提高pH,降低重金屬有效態(tài)含量。采取生物或化學(xué)措施對(duì)礦區(qū)土地進(jìn)行基質(zhì)改良,如施加化學(xué)肥料(150—400 kg/hm2)、石灰(10—20 t/hm2)、有機(jī)肥(10—30 t/hm2)等材料,培肥地力,中和土壤酸度,降低土壤重金屬的生物有效性,以保證一個(gè)良好的植被重建效果。③恢復(fù)植被。選擇耐貧瘠、耐干旱的重金屬耐性植物種類：如泡桐、桉樹(shù)、苧麻、類蘆、高羊茅和狗牙根等,在群落結(jié)構(gòu)配置上以草灌植物為主,適當(dāng)客土移植耐性較強(qiáng)的喬木。喬木、灌木和草本的配置面積比例分別是10%—20%、20%—30%、50%—70%。為了提高植物的成活率,植物采用營(yíng)養(yǎng)袋繁殖,將營(yíng)養(yǎng)袋直接放置到植穴或種植帶中,草種采用撒播、條播或點(diǎn)播法,然后蓋一層稻草保證種子順利萌發(fā)。

(2)Ⅱ區(qū)核心制約因素為水質(zhì),阻力值在該區(qū)域較小,其次是酸堿度和土壤重金屬污染。①河流治理。建設(shè)礦區(qū)外排水處理廠,實(shí)現(xiàn)污染減排的目的；繼續(xù)加大實(shí)施攔泥庫(kù)清淤工程力度,增強(qiáng)攔泥庫(kù)的廢水調(diào)節(jié)能力,同時(shí)加大河道清理和河堤工程建設(shè)等。②治土。根據(jù)污染程度對(duì)礦區(qū)周邊耕地進(jìn)行分類污染治理。對(duì)于污染嚴(yán)重的農(nóng)田,宜采用植物穩(wěn)定技術(shù)。如能源作物麻瘋樹(shù)+土壤改良(0.5%石灰石和2%粉煤灰)組合。對(duì)于污染程度中等的農(nóng)田土壤,適宜采用植物提取技術(shù)進(jìn)行修復(fù),即通過(guò)種植重金屬超富集植物(東南景天、籽粒莧),將土壤中的重金屬含量逐步降低?？刹捎脰|南景天+低累積玉米套種+土壤改良(0.2%石灰石)或刈割處理籽粒莧模式,來(lái)恢復(fù)土壤的農(nóng)業(yè)使用用途。

(3)Ⅲ、Ⅳ區(qū)為輕度污染農(nóng)田區(qū),核心在于治土。采用植物阻隔技術(shù)進(jìn)行修復(fù)。低積累玉米、低積累水稻、豆角、花生等作物均適合在適當(dāng)?shù)母牧蓟A(chǔ)上(輕度污染旱地為白云石,輕度污染水田為粉煤灰),作為植物阻隔技術(shù)的實(shí)施材料,在逐步改善土壤條件的同時(shí),又可為農(nóng)民帶來(lái)經(jīng)濟(jì)收益。

6 討論與結(jié)論

6.1 討論

山水林田湖草生態(tài)保護(hù)修復(fù)強(qiáng)調(diào)綜合治理和整體修復(fù),隨著試點(diǎn)工作逐步開(kāi)展,試點(diǎn)過(guò)程中科學(xué)研究已成為山水林田湖草生態(tài)保護(hù)修復(fù)工作的瓶頸,項(xiàng)目立項(xiàng)和不同階段實(shí)施片區(qū)的確定具有很大盲目性,導(dǎo)致區(qū)域核心生態(tài)問(wèn)題把脈不準(zhǔn),項(xiàng)目落地時(shí)空不合理,治理缺乏針對(duì)性等問(wèn)題。本文結(jié)合相關(guān)研究成果,提出了大尺度區(qū)域優(yōu)先級(jí)分區(qū)治理的模式,具有一定的理論和現(xiàn)實(shí)意義,可以為修復(fù)項(xiàng)目在時(shí)空尺度上落地和科學(xué)治理提供科學(xué)依據(jù)。

早在2001年,陳述彭院士就提出了“統(tǒng)一規(guī)劃,分類指導(dǎo),分區(qū)推進(jìn),嚴(yán)格監(jiān)督”的生態(tài)環(huán)境的指導(dǎo)思想和基本原則[36]。為改變單要素治理的弊端,更好的助力于山水林田湖草生態(tài)環(huán)境治理,需要考慮由于風(fēng)險(xiǎn)源所表現(xiàn)出的內(nèi)外空間分布差異,針對(duì)不同區(qū)域分析多種因素間的綜合作用。本文結(jié)合研究區(qū)自然條件和礦業(yè)開(kāi)采的歷史遺留問(wèn)題,系統(tǒng)分析了由于風(fēng)險(xiǎn)源所引起的一系列生態(tài)修復(fù)阻力因素。評(píng)價(jià)指標(biāo)體系符合吳濃娣和余新曉等人提出的全面分析、系統(tǒng)治理的思想體系,遵循山水林田湖草全要素系統(tǒng)理論的要求[37-38]。所提出的分區(qū)治理思路與王英等人的研究成果相似,進(jìn)一步論證了優(yōu)先級(jí)分區(qū)治理模式的合理性[39]。制定的分區(qū)治理重點(diǎn)與對(duì)策,符合研究區(qū)實(shí)際情況,與《廣東粵北南嶺山區(qū)山水林田湖草生態(tài)保護(hù)修復(fù)試點(diǎn)實(shí)施方案》相契合。山水林田湖草生態(tài)環(huán)境修復(fù)是一項(xiàng)長(zhǎng)期工程,需要全方位的考慮,因此在后續(xù)的研究過(guò)程中應(yīng)針對(duì)不同優(yōu)先級(jí)區(qū)域進(jìn)行更為深入的專題研究,以便更加科學(xué)地制定出修復(fù)方案。

6.2 結(jié)論

以粵北南嶺大寶山礦流域?yàn)檠芯繀^(qū),采用最小累積阻力模型,構(gòu)建區(qū)域生態(tài)阻力面,根據(jù)生態(tài)阻力對(duì)研究區(qū)進(jìn)行修復(fù)優(yōu)先級(jí)分析,并提出分區(qū)治理對(duì)策。

(1)不同評(píng)價(jià)指標(biāo)阻力值均呈現(xiàn)出東南部大,西北部小的空間格局。長(zhǎng)期的礦業(yè)開(kāi)采活動(dòng)污染了礦區(qū)流域水質(zhì)和土壤,造成了區(qū)域重金屬富集、土壤酸化、土地破壞等問(wèn)題,并且距離礦區(qū)越近的區(qū)域生態(tài)修復(fù)的優(yōu)先性越高。

(2)阻力值越低,生態(tài)修復(fù)越具有緊迫性和優(yōu)先性。研究區(qū)共分為4個(gè)優(yōu)先級(jí)區(qū),研究區(qū)修復(fù)優(yōu)先性依次為Ⅰ區(qū)>Ⅱ區(qū)>Ⅲ區(qū)>Ⅳ區(qū)。

(3)Ⅰ區(qū)重點(diǎn)在于治水,降低pH,降低重金屬有效態(tài)含量,提高植被覆蓋度,主要措施包括清污分流、土壤改良和植被恢復(fù)；Ⅱ區(qū)重點(diǎn)在于河流治理與治土,主要措施包括建立污水處理廠、植物穩(wěn)定和植物提??；Ⅲ區(qū)、Ⅳ區(qū)核心在治土,主要措施為植物阻隔。