張夢(mèng)宇,張 黎,何洪林,任小麗,牛忠恩,呂 妍,徐 茜,常清青,劉衛(wèi)華,李 攀

1 中國(guó)科學(xué)院地理科學(xué)與資源研究所 生態(tài)系統(tǒng)網(wǎng)絡(luò)觀測(cè)與模擬重點(diǎn)實(shí)驗(yàn)室,北京 100101 2 國(guó)家生態(tài)科學(xué)數(shù)據(jù)中心,北京 100101 3 中國(guó)科學(xué)院大學(xué),北京 100049 4 中國(guó)科學(xué)院大學(xué)資源與環(huán)境學(xué)院,北京 100190 5 天津大學(xué) 地球系統(tǒng)科學(xué)學(xué)院,天津 300072

生態(tài)系統(tǒng)質(zhì)量是診斷生態(tài)系統(tǒng)是否健康的指標(biāo),是指在一個(gè)具體的時(shí)間和空間范圍內(nèi)生態(tài)系統(tǒng)的總體或部分生命組分的質(zhì)量,主要表現(xiàn)在其生產(chǎn)能力、自我維持能力、受到外界干擾后的抗干擾能力以及對(duì)人類(lèi)的生存及社會(huì)經(jīng)濟(jì)持續(xù)發(fā)展的影響[1-4]。一個(gè)高質(zhì)量的生態(tài)系統(tǒng)能夠在壓力下維持自身組織結(jié)構(gòu)和穩(wěn)定,具有自我恢復(fù)的能力[1]。改革開(kāi)放后中國(guó)人口與經(jīng)濟(jì)快速增長(zhǎng),耕地不斷擴(kuò)張[5],森林和草地等生態(tài)系統(tǒng)發(fā)生了退化[6-7],帶來(lái)自然棲息地破壞和生態(tài)系統(tǒng)功能下降等生態(tài)問(wèn)題[8-10]。但自1999年以來(lái)中國(guó)實(shí)施了草地保護(hù)、天然林保護(hù)、退耕還林以及石漠化綜合治理等工程,退化的生態(tài)系統(tǒng)得到一定的恢復(fù)[11-14]。黨中央在十九大報(bào)告中將生態(tài)文明提升到戰(zhàn)略地位,要通過(guò)加大生態(tài)系統(tǒng)保護(hù)力度,實(shí)施重要生態(tài)系統(tǒng)保護(hù)和修復(fù)重大工程,以提升生態(tài)系統(tǒng)質(zhì)量和穩(wěn)定性。在當(dāng)前我國(guó)土地利用需求不斷增加和生態(tài)工程大規(guī)模實(shí)施的背景下,迫切需要對(duì)我國(guó)生態(tài)系統(tǒng)質(zhì)量狀況進(jìn)行評(píng)估,為土地利用和生態(tài)系統(tǒng)可持續(xù)管理提供參考[15]。

我國(guó)已開(kāi)展一系列與生態(tài)系統(tǒng)質(zhì)量有關(guān)的評(píng)估工作,主要包括生態(tài)系統(tǒng)功能評(píng)估[16-17]、生態(tài)系統(tǒng)健康評(píng)估[18]和生態(tài)環(huán)境質(zhì)量評(píng)估[19]等方面,為國(guó)家和地方生態(tài)保護(hù)政策的制定提供了科學(xué)依據(jù)。在全國(guó)層面,根據(jù)“全國(guó)生態(tài)環(huán)境十年變化(2000—2010年)調(diào)查評(píng)估”項(xiàng)目研究結(jié)果,2000—2010年間食物生產(chǎn)、水源涵養(yǎng)、土壤保持、防風(fēng)固沙、洪水調(diào)蓄、固碳這6項(xiàng)生態(tài)系統(tǒng)服務(wù)得到明顯改善,而生物多樣性維護(hù)功能有所下降[17]；中國(guó)生態(tài)系統(tǒng)健康水平呈現(xiàn)從西北到東南上升的格局,水分指數(shù)和土地利用強(qiáng)度的貢獻(xiàn)分別為24.5%和20.7%[18]。我國(guó)生態(tài)環(huán)境部依據(jù)生態(tài)環(huán)境狀況評(píng)價(jià)技術(shù)規(guī)范(HJ 192—2015),定期發(fā)布《全國(guó)生態(tài)環(huán)境質(zhì)量報(bào)告》,公布全國(guó)省域和縣域生態(tài)環(huán)境狀況及變化評(píng)估結(jié)果。此外,國(guó)內(nèi)學(xué)者也開(kāi)展了一些典型區(qū)域的生態(tài)系統(tǒng)質(zhì)量評(píng)價(jià)方法和應(yīng)用研究。對(duì)于自然生態(tài)系統(tǒng),主要從生產(chǎn)能力、功能穩(wěn)定性和人口壓力[20],植被覆蓋度和葉面積指數(shù)[21],生態(tài)系統(tǒng)生產(chǎn)能力、穩(wěn)定性和承載力[22-23]等方面評(píng)估生態(tài)系統(tǒng)質(zhì)量。對(duì)于城市生態(tài)系統(tǒng),重點(diǎn)從經(jīng)濟(jì)、自然和社會(huì)環(huán)境質(zhì)量[24],結(jié)構(gòu)、功能和脅迫[25]等角度進(jìn)行評(píng)估。對(duì)于干旱區(qū)流域水文生態(tài)系統(tǒng),綜合植被覆蓋度、墾區(qū)土壤鹽漬化率、沙漠化面積比率、水土流失強(qiáng)度和生態(tài)用水比率等因素建立生態(tài)系統(tǒng)質(zhì)量框架[26]。計(jì)算生態(tài)系統(tǒng)質(zhì)量或健康綜合指數(shù)的方法主要包括求積法[27-28]、PSR模型[2,29]、主成分分析法[30]、層次分析法[31]等。然而,以往有關(guān)生態(tài)系統(tǒng)功能、健康和質(zhì)量的評(píng)估往往聚焦于生態(tài)系統(tǒng)指標(biāo)的現(xiàn)實(shí)值,其數(shù)值大小在很大程度上受氣候因子(如降水、溫度、總輻射)和地理本底條件的差異主導(dǎo)[32-33],導(dǎo)致不同地區(qū)評(píng)價(jià)結(jié)果的可比性較低,并且對(duì)于全國(guó)范圍生態(tài)系統(tǒng)質(zhì)量及其變化的認(rèn)識(shí)仍較為缺乏,尚不足以滿(mǎn)足國(guó)家生態(tài)系統(tǒng)管理的需求[15]。

根據(jù)生態(tài)系統(tǒng)完整性[4,34-36]和穩(wěn)定性[37]理論,生態(tài)系統(tǒng)質(zhì)量評(píng)估不僅應(yīng)綜合考慮一定時(shí)間、空間范圍內(nèi)生態(tài)系統(tǒng)結(jié)構(gòu)、功能、穩(wěn)定性特征,以及人類(lèi)脅迫或干擾狀況,從而反映生態(tài)系統(tǒng)的整體質(zhì)量狀態(tài)；還應(yīng)分區(qū)、分類(lèi)型地開(kāi)展評(píng)估參照生態(tài)系統(tǒng)(簡(jiǎn)稱(chēng)參照系)的設(shè)定,使得評(píng)估結(jié)果在不同時(shí)間和空間上具有較好的可比性[15]。由于現(xiàn)實(shí)世界中很難找到完全不受人類(lèi)影響的理想?yún)⒄障?將與待評(píng)估生態(tài)系統(tǒng)具有相似植被、土壤、氣候條件且受人類(lèi)活動(dòng)影響相對(duì)較小的自然生態(tài)系統(tǒng)作為參照系具有較好的合理性與可行性。參照生態(tài)系統(tǒng)主要包括：(1)極小干擾參照狀態(tài),生態(tài)系統(tǒng)最大程度地保留了完整性,幾乎沒(méi)有人類(lèi)活動(dòng)；(2)歷史條件參照狀態(tài),在人類(lèi)干擾之前的生態(tài)系統(tǒng)狀態(tài)；(3)最佳可達(dá)狀態(tài),在某一特定地理和氣候本底條件下,生態(tài)系統(tǒng)所能達(dá)到的最佳狀態(tài),可以根據(jù)分位數(shù)法選取同一條件區(qū)域內(nèi)一組參照點(diǎn)作為參照系[38-40]。針對(duì)我國(guó)陸地生態(tài)系統(tǒng)特點(diǎn)和數(shù)據(jù)可獲取性狀況,將我國(guó)國(guó)家自然保護(hù)區(qū)核心區(qū)、中國(guó)生態(tài)系統(tǒng)研究網(wǎng)絡(luò)(CERN)野外觀測(cè)臺(tái)站永久觀測(cè)樣地等較少受人類(lèi)活動(dòng)干擾且狀態(tài)優(yōu)良的生態(tài)系統(tǒng)作為參照系,是當(dāng)前我國(guó)陸地生態(tài)系統(tǒng)質(zhì)量評(píng)估的一種新思路[15]。

本文利用基于參照系的生態(tài)系統(tǒng)質(zhì)量評(píng)估方法,以中國(guó)陸地生態(tài)系統(tǒng)為研究對(duì)象,綜合土地覆蓋數(shù)據(jù),降水、氣溫、總輻射的年度氣候數(shù)據(jù)和生態(tài)系統(tǒng)生產(chǎn)力、固碳、水源涵養(yǎng)和土壤保持四項(xiàng)生態(tài)系統(tǒng)功能數(shù)據(jù),評(píng)估1990、2000、2010、2015年四期生態(tài)系統(tǒng)質(zhì)量,重點(diǎn)分析參照生態(tài)系統(tǒng)的質(zhì)量變化特征,基于參照系的中國(guó)陸地生態(tài)系統(tǒng)質(zhì)量的空間格局,以及生態(tài)工程大規(guī)模實(shí)施前后的兩個(gè)時(shí)段(即1990—2000年和2000—2015年)中國(guó)陸地生態(tài)系統(tǒng)質(zhì)量變化特征及其歸因,以期為我國(guó)生態(tài)工程的成效評(píng)估提供依據(jù),為土地利用和生態(tài)系統(tǒng)管理提供支撐。

1 數(shù)據(jù)和方法

1.1 數(shù)據(jù)

1.1.1土地覆蓋數(shù)據(jù)

采用1990年、2000年、2010年、2015年空間分辨率為1 km的土地覆蓋數(shù)據(jù)(ChinaCover)[41]計(jì)算生態(tài)系統(tǒng)結(jié)構(gòu)指數(shù)。該數(shù)據(jù)是基于國(guó)產(chǎn)環(huán)境災(zāi)害衛(wèi)星(HJ-1A/B)和美國(guó)陸地衛(wèi)星(Landsat OLI)數(shù)據(jù),采用面向?qū)ο蟮亩喑叨确指钆c變化檢測(cè)的分類(lèi)方法進(jìn)行制圖,一級(jí)地類(lèi)的精度達(dá)94%。

1.1.2環(huán)境要素?cái)?shù)據(jù)

采用1990—2015年1 km分辨率的氣候要素(年均溫、年降水、年總輻射)數(shù)據(jù)分析氣候因子對(duì)生態(tài)系統(tǒng)質(zhì)量變化的影響。其中,氣溫和降水?dāng)?shù)據(jù)基于中國(guó)地面氣候日值數(shù)據(jù)集和全球歷史氣候網(wǎng)格數(shù)據(jù)集采用ANUSPLIN軟件插值計(jì)算生成[42],總輻射數(shù)據(jù)由日照時(shí)數(shù)計(jì)算生成[43]。根據(jù)《1∶1000000 中國(guó)植被圖集》和中國(guó)生態(tài)地理分區(qū)生成168個(gè)評(píng)估單元,數(shù)據(jù)源自中國(guó)科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心(http://www.resdc.cn/)。

1.1.3生態(tài)系統(tǒng)功能數(shù)據(jù)

1990—2015年生態(tài)系統(tǒng)功能指標(biāo)數(shù)據(jù)由生態(tài)系統(tǒng)過(guò)程模型(Carbon Exchange between Vegetation,Soil,and Atmosphere—Ecosystem services,CEVSA-ES[44])模擬得到。CEVSA-ES模型模擬的NPP、NEP、土壤含水量和土壤保持量與觀測(cè)數(shù)據(jù)具有較好的一致性,二者線(xiàn)性擬合的決定系數(shù)分別為0.96、0.81、0.47、0.52[44]。模擬結(jié)果的空間分辨率為0.05°。

1.1.4其他數(shù)據(jù)

1990—2015年GLASS LAI數(shù)據(jù)[45-46]用于分析植被生長(zhǎng)狀況,數(shù)據(jù)源自在國(guó)家地球系統(tǒng)科學(xué)數(shù)據(jù)中心,其空間分辨率為1 km,時(shí)間分辨率為8 d。本文將其計(jì)算為年最大LAI(LAImax)。國(guó)家自然保護(hù)區(qū)空間分布數(shù)據(jù)(截至2013年9月)用于選取參照系,數(shù)據(jù)來(lái)自于中國(guó)科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心(http://www.resdc.cn/)。CERN野外觀測(cè)臺(tái)站地理位置數(shù)據(jù)從國(guó)家生態(tài)科學(xué)數(shù)據(jù)中心(http://rs.cern.ac.cn/data/initDRsearch?cid=STA_DP2011)獲取。采用1990—2015年1 km分辨率的植被覆蓋度數(shù)據(jù)計(jì)算生態(tài)系統(tǒng)結(jié)構(gòu)指數(shù),該數(shù)據(jù)是以基于MODIS NDVI數(shù)據(jù),利用像元二分模型計(jì)算得到,在計(jì)算過(guò)程中利用ChinaCover產(chǎn)品對(duì)純土壤、植被像元的NDVI進(jìn)行了改進(jìn)[41]。

1.2 基于參照系的國(guó)家生態(tài)系統(tǒng)質(zhì)量評(píng)價(jià)方法

根據(jù)生態(tài)系統(tǒng)完整性和穩(wěn)定性理論[4,34-37],本文從生態(tài)系統(tǒng)結(jié)構(gòu)、功能、穩(wěn)定性和人類(lèi)脅迫四個(gè)方面構(gòu)建生態(tài)系統(tǒng)質(zhì)量評(píng)價(jià)指標(biāo)(表1),將生態(tài)系統(tǒng)質(zhì)量指數(shù)(Ecosystem quality Index,ESQI)表示為生態(tài)系結(jié)構(gòu)指數(shù)、功能指數(shù)、穩(wěn)定性指數(shù)統(tǒng)和脅迫指數(shù)的加權(quán)和(公式1)。ESQI高表示人類(lèi)脅迫小,生態(tài)系統(tǒng)景觀結(jié)構(gòu)優(yōu)良,并具有良好且穩(wěn)定的功能。

表1 國(guó)家生態(tài)系統(tǒng)質(zhì)量指標(biāo)體系Table 1 Indicator framework for national ecosystem quality

(1)

式中,xi為一級(jí)指標(biāo),包括標(biāo)準(zhǔn)化結(jié)構(gòu)指數(shù)、標(biāo)準(zhǔn)化功能指數(shù)、標(biāo)準(zhǔn)化穩(wěn)定性指數(shù)和標(biāo)準(zhǔn)化脅迫指數(shù)。ωi為對(duì)應(yīng)的權(quán)重,根據(jù)熵權(quán)法計(jì)算[47-48],分別為0.299、0.335、0.317和0.048。

生態(tài)系統(tǒng)結(jié)構(gòu)指標(biāo)用生態(tài)系統(tǒng)結(jié)構(gòu)指數(shù)(Ecosystem Structure Index,EStrI)表示,由不同土地覆蓋類(lèi)型覆蓋率的加權(quán)和表示,參照生態(tài)環(huán)境部2015年發(fā)布的生態(tài)環(huán)境狀況評(píng)價(jià)技術(shù)規(guī)范,森林、草地、水域濕地、耕地、建設(shè)用地、未利用地的權(quán)重分別設(shè)為0.35、0.21、0.28、0.11、0.04和0.01。

(2)

式中,flc為土地覆蓋類(lèi)型lc的覆蓋率,ωlc為單位面積的土地覆蓋類(lèi)型lc的權(quán)重。

生態(tài)系統(tǒng)功能指標(biāo)由生態(tài)系統(tǒng)功能指數(shù)(Ecosystem Function Index,EFunI)表示,設(shè)定為四項(xiàng)重要生態(tài)系統(tǒng)功能(包括物質(zhì)生產(chǎn)、固碳、水源涵養(yǎng)和土壤保持功能)指標(biāo)的加權(quán)和。分別采用凈初級(jí)生產(chǎn)力(Net primary productivity,NPP)、凈生態(tài)系統(tǒng)生產(chǎn)力(Net ecosystem productivity,NEP)、土壤含水量及土壤保持量表征物質(zhì)生產(chǎn)、固碳、水源涵養(yǎng)和土壤保持功能。NPP、NEP和土壤含水量采用min—max標(biāo)準(zhǔn)化到0—1之間。由于土壤保持量呈現(xiàn)明顯的右偏態(tài)分布,本文對(duì)其數(shù)值取對(duì)數(shù)后再采用min—max標(biāo)準(zhǔn)化到0—1之間。最后利用熵權(quán)法確定四個(gè)功能指標(biāo)的權(quán)重。

生態(tài)系統(tǒng)穩(wěn)定性指數(shù)(Ecosystem Stability Index,EStaI)用待評(píng)估時(shí)間點(diǎn)的前十年NPP的時(shí)間穩(wěn)定性指數(shù)表示。該指數(shù)假設(shè)生態(tài)系統(tǒng)越接近長(zhǎng)時(shí)間序列的平均狀態(tài)就越穩(wěn)定,能較好地綜合反映生態(tài)系統(tǒng)的抵抗力和恢復(fù)力[49]?？紤]到中國(guó)大規(guī)模實(shí)施了生態(tài)工程并具有城市擴(kuò)張的背景,生態(tài)系統(tǒng)具有一定的恢復(fù)或退化的趨勢(shì),本研究采用去除時(shí)間趨勢(shì)后的生產(chǎn)力時(shí)間穩(wěn)定性(公式3)來(lái)表征中國(guó)生態(tài)系統(tǒng)的穩(wěn)定性指數(shù),并采用min—max方法將其標(biāo)準(zhǔn)化到0—1之間。

(3)

式中,meanfun為生態(tài)系統(tǒng)功能的多年平均值,stddetrend_fun為去趨勢(shì)后生態(tài)系統(tǒng)功能的標(biāo)準(zhǔn)差。

人類(lèi)脅迫指標(biāo)用人類(lèi)干擾指數(shù)(Human Interference Index,HII)表示,考慮每種土地覆蓋類(lèi)型的面積(覆蓋率)及其對(duì)生態(tài)系統(tǒng)的干擾強(qiáng)度[50-51],同時(shí)考慮周?chē)恋仡?lèi)型對(duì)其干擾的距離效應(yīng)[52]。采用min—max標(biāo)準(zhǔn)化方法將其負(fù)向標(biāo)準(zhǔn)化到0—1之間。

(4)

表2 土地覆蓋類(lèi)型的人類(lèi)干擾度范圍Table 2 The range of human interference index of different land cover types

為盡可能消除評(píng)價(jià)結(jié)果空間不可比問(wèn)題,依據(jù)每個(gè)評(píng)估單元內(nèi)較少受人類(lèi)活動(dòng)明顯干擾且優(yōu)良的參照系,計(jì)算基于參照系的生態(tài)系統(tǒng)質(zhì)量指數(shù)(ESQI′；公式5)。

ESQI′=ESQI/ESQIref

(5)

ESQIref為生態(tài)系統(tǒng)質(zhì)量的參照值,ESQI為生態(tài)系統(tǒng)質(zhì)量的現(xiàn)實(shí)值。當(dāng)ESQI′≥1時(shí),表示生態(tài)系統(tǒng)質(zhì)量?jī)?yōu)于參照系；反之生態(tài)系統(tǒng)質(zhì)量劣于參照系。根據(jù)分位數(shù)將ESQI′劃分為優(yōu)良(80%—100%)、較好(60%—80%)、中等(40%—60%)、較差(20%—40%)、很差(0%—20%)五個(gè)等級(jí)。生態(tài)系統(tǒng)質(zhì)量的參照值是指某一時(shí)期某個(gè)評(píng)估單元內(nèi)的參照系的質(zhì)量的平均值。本研究中針對(duì)每個(gè)具有相似地理和氣候條件的評(píng)估單元,按照下面方法的優(yōu)先次序篩選出評(píng)估單元內(nèi)適合作為參照生態(tài)系統(tǒng)的備選區(qū)域：(1)具有CERN野外觀測(cè)臺(tái)站或國(guó)家自然保護(hù)區(qū)核心區(qū)等較少受人類(lèi)影響的生態(tài)系統(tǒng)；(2)凈初級(jí)生產(chǎn)力主要受氣候因子控制的生態(tài)系統(tǒng)；(3)評(píng)估單元內(nèi)生態(tài)系統(tǒng)質(zhì)量指數(shù)為優(yōu)良的生態(tài)系統(tǒng)。在這些備選區(qū)域內(nèi)選擇未發(fā)生土地利用變化且生態(tài)系統(tǒng)質(zhì)量相對(duì)優(yōu)良的區(qū)域作為參照系,下文分別簡(jiǎn)稱(chēng)為CERN參照、國(guó)家自然保護(hù)區(qū)參照、最受氣候影響區(qū)參照和評(píng)估單元參照。

2 結(jié)果

2.1 參照系的生態(tài)系統(tǒng)質(zhì)量變化特征

對(duì)1990年、2000年、2010年、2015年不同類(lèi)型參照系的生態(tài)系統(tǒng)質(zhì)量指標(biāo)參照值進(jìn)行統(tǒng)計(jì),其占評(píng)估單元總體的分位數(shù)平均值如表3所示。總體上,參照系的生態(tài)系統(tǒng)質(zhì)量指數(shù)數(shù)值位于評(píng)估單元89%分位數(shù),其中脅迫標(biāo)準(zhǔn)化指數(shù)、結(jié)構(gòu)標(biāo)準(zhǔn)化指數(shù)、功能標(biāo)準(zhǔn)化指數(shù)、穩(wěn)定性標(biāo)準(zhǔn)化指數(shù)四個(gè)分量的參照值分別位于76%、80%、75%和69%。按照不同類(lèi)型的參照系來(lái)看,以國(guó)家自然保護(hù)區(qū)為參照的評(píng)估單元數(shù)最多,參照系的人類(lèi)干擾標(biāo)準(zhǔn)化指數(shù)位于評(píng)估單元的90.8%分位數(shù),在四類(lèi)參照系中的受干擾程度最低。以CERN永久樣地為基準(zhǔn)的參照生態(tài)系統(tǒng)具有相對(duì)較好的結(jié)構(gòu)、功能和穩(wěn)定性,均位于評(píng)估單元分位數(shù)的60%以上。對(duì)于最受氣候影響的參照系,其結(jié)構(gòu)指數(shù)所處于評(píng)估單元的分位數(shù)水平(77%)>功能指數(shù)(69.7%)>穩(wěn)定性指數(shù)(68.3%)。以評(píng)估單元優(yōu)良值為基準(zhǔn)的參照生態(tài)系統(tǒng)的質(zhì)量位于評(píng)估單元的91.6%分位數(shù)。可以看出,在具有相似氣候和地理?xiàng)l件的評(píng)估單元內(nèi),參照系具有較低的人類(lèi)脅迫和良好的生態(tài)系統(tǒng)結(jié)構(gòu)、功能和穩(wěn)定性,說(shuō)明本研究根據(jù)CERN生態(tài)系統(tǒng)永久樣地、國(guó)家自然保護(hù)區(qū)核心區(qū)、最受氣候影響區(qū)來(lái)選取受人類(lèi)干擾較少且優(yōu)良的生態(tài)系統(tǒng)作為參照系是合理可行的。

表3 不同類(lèi)型參照系的質(zhì)量參照值位于相應(yīng)評(píng)估單元分位數(shù)的平均值Table 3 The average of quantile of the reference value of different indicator in the corresponding assessment unit for different reference system

ESQref的四期平均值總體上呈現(xiàn)東南高(0.9)、西北低(0.07),濕潤(rùn)區(qū)高(0.62以上)、干旱區(qū)低(0.35以下)的格局(圖2)。全國(guó)參照生態(tài)系統(tǒng)的質(zhì)量平均值在生態(tài)工程大規(guī)模實(shí)施前后的兩個(gè)時(shí)段(即1990—2000年和2000—2015年)分別呈下降和上升的趨勢(shì)(圖1)。1990—2000年全國(guó)生態(tài)系統(tǒng)質(zhì)量的參照值下降了0.0013,ESQref下降的參照系主要集中在兩個(gè)區(qū)域,一是在暖溫帶亞濕潤(rùn)地區(qū)、中溫帶半干旱地區(qū)、北亞熱帶濕潤(rùn)地區(qū)以及中國(guó)中部多種生態(tài)地理分區(qū)的交界處(圖1),該區(qū)域年降水量減少123 mm,ESQref平均下降了0.06；二是中溫帶亞濕潤(rùn)地區(qū)和寒溫帶濕潤(rùn)地區(qū),該區(qū)域年降水量下降162.40 mm,年平均氣溫下降1.20℃,ESQref平均下降0.013。而ESQ增加的區(qū)域主要分布于三個(gè)區(qū)域,其中中亞熱帶濕潤(rùn)地區(qū)年降水量平均增加212 mm,總輻射平均增加35.03 MJ/m2,ESQref平均增加0.053；青藏高原溫帶干旱地區(qū)和青藏高原亞寒帶半干旱地區(qū)的東部地區(qū)年降水量增加46 mm,年平均氣溫上升0.57℃,ESQref平均增加了0.018；而暖溫帶干旱地區(qū)和中溫帶干旱地區(qū)的西部區(qū)域受年降水量增加(+22 mm)影響,ESQref平均增加0.021(圖1)。

2000—2015年全國(guó)參照生態(tài)系統(tǒng)的質(zhì)量值平均增加了0.0306,ESQref增加的參照系集中分布在中溫帶亞濕潤(rùn)地區(qū)、中溫帶濕潤(rùn)地區(qū),中溫帶半干旱地區(qū)、暖溫帶亞濕潤(rùn)地區(qū)以及北亞熱帶濕潤(rùn)地區(qū)(圖1)。受氣溫上升影響,中溫帶亞濕潤(rùn)地區(qū)和中溫帶濕潤(rùn)地區(qū)的ESQref增加了0.117；中溫帶半干旱地區(qū)和暖溫帶亞濕潤(rùn)地區(qū)ESQref平均增加0.136。而北亞熱帶濕潤(rùn)地區(qū)ESQref的增加(+0.123)主要與該區(qū)域年降水量增加(+380 mm)有關(guān)。該時(shí)段青藏高原亞寒帶半干旱地區(qū)東部地區(qū)受年降水量下降(-293 mm)影響,ESQref平均下降0.003(圖1)。

圖1 參照生態(tài)系統(tǒng)的質(zhì)量變化及其主要?dú)夂蛴绊懸蜃覨ig.1 Changes in the quality of reference ecosystems and their main climate influencing factors

2.2 基于參照系的中國(guó)陸地生態(tài)系統(tǒng)質(zhì)量的空間分布特征

生態(tài)系統(tǒng)質(zhì)量現(xiàn)實(shí)值(ESQ)和基于參照系的生態(tài)系統(tǒng)質(zhì)量(ESQ′)四期均值的空間分布如圖2所示。ESQ的全國(guó)平均值為0.305,與ESQref具有一致的空間分異規(guī)律。ESQ′平均值為0.72,表明在全國(guó)水平上,生態(tài)系統(tǒng)質(zhì)量現(xiàn)實(shí)值與參照值的平均差距為28%。ESQ′優(yōu)良的地區(qū)主要分布在青藏高原亞寒帶、溫度半干旱地區(qū)、亞寒帶亞濕潤(rùn)地區(qū)草地生態(tài)系統(tǒng),暖溫帶、中溫帶干旱地區(qū)森林生態(tài)系統(tǒng),南亞熱帶、中亞熱帶、北亞熱帶濕潤(rùn)地區(qū)森林生態(tài)系統(tǒng)的東部和北部,中溫帶、寒溫帶濕潤(rùn)地區(qū)和中溫帶亞濕潤(rùn)地區(qū)的森林生態(tài)系統(tǒng)。以西北地區(qū)為例,生態(tài)系統(tǒng)質(zhì)量現(xiàn)實(shí)值評(píng)估結(jié)果“較好”和“良好”水平的區(qū)域僅占17.72%,而基于參照系的生態(tài)系統(tǒng)質(zhì)量評(píng)估結(jié)果中,二者比例占41.16%(圖2)。例如,三江源草原草甸濕地生態(tài)功能區(qū)雖然生態(tài)系統(tǒng)質(zhì)量現(xiàn)實(shí)值較低,但該地區(qū)的生態(tài)系統(tǒng)功能指數(shù)優(yōu)于參照系(生態(tài)系統(tǒng)功能指數(shù)現(xiàn)實(shí)值/參照系=1.25)且生態(tài)系統(tǒng)穩(wěn)定性指數(shù)接近參照系(生態(tài)系統(tǒng)穩(wěn)定性指數(shù)現(xiàn)實(shí)值/參照系=0.92),因此ESQ接近參照系的質(zhì)量。這表明基于參照系的生態(tài)系統(tǒng)質(zhì)量在一定程度上克服了質(zhì)量好壞受氣候、自然地理?xiàng)l件主導(dǎo)的問(wèn)題。

圖2 生態(tài)系統(tǒng)質(zhì)量的空間分布和分區(qū)統(tǒng)計(jì)Fig.2 Spatial distribution and regional statistics of ecosystem quality

2.3 中國(guó)陸地生態(tài)系統(tǒng)質(zhì)量變化原因分析

圖3為1990—2000年和2000—2015年兩個(gè)時(shí)段ESQ′變化的空間分布。相對(duì)于參照系ESQ(ESQref)的變化而言,1990—2000年中國(guó)陸地生態(tài)系統(tǒng)質(zhì)量總體下降了1.2%,而2000—2015年生態(tài)系統(tǒng)質(zhì)量上升了3.69%。其中,ESQ′轉(zhuǎn)變比較明顯的區(qū)域主要分布在陜西-甘肅-寧夏(簡(jiǎn)稱(chēng)陜-甘-寧)、東北和西南石漠化治理工程區(qū)。

圖3 基于參照系的生態(tài)系統(tǒng)質(zhì)量(ESQ′)的變化Fig.3 Changes in ecosystem quality based on reference system

對(duì)于陜-甘-寧地區(qū),1990—2000年因總輻射量增加(+101.74 MJ/m2)會(huì)提升生態(tài)系統(tǒng)功能,參照值相應(yīng)增加0.01,而農(nóng)田擴(kuò)張導(dǎo)致該區(qū)域受人類(lèi)干擾加劇,生態(tài)系統(tǒng)結(jié)構(gòu)指數(shù)下降,同時(shí)引起葉面積指數(shù)下降(-0.13)和生態(tài)系統(tǒng)穩(wěn)定性的下降(-0.11),從而引起該區(qū)域質(zhì)量變化偏離參照系質(zhì)量變化。2000—2015年該區(qū)域ESQ′平均增加0.133,特別是黃土高原丘陵溝壑水土保持生態(tài)功能區(qū)(+0.203)。在ESQref平均增加0.041的背景下,退耕還林還草工程的大規(guī)模實(shí)施增加了該地區(qū)的林草地比例,這促使生態(tài)系統(tǒng)結(jié)構(gòu)得以改善(生態(tài)系統(tǒng)結(jié)構(gòu)指數(shù)增加0.089；LAI增加0.429),生態(tài)系統(tǒng)功能指數(shù)和生態(tài)系統(tǒng)穩(wěn)定性指數(shù)分別增加了0.139和0.186,以致生態(tài)系統(tǒng)質(zhì)量指數(shù)的增加量反而大于參照生態(tài)系統(tǒng)的質(zhì)量變化(表5)。

對(duì)于中國(guó)東北地區(qū),1990—2000年?yáng)|北平原的ESQ′下降受NPP(-60.89 gC m-2a-1)、NEP(-139.63 gC m-2a-1)和水源涵養(yǎng)功能下降(-49.59 mm)主導(dǎo),可能與降水的減少(-210.34 mm)、溫度的下降(-0.79℃)和土地利用變化有關(guān)；而長(zhǎng)白山地區(qū)受年均溫顯著下降(-1.25℃)影響,參照系質(zhì)量下降0.004,而該區(qū)域由6750 km2森林轉(zhuǎn)為農(nóng)田引起的穩(wěn)定性下降(-0.16)則主導(dǎo)了整個(gè)區(qū)域質(zhì)量下降(表4)。2000—2015年,雖然氣候變化的不利影響導(dǎo)致東北地區(qū)參照值下降了0.051,但長(zhǎng)白山森林生態(tài)功能區(qū)和大小興安嶺森林生態(tài)功能區(qū)表現(xiàn)出大面積的質(zhì)量恢復(fù),ESQ′增加主要由生態(tài)系統(tǒng)功能和穩(wěn)定性增加引起,這在很大程度上與三北防護(hù)林工程[53]有關(guān),該時(shí)段有6450 km2農(nóng)田轉(zhuǎn)換為森林。

表4 1990—2000年ESQ′轉(zhuǎn)變明顯區(qū)域的生態(tài)系統(tǒng)質(zhì)量、氣候因子和LAI的變化的平均值,以及主要的土地利用變化Table 4 The average value of changes in ecosystem quality,climate factors and LAI,as well as main land-use changes during 1990—2000 in the regions with obvious ESQ′ transformations

對(duì)于西南石漠化治理工程實(shí)施的縣域,1990—2000年ESQ增加的幅度低于參照值變化。而2000—2015年該地區(qū)降水和總輻射分別增加了92.77 mm和195.33 MJ/m2,參照生態(tài)系統(tǒng)的質(zhì)量提升了0.004；這期間由于石漠化綜合治理工程的實(shí)施,該地區(qū)分別有13375 km2和5150 km2的農(nóng)田轉(zhuǎn)為森林和草地,同時(shí)又有5800 km2的草地轉(zhuǎn)換成森林,生態(tài)系統(tǒng)結(jié)構(gòu)指數(shù)、生態(tài)系統(tǒng)功能指數(shù)和生態(tài)系統(tǒng)穩(wěn)定性指數(shù)分別增加了0.043、0.027和0.090,尤其是生產(chǎn)力和固碳功能分別增加了150.03 gC m-2a-1和 62.19 gC m-2a-1,最終促使整個(gè)區(qū)域ESQ的增幅(+0.051)遠(yuǎn)高于參照系。

3 討論

2000年前后中國(guó)實(shí)施了天然林資源保護(hù)、退耕還林還草、京津風(fēng)沙源治理等一系列重大生態(tài)建設(shè)工程,這些大規(guī)模的土地利用管理引領(lǐng)了全球陸地變綠,2000—2017年中國(guó)植被葉面積的增加達(dá)到全球綠化面積增加的25%[13]。在氣候變化和人類(lèi)活動(dòng)的共同作用下,2000—2010年中國(guó)多項(xiàng)生態(tài)系統(tǒng)服務(wù)(包括食物生產(chǎn)、水源涵養(yǎng)、土壤保持、防風(fēng)固沙、洪水調(diào)蓄、固碳)得到明顯改善[17]。本研究基于參照系的生態(tài)系統(tǒng)質(zhì)量(ESQ′)評(píng)估結(jié)果進(jìn)一步表明,因人類(lèi)活動(dòng)引起的生態(tài)系統(tǒng)質(zhì)量變化在2000—2015年間增加了3.69%,尤其在黃土高原丘陵溝壑水土保持生態(tài)功能區(qū)、青海省東部、長(zhǎng)白山地區(qū)質(zhì)量增加明顯,這與Ouyang等研究中2000—2010年多項(xiàng)生態(tài)系統(tǒng)服務(wù)增加的區(qū)域一致[17]。同時(shí),本研究還表明西南石漠化工程治理縣域和長(zhǎng)白山地區(qū)這兩個(gè)LAI“變綠”明顯的區(qū)域[13],ESQ′也存在明顯增加。統(tǒng)計(jì)結(jié)果表明,這兩個(gè)地區(qū)實(shí)施了大量的造林工程(退耕還林工程[5,54],三北防護(hù)林工程[53]),2000—2015年分別有13375 km2和6450 km2的農(nóng)田轉(zhuǎn)為森林(表5)。此外,在中國(guó)西南喀斯特區(qū)域,隨著退耕還林還草、石漠化綜合治理等工程的實(shí)施,2000—2015年貴州和廣西省的ESQ′增加明顯,這與Tong等NDVI變化(2000—2011年與1982—2000年NDVI變化趨勢(shì)之差)的空間格局一致[55]。

表5 2000—2015年ESQ′轉(zhuǎn)變明顯區(qū)域的生態(tài)系統(tǒng)質(zhì)量、氣候因子和LAI的變化的平均值,以及主要的土地利用變化Table 5 The average value of changes in ecosystem quality,climate factors and LAI,as well as main land-use changes during 2000—2015 in the regions with obvious ESQ′ transformations

本研究是從人類(lèi)脅迫、生態(tài)系統(tǒng)結(jié)構(gòu)、功能和穩(wěn)定性四個(gè)方面評(píng)估生態(tài)系統(tǒng)質(zhì)量。造林工程能直接增加植被覆蓋率,改善生態(tài)系統(tǒng)結(jié)構(gòu)[56],而退耕還林還草工程能促使生態(tài)系統(tǒng)由農(nóng)田的高人類(lèi)脅迫轉(zhuǎn)為草地或森林的低人類(lèi)脅迫[52]。由于森林的抗旱能力高于農(nóng)田[57],退耕還林工程同時(shí)也會(huì)增加生態(tài)系統(tǒng)的穩(wěn)定性。生態(tài)工程對(duì)生態(tài)系統(tǒng)功能的綜合效益因氣候背景而異,在降水充足的地區(qū),造林無(wú)疑會(huì)增加生態(tài)系統(tǒng)多項(xiàng)功能,如在西南喀斯特地區(qū)生態(tài)系統(tǒng)功能之間的關(guān)系以協(xié)同為主[58]；然而,在降水限制區(qū)生態(tài)工程對(duì)生態(tài)系統(tǒng)功能的影響還存在較大爭(zhēng)議,如三北防護(hù)林工程京津風(fēng)沙源治理工程增加了半干旱地區(qū)防風(fēng)固沙和土壤保持功能,但減少了地表徑流和加速了地下水的消耗[59]。綜合這些研究可以看出,生態(tài)工程在短期內(nèi)會(huì)通過(guò)減小人類(lèi)脅迫、改善生態(tài)系統(tǒng)結(jié)構(gòu)、功能,提升穩(wěn)定性來(lái)增加生態(tài)系統(tǒng)質(zhì)量。此外,人工林會(huì)隨時(shí)間變化生長(zhǎng)而導(dǎo)致需水量增加,由此而產(chǎn)生的與生產(chǎn)力、固碳和土壤保持的權(quán)衡關(guān)系的增強(qiáng)(尤其是在半干旱區(qū)和干旱區(qū))在長(zhǎng)期內(nèi)如何影響生態(tài)系統(tǒng)質(zhì)量還需要進(jìn)一步的考量。

根據(jù)《西部地區(qū)重點(diǎn)生態(tài)區(qū)綜合治理規(guī)劃綱要(2012-2020年)》,在西北草原荒漠化防治區(qū)、黃土高原水土保持區(qū)、青藏高原江河水源涵養(yǎng)區(qū)、西南石漠化防治區(qū)、重要森林生態(tài)功能區(qū)等五大重點(diǎn)生態(tài)區(qū),退耕還林、退牧還草、自然保護(hù)區(qū)建設(shè)和水土流失治理等生態(tài)治理區(qū)域面積到2020年將分別達(dá)到151.1×104km2、27.6×104km2、88×104km2、26.7×104km2、106.2×104km2的面積。同時(shí),根據(jù)石漠化綜合治理工程“十三五”規(guī)劃,到2020年巖溶土地和石漠化土地治理面積分別達(dá)到5×104km2和2×104km2以上,建設(shè)和保護(hù)林草植被195×104hm2[60]。我國(guó)的森林覆蓋率從1998年的13.92%提升至2015年的21.63%[61],預(yù)期到2025年將增加至24.1%[62]。我國(guó)《十四五規(guī)劃和2035年遠(yuǎn)景目標(biāo)綱要》指出,要堅(jiān)持山水林田湖草治理,通過(guò)在重點(diǎn)生態(tài)功能區(qū)、國(guó)家級(jí)自然保護(hù)區(qū)等地實(shí)施生態(tài)系統(tǒng)保護(hù)和修復(fù)工程,改善生態(tài)系統(tǒng)質(zhì)量。這些工程將布局在“三生態(tài)區(qū)-四森林帶”,預(yù)期在青藏高原生態(tài)屏障區(qū)(如三江源、祁連山)治理退化草原320×104hm2,在黃河重點(diǎn)生態(tài)區(qū)(如黃土高原、秦嶺)治理水土流失和沙地280×104hm2,保護(hù)林草地80×104hm2,在長(zhǎng)江重點(diǎn)生態(tài)區(qū)造林110×104hm2,治理水土流失和石漠化面積共600×104hm2。同時(shí)將在東北森林帶、北方防沙帶和南方丘陵帶造林約300×104hm2,并且預(yù)期將在北方防沙帶治理沙化土地750×104hm2。這些積極的生態(tài)保護(hù)工程預(yù)估將進(jìn)一步增強(qiáng)生態(tài)脆弱區(qū)和重點(diǎn)生態(tài)功能區(qū)的生態(tài)系統(tǒng)結(jié)構(gòu)與穩(wěn)定性,提升生態(tài)系統(tǒng)功能和質(zhì)量。受全球變暖影響,預(yù)計(jì)未來(lái)三十年(2021—2050)中國(guó)大部分地區(qū)將遭受更嚴(yán)重、更頻繁和持續(xù)時(shí)間更長(zhǎng)的干旱,尤其在西北地區(qū)[63],而2040年后中國(guó)降水整體呈增加趨勢(shì)(集中在北方地區(qū))[64]。因此未來(lái)幾十年降水能否承載起生態(tài)工程導(dǎo)致的森林增加和生長(zhǎng)的需求是個(gè)很大的挑戰(zhàn),未來(lái)我國(guó)半干旱半濕潤(rùn)區(qū)的生態(tài)工程的布設(shè)應(yīng)關(guān)注具有變濕趨勢(shì)的地區(qū)。

本研究通過(guò)考慮參照系的生態(tài)系統(tǒng)質(zhì)量變化,在一定程度上消除不同區(qū)域氣候波動(dòng)差異對(duì)質(zhì)量評(píng)估的影響,為量化保護(hù)、治理等生態(tài)工程對(duì)生態(tài)系統(tǒng)的影響提供了一種新途徑。例如長(zhǎng)白山參照系無(wú)土地利用變化,其ESQ可能會(huì)隨氣候限制因子的消極變化而下降(-0.051),而整個(gè)區(qū)域的ESQ能抵消氣候變化的不利影響進(jìn)而使生態(tài)系統(tǒng)質(zhì)量有較大的提升(+0.110)。ESQ′在量化生態(tài)工程等人為影響方面具有較大潛力,可以作為政府部門(mén)績(jī)效評(píng)估的依據(jù)之一,也能為生態(tài)系統(tǒng)和土地利用管理提供參考。

參照系的合理選取是評(píng)價(jià)生態(tài)系統(tǒng)質(zhì)量變化的重要前提。本研究中大部分評(píng)估單元的參照系符合優(yōu)良生態(tài)系統(tǒng)的假設(shè),然而受到數(shù)據(jù)獲取的限制,參照系選取主要基于CERN定位觀測(cè)臺(tái)站永久樣地和國(guó)家自然保護(hù)區(qū)核心區(qū),尚不能充分利用國(guó)家林草局、國(guó)家氣象局、生態(tài)環(huán)境部、高校等其他部門(mén)野外觀測(cè)臺(tái)站的永久樣地?cái)?shù)據(jù)。因此,僅考慮極小干擾參照系很難覆蓋全部的評(píng)估單元,本文根據(jù)數(shù)據(jù)的可獲取性選擇最受氣候影響區(qū)和評(píng)估單元優(yōu)良值基準(zhǔn)作為剩余評(píng)估單元的參照系,基于多種基準(zhǔn)參照系的基礎(chǔ)上對(duì)國(guó)家尺度的生態(tài)系統(tǒng)質(zhì)量進(jìn)行評(píng)估；同時(shí)由于數(shù)據(jù)觀測(cè)和處理、模型模擬過(guò)程中存在的隨機(jī)誤差和系統(tǒng)誤差[65-67],導(dǎo)致部分評(píng)估單元的參照系并不完全符合優(yōu)良生態(tài)系統(tǒng)的假設(shè),一些功能處于評(píng)估單元內(nèi)的中等水平。未來(lái)國(guó)家不同部門(mén)野外觀測(cè)研究站網(wǎng)規(guī)模的擴(kuò)大和數(shù)據(jù)共享,以及中國(guó)生態(tài)保護(hù)紅線(xiàn)地面觀測(cè)網(wǎng)絡(luò)的建立,形成一體化、長(zhǎng)期、連續(xù)的生態(tài)系統(tǒng)質(zhì)量監(jiān)測(cè)體系,將為我國(guó)生態(tài)系統(tǒng)質(zhì)量評(píng)估提供更多可靠的參照系數(shù)據(jù)、模型驗(yàn)證和校準(zhǔn)數(shù)據(jù),使基于參照系的生態(tài)系統(tǒng)質(zhì)量評(píng)估結(jié)果更具可比性和科學(xué)性,從而更好地服務(wù)于我國(guó)生態(tài)文明建設(shè)、重點(diǎn)生態(tài)區(qū)綜合治理、石漠化綜合治理、“山水林田湖草”等工程的治理成效評(píng)估和績(jī)效考核。

4 結(jié)論

本文基于1990—2015年中國(guó)土地覆蓋、植被覆蓋度、生態(tài)系統(tǒng)功能和氣候等數(shù)據(jù),計(jì)算了生態(tài)系統(tǒng)結(jié)構(gòu)、功能、穩(wěn)定性和人類(lèi)脅迫指數(shù)；根據(jù)生態(tài)地理分區(qū)和植被類(lèi)型劃分的評(píng)估單元,基于CERN臺(tái)站、國(guó)家自然保護(hù)區(qū)和最受氣候影響區(qū)等參照生態(tài)系統(tǒng),結(jié)合生態(tài)系統(tǒng)質(zhì)量評(píng)估模型,評(píng)估了基于參照系的中國(guó)陸地生態(tài)系統(tǒng)質(zhì)量及其變化情況。主要研究結(jié)論如下：

(1)本研究選用的生態(tài)系統(tǒng)質(zhì)量評(píng)估參照系具有較低的人類(lèi)脅迫和良好的生態(tài)系統(tǒng)功能和穩(wěn)定性,參照系的生態(tài)系統(tǒng)質(zhì)量指數(shù)總體位于評(píng)估單元89%分位數(shù)。1990—2000年和2000—2015年參照系ESQI的全國(guó)平均值分別呈下降和上升趨勢(shì)。與參照系相比,中國(guó)陸地生態(tài)系統(tǒng)質(zhì)量指數(shù)的平均值偏低28%。

(2)1990—2000年和2000—2015年ESQ′分別下降1.2%和上升3.69%。ESQ′轉(zhuǎn)變比較明顯的區(qū)域主要分布在“陜-甘-寧”地區(qū)、長(zhǎng)白山地區(qū)和西南石漠化治理工程實(shí)施縣域,在這些區(qū)域?qū)嵤┑耐烁€林還草工程、三北防護(hù)林工程和石漠化綜合治理工程發(fā)揮了重要作用。未來(lái)我國(guó)計(jì)劃進(jìn)一步實(shí)施生態(tài)保護(hù)和修復(fù)工程,加強(qiáng)山水林田湖草沙治理,這些措施將進(jìn)一步提升我國(guó)生態(tài)系統(tǒng)質(zhì)量。同時(shí)需要關(guān)注未來(lái)氣候變化以及生態(tài)系統(tǒng)功能之間的權(quán)衡關(guān)系對(duì)生態(tài)系統(tǒng)質(zhì)量的影響。

(3)針對(duì)傳統(tǒng)評(píng)估方法因包含了地理、氣候本底條件差異而存在評(píng)估結(jié)果在空間上不可比的問(wèn)題,基于參照系的生態(tài)系統(tǒng)質(zhì)量評(píng)估方法通過(guò)考慮參照系的生態(tài)系統(tǒng)質(zhì)量變化,可以在一定程度上消除不同區(qū)域氣候波動(dòng)差異對(duì)質(zhì)量評(píng)估的影響,為科學(xué)量化保護(hù)、治理等生態(tài)工程對(duì)生態(tài)系統(tǒng)的影響提供了一種新途徑。未來(lái)應(yīng)加強(qiáng)全國(guó)不同區(qū)域、不同類(lèi)型的參照生態(tài)系統(tǒng)聯(lián)網(wǎng)觀測(cè),為提升我國(guó)生態(tài)系統(tǒng)質(zhì)量評(píng)估能力提供科學(xué)數(shù)據(jù)支撐。