席浩郡 袁一斌 昝曉輝 李天宏,?
1.北京大學環(huán)境科學與工程學院, 北京 100871; 2.國家環(huán)境保護河流全物質通量重點實驗室, 北京 100871;3.成都市環(huán)境保護科學研究院 成都 610072; ? 通信作者, E-mail: lth@pku.edu.cn
河流是地球上物質循環(huán)和能量遷移的主要通道, 河流的健康狀況對人類健康和社會經(jīng)濟發(fā)展有重要影響[1]。河流系統(tǒng)的結構和功能會因自然與人類活動的作用而不斷發(fā)生時空變化, 一旦這些變化超出河流系統(tǒng)自身負反饋調節(jié)能力, 河流健康就會受到損害[2]。因此, 如何以維護河流健康為目標對河流進行科學管理, 是一個亟待解決的關鍵問題。“河流健康”一詞最早由 Scrimgeour 等[3]提出, 認為河流健康管理中應該關注生態(tài)系統(tǒng)的完整性。健康的河流應該是穩(wěn)定、完整、充滿活力和有彈性的[4]。隨著社會的發(fā)展, 河流健康管理關注的對象由自然生態(tài)系統(tǒng)轉到受人類主導的社會–經(jīng)濟–自然復合生態(tài)系統(tǒng)[5], 認為健康的河流除要符合所在流域特征的結構、功能和生物多樣性, 對長期或突發(fā)的擾動能保持彈性、穩(wěn)定性以及一定的自我恢復和適應能力外, 還能持續(xù)地為人類提供社會經(jīng)濟服務功能,包括供水、通航、行洪、輸沙、科研文化和景觀娛樂等[6]。
河流健康評價是河流健康管理的重要依據(jù), 通過建立指標體系, 對河流現(xiàn)狀和問題進行評價和診斷, 從而確立河流生態(tài)修復的目標, 并對修復過程進行監(jiān)控, 為河流管理提供信息支撐。近年來, 河流健康評估方法的研究已經(jīng)取得很大的進展, 學者們提出許多方法進行河流健康評價[7], 包括歐盟水框架指令(WFD)、美國快速生物評價規(guī)程(RBPs)、澳大利亞溪流狀況指數(shù)(ISC)、英國河流無脊椎動物預測分類系統(tǒng)(RIVPACS)和南非河流健康計劃(RHP)等。河流健康評價的方法一般分為生物指標評價法和多指標綜合評價法。生物指標主要包括浮游植物、浮游動物、底棲無脊椎動物和魚類等水生生物, 能夠較好地反映水體健康的物理、化學和生物效應[8]。在河流水生生物的監(jiān)測中, 環(huán)境 DNA(eDNA)分析是一種快速且成本低效益高的生物多樣性監(jiān)測工具, 克服了傳統(tǒng)取樣的局限性, 采樣簡便, 分析準確度高, 在水生生物保護和管理方面有巨大的應用潛力[9]。例如, Li 等[10]采用環(huán)境 DNA宏條形碼技術, 揭示長江與太湖流域河流水質與生態(tài)健康。然而, 由于社會–經(jīng)濟–自然復合生態(tài)系統(tǒng)的復雜性, 在空間和時間上均具有多樣的變量函數(shù),因此需要采用多指標綜合評價的方法為河流生態(tài)系統(tǒng)提供更加全面的評價。多指標綜合評價指標體系通常包括物理形態(tài)、水質、水文、水生生物和河岸帶等指標, 所選指標的原始數(shù)據(jù)可以通過歷史資料分析、實地調研監(jiān)測和遙感、GIS 等手段獲得, 再通過層次分析法、模糊綜合評價法、數(shù)據(jù)包絡分析、灰色關聯(lián)分析法和人工神經(jīng)網(wǎng)絡等[11]多種數(shù)學方法進行評價。
成都平原岷江內江水系是自然條件和人類活動綜合作用的結果。近幾十年來, 快速城市化和經(jīng)濟發(fā)展過程中, 人類活動對河流水系的影響尤其明顯,這些河流水系屬于典型的社會–經(jīng)濟–自然復合生態(tài)系統(tǒng)。當前, 流域面臨斷面水質達標不穩(wěn)定、水資源開發(fā)強度大、水生態(tài)系統(tǒng)整體性需提升等問題。因此, 本文以流域中典型河流柏條河–府河為研究對象, 采用多指標綜合評價指標體系, 建立滿足河流生態(tài)系統(tǒng)完整性和人類服務需求的河流健康評價方法, 分析不同河段河流健康水平和空間變化特征,以期為成都平原河流的可持續(xù)管理提供科學依據(jù),并為其他深受人類活動影響河流系統(tǒng)的健康評價探索新方法。
錦江流域, 又稱岷江內江流域, 位于四川省成都市中部, 地勢西北高、東南低, 總面積約為 2380 km2。岷江是長江上游的一條重要支流, 經(jīng)都江堰水利工程分為內、外二江, 內江通過寶瓶口流入成都平原, 從天府新區(qū)黃龍溪出境, 流域概況如圖1所示。研究區(qū)屬于亞熱帶濕潤季風氣候區(qū),年平均氣溫為 17oC,年降水量約為 860 mm, 季節(jié)分配不均, 主要集中在 6—9月。錦江流域是成都市社會經(jīng)濟最集中的區(qū)域, 包括都江堰市、溫江區(qū)、郫都區(qū)、雙流區(qū)和天府新區(qū)成都直管區(qū)等 11個區(qū)。流域內人口密集, 常住人口約占全市總人口的 45%,全 流 域 GDP 為 1.03×104億元 , 約占全市 GDP 的60%。研究區(qū)域內土地類型多樣, 以耕地和建設用地為主, 分別占 52.4%和 41.5%。
圖1 研究區(qū)位置Fig.1 Location of study area
柏條河是都江堰內江四大干渠之一, 是成都市自來水六廠、七廠的主要水源河道。府河是古代成都通往中國沿海的黃金水道, 根據(jù)《成都港總體規(guī)劃》(草案), 成都港將成為四川省的重要港口, 主要發(fā)展旅游客運, 兼顧發(fā)展水上公交。府河亦匯合了成都市境內其他的河流, 包括南河、沙河和江安河等。
河流健康評價指標體系應遵循系統(tǒng)性、完備性和可操作性相結合的原則[12]。本研究將河流健康評價設計成目標層、準則層和指標層的層次結構,目標層為“河流健康”, 基于健康河流應具備的生態(tài)環(huán)境功能和社會經(jīng)濟服務功能, 選取生境特征、生物群落和社會服務為準則層。每個準則層由若干狀態(tài)層構成, 每個狀態(tài)層由一個或若干具體指標進行表征。本文中指標參考水利部《河湖健康評估技術導則》[13]和文獻[11–12,14–21], 結合流域特征, 共選取 16個評價指標, 建立河流健康評價體系。表1給出各個指標的計算依據(jù)和量化法。
表1 河流健康評價指標體系Table 1 Index system of river health assessment
2.1.1 生境特征
生境特征主要表現(xiàn)在河流水文、水質和物理結構上, 其中物理結構包括河流本身的形態(tài)、河岸帶的形態(tài)以及河流周邊的土地利用。生態(tài)需水滿足程度反映河流水資源配置時生態(tài)用水的保障程度, 決定著河流系統(tǒng)的穩(wěn)定性和多樣性, 是河流生態(tài)健康的核心[22]。因此, 選擇生態(tài)基流滿足程度作為水文評價指標。由于錦江流域河流的流量受人類影響劇烈, 大部分水文測站已不能反映天然徑流的狀況,且不能有效地進行還原計算, 故水文變異值不納入評價指標。水質指標中采用流域主要污染物的水質綜合指數(shù)反映河流的污染程度, 對內江流域地表水的水質指標進行篩選, 識別出主要污染物為氨氮、總磷和化學需氧量。此外, 采用溶解氧含量來反映河流生態(tài)系統(tǒng)的自凈功能。水質監(jiān)測斷面位置如圖1 所示。
選取河流縱向連通性作為河流形態(tài)指標, 用來表示生物、物質、能量在河流縱向上運移的通暢程度, 使用 Google Earth 軟件計算單位河長內影響河流連通性的建筑物或設施數(shù)量[16]。河流生態(tài)緩沖帶是河流生態(tài)系統(tǒng)的重要載體, 緩沖帶生態(tài)修復是《“十四五”重點流域水生態(tài)環(huán)境保護規(guī)劃》的重要工作。因此, 選取河岸帶穩(wěn)定性、河岸帶植被覆蓋率和河岸緩沖帶寬度作為河岸緩沖帶指標。人類活動對河流生態(tài)系統(tǒng)的影響越來越深刻, 而土地利用強度是河流生態(tài)系統(tǒng)健康的第一位驅動因素[23]。因此, 利用 ArcGIS 10.6 計算土地利用強度, 表征由自然因素和人為因素對河流健康雙重作用的結果。
2.1.2 生物群落
生物群落與生境的統(tǒng)一性是生態(tài)系統(tǒng)完整性的基本特征, 生物群落是生態(tài)系統(tǒng)的主體, 生境是生物群落的生存條件[24]。藻類、底棲生物和魚類都是河流健康評價的重要指示物種, 對水生生態(tài)系統(tǒng)的物質循環(huán)、能量流動及維持生態(tài)系統(tǒng)平衡起重要作用[25]。Shannon-Wiener 指數(shù)綜合考慮群落的豐富度和均勻度, 數(shù)值越大, 表明群落結構的復雜程度和穩(wěn)定性越高, 生物多樣性越高[8]。本研究依托北京大學構建的以環(huán)境樣品遺傳物質 eDNA 為核心的水生生物監(jiān)測技術, 獲取研究區(qū)域內藻類以及底棲動物和魚類生物群落的多樣性與結構特征, 并采用Shannon-Wiener指數(shù)表述藻類、底棲生物和魚類多樣性狀況。采樣點位置如圖1 所示。
2.1.3 社會服務
本文采用防洪工程措施達標率來表征河流的輸水泄洪功能, 采用水資源開發(fā)利用率和萬元 GDP用水量來表征河流的供水功能。水資源開發(fā)利用率是判斷人類開發(fā)利用的水資源是否擠占生態(tài)用水的重要標準, 萬元 GDP 用水量反映人類對水資源的利用效率[23], 這兩個指標從不同方面間接地反映供水功能及其所受人類社會活動的影響。采用公眾滿意度和水文化傳承載體數(shù)量分別反映河流的景觀娛樂功能和文化傳承功能。
根據(jù)河流健康評價體系, 確定各單項指標的現(xiàn)狀值。根據(jù)河流健康評價指標評價標準, 將各指標現(xiàn)狀值進行量化打分, 優(yōu)、良、中、差和劣 5個級別分別對應量化分值 5, 4, 3, 2 和 1。各指標劃分方法包括: 1) 采取現(xiàn)有國家、地方、行業(yè)標準或國際標準; 2) 結合相關文獻中相似研究區(qū)域的指標分級情況; 3) 專家咨詢。表2 給出各指標賦值標準。
表2 河流健康評價指標賦值標準Table 2 Evaluation standard for river health
生境特征、生物群落和社會服務均是河流健康的主要組成部分, 因此, 本研究采用等權的綜合指數(shù)法進行健康評價, 認為生境特征、生物群落和社會服務的權重相同, 每個準則層內的健康指標權重相同。采用綜合評價模型計算河流健康評價指數(shù)(River Health assessment Index, RHI), 計算公式為
式中,Ii為第i個指標的得分值,wi為第i個指標的權重。
將河流健康評價結果劃分為健康、亞健康、中等、亞病態(tài)和病態(tài) 5個等級, 見表3。
表3 河流健康狀況等級劃分表Table 3 Classification of river health status
在河流健康評價中, 對于大中型河流, 由于涉及面積廣, 自然條件和社會經(jīng)濟背景差異大, 需要對河流進行分段評價[28]。結合流域地貌形態(tài)、水文分區(qū)、水功能分區(qū)、土地利用狀況和行政區(qū)劃等, 將錦江流域的柏條河–府河河流分為 4 段(圖1),從上游到下游依次為柏條河(I)、府河上游(II)、府河中游(III)和府河下游(IV)。
從圖2 可以看出, 不同河段間生境特征指標和社會服務指標差別較大, 生物群落指標差別較小。柏條河 N4, N10 和 N11 的賦值小于 3, 府河上游 N4,N12 和 N13 的賦值小于 3, 府河中游 N4, N8, N11 和N13 的賦值小于 3, 府河下游 N4, N10, N11, N12 和N13 的賦值小于 3。
圖2 指標層得分Fig.2 Health score of index layer
對評價河段的健康評價指標進行計算, 不同河段指標賦值情況如圖2 所示。
3.2.1 生境特征
利用綜合評價模型計算各河段準則層的健康指數(shù), 結果如圖3 所示, 可以看出, 柏條河–府河生境特征的健康指數(shù)取值范圍為 3.38~4.25。從河流健康分級情況(圖4)可以看出, 除柏條河的生境特征評價等級為健康以外, 其余河段的生境特征評價等級為亞健康。
圖3 柏條河-府河河流健康評價指數(shù)Fig.3 RHI of Baitiao-Fu River
圖4 柏條河-府河河流健康狀況分級Fig.4 Health status classification of Baitiao-Fu River
從水文指標看, 各評價河段的生態(tài)基流滿足程度均大于 90%, 評價等級為優(yōu)。從水質指標看, 柏條河–府河河流斷面水質綜合指數(shù)范圍為 0.43~1.20,府河中下游斷面的水質最差, 主要超標污染物為氨氮和總磷。出現(xiàn)這種結果的主要原因是府河中下游流經(jīng)全市人口最為密集、商業(yè)最集中的中心城區(qū),南河、沙河、江安河和鹿溪河等支流沿程依次匯入。由于流域產(chǎn)污量大以及周邊區(qū)縣的污水收集處理率不高, 導致城鎮(zhèn)生活直排污水成為污染物的主要來源。此外, 受到降水影響, 雨季時非點源污染負荷增加, 導致河流污染物通量增加。府河設置較多的水閘和橡膠壩, 河流縱向連通性破壞嚴重, 對洄游魚類造成嚴重影響[29]。唐家璇等[16]對長江流域河流縱向連通性特征進行分析, 結果表明, 近 60年來河流縱向連通性呈現(xiàn)明顯下降趨勢, 評價等級為劣的河流主要為岷沱江、嘉陵江和洞庭湖水系,與本文對錦江流域中柏條河–府河的評價結果一致。結合實地勘查結果與遙感影像目視解譯, 對河岸帶穩(wěn)定性、護岸類型、河岸帶植被覆蓋率以及河岸緩沖帶寬度進行分析。
柏條河以自然土質岸坡為主, 河岸穩(wěn)定性評價結果為“良”, 植被覆蓋率大于 75%, 緩沖帶寬度 30 m。府河上游以自然土質岸坡和生態(tài)護岸為主, 穩(wěn)定性較好, 植被覆蓋率為50%~75%, 緩沖帶寬度范圍為 20~50 m。府河中游位于中心城區(qū), 河流護岸以直立式漿砌石護岸為主, 河道渠化明顯, 穩(wěn)定性強, 植被覆蓋率為 50%~75%, 河流臨近道路, 緩沖帶寬度范圍為 2~15 m。府河下游以自然土質岸坡為主, 植被覆蓋率大于 75%, 左右兩岸的緩沖帶寬度均在 20 m 以上。但是, 受降水沖刷和水土流失的影響, 河岸穩(wěn)定性為河段中最差。
沿河流的兩岸向外延展 1 km[11]作為緩沖區(qū), 通過對遙感圖像進行解譯, 獲得各土地利用類型的面積數(shù)據(jù), 各河段土地利用強度指數(shù)范圍為 3.31~4.88, 府河中游土地利用類型以城鎮(zhèn)用地為主(93.1%), 土地利用強度最大。柏條河和府河下游土地利用類型以耕地和農(nóng)村居民點為主, 土地利用強度較小。
3.2.2 生物群落
采用環(huán)境 DNA (environmental DNA, eDNA)技術, 對柏條河–府河采樣點的藻類、底棲動物和魚類的生物多樣性進行分析。在流域中檢出藻類 15門 35 綱 76 目 127 科 197 屬 329種 , 硅藻門(Bacilla
riophyta)、綠藻門(Chlorophyta)和隱藻門(Cryptophyta)在空間分布中占主要優(yōu)勢。流域中有底棲動物共 6 門 8 綱 22 目 66 科 139 屬 192種, 優(yōu)勢科主要為搖蚊科(Chironomidae)、水???Hydridae)和水虻科(Stratiomyidae)等。流域中的魚類都屬于脊索動物門條鰭魚綱, 共 7 目 20 科 28 屬 33種, 優(yōu)勢屬主要為鯉屬(Cyprinus)、頜須鮈屬(Gnathopogon)和云南鰍屬(Yunnanilus)等。
從圖3 可以看出, 柏條河–府河生物群落的健康指數(shù)取值范圍為 2.67~4.00, 府河上游的生物群落評價等級為健康, 柏條河和府河中游的生物群落評價等級為亞健康, 府河下游的生物群落評價等級為中等。從各斷面的環(huán)境 DNA 檢測結果看, 所有點位的藻類多樣性等級都在“中等”以上, 66.7%點位的底棲動物多樣性等級在“中等”以上, 16.7%點位的魚類多樣性等級為“中等”, 其余點位魚類多樣性等級為“差”??偟膩碚f, 多樣性指數(shù)較低的點位多位于府河中下游, 且多樣性指數(shù)有隨距離增加而衰減的趨勢, 但變化幅度不大。
與文獻[30]對比, eDNA 技術相較于傳統(tǒng)形態(tài)學可以識別出更多的物種, 促進研究人員對復合生態(tài)系統(tǒng)的理解。此外, 要真正了解流域內物理因子(溫度和流速等)、化學因子(溶解氧、氮磷和重金屬等)、生物因子(種內種間的競爭和捕食)和人為干擾(城市化以及閘壩建設)等對物種多樣性和功能完整性的影響機制, 需要與傳統(tǒng)形態(tài)學方法相結合, 獲得更全面的生物信息, 并需要進行多次采樣調查分析。
3.2.3 社會服務
從評價結果看, 柏條河–府河社會服務的健康指數(shù)取值范圍為 3.40~4.20, 柏條河和府河中游的社會服務評價等級為健康, 府河上游和下游的社會服務評價等級為亞健康。
對于防洪指標, 中心城區(qū)的防洪措施達標率已達到 100%, 府河上游和下游區(qū)域河段的防洪措施達標率均小于 60%, 堤防建設長度仍有待加強。對于供水指標, 根據(jù)《成都市水資源公報》計算, 柏條河、府河上中下游的水資源開發(fā)利用率分別為35%, 143%, 332%和 93%, 府河的水資源開發(fā)利用率超過國際公認的 40%的開發(fā)利用紅線[31]。柏條河以及府河上、中、下游的萬元 GDP 用水量分別為 52.2, 46.5, 9.3 和 39.9 m3, 除柏條河的評價等級為良外, 其他河段的評價等級為優(yōu)。流域內水資源可利用量包括本地水和岷江過境水, 而本地水資源時空分布不均: 在時間分配上, 豐水期和枯水期懸殊較大, 季節(jié)性缺水嚴重; 在空間分配上, 徑流深由西北向東南遞減, 缺少調蓄設施調節(jié)天然徑流[32]。岷江過境水資源豐富, 彌補了本地水資源的不足。但是, 受全球氣候變化和人類活動的影響, 河流水文條件發(fā)生較大的變化, 岷江上游來水減少[33]。
對于文化傳承指標, 根據(jù)《成都市城市總體規(guī)劃(2016—2035)》[34], 錦江流域河流具有豐富多樣的水文化。因此, 評價河段的水文化傳承體數(shù)量等級均為優(yōu)。受數(shù)據(jù)限制, 僅從水文化傳承載體數(shù)量對河流文化功能進行評價具有一定局限性, 從水文化傳承、發(fā)展和創(chuàng)新等方面對河流文化功能進行量化是未來研究的發(fā)展方向。對于景觀娛樂指標, 通過分發(fā)問卷, 對柏條河-府河的公眾滿意度進行調查。調查內容包括水質、水生態(tài)、岸線衛(wèi)生、植被綠化和水文化宣傳等, 共收回 51 份有效問卷, 公眾滿意度范圍為 76.0%~98.0%。其中府河上游平均公眾滿意度最高(96.0%), 府河中游的公眾滿意度最低(81.7%)。在社會–經(jīng)濟–自然復合生態(tài)系統(tǒng)中, 河流景觀娛樂功能受多種因素影響。后續(xù)研究中應從景觀格局異質性和連通性、舒適宜居設施以及自然景觀保護設施等方面進一步完善河流景觀娛樂功能評價指標。
本研究充分考慮河流健康內涵和特征, 構建基于生態(tài)完整性指標和社會服務功能的河流健康評價指標體系, 較全面地反映河流系統(tǒng)自身運行狀況和社會服務水平。所選各項指標的原始數(shù)據(jù)可以通過實地調研監(jiān)測、遙感、GIS 和環(huán)境 DNA 等手段獲得, 具有簡潔明了、計算方便的特點。河流健康綜合評價結果(圖3 和 4)表明, 各評價河段的健康指數(shù) RHI 范圍為 3.23~3.82, 根據(jù)各河段長度占總河長的比例, 確定權重, 得到柏條河–府河的 RHI 為3.55, 整體上處于亞健康水平, 尚未實現(xiàn)“人水和諧”的目標。柏條河的河流健康狀況最好, 接近健康水平, 生境特征、生物群落和社會服務的健康狀況均在“良”以上。由于柏條河是成都市的重要飲用水源地, 水資源保護、水生態(tài)保護等基礎設施均優(yōu)于其他河段。府河上游和中游健康狀況為“亞健康”, 府河下游的健康狀況最差, 接近中等水平, 其中生物群落的評分最低(2.67)。
總體來看, RHI 從上游到下游呈現(xiàn)沿程降低的趨勢, 上游的健康狀況明顯好于中下游, 反映不同河段的河流健康水平存在空間異質性, 表明人類活動對河流系統(tǒng)健康狀況造成影響。歐陽莉莉等[35]選取水質、生境和生物指標對 2016年成都市河流健康狀況進行評價, 結果表明岷江上游河流的生態(tài)環(huán)境狀況良好, 中下游河流生態(tài)健康狀況等級為差或很差。評價結果與本文結果的趨勢大體上一致,也說明幾年來這一趨勢仍未改變。此外, 本文在生態(tài)完整性的基礎上增加了社會服務功能指標, 從復合生態(tài)系統(tǒng)的角度說明河流的健康水平, 本文的評價結果更接近實際情況, 具有較強的適用性與可靠性。
本文基于河流健康的內涵, 構建基于生境特征、生物群落和社會服務的河流健康評價體系, 從定量的角度評價錦江流域典型河流柏條河–府河的健康水平, 主要結論如下。
1) 不同河段間各生境特征指標和社會服務指標差別較大, 生物群落指標差別較小。各河段生境特征、生物群落和社會服務的健康指數(shù)范圍分別為3.38~4.25, 2.67~4.00 和 3.40~4.20, 其中府河下游有5個指標處于“差”或“劣”的狀態(tài)。
2) 河流健康綜合評價指數(shù)范圍為 3.23~3.82,健康狀況處于亞健康水平。河流健康評價指數(shù)呈現(xiàn)從上游到下游沿程降低的趨勢, 反映不同河段的河流健康水平存在空間異質性, 表明人類活動對河流系統(tǒng)健康造成干擾。
3) 本文提出的河流健康評價框架采用遙感、GIS 和環(huán)境 DNA 等新型技術手段, 評價結果具有較強的適用性與可靠性, 可為柏條河–府河以及該地區(qū)類似河流的健康管理提供科學工具。
河流這一復合生態(tài)系統(tǒng)具有復雜的分層等級結構, 評價尺度的不同會導致對河流各組成因素間相互作用規(guī)律產(chǎn)生不同的認識。目前, 對河流健康評價的研究大多數(shù)以單條河流為主, 因此, 進行河段–河流–流域的多尺度健康評價是未來研究的發(fā)展方向。
人類活動是損害河流生態(tài)系統(tǒng)健康的主要原因, 基于河流健康評價模型識別影響河流健康的關鍵因素, 量化人類活動對河流生態(tài)系統(tǒng)健康的貢獻,揭示流域內人類活動影響河流生態(tài)健康的機制, 進而采取針對性的河流保護和修復措施, 也是未來研究的重要內容。
本研究提出的指標體系是一個開放的體系, 在實際應用時, 可根據(jù)河流功能、管理需求和人類活動方式等, 對評價指標和權重進行適當?shù)恼{整與完善。為全面掌握河流健康的變化, 進一步探索河流健康的驅動機制, 需要開展河流全物質通量監(jiān)測和大數(shù)據(jù)庫建設, 進行長時間序列的河流系統(tǒng)健康動態(tài)評價。