刘凤莲, 杨人懿

(1.云南财经大学 国土资源与持续发展研究所, 昆明 650221; 2.云南财经大学 精准扶贫与发展研究院, 昆明 650221)

土地利用变化是全球变化研究的热点问题[1]，土地利用变化直接影响着全球环境和气候变化[2-3]。土地利用变化影响生态系统的结构、格局及功能[4]，是影响生态系统服务功能变化的重要因素[5-6]。生态系统服务功能是评估区域可持续发展情形的重要视角[6-7]，生态系统服务功能(ecosystem service values,ESV)衰退和丧失也会大大影响土地利用系统的结构和效率，影响人类的可持续发展[8]。研究土地利用变化对生态系统服务价值的影响，对于剖析区域生态环境状况变化，提高地区土地利用水平具有重要的现实意义[9]。近年来遥感技术(Remote Sensing,RS)和地理信息系统(Geographical Information System,GIS)技术的快速发展，RS和GIS技术在土地利用变化研究中的应用也逐渐增多[10-11]，应用RS和GIS技术研究区域土地利用变化，并进一步估算生态系统服务价值的变化，可为区域科学制定环境保护法规和生态补偿方案以及区域的可持续发展提供决策参考[12]。

1997年Costanza等提出了各类生态系统服务价值系数，定量研究了生态系统服务价值，奠定了ESV研究基础[7]。谢高地等展开了进一步研究，估算了我国生态系统服务价值[13-15]。国内其他学者选择了典型区域研究了土地利用结构和格局变化对生态系统服务价值的影响[6,11,16-18]。生态系统服务价值估算方法主要有两种，功能价值法和当量因子法[19]。功能价值法计算复杂，有大量计算参数，且受主观因素影响较大，而当量因子法由于具有操作简单，综合评价，适应性强的优点，被普遍采用[20]。

武汉市地处长江中游，江汉平原东部，是我国中部地区的国家中心城市，无论在自然地理位置，还是在经济地理圈层，武汉市都处于优越的中心位置，区位优势和交通优势显著，有“中国经济地理中心”之称。近年来，武汉市经济飞速发展，加之武汉市是两型社会建设的核心城市，拥有得天独厚的资源条件，是一座极具特色的滨江、滨湖城市。在经济飞速发展的同时，武汉市土地利用格局也发生了巨大改变，也必然会带动生态服务功能的变化。在前人研究的基础上，本研究选取武汉市近25 a来的六期遥感影像，综合利用RS和GIS技术，运用遥感数据解译、空间数据处理、海量数据建库、集成分析等方法[21-22]，对武汉市近25 a来的土地利用变化的时空演变特征进行监测分析，并根据当地实际情况修正各生态系统服务价值当量，研究区域生态系统服务价值变化，为区域国土资源优化配置和生态文明建设提供决策参考。

1 研究区概况

武汉市地处长江中游，江汉平原东部，包括江岸区、江汉区、硚口区、汉阳区、武昌区、洪山区、青山区7个中心城区和东西湖区、蔡甸区、江夏区、黄陂区、新洲区、汉南区6个新城区，地理位置介于113°41′—115°05′E，29°58′—31°22′N，土地总面积8 549.09 km2。武汉市中部平整，北面丘陵林立，南面多岗地，境内水网密集，水资源丰富，水域面积约占全市面积的1/4，拥有优越的自然条件。气候属亚热带季风性湿润气候，雨量充沛、日照充足、四季分明，夏季湿热，冬季凉爽湿润，夏季长达135 d，初夏梅雨季节雨量较集中。2018年末常住人口1 108.1万人，城镇人口889.69万人，城镇化率达到80.29%，地区生产总值为14 847.29亿元，近年来武汉市经济快速发展，土地利用也发生了很大变化。一方面建设用地快速扩张，为经济发展提供了有利支撑，另一方面也出现了一些不容忽视的环境问题，如土地利用综合效益偏低[23-24]，生态环境下降[25-27]，土地资源供需失衡等[28-29]。因此，提高土地集约利用水平，调整土地利用结构，优化土地利用布局，成为武汉市亟待解决的问题。

2 数据来源与研究方法

2.1 数据来源与处理

本研究选取武汉市1990年、1995年、2000年、2005年、2010年、2015年的六期Landsat TM/OLI卫星影像为数据源，影像从美国地质调查局(USGS)网站(https:∥www.usgs.gov/)和地理空间数据云(https:∥www.gscloud.cn/)获取，空间分辨率30 m，影像云量小于3%，便于解译。覆盖武汉市的每期遥感影像均采用3景遥感影像镶嵌而成，3景影像的轨道号分别为123—038，123—039,122—039。借助ENVI5.3，以研究区1∶10 000地形图为基准并结合GPS野外调查控制点对六期遥感影像进行多波段融合、几何校正、影像镶嵌和裁减等预处理工作。本文采用432的波段组合为标准假彩色图像，在这种组合中，所有的植被都显示为红色，在植被、农作物、土地利用和湿地分析的遥感方面，这是最常用的波段组合。

利用ENVI5.3和ArcGIS 10.5软件，结合研究实际需求，首先建立土地利用信息提取标志(表1)，在整个研究区均匀选取训练样本，并对训练样本进行评价，然后选取“支持向量机”(Support Vector Machine Classification，SVM)的分类器进行监督分类，最后对分类结果进行精度评价。各期遥感影像的的总分类精度分别为88.24%，87.56%，89.73%，85.48%，84.18%，88.52%，Kappa系数分别为0.83,0.82,0.86,0.84,0.79,0.85，满足研究精度要求。

2.2 研究方法

2.2.1 土地利用变化分析

(1) 净变化量是表征土地利用变化幅度的指标，用来计算某种地类在研究期内的变化量。其计算公式为：

Ci=Uib-Uia

式中:Ci为研究期内第i种地类的净变化量；Uia和Uib分别为研究期初和研究期末第i种地类的面积。

表1 研究区假彩色影像土地利用信息提取标志

(2) 土地利用动态度是研究土地利用变化的经典范式，其表征的是研究区内特定时间段内不同地类的变化速率。其计算公式为：

式中:Di为第i种地类的土地利用动态度;Ci为研究期内第i种地类的净变化量；Uia为研究期初第i种地类的面积；T为研究时长。综合土地利用动态度为各种地类Di之和。

(3) 土地利用转移矩阵是研究不同土地利用类型转移方向和数量变化的经典方法，可以直观展现各地类格局的演化过程。其数学表达式为：

式中:Aij为研究期初为第i种土地利用类型在研究期末转换成第j种土地利用类型的土地面积；i(i=1,2,…，n),j(j=1,2,…，n)分别表示研究期初和研究期末的地类；n代表土地利用类型数。

2.2.2 生态系统服务价值

(1) 生态系统服务价值评估。参照Costanza[30]和谢高地等[31]提出的生态系统服务价值系数，结合谢高地[13]2015年对中国陆地生态系统服务当量因子的改进，为研究区各土地利用类型各项生态系统服务价值当量赋值(表2)，建设用地的价值当量为0，并通过以下公式计算武汉市生态系统服务价值。

式中:ESV是研究区生态系统服务价值;VCi是第i种地类的单位面积生态系统服务价值；Ai是第i种地类的面积;i为土地利用类型;ESVf是第f项生态系统服务价值;VCfi是第i种地类的单位面积的第f项生态系统服务价值;f为生态系统服务价值的项数;ECf是某种地类第f项生态系统服务的价值当量;Ea为1个标准当量生态系统服务价值。谢高地等[13-14]确定的中国1个陆地生态系统价值当量的经济价值为3 406.50元/hm2,本文选取研究区2015年单位面积粮食产量和当年全国单位面积粮食产量的比值作为修正系数,计算出研究区生态系统服务经济价值当量因子为3 613.36元/hm2。

表2 武汉市单位面积生态服务价值当量

(2) 生态贡献率。生态贡献率是研究周期内不同地类的生态系统服务价值变化量占整个研究区总生态服务价值变化量的比例，能够揭示影响研究区生态服务价值变化的主要贡献因子和敏感因素。其数学表达式为：

式中:Si为第i类土地利用类型在研究期内的生态服务贡献率；ΔESVi是第i类土地利用类型在研究期内的生态服务价值变化量；i为土地利用类型。

3 结果与分析

3.1 武汉市土地利用时空变化

3.1.1 土地利用结构变化 根据遥感解译分类结果(附图1)可以看出，武汉市主要地类为耕地，占土地总面积的一半以上，主要分布在在建设用地的外围，在远城区连绵成片，随着建设用地的扩张，耕地分布趋于破碎；其次为水域，占总面积的1/5左右，以长江和汉水以及汤逊湖、东湖和严西湖等数百个湖泊为主；林地占9%以上，主要分布在黄陂区和新洲区；建设用地占6%以上，主要分布在长江两侧的滨水平原；其他土地利用类型占比较少。

近25 a间，武汉市面积缩减的地类包括耕地、林地和未利用地，其中耕地减少最多，净减少量为811.02 km2，林地和未利用地分别减少了22.12，22.93 km2。建设用地净增加量最大，净增加量为720.08 km2，建设用地比例由6.45%上升到14.87%，其次为水域，增加了138.41 km2。分阶段来看，2005—2010年建设用地、耕地、水域、林地的净变化量均比其他时段大，显示武汉市在这个时段土地利用发生剧烈变化。从具体地类来看，林地1995—2000年增加，其他时段减少，这和武汉市1998年洪灾之后政府全面实施封山育林、退耕还林有关。

研究期内武汉市各地类具有明显的时间变化特征(图1)。具体表现为：耕地面积大幅度缩减，尤其是2005—2010年耕地面积急剧减少，净减少量为347.66 km2，是整个研究周期减少速率最快的时段；建设用地面积不断扩大，且增长速度越来越快，研究期内净增加量为720.08 km2，建设用地的增加量和耕地的减少量基本相当；林地面积总体上呈减少趋势，但是2010—2015年变化趋于平缓；未利用地和草地变化率较大，但面积的绝对变化量较小；水域面积有所增加。

图1 1990-2015年武汉市各类土地利用面积变化

利用ArcGIS的空间分析功能对两个不同年份的土地利用数据叠加分析，得到土地利用转移矩阵(表4)。

表3 1990-2015年武汉市土地利用变化

表4 1990-2015年武汉市土地利用转移矩阵 km2

分析可知，武汉市土地利用的主要转移发生在耕地向建设用地的转移。具体来看，耕地总转化量(1 413.09 km2)的一半(730.71 km2))转化为建设用地，这主要发生在城市边缘区以及原有的农村居民点附近，建设用地扩张大量侵占耕地；其次耕地的重要流向为水域(477.26 km2,33.77%)，这主要发生在湖泊周围，得益于武汉市大力推行退耕还湖政策；耕地向林地转移了154.84 km2，约占耕地总转移量的10.96%，主要发生在山丘与平原的交错地带。水域的总转移量里也有相当大的部分转化为耕地(303.09 km2,61.49%)和建设用地(131.63 km2,26.7%)。耕地和水域的双向转移主要发生在坑塘和湖泊周围。这表明武汉市的城市建设由于山水相间的自然地形格局的约束，城市建设用地扩张主要通过占用耕地和填埋水域来完成。

3.1.2 土地利用动态度 土地利用动态度表征地类变化速度的快慢，代表着某个时期内人类活动对土地利用覆被变化影响的大小。根据武汉市土地利用数据，计算土地利用动态度(表5)。参考全国和武汉市的土地利用变化数据，可以看出如下特点：(1) 从总体情况来看，武汉市1990—2015年土地利用综合动态度为8.84%，为同时期全国水平(1.90%)的4倍多，说明武汉市作为中国中部地区最大都市，社会经济快速发展引起的土地利用覆被变化活跃程度和变化速度远超过全国平均水平。(2) 从各种土地利用类型来看，建设用地的动态度远高于其他用地类型，为5.23%，排在第二位的未利用地为2.04%，耕地的动态度是0.61%，说明近25 a间，武汉市城市化建设活动处于活跃期。(3) 从分阶段来看，武汉市2005—2010年土地利用动态度最高(12.91%)，是研究期内的低值(1990—1995年，5.35%)的2倍多。(4) 分阶段分地类的动态度分析结果显示，研究期内武汉市的城市建设活动处于增速，尤其是2005—2010年，城市建设处于急剧加速时期。从城镇化率来看，2015年末武汉市城镇化率为79.77%，远高于同期全国城镇化率水平56.1%，表明武汉市城市化发展迅速。

表5 1990-2015年武汉市土地利用动态度 %

3.2 武汉市生态系统服务价值变化

3.2.1 生态系统服务价值变化 近25 a来，武汉市生态系统服务价值有所增加，从1990年的1 006.00亿元增加到2015年的1 040.32亿元，增加了3.41%(表6)。水域是武汉市生态系统服务价值构成的主体，占了全市ESV的3/4以上，其次为耕地和林地。总体上来看，水域和草地的ESV有所增加，耕地、林地和未利用地的ESV减少，这和各土地利用类型面积的变化是一致的。从近25 a间武汉市各地类生态系统服务价值的变化可以看出，耕地面积减少引起的耕地生态系统服务价值减少是武汉市生态系统服务价值减少的主要原因，水域面积增加引起的水域生态系统服务价值增加是武汉市生态系统服务价值增加的主要原因，耕地和水域不仅是武汉市土地利用构成的主要部分，也是影响武汉市生态系统服务价值变化的重要原因。

表6 1990-2015年武汉市各地类生态系统服务价值

武汉市各类生态系统服务价值中(表7)，以水文调节和水资源供给为主，其次为气候调节和净化环境等。近25 a间，单项ESV变化最大的是水文调节，增加42.81亿元，增长了6.26%；其次是供给服务中的食物生产，减少3.60亿元，变化率达到了-11.24%；水资源供给减少3.59亿元，变化率为-3.52%；调节服务中的气体调节减少3.03亿元，变化率为-9.20%；净化环境功能增加2.11亿元，变化率为4.93%。水域的ESV系数中水文调节、水资源供给和净化环境等系数较大，水域面积在研究期内增加，直接影响了研究区总体生态系统服务价值的增加；耕地的各项ESV系数中水文调节、水资源供给、食物生产和气体调节等系数较大，耕地面积在研究期内的大量减少，是研究区生态系统服务价值减少的重要原因。由于武汉市对水域的良好保护，使得水域面积不减反增，这对于武汉市维持生态系统服务价值稳定起到了积极的正向作用。

表7 1990-2015年武汉市土地生态系统单项服务价值

水域、耕地和林地是武汉市生态系统服务价值的最重要的组成部分，未来的城市建设工作，应重点加强水域、林地和耕地等生态用地的保护，加快推进区域生态保护红线划定工作，重视城市的合理开发和有序建设，盘活城市的存量用地，避免“摊大饼”式的过度膨胀导致的生态系统失衡甚至崩溃。

3.2.2 各地类对生态系统服务价值变化的生态贡献率 从整体上看，在近25 a的研究期内，水域对生态系统服务价值变化的贡献率最多，其生态贡献率为68.41%；耕地次之，耕地的生态贡献率为29.19%；林地的生态贡献率为1.99%，草地(0.27%)和未利用地(0.13%)较小。分阶段来看，在2000—2005年耕地的生态贡献率最高，达到了59.49%，这个时间段也是林地在整个研究周期生态贡献率最高的阶段，林地的生态贡献率为6.13%。由于水域和耕地在地类转移里双向流动较多，所以水域和耕地的生态系统服务价值变化较大，其生态贡献率变化也较大，最大变化分别是44.49%，47.33%。水域、耕地和林地三者生态贡献率之和在98%以上，水域、耕地和林地的生态系统服务价值变化量对研究区总生态系统服务价值量变化影响较大，因此是主要贡献因子和敏感因子(表8)。

表8 1990-2015年武汉市各地类对ESV变化的生态贡献率 %

4 结 论

(1) 1990—2015年，耕地和水域是研究区主体土地利用类型，土地利用结构发生了显著的变化，建设用地大量增加，增加量为720.08 km2；耕地显著减少，减少量为811.02 km2；水域增加了138.41 km2，林地和未利用地变化不大。土地利用转移主要发生在耕地向建设用地的转移，建设用地扩张主要通过占用耕地和填埋水域来实现。

(2) 1990—2015年，研究区土地利用综合动态度是同时期全国水平的4倍多，土地利用覆被变化活跃程度和变化速度远超全国平均水平。建设用地的动态度为5.23%，远超其他地类，城市建设活跃。尤其是2005—2010年，土地利用动态度最高，土地利用变化加快，2010—2015年土地利用动态度变缓，土地利用趋于稳定。

(3) 1990—2015年，研究区生态系统服务价值呈增加趋势，增加了34.32亿元，3.41%。其中水域、耕地和林地是武汉市生态系统服务价值构成的主体。单项生态服务价值以水文调节和水资源供给为主，其次是气候调节和净化环境，单项ESV变化最大的是水文调节，增加42.81亿元，增长了6.26%。武汉市水域面积在研究期内增加，使得水文调节和水资源供给等单项生态系统服务价值增加，这也是研究区总生态系统服务价值增加的主要原因。

(4) 1990—2015年，水域对ESV变化的贡献率最大，其生态贡献率为68.41%，耕地以29.19%的生态贡献率居第二位，第三是林地。水域、耕地和林地的生态贡献率之和超过了98%，是ESV变化的主要贡献因子和敏感因子。

研究期内武汉市建设用地大量增加，耕地大量减少，土地利用发生变化活跃。将来的土地利用规划中，应重点加强水域、林地和耕地等生态用地的保护，注重土地利用空间格局的优化调控，及时划定生态保护红线，保证生态用地维持在一定数量。区域协调快速城市化发展和生态环境的可持续发展，是未来值得重视的问题。