2001—2020 年成渝地区城市扩展及建成区植被变化

薛飞阳, 王鑫, 青玲萱, 苏子昕, 杨存建*

1. 西南土地资源评价与监测重点实验室, 四川 成都 610068；

2. 四川师范大学地理与资源科学学院, 四川 成都 610068

城市化是现代化进程的重要结构性变迁，对经济社会的发展具有多重作用，是以农业为主的传统乡村型社会向非农产业为主的现代城市社会逐渐转变的历史过程，是人口向城市地区聚集及农村地区或者自然区域转变为城市地区的过程[1]。近年来我国城镇化发展迅速，截至2021 年末我国常住人口城镇化率为64.72%。随着经济飞速发展，城市化进程不断加速，城市的发展也受到越来越多的重视[2]。

城市化的飞速发展，大量人口涌入城市，促使一些城市对外进行扩张，而城市建成区的不断扩大，必然影响到建成区内外植被的生长、分布情况[3]。赵安周等[4]分析了2000—2015 年京津冀13 个城市生长季及不同季节城市扩张对植被和地表城市热岛影响的时间变化趋势以及二者之间的关系，结果表明，京津冀13 个城市主城区生长季及不同季节的城市扩张对植被均存在消极影响。董晨炜等[5]基于DMSP/OLS 夜间灯光数据和MODIS NDVI 数据，研究环杭州湾地区城市扩张和建成区植被变化，结果表明2000—2013 年，城市扩张对建成区植被均产生了不利影响。王静等[6]定量研究了2000—2010 年，京津冀地区植被覆盖度及其景观格局的动态变化，揭示了城市化进程对植被景观的干扰过程及生态质量的影响。

夜间灯光遥感数据可以直观反映夜间人类活动差异，具有覆盖范围大、时效快和易获取等优势，可以广泛应用于多尺度长时序的城市问题研究[7]。同时，传统遥感数据在城市社会经济特征的监测中没有有效的直接手段，存在估算准确率不高的问题[8]。成渝地区位于我国西南地区，城市化水平远低于东部沿海城市，对成渝地区的城市发展可以构建科学合理的城市体系，同时对促进整个西南地区城市质量的提升和城市持续发展提供现实意义。因此，本文将采用夜间灯光数据对成渝地区城市扩展及建成区植被变化进行分析。

1 数据与方法

1.1 研究区概况

成渝地区位于我国四川盆地腹地，是我国西部最发达的地区，也是我国重要的城市群之一，成渝地区包括四川省的成都、绵阳、德阳等十五个地级市以及重庆市。该地区总面积为20 万km2，人口8 724 万，分别占中国西南地区的6.6%和35%。成渝地区以亚热带季风气候为主，地形多山地丘陵，河流众多，气候适宜。由于毗邻区域合作不断深化，川渝合作进程逐步加快，一体化发展的趋势日益明显。成渝地区在中国经济区中人口最多，也是西部人口密集度最高的区域，同时也拥有成渝两座国家中心城市和数座大中型城市（见图1）。

图1 研究区地理位置Fig. 1 Geographical location of the study area

1.2 数据源与处理

数据主要为2001、2006、2010 年的DMSP/OLS夜间灯光数据和2015、2020 年的NPP/VIRS 夜间灯光数据（美国国家海洋和大气管理局https://www.noaa.gov/）及其城市建成区数据（来自《中国城市统计年鉴》）；SPOT/VEGETATION NDVI 年度合成产品（分辨率1 km，用最大值合成法由月度数据生成），全国1∶400 万县级行政界限（均来自中国科学院地理科学与资源研究所https://www.resdc.cn/）

1.3 研究方法

1.3.1 城市扩展分析

分别选取城市扩展速度和扩展动态度两个指标来量化城市扩张的变化程度，城市扩张速度是在一定时期内城市用地面积变化的快慢，城市扩张动态度指数表示单位时间内城市用地面积变化的幅度[9]，城市扩展速度与扩展动态度计算公式如下：

式中，V 和K 表示城市扩展速度和扩展动态度，Ai和An+i分别为第i 年和第n+i 年的城市建成区斑块面积，n 为时间跨度。

紧凑度可以反映城市外围轮廓在空间上的离散程度[10]，以紧凑度的数值变化来分析城市扩展的空间变化特征[11]。形态紧凑度计算公式：

式中，C 表示紧凑度， Ai为第i 年城市建成区斑块面积，P 为第i 年城市建成区斑块外围周长。值介于0-1 之间，值越大，表示城市外部的形状越接近圆形，值越小，表示城市外部形状越狭长。

城市重心转移的变化趋势可以表示城市在演变中的空间变化特征[12]。重心的经纬度坐标计算公式如下：

式中： wi为第i 个离散目标物权重，由各个离散目标的面积与总面积的比例来确定， xi和 yi分别是第i 个离散目标的经、纬度坐标。

式中：Dn+1为重心偏移距离，（xn、yn）和（xn+1、yn+1）分别是第n、n+1 年建成区区域的重心坐标。

1.3.2 NDVI 植被分析

基于提取到的城市建成区的边界，采用NDVI 均值对城市的植被演变特征进行表示。

式中，CNDVI为NDVI 均值，pi为每个像元的NDVI值，n 为像元个数。由于得出的NDVI 值介于0-1 之间，为了便于计算，在不影响精度的情况下，对得出的值乘以10 000 进行分析。

对建成区与背景区的NDVI 均值进行比值分析，来衡量在城市演变过程中人为因素对植被变化的影响。表达式为：

式中，Ri为第i 年NDVI 比值，CNDVi为第i 年建成区NDVI 均值，BNDVi为第i 年背景区NDVI 均值。

趋势分析和变化百分率2 个方面来分析城市扩展过程中植被变化趋势[13]，表达式为：

式中，slope 为趋势线的斜率，NDVIi为第i 年的NDVI 值，n 为5，时间增加，NDVI 为上升趋势，反之，则为下降趋势。NDVIC为变化百分率，Mean为5 年的平均NDVI。

变异系数可以描述数据的离散程度，可以用于评估NDVI 在时间序列上的稳定性[14]。表达式为：

式中，CV 为变异系数，ST 和Mean 为NDVI 标准差和5 年的平均NDVI 值。若CV≤0.1,表示NDVI 值很稳定，0.1＜CV≤0.2 为稳定，0.2＜CV≤0.3 为不稳定，＞0.3 则为很不稳定。

2 结果与分析

2.1 城市建成区提取

采用统计数据比较法[15]，根据阈值分别提取各时期相应城市建成区范围。首先，分别设定各阈值城市市区用地信息，再将提取到的城市建成区面积与统计年鉴进行比较，修改阈值至二者差距最小为止，最后基于最佳阈值（见图2）提取出16 个城市的城市建成区空间范围。对比提取结果与统计结果，提取误差均小于2%（见表1），提取精度较高，可用于后续分析。

表1 基于参考比较法的成渝地区用地面积Tab. 1 Land use area in Chengdu-Chongqing area based on reference comparison method

图2 2001—2020 年各城市不同时期建成区提取阈值Fig. 2 Extraction threshold of built-up areas in different periods of cities from 2001 to 2020

2.2 城市建成区时空动态特征

2.2.1 时空分布特征

2001、2006、2010、2015 和2020 年成渝地区各城市的夜间灯光值呈现出明显的空间分布差异（见图3）。根据提取结果，2001—2020 年间研究区各地级市城市建成区面积处于不断增长中，城市规模比较集中分布在成都和重庆，导致大城市的规模比较大，而中小城市规模较小。整体上，2001—2020 年中小城市的数量和规模持续增长，成都、重庆城市规模仍保持较快速度增长。

图3 2001—2020 年城市建成区提取结果叠加图Fig. 3 Stacked map of extraction results of urban built-up areas from 2001 to 2020

2.2.2 速度变化特征及动态度变化特征

由表2 可以看出，2001—2006 年重庆扩展速度最快，成都次之，乐山为负增长，其他城市之间差距较小；2006 年—2010 年，重庆快速扩展，成都增速趋缓，雅安最慢；2010—2015 年，重庆最快，成都开始快速发展，资阳最慢；2015—2020 年，成都最快，重庆速度减缓，乐山最慢。从年际变化来看，成都城市扩展速度在2006—2010 年（14.66 km2·a-1）有所减少，但在2010—2020 年（32.03 km2·a-1）开始快速发展，到2015—2020 年（72.28 km2·a-1）的速度最快；重庆在2010—2020 年间城市扩展速度趋缓，由2006—2010 年的81.46 km2·a-1减少到2015—2020 年的43.66 km2·a-1；南充和宜宾的城市扩展速度在2001—2020 年持续增加。总的来看，2001—2020 年，成渝地区城市扩展速度整体稳步上升。

表2 2001—2020 年成渝地区扩展速度及扩展动态度变化Tab. 2 Changes in expansion speed and dynamic degree of expansion in Chengdu-Chongqing area from 2001 to 2020

从城市扩展动态度来看（见表2），2001—2006年，除德阳、乐山、广安和遂宁的城市扩展速度停滞或萎缩外，其余城市都有不同程度的扩展，其中，重庆和成都都在快速扩张，南充城市扩展动态度最大（16%）而广安最小（-8.14%）；2006—2010年，重庆与成都呈中高速扩展，其余城市低速扩展，其中，达州扩展最快；2010—2015 年，在扩展速度上，成都快速扩张，而重庆扩展速度减缓，在扩展动态度，重庆和成都间差距并不明显，其中，内江扩展动态度最大；2015—2020 年，成都城市扩展速度和扩展动态度超过重庆，其中，宜宾和重庆扩展动态度最大。

2.2.3 紧凑度

除重庆外，研究区其他城市的紧凑度指数较高（见图4）。由于重庆市城区位于四川盆地东部丘陵地区，受南北山地的地形所限，建成区多沿山间盆地呈狭长形分布，紧凑度低。2001—2020 年，乐山和资阳的紧凑度指数呈增加趋势，表示其城市发展由相对无序向相对有序发展；成都、绵阳、南充、宜宾和重庆的紧凑度指数持续递减，其中成都的紧凑度指数变化最大；而乐山、内江、遂宁和资阳紧凑度指数相对稳定，表明城市发展有序。

图4 成渝地区不同时期紧凑度特征Fig. 4 Compactness characteristics of different periods in Chengdu-Chongqing area

2.2.4 重心转移

2001 年成渝地区的建成区重心在资阳市（见图5），2001—2006 年重心向东北方向移动了15.05 km，但仍在资阳市内，2006—2010 年重心向东南方向移动了13.62 km，由资阳转移到了重庆，这一时期，重庆的建成区面积超过其他地区，2010—2015 年重心向东南方向迁移了1.10 km，仍在重庆市内，2015—2020 年重心向西北方向移动了11.63 km，又由重庆转移到资阳。

图5 成渝地区城市建成区重心变化Fig. 5 Changes of gravity center of urban built-up areas in Chengdu-Chongqing area

2.3 植被变化

2.3.1 植被时间变化特征

从建成区面积总和与NDVI 总和来看（见图6）。2001—2006 年建成区面积增加了586.03 km2，建成区NDVI 总和减少了1 690.64；2006—2010 年建成区面积与建成区NDVI 总和分别增加了616.58 km2、8 125.12；2010—2015 与2015—2020 年间建成区面积增加了898.51 km2、978.74 km2，建成区NDVI 总和减少了6 809.74、6 688.66。总的来看，随着城市的扩展，城市建成区面积总和持续增加，在2010-2020 年间城市发展速度呈骤增趋势，但城市建成区NDVI 总和呈先上升后下降的趋势。

图6 成渝地区建成区面积总和与NDVI 总和Fig. 6 The sum of built-up area and NDVI in Chengdu-Chongqing area

2.3.2 NDVI 均值

从建成区年均NDVI 变化情况来看（见图7），2001—2006 近一半城市NDVI 均值在上升，其中成都和重庆植被情况有所改善，其中绵阳植被增加量最大；2006—2010 年有近1/3 的城市植被均值持续下降，植被变差，其中达州植被退化最为明显；2010—2015 年近1/3 的城市植被均值增加，成都NDVI 均值增加量最大，而重庆植被退化明显；2015—2020年仅有1/4 的城市植被均值增加，成都和重庆NDVI均值有不同程度的增加，表明城市建成区植被情况有所改善。结合表3 可以看出，2001—2020 年近3/4的城市建成区植被均值持续降低。

表3 2001—2020 年成渝地区城市建成区内NDVI 比值及变化量Tab. 3 NDVI ratio and its variation in urban built-up areas in Chengdu-Chongqing area from 2001 to 2020

表4 2001—2020 年成渝地区NDVI 年际变化率和变化百分率Tab. 4 Annual change rate and percentage of NDVI in Chengdu-Chongqing area from 2001 to 2020

表5 成渝地区城市建成区NDVI 变异系数Tab. 5 Variation coefficient of NDVI in urban built-up areas in Chengdu-Chongqing area

图7 成渝地区城市建成区植被变化情况Fig. 7 Vegetation changes in urban built-up areas in Chengdu-Chongqing area

2.3.3 比值及变化量

2001—2020 年，各城市建成区NDVI 比值均小于1，说明城市扩张对建成区植被均产生不利影响。2001—2006 年，成都、乐山、泸州、绵阳、内江和重庆市建成区NDVI 比值变大，建成区植被有所改善，其他城市扩张对建成区的不利影响增大；2006—2010 年，眉山、绵阳和遂宁市建成区NDVI 比值减小，对城市建成区植被不利影响增大，而其他城市的扩张对植被影响减小；2010—2015 年，成都、德阳、眉山和自贡市建成区NDVI 比值增加，建成区植被改善，而其他城市扩张都对建成区植被不利影响增大；2015—2020 年，南充、宜宾、重庆和资阳市的城市扩张对建成区植被的不利影响减小，其他城市的扩张对建成区的植被不利影响增大。总体来看，2001—2020 年，除成都、重庆、绵阳、内江和泸州市的城市扩张对建成区植被不利影响下降，其他近2/3 的城市扩张对建成区植被状况的不利影响明显增加。

2.3.4 趋势分析和变化百分率

从NDVI 年际变化率来看，城市建成区有近2/3的城市自然植被呈下降趋势，其中眉山和遂宁下降显著，建成区内NDVI 变化百分率处于负增长，说明成渝地区城市建成区内植被趋于退化；而背景区有近2/3 的城市自然植被上升趋势，说明该区域植被处于良好发展。总体来看，城市建成区植被处于退化趋势。

2.3.5 植被稳定性

从城市建成区NDVI 变异系数来看，其中乐山变异系数最大，是植被稳定性最差的城市。2001—2006 年，有近1/3 的城市变异系数呈增加趋势；2006—2010 年，近1/2 的城市变异系数增加，呈上升趋势；2010—2015 年，近4/5 的城市变异系数呈增加趋势；2015—2020 年，近2/3 的城市变异系数增加。总体来看，在2001—2020 年变异系数均有不同程度增加，说明该时段内研究区植被比较不稳定。

3 讨论

城市的扩张是城市发展的必然结果，受到社会经济因素和自然因素的共同影响，尤其是政策、经济与人口的发展[16]。成渝地区城市在2001—2010 年间城市发展较为缓慢，但在2010 年后保持可观的速度扩展，社会经济活动形势向好。2016 年《成渝城市群发展规划》的发布，成渝城市群开始快速发展，目前城市群正处于经济社会发展上升期，但是生态保护与经济发展的平衡尚未建立[17]。欧阳晓等[18]在中国城市群城市用地扩张时空动态特征研究中发现东部地区城市在2005 年以前扩张强度达到极大值，2005 年之后，扩张强度呈现下降趋势。与本文研究结果恰好相反，说明我国城市群呈现不平衡的发展趋势，东部城市发展明显高于中西部城市。

随着建成区的扩展，2001—2020 年间城市扩张对城市建成区植被产生不利影响。研究区位于四川盆地，植被覆盖的变化不仅与气候因子密切相关，而且也受到人为因素的影响[19]。宋梦来等[20]以天津市为例，得出植被变化受到气候和人类活动因素的共同影响，人类活动因素对植被覆盖度的影响远远超过气候因素[21]。研究发现，成渝地区城市建成区植被生长状况较自然背景区植被较差，与余晨渝等[22]研究得出人类活动对植被覆盖负向影响区域主要集中于人类活动密集、城市化进程较快的城区基本一致。随着建成区的扩展，即使城市内有着完善的植被保护措施，但还是会不可避免的占用或破坏天然林草地，导致植被覆盖度的下降。

研究还存在一些不足，有待进一步改进。一方面，虽然将DMSP/OLS 和NPP/VIRS 两种夜间灯光进行拟合，但仍存在误差，而且由于DMSP/OLS 的分辨率为1km，空间分辨率过大，不利于进一步研究城市建成区内植被的变化过程。另一方面，比较依赖于统计数据来提取城市建成区的边界，统计年鉴中数据的误差或缺失会对研究结果的准确性产生影响，未来可引入珞珈一号等高分辨率的夜间灯光数据或土地资料数据提取建成区。

4 结论

基于2001—2020 年整合夜间灯光数据，运用统计数据比较法提取成渝地区城市建成区边界，对城市扩展特征进行分析，并引入SPOT/VEGETATION NDVI 数据，以NDVI 总和、比值及变化量和稳定性等对城市建成区植被进行分析，得出以下结论。

（1）成渝地区城市扩展年均速呈上升趋势。2001—2006、2006—2010、2010—2015、2015—2020年城市扩展速度为7.33 km2·a-1、9.63 km2·a-1、11.23 km2·a-1、12.23 km2·a-1。

（2）成渝地区除重庆外，研究区其他城市的紧凑度指数较高。由于重庆市城区受南北山地的地形所限，建成区多沿山间盆地呈狭长形分布，紧凑度较低。

（3）成渝地区城市重心在2001—2006 年重心向东北方向移动，2006—2010 年重心向东南方向移动，2010—2015 年重心向东南方向迁移，2015—2020年重心向西北方向迁移。

（4）成渝地区城市建成区植被在时空分布上总体呈退化趋势。NDVI 总和在2001—2020 年间呈先上升后下降趋势，在2006—2010 年间城市建成区植被NDVI 总和随着城市扩展不断上升，但在2010—2020年间城市建成区植被随着城市建成区扩展不断退化。

（5）成渝地区城市建成区植被普遍较背景区差。2001—2020 年间半数以上的城市建成区植被的NDVI 比值及变化量呈下降趋势，而背景区植被NDVI 均值整体超过建成区植被，且半数以上城市呈增加趋势。

（6）成渝地区城市植被稳定性变差。2001—2020 年超过半数城市的变异系数均有不同程度的增加，表示该时段城市植被比较不稳定。