基于土地利用的三峡库区近30年生态系统服务价值时空变化特征

李 辉, 张晓媛, 国洪磊

(1.生态环境空间信息数据挖掘与大数据集成重庆市重点实验室, 重庆 401320;2.重庆财经学院, 重庆 401320; 3.南京师范大学 虚拟地理环境教育部重点实验室, 南京 210023)

土地资源是人类赖以生存和发展的基础[1-6]。土地利用/覆被变化(LUCC)是全球气候变化的主要原因并且与人类活动有着非常紧密的联系，因此，对土地利用或覆被变化的研究是全球十分热点的研究内容之一[7-8]。随着我国发展的需要，造成人地矛盾越来越突出，越来越重视对土地利用的研究。根据土地利用的研究，我们可以对人类活动对全球环境变化产生的影响进行定性和定量的评价，从而为资源环境的可持续开发利用提供依据[9]。生态系统服务是指通过生态系统的结构、过程和功能直接或间接得到的生命支持产品和服务，是人类社会的关键自然资本，维持和保育生态系统服务是实现可持续发展的基础[10]。土地利用/覆盖变化(LUCC)能改变生态系统的结构、过程和功能，进而影响生态系统服务[11]。当前全球所面临的生态环境危机的本质是由于生态系统服务功能被破坏的后果。目前，国内外在生态系统服务价值的评估方法上都开展过大量的研究，但统一的生态系统服务价值的评估体系尚未形成[12]。总体来说，生态系统服务价值的评估可以归纳为两大类，一是基于单位服务功能价格方法，二是基于单位面积价值当量因子方法，前者所采用的评估方法涉及的参数较多且计算繁琐，同时标准无法统一，因此在具体理论及方法上尚未形成突破[13]；后者在评估方法的选取上充分借鉴了Costanza相关研究的成果，相对于前者而言，当量因子法相对直观、数据量少、评估比较全面且可比性高，特别适用于区域和全球尺度生态系统服务价值评估。我国学者谢高地等[14]在Costanza评价模型基础上改进并计算得到中国生态系统服务价值的基础当量，被大多数学者和研究人员借鉴并广泛应用于评估国家、省市县、城市群、流域等不同空间尺度的生态系统服务价值。目前，对于土地利用/覆被变化和生态系统服务价值之间关系的研究相对较少，如韩增林等[15]运用地统计学方法及地理加权回归模型对辽东湾北部区近30 a的土地利用强度与生态系统服务价值的时空关系进行了分析；但从长时间序列和土地利用角度对三峡库区生态系统服务价值进行全面研究还相对较少。

三峡库区拥有十分独特而脆弱的生态系统，是长江经济带乃至全国重要的生态宝库，在当前及今后相当长的时期内三峡库区的生态环境在长江经济带的发展中将始终占据着十分重要的地位[16]。随着长江上游地区经济快速发展和三峡大坝的建设及蓄水，使得库区土地利用格局发生了较大变化，大量生态用地被淹没，建设用地规模不断扩大，使得区域的生态系统服务价值也发生了一定的变化。本文利用RS和GIS技术，提取三峡库区1986年、1995年、2000年、2007年、2010年和2018年的各类土地利用的数量以及空间演变情况，进而对研究区1986—2018年的土地利用变化规律进行时空演变的分析，并通过相关模型对研究区的生态系统服务价值的时空演变特征进行探究和分析，为库区的可持续发展提供参考。

1 研究区概况

三峡库区的地理位置位于东经106°20′—110°30′，北纬29°00′—31°50′(长江中上游地区)，并与四川盆地的相互连接，南北两面分别为川鄂高原和大巴山，其包括巴南区、江津区、长寿区、石柱县、巫山县、渝中区、武隆县、巫溪县、北碚区、丰都县、奉节县、沙坪坝区、南岸区、云阳县、渝北区、涪陵区、开县、九龙坡区、大渡口区、江北区、万州区和忠县22个重庆市所辖的区县以及湖北省的4个区县(夷陵区、秭归县、兴山县和巴东县)(图1)。亚热带季风气候是三峡库区的主要气候类型，每年各个地区的降水量几乎相同，总量为1 000～1 800 mm。在库区内水土流失严重，滑坡、泥石流、崩塌以及地震经常发生。随着三峡工程的修建，库区蓄水，从而改变了库区的气候环境，扩大了水域面积，减小了日夜温差。因此，对三峡库区土地利用变化特征的研究具有重大意义。三峡库区幅员面积大，其土地资源的总量大而耕地少，人均耕地严重不足。研究区以紫色土、粗骨土、石灰土以及水稻土作为主要的土壤类型，植被种类多样，自然资源储存量丰富。三峡库区土地类型相对单一，主要分类类型为山地及丘陵，山地面积占库区总面积的76%，丘陵其次，占库区总面积的15%，地形起伏大、地貌类型多样，林地在农用地中占比较大，耕地以坡耕地为主，零散分布，旱地居多；可利用率的建设用地很低，基础设施的建设成本高、投资大，各土地利用类型分布特点各异，土地利用水平空间上的差异较大，是典型的大库区、大山区。

图1 研究区概况及评价单元格划分

2 数据源和研究方法

2.1 数据源

文章以轨道号为125038，125039，126038，126039，126040，127038，127039，127040，128039，128040的三峡库区1986年、1995年、2000年、2007年、2010年和2018年的Landsat TM/Landsat 8 OLI遥感影像作为原始数据源(从影像的获取质量选取年份)。对于所获取的所有遥感图像均以ALBERS投影进行投影，使用中央经线为东经110°，分别使用北纬25°和47°作为双标准纬线，其中获取的三峡库区TM影像数据的1—5及7波段的空间分辨率为30 m，第6波段的空间分辨率为120 m。文中所涉及到的地理空间矢量行政区数据的来源为《地理国情监测云平台》(http:∥www.dsac.cn/)，其他相关数据源自研究区历年统计公报，如重庆市统计年鉴数据和湖北省统计年鉴数据和各年全国农产品成本收益资料汇编数据。

2.2 数据处理

由于遥感传感器自身就具有一定的误差，再加上辐射受大气的影响，所以对获取到存在几何畸变和辐射畸变等问题的遥感数据进行了一定的遥感技术处理，如对遥感影像进行判读前，需要使用Erdas 9.1对遥感影像进行适当的校正(辐射和几何)以及镶嵌处理，并对影像进行增强处理，使影像更为清晰，解译标志更为明显[17-18]。而后对遥感影像进行解译，得到需要的土地利用现状数据。本研究将计算机图像处理和目视解译的方法结合起来对1986年、1995年、2000年、2007年、2010年、2018年的TM遥感影像图进行解译，从而获得对应年份的土地利用现状数据。

2.3 研究方法

2.3.1 土地利用数量变化研究 本文从研究区各类土地利用的数量变化和空间分布状况的演变以及土地利用程度的变化等情况进行量化分析，以揭示区内土地利用特征[19-22]。引入GIS等手段，运用软件空间统计分析功能、制图功能等将三峡库区1986年、1995年、2000年、2007年、2010年、2018年的土地利用现状数据进行分类统计，并对在时间与空间上对其进行分析。

2.3.2 土地利用动态指数模型 文章引入单一土地利用动态度和综合土地利用动态度，并运用土地利用动态指数模型，分别对研究区各类土地利用类型的数量变化和速度的变化进行分析，同时对整个研究区各类土地利用的数量和速度变化进行研究分析，土地利用动态指数模型如下[23-24]：

(1) 单一土地利用动态度

(1)

式中:K为某一种土地利用类型在研究年限内的动态度;ua,ub为要研究的起始期以及研究的结束期某一种土地利用类型的数量;T为研究的时间跨度(单位为a);K为某一土地利用类型在研究时间跨度内的变化速率[25]。

(2) 综合土地利用动态度

(2)

式中:Lui为起始期第i类土地利用类型的面积;ΔLui-j为研究期内第i类土地利用类型转化为非i类土地利用类型的面积的绝对值;T为研究的时间跨度(单位为a);Lc为整个研究区的土地利用综合变化速率。

2.3.3 土地利用程度指数模型 土地利用的广度和深度主要以土地利用程度进行测算[26]，体现出土地收到其自然属性以及人类活动的双向影响程度。本文引入刘纪远等[27]提出的土地利用程度综合分析方法，并充分结合三峡库区的现实情况，将研究区的土地利用程度划分为4个等级，土地利用程度的分级指数见表1，可以通过计算得到研究区的土地利用程度指数，具体计算公式如下:

(3)

式中:L为整个研究区的土地利用程度指数;Ai为研究区内第i级土地的土地利用程度分级指数;Ci为研究区内第i级土地面积占整个研究区面积的百分比;n为土地利用程度的分级数。

Ci=CCi/HJ

(4)

式中:CCi为第i级土地利用类型的总面积;HJ为整个研究区土地的总面积。

(5)

式中:Ia,Ib分别为研究初始期及研究终止期的土地利用程度指数;Ai为第i级土地利用程度分级指数;Cib,Cia分别为研究初始期及研究终止期第i级土地利用程度的面积。

表1 土地利用程度指数

2.3.4 生态系统服务价值核算模型 本文借鉴谢高地等[14]2007年所修订的生态系统生态服务价值当量核算因子表，标准当量为单位(1 hm2)面积农田生态系统粮食生产的净利润作为1个标准当量因子的生态系统服务价值量。研究区以旱地和水田作为主要的耕地类型，大部分种植水稻及玉米；有林地、灌木林地和针阔叶混交林为主要的林地类型；盐碱地和滩涂为未利用地主要类型。本文研究尺度相对较小，区域自然地理特征相对一致，且农业生产方式无大差异，生态系统属性基本一致，能够运用当量赋值的方法对区域生态系统服务价值进行核算。在借鉴前人研究成果的基础上，充分结合三峡库区真实状况对标准当量的系数进行适当的修订，具体修正模型如下：

ESVk=∑Ak×VCk

(6)

(7)

式中: ESVk和ESV分别为研究区第k类型的服务价值和总服务价值;Ak为第k类型的总面积;VCk为价值系数,是第k类型单位面积的生态系统服务价值。

本文所运用的基础数据来源于中国统计年鉴数据和重庆市统计年鉴以及湖北省统计年鉴，并参考了全国各年农产品成本收益相关资料的汇编，部分数据在参看数据源的基础上，通过一定的计算获得。根据三峡库区的真实状况，借鉴和参考谢高地等[28]的相关研究成果，修正了研究区各类型的生态系统的服务价值系数，以三峡库区1986年、1995年、2000年、2007年、2010年以及2018年6年的粮食平均产量(4 468 kg/hm2)以及平均价格(1.46元/kg)为研究区基准单产和价格基准，核算出了研究区生态系统服务价值当量因子，具体数据为931.90元/hm2，利用公式(6)及公式(7)进行计算，得到三峡库区1986—2018年的生态系统单位面积生态服务价值的平均值(表2)。

表2 1986-2018年三峡库区生态系统单位面积生态服务价值平均值 元/hm2

3 结果与分析

3.1 三峡库区土地利用类型结构分析

在RS和GIS的技术支撑下，对1986年、1995年、2000年、2007年、2010年和2018年三峡库区土地利用现状数据进行分类统计，进而得到各年对应各类土地利用类型面积和比例，结果见表3；运用ArcGIS 10.1软件的制图功能，得到了三峡库区1986年、1995年、2000年、2007年、2010年和2018年6期的土地利用类型空间分异图(图2)。

表3 1986-2018年三峡库区各土地利用类型面积及占比

图2 1986-2018年三峡库区土地利用类型空间分异

通过表3和图2我们可以看出三峡库区1986—2018年土地利用的情况，总体来看，三峡库区土地利用的主要类型是林地和耕地，耕地在库区的中上游段分布居多。库区的下游段和重庆都市区的4条平行山脉(缙云山、中梁山、铜锣山和明月山)以林地为主要的土地利用方式；草地零星地分布在三峡库区的部分区域；建设用地规模随着时间地变化不断变大。在研究年限内，水域、建设用地和未利用地面积都呈现增多的状态；草地、耕地和林地呈现减少的状态，三峡库区用地类型以林地与耕地为主。在近30 a的时间里，建设用地是三峡库区所有土地利用类型中变化幅度最大的，其面积由234.59 km2增长到了1 861.99 km2，增加了2.84%，并呈现不断增加的趋势；其次是草地，面积由2 306.04 km2减少到1 396.06 km2，减少了1.59%；第三是耕地，面积由22 141.89 km2减少到21 340.44 km2，减少了1.4%，并呈现逐年减少的趋势；水域面积从764.86 km2增加到了1 513.67 km2，这是由于三峡大坝工程的建设和使用，导致库区水面面积不断增加；林地面积由31 879.08 km2减少到31 206.16 km2；未利用地占三峡库区总面积的比例本身比较小，变化不明显，但在总量上有一定的增幅，其具体变化过程分为4个阶段：先增加后减少，再次增加，再次减少。鉴于此，三峡库区的国土管理部门需要有针对性的制定相关耕地和林地的保护政策并切实推行，以保障库区的粮食和生态安全。

3.2 三峡库区土地利用变化速度分析

引入单一土地利用动态度及综合土地利用动态度对三峡库区土地利用变化速率进行分析，在利用1986年、1995年、2000年、2007年、2010年和2018年三峡库区土地利用现状分类数据的基础上，采用土地利用动态度的计算公式获得各土地利用类型的单一土地利用动态度及综合土地利用动态度，见表4。

从表4可以看出，三峡库区近30 a来建设用地的土地利用变化速度最快，在1995—2000年时段变化率最明显，增长率达到了38.596%；耕地在2010—2018年时段变化最明显，减少率为0.259%；林地、水域、草地和未利用地在1995—2000年时段变化最明显，林地减少率、水域增加率、草地减少率和未利用地增加率分别为0.157%，7.445%，7.674%，18.05%。

研究区土地利用综合动态度最大的是1995—2000年时段，综合动态度为0.421，由此可知，在整个研究期内，三峡库区土地利用的变化速度最快的时期在1995—2000年时段，变化率为0.421%；在这一时段水域、建设用地和未利用地都以较快的速度增长，而耕地、林地以及草地都以不同程度的变化速度减少，这是由于1997年重庆直辖，加快了三峡库区的发展。三峡库区在2010—2018年时段土地利用变化速度相对较快，变化率为0.175%，由此可见，建设用地和水域的增长速度加快，而其余类型的地类都以不同程度的速度减少。三峡库区在2007—2010年时段和在1986—1995年时段土地利用变化率分别为0.099%，0.077%，水域、建设用地和未利用地都以不同程度的变化速度增长，耕地、林地以及草地都以不同程度的变化速度减少；三峡库区在2000—2007年时段的土地变化率为0.052%，水域和建设用地都以不同程度的变化速度增长，耕地、林地、草地以及未利用地都以不同程度的变化速度减少。

综上分析，水域和建设用地的土地利用变化率在5个时段均为正值，表明研究区近30 a来水域和建设用地的面积不断增加；耕地、林地和草地的土地利用变化率在5个时段均为负值，表明其耕地、林地和草地的面积不断减少；而未利用地的面积变化在这近30 a的5个时段中呈先增加后减少再增加又减少的趋势。在1995—2000年时段三峡库区土地利用变化率达到峰值，说明在该时段三峡库区城乡规划建设力度较大。

3.3 三峡库区土地利用程度分析

通过对研究区土地利用程度的计算和分析，能够进一步认识土地利用变化的发展程度[11]，分析发现研究区土地利用是否处于发展状态。利用1986年、1995年、2000年、2007年、2010年和2018年三峡库区土地利用现状数据和土地利用程度的核算公式，本文得到了研究区各土地利用类型的土地利用程度和三峡库区近30 a的土地利用程度综合指数以及三峡库区近30 a来土地利用程度变化量，见表5。

表5 三峡库区1986-2018年各类土地利用程度综合指数

三峡库区在1986—2018年里耕地的土地利用程度是最高的，未利用地的土地利用程度是最低的；水域和建设用地的土地利用程度在1986—2018年近30 a里呈现逐年上升的趋势，而耕地、林地和草地的土地利用程度在1986—2018年近30 a里呈现逐年下降的趋势。从总体来看，未利用地的土地利用程度在1986—2018年期间呈现上升的趋势，但其实际变化不明显，这是由于未利用地面积本身基数较小造成的。

研究区在1986—1995年的土地利用程度变化量为0.003 1，1995—2000年的土地利用程度变化量为0.027 7，2000—2007年的土地利用程度变化量为0.001，2007—2010年的土地利用程度变化量为0.001 8，2010—2018年的土地利用程度变化量为0.009 1，由于该5个时段的土地利用程度变化量都大于0，因而在这5个时段研究区都处于发展期。由于研究区1995—2000年的土地利用程度变化量是最大的，在一定程度上可以说是发展速度最快，这是由于1997年重庆直辖加速了三峡库区的发展。

3.4 生态系统服务价值时空演变分析

在GIS技术的支持下，结合三峡库区生态系统单位面积生态服务价值表，得到了1986—2018年三峡库区生态系统服务价值核算结果以及结构变化情况表(表6)。通过对表6分析可知，研究区1986—2018年生态系统服务价值总量由1 639.56亿元持续增长到1 662.01亿元，增长率1.37%；单位面积生态系统服务价值由28 596.02元/hm2增长到28 987.57元/hm2。研究区各土地利用类型的生态系统服务价值占比大小依次为:林地>耕地>水域>草地>未利用地>建设用地；研究期内耕地、林地和草地的生态系统服务价值持续减少，变化率分别为-3.62%，-2.11%，-39.45%；水域和未利用地生态系统服务价值持续增加，变化率为97.91%，0.01%；耕地生态系统服务价值减速逐渐降低；林地、草地和水域生态服务价值呈现先增后减趋势；未利用地生态服务价值逐年增加。

表6 1986-2018年三峡库区各地类生态系统服务价值结构及变化

本文通过计算评价单元格单位面积生态系统服务价值，引入中心点赋值法，运用GIS软件进行生态系统服务价值空间插值，结合几何间隔分类，采用GIS自然间断点法，并充分考虑研究区的真实状况，将研究区生态系统服务价值评价单元的价值划分为低[0，2000)、较低[2000，4000)、中[4000，6000)、较高[6000，8000)、高[8000，∞)共5个价值等级，并得到研究区6期单位面积生态系统服务价值空间格局分布图(图3)。

图3 1986-2018年三峡库区生态系统服务价值空间格局演变

从图3可以看出，三峡库区单位面积生态系统服务价值主要以较高和高两个等级为主，其面积平均占比为53%，32.53%；高和中等级面积比例呈现减少趋势；低、较低和较高生态服务价值面积占比呈现增加趋势，所有等级的面积变化速率为先增后减；东北部、东南部山地以及长江沿岸区域单位面积生态系统服务价值较高，主要为高级和较高等级分布区域，江津区、长寿区以及忠县是中等级生态系统价值的主要分布区，重庆主城区和渝西地区是较低和低级生态系统服务的集中区域。

3.5 生态系统服务价值空间自相关分析

从表7空间自相关分析结果显示，研究区1986-2018年的6个时期的生态系统服务价值全局Moran′sI值全部大于0，且p值全部低于0.001，说明三峡库区生态系统服务价值的空间分布呈现出一定的正向自相关性，同时在空间上呈现出显著的聚集性，其中高值区聚集明显，中低值区分区域相邻布局；研究区1986-2018年6个时期的全局Moran′sI指数总体上呈现上升的趋势，在1986年的值最低，为0.012 8，表明生态系统服务价值的空间自相关性持续增强，由于库区的建设与区域发展等因素促使三峡库区生态系统服务价值的空间聚集性持续加强。

表7 1986-2018年三峡库区ESV空间自相关性

3.6 土地利用程度对不同生态系统服务类型的影响分析

为了定量验证并进一步分析土地利用与生态系统服务价值的相关关系，本文计算出了研究区1986年、1995年、2000年、2007年、2010年以及2018年6个时期的土地利用程度综合指数分别为2.394 2，2.397 3，2.425 1，2.426 1，2.427 8，2.436 9，再核算研究区土地利用程度综合指数与生态系统服务价值的相关系数，同时在Origin 2017软件中生成二者关系的散点图，并进行线性拟合，结果见图4。

图4 三峡库区土地利用程度与生态系统服务价值相关关系

研究区土地利用程度综合指数与生态系统服务价值的相关系数为0.815 7，表示二者关系为强相关，且拟合曲线显示研究区生态系统服务价值随土地利用程度的增加呈上升趋势(图4)，说明二者存在显著的正相关关系。

土地利用程度的加深表明人类对土地的劳动及资本投入增加，同时表明人类对生态系统干扰也越多。研究区土地利用程度的增加主要是由于库区社会经济的快速发展，催生了对大量建设用地的需求，而建设用地的主要来源是耕地和林地，使得这两种地类有所减少，但由于“天保工程”，“退耕还林”等生态保护政策的有效实施，又使得林地和耕地在大量转化为建设用地的同时有一定的补偿来源。另一方面，由于库区的三次蓄水，使得库区水平面有所上升，水域面积增加较为显著。在生态系统服务价值的构成中，林地和耕地的价值始终低于水域的价值，林地和耕地的减少导致生态系统服务价值的降低要低于水域增加的价值，最终使得研究区生态系统服务价值增加。

4 讨论与结论

4.1 讨 论

(1) 本文研究方法存在一定的局限性。生态系统服务价值的评估存在许多主观因素，同一区域在不同的评价方法、研究视角以及参数选择等方面都会存在差异，从而影响评估结果，但生态系统服务价值的变化趋势不会有较大差异，不会对研究结果的比较产生影响。如王大菊等[29]采用替代花费法、影子工程法、碳税法、市场价值法和水量平衡法对三峡库区的生态系统服务价值进行量化评估的结果为2 027.59亿元(2015年)；齐静等[30]采用影响生态系统质量状况的生态参数修订当量因子并建立生态系统服务价值(ESV)估算模型对三峡库区2000—2015年ESV进行了核算，其结果为和3 781.16亿元。本文采用修订后的单位面积价值当量因子法评估了三峡库区1986—2018年近30 a的生态服务价值，得到的结果为2 106.6～2 185.7亿元，与王大菊等[29]的评估结果相接近，但与齐静等[30]相差较大。因此，参数修订当量因子法也存在一定的不足，相同研究区域修正方法和参数的不同都会使评价结果产生差异，而且局部的差异也难以用统一的参数来衡量。因此，迫切需要有精准性更好地评估模型来提高生态系统服务价值量化评估的准确性。

(2) 土地利用变化对生态系统服务价值影响。三峡库区作为我国重要的水源和生态涵养区，其30 a来土地利用发生剧烈变化。林地是研究区最主要的生态系统类型，占50%以上，林地面积逐年降低，在研究期始末降至54.74%；建设用地面积增加迅猛，由1986年的234.59 km2增至2018年的1 861.99 km2，增长率高达693.76%。土地利用主要表现为耕地、林地转为建设用地和水域的变化特征，这也使得生态系统服务价值的协同改变。

(3) 通过对土地利用程度与生态系统服务价值的时空耦合关系分析可知，土地利用程度的增加在短期内会带来显著的经济效益，由于本文研究年限跨越三峡库区的蓄水，才使得生态系统服务总价值未呈现减少，但长此以往，土地利用程度的无限增加会使得生态系统服务价值减少，同时造成的生态经济损失将难以估量。鉴于此，在对土地进行投资和使用时，慎重决定投资方式及建设力度，协调好各种土地利用类型的比例，构建完善的土地利用类型结构，尤其是研究区的渝西区域和湖北4县地区。

(4) 生态系统服务价值的变化研究是区域可持续发展相关研究的热点内容之一，本文以三峡库区为研究对象，对其30 a来土地利用格局时空变化特征和生态系统服务价值进行了评估和分析，研究可为三峡库区国土开发政策制定和核算生态系统服务价值提供参考。但由于生态系统服务价值研究还处于发展阶段，相关理论体系和研究方法有待进一步完善和提升，如何更准确地评估区域生态系统服务价值以及分析其与土地利用之间的内在联系和变化机理将有待于进一步探讨。

4.2 结 论

(1) 林地和耕地是研究区最主要的土地利用类型，土地利用类型变化速度最快的是建设用地，最慢的是林地；建设用地、未利用地和水域呈现增加的状态，草地、耕地和林地正好相反，呈现减少的状态，建设用地、未利用地和水域以及草地、耕地和林地的变化程度和变化过程又存在着明显的差异。因此，继续保持耕地和林地作为研究区主要地类的优势地位依然重要。

(2) 三峡库区土地变化速率均为正值的是水域以及建设用地，建设用地和水域用地面积持续增加；草地、耕地和林地正好相反，它们的变化速率为负值，耕地、林地和草地用地面积持续减少，这需要引起相关部门的重视，及时制定有效政策缓解生态用地的减少。

(3) 1986—2018年，研究区单位面积生态系统服务价值主要以较高和高两个等级为主，空间分异特征比较明显，其中库区的东北部、东南部山地以及长江沿岸区域单位面积生态系统服务价值较高，重庆主城区和渝西地区是较低和低级生态系统服务的集中区域。因此，需要国土空间规划和生态环境保护等部门加强城镇建成区和扩展新区生态服务功能的补偿与维护。

(4) 1986—2018年，研究区生态系统服务价值的总量有所增加，生态安全状况总体上有一定的改善，但局部地区还需特别加强保护，从土地利用角度制定科学合理的生态和环境保护政策。

(5) 1986—2018年，研究区6个时期的生态系统服务价值的全局Moran′sI值全部大于0，体现出明显的空间正向自相关性，其中高值区聚集明显，中低值区分区域相邻布局；同时在时间上体现出差异性，主要是由于研究区特殊的地貌条件使得其优势地类的主要地位未发生变化。