曹秉帥, 單 楠, 顧羊羊, 敖 文, 龐 博, 竇華山, 王文林, 鄒長新*
1.生態(tài)環(huán)境部南京環(huán)境科學研究所, 江蘇 南京 210042 2.呼倫貝爾市北方寒冷干旱地區(qū)內(nèi)陸湖泊研究院, 內(nèi)蒙古 呼倫貝爾 021000 3.生態(tài)環(huán)境部南京環(huán)境科學研究所呼倫湖研究中心, 江蘇 南京 210042
呼倫湖是我國北方第一大湖泊,也是我國第五大湖泊,歷史上水域面積曾達到 2 339 km2,具有涵養(yǎng)水源、生物多樣性維護、氣候調(diào)節(jié)等重要生態(tài)系統(tǒng)功能[1],且其生態(tài)系統(tǒng)所能提供的生態(tài)價值遠大于經(jīng)濟價值[2]. 然而,近20年來,隨著氣候暖干化加劇,呼倫湖流域生態(tài)系統(tǒng)表現(xiàn)出湖泊濕地面積萎縮[3-4]、草地退化[5-9]、土地沙化[10-12]、魚類小型化[13]等生態(tài)問題,對我國北方生態(tài)安全屏障造成威脅,不容忽視. 因此,研究呼倫湖生態(tài)安全狀況,對于補齊短板、精準施策乃至維持區(qū)域可持續(xù)發(fā)展具有重要意義.
生態(tài)安全評價是生態(tài)安全研究的基礎和前提[14],準確可靠的生態(tài)安全評價結(jié)果能夠反映生態(tài)環(huán)境所存在的安全問題及隱患[15]. 目前,國內(nèi)關(guān)于湖泊的生態(tài)安全評價研究多集中于太湖[16]、洞庭湖[17-18]、巢湖[19]、滇池[20]等南方受人類活動影響較大的湖泊,評價關(guān)注的重點也更多在水環(huán)境〔如化學需氧量(COD)、總磷(TP)、總氮(TN)等〕、水生態(tài)(如生物多樣性指數(shù)、水華指數(shù)等)以及社會經(jīng)濟類指標. 也有學者對東北湖區(qū)(山口湖流域)[21]、新疆博斯騰湖[22]等開展了生態(tài)安全評價研究,但評價指標依然采用南方湖泊評價指標體系,如水質(zhì)、人口、經(jīng)濟、社會指標等,過于強調(diào)人為因素影響. 然而,以呼倫湖為代表的北方寒冷干旱區(qū)內(nèi)陸湖泊,受人類活動影響有限,其水體生態(tài)環(huán)境質(zhì)量受背景值影響較大,水質(zhì)變化主要受自然因素影響,且不具備飲用水源等功能. 有研究[23]表明,呼倫湖水體COD與總有機碳濃度呈顯著正相關(guān),表明水體有機物以含碳有機物為主,含碳有機物主要為難降解腐殖質(zhì),主要來源于本地及河流輸入草原難分解有機物. 因此,在評價這一類湖泊的生態(tài)安全狀況時應“陸-湖”統(tǒng)籌,充分考慮陸域生態(tài)要素的變化對湖泊、河流生態(tài)環(huán)境的影響,而不僅是關(guān)注水體本身.
該文以“山水林田湖草沙”系統(tǒng)觀為指導思想,通過壓力-狀態(tài)-響應(PSR)模型建立涵蓋水資源、水環(huán)境、水生態(tài)及陸域生態(tài)等多要素的呼倫湖生態(tài)安全綜合評價指標體系,并基于遙感和GIS技術(shù)分別從湖泊和流域等多尺度對呼倫湖1990—2018年的生態(tài)安全格局進行評價,以期為科學指導呼倫湖綜合治理、保障呼倫湖生態(tài)安全底線提供技術(shù)支撐.
圖1 呼倫湖流域生態(tài)系統(tǒng)類型Fig.1 Ecosystem type map of Hulun Lake Basin
研究區(qū)位于位于內(nèi)蒙古自治區(qū)呼倫貝爾草原西部的呼倫湖流域,117°00′10″E~117°41′40″E、48°30′40″N~49°20′40″N,是我國中高緯度溫帶半干旱區(qū)數(shù)千公里內(nèi)唯一的大澤,其流域呈放射狀分布,包括呼倫湖、貝爾湖局部、烏爾遜河、克魯倫河、額爾古納河、海拉爾河、新開河、輝河、伊敏河等主要河流及周邊沼澤地,面積約7.9×104km2. 研究區(qū)內(nèi)生態(tài)系統(tǒng)類型以草原、森林和濕地為主,三者分別占研究區(qū)總面積的65%、19%和9%.
圖2 呼倫湖流域生態(tài)安全影響因子及其作用機制示意[24]Fig.2 Schematic diagram of ecological security influencing factors and their mechanism in Hulun Lake Basin[24]
1.2.1指標篩選原則及體系構(gòu)建
筆者通過前期調(diào)查研究初步概括了呼倫湖流域各生態(tài)安全影響因子及其相互作用(見圖2)[24],并對呼倫湖生態(tài)安全問題及其成因進行了全面診斷,研究發(fā)現(xiàn):①2000年以來,呼倫湖流域氣候暖干化趨勢加劇,在長期超載過牧等人類活動壓力下,草地退化顯著,土地沙化、鹽堿化面積擴大;②受氣候暖干化和人類活動的交織作用,濕地面積萎縮,水生態(tài)系統(tǒng)紊亂,藍藻水華頻發(fā)[24]. 在此基礎上,結(jié)合指標構(gòu)建的可操作性、系統(tǒng)性、科學性、針對性以及數(shù)據(jù)的可獲得性等原則,該研究采用通用的壓力-狀態(tài)-響應(Pressure-State-Response,PSR)模型[25],以氣候條件、人口、水資源及放牧強度等因素構(gòu)成壓力子系統(tǒng);以水環(huán)境、水生態(tài)、陸域生態(tài)、生態(tài)系統(tǒng)服務功能與敏感性以及景觀格局等因素組成狀態(tài)子系統(tǒng),以生態(tài)恢復措施和生態(tài)治理投入作為響應子系統(tǒng),共同組建呼倫湖生態(tài)安全評價指標體系(見表1).
1.2.2指標權(quán)重確定原則與方法
層次分析法作為指標賦權(quán)的常用方法,在處理難以定量分析的復雜問題時能夠充分利用專家經(jīng)驗,具有一定權(quán)威性[31],已被廣泛應用于生態(tài)安全評價指標權(quán)重的計算中. 該研究緊密圍繞呼倫湖水生態(tài)-水環(huán)境-水資源“三水”共生目標,基于生態(tài)安全問題及其成因分析結(jié)果,結(jié)合專家先驗知識,采用專家打分法與層次分析法相結(jié)合的方式計算指標權(quán)重,結(jié)果如表1所示.
表1 呼倫湖流域生態(tài)安全評價指標體系及指標權(quán)重和標準值Table 1 Index system, index weight and standard value of ecological security assessment in Hulun Lake Basin
1.2.3指標標準值(閾值)確定原則
為了能夠突出呼倫湖的地域性特征,該文指標標準值的獲取原則是:①優(yōu)先采用實地調(diào)查法;②有國家標準的實行國家標準;③沒有國家標準的實行地方標準;④采用參考文獻中使用的標準;⑤采用往年平均值或相似區(qū)域的平均值;⑥參考類似發(fā)達地區(qū)的水平標準或采用科學研究已判定的標準. 指標體系標準值及其依據(jù)如表1所示.
1.2.3.1氣候條件
呼倫湖流域草地植被的生長受降水、溫度與人類活動的共同影響. 干旱氣候條件下牧草產(chǎn)量不足豐水年的70%,且隨著干旱強度和頻次的增加,土地沙化、裸露逐漸加重,草地植被有進一步退化的風險. 而降水少、氣溫高是干旱氣候條件的典型特征,且隨干旱程度的加大而更加明顯. 通過部門訪談與文獻資料調(diào)研[8,32-34]得知,呼倫湖歷史上高強度、高頻次干旱發(fā)生年份為1987年、1994年、2000—2012年,而受氣候暖干化影響,呼倫湖降水量2000—2012年急劇下降(見圖3). 將歷史上高頻干旱發(fā)生年份與降雨量急劇下降的年份進行疊加,最終確定2000—2012年的氣候條件為對呼倫湖生態(tài)安全影響最大的時期. 因此,該研究將年降雨量、年均氣溫和年蒸發(fā)量在2000—2012年間的平均值作為其各自標準值.
1.2.3.2水資源
圖3 1981—2014年呼倫湖流域降水量和年均氣溫的變化Fig.3 The variation of precipitation and temperature in the Hulun Lake Basin from 1981 to 2014
有研究表明生態(tài)水面法能夠較為客觀地確定濕地最小生態(tài)水面面積[27]. 生態(tài)水面法認為,對于高頻率出現(xiàn)的水面面積,是濕地適應自然的穩(wěn)態(tài)表現(xiàn),一般其生態(tài)系統(tǒng)健康狀況較為良好;而濕地最小生態(tài)水面面積是維持濕地自身存在的水面積,必然出現(xiàn)在低頻率的年份中[27]. 該文通過分析呼倫湖長時間序列的水面面積數(shù)據(jù),得到不同水面面積的頻率直方圖(見圖4). 可以看出,近60年呼倫湖湖面面積低于 2 036 km2的次數(shù)占比不足30%,即70%以上的時間湖面面積都維持在 2 036 km2以上,如果時間尺度擴大,該比例會更高. 因此,以 2 036 km2水面面積作為該指標的標準值. 采用同樣的方法,確定入湖徑流量指標的標準值為3.1×108m3.
圖4 1961—2019年呼倫湖湖面面積頻率直方圖Fig.4 Frequency histogram of Hulun Lake area from 1961 to 2019
1.3.1氣象數(shù)據(jù)
該文使用的降雨量和氣溫數(shù)據(jù)來自英國東英格利亞大學氣候研究所(Climatic Research Unit, CRU)全球氣象數(shù)據(jù)產(chǎn)品,該產(chǎn)品是基于中國氣象局和其他外部代理商每天或次日數(shù)據(jù)計算得出的每月觀測數(shù)據(jù)的每月網(wǎng)格字段,數(shù)據(jù)分辨率為0.5°×0.5°. 該研究通過獲取1990年、2000年、2010年、2015年和2018年5個時間節(jié)點的流域逐月降水量和平均氣溫,在ArcGIS中使用柵格計算器得到逐年氣象數(shù)據(jù)(年降水量、年均氣溫),采用IDW反距離權(quán)重插值重采樣成250 m×250 m分辨率,并根據(jù)呼倫湖流域氣象站點實測數(shù)據(jù)進行驗證. 該文使用的蒸發(fā)量數(shù)據(jù)采用中國氣象科學數(shù)據(jù)共享服務網(wǎng)(http://data.cma.cn)提供的1990年、2000年、2010年、2015年和2018年滿洲里(117°26′E、49°34′N)、新巴爾虎右旗(116°49′E、48°40′N)及新巴爾虎左旗(118°16′E、48°13′N)等氣象站點的月均值.
1.3.2水資源、水環(huán)境、水生態(tài)數(shù)據(jù)
該文使用的水資源數(shù)據(jù)(包括湖泊面積和入湖量)主要來自呼倫湖水文站1961—2019年的長期監(jiān)測數(shù)據(jù),由當?shù)厮块T提供;水環(huán)境數(shù)據(jù)(包括COD、TP、TN濃度)來源于呼倫湖流域小河口、甘珠花、拴馬樁等18個“十三五”國控重點流域河流水質(zhì)監(jiān)測斷面的長時間序列(1987—2018年)水質(zhì)監(jiān)測數(shù)據(jù);水生態(tài)指標中藍藻水華面積參考HJ 1098—2020《水華遙感與地面監(jiān)測評價技術(shù)規(guī)范》計算,水鳥數(shù)量歷年數(shù)據(jù)參照《生物多樣性觀測技術(shù)導則 鳥類》標準,分別在小河口、五一、拴馬樁、嘎拉達白辛、烏蘭諾爾、烏都魯、甘珠花等地觀測獲得,浮游動物生物量數(shù)據(jù)根據(jù)《淡水浮游生物研究方法》[35]采用25#的浮游生物網(wǎng)(網(wǎng)孔直徑0.064 mm)采集獲取.
1.3.3放牧強度數(shù)據(jù)
草地載畜量數(shù)據(jù)采用《呼倫貝爾市統(tǒng)計年鑒2020》[36]中的呼倫貝爾市歷年(1949—2017年)大牲畜和羊的統(tǒng)計數(shù)據(jù),并基于NY/T 635—2015《天然草地合理載畜量的計算》中的成年家畜折算系數(shù)將牲畜統(tǒng)計數(shù)據(jù)統(tǒng)一折算為標準家畜單位,即綿羊單位. 草地載畜強度是指單位面積的羊單位數(shù)量.
1.3.4陸域生態(tài)數(shù)據(jù)
土地沙化與鹽堿化面積是利用1990年、2000年、2010年、2015年和2018年的高分一號、高分二號等高空間分辨率衛(wèi)星遙感影像以及Landsat TM/ETM等中分辨率遙感影像,采用人機交互和面向?qū)ο蠓诸惙椒?,通過地面核查進行復核、完善提取得到,空間分辨率達5 m. 植被覆蓋度和地上部生物量數(shù)據(jù)來源于MODIS植被指數(shù)產(chǎn)品中的MOD13Q1產(chǎn)品,空間分辨率為250 m,時間分辨率為16 d,時間范圍為1990年、2000年、2010年、2015年和2018年.
1.3.5景觀格局數(shù)據(jù)
該文參考陳昆侖等[37]通過計算景觀指數(shù)來概括景觀格局的方法,選取景觀類型斑塊密度、平均斑塊形狀指數(shù)、聚集度指數(shù)和斑塊所占景觀面積比例為主要指標,建立景觀指數(shù)與景觀生態(tài)安全度的定量化描述. 針對呼倫湖區(qū)域1990年、2000年、2010年、2015年和2018年共5個時期的Landsat TM/ETM中等分辨率衛(wèi)星遙感影像,通過影像校正、影像裁剪、影像拼接等數(shù)據(jù)預處理后,經(jīng)ENVI軟件計算分別獲取草地、濕地、林地等景觀類型的斑塊密度、平均斑塊形狀指數(shù)、聚集度指數(shù)和斑塊所占景觀面積比例.
1.4.1評價尺度
為了能夠充分反映呼倫湖流域生態(tài)安全在時空尺度上的變化特征,考慮到指標數(shù)據(jù)的易獲得性和可靠性,該研究以1990年、2000年、2010年、2015年和2018年為評價時段,空間分辨率為250 m.
1.4.2指標歸一化方法
多指標體系中,由于各指標的單位、量綱、數(shù)量級不同,不具備可比性,為了消除這一差異,應對指標進行無量綱處理. 該文運用閾值標準化法對指標原始數(shù)據(jù)進行標準化處理.
對于正向型指標(值越大越安全的指標),若以安全值為標準值:
(1)
若以不安全值為標準值:
(2)
對于負向型指標(值越小越安全的指標),若以安全值為標準值:
(3)
若以不安全值為標準值:
(4)
式中,Xij和Xij′分別為第i年第j項指標的原始值和標準化指,Xb為第j項指標的標準值.
1.4.3生態(tài)安全評價綜合指數(shù)計算
指標無量綱化后,采用式(5)計算第i年生態(tài)安全指數(shù)(ESI):
(5)
式中,Xj為第j個指標的標準化值,wj為第j個指標權(quán)重,n為指標個數(shù).
1.4.4生態(tài)安全等級劃分
該文采用均分法將生態(tài)安全狀態(tài)等比例劃分為很差、較差、預警、良好、理想等5個等級,生態(tài)安全指數(shù)得分越高,表明生態(tài)安全狀況越好,具體分級標準如表2所示.
表2 生態(tài)安全綜合評估分級標準Table 2 Classification standard of comprehensive ecological security assessment
呼倫湖及其流域1990—2018年生態(tài)安全指數(shù)平均值如圖5所示. 結(jié)果顯示,1990—2018年,呼倫湖與其流域生態(tài)安全變化趨勢均呈現(xiàn)先降低后逐步向好的態(tài)勢,且流域的生態(tài)安全水平提升更加明顯,2018年其生態(tài)安全指數(shù)分別為0.495和0.774,分別處于預警和良好狀態(tài). 根據(jù)不同時期呼倫湖及其流域生態(tài)安全狀況,可進一步將其生態(tài)安全變化趨勢劃分為3個時期:①1990—2000年為驟降期,呼倫湖生態(tài)安全指數(shù)由1990年的0.548降至2000年的0.472,降幅達到13.9%,而其流域生態(tài)安全指數(shù)也下降了2.2%;②2000—2015年為穩(wěn)定期,這一時期呼倫湖及其流域生態(tài)安全狀況基本保持穩(wěn)定,生態(tài)安全指數(shù)分別維持在0.470和0.718左右;③2015—2018年為恢復期,呼倫湖生態(tài)安全指數(shù)由0.478提升至0.495,增幅約3.6%,而其流域生態(tài)安全水平更是提高了8.1%,恢復態(tài)勢明顯.
圖5 1990—2018年呼倫湖及其流域范圍 生態(tài)安全指數(shù)平均值Fig.5 The mean value of ecological security index of Hulun Lake and its watershed from 1990 to 2018
圖6 2018年呼倫湖流域生態(tài)安全格局Fig.6 The ecological security pattern of the Hulun Lake Basin in 2018
根據(jù)生態(tài)安全等級劃分標準,將2018年呼倫湖流域生態(tài)安全水平進行等級劃分,結(jié)果如圖6所示. 從圖6可以看出,呼倫湖流域生態(tài)安全呈“一林一草一湖”分布格局,即東部大興安嶺森林生態(tài)功能區(qū)—中部呼倫貝爾草原草甸生態(tài)功能區(qū)—西部呼倫湖湖區(qū)等“三區(qū)”分布格局. 從生態(tài)安全等級來看,東部林區(qū)生態(tài)安全水平最好,為Ⅰ級理想狀態(tài),占評價區(qū)面積的33.7%,生態(tài)系統(tǒng)服務功能和結(jié)構(gòu)完整,生態(tài)環(huán)境基本沒有受到破壞;中部草原區(qū)生態(tài)安全水平適中,為Ⅱ級良好狀態(tài),占評價區(qū)面積的63.5%;西部湖區(qū)生態(tài)安全水平最低,為Ⅲ級預警狀態(tài),僅占評價區(qū)面積的2.8%,生態(tài)系統(tǒng)服務功能已有退化,生態(tài)系統(tǒng)結(jié)構(gòu)有變化,但尚可維持基本功能,生態(tài)問題顯現(xiàn),受干擾后易惡化.
呼倫湖1990—2018年生態(tài)安全時空變化如圖7所示. 由圖7可見,呼倫湖生態(tài)安全經(jīng)歷了驟降(1990—2000年)、穩(wěn)定(2000—2015年)和恢復(2015—2018年)等3個時期. 其中,最低生態(tài)安全水平發(fā)生在2010年,95.8%的湖體生態(tài)安全指數(shù)低于0.48,較差區(qū)域分布在新開河入湖口、湖西岸大部分區(qū)域、湖中心以及湖東南方向灣口周圍;而較好區(qū)域主要分布在烏爾遜河入湖口. 2010—2015年,呼倫湖生態(tài)安全狀況有所恢復,但沿湖濱岸帶以及呼倫溝周圍生態(tài)安全水平仍較低. 到2018年,全湖區(qū)生態(tài)安全已有明顯回升,生態(tài)安全指數(shù)基本維持在0.48~0.54之間,僅有東南部灣區(qū)少量區(qū)域(不足3.1%)生態(tài)安全水平(<0.48)較低.
呼倫湖生態(tài)安全的上述變化規(guī)律在其流域范圍內(nèi)則表現(xiàn)得更為明顯,如圖8所示,根據(jù)生態(tài)安全變化分布格局可將呼倫湖流域分為沿海拉爾河上游—免渡河流域的東部林區(qū)、新巴爾虎左旗境內(nèi)的中部草原區(qū)(沙帶)及位于新巴爾虎右旗境內(nèi)沿克魯倫河兩岸的西部草原區(qū)等3個區(qū)域. 其中,東部林區(qū)生態(tài)安全水平較高,基本維持在0.74以上;西部草原區(qū)次之,絕大部分區(qū)域生態(tài)安全指數(shù)在0.70以上,少部分低于0.70的區(qū)域零星分布在沿克魯倫河兩岸;而生態(tài)安全水平下降最為嚴重的區(qū)域主要位于新左旗中部、海拉爾河流域以及呼倫溝等地,尤其是沿烏爾遜河上游東側(cè)地帶,生態(tài)安全指數(shù)低于0.68.
圖8 1990—2018年呼倫湖流域生態(tài)安全時空分布格局變化趨勢Fig.8 The spatial and temporal distribution pattern of ecological security in Hulun Lake Basin from 1990 to 2018
為指導完善呼倫湖生態(tài)安全體系建設,找準“痛點”問題,挖掘呼倫湖流域生態(tài)安全水平提升途徑,計算得到1990—2018年呼倫湖流域生態(tài)安全影響要素的變化趨勢如圖9所示. 從圖9可以看出:①水資源、水生態(tài)和生態(tài)系統(tǒng)服務功能與敏感性的生態(tài)安全變化趨勢與流域生態(tài)安全變化趨勢相似,均呈現(xiàn)出先驟降、再穩(wěn)定、后恢復的變化趨勢. 通過進一步分析指標層生態(tài)安全指數(shù)變化可知,對上述3項要素生態(tài)安全影響最大的指標分別為入湖徑流量、藍藻水華面積占比和水源涵養(yǎng)量,說明水資源保障與藍藻水華治理對提升呼倫湖生態(tài)安全水平具有重要作用. ②放牧強度、水環(huán)境和景觀格局生態(tài)安全呈下降趨勢,亟待采取有效措施加以控制和恢復. 這主要與超載過牧、水質(zhì)惡化及濕地景觀生態(tài)安全度下降有關(guān). 據(jù)調(diào)查,2000—2010年,位于呼倫湖流域境內(nèi)的呼倫貝爾市牧業(yè)四旗(包括新巴爾虎左旗、新巴爾虎右旗、陳巴爾虎旗和鄂溫克旗),單位面積放牧強度增加了176%;呼倫湖COD濃度由1990年的81.6 mg/L升至2018年的102.4 mg/L,而2015年出現(xiàn)短暫好轉(zhuǎn)主要與當年入湖徑流量增大從而導致入湖污染物稀釋有關(guān)(2015年入湖徑流量比2018年多45.6%);1998—2011年,隨著水位的快速降低,濕地面積由140.72 km2萎縮到5.49 km2,隨著海拉爾河引水工程的實施,2013年濕地面積才恢復至21.58 km2[24].
圖9 1990—2018年呼倫湖流域生態(tài)安全影響要素變化趨勢Fig.9 The variation trend of ecological security impact factors in Hulun Lake Basin from 1990 to 2018
a) 近30年來,呼倫湖流域生態(tài)安全狀況經(jīng)歷了1990—2000年的驟降期、2000—2015年的穩(wěn)定期和2015—2018年的恢復期等3個時期,目前正處于穩(wěn)中向好的恢復期,呈現(xiàn)“一林一草一湖”的生態(tài)安全分布格局.
b) 呼倫湖生態(tài)安全水平最低值出現(xiàn)在2010年,95.8%的湖體生態(tài)安全指數(shù)低于0.48,主要分布在新開河入湖口、湖西岸大部分區(qū)域、湖中心以及湖東南方向灣口周圍;流域生態(tài)安全水平最低值出現(xiàn)2015年,較差區(qū)域主要位于新左旗中部、海拉爾河流域以及呼倫溝等地,尤其是沿烏爾遜河上游東側(cè)地帶.
c) 影響呼倫湖流域生態(tài)安全水平的主要因素為入湖徑流量、藍藻水華面積占比和水源涵養(yǎng)量,而長期超載過牧、水體污染物濃縮效應以及濕地面積萎縮是限制生態(tài)安全水平進一步提升的另一重要因素. 因此,通過保障水資源供給、建立藍藻水華風險防控體系,堅持以水定牧、以草定牧、草畜平衡,以資源環(huán)境承載力為基礎,合理核定載畜量,開展蘆葦濕地生態(tài)保護與修復等一系列措施,對于有效提升呼倫湖生態(tài)安全恢復潛力具有重要作用.