生态系统服务指生态系统与生态过程形成及维持的人类赖以生存的自然环境条件与效用, 分为供给服务、调节服务、支持服务和文化服务等^[1-3]。土地作为人类生产、生活的载体, 一方面土地利用模式直接影响生态系统提供的服务种类与供给程度; 另一方面, 土地覆盖变化也会驱动生态系统服务功能的改变^[4-6]。有关土地利用与生态系统服务关系的研究主要聚焦于两方面：其一为利用土地覆盖类型量化生态系统服务价值(ESV, ecosystem services value), 欧阳志云^[7]、谢高地^[8-10]以Costanza^[3]等人研究成果为基础, 构建了生态系统服务评估单价体系, 之后学者引入月度净初级生产力、物价修正系数、生态干扰度等指标, 以土地利用为基础的生态系统服务价值量化框架不断被完善^[11-15]; 其二为评估土地利用变化与生态系统服务的动态响应, 如土地利用/土地覆盖变化与气候调节^[16-17]、固碳释氧^[18-19]、水土保持^[20-23], 以及娱乐和美学价值^[24]等生态服务功能存在高度相关性, 依据土地退化面积、土地利用空间格局的改变估算生态系统服务变化量, 为土地管理与生态空间管控提供指导^{[6, 25-30]}。

探究土地利用变化对生态系统服务的影响可以作为连接生态过程和城市化等人类活动进程的手段^{[27-28, 30-32]}, 生态恢复效果对区域生态系统服务空间格局的评价与预判起重要反馈作用。当前有关生态系统服务功能评价的研究往往只利用截面数据, 生态系统服务价值的改善仅作为土地管理政策的提出依据独立于评价方法, 形成一些生态管控措施与评价结果不匹配的现象。基于此, 本研究以内蒙古自治区为研究区, 提出一种将生态系统服务价值改善效果纳入生态服务功能评价框架的分级方法, 使评价结果可以更为有效地引导和协调土地资源合理配置与生态空间管控。具体而言, 本研究的主要目标包括：(1)估算研究区土地利用变化对静态生态系统服务价值的影响, 对2005—2015年内蒙古各类型生态系统服务的动态改善效果以及各盟市生态系统服务价值的空间差异进行定量分析; (2)引入生态系统服务价值相对变化率重构生态系统服务价值空间格局, 依据分级结果识别生态改善效果与生态系统服务高、低值区对应关系, 提出区域生态改善建议。

1 研究区概况

本文以内蒙古自治区为研究对象, 从省域与市域尺度评价生态脆弱区的生态系统服务价值。内蒙古自治区地域辽阔, 位于37°24′—53°23′N, 97°12′—126°04′E之间, 呈狭长状横贯我国东北、华北和西北地区, 土地面积约118.3万km²; 地势由南向北、西向东倾斜, 大部分地区的海拔位于1000—1500 m(图 1); 气候以温带大陆性季风气候为主, 从东至西依次涵盖温带湿润区、半湿润区、半干旱区、干旱区和极端干旱区; 年降水量约30—500 mm, 由东北向西南递减^[24]; 地貌类型呈现出高原、山地、平原、高原带状分布的格局; 土地利用类型存在地域分异, 东部地区以林地与耕地为主, 中部与南部地区分布着草地, 西部则以裸地为主。

作为我国北方典型的农牧交错带区域, 内蒙古土地利用政策从大面积开垦草地和林地以满足经济发展需求逐步过渡为现阶段“生态优先”的发展路径^[33], 土地利用与生态系统格局均发生了较大改变, 通过退耕休牧、还林还草、风沙源治理、生态移民、封育轮牧等生态管控措施, 内蒙古草地退化与土壤侵蚀有了一定程度的改善^{[24, 33-35]}。

2 数据与研究方法 2.1 数据来源

本研究使用的2005年、2010年与2015年内蒙古土地利用数据来自中国科学院地理科学与资源研究所, 地类分为耕地、林地、草地、水域、城乡工矿居民用地、未利用土地6个一级类型以及25个二级类型, 空间分辨率为30 m×30 m, 精度达到90%以上; 高程信息来自美国航天飞机雷达地形测绘数据(SRTM, http://srtm.csi.cgiar.org); 用于生态系统服务价值核算的农作物播种面积、产量数据和CPI指数来自《内蒙古统计年鉴(2016)》; 内蒙古农作物净利润指标来自《全国农产品成本收益资料汇编》。

2.2 研究方法 2.2.1 静态生态系统服务价值

采用当量因子法评估内蒙古及各盟市的ESV, 步骤为标准当量价值量的确定、不同土地类型下单位面积价值表的构建和静态价值计算。

(1) 标准当量的确定与修正：1单位标准当量的价值量指1 hm²平均产量的耕地每年粮食产量的货币价值^[8-10], 一般将单位面积耕地稻谷、小麦和玉米的净利润作为1个标准当量因子的生态系统服务价值量^[10]。以往研究对ESV的量化仅使用静态指标, 考虑到年际影响与融入ESV改善时存在的价格变动, 本研究使用CPI指数对多年净利润平均值确定的标准当量进行修正, 方法如下：

(1)

式中, D为1个标准当量因子的生态系统服务价值量(元/hm²); S_ri、S_wi和S_ci为第i年稻谷、小麦和玉米的播种面积占三种作物总面积百分比; F_ri、F_wi和F_ci为第i年全国稻谷、小麦和玉米的单位面积平均净利润(元/hm²); α_i为第i年相较于基年的CPI修正系数。经计算, 内蒙古自治区2005—2015年平均价值量为2725.19元/hm², 物价修正后D值为1587.01元/hm²。

(2) 不同土地利用类型的静态ESV：以研究区2005年、2010年和2015年3期土地利用数据为基础计算ESV变化量。将内蒙古生态系统服务分为供给服务、调节服务、支持服务和文化服务共4类11子类。根据研究区自然环境, 结合谢高地^[10]制定的中国陆地生态系统当量表对生态系统服务价值系数进行修正。其中, 耕地使用农田生态系统当量因子; 研究区有林地、疏林地和其他林地处于欧亚针叶林植物区, 同时受欧亚草原植物区和东亚阔叶林植物区的影响^[36], 取针叶、针阔混交与阔叶生态系统当量因子的平均值; 草地方面, 高覆盖度草地分布于大兴安岭东西两麓, 以林缘草甸与草甸草原为主, 选择草甸草地当量因子, 中、低覆盖度草地选择草原草地当量因子^[37-38]; 土地利用为河渠、湖泊和水库坑塘的选择水域当量因子, 滩地与沼泽地选择湿地当量因子; 沙地、戈壁、盐碱地和其他用地的选择荒漠当量因子, 裸土地和裸岩石质地选择裸地当量因子; 城镇用地、农村居民点与其他建设用地的当量因子设为0。上述土地利用类型下的生态系统服务当量因子与D值相乘得到内蒙古自治区单位面积生态系统服务价值系数表(表 1)。计算研究区的生态系统服务价值如下式：

(2)

式中, A_k为土地利用类型为k的面积; VC_k为单位面积生态系统服务价值系数; ESV^(f)为第f项生态系统服务功能价值; VC_fk为第f项生态系统服务功能单位面积生态系统服务价值系数。

2.2.2 生态系统服务价值相对变化率

本研究提出双维度下的ESV相对变化率计算方法：其一为区域维度上的相对变化率R_i; 其二为单项生态服务功能维度上的相对变化率R^(f), 两指标用于分析空间和功能尺度上ESV的增减趋势, 以反映生态改善效果, 公式如下：

(3)

式中, R_i为区域i的ESV相对变化率, L_t+j、L_t分别为区域i在t+j与t期的局部ESV, C_t+j、C_t分别为t+j与t期的总体ESV。R^(f)为第f项ESV相对变化率, 其他指标解释同式(2)。|R|＞1表示区域i或第f项生态服务功能的ESV变动幅度大于研究区整体变化幅度; |R|≤1表示区域i或第f项生态服务功能的ESV变动幅度小于等于研究区。R值的正负表示区域i或第f项ESV变化方向与研究区变化方向的一致性。

2.2.3 生态服务功能分级评价

研究区生态系统服务价值格局的研究单元应反映该单元多类生态系统的ESV, 采用ArcGIS块统计对90 m×90 m内的各单项ESV进行求和, 基于R^(f)对研究区各单项ESV的修正, 得到各研究单元ESV改善情况, 公式如下：

(4)

式中, ESV_i^f为第i个30 m×30 m土地栅格内土地利用类型对应的生态系统服务价值; ESV^*为此90 m×90 m研究单元在考虑相对变化率后的生态系统服务价值。

经R^(f)修正后, 高ESV^*值表明该地区ESV改善效率较高, 当前生态管控措施下生态恢复效果显著; 低ESV^*值表明该地区生态系统服务的提升能力较差, ESV改善效果甚微, 需调整地区的生态管控措施。因此, 将ESV^*指标纳入生态服务功能评价框架可实现对部分高ESV且生态改善趋势明显, 以及低ESV且生态不断恶化地区的关注, 有利于决策者协调各地区的生态管控强度, 实现对同一ESV地区内等级的细分。具体方法为利用分位数断点对供给、调节、支持、文化与ESV改善5项数据分级赋值, 形成低、较低、一般、较高、高五个等级。

(5)

由上述5项分级结果, 利用式(5)进行综合评价。Score_i为i单元ESV等级值; k_i^(p)、k_i^(r)、k_i^(s)、k_i^(c)、k_i^(im)分别为供给、调节、支持、文化与生态恢复效果等级值。利用ArcGIS空间分析, 将Score_i重分类为低与较低区间内的土地单元归类为低生态价值区, 高与较高区间内的归类为高生态价值区, 形成基于不同生态服务功能与生态恢复效果的综合分区。

3 实证结果 3.1 生态系统服务价值时空变化

内蒙古ESV的绝对量方面, 气候调节与生物多样性提供的生态价值较高, 水资源供给与维持养分循环相对较小。2005—2015年间, 内蒙古生态系统服务价值实现增长, 从2005年的17069亿元增长至2015年的18667亿元, ESV变化率达9.36%(表 2)。其中, 2005—2010年生态系统服务供给能力轻微减少, 体现在净化环境功能、水文调节功能提供程度的下降; 2010年以来, 除气候调节功能外, 其余生态系统服务价值呈增加趋势, 水资源供给、水文调节、生物多样性维护和美学景观功能改善明显, ESV的增加率达到15%以上。

通过土地利用转移矩阵揭示ESV变动的影响机制(表 3), 结果表明：研究区的土地退化主要体现为草地面积的减少与其他用地面积的增加。草地与未利用地的相互转化、草地向林地与水田的转化面积均占总转移面积的10%以上, 4项累计占比达60.8%。其中, 草地退化为荒漠与裸地的净面积达到25239.93 km², 土地沙漠化成为制约研究区气候调节与土壤保持的主要因素。另一方面, 研究区的生态恢复体现为水域、滩地、沼泽地面积的增加, 与旱地比重的下降, 这些土地利用变化使水资源的供给和养护能力大幅提升。城镇空间与生态空间的协调方面, 除耕地与草地外, 各类其他用地转换为城镇用地的占比均小于1%, 城镇用地的净转入面积也在所有用地类型中为次小项。结合区内生态政策可知, 随着退耕还林、退草还林等生态过程对森林保有量与生物多样性的增加, 林地与水域的生态恢复有序推进, 研究区生态系统服务价值的提升源于生态空间的增加和用地结构的优化, 城镇空间对生态空间的侵占得到较好控制, 但草地退化仍是自治区当前的主要生态问题。

3.2 生态系统服务价值改善效果分析

图 2与图 3分别显示内蒙古自治区12个盟市的R_i值, 与各单项生态服务功能的R^(f)值。就空间差异性而言, 2005—2015年内蒙古自治区的ESV变化呈现空间分布不均衡的格局。呼伦贝尔市、兴安盟与赤峰市的ESV变化与研究区整体水平一致。其中, 呼伦贝尔市ESV增长率远高于平均值。上述盟市位于大兴安岭一带, 生态空间以森林生态系统为主, 森林景观的生态恢复效果显著, 该地区的高生态服务价值带动了研究区整体生态服务功能的增强。另一方面, 内蒙古其他盟市的ESV均呈现不同程度的下降, 其中乌海市、阿拉善盟与锡林郭勒盟的R_i小于-1, 生态恶化程度大于平均值。阿拉善盟位于研究区最西部, 干旱少雨, 植被覆盖度低, 土地沙漠化现象严重, 导致草地与林地退化为裸地, 区内生态系统服务的提供能力处于下降阶段; 乌海市为高工业化水平的城镇聚集区, 且沙丘、裸岩面积较大, 土地利用体现为城镇用地对耕地、林地与草地的侵占; 锡林郭勒盟位于内蒙古自治区中部, 草场面积广大, 受草地退化影响, ESV随草场面积, 尤其是高覆盖度草地面积的减少而下降。

就功能差异性而言, 除气候调节功能外, 其他各生态服务功能的变化与总体生态服务功能变化趋势一致。究其原因, 气候调节功能受草地退化影响小幅下降, 同时土壤保持、原料生产与气体调节功能的R^(f)小于0.2, 滞后于研究区平均水平, 进一步印证生态脆弱区土地退化与水土流失的严重性。研究区水源涵养类生态服务功能如水资源供给、水文调节的ESV变化幅度较大, 表明水域的恢复和水质的养护能力提升, 是内蒙古整体生态系统服务价值改善的重要驱动力。

3.3 生态系统服务价值分级结果

2015年内蒙古单位面积生态系统服务价值的变异系数为0.69, 具有极差大、离散程度强的特征。各单项生态服务功能与生态恢复效果大致遵从由北至南, 评价等级由高至低的空间分布规律(图 4)。供给功能在研究区南部边缘形成主要的低价值区, 该区域城镇空间密集, 且位于鄂尔多斯高原与黄土高原过渡地带, 旱地、沙地交错分布, 水资源供给不足, 土地生态系统向人类提供的物质基础与服务能力较低, 为主要的低价值区; 调节功能由东北至西南递减, 主要源于东北部季风气候区植被覆盖度较高, 降雨充沛, 其净化环境、水文调节服务功能优于大陆性气候区; 西部阿拉善盟地区以沙地、裸地为主, 土壤保持能力较差, 沿黄河地带受风水混合侵蚀水土流失严重, 导致该区为支持功能的低价值区; 文化功能的高价值区集中在植被覆盖度高、自然景观分布广泛的森林区与锡林郭勒盟东部的高覆盖草场。

ESV改善效果方面, 研究区北部ESV绝对量大且改善幅度明显, 蒙西及内蒙古南部边缘等生态服务功能提供能力较差的地区改善幅度不显著。不同生态服务功能与生态恢复的空间格局具有相似性, 供给、调节、支持和文化服务功能分级与ESV改善分级的关联度较高, 即高服务功能区的生态价值对应较强的ESV改善能力, 生态价值不断巩固强化; 而生态短板地区对应低ESV改善能力, 区域生态恢复不足且不断恶化, 呈现不同区域ESV差距加大的极化结果。

基于单项生态服务功能与ESV改善分级, 形成内蒙古自治区生态系统服务价值综合评价结果(图 5)。根据空间分布特征, 划分出北部大兴安岭核心区、中部高覆盖草原区2个高价值区; 以及西部沙地区、南部沿黄河地区和东部科尔沁地区3个低价值区。高值区约占研究区面积26.3%, 低值区约占36.2%, 各地区的具体特征见表 4。

4 结论与讨论

本研究结合内蒙古自治区实际生态环境, 利用当量因子法评估了不同土地利用类型的静态生态系统服务价值; 同时, 提出将ESV改善效果纳入生态服务功能综合评价框架的分级方法。由ESV相对变化率可知, 研究区ESV变化趋势具有极差大, 高、低值区显著分化的特征, 生态服务功能空间分布不均衡。其中, 水资源养护类服务功能的增强速率远高于其他指标, 是区域整体生态环境优化的驱动力。结果显示, 虽然2005—2015年间内蒙古ESV有所增加, 水域、林地的恢复速率提高, 城镇用地扩张得到有效控制, 但北部与中部的森林与草地保育区, 以及蒙西低生态价值区的生态改善效果具有明显差异, 生态恢复在空间分布上呈现极化效应, 地区间生态环境质量的差距不断扩大, 生态保护的协调性面临较大挑战。因此, 区域生态管控措施宜进一步加大沙地区、黄河流域与东部科尔沁的生态恢复力度, 尤其调整、加强区域的草地恢复力度, 改善低植被覆盖带来的土地沙漠化问题, 缩小区域南北向生态价值的差距, 保证区域生态功能与生态安全格局的平衡。

此外, 本研究使用的生态系统服务价值评价方法仍存在缺点与改进空间, 包括：(1)ESV估计结果的准确性取决于当量因子的可靠性, 而11类土地利用类型与生态系统的对应只是实际生态过程的简化; (2)ESV分级评价中尚未考虑社会经济活动与人类资本对ESV的相互作用, 这需要对区域居民支付意愿、社会文化与景观结构等进行更全面、深入的研究。然而, 改进后的生态系统服务价值评价可以用一种简化、通用的方法体现不同生态服务功能与生态恢复效果空间格局的一致性, 在反映当前ESV高低分布特征的同时, 揭示那些低ESV地区是处于生态改善进程或是持续恶化阶段, 对ESV相同的区域进行了细分, 加强了生态服务功能评价与土地利用管理政策之间的衔接, 可以为生态管控措施的调整提供决策信息。