晉江流域生態(tài)環(huán)境脆弱性評(píng)價(jià)

李劍鋒,朱玉晨,劉春雷,李亞松,郝奇琛,李政紅

(1.中國(guó)地質(zhì)科學(xué)院水文地質(zhì)環(huán)境地質(zhì)研究所,河北 石家莊 050061;2.自然資源部地下水科學(xué)與工程重點(diǎn)實(shí)驗(yàn)室,河北 石家莊 050061)

近年來,我國(guó)經(jīng)濟(jì)經(jīng)歷了飛速的發(fā)展,在發(fā)展經(jīng)濟(jì)的同時(shí),引發(fā)的水土流失、環(huán)境污染、植被退化等一系列生態(tài)問題威脅著人類賴以生存的環(huán)境[1]。在這種環(huán)境下,學(xué)者從多角度開展了生態(tài)環(huán)境脆弱性研究。從研究尺度上看,省、市以及重點(diǎn)地區(qū)、城市群均有學(xué)者開展研究:陸海燕[2]、任源鑫等[3]以省級(jí)單位為評(píng)價(jià)單元,分別開展了新疆維吾爾族自治區(qū)、陜西省的生態(tài)脆弱性研究;王貝貝等[4]以市為研究單元,對(duì)南昌市開展脆弱性評(píng)價(jià)研究;張慧琳等[5]、周玲美等[6]以地區(qū)作為評(píng)價(jià)單元,分別對(duì)山西五臺(tái)山地區(qū)、內(nèi)蒙古呼日查干淖爾地區(qū)開展脆弱性評(píng)價(jià)研究。從評(píng)價(jià)體系看,多采用SRP概念模型[2-5]來搭建評(píng)價(jià)體系,主成分分析法[4]、綜合復(fù)權(quán)法[3]、層次分析法[6]等也有學(xué)者采用。POLSKY C等[7]受美國(guó)公共空間計(jì)劃整合框架的啟示,發(fā)展了基于“暴露-敏感-適應(yīng)”的VSD(Vulnerability Scoping Diagram)評(píng)價(jià)整合模型。

晉江流域是典型的南方丘陵地區(qū),多山地、丘陵。素有“中國(guó)茶葉之都”之稱的安溪縣就位于晉江流域的上游。然而,在茶葉產(chǎn)業(yè)發(fā)展的過程中,由于土地利用不當(dāng),茶園過度開發(fā),出現(xiàn)了大面積紅壤嚴(yán)重退化、水土流失、滑坡侵蝕等諸多環(huán)境問題[8]。本文依托地質(zhì)調(diào)查項(xiàng)目“廈漳泉同城化地區(qū)綜合地質(zhì)調(diào)查”,以福建省晉江流域作為評(píng)價(jià)區(qū),基于“暴露性-易損性-恢復(fù)力”三個(gè)維度開展評(píng)價(jià),應(yīng)用Getis-Ord Gi*統(tǒng)計(jì)模型、Calinski-Harabasz指數(shù)評(píng)價(jià)聚類模型等評(píng)價(jià)方法進(jìn)行分析,為流域尺度上開展脆弱性評(píng)價(jià)提供了新思路。

1 研究區(qū)概況

晉江流域主要位于福建省泉州市(圖1),地勢(shì)西北高、東南低，既有低山丘陵,又有高峰峽谷,河流盆地穿插其中。有東溪和西溪兩大支流,于南安市雙溪口匯入晉江干流,水系具有格子狀扇形水系的特征,流域面積為5 600 km2。晉江流域地處低緯度地區(qū),東臨海洋,屬于中亞熱帶海洋性季風(fēng)氣候,冬短夏長(zhǎng),光熱豐富,降雨充沛,降雨量年際間變化率大,季風(fēng)氣候顯著,干、濕季分明。土壤主要類型以南園屬紅壤、河田屬紅壤、魁岐屬紅壤為主,其他為老紅壤、赤紅壤、黃紅壤、潮土、鹽土、黃赤土等土壤類型。植被主要有亞熱帶季風(fēng)雨林和中亞熱帶常綠闊葉林,以次生林和人工林為主,杉木林、馬尾松林、建柏林等人工林也較為常見。

1.水系;2.公路;3.縣區(qū);4.流域范圍

2 評(píng)價(jià)指標(biāo)體系構(gòu)建

2.1 影響因素分析

生態(tài)系統(tǒng)脆弱性是由多種因素綜合作用引起的,從外界威脅的程度、自我保護(hù)的能力與受到破壞后的恢復(fù)力3方面考慮,可將脆弱性劃分為暴露性、易損性和恢復(fù)力[9-10]。

暴露性反映了生態(tài)系統(tǒng)遭受外界的影響和破壞程度,主要包含自然破壞與人為破壞。自然破壞方面,降雨量是引發(fā)水土流失等生態(tài)地質(zhì)問題的首要條件,因此降雨量越大,生態(tài)地質(zhì)問題發(fā)生的可能性就越大。人為破壞分為2類:一類是由于人類生活導(dǎo)致的生態(tài)系統(tǒng)破壞,從空間分布角度看,海拔200～500 m、起伏度50～200 m為丘陵地帶,該區(qū)域是人類活動(dòng)最頻繁的區(qū)域;另一類是人類生產(chǎn)行為對(duì)自然的破壞,晉江流域是茶葉的重要產(chǎn)地,由于不合理的開發(fā)導(dǎo)致茶園易發(fā)生生態(tài)地質(zhì)問題,茶園種植面積越大,該地區(qū)受到的生態(tài)風(fēng)險(xiǎn)也越大。從土地利用類型來看,建設(shè)用地是人類改造自然的區(qū)域,這兩類區(qū)域都是生態(tài)地質(zhì)問題易發(fā)的區(qū)域。

易損性是生態(tài)系統(tǒng)受到外界影響破壞時(shí)自身的抗干擾能力,取決于環(huán)境自身的屬性特征。坡度越陡的地區(qū)越容易發(fā)生地質(zhì)災(zāi)害,易損性越高。不同土壤類型對(duì)水土保持、植被發(fā)育、生態(tài)適宜性也有所區(qū)別,適宜性差的區(qū)域,受到外界脅迫時(shí)更易受損。土壤有機(jī)質(zhì)不同的區(qū)域土壤保肥性不同,有機(jī)質(zhì)越高的區(qū)域植被生長(zhǎng)越好,生態(tài)條件越穩(wěn)定，易損性越低。

恢復(fù)力指的是生態(tài)系統(tǒng)在遭受外界破壞后的自我恢復(fù)能力和調(diào)節(jié)能力,植被初級(jí)生產(chǎn)力是植物通過光合作用生產(chǎn)有機(jī)物質(zhì)的總量,是估算地球支持能力和評(píng)價(jià)陸地生態(tài)系統(tǒng)可持續(xù)發(fā)展的重要指標(biāo),植被初級(jí)生產(chǎn)力越高,生態(tài)系統(tǒng)恢復(fù)的可能性越大,恢復(fù)的也越快。

綜上所述,選擇降雨量、丘陵面積比、茶園種植面積比、建設(shè)用地面積比、坡度、土壤類型、土壤有機(jī)質(zhì)與植被初級(jí)生產(chǎn)力8個(gè)指標(biāo)建立脆弱性評(píng)價(jià)指標(biāo)體系。

2.2 評(píng)價(jià)指標(biāo)權(quán)重

晉江流域生態(tài)系統(tǒng)脆弱性評(píng)價(jià)指標(biāo)體系中各個(gè)指標(biāo)的權(quán)重采用層次分析法確定[11]。根據(jù)模型選定的評(píng)價(jià)指標(biāo),將評(píng)價(jià)指標(biāo)兩兩比對(duì),建立判斷矩陣(表1)。首先,對(duì)判斷矩陣進(jìn)行一致性檢驗(yàn),判斷矩陣是否具有一致性,通常采用計(jì)算CI及CR值的方法,當(dāng)CR≤0.1時(shí)具有可接受的一致性。經(jīng)計(jì)算λmax=8.05,CI=(λmax-10)/9=0.007 2,查詢8階判斷矩陣RI=1.41,得到CR=CI/RI=0.005 1,滿足CR值≤0.1,滿足一致性要求。由此得到各指標(biāo)圖層及綜合權(quán)重值(表2)。

表1 層次分析判斷矩陣

表2 晉江流域生態(tài)脆弱性評(píng)價(jià)指標(biāo)及權(quán)重

3 脆弱性評(píng)價(jià)

3.1 數(shù)據(jù)收集與處理

系統(tǒng)收集研究區(qū)用于脆弱性評(píng)價(jià)的數(shù)據(jù),包括多年平均降雨量、研究區(qū)高程、植被初級(jí)生產(chǎn)力、遙感影像、土壤類型、土壤有機(jī)質(zhì)等數(shù)據(jù)(表3)。對(duì)遙感影像數(shù)據(jù)進(jìn)行解譯,獲取了研究區(qū)茶園種植區(qū)分布圖與建設(shè)用地分布圖,進(jìn)而統(tǒng)計(jì)各類面積比值。

表3 數(shù)據(jù)來源

3.2 脆弱性評(píng)價(jià)分區(qū)

在晉江流域脆弱性指數(shù)的基礎(chǔ)上,將晉江流域生態(tài)環(huán)境按照生態(tài)環(huán)境脆弱性指數(shù)的高低分為脆弱性較高分區(qū)、脆弱性中等分區(qū)、脆弱性較低分區(qū)3個(gè)級(jí)別(圖2)。本文采用自然斷點(diǎn)分級(jí)法,該方法在考慮了數(shù)據(jù)聚類特征的基礎(chǔ)上,進(jìn)一步兼顧了數(shù)據(jù)的數(shù)量與范圍。

1.脆弱性較高分區(qū);2.脆弱性中等分區(qū);3.脆弱性較低分區(qū);4.河流;5.縣市

通過脆弱性分析可知,安溪縣感德鎮(zhèn)、長(zhǎng)坑鄉(xiāng)、蓬萊鎮(zhèn)所在的子流域脆弱性最高,南安市南部到泉州市的平原區(qū)脆弱性較低,其他地區(qū)脆弱性一般。整體來看,西溪比東溪脆弱,上游比下游脆弱。

安溪縣幾個(gè)鄉(xiāng)鎮(zhèn)脆弱性較高的主要原因是受茶園種植以及土壤性質(zhì)的影響。這幾個(gè)鄉(xiāng)鎮(zhèn)的主要土壤類型為河田屬紅壤,土層中富含白色石英砂,茶園墾殖過程中常年鏟除雜草,雨后細(xì)土隨地表徑流流失,土壤逐漸砂化。由于茶樹的聚鋁性特征[17],導(dǎo)致茶園表土迅速酸化,進(jìn)而加劇土壤表層肥力的退化。茶園的開采對(duì)區(qū)域的原生植被有破壞作用,植被覆蓋信息與初級(jí)生產(chǎn)力指標(biāo)均表明這幾個(gè)鄉(xiāng)鎮(zhèn)恢復(fù)力一般,特別是感德鎮(zhèn)與長(zhǎng)坑鄉(xiāng),初級(jí)生產(chǎn)力與植被覆蓋指標(biāo)均較差,更加凸顯了脆弱性。

4 生態(tài)環(huán)境脆弱性評(píng)價(jià)

4.1 基于多維偽F統(tǒng)計(jì)量動(dòng)態(tài)K-均值聚類分析算法的脆弱性主導(dǎo)因素識(shí)別

基于上述評(píng)價(jià)模型可知,影響晉江流域生態(tài)環(huán)境脆弱性的二級(jí)指標(biāo)共8種,為了分析影響研究區(qū)脆弱性在不同子流域的特點(diǎn),對(duì)這些指標(biāo)進(jìn)行聚類分析[18],以便深入了解研究區(qū)不同位置的脆弱性主導(dǎo)因素。

研究采用基于多維偽F統(tǒng)計(jì)量動(dòng)態(tài)K-均值聚類分析算法。K-均值聚類分析算法是一種廣泛應(yīng)用于數(shù)據(jù)聚類分析的迭代變換算法,通過指定類別數(shù)K,基于給定的聚類目標(biāo)函數(shù),采用迭代更新的方法,使得最終聚類結(jié)果目標(biāo)函數(shù)值為極小值,達(dá)到較優(yōu)的聚類效果[19-20]。類別數(shù)K通過Calinski-Harabasz指數(shù)評(píng)價(jià)聚類模型[21]多維偽F統(tǒng)計(jì)量來確定(表4)。

表4 偽F統(tǒng)計(jì)量匯總表

對(duì)聚類模型來說,通常聚類結(jié)果為同一類的數(shù)據(jù)距離越近越好,而不同類別的數(shù)據(jù)距離越遠(yuǎn)越好,偽F統(tǒng)計(jì)量反映的是聚類間方差和聚類內(nèi)方差的比率,即反映組內(nèi)相似性和組間差異的比率。該值越大,代表組內(nèi)越緊密,不同組之間的差異越大,即有更優(yōu)的聚類結(jié)果。二級(jí)指標(biāo)的偽F統(tǒng)計(jì)值結(jié)果顯示,隨著類別數(shù)量增加,偽F統(tǒng)計(jì)值呈逐漸降低的趨勢(shì),即將研究目標(biāo)分為3類時(shí),分類效果最好。由此選定K值為3,進(jìn)行K-均值聚類分析,得到聚類分布圖(圖3)。

1.組1;2.組2;3.組3;4.河流;5.縣市

組1為研究區(qū)中部大部分地區(qū),面積3 862 km2,占研究區(qū)面積的68.88 %;組2為東南部城市聚集區(qū),面積915 km2,占研究區(qū)面積的16.33 %;組3為研究區(qū)西部劍斗鎮(zhèn)、長(zhǎng)坑鄉(xiāng)、西坪鎮(zhèn)等茶葉種植區(qū),面積829 km2,占研究區(qū)面積的14.79 %。創(chuàng)建基于組匯總的變量統(tǒng)計(jì)值見表5。

基于統(tǒng)計(jì)表反映出,組2的主導(dǎo)因素前三名分別為建設(shè)用地面積比、初級(jí)生產(chǎn)力與土壤有機(jī)質(zhì),坡度與降雨量的影響很小,充分反映了組2作為城市聚集區(qū)的特點(diǎn):建設(shè)用地面積比非常高,同時(shí)由于城市化對(duì)于土地的破壞,初級(jí)生產(chǎn)力與土壤有機(jī)質(zhì)均非常低,在城市聚集區(qū),脆弱性主要體現(xiàn)為居民活動(dòng)所造成的影響。因此,改善城市聚集區(qū)的脆弱性應(yīng)當(dāng)從以下方面入手:減緩?fù)恋氐恼加?、植樹造林增加植被覆蓋以及土壤修復(fù)。組3的主導(dǎo)因素為茶園,相對(duì)于其他因素非常突出,其他因素則相差不大,組3所反映的區(qū)域恰為研究區(qū)內(nèi)茶葉種植密集的區(qū)域,說明茶葉的過度種植對(duì)生態(tài)環(huán)境脆弱性的影響非常大,需要合理規(guī)劃,適度科學(xué)種植,避免過度種植引發(fā)生態(tài)破壞與災(zāi)害[22]。組1為研究區(qū)的中部大部分地區(qū),從統(tǒng)計(jì)指標(biāo)來看,各項(xiàng)要素較為均衡,差別較小,沒有出現(xiàn)特別突出的主導(dǎo)因素,說明在研究區(qū)的大部分區(qū)域,各項(xiàng)二級(jí)指標(biāo)對(duì)脆弱性共同產(chǎn)生作用,在各方面均需要適度注意。

4.2 基于Getis-Ord Gi*統(tǒng)計(jì)模型顯著性熱點(diǎn)分析的脆弱性聚集評(píng)價(jià)

應(yīng)用Getis-Ord Gi*統(tǒng)計(jì)模型對(duì)研究區(qū)評(píng)價(jià)結(jié)果進(jìn)行統(tǒng)計(jì)顯著性分析,開展晉江流域脆弱性較高區(qū)域與較低區(qū)域的聚集特征評(píng)價(jià)(圖4)。Getis-Ord Gi*是一種用于識(shí)別具有統(tǒng)計(jì)顯著性的高值(熱點(diǎn))和低值(冷點(diǎn))的空間聚類統(tǒng)計(jì)方法,它通過自動(dòng)聚合數(shù)據(jù),識(shí)別適當(dāng)?shù)姆治龇秶?并糾正多重測(cè)試和空間依賴性,對(duì)數(shù)據(jù)進(jìn)行查詢,以確定可優(yōu)化熱點(diǎn)分析結(jié)果的設(shè)置[23]。對(duì)本研究來說,生態(tài)環(huán)境脆弱性指數(shù)作為評(píng)價(jià)數(shù)據(jù),評(píng)價(jià)模型結(jié)果的高值區(qū),即熱點(diǎn)區(qū)域,代表了脆弱性相對(duì)較高的聚集區(qū)域。低值區(qū)，即冷點(diǎn)區(qū)域,代表了脆弱性相對(duì)較低的聚集區(qū)域。

1.熱點(diǎn)區(qū)間(99%置信);2.熱點(diǎn)區(qū)間(95%置信);3.熱點(diǎn)區(qū)間(90%置信);4.無明顯特征區(qū)間;5.冷點(diǎn)區(qū)間(90%置信);6.冷點(diǎn)區(qū)間(95%置信);7.冷點(diǎn)區(qū)間(99%置信)

晉江流域統(tǒng)計(jì)顯著性的聚類分析結(jié)果表明,脆弱性高值聚集區(qū)分布在感德鎮(zhèn)—蓬萊鎮(zhèn)一帶以及永春縣湖洋鎮(zhèn)一帶。其中,感德鎮(zhèn)是最高的脆弱性高值聚集區(qū),該區(qū)域即為晉江流域脆弱性最為聚集的區(qū)域。

5 茶園種植優(yōu)化建議

綜上所述,本文對(duì)晉江流域茶園種植分布現(xiàn)狀合理性進(jìn)行了評(píng)估。其中,生態(tài)地質(zhì)脆弱性高、坡度高、植被初級(jí)生產(chǎn)力、土壤有機(jī)質(zhì)含量與高程較低的茶園種植區(qū),不宜進(jìn)行茶葉種植,建議作為茶葉種植調(diào)整優(yōu)化區(qū)。位于生態(tài)地質(zhì)脆弱性中等的茶園種植區(qū),對(duì)生態(tài)地質(zhì)系統(tǒng)影響不突出,建議作為維持現(xiàn)狀區(qū)。而對(duì)于晉江流域中生態(tài)地質(zhì)脆弱性低、坡度較低、植被初級(jí)生產(chǎn)力與土壤有機(jī)質(zhì)含量較高的區(qū)域,生態(tài)恢復(fù)力高,脆弱性低,具有較好的種植條件,作為適度擴(kuò)大種植區(qū)。根據(jù)以上原則,分析得到晉江流域茶園種植建議面積統(tǒng)計(jì)表(表6)和優(yōu)化建議圖(圖5)。

表6 茶園種植建議面積統(tǒng)計(jì)表

建議茶園種植調(diào)整優(yōu)化區(qū)面積53.18 km2,主要分布在感德鎮(zhèn)、金谷鎮(zhèn)、長(zhǎng)坑鄉(xiāng)、金谷鎮(zhèn)四個(gè)鄉(xiāng)鎮(zhèn)與劍斗鎮(zhèn)南部;維持現(xiàn)狀區(qū)面積168.00 km2,主要位于桃舟鄉(xiāng)、劍斗鎮(zhèn)北側(cè)、玉斗鎮(zhèn)、西坪鎮(zhèn)、虎邱鎮(zhèn)等地區(qū);適度擴(kuò)大種植區(qū)主要分布在一都鎮(zhèn)、下洋鎮(zhèn)、桂陽(yáng)鎮(zhèn)、以及藍(lán)田鄉(xiāng),總面積約為87.06 km2,面積較大,即使擴(kuò)大種植業(yè)仍需要避免過度集中。

6 結(jié)論

(1)晉江流域脆弱性分級(jí)結(jié)果顯示,脆弱性較高的區(qū)域?yàn)楦械骆?zhèn)、長(zhǎng)坑鄉(xiāng)、蓬萊鎮(zhèn),南安市南部至泉州市的平原區(qū)脆弱性較低,整體來看,西溪比東溪脆弱,上游比下游脆弱。

(2)晉江流域的脆弱性條件在不同的地區(qū)有著不同的主導(dǎo)要素,可分為3種類型:東南部城市聚集區(qū),影響脆弱性的主導(dǎo)因素分別為建筑用地、初級(jí)生產(chǎn)力與土壤有機(jī)質(zhì),脆弱性主要體現(xiàn)為高頻率居民活動(dòng)所造成的影響;在晉江流域中西部,影響脆弱性的主導(dǎo)因素為茶葉的過度種植;其他區(qū)域脆弱性主導(dǎo)要素不明顯,多種要素共同作用。

(3)為了改善晉江流域的脆弱性,建議在城市聚集區(qū)減緩?fù)恋氐恼加?植樹造林增加植被覆蓋以及土壤修復(fù);在茶葉種植密集區(qū),合理規(guī)劃,適度科學(xué)種植,建議感德鎮(zhèn)、金谷鎮(zhèn)、長(zhǎng)坑鄉(xiāng)、金谷鎮(zhèn)4個(gè)鄉(xiāng)鎮(zhèn)與劍斗鎮(zhèn)南部不宜種植茶葉,進(jìn)行茶葉種植調(diào)整優(yōu)化;在一都鎮(zhèn)、下洋鎮(zhèn)、桂陽(yáng)鎮(zhèn)、以及藍(lán)田鄉(xiāng)適度擴(kuò)展種植茶葉。