黃良美,陳 蓓,田 艷,黃 寧,黎 寧

廣西壯族自治區(qū)環(huán)境監(jiān)測(cè)中心站, 南寧 530028

岸邊帶生態(tài)系統(tǒng)研究對(duì)飲用水資源可持續(xù)的合理開(kāi)發(fā)利用與水環(huán)境規(guī)范、精確管理具有重要意義[1]。國(guó)內(nèi)外飲用水源相關(guān)研究多見(jiàn)報(bào)道,如飲用水源的水質(zhì)狀況、污染成因、生態(tài)安全、健康風(fēng)險(xiǎn)、等級(jí)分類(lèi)及其演變規(guī)律和評(píng)價(jià)方法等[2- 7],對(duì)水源地水體特征污染物、有毒有害物質(zhì)、微生物和沉積物輻射[3,5,8- 13]也進(jìn)行了深入系統(tǒng)檢測(cè)分析,另外,基于水質(zhì)模型、GIS技術(shù)、多元分析方法、信息熵、同位素示蹤法對(duì)水質(zhì)進(jìn)行評(píng)價(jià)、診斷、預(yù)測(cè)和預(yù)報(bào)[2,6-10,14]的研究受到重點(diǎn)關(guān)注。受我國(guó)飲用水源水質(zhì)達(dá)標(biāo)化管理和環(huán)境質(zhì)量目標(biāo)改善的環(huán)境保護(hù)政策引導(dǎo),流域復(fù)合生態(tài)系統(tǒng)有機(jī)整合、主體功能區(qū)規(guī)劃、生態(tài)紅線劃分、國(guó)土空間格局重構(gòu)及其土地利用方式優(yōu)化、景觀格局驅(qū)動(dòng)、水質(zhì)多媒介關(guān)聯(lián)、生態(tài)供需水量平衡、水源涵養(yǎng)分區(qū)及過(guò)程管制等研究[15-18]正成為新的興趣熱點(diǎn)。研究對(duì)當(dāng)前飲用水安全問(wèn)題的解決提供了有益探索和實(shí)踐。

目前,對(duì)飲用水源涵養(yǎng)功能與景觀格局、水資源、水環(huán)境質(zhì)量、污染物排放間同時(shí)進(jìn)行綜合定量分析及其關(guān)系耦合的研究報(bào)道較少。而我國(guó)飲用水源存在水源單一、水質(zhì)超標(biāo)、清潔供水不足、生態(tài)破壞、標(biāo)準(zhǔn)化建設(shè)與規(guī)范化管理落后等問(wèn)題,水環(huán)境保護(hù)、水污染防治和水供需不平衡等壓力很大,建立健康持續(xù)的飲用水安全保障體系尤為迫切。飲用水安全保障體系的建立需要長(zhǎng)遠(yuǎn)的飲用水源功能規(guī)劃、嚴(yán)格的保護(hù)區(qū)劃定、精準(zhǔn)的水量水質(zhì)管理和科學(xué)的水源地結(jié)構(gòu)功能調(diào)整與優(yōu)化配置,從而為保障水源涵養(yǎng)功能和實(shí)現(xiàn)水質(zhì)達(dá)標(biāo)奠定扎實(shí)基礎(chǔ)。我國(guó)《飲用水水源保護(hù)區(qū)劃分技術(shù)規(guī)范》(HJ T338—2007)、《集中式飲用水水源地規(guī)范化建設(shè)環(huán)境保護(hù)技術(shù)要求》(HJ773—2015)和《水源涵養(yǎng)林建設(shè)規(guī)范》(GB/T26903—2011)實(shí)施對(duì)飲用水源水質(zhì)安全和水源涵養(yǎng)起到了積極的規(guī)范化約束和標(biāo)準(zhǔn)化參考作用,但在實(shí)踐過(guò)程中存在“過(guò)度保護(hù)”和“保護(hù)不足”的問(wèn)題；另外,在森林水源涵養(yǎng)功能的邊界界定和計(jì)量方法方面也有不足之處。解決這些問(wèn)題一方面需要全面掌握當(dāng)?shù)氐牡乩砦恢?、水文、氣象、地質(zhì)特征、水動(dòng)力特性、水域污染類(lèi)型、污染特征、污染源分布、排水區(qū)分布、水源地規(guī)模、水量需求等情況,另一方面需要精確計(jì)算水源涵養(yǎng)功能的植被類(lèi)型、持水量大小及其最小覆蓋面積,以耦合優(yōu)化飲用水源景觀結(jié)構(gòu)與涵養(yǎng)功能的數(shù)量關(guān)系。

南寧市雖然水系發(fā)達(dá),但飲用水資源利用單一性特點(diǎn)突出,嚴(yán)重依賴河流型水源,湖庫(kù)型水源的標(biāo)準(zhǔn)化建設(shè)和規(guī)范化管理還不夠完善,存在生態(tài)破壞及工業(yè)、交通、養(yǎng)殖、生活污染等問(wèn)題,歷史上發(fā)生過(guò)2起水庫(kù)環(huán)境污染事故。為有效管理飲用水資源利用、污染防治和生態(tài)修復(fù),保障水質(zhì)安全,開(kāi)展南寧市湖庫(kù)型飲用水源景觀格局和水源涵養(yǎng)功能關(guān)系耦合優(yōu)化的基礎(chǔ)性研究和調(diào)查極為必要和迫切。

1 研究對(duì)象及方法

1.1 研究對(duì)象

南寧市位于107°45′—108°51′E、22°13′—23°32′N(xiāo)之間,南北長(zhǎng)201 km,東西寬234 km,建成區(qū)土地面積約190 km2。南寧地勢(shì)東北高西南低,以丘陵為主,平均海拔81 m,最高處為496 m；地形呈以邕江河谷為中心的盆地形態(tài),盆地向東開(kāi)口,南、北、西三面均為山地圍繞；地貌分平地、低山、石山、丘陵、臺(tái)地5種類(lèi)型。全市共有大、中、小型水庫(kù)779座,其中：庫(kù)容1×108m3以上的大型水庫(kù)5座,1×107m3以上的中型水庫(kù)25座,小型水庫(kù)749座,總庫(kù)容38×108m3,有效庫(kù)容16.51×108m3,現(xiàn)狀供水能力23.34×108m3。引水工程3154處,總引水規(guī)模達(dá)113.52 m3/s,現(xiàn)狀供水能力7.16×108m3。提水工程3950處,總提水規(guī)模為155.33 m3/s,現(xiàn)狀供水能力11.35×108m3。本研究湖庫(kù)型飲用水源主要包括大王灘水庫(kù)、鳳亭水庫(kù)、屯六水庫(kù)、西云江水庫(kù)、龍?zhí)端畮?kù)、天雹水庫(kù)、東山水庫(kù)、峙村河水庫(kù)和老虎嶺水庫(kù)等9個(gè)在用或備用水源地,其空間分布范圍見(jiàn)圖1。

圖1 南寧市9個(gè)湖庫(kù)型飲用水源分布范圍Fig.1 Study area location of the 9 reservoirs in Nanning city

1.2 研究方法

在評(píng)估分析南寧市湖庫(kù)型飲用水源景觀格局和水源涵養(yǎng)功能的基礎(chǔ)上,通過(guò)建立湖庫(kù)型飲用水源取水量與水源涵養(yǎng)持水量的耦合關(guān)系,對(duì)湖庫(kù)岸邊帶不同林份水源涵養(yǎng)所需最小面積進(jìn)行分析和優(yōu)化,以探討水源地保護(hù)區(qū)、集雨區(qū)和緩沖區(qū)劃分的適宜性；應(yīng)用多元逐步回歸分析方法擬合水源涵養(yǎng)功能與景觀格局、水資源、水質(zhì)、污染物排放間的數(shù)量關(guān)系；用多維尺度分析技術(shù)診斷飲用水源景觀結(jié)構(gòu)、水源涵養(yǎng)功能及其水環(huán)境質(zhì)量的風(fēng)險(xiǎn)因子。

1.2.1研究區(qū)邊界劃分

根據(jù)《飲用水水源保護(hù)區(qū)劃分技術(shù)規(guī)范》(HJ T338—2007)及《南寧市市區(qū)飲用水水源保護(hù)區(qū)劃定方案》批復(fù)文件有關(guān)要求,基于2012年9月—12月美國(guó)Worldview和德國(guó)Rapideye遙感影像數(shù)據(jù),按照“水體—水位線-陸域”的自然復(fù)合生態(tài)系統(tǒng)和景觀空間格局理念,參考當(dāng)前岸邊帶、緩沖區(qū)、生態(tài)交錯(cuò)帶、生態(tài)屏障等概念[1,18],并考慮后續(xù)數(shù)據(jù)統(tǒng)計(jì)處理與比較上的便利性,以水位線為界,劃分500 m和3 km的2個(gè)緩沖區(qū)對(duì)南寧市9個(gè)湖庫(kù)型飲用水源的景觀格局特點(diǎn)及水源涵養(yǎng)功能進(jìn)行研究。其中水體統(tǒng)一包含于2個(gè)緩沖區(qū),以避免《飲用水水源保護(hù)區(qū)劃分技術(shù)規(guī)范》中水體的約束性分割和多重性劃定,從而更符合景觀結(jié)構(gòu)整體性觀點(diǎn)。

1.2.2數(shù)據(jù)源及其處理技術(shù)方法

(1)遙感影像預(yù)處理

Rapideye和Worldview遙感影像的平面坐標(biāo)系統(tǒng)采用1980西安坐標(biāo)系,橫坐標(biāo)加帶號(hào)；在處理軟件中統(tǒng)一規(guī)定橢球和基準(zhǔn)面標(biāo)準(zhǔn),橢球長(zhǎng)、短半軸長(zhǎng)分別定義為6378140.0 m、6356755.288 m。投影方式采用高斯-克呂格投影；分帶方式按3°分帶；中央經(jīng)線經(jīng)度為108°。高程系統(tǒng)采用1985國(guó)家高程基準(zhǔn)。

參考底圖主要采用源于Google地圖的高分辨率衛(wèi)星影像,投影坐標(biāo)統(tǒng)一轉(zhuǎn)換至高斯-克呂格投影坐標(biāo)系統(tǒng)。正射糾正基于ERDAS遙感圖像處理軟件,利用已有的參考底圖和外業(yè)調(diào)查GPS靜態(tài)實(shí)測(cè)高精度控制點(diǎn)作為控制資料,引入數(shù)字高程模型數(shù)據(jù)(DEM),對(duì)所獲取的高分辨率衛(wèi)星影像進(jìn)行正射糾正、配準(zhǔn)、融合等處理,建立覆蓋全部調(diào)查范圍的數(shù)字正射影像。當(dāng)影像與基礎(chǔ)底圖采樣間隔均≤1 m時(shí),糾正圖像相對(duì)誤差控制在2.0—8.0倍；當(dāng)至少有一類(lèi)>1 m時(shí),其相對(duì)誤差限差控制在2.0—4.0倍。

(2)遙感影像解譯及其野外核查

采集南寧市區(qū)土地利用現(xiàn)狀圖件,并將其掃描、配準(zhǔn)、矢量化建庫(kù),為生態(tài)遙感解譯及其調(diào)查提供輔助信息。生態(tài)遙感解譯分類(lèi)系統(tǒng)主要參考2001年國(guó)土資源部頒布的《全國(guó)土地分類(lèi)(試行)》分類(lèi)系統(tǒng)。通過(guò)建立南寧市土地利用分類(lèi)的主要遙感影像解譯標(biāo)志,設(shè)置解譯判讀圖班大小規(guī)則(Worldview：2 m×5 m,Repideye：20 m×20 m)和精度控制要求,采用ArcGIS軟件完成遙感圖像的內(nèi)業(yè)解譯。

在龍?zhí)端畮?kù)和大王灘水庫(kù)選取了24個(gè)點(diǎn)位進(jìn)行核查,一級(jí)分類(lèi)的判讀精度為91.7%,二級(jí)分類(lèi)的判讀精度相對(duì)較低,為83.3%。

(3)水質(zhì)監(jiān)測(cè)

由于僅天雹水庫(kù)、峙村河水庫(kù)、老虎嶺水庫(kù)、龍?zhí)端畮?kù)、大王灘水庫(kù)有歷史性監(jiān)測(cè)數(shù)據(jù),且其監(jiān)測(cè)項(xiàng)目?jī)H為常規(guī)24項(xiàng)水質(zhì)指標(biāo),為了更加全面地掌握9個(gè)水庫(kù)的水質(zhì)和作為備用飲用水源的條件,在9個(gè)飲用水源取水口上游100 m各設(shè)置1個(gè)水質(zhì)監(jiān)測(cè)斷面,于2013年4月對(duì)9個(gè)水庫(kù)開(kāi)展水質(zhì)全分析補(bǔ)充監(jiān)測(cè),監(jiān)測(cè)項(xiàng)目為《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838—2002)中表1、表2、表3的108項(xiàng)指標(biāo),及透明度、葉綠素,共110項(xiàng)。

(4)數(shù)據(jù)處理方法

主要對(duì)遙感解譯數(shù)據(jù)、資料收集與野外調(diào)查數(shù)據(jù)和水質(zhì)監(jiān)測(cè)數(shù)據(jù)進(jìn)行分類(lèi)、集成與融合等信息綜合與數(shù)理統(tǒng)計(jì)分析。

1)水源涵養(yǎng)功能指數(shù)和景觀格局指數(shù)計(jì)算

基于遙感解譯的斑塊面積和數(shù)量,根據(jù)《生態(tài)環(huán)境狀況評(píng)價(jià)技術(shù)規(guī)范》(HJ/T192—2014)水源涵養(yǎng)指數(shù)的計(jì)算模式測(cè)算9個(gè)湖庫(kù)型飲用水源地在500 m和3 km 2個(gè)緩沖帶內(nèi)的綜合水源涵養(yǎng)功能大小。將林地、草地、水體和未利用地歸為自然用地類(lèi)型,耕地、建設(shè)用地歸為人工用地類(lèi)型,計(jì)算2種緩沖區(qū)和保護(hù)區(qū)的有關(guān)自然/人工景觀格局指數(shù)。

水源涵養(yǎng)功能指數(shù)=Nc×(0.45×(0.1×河流面積 + 0.3×湖庫(kù)面積 + 0.6×沼澤灘涂面積)+ 0.35×(0.6×有林地面積 + 0.25×灌木林地面積 + 0.15×疏林地和其他林地面積)+ 0.20×(0.6×高蓋度草地面積 + 0.3×中蓋度草地面積 + 0.1×低蓋度草地面積))/總面積,Nc為歸一化系數(shù),Nc=526.7。

自然/人工比(斑塊面積)=(林地+草地+水體+未利用地)÷(耕地+建設(shè)用地)

自然/人工比(斑塊個(gè)數(shù))=(林地+草地+水體+未利用地)÷(耕地+建設(shè)用地)

自然/人工比(綜合)=自然/人工比(面積)+ 自然/人工比(斑塊)

斑塊密度=斑塊數(shù)量÷斑塊面積

斑塊比例(數(shù)量)=類(lèi)型斑塊數(shù)量÷總斑塊數(shù)量

斑塊比例(面積)=類(lèi)型斑塊面積÷總斑塊面積

2)取水量與水源涵養(yǎng)功能關(guān)系耦合：根據(jù)生態(tài)學(xué)群落樣方、生物多樣性、景觀破碎度的最小面積原理[14-16],假設(shè)南寧市湖庫(kù)型飲用水源水資源為滿足南寧市居民生產(chǎn)生活飲用水需求的取水量與森林生態(tài)系統(tǒng)對(duì)水源涵養(yǎng)的持水量是一種等式關(guān)系[17-18],根據(jù)廣西森林生態(tài)系統(tǒng)水源涵養(yǎng)功能的測(cè)算結(jié)果,和南寧市飲用水源規(guī)劃與劃定方案中確定的取水量,測(cè)算出每個(gè)水庫(kù)不同森林生態(tài)系統(tǒng)類(lèi)型結(jié)構(gòu)與功能滿足假設(shè)等式關(guān)系的最小面積。

根據(jù)《森林生態(tài)系統(tǒng)服務(wù)功能評(píng)估規(guī)范》(LY/T1721—2008),森林生態(tài)系統(tǒng)涵養(yǎng)水量公式：

G調(diào)= 10A(P-E-C)

式中：G調(diào)為林分涵養(yǎng)水量(m3/a)；P為林外降水量,E為林分蒸散量,C為地表快速?gòu)搅髁?mm/a)；A為林分面積(hm2)。

3)多元回歸分析：設(shè)水源涵養(yǎng)功能指數(shù)為因變量y,分層設(shè)景觀格局、水資源、水質(zhì)和污染物排放等4種類(lèi)型的影響因子為自變量(x1,x2…xn),共25個(gè),開(kāi)展多元逐步回歸分析,計(jì)算矩陣為9×26,計(jì)算軟件為IBM SPSS 20.0。其中,水質(zhì)自變量先用主成分分析法提取2個(gè)主成分,計(jì)算矩陣9×110,即9個(gè)湖庫(kù)《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB3838—2002)表1、表2、表3的108項(xiàng)指標(biāo)和透明度、葉綠素；另外,用《地表水環(huán)境質(zhì)量評(píng)價(jià)辦法》(試行—2012)計(jì)算出水質(zhì)綜合污染指數(shù)。

4)主成分分類(lèi)分析：為推斷飲用水源空間結(jié)構(gòu)與屬性功能耦合關(guān)系,找出內(nèi)在規(guī)律性特點(diǎn),基于矩陣9×26進(jìn)行主成分分類(lèi)分析,生成空間結(jié)構(gòu)與指標(biāo)屬性間的耦合模型及其成果圖件。主成分分類(lèi)分析方法本質(zhì)上是一種最優(yōu)尺度上的同質(zhì)性分析技術(shù)[3,8,12],一方面可以起到因子分析和多維尺度分析的降維效果,另一方面又能集成對(duì)應(yīng)分析探討品質(zhì)類(lèi)型多變量關(guān)系的特點(diǎn),實(shí)現(xiàn)多變量多品質(zhì)在維度與空間分布上的最優(yōu)擬合過(guò)程。

2 研究結(jié)果

2.1 景觀格局

(1)基質(zhì)特征。表1結(jié)果表明,湖庫(kù)景觀的林地、草地、水體、耕地、建設(shè)用地等要素齊全,湖庫(kù)岸邊帶主要是水體與林地,其次草地,沿岸的農(nóng)業(yè)用地占一定斑塊數(shù)量,但各湖庫(kù)略有差異。9個(gè)湖庫(kù)主要體現(xiàn)以林地為主要景觀基質(zhì)的特征,其面積比例42.6%—87.3%。其中鳳亭河水庫(kù)500 m緩沖區(qū)林地基質(zhì)面積比例達(dá)74.3%,水體面積比例19.9%,耕地面積比例4.1%,建設(shè)用地面積比例0.2%,生態(tài)系統(tǒng)呈現(xiàn)自然穩(wěn)定性較人工干擾性略高的特點(diǎn),兩者比例為95.8∶4.2。3 km緩沖區(qū)用地類(lèi)型較500 m有所提高,表明景觀多樣性有所上升,異質(zhì)性提高。龍?zhí)端畮?kù)在3 km緩沖帶是唯一呈現(xiàn)以耕地為主要基質(zhì)的水庫(kù),面積比例達(dá)52.5%,其自然穩(wěn)定性較人工干擾性低,兩者比例為45.8∶54.2。峙村河水庫(kù)500 m緩沖區(qū)是唯一一個(gè)呈全自然景觀格局的生態(tài)系統(tǒng),耕地和建設(shè)用地比例極低；該水庫(kù)為降低人工干擾性,提高水庫(kù)的水源涵養(yǎng)功能,庫(kù)區(qū)實(shí)施封閉管理和移民搬遷。

表1 南寧市9個(gè)湖庫(kù)2種緩沖區(qū)景觀要素面積百分比例比較

(2)斑塊特征。比較湖庫(kù)型500 m緩沖區(qū)與3 km緩沖區(qū)的斑塊數(shù)量(表2),多呈現(xiàn)水體占優(yōu)的情況。其中,鳳亭河水庫(kù)500 m緩沖區(qū)水體斑塊數(shù)占3 km緩沖區(qū)水體斑塊數(shù)的85.8%,同理,工業(yè)用地占80.7%,林地占70.7%,草地占43.8%,農(nóng)業(yè)種植用地占35.6%,居民用地占12.5%。兩者用地類(lèi)型的斑塊數(shù)量比例表明,500 m岸邊帶一方面林地、水體和草地被開(kāi)發(fā)占用,另一方面各用地類(lèi)型較多,存在多種森林結(jié)構(gòu)如闊葉、針葉和經(jīng)濟(jì)林,灌木和草本。大王灘水庫(kù)岸邊帶各種用地破碎化明顯,隨著往陸域地帶縱深,破碎化下降,這種景觀格局表明大王灘水庫(kù)岸邊帶人為活動(dòng)干擾較強(qiáng)烈。西云江水庫(kù)500 m向3 km陸域縱深后水體斑塊基本不再發(fā)生變化,林地變化的空間也較小,居民點(diǎn)分布也變稀,但草地和農(nóng)業(yè)種植用地斑塊數(shù)均有較大的增加。天雹水庫(kù)500 m緩沖帶未見(jiàn)有居民點(diǎn)分布,但往3 km縱深后,居地用地增加了11個(gè)；同理,工業(yè)用地由1個(gè)增加到3個(gè),由此表明天雹水庫(kù)景觀破碎化以工業(yè)用地最為明顯。龍?zhí)端畮?kù)500 m岸邊帶未見(jiàn)有居民點(diǎn)分布,往3 km縱深后,居地用地增加了12個(gè)。峙村河水庫(kù)500 m岸邊帶未見(jiàn)有居民點(diǎn)分布、農(nóng)業(yè)種植用地和工業(yè)用地,往3 km縱深后,居民用地和工業(yè)用地各自增加了5個(gè),農(nóng)業(yè)種植用地增加了33個(gè)。

表2 南寧市9個(gè)水庫(kù)2種緩沖區(qū)斑塊數(shù)量比較

2.2 水源涵養(yǎng)功能

由圖2可知,500 m緩沖帶景觀格局對(duì)水源涵養(yǎng)的情況是鳳亭水庫(kù)(100)>峙村河水庫(kù)(86.1)>屯六水庫(kù)(85.3)>西云江水庫(kù)(85.2)>老虎嶺水庫(kù)(78.4)>天雹水庫(kù)(75.1)>東山水庫(kù)(53.1)>大王灘水庫(kù)(47.4)>龍?zhí)端畮?kù)(38.9)。測(cè)算結(jié)果與實(shí)際岸邊帶調(diào)查情況基本一致,鳳亭、屯六和西云江3個(gè)水庫(kù)均蓄水量豐富,周邊環(huán)境人為干擾較小,因此,水源涵養(yǎng)功能明顯。大王灘水庫(kù)由于曾經(jīng)作為游覽景區(qū)開(kāi)發(fā),岸邊帶建設(shè)了較多的各類(lèi)旅游景觀設(shè)施,對(duì)水環(huán)境有負(fù)面影響。龍?zhí)端畮?kù)岸邊帶目前存在較多環(huán)境問(wèn)題,種植和養(yǎng)殖排放對(duì)水環(huán)境造成明顯水質(zhì)惡化。

9個(gè)水庫(kù)的3 km緩沖帶中鳳亭河水庫(kù)和屯六水庫(kù)水源涵養(yǎng)功能大小變化趨勢(shì)與500 m相近,但是峙村河水庫(kù)、老虎嶺水庫(kù)、天雹水庫(kù)和龍?zhí)端畮?kù)的水源涵養(yǎng)功能較500 m緩沖帶則下降較大,大王灘水庫(kù)在3 km緩沖區(qū)尺度上水源涵養(yǎng)功能略有升高。鳳亭、屯六2個(gè)水庫(kù)水源涵養(yǎng)指數(shù)基本沒(méi)什么變化,說(shuō)明這2個(gè)水庫(kù)在更大尺度范圍的生態(tài)環(huán)境變化不明顯,呈林地>耕地>建設(shè)用地的景觀格局。大王灘水庫(kù)3 km緩沖帶水源涵養(yǎng)功能還略有提升,說(shuō)明大王灘目前主要是岸邊帶人為活動(dòng)對(duì)水環(huán)境有較大影響。峙村河水庫(kù)、老虎嶺水庫(kù)、天雹水庫(kù)、龍?zhí)端畮?kù)與西云江水庫(kù)3 km緩沖帶水源涵養(yǎng)功能均下降較大,表明水庫(kù)受周邊環(huán)境影響較大,目前這些湖庫(kù)岸邊帶的500 m水源地保護(hù)帶并不足于抵抗周邊環(huán)境開(kāi)發(fā)帶來(lái)的壓力,因此,需對(duì)水庫(kù)水質(zhì)的面源污染做詳盡調(diào)查與排查,預(yù)防和控制面源上的影響。

圖2 南寧市9個(gè)湖庫(kù)2種緩沖區(qū)水源涵養(yǎng)功能指數(shù)比較Fig.2 The water conservation function index comparation between the two types of buffer area for the 9 reservoirs,Nanning City

2.3 飲用水源水質(zhì)

飲用水源水質(zhì)監(jiān)測(cè)結(jié)果分析表明,飲用水源生態(tài)環(huán)境好,景觀系統(tǒng)結(jié)構(gòu)穩(wěn)定,自然/人工百分比例高,水源涵養(yǎng)功能指數(shù)高,飲用水源水質(zhì)就較好。如鳳亭河水庫(kù)、屯六水庫(kù)、西云江水庫(kù)水源涵養(yǎng)指數(shù)排前3位,而水質(zhì)監(jiān)測(cè)結(jié)果均滿足Ⅲ類(lèi)水質(zhì)標(biāo)準(zhǔn),水質(zhì)良好。而人為干擾性強(qiáng)的水庫(kù)則出現(xiàn)水質(zhì)超標(biāo)的情況。如大王灘水庫(kù)的主要污染因子有溶解氧、石油類(lèi)、汞、總磷、硒、高錳酸指數(shù)、五日生化需要量、氨氮、化學(xué)需氧量、總氮,其中化學(xué)需氧量、氨氮污染程度較大；龍?zhí)端畮?kù)的主要污染因子為溶解氧、化學(xué)需氧量?？偟鳛閰⒃u(píng)指標(biāo)因子,污染程度較大,多次出現(xiàn)劣Ⅴ類(lèi)水質(zhì)。

2.4 水源涵養(yǎng)功能優(yōu)化及其影響因子關(guān)系耦合

2.4.1飲用水源保護(hù)區(qū)范圍與水源涵養(yǎng)功能關(guān)系

根據(jù)《南寧市飲用水水源保護(hù)區(qū)劃定方案》,湖庫(kù)型飲用水源保護(hù)區(qū)分為一級(jí)保護(hù)區(qū)和二級(jí)保護(hù)區(qū),研究中9個(gè)在用或備用水庫(kù)型水源的保護(hù)區(qū)面積、集雨區(qū)面積與500 m緩沖區(qū)面積、3 km緩沖區(qū)面積見(jiàn)表3。

鳳亭河水庫(kù)由于跨市域,在南寧市轄區(qū)的保護(hù)區(qū)面積僅14.68 km2,與緩沖區(qū)范圍比較差異很大。另外,大王灘水庫(kù)、老虎嶺水庫(kù)和峙村河水庫(kù)保護(hù)區(qū)邊界部分在500 m緩沖區(qū)以外,但差異不明顯。大王灘水庫(kù)、龍?zhí)端畮?kù)、西云江水庫(kù)和鳳亭河水庫(kù)集雨區(qū)面積較3 km緩沖區(qū)面積略偏大?？傮w而言,除鳳亭河水庫(kù)跨界保護(hù)例外,湖庫(kù)型飲用水源保護(hù)區(qū)邊界劃定范圍與2種緩沖帶范圍不存在明顯沖突,基本能包括在2種緩沖帶范圍內(nèi)。因此,2種緩沖區(qū)景觀格局和水源涵養(yǎng)功能測(cè)算與分析基本能涵蓋保護(hù)區(qū)邊界劃定的要求。

表3 南寧市9個(gè)湖庫(kù)型飲用水源的保護(hù)區(qū)、集雨區(qū)及緩沖區(qū)面積

2.4.2飲用水源涵養(yǎng)功能最小面積優(yōu)化

廣西森林生態(tài)系統(tǒng)水源涵養(yǎng)功能位于0.00—68.53×108m3/a,不同森林類(lèi)型水源涵養(yǎng)功能調(diào)節(jié)水量：軟闊林5359.35 t hm-2a-1,松樹(shù)林2778.86 t hm-2a-1,桉樹(shù)林3190.49 t hm-2a-1,杉木林3787.38 t hm-2a-1,石山灌木林1140.00 t hm-2a-1,灌木經(jīng)濟(jì)林2727.22 t hm-2a-1,櫟林4746.52 t hm-2a-1,喬木經(jīng)濟(jì)林2598.17 t hm-2a-1,竹林4009.08 t hm-2a-1,硬闊葉林4574.81 t hm-2a-1,土山灌木林2645.76 t hm-2a-1。各森林類(lèi)型所占比例依次是軟闊林38.4%,松樹(shù)林14.7%,桉樹(shù)林11.3%,杉木林10.5%,石山灌木林5.8%,灌木經(jīng)濟(jì)林4.4%,櫟林3.6%,喬木經(jīng)濟(jì)林3.2%,竹林3.0%,硬闊葉林3.0%,土山灌木林2.2%。

據(jù)此測(cè)算的9個(gè)湖庫(kù)每種森林類(lèi)型水源涵養(yǎng)所需的最小面積(表4)表明當(dāng)前劃定的南寧飲用水源保護(hù)區(qū)林地面積只有東山水庫(kù)水源涵養(yǎng)功能能滿足年度取水量需求間的關(guān)系；老虎嶺水庫(kù)需要將所有林地保證為軟闊葉林和櫟林的前提下,方可滿足取水量需求間的關(guān)系；峙村河水庫(kù)則需要林地均為軟闊葉林、櫟林、硬闊葉林、竹林或杉木林。其他6個(gè)水庫(kù)在現(xiàn)狀前提下均無(wú)法滿足水源涵養(yǎng)功能與取水量間的關(guān)系,缺口最大的是鳳亭河和屯六水庫(kù),這兩個(gè)水庫(kù)為未來(lái)發(fā)展戰(zhàn)備水庫(kù),規(guī)劃的取水量為8789.5—17579.5萬(wàn)m3/a,而保護(hù)區(qū)面積劃定并不明確,測(cè)算中以劃定方案中確定的水庫(kù)正常水位線以外徑向距離為200 m的陸域。其次缺口較大的水庫(kù)還有大王灘水庫(kù)、西云江水庫(kù)、天雹水庫(kù)和龍?zhí)端畮?kù)。500 m緩沖區(qū)與保護(hù)區(qū)面積情況基本相似。緩沖到3 km后,鳳亭河與屯六水庫(kù)仍然因?yàn)槿∷烤薮?無(wú)法滿足水源涵養(yǎng)功能與取水量間的關(guān)系；另外龍?zhí)端畮?kù)和天雹水庫(kù)也無(wú)法滿足水源涵養(yǎng)功能與取水量間的關(guān)系。

根據(jù)現(xiàn)場(chǎng)調(diào)查與資料收集,9個(gè)水庫(kù)以軟闊葉林、櫟林、硬闊葉林、竹林、杉林、桉林、松林、灌木經(jīng)濟(jì)林、土山灌木林、喬木經(jīng)濟(jì)林的平均水源涵養(yǎng)功能為參照系,大王灘需要105.51 km2的水源涵養(yǎng)區(qū),東山水庫(kù)需要3.16 km2,鳳亭河與屯六水庫(kù)均需要257.61 km2、老虎嶺水庫(kù)和峙村河水庫(kù)需要15.83 km2,龍?zhí)端畮?kù)和天雹水庫(kù)需要21.10 km2,西云江水庫(kù)需要31.65 km2?，F(xiàn)劃定飲用水源保護(hù)區(qū)范圍,東山水庫(kù)和峙村河水庫(kù)基本符合要求,大王灘和峙村河水庫(kù)需保護(hù)區(qū)均為林地,而不得存在其他商業(yè)類(lèi)建設(shè)用地和農(nóng)田種植情況下,方可保證水源涵養(yǎng)功能與取水量間的關(guān)系。鳳亭河水庫(kù)、屯六水庫(kù)和、老虎嶺水庫(kù)、天雹水庫(kù)、龍?zhí)端畮?kù)和西云江水庫(kù)均無(wú)法滿足要求。

根據(jù)解譯數(shù)據(jù)結(jié)果及現(xiàn)場(chǎng)調(diào)查咨詢情況,部分水庫(kù)周邊林地多以桉樹(shù)林為主,其中天雹水庫(kù)、峙村河水庫(kù)、東山水庫(kù)、西云江水庫(kù)、老虎嶺水庫(kù)桉樹(shù)林面積占50%—75%。因此以桉樹(shù)類(lèi)型測(cè)算的最小水源涵養(yǎng)面積為參照,僅有東山水庫(kù)能滿足取水量與水源涵養(yǎng)功能間的關(guān)系。

表4 南寧市9個(gè)湖庫(kù)型飲用水源涵養(yǎng)水源最小森林面積測(cè)算

2.4.3飲用水源涵養(yǎng)功能影響因子多元關(guān)系耦合

多元回歸分析表明,水源涵養(yǎng)功能指數(shù)與水資源量、自然/人工的景觀面積比、景觀斑塊比和綜合比、水質(zhì)主成分因子1、水質(zhì)主成分因子2、水質(zhì)綜合污染指數(shù)、污水排放量、化學(xué)需氧量、氨氮、總磷、總氮和生活垃圾等指標(biāo)可建立起若干個(gè)逐步回歸方程式(表5)。

表5 飲用水源涵養(yǎng)功能及其影響因子的多元回歸分析

y水源涵養(yǎng)功能指數(shù) Water conservation function index；景觀格局x1—x3：x1自然/人工斑塊面積比 Area ratio of the natural and artifical patches、x2綜合比 Total ratio、x3斑塊個(gè)數(shù)比 Patches No. ratio of the natural and artifical patches；水資源x1—x4：x1水體面積 Water area、x2水位 Water level、x3集雨區(qū)面積 Catchment area、x4庫(kù)容量Reserviors capcity；水質(zhì)x1—x3：x1綜合污染指數(shù) Comprehensive pollutant index、x2主成分因子2 Principal component factor 2、x3主成分因子1 Principal component factor 1；污染物排放x1—x5：x1總氮 Total nitrogen、x2垃圾 Refuse、x3總磷 Total phosphorus、x4廢水 Waste water、x5污染源距離 Distance to the pollutant sources

水源涵養(yǎng)功能與水資源量間,影響水源涵養(yǎng)指數(shù)較大的因子有水體面積和水位,而這兩個(gè)因子的乘積就是水資源量；集雨區(qū)面積對(duì)水源涵養(yǎng)有弱正相關(guān),而庫(kù)容量則對(duì)水源涵養(yǎng)呈弱負(fù)相關(guān)。水源涵養(yǎng)功能與水質(zhì)間,最大污染指數(shù)影響較大,主成分1次之,而主成分因子1載荷量較大的水質(zhì)指標(biāo)有高錳酸鹽指數(shù)、氨氮、總磷、陰離子表面活性劑、化學(xué)需氧量、五日生化需氧量、總氮、石油類(lèi)、硫化物、揮發(fā)酚,這些指標(biāo)均為水質(zhì)監(jiān)測(cè)指標(biāo)中易檢出或超標(biāo)的因子。水源涵養(yǎng)功能與景觀指標(biāo)間,面積比例和斑塊比例均影響明顯,但方向相反,即景觀破碎化不利于水源涵養(yǎng)功能的保持。水源涵養(yǎng)功能與污染物排放量間,主要影響指標(biāo)是總氮、垃圾、總磷、廢水排放和污染源距離。

對(duì)水源涵養(yǎng)功能與各影響因子進(jìn)行主成分分類(lèi)分析(圖3),指標(biāo)屬性方面：集雨區(qū)面積、緩沖區(qū)面積、林地面積、水庫(kù)容量、污染物排放量等一般性指標(biāo)在第一象限；水源涵養(yǎng)指數(shù)、景觀格局指數(shù)、水體面積、水位主要分布在第二象限；景觀格局指數(shù)的緩沖區(qū)自然/人工面積比、污染源距離在第三象限；水質(zhì)主成分1和主成分2、綜合污染指數(shù)等水質(zhì)指標(biāo),及斑塊密度分布在第四象限。與之對(duì)應(yīng)的大王灘水庫(kù)在第一象限；鳳亭水庫(kù)和屯六水庫(kù)、西云江水庫(kù)分布在第二象限；東山水庫(kù)老虎嶺水庫(kù)、天雹水庫(kù)、峙村河水庫(kù)分布在第三象限；龍?zhí)端畮?kù)分布在第四象限。

圖3 湖庫(kù)型飲用水源結(jié)構(gòu)與功能多維尺度分析Fig.3 The multidemension scale analysis on the structure and the function of the drinking water sources for the reservoirsWCindex：水源涵養(yǎng)功能指數(shù) Water conservation function index；Landscape x1—x7：x1緩沖區(qū)自然/人工斑塊面積比 The natural and artifical patches area ratio of the buffer zone、x2緩沖區(qū)自然/人工斑塊面積—個(gè)數(shù)綜合比 The natural and artifical patches total ratio of the buffer zone、x3緩沖區(qū)自然/人工斑塊個(gè)數(shù)比The natural and artifical patches No. ratio of the buffer zone、x4保護(hù)區(qū)自然/人工斑塊密度 The natural and artifical patches density of the protected zone、x5保護(hù)區(qū)自然/人工斑塊個(gè)數(shù)比 The natural and artifical patches density of the protected zone、x6保護(hù)區(qū)自然/人工斑塊面積比The natural and artifical patches total ratio of the protected zone、x7保護(hù)區(qū)自然/人工面積—個(gè)數(shù)綜合比The natural and artifical patches No. ratio of the protected zone；WS x1—x3：x1水體面積 Water area、x2 水位 Water level、x3庫(kù)容量Reserviors capcity；Area x1—x5：x1保護(hù)區(qū) Protected area、x2 500 m 緩沖區(qū) 500 m buffer zone、x3 3 km 緩沖區(qū) 3 km buffer zone、x4集雨區(qū)Catchment、x5保護(hù)區(qū)林地 The forest of protected area；WFactor 1：水質(zhì)主成分因子1 the PCA factor 1 of water quality；WFactor 2：水質(zhì)主成分因子2 the PCA factor 2 of water quality；WPCindex：水質(zhì)綜合污染指數(shù)Pollutional comprehensive index of the water quality；Pollution x1—x7：x1廢水 Waste water、x2化學(xué)需氧量Chemical oxygen demand、x3氨氮 Ammonia nitrogen、x4總磷 Total phosphorus、x5 總氮Total nitrogen、x6垃圾 Refuse、x7污染源距離 Distance to the pollutant sources;○1—9：水庫(kù)點(diǎn)位代碼 Reservoirs code,1大王灘水庫(kù)、2東山水庫(kù)、3鳳亭河水庫(kù)、4老虎嶺水庫(kù)、5龍?zhí)端畮?kù)、6天雹水庫(kù)、7屯六水庫(kù)、8西云江水庫(kù)、9峙村河水庫(kù)

大王灘和龍?zhí)端畮?kù)均分布在右邊象限,影響因子是一般性指標(biāo)和水質(zhì)指標(biāo),呈人工干擾強(qiáng)烈的指標(biāo)特性,但大王灘目前整治效果較明顯,因此結(jié)構(gòu)自組織的內(nèi)斂機(jī)制將其排列到上方的第一象限。鳳亭、屯六和西云江水庫(kù)均自然度高,故分布在第二象限,與龍?zhí)端畮?kù)呈相反特性。東山水庫(kù)、老虎嶺水庫(kù)、天雹水庫(kù)、峙村河水庫(kù)景觀本底上是自然度高于人工干擾性,但污染源距離、緩沖區(qū)自然/人工面積比2個(gè)自然與人工特性指標(biāo)的偶然性、靈敏性較強(qiáng),其不確定影響可能改變本底特性,引起突變事故,即生態(tài)系統(tǒng)自然穩(wěn)定度與調(diào)節(jié)功能很脆弱,抗性低,而易系統(tǒng)崩潰。龍?zhí)端畮?kù)作為水質(zhì)影響最明顯的案例,其自然度最低,污染較嚴(yán)重,水質(zhì)呈富營(yíng)養(yǎng)狀態(tài)。

耦合模型總變量解釋度99%,功能屬性變量解釋度89.2%,空間結(jié)構(gòu)變量解釋度96.6%,為一理想的主分量分類(lèi)耦合模型。

3 討論

3.1 景觀分類(lèi)與尺度分析

景觀分析始于景觀分類(lèi),常見(jiàn)有土地利用/下墊面、自然/人工、源-匯、生態(tài)界面等分類(lèi)系統(tǒng)；不同景觀分類(lèi)延伸的科學(xué)內(nèi)涵及其意義各有特點(diǎn),又存共性。岳雋等應(yīng)用生態(tài)學(xué)“源-匯理論”分類(lèi)景觀要素,判定深圳市西部庫(kù)區(qū)流域污染物輸入“匯”景觀格局空間變化對(duì)水體水質(zhì)有密切關(guān)聯(lián)[13]。趙志軒等[14]探討了自然和人類(lèi)活動(dòng)對(duì)流域景觀類(lèi)型、規(guī)模、組成、布局發(fā)生改變的二元驅(qū)動(dòng)機(jī)制,發(fā)現(xiàn)景觀空間梯度變化在過(guò)渡地帶波動(dòng)幅度相對(duì)較大；張華兵等[19]將鹽城國(guó)家級(jí)自然保護(hù)區(qū)核心區(qū)劃分為人工管理區(qū)和自然濕地區(qū)兩種模式,對(duì)濕地景觀格局驅(qū)動(dòng)力變化進(jìn)行分析表明：自然過(guò)程如地貌過(guò)程和植物群落演替過(guò)程等呈連續(xù)性變化,從而導(dǎo)致景觀演變呈帶狀連續(xù)性發(fā)展,人為管理往往使生態(tài)過(guò)程的連續(xù)性發(fā)生突變,從而使景觀演變呈多向性特征。為了生物遷移、物質(zhì)流、能量流和信息流在邊界判定和空間動(dòng)態(tài)數(shù)學(xué)模擬不受系統(tǒng)和尺度限制,Calnasso等提出生態(tài)界面的科學(xué)概念[20],由于生態(tài)界面研究?jī)?nèi)容上注重結(jié)構(gòu)、功能與過(guò)程的變異分析和動(dòng)態(tài)模擬,方法上集成數(shù)學(xué)、統(tǒng)計(jì)學(xué)和工程模型的綜合應(yīng)用,正成為引領(lǐng)景觀生態(tài)研究中的興趣熱點(diǎn),其同義表達(dá)則可追溯為邊緣、交錯(cuò)帶、界面層、梯度、漸變?nèi)骸⑦^(guò)渡帶、交界面、邊界和界面等[21]。本研究采用國(guó)家規(guī)范標(biāo)準(zhǔn)的土地利用、自然/人工的景觀分類(lèi)方法,與“源-匯”景觀分類(lèi)有異曲同工之處,而湖庫(kù)岸邊帶又是生態(tài)交錯(cuò)區(qū)域,具有生態(tài)界面的科學(xué)內(nèi)涵；500 m 和3 km 兩種緩沖帶的自然/人工景觀指數(shù)結(jié)果表明湖庫(kù)型飲用水源人工干擾性增強(qiáng),景觀破碎度就越高,從而影響水體水質(zhì)和水源涵養(yǎng)；部分湖庫(kù)應(yīng)根據(jù)“集中與分散相結(jié)合”的最優(yōu)景觀格局原理,在關(guān)鍵地段對(duì)“源匯”景觀或自然/人工景觀予以合理配置和協(xié)調(diào)[13-14,22-23]。

尺度是正確理解景觀格局和過(guò)程變化的重要范式,極具復(fù)雜性和多樣性,尺度效應(yīng)及其適宜分析是景觀生態(tài)科學(xué)研究的重要基礎(chǔ)內(nèi)容。尺度常見(jiàn)有空間尺度、時(shí)間尺度和方法尺度等類(lèi)型,如蔡博峰等[22]用半方差分析、尺度方差、小波分析和孔隙度指數(shù)分析等4種尺度算法判定三北防護(hù)林地區(qū)的林地景觀空間格局特征尺度約為510—530 km。王小平等[23]用1、2、3、4、5 km等5個(gè)緩沖區(qū)尺度對(duì)艾比湖區(qū)域景觀格局與水質(zhì)進(jìn)行關(guān)聯(lián)分析,基于去趨勢(shì)分析方法的方差解釋量比較,結(jié)合Bootstrip度量梯度變異的關(guān)鍵區(qū)間,認(rèn)為4 km緩沖區(qū)景觀格局具有最大的水質(zhì)空間分異解釋能力。基于國(guó)家飲用水源保護(hù)區(qū)劃分規(guī)范標(biāo)準(zhǔn),本研究采用500 m和3 km 兩種緩沖區(qū)分析了湖庫(kù)型飲用水源景觀格局,但從景觀整體性角度考慮,水體統(tǒng)一包含于緩沖區(qū)內(nèi),則可區(qū)別于規(guī)范標(biāo)準(zhǔn)劃定的一級(jí)保護(hù)區(qū)和二級(jí)保護(hù)區(qū)相互間的約束性分割和差別化交錯(cuò),而景觀格局分析和水源涵養(yǎng)功能解釋又能滿足于兩種保護(hù)區(qū)劃定范圍。劉文竹等[13]對(duì)密云水庫(kù)小流域不同尺度景觀格局分析認(rèn)為同一河段不同尺度下,100 m×500 m尺度下的破碎化程度高于100 m×1 000 m尺度,即較長(zhǎng)尺度的破碎程度越低,而本研究結(jié)果是500 m較3 km的緩沖區(qū)景觀破碎程度低；結(jié)果相反的原因是一個(gè)面向河流,一個(gè)面向湖庫(kù),兩者的岸邊帶人類(lèi)活動(dòng)對(duì)水系依存度不同,河流高于湖庫(kù)；另外是南北方氣候和土地利用方式有差異,北方土地耕作強(qiáng)度高,而本研究中500 m緩沖區(qū)均屬于保護(hù)區(qū)范圍,僅有一些養(yǎng)殖活動(dòng)和岸邊帶設(shè)施工程,耕作強(qiáng)度較高者僅在龍?zhí)端畮?kù)周邊。

3.2 耦合模式與特征表達(dá)

構(gòu)建耦合模型是研究格局-過(guò)程相互關(guān)系的重要途徑[24],模型可以充分利用實(shí)驗(yàn)和觀測(cè)數(shù)據(jù)并綜合不同時(shí)間和空間尺度上的信息提煉規(guī)律或揭示內(nèi)在機(jī)制,反演景觀格局特征、生態(tài)過(guò)程發(fā)生/發(fā)育狀況和關(guān)鍵環(huán)節(jié),預(yù)測(cè)發(fā)展趨勢(shì)。通過(guò)擬合景觀格局指數(shù)與環(huán)境因子在梯度分析上的最適排序,可以探討景觀梯度格局中驅(qū)動(dòng)機(jī)制、尺度效應(yīng)和影響因子[14]。借鑒水質(zhì)-水量[2]、景觀-水質(zhì)[13- 16,25]的概念模式與方法,本研究基于多元逐步回歸分析和主分量分類(lèi)分析方法,定性與定量的耦合了水源涵養(yǎng)功能與景觀格局空間結(jié)構(gòu)、水質(zhì)、水資源和污染物排放等影響因子間模型關(guān)系,擬合/耦合模型集方程解釋、圖形表達(dá)、空間與屬性相關(guān)/關(guān)聯(lián)于一體,其中主分量分類(lèi)分析的變量方差總解釋度[23,25]達(dá)99.8%,空間解釋度達(dá)96.6%,指標(biāo)屬性解釋度達(dá)89.2%,耦合模型理想。其分析結(jié)果可推斷出天雹水庫(kù)、東山水庫(kù)和老虎嶺水庫(kù)極易由這些相關(guān)性因子引發(fā)生態(tài)系統(tǒng)的脆弱性,降低其水源涵養(yǎng)調(diào)節(jié)功能的穩(wěn)定作用；龍?zhí)端畮?kù)一方面現(xiàn)有水源地保護(hù)區(qū)生態(tài)系統(tǒng)結(jié)構(gòu)不足以承載水源涵養(yǎng)功能,且景觀格局強(qiáng)烈破碎化,幾無(wú)外向型結(jié)構(gòu)整合優(yōu)化空間,另一方面養(yǎng)殖和農(nóng)業(yè)污染面源對(duì)水庫(kù)水質(zhì)影響嚴(yán)重；因此需要優(yōu)化生態(tài)空間[23,26],加強(qiáng)農(nóng)業(yè)污染面源的整治,限制岸邊帶養(yǎng)殖；鳳亭河水庫(kù)與屯六水庫(kù)功能定位不明確將影響未來(lái)飲用水安全保障問(wèn)題,演生較多的水資源利用與生態(tài)功能錯(cuò)位情況,尤其是部分地區(qū)水資源先天充沛卻面臨無(wú)水可用的情況,出現(xiàn)人為“搶水”行為導(dǎo)致空間性缺水問(wèn)題[2]需要引起政府及有關(guān)方面的關(guān)注。另外,部分湖庫(kù)型水源地周邊桉樹(shù)林所占面積比例較高,需要更大區(qū)域的林地保持量才能保障水源涵養(yǎng)的功能；然而,在人地矛盾突出,環(huán)境保護(hù)與經(jīng)濟(jì)發(fā)展雙重壓力下,劃出更大區(qū)域的做法無(wú)疑面臨較大挑戰(zhàn)。

3.3 景觀結(jié)構(gòu)與功能優(yōu)化

通過(guò)科學(xué)構(gòu)建生態(tài)源地、生態(tài)廊道和生態(tài)節(jié)點(diǎn)等景觀組分,可優(yōu)化景觀格局/結(jié)構(gòu),有效提高景觀效益和穩(wěn)定性,而其關(guān)鍵在于保護(hù)核心斑塊源地、提高景觀連通度,如孫賢斌等[26]基于生態(tài)功能評(píng)價(jià)及其效應(yīng)分析,用最小景觀耗費(fèi)強(qiáng)度-景觀連通指數(shù)的方法對(duì)鹽城濕地生態(tài)景觀格局和服務(wù)功能進(jìn)行了探索性優(yōu)化。基于水量平衡理論,本研究通過(guò)對(duì)11種森林類(lèi)型水源涵養(yǎng)最小面積進(jìn)行測(cè)算以判定和確定現(xiàn)有水源保護(hù)區(qū)邊界的適宜度和調(diào)整優(yōu)化空間,提出優(yōu)化現(xiàn)有飲用水源保護(hù)區(qū)內(nèi)森林生態(tài)系統(tǒng)的林分類(lèi)型,提高適生軟闊葉林所占比重,以及在湖庫(kù)部分水體區(qū)域營(yíng)造濕地生態(tài)系統(tǒng),從而提高生態(tài)系統(tǒng)、群落和物種的多樣性[1-2,18],增強(qiáng)水質(zhì)凈化能力和涵養(yǎng)功能。從森林生態(tài)學(xué)與造林學(xué)角度看,影響森林水源涵養(yǎng)功能的因素還包括林木水平分布格局、樹(shù)種混交、種內(nèi)及種間競(jìng)爭(zhēng)、多層次垂直結(jié)構(gòu)等因素,因此,有森林生態(tài)學(xué)研究人員基于系統(tǒng)演替理論,提出結(jié)構(gòu)化經(jīng)營(yíng)森林的水源涵養(yǎng)林空間結(jié)構(gòu)多目標(biāo)優(yōu)化模型[27],通過(guò)目標(biāo)樹(shù)種空間位置更新,慣性因子(林隙)調(diào)整,適應(yīng)度函數(shù)設(shè)置(主要選用混交度、競(jìng)爭(zhēng)指數(shù)、角尺度、林層指數(shù)、空間密度指數(shù)、開(kāi)闊比數(shù)作為林木調(diào)整的目標(biāo)函數(shù)),借助智能粒子群優(yōu)化(PSO)運(yùn)算,得出目標(biāo)樹(shù)種及其伴生樹(shù)種達(dá)到理想頂級(jí)群落狀態(tài)下的擇伐、補(bǔ)植空間結(jié)構(gòu)單元,使目標(biāo)樹(shù)均質(zhì)性指數(shù)增加,其空間布局趨向合理,從而整體上優(yōu)化林分結(jié)構(gòu),促進(jìn)生態(tài)系統(tǒng)水源涵養(yǎng)功能的健康和穩(wěn)定[18,23,28]。

4 結(jié)論

基于多源數(shù)據(jù)、遙感技術(shù)和多元分析方法,對(duì)南寧市湖庫(kù)型飲用水源景觀結(jié)構(gòu)與涵養(yǎng)功能關(guān)系耦合優(yōu)化進(jìn)行了探索性研究,結(jié)論如下：

(1)南寧市湖庫(kù)型飲用水源隨著緩沖區(qū)域由500 m向3 km延伸擴(kuò)展,受人類(lèi)活動(dòng)干擾加強(qiáng),景觀破碎度有所增加,水源涵養(yǎng)功能減弱,形成景觀格局梯度。

(2)南寧市湖庫(kù)型水源涵養(yǎng)功能受景觀格局、水資源、水質(zhì)、污染物排放等因子影響,相互之間存在數(shù)量相關(guān)性或關(guān)聯(lián)性,可建立起多重?cái)M合優(yōu)化模型方程；基于主分量分類(lèi)法的多維尺度耦合模型可揭示當(dāng)前南寧市湖庫(kù)型飲用水源存在問(wèn)題和面臨壓力。

(3)“最小面積-多維空間耦合”方法的綜合應(yīng)用可優(yōu)化南寧市湖庫(kù)型飲用水源景觀結(jié)構(gòu)與水源涵養(yǎng)功能,診斷飲用水源在生態(tài)系統(tǒng)穩(wěn)定性、水資源量和水質(zhì)方面的潛在生態(tài)安全與風(fēng)險(xiǎn)問(wèn)題,對(duì)當(dāng)前飲用水源生態(tài)保護(hù)與建設(shè)、資源合理開(kāi)發(fā)與保護(hù)、經(jīng)濟(jì)發(fā)展方式轉(zhuǎn)變、區(qū)域生態(tài)承載能力增強(qiáng)等工作具有借鑒參考意義。

致謝：南寧市規(guī)劃局、南寧市環(huán)境保護(hù)局、廣西壯族自治區(qū)林業(yè)勘查設(shè)計(jì)院提供相關(guān)研究資料，特此致謝。