改革开放以来, 快速城镇化已成为我国经济和社会发展的重要特征。城镇化水平由1978年17.9%提高到2015年的56.1%, 预测至2030年, 中国的城镇化水平将达到70%^[1]。城镇空间扩展是城镇化作用于地理空间的直接表现形式, 其引起的城镇用地增长及土地利用变化是土地利用/覆被变化(Land Use and Land Cover Change, LUCC)的重要研究内容^[2]。然而, 伴随着城镇化进程的持续推进, 城镇用地面积大幅增加, 导致城镇周边大量具有重要生态功能用地(如耕地、林地等)被侵占^[3], 由此带来生境退化、湿地减少、生态系统服务价值损失、生物多样性降低等一系列生态环境安全问题^[4-5], 已经严重影响到人类社会经济的可持续发展。因此, 如何实现城镇化和区域生态安全的“双赢”？一直以来是地理学研究者关注的热点话题。

生境质量是指生态系统提供给个体或种群持续发展生存条件的能力, 被视为区域生物多样性和生态服务水平的重要表征^[6]。城镇化过程形成的土地利用剧烈变化, 通过改变生境结构与组成要素, 使生境的生产能力和服务功能逐渐下降, 进而影响区域气候、植被和生物栖息地等生态环境^[7]。近年来, 国内外学者围绕城镇发展与生态环境之间的关系, 主要开展了以下几个方面研究：(1)城镇空间扩展对某种生态环境要素的影响研究, 包括对气候^[8]、水文^[9]、土壤性质的影响^[10]等。(2)城镇空间扩展对生态环境的响应研究, 如对自然生境的作用^[11]、生态系统服务价值的变化^[12-13]、生态环境质量的评价^[14]等。(3)通过对影响生态敏感性、生态适宜性和生态风险性等因子进行评价^[15-16], 构建生态安全格局, 来倒逼城镇空间增长界限, 优化城镇空间布局^[17]。(4)采用城市演化动态模型对城镇土地利用变化过程进行模拟, 预测城镇空间扩展过程中的生态环境效应^[18-19]。总体来看, 已有研究对城镇空间扩展与生态环境影响的单向评价较多, 而将城镇发展与生态环境保护有机结合起来, 对生态环境质量约束下的城镇自组织增长规律研究较少, 缺乏将生境质量融入到城镇空间扩展中的前置性应用研究。

黑河流域是我国西北干旱区典型的内陆河流域, 地处河西走廊中部, 即是国家“一带一路”战略布局的重要节点, 也是甘肃省祁连山内陆河生态安全屏障的关键区域。黑河中游地区是全流域绿洲分布最为集中、开发条件较好的区域。随着水土资源利用和城镇化进程的快速推进, 不断加剧的人类活动导致流域生态环境质量持续退化^[20]。如何在推进城镇化进程中, 确保绿洲生态系统的稳定性, 是干旱区内陆河流域面临的重大问题。本文以黑河中游地区为研究区, 以该地区绿洲城镇为研究对象。在对研究区生境质量及其退化状况进行评估的基础上, 将其评估结果融入到城镇扩展模型中, 来模拟适宜生境阈值下的城镇规模, 并从城镇扩展占用地类和景观格局指数两方面深入探讨有(无)生境质量约束下的城镇扩展差异性, 进而揭示不同等级城镇空间扩展的基本规律和格局特征, 为绿洲城镇发展与生态环境保护提供科学依据。

1 研究区概况与数据来源 1.1 研究区概况

黑河中游是指从莺落峡至正义峡之间的河段, 长约185 km, 地理位置介于98°57′—100°52′E和38°39′—39°59′N, 南临祁连山, 北依合黎山。行政边界包括甘州区、临泽县和高台县, 总面积为10737.71 km²。黑河贯穿全境, 兰新铁路、G30连霍高速、G312等主要交通干线经过。属温带大陆性干旱气候, 年降雨量为50—250 mm, 年蒸发量为2000—3500 mm。土壤类型以灰棕漠土和灰漠土为主, 植被为小灌木、半灌木等荒漠植被。地势南高北低, 东高西低, 自南向北分为山前洪积扇区、绿洲平原区和戈壁荒漠区。该地区现有36个城镇(图 1), 主要分布在黑河沿岸和绿洲核心区, 且沿黑河干流和交通道路呈串珠状展布, 表现出典型的“亲水性”和“亲路性”特征。1992—2015年城镇建设用地面积从22.33 km²增加至82.53 km², 年均增长率为5.85%。

1.2 数据来源

研究数据主要来源于1992年和2000年Landsat5 TM、2008年Landsat7 ETM+、2015年Landsat8 OLI影像(行带号33/133、33/134)、ASTER GDEM(30 m分辨率), 下载自地理空间数据云网站(http://www.gscloud.cn/)。应用ENVI 4.7和ArcGIS 10.2软件对影像进行几何校正、波段融合、图像增强、拼接等处理。解译时以研究区2015年土地利用变更调查数据为基准, 首先对2015年遥感影像进行人工目视解译, 再参照2015年解译结果对1992年、2000年和2008年遥感影像进行解译, 统一采用WGS_1984坐标系, 最后通过野外抽样验证和谷歌地球高清影像验证, 四期分类数据判读精度高于92%。根据研究需要, 参考土地利用现状分类标准(GB/T21010—2007)及研究区土地利用特点, 建立二级土地利用类型分类体系(表 1)。

2 研究方法 2.1 生境质量评估模型

InVEST模型的生境质量模块(Habitat Quality)常被用于区域生境质量的定量评估^[21-22]。该模型将土地利用/覆被类型与威胁源建立联系, 通过对各威胁因子的敏感性和外界威胁强度来模拟评估生境质量。其建模思路是基于人类活动影响因子进行评价, 通过生境适应度、生境对各威胁因子的敏感性、威胁因子的影响距离和权重、土地受法律保护的程度四个因素来计算生境退化度和生境质量。首先计算生境退化度：

(1)

(2)

(3)

式中, D_xj为土地利用/覆被类型j中x栅格的生境退化度；R为威胁源个数；w_r为威胁源r的权重；Y_r为威胁源的栅格数；r_y为栅格y的胁迫值；i_rxy为栅格y的胁迫值r_y对栅格x的胁迫水平；β_x为威胁源对栅格x的可达性(受法律保护的区域为0, 其余区域为1)；S_jr为土地利用/覆被类型j对威胁源r的敏感度；d_xy为栅格x与栅格y的直线距离；d_rmax为威胁源r的最大胁迫距离。

在此基础上计算生境质量, 其计算公式为：

(4)

式中, Q_xj为土地利用/覆被类型j中x栅格的生境质量指数；H_j为土地利用/覆被类型j的生境适宜度, 值域为[0, 1]；D_xj为土地利用/覆被类型j中x栅格的生境退化度；k为半饱和常数, 取D_xj最大值的一半；z为归一化常量, 通常设为2.5。

本文以2015年土地利用/覆被数据为基础, 假定土地的法律保护处于有效状态, 将耕地、林地、草地和水域等作为生境, 选取人类活动较强的城镇、农村居民点、交通用地和耕地作为生境威胁源。参考模型说明^[23]和研究区实际情况, 对威胁源的最大影响距离及权重、不同类型生境的适宜度及各生境对威胁因子的敏感性进行设置(表 2、表 3)。

2.2 城镇扩展模拟模型

SLEUTH模型能从中观到宏观尺度上模拟城市增长, 且可以考虑政策规划、生态保护等不同城市发展理念, 具有较强的普适性和兼容性, 在国内外一些城市广泛应用^[24-26]。该模型假设历史增长趋势是持续的, 未来的现象可以由过去真实演变数据模拟得到, 通过对历史时期城市范围数据反复训练得到的五个系数(散布系数、繁殖系数、扩展系数、坡度系数、道路引力系数)模拟城市的4种增长方式(自发增长、新扩展中心增长、边缘增长与道路引力增长)。模型的运行分为数据准备阶段、校准阶段、预测阶段。

模型数据准备阶段要求至少一个坡度层(slope)、两个土地利用层(landuse)、一个排除层(exclude)、四个城市范围层(urban)、两个交通道路层(transport)和一个山体阴影层(hillshade)。坡度层和阴影层通过对研究区DEM边界裁剪后, 在ArcGIS 10.2中经空间分析后所得；土地利用层采用1992年和2015土地利用现状；从1992、2000、2008、2015年土地利用中提取城镇建设用地作为城镇范围层；根据遥感影像和张掖市综合交通现状图, 经几何校正后数字化, 提取高速公路、铁路、国省道、县乡道路以及城镇主干道、次干道形成1992、2015年交通道路层；将保护区和河流、水库等禁止城市化区域作为排除层。所有数据均转换为模型需要的GIF格式栅格数据, 栅格大小30 m×30 m, 且保持相同的空间范围和统一的投影坐标系(图 2)。

模型校准阶段是模拟历史年份城市增长过程, 得到一套反映城市演变规律的最优参数。模型采用蒙特卡洛迭代方法, 分粗校准、精校准、终校准和预测参数获取四个阶段进行。在校准的每个阶段, 通过计算Lee-Salee指数或OSM指数对系数组合进行排序, 逐步缩小5个增长系数的取值范围。最后取步长为1, 进行100次蒙特卡罗迭代, 得到研究区最佳增长系数组合, 模型校准的具体参数设置见表 4。本文通过加载土地利用层同时校准UGM模型和LCD模型, 将1992年数据图层作为模型校准的初始图层, 以2000、2008、2015年3个时期的数据图层作为校准图层, 导入模型进行参数校准, 并选用OSM指数^[23]作为参数区间缩小的主要依据, 经过四个阶段校准后, 最终得到5个最佳系数值(表 4)。利用最佳系数初始化模型的预测模块, 反演研究区1992—2015年城镇土地利用变化过程, 运用ROC曲线对2015年模拟图与实际图的一致性进行精度验证。ROC曲线根据其下方面积(Area Under Curve, AUC)来判断各用地类型的拟合度, AUC的取值范围在0.5—1之间, 值大于0.7时, 模型拟合效果较好^[27]。

模型预测阶段通过修改增长系数、调整排除图层来设置不同情景方案, 模拟未来城镇土地利用变化。为探究生境质量对城镇空间扩展的影响, 通过修改排除层分别设置有(无)生境质量约束两种情景(图 2), 来模拟2015—2030年城镇扩展的土地利用变化格局。排除层取值在0—100之间, 值越大, 城市化概率越小。在不考虑城镇扩展中的生境质量约束情景下(情景一), 仅将河流、水库、黑河湿地保护区和东大山自然保护区等禁止建设区域值设为100形成排除层；在考虑城镇扩展中的生境质量约束情景下(情景二), 借助ArcGIS 10.2空间分析的栅格计算器工具, 将InVEST模型输出的生境质量图中各地块单元生境值(0—1)转化为0—100之间, 并叠加黑河湿地保护区和东大山自然保护区(值设为100)形成排除层, 重叠区域取像元最大值。

2.3 景观格局指数

景观格局指数通过高度浓缩景观格局信息, 来定量描述斑块、斑块类别间以及整个景观的空间特征和异质性^[28]。研究以体现景观空间的破碎度、连接性与复杂性为依据来选择景观格局指数, 从景观类型和景观水平两方面对比分析有(无)生境质量约束情景下城镇扩展特征。选择的景观格局指数有斑块数目、斑块密度、平均斑块面积、面积加权平均分维数、平均邻近指数、聚集度。具体计算由Fragstats4.2软件完成。各指数含义如表 5所示。

3 结果分析 3.1 生境质量整体较低, 质量阶序和退化程度分异明显

运行InVEST模型, 输出黑河中游地区2015年生境质量和生境退化度图(图 3), 通过计算, 平均生境质量指数为0.32, 整体生境质量值较低。为便于比较生境质量空间分布状况, 采用ArcGIS 10.2软件的自然断裂法将研究区生境质量划分为高[0.8—1.0]、中[0.4—0.8)、低[0—0.4)3个等级。发现绿洲平原区和山前洪积扇生境质量存在明显的空间分异, 且生境退化程度由绿洲边缘区向核心区递增。高等级生境质量集中分布在黑河干流沿线和植被较好的区域, 以水域和林地、草地为主, 仅占总面积的5.69%, 但出现了一定程度的退化；而在绿洲核心区, 受长期农业活动和城镇村建设等威胁源干扰, 生境出现了较为严重的退化, 尤其是在绿洲腹地的耕作区和城镇周边, 生境质量值为中等, 占总面积的32.77%；低等级生境质量区域主要分布在城镇建设区和绿洲外围的戈壁荒漠区, 占总面积的61.54%, 以裸地、荒草地和沙地为主, 其退化程度低。由此可见, 研究区生境质量和退化程度在空间上不匹配, 城镇发育区域生境质量较好, 但退化程度也高, 城镇发展对生境产生了较强的胁迫性。因此, 为了保证绿洲区域生境质量状况, 必须扼制城镇建设用地的无序蔓延。

3.2 有(无)生境质量约束情景下, 城镇扩展差异很大 3.2.1 SLEUTH模型校准结果与精度验证

由表 4可见, 模型校准最终得到的5个最佳系数值中, 扩展系数(81)和道路引力系数(75)较大, 说明城镇主要以边缘扩展和道路引力扩展为主, 受已有城镇斑块及交通牵引带动作用较强；坡度系数仅为5, 说明城镇扩展较少受地形、坡度因素的制约, 这与多数城镇位于绿洲平原区有关；散步系数控制着城镇的自发增长程度, 其值为55, 表明出现新的城镇斑块的随机性较大；繁殖系数为47, 说明自发增长形成的城镇斑块有可能成为新的城镇增长中心。使用ROC曲线对2015年土地利用模拟图与现状图进行验证, 结果显示, 不同土地利用类型的AUC值中, 林地(0.88)＞未利用地(0.86)＞城镇建设用地(0.83)＞水域(0.80)＞耕地(0.80)＞其他建设用地(0.75)＞草地(0.74), 除其他建设用地和草地的AUC值较低外, 其他用地类型均达到0.8以上, 说明该模型能较好地用于模拟研究区城镇土地利用变化过程。

3.2.2 城镇扩展占用地类差异

以2015年土地利用现状图为基础, 输入各数据图层及最佳增长系数, 调用SLEUTH模型预测模块, 得到2030年城镇扩展的土地利用变化模拟图(图 4), 可见, 有生境质量约束情景下, 城镇扩展规模较小, 空间形态较为紧凑, 有效地抑制了城镇建设用地增长。统计两种情景下城镇扩展占用地类情况(表 6), 发现无生境质量约束情景下, 城镇建设用地面积为181.82 km², 年均增长率为5.41%, 与历史时期基本持平, 其中占用耕地、林地、草地和水域面积分别为47.53 km²、3.62 km²、2.64 km²和1.83 km², 使具有生态功能的用地被较多侵占；有生境质量约束情景下, 城镇建设用地面积为132.58 km², 年均增长率为3.21%, 较历史时期降低了2.64%, 相比较无生境质量约束情景, 占用了较少的生态用地, 分别减少了61.58%、68.89%、95.83%和57.38%。由此可见, 生境质量约束情景下, 能够有效地控制城镇建设用地扩展的数量与速度, 避免了对绿洲生态用地的大量占用, 有利于研究区生态环境保护。

3.2.3 城镇建设用地扩展模式差异

采用景观格局指数对预测年份城镇斑块进行分析, 定量对比有(无)生境质量约束情景下的城镇建设用地扩展特征(表 7)。无生境质量约束情景的斑块数目NP是有生境质量约束情景的1.3倍, 但平均斑块面积AREA_MN仅为有生境质量约束情景的0.6倍, 且边界密度ED是有生境质量约束情景的1.4倍, 表明无生境质量约束情景下的城镇斑块增长速度较快, 具有一定程度的破碎性。比较两种情景下的平均邻近度MPI和聚集度AI, 发现有生境质量约束情景的值都要大于无生境质量约束情景, 说明其相邻城镇斑块之间距离较近, 斑块间的连接性较好, 新增城镇建设用地集聚程度较高, 以“填充式”扩展为主, 有利于保护绿洲生态空间；无生境质量约束情景的FRAC_AM值也较大, 说明城镇斑块形状更加复杂与不规则, 新增城镇建设用地的随机性较强, 以“蔓延式”扩展为主, 对绿洲生境质量产生较强的胁迫性, 是一种不可持续的扩展方式。

3.3 生境质量对研究区城镇扩展的影响

以有生境质量约束情景下的2030年城镇土地利用模拟图为对象, 统计不同等级城镇新增城镇建设用地面积(表 8)。发现预测年份一级城镇新增面积为30.86 km², 较2015年增长了59.47%；二级城镇新增面积为14.66 km², 较2015年增长了77.40%；重点镇新增面积为1.01 km², 较2015年增长了24.4%；一般集镇新增面积为3.52 km², 较2015年增长了46.56%。将不同等级城镇新增城镇建设用地与生境质量分布图进行叠置分析, 发现占用高等级生境质量面积为0.88 km², 仅占新增城镇建设用地的1.76%, 而占用中、低等级生境质量面积为17.46 km²、31.71 km², 分别占新增城镇建设用地的34.88%、63.36%(表 9)。具体来看, 一级城镇主要占用中(33.54%)、低(65.68%)等级生境质量；二级城镇主要占用低等级生境质量(74.56%)；重点镇(71.29%)和一般集镇(84.38%)主要占用中等级生境质量, 但是增长面积有限, 仅占研究区新增城镇建设用地总面积的7.37%。由此可见, 城镇空间扩展总体体现了生境质量的要求, 预测年份新增城镇建设用地以占用中、低等级生境质量为主, 并且低等级生境质量是城镇空间扩展的主要区域。

从新增城镇建设用地空间分布来看(表 8), 绿洲平原区生境质量约束性较强, 城镇空间扩展速度缓慢, 其中：绿洲核心区新增城镇建设用地面积为30.55 km², 较2015年增长了52.01%, 黑河干流沿岸新增城镇建设用地面积为4.89 km², 较2015年仅增长了48.8%；而位于荒漠-绿洲过渡带的绿洲边缘区, 生境质量约束性较低, 且受过境交通牵引力作用, 城镇空间扩展速度较快, 主要体现在城镇产业园区的增长, 新增城镇建设用地面积为14.18 km², 较2015年增长了110.18%；山前洪积扇区城镇距绿洲核心区较远, 且并未处于主要交通沿线, 城镇发展相对孤立, 在局地生境约束下, 新增城镇建设用地面积仅为0.43 km², 较2015年增长了48.31%。

4 讨论

研究采用InVEST模型的生境质量模块, 定量评估了受人类活动影响的生境质量状况, 突破了以往对区域生态空间质量的传统经验主义认识。同时, 在设置城镇空间扩展的情景模拟时将其转化为SLEUTH模型排除层的概率值, 改变了以往人为主观设置的随意性, 也提高了模拟结果的准确性。研究所得的生境质量评估值及不同等级城镇空间扩展特征, 将为研究区生态环境保护政策和城镇发展路径提供重要依据。

干旱区水资源短缺决定了绿洲脆弱的生态环境状况, 为解决绿洲城镇发展过程中的生态环境问题, 可通过整合区域性的生产、生活和生态空间, 实施差异化的空间管制策略, 以寻求可持续的绿洲城镇发展模式。即引导大规模、高强度的生产性活动向生境质量约束力较低的区域转移, 降低高、中等级生境质量区域的城镇开发强度, 使其主要承担生活服务功能和生态保育功能。研究区的绿洲边缘地带, 有生境质量相对较低的戈壁、荒滩等未利用地, 引导城镇产业优先向满足建设条件的该区域布局, 既可以解决城镇发展中的用地问题, 也可以将其对生态环境的负面影响降到最低。而在绿洲平原区应充分发挥城镇集聚效应, 促进区域人口、经济和服务等生活性要素向城镇聚集, 以减轻人类活动对区域生态环境造成的压力。同时, 必须实施集约化的城镇扩展模式, 避免低效率的零散式扩展, 引导城镇空间形态呈团聚状分布, 减少其对绿洲生境的影响。

5 结论

本文以典型的干旱区内陆河流域黑河中游地区为例, 采用InVEST模型对研究区生境质量及其退化状况进行定量评估, 并将评估结果融入到SLEUTH模型中, 对未来城镇空间扩展的土地利用变化格局进行模拟, 分析生境质量对绿洲城镇空间扩展的影响。主要结论如下：

(1) 根据生境质量评估结果, 研究区整体生境质量较低, 高等级生境质量仅占5.69%, 主要分布在黑河干流和植被较好的区域；中等级生境质量占32.77%, 主要分布在绿洲核心区, 退化程度严重；低等级生境质量占61.54%, 主要分布在城镇建设区和绿洲外围的戈壁荒漠区。

(2) 对比两种情景下城镇空间扩展模拟结果, 发现无生境质量约束, 城镇扩展速度较快, 斑块破碎化显著, 以“蔓延式”为主。而在有生境质量约束下, 城镇扩展速度减缓, 城镇建设用地面积和占用生态用地总面积分别减少了27.08%、63.61%, 并且, 城镇空间形态紧凑, 城镇斑块间连接性较好、以“填充式”扩展为主。

(3) 城镇空间扩展总体体现了生境质量的要求, 新增城镇建设用地占用中、低等级生境质量。绿洲平原区和山前洪积扇区生境质量约束性强, 城镇扩展速度缓慢, 新增城镇建设用地面积分别为35.44 km²、0.43 km²；绿洲边缘区生境质量约束性较低, 城镇扩展速度较快, 新增城镇建设用地面积为14.18 km²。

快速城镇化进程中, 社会经济发展往往会凌驾于生态环境保护之上, 城镇空间扩展愈发强烈。后续研究中需要进一步探究生境质量的演化过程及其影响机制, 以不断调整和优化生境结构与组成要素, 科学判断生态环境所能承载的最大城镇规模, 划定城镇空间增长边界, 有效地控制城镇扩展规模。此外, 本文是基于“生态优先”理念, 将生境质量结果融入到城镇扩展模拟模型中, 来研究未来城镇空间自组织增长规律。但是, 城镇空间扩展受社会、经济、政策等多因素和多微观主体博弈的综合作用过程, 存在许多不确定性。因此, 综合考虑影响城镇发展各类因子, 建立融合多因子的地理模拟系统, 模拟未来城镇发展趋势, 将是未来研究的重要议题。