张梦宇,张 黎,何洪林,任小丽,牛忠恩,吕 妍,徐 茜,常清青,刘卫华,李 攀

1 中国科学院地理科学与资源研究所 生态系统网络观测与模拟重点实验室,北京 100101 2 国家生态科学数据中心,北京 100101 3 中国科学院大学,北京 100049 4 中国科学院大学资源与环境学院,北京 100190 5 天津大学 地球系统科学学院,天津 300072

生态系统质量是诊断生态系统是否健康的指标,是指在一个具体的时间和空间范围内生态系统的总体或部分生命组分的质量,主要表现在其生产能力、自我维持能力、受到外界干扰后的抗干扰能力以及对人类的生存及社会经济持续发展的影响[1-4]。一个高质量的生态系统能够在压力下维持自身组织结构和稳定,具有自我恢复的能力[1]。改革开放后中国人口与经济快速增长,耕地不断扩张[5],森林和草地等生态系统发生了退化[6-7],带来自然栖息地破坏和生态系统功能下降等生态问题[8-10]。但自1999年以来中国实施了草地保护、天然林保护、退耕还林以及石漠化综合治理等工程,退化的生态系统得到一定的恢复[11-14]。党中央在十九大报告中将生态文明提升到战略地位,要通过加大生态系统保护力度,实施重要生态系统保护和修复重大工程,以提升生态系统质量和稳定性。在当前我国土地利用需求不断增加和生态工程大规模实施的背景下,迫切需要对我国生态系统质量状况进行评估,为土地利用和生态系统可持续管理提供参考[15]。

我国已开展一系列与生态系统质量有关的评估工作,主要包括生态系统功能评估[16-17]、生态系统健康评估[18]和生态环境质量评估[19]等方面,为国家和地方生态保护政策的制定提供了科学依据。在全国层面,根据“全国生态环境十年变化(2000—2010年)调查评估”项目研究结果,2000—2010年间食物生产、水源涵养、土壤保持、防风固沙、洪水调蓄、固碳这6项生态系统服务得到明显改善,而生物多样性维护功能有所下降[17]；中国生态系统健康水平呈现从西北到东南上升的格局,水分指数和土地利用强度的贡献分别为24.5%和20.7%[18]。我国生态环境部依据生态环境状况评价技术规范(HJ 192—2015),定期发布《全国生态环境质量报告》,公布全国省域和县域生态环境状况及变化评估结果。此外,国内学者也开展了一些典型区域的生态系统质量评价方法和应用研究。对于自然生态系统,主要从生产能力、功能稳定性和人口压力[20],植被覆盖度和叶面积指数[21],生态系统生产能力、稳定性和承载力[22-23]等方面评估生态系统质量。对于城市生态系统,重点从经济、自然和社会环境质量[24],结构、功能和胁迫[25]等角度进行评估。对于干旱区流域水文生态系统,综合植被覆盖度、垦区土壤盐渍化率、沙漠化面积比率、水土流失强度和生态用水比率等因素建立生态系统质量框架[26]。计算生态系统质量或健康综合指数的方法主要包括求积法[27-28]、PSR模型[2,29]、主成分分析法[30]、层次分析法[31]等。然而,以往有关生态系统功能、健康和质量的评估往往聚焦于生态系统指标的现实值,其数值大小在很大程度上受气候因子(如降水、温度、总辐射)和地理本底条件的差异主导[32-33],导致不同地区评价结果的可比性较低,并且对于全国范围生态系统质量及其变化的认识仍较为缺乏,尚不足以满足国家生态系统管理的需求[15]。

根据生态系统完整性[4,34-36]和稳定性[37]理论,生态系统质量评估不仅应综合考虑一定时间、空间范围内生态系统结构、功能、稳定性特征,以及人类胁迫或干扰状况,从而反映生态系统的整体质量状态；还应分区、分类型地开展评估参照生态系统(简称参照系)的设定,使得评估结果在不同时间和空间上具有较好的可比性[15]。由于现实世界中很难找到完全不受人类影响的理想参照系,将与待评估生态系统具有相似植被、土壤、气候条件且受人类活动影响相对较小的自然生态系统作为参照系具有较好的合理性与可行性。参照生态系统主要包括：(1)极小干扰参照状态,生态系统最大程度地保留了完整性,几乎没有人类活动；(2)历史条件参照状态,在人类干扰之前的生态系统状态；(3)最佳可达状态,在某一特定地理和气候本底条件下,生态系统所能达到的最佳状态,可以根据分位数法选取同一条件区域内一组参照点作为参照系[38-40]。针对我国陆地生态系统特点和数据可获取性状况,将我国国家自然保护区核心区、中国生态系统研究网络(CERN)野外观测台站永久观测样地等较少受人类活动干扰且状态优良的生态系统作为参照系,是当前我国陆地生态系统质量评估的一种新思路[15]。

本文利用基于参照系的生态系统质量评估方法,以中国陆地生态系统为研究对象,综合土地覆盖数据,降水、气温、总辐射的年度气候数据和生态系统生产力、固碳、水源涵养和土壤保持四项生态系统功能数据,评估1990、2000、2010、2015年四期生态系统质量,重点分析参照生态系统的质量变化特征,基于参照系的中国陆地生态系统质量的空间格局,以及生态工程大规模实施前后的两个时段(即1990—2000年和2000—2015年)中国陆地生态系统质量变化特征及其归因,以期为我国生态工程的成效评估提供依据,为土地利用和生态系统管理提供支撑。

1 数据和方法

1.1 数据

1.1.1土地覆盖数据

采用1990年、2000年、2010年、2015年空间分辨率为1 km的土地覆盖数据(ChinaCover)[41]计算生态系统结构指数。该数据是基于国产环境灾害卫星(HJ-1A/B)和美国陆地卫星(Landsat OLI)数据,采用面向对象的多尺度分割与变化检测的分类方法进行制图,一级地类的精度达94%。

1.1.2环境要素数据

采用1990—2015年1 km分辨率的气候要素(年均温、年降水、年总辐射)数据分析气候因子对生态系统质量变化的影响。其中,气温和降水数据基于中国地面气候日值数据集和全球历史气候网格数据集采用ANUSPLIN软件插值计算生成[42],总辐射数据由日照时数计算生成[43]。根据《1∶1000000 中国植被图集》和中国生态地理分区生成168个评估单元,数据源自中国科学院资源环境科学数据中心(http://www.resdc.cn/)。

1.1.3生态系统功能数据

1990—2015年生态系统功能指标数据由生态系统过程模型(Carbon Exchange between Vegetation,Soil,and Atmosphere—Ecosystem services,CEVSA-ES[44])模拟得到。CEVSA-ES模型模拟的NPP、NEP、土壤含水量和土壤保持量与观测数据具有较好的一致性,二者线性拟合的决定系数分别为0.96、0.81、0.47、0.52[44]。模拟结果的空间分辨率为0.05°。

1.1.4其他数据

1990—2015年GLASS LAI数据[45-46]用于分析植被生长状况,数据源自在国家地球系统科学数据中心,其空间分辨率为1 km,时间分辨率为8 d。本文将其计算为年最大LAI(LAImax)。国家自然保护区空间分布数据(截至2013年9月)用于选取参照系,数据来自于中国科学院资源环境科学数据中心(http://www.resdc.cn/)。CERN野外观测台站地理位置数据从国家生态科学数据中心(http://rs.cern.ac.cn/data/initDRsearch?cid=STA_DP2011)获取。采用1990—2015年1 km分辨率的植被覆盖度数据计算生态系统结构指数,该数据是以基于MODIS NDVI数据,利用像元二分模型计算得到,在计算过程中利用ChinaCover产品对纯土壤、植被像元的NDVI进行了改进[41]。

1.2 基于参照系的国家生态系统质量评价方法

根据生态系统完整性和稳定性理论[4,34-37],本文从生态系统结构、功能、稳定性和人类胁迫四个方面构建生态系统质量评价指标(表1),将生态系统质量指数(Ecosystem quality Index,ESQI)表示为生态系结构指数、功能指数、稳定性指数统和胁迫指数的加权和(公式1)。ESQI高表示人类胁迫小,生态系统景观结构优良,并具有良好且稳定的功能。

表1 国家生态系统质量指标体系Table 1 Indicator framework for national ecosystem quality

(1)

式中,xi为一级指标,包括标准化结构指数、标准化功能指数、标准化稳定性指数和标准化胁迫指数。ωi为对应的权重,根据熵权法计算[47-48],分别为0.299、0.335、0.317和0.048。

生态系统结构指标用生态系统结构指数(Ecosystem Structure Index,EStrI)表示,由不同土地覆盖类型覆盖率的加权和表示,参照生态环境部2015年发布的生态环境状况评价技术规范,森林、草地、水域湿地、耕地、建设用地、未利用地的权重分别设为0.35、0.21、0.28、0.11、0.04和0.01。

(2)

式中,flc为土地覆盖类型lc的覆盖率,ωlc为单位面积的土地覆盖类型lc的权重。

生态系统功能指标由生态系统功能指数(Ecosystem Function Index,EFunI)表示,设定为四项重要生态系统功能(包括物质生产、固碳、水源涵养和土壤保持功能)指标的加权和。分别采用净初级生产力(Net primary productivity,NPP)、净生态系统生产力(Net ecosystem productivity,NEP)、土壤含水量及土壤保持量表征物质生产、固碳、水源涵养和土壤保持功能。NPP、NEP和土壤含水量采用min—max标准化到0—1之间。由于土壤保持量呈现明显的右偏态分布,本文对其数值取对数后再采用min—max标准化到0—1之间。最后利用熵权法确定四个功能指标的权重。

生态系统稳定性指数(Ecosystem Stability Index,EStaI)用待评估时间点的前十年NPP的时间稳定性指数表示。该指数假设生态系统越接近长时间序列的平均状态就越稳定,能较好地综合反映生态系统的抵抗力和恢复力[49]。考虑到中国大规模实施了生态工程并具有城市扩张的背景,生态系统具有一定的恢复或退化的趋势,本研究采用去除时间趋势后的生产力时间稳定性(公式3)来表征中国生态系统的稳定性指数,并采用min—max方法将其标准化到0—1之间。

(3)

式中,meanfun为生态系统功能的多年平均值,stddetrend_fun为去趋势后生态系统功能的标准差。

人类胁迫指标用人类干扰指数(Human Interference Index,HII)表示,考虑每种土地覆盖类型的面积(覆盖率)及其对生态系统的干扰强度[50-51],同时考虑周围土地类型对其干扰的距离效应[52]。采用min—max标准化方法将其负向标准化到0—1之间。

(4)

表2 土地覆盖类型的人类干扰度范围Table 2 The range of human interference index of different land cover types

为尽可能消除评价结果空间不可比问题,依据每个评估单元内较少受人类活动明显干扰且优良的参照系,计算基于参照系的生态系统质量指数(ESQI′；公式5)。

ESQI′=ESQI/ESQIref

(5)

ESQIref为生态系统质量的参照值,ESQI为生态系统质量的现实值。当ESQI′≥1时,表示生态系统质量优于参照系；反之生态系统质量劣于参照系。根据分位数将ESQI′划分为优良(80%—100%)、较好(60%—80%)、中等(40%—60%)、较差(20%—40%)、很差(0%—20%)五个等级。生态系统质量的参照值是指某一时期某个评估单元内的参照系的质量的平均值。本研究中针对每个具有相似地理和气候条件的评估单元,按照下面方法的优先次序筛选出评估单元内适合作为参照生态系统的备选区域：(1)具有CERN野外观测台站或国家自然保护区核心区等较少受人类影响的生态系统；(2)净初级生产力主要受气候因子控制的生态系统；(3)评估单元内生态系统质量指数为优良的生态系统。在这些备选区域内选择未发生土地利用变化且生态系统质量相对优良的区域作为参照系,下文分别简称为CERN参照、国家自然保护区参照、最受气候影响区参照和评估单元参照。

2 结果

2.1 参照系的生态系统质量变化特征

对1990年、2000年、2010年、2015年不同类型参照系的生态系统质量指标参照值进行统计,其占评估单元总体的分位数平均值如表3所示。总体上,参照系的生态系统质量指数数值位于评估单元89%分位数,其中胁迫标准化指数、结构标准化指数、功能标准化指数、稳定性标准化指数四个分量的参照值分别位于76%、80%、75%和69%。按照不同类型的参照系来看,以国家自然保护区为参照的评估单元数最多,参照系的人类干扰标准化指数位于评估单元的90.8%分位数,在四类参照系中的受干扰程度最低。以CERN永久样地为基准的参照生态系统具有相对较好的结构、功能和稳定性,均位于评估单元分位数的60%以上。对于最受气候影响的参照系,其结构指数所处于评估单元的分位数水平(77%)>功能指数(69.7%)>稳定性指数(68.3%)。以评估单元优良值为基准的参照生态系统的质量位于评估单元的91.6%分位数。可以看出,在具有相似气候和地理条件的评估单元内,参照系具有较低的人类胁迫和良好的生态系统结构、功能和稳定性,说明本研究根据CERN生态系统永久样地、国家自然保护区核心区、最受气候影响区来选取受人类干扰较少且优良的生态系统作为参照系是合理可行的。

表3 不同类型参照系的质量参照值位于相应评估单元分位数的平均值Table 3 The average of quantile of the reference value of different indicator in the corresponding assessment unit for different reference system

ESQref的四期平均值总体上呈现东南高(0.9)、西北低(0.07),湿润区高(0.62以上)、干旱区低(0.35以下)的格局(图2)。全国参照生态系统的质量平均值在生态工程大规模实施前后的两个时段(即1990—2000年和2000—2015年)分别呈下降和上升的趋势(图1)。1990—2000年全国生态系统质量的参照值下降了0.0013,ESQref下降的参照系主要集中在两个区域,一是在暖温带亚湿润地区、中温带半干旱地区、北亚热带湿润地区以及中国中部多种生态地理分区的交界处(图1),该区域年降水量减少123 mm,ESQref平均下降了0.06；二是中温带亚湿润地区和寒温带湿润地区,该区域年降水量下降162.40 mm,年平均气温下降1.20℃,ESQref平均下降0.013。而ESQ增加的区域主要分布于三个区域,其中中亚热带湿润地区年降水量平均增加212 mm,总辐射平均增加35.03 MJ/m2,ESQref平均增加0.053；青藏高原温带干旱地区和青藏高原亚寒带半干旱地区的东部地区年降水量增加46 mm,年平均气温上升0.57℃,ESQref平均增加了0.018；而暖温带干旱地区和中温带干旱地区的西部区域受年降水量增加(+22 mm)影响,ESQref平均增加0.021(图1)。

2000—2015年全国参照生态系统的质量值平均增加了0.0306,ESQref增加的参照系集中分布在中温带亚湿润地区、中温带湿润地区,中温带半干旱地区、暖温带亚湿润地区以及北亚热带湿润地区(图1)。受气温上升影响,中温带亚湿润地区和中温带湿润地区的ESQref增加了0.117；中温带半干旱地区和暖温带亚湿润地区ESQref平均增加0.136。而北亚热带湿润地区ESQref的增加(+0.123)主要与该区域年降水量增加(+380 mm)有关。该时段青藏高原亚寒带半干旱地区东部地区受年降水量下降(-293 mm)影响,ESQref平均下降0.003(图1)。

图1 参照生态系统的质量变化及其主要气候影响因子Fig.1 Changes in the quality of reference ecosystems and their main climate influencing factors

2.2 基于参照系的中国陆地生态系统质量的空间分布特征

生态系统质量现实值(ESQ)和基于参照系的生态系统质量(ESQ′)四期均值的空间分布如图2所示。ESQ的全国平均值为0.305,与ESQref具有一致的空间分异规律。ESQ′平均值为0.72,表明在全国水平上,生态系统质量现实值与参照值的平均差距为28%。ESQ′优良的地区主要分布在青藏高原亚寒带、温度半干旱地区、亚寒带亚湿润地区草地生态系统,暖温带、中温带干旱地区森林生态系统,南亚热带、中亚热带、北亚热带湿润地区森林生态系统的东部和北部,中温带、寒温带湿润地区和中温带亚湿润地区的森林生态系统。以西北地区为例,生态系统质量现实值评估结果“较好”和“良好”水平的区域仅占17.72%,而基于参照系的生态系统质量评估结果中,二者比例占41.16%(图2)。例如,三江源草原草甸湿地生态功能区虽然生态系统质量现实值较低,但该地区的生态系统功能指数优于参照系(生态系统功能指数现实值/参照系=1.25)且生态系统稳定性指数接近参照系(生态系统稳定性指数现实值/参照系=0.92),因此ESQ接近参照系的质量。这表明基于参照系的生态系统质量在一定程度上克服了质量好坏受气候、自然地理条件主导的问题。

图2 生态系统质量的空间分布和分区统计Fig.2 Spatial distribution and regional statistics of ecosystem quality

2.3 中国陆地生态系统质量变化原因分析

图3为1990—2000年和2000—2015年两个时段ESQ′变化的空间分布。相对于参照系ESQ(ESQref)的变化而言,1990—2000年中国陆地生态系统质量总体下降了1.2%,而2000—2015年生态系统质量上升了3.69%。其中,ESQ′转变比较明显的区域主要分布在陕西-甘肃-宁夏(简称陕-甘-宁)、东北和西南石漠化治理工程区。

图3 基于参照系的生态系统质量(ESQ′)的变化Fig.3 Changes in ecosystem quality based on reference system

对于陕-甘-宁地区,1990—2000年因总辐射量增加(+101.74 MJ/m2)会提升生态系统功能,参照值相应增加0.01,而农田扩张导致该区域受人类干扰加剧,生态系统结构指数下降,同时引起叶面积指数下降(-0.13)和生态系统稳定性的下降(-0.11),从而引起该区域质量变化偏离参照系质量变化。2000—2015年该区域ESQ′平均增加0.133,特别是黄土高原丘陵沟壑水土保持生态功能区(+0.203)。在ESQref平均增加0.041的背景下,退耕还林还草工程的大规模实施增加了该地区的林草地比例,这促使生态系统结构得以改善(生态系统结构指数增加0.089；LAI增加0.429),生态系统功能指数和生态系统稳定性指数分别增加了0.139和0.186,以致生态系统质量指数的增加量反而大于参照生态系统的质量变化(表5)。

对于中国东北地区,1990—2000年东北平原的ESQ′下降受NPP(-60.89 gC m-2a-1)、NEP(-139.63 gC m-2a-1)和水源涵养功能下降(-49.59 mm)主导,可能与降水的减少(-210.34 mm)、温度的下降(-0.79℃)和土地利用变化有关；而长白山地区受年均温显著下降(-1.25℃)影响,参照系质量下降0.004,而该区域由6750 km2森林转为农田引起的稳定性下降(-0.16)则主导了整个区域质量下降(表4)。2000—2015年,虽然气候变化的不利影响导致东北地区参照值下降了0.051,但长白山森林生态功能区和大小兴安岭森林生态功能区表现出大面积的质量恢复,ESQ′增加主要由生态系统功能和稳定性增加引起,这在很大程度上与三北防护林工程[53]有关,该时段有6450 km2农田转换为森林。

表4 1990—2000年ESQ′转变明显区域的生态系统质量、气候因子和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

对于西南石漠化治理工程实施的县域,1990—2000年ESQ增加的幅度低于参照值变化。而2000—2015年该地区降水和总辐射分别增加了92.77 mm和195.33 MJ/m2,参照生态系统的质量提升了0.004；这期间由于石漠化综合治理工程的实施,该地区分别有13375 km2和5150 km2的农田转为森林和草地,同时又有5800 km2的草地转换成森林,生态系统结构指数、生态系统功能指数和生态系统稳定性指数分别增加了0.043、0.027和0.090,尤其是生产力和固碳功能分别增加了150.03 gC m-2a-1和 62.19 gC m-2a-1,最终促使整个区域ESQ的增幅(+0.051)远高于参照系。

3 讨论

2000年前后中国实施了天然林资源保护、退耕还林还草、京津风沙源治理等一系列重大生态建设工程,这些大规模的土地利用管理引领了全球陆地变绿,2000—2017年中国植被叶面积的增加达到全球绿化面积增加的25%[13]。在气候变化和人类活动的共同作用下,2000—2010年中国多项生态系统服务(包括食物生产、水源涵养、土壤保持、防风固沙、洪水调蓄、固碳)得到明显改善[17]。本研究基于参照系的生态系统质量(ESQ′)评估结果进一步表明,因人类活动引起的生态系统质量变化在2000—2015年间增加了3.69%,尤其在黄土高原丘陵沟壑水土保持生态功能区、青海省东部、长白山地区质量增加明显,这与Ouyang等研究中2000—2010年多项生态系统服务增加的区域一致[17]。同时,本研究还表明西南石漠化工程治理县域和长白山地区这两个LAI“变绿”明显的区域[13],ESQ′也存在明显增加。统计结果表明,这两个地区实施了大量的造林工程(退耕还林工程[5,54],三北防护林工程[53]),2000—2015年分别有13375 km2和6450 km2的农田转为森林(表5)。此外,在中国西南喀斯特区域,随着退耕还林还草、石漠化综合治理等工程的实施,2000—2015年贵州和广西省的ESQ′增加明显,这与Tong等NDVI变化(2000—2011年与1982—2000年NDVI变化趋势之差)的空间格局一致[55]。

表5 2000—2015年ESQ′转变明显区域的生态系统质量、气候因子和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

本研究是从人类胁迫、生态系统结构、功能和稳定性四个方面评估生态系统质量。造林工程能直接增加植被覆盖率,改善生态系统结构[56],而退耕还林还草工程能促使生态系统由农田的高人类胁迫转为草地或森林的低人类胁迫[52]。由于森林的抗旱能力高于农田[57],退耕还林工程同时也会增加生态系统的稳定性。生态工程对生态系统功能的综合效益因气候背景而异,在降水充足的地区,造林无疑会增加生态系统多项功能,如在西南喀斯特地区生态系统功能之间的关系以协同为主[58]；然而,在降水限制区生态工程对生态系统功能的影响还存在较大争议,如三北防护林工程京津风沙源治理工程增加了半干旱地区防风固沙和土壤保持功能,但减少了地表径流和加速了地下水的消耗[59]。综合这些研究可以看出,生态工程在短期内会通过减小人类胁迫、改善生态系统结构、功能,提升稳定性来增加生态系统质量。此外,人工林会随时间变化生长而导致需水量增加,由此而产生的与生产力、固碳和土壤保持的权衡关系的增强(尤其是在半干旱区和干旱区)在长期内如何影响生态系统质量还需要进一步的考量。

根据《西部地区重点生态区综合治理规划纲要(2012-2020年)》,在西北草原荒漠化防治区、黄土高原水土保持区、青藏高原江河水源涵养区、西南石漠化防治区、重要森林生态功能区等五大重点生态区,退耕还林、退牧还草、自然保护区建设和水土流失治理等生态治理区域面积到2020年将分别达到151.1×104km2、27.6×104km2、88×104km2、26.7×104km2、106.2×104km2的面积。同时,根据石漠化综合治理工程“十三五”规划,到2020年岩溶土地和石漠化土地治理面积分别达到5×104km2和2×104km2以上,建设和保护林草植被195×104hm2[60]。我国的森林覆盖率从1998年的13.92%提升至2015年的21.63%[61],预期到2025年将增加至24.1%[62]。我国《十四五规划和2035年远景目标纲要》指出,要坚持山水林田湖草治理,通过在重点生态功能区、国家级自然保护区等地实施生态系统保护和修复工程,改善生态系统质量。这些工程将布局在“三生态区-四森林带”,预期在青藏高原生态屏障区(如三江源、祁连山)治理退化草原320×104hm2,在黄河重点生态区(如黄土高原、秦岭)治理水土流失和沙地280×104hm2,保护林草地80×104hm2,在长江重点生态区造林110×104hm2,治理水土流失和石漠化面积共600×104hm2。同时将在东北森林带、北方防沙带和南方丘陵带造林约300×104hm2,并且预期将在北方防沙带治理沙化土地750×104hm2。这些积极的生态保护工程预估将进一步增强生态脆弱区和重点生态功能区的生态系统结构与稳定性,提升生态系统功能和质量。受全球变暖影响,预计未来三十年(2021—2050)中国大部分地区将遭受更严重、更频繁和持续时间更长的干旱,尤其在西北地区[63],而2040年后中国降水整体呈增加趋势(集中在北方地区)[64]。因此未来几十年降水能否承载起生态工程导致的森林增加和生长的需求是个很大的挑战,未来我国半干旱半湿润区的生态工程的布设应关注具有变湿趋势的地区。

本研究通过考虑参照系的生态系统质量变化,在一定程度上消除不同区域气候波动差异对质量评估的影响,为量化保护、治理等生态工程对生态系统的影响提供了一种新途径。例如长白山参照系无土地利用变化,其ESQ可能会随气候限制因子的消极变化而下降(-0.051),而整个区域的ESQ能抵消气候变化的不利影响进而使生态系统质量有较大的提升(+0.110)。ESQ′在量化生态工程等人为影响方面具有较大潜力,可以作为政府部门绩效评估的依据之一,也能为生态系统和土地利用管理提供参考。

参照系的合理选取是评价生态系统质量变化的重要前提。本研究中大部分评估单元的参照系符合优良生态系统的假设,然而受到数据获取的限制,参照系选取主要基于CERN定位观测台站永久样地和国家自然保护区核心区,尚不能充分利用国家林草局、国家气象局、生态环境部、高校等其他部门野外观测台站的永久样地数据。因此,仅考虑极小干扰参照系很难覆盖全部的评估单元,本文根据数据的可获取性选择最受气候影响区和评估单元优良值基准作为剩余评估单元的参照系,基于多种基准参照系的基础上对国家尺度的生态系统质量进行评估；同时由于数据观测和处理、模型模拟过程中存在的随机误差和系统误差[65-67],导致部分评估单元的参照系并不完全符合优良生态系统的假设,一些功能处于评估单元内的中等水平。未来国家不同部门野外观测研究站网规模的扩大和数据共享,以及中国生态保护红线地面观测网络的建立,形成一体化、长期、连续的生态系统质量监测体系,将为我国生态系统质量评估提供更多可靠的参照系数据、模型验证和校准数据,使基于参照系的生态系统质量评估结果更具可比性和科学性,从而更好地服务于我国生态文明建设、重点生态区综合治理、石漠化综合治理、“山水林田湖草”等工程的治理成效评估和绩效考核。

4 结论

本文基于1990—2015年中国土地覆盖、植被覆盖度、生态系统功能和气候等数据,计算了生态系统结构、功能、稳定性和人类胁迫指数；根据生态地理分区和植被类型划分的评估单元,基于CERN台站、国家自然保护区和最受气候影响区等参照生态系统,结合生态系统质量评估模型,评估了基于参照系的中国陆地生态系统质量及其变化情况。主要研究结论如下：

(1)本研究选用的生态系统质量评估参照系具有较低的人类胁迫和良好的生态系统功能和稳定性,参照系的生态系统质量指数总体位于评估单元89%分位数。1990—2000年和2000—2015年参照系ESQI的全国平均值分别呈下降和上升趋势。与参照系相比,中国陆地生态系统质量指数的平均值偏低28%。

(2)1990—2000年和2000—2015年ESQ′分别下降1.2%和上升3.69%。ESQ′转变比较明显的区域主要分布在“陕-甘-宁”地区、长白山地区和西南石漠化治理工程实施县域,在这些区域实施的退耕还林还草工程、三北防护林工程和石漠化综合治理工程发挥了重要作用。未来我国计划进一步实施生态保护和修复工程,加强山水林田湖草沙治理,这些措施将进一步提升我国生态系统质量。同时需要关注未来气候变化以及生态系统功能之间的权衡关系对生态系统质量的影响。

(3)针对传统评估方法因包含了地理、气候本底条件差异而存在评估结果在空间上不可比的问题,基于参照系的生态系统质量评估方法通过考虑参照系的生态系统质量变化,可以在一定程度上消除不同区域气候波动差异对质量评估的影响,为科学量化保护、治理等生态工程对生态系统的影响提供了一种新途径。未来应加强全国不同区域、不同类型的参照生态系统联网观测,为提升我国生态系统质量评估能力提供科学数据支撑。