地球迈入“人类世”后, 认知人类活动与地表环境系统的影响机理进而调控和管理人类行为, 成为国际“Future Earth”研究计划、IPCC项目、生物多样性和生态系统服务政府间科学-政策平台(IPBES)等重大研究计划倡导的前沿领域与优先级主题^[1-2]。自20世纪80年代, 国内地理学者也先后强调要将自然环境过程与人文社会经济现象的互动机理研究作为当代地理学发展的重要方向之一^[3-4]。近年来, 生态环境与城镇化两大系统的耦合与胁迫机理调控的研究框架强调自然与人文两大系统之间的动态作用, 也要认知系统内部各要素之间对另一系统要素的动态效应^[5-6]。全球范围的城镇化引发的城镇土地扩张, 是人类活动影响自然环境系统改变的直观体现^[7]。城镇建设用地, 是陆地表层系统演化过程中人地矛盾最突出、最不可逆转的地域, 也是城镇化复杂系统内部自组织和动态演化效应对生态环境系统影响最深、作用最直接的要素之一^[8-9]。由城镇建设用地扩展引发的区域生物地球化学循环^[10]、气温变化^[11]、生物多样性减少^[12]、不透水表面增加^[13]等相继引起科学家们的重视。因此, 在对城镇体系规划、资源环境管理和调控、市场分析和服务分配等领域宏观决策前, 探明历史阶段和现状实情的城镇建设用地格局与扩展过程, 实证解析城镇建设用地扩展变化的主控要素及影响机理, 是区域可持续发展的一项必要性工作^[14]。

对已有城镇建设用地扩展研究作以归纳和分析, 在用地数据源方面, 美国国家地理信息中心(NGDC)提供的夜间灯光影像适宜全球大、中尺度城镇信息提取^[15], Landsat等多源遥感数据与统计数据^[16-17]多用于区域中、小尺度城镇信息提取, 我国搭建了LUCC时空平台提供上世纪末至今的LUCC数据集^[18]。城镇建设用地时空特征模拟方面, 除传统土地利用动态度等表征指标外, 近年来相继应用自组织映射神经网络^[19]、Dyna-CLUE^[20]、FLUS^[21]、元胞自动机CA(Cellular Automata)与多智能体ABM(Agent-Based Model)^[22]等新模型。在识别驱动因子的方法和技术手段上, 相关或回归分析^[23]、SD系统动力学^[24]、协整分析^[25]等及基于空间分层异质性原理开发的地理探测器(Geodetector)新方法^[26-28]先后被有效应用。已有文献结论^[27-29]证实自然要素中的地形、海拔、到市中心距离等和人文要素中的GDP、固定资产投资、人口密度、二三产业产值、城镇化率等是影响城镇建设用地扩展的主要因素。综合来看, 遥感监测数据与新模型方法的结合、自然条件和社会经济条件综合考虑探寻演变机理, 是国内外城镇建设用地扩展研究的新近趋势。

2014年5月二次中央新疆工作会议, 新疆被确立为“丝绸之路经济带核心区”, 2015年4月, 丝绸之路经济带被全国两会列为中国新时期发展的三条战略支撑带之一, 2016年3月, 丝绸之路经济带核心区建设作为“十三五”规划国家完善对外开放布局的重要部署。城镇建设用地作为丝绸之路经济带核心区生产与生活的重要场所, 也是丝绸之路经济带国内段边疆支点与走廊门户建设水准的重要体现, 把控好丝绸之路经济带核心区新疆今天与未来城镇建设用地布局的现实意义深远。然而, 丝绸之路经济带核心区新疆建设用地时空格局及其影响因子的目前研究^[30-33]尚存在一些不足之处：①在城镇建设用地时空格局的研究单元上, 主要以个别区域或个别城市为案例区, 全疆尺度的完整分析较少；②研究时段上, 以近10—20年时间序列居多；③驱动要素探索上, 多局限于分析人文经济影响因素, 而鲜少见自然环境本底因素和人文经济因素共同解释用地演变的时空分异程度及各因素交互作用影响的着眼点。鉴于此, 本研究通过提取1980—2015年丝绸之路经济带核心区新疆城镇建设用地变化信息, 从总量趋势、分地州市级区域、分用地类型及特殊性全面探析扩展的时空演变特征, 在县域尺度与各影响要素进行空间耦合匹配分析, 定量判别影响城镇建设用地扩展的人文主控要素和自然主控要素, 探寻其作用机理, 并以核心区重点发展区域天山北坡城市群和喀什都市圈的区域分异特征作分析和对比, 以期对丝绸之路经济带核心区新疆城镇建设用地有序调控和可持续管理提供科学认知及决策依据。

1 数据与研究方法 1.1 研究区概况

丝绸之路经济带核心区新疆, 位于欧亚大陆中部, 横跨东经73°40′—96°23′, 纵踞北纬34°25′—49°10′。疆域面积166 km²位列中国省区最广, 边境线5600 km位列中国省区最长, 接壤邻境8国位列中国省区最多, 是“古丝绸之路”的重要经商通道、民族迁徙走廊和多元文化交融地, 是“新丝绸之路”北、中、南三条大通道的必经门户。疆域典型地貌构造单元由高纬至低纬依次是：阿尔泰山、准噶尔盆地、天山、塔里木盆地、昆仑山, 地形起伏大, 海拔高程-155 m至8611 m, 地物分布的垂直地带性明显, 温带大陆性干旱气候区, 年均降水约155 mm, 现有绿洲面积约14.3万km²。2015年土地覆被类型按分布面积占比值大小依次为：未利用地(60.8%)、草地(28.6%)、耕地(4.7%)、水域(3.2%)、林地(2.3%)、城乡建设用地(0.4%), 见图 1。20世纪80年代至2015年, 丝绸之路经济带核心区新疆全域城镇化率由29.0%提升至47.2%, 人口密度7.7%增至14.2%, 年末总人口达2360万人, 约占全国总人口的1.7%；GDP由53亿元增至9325亿元, 约占全国总GDP的1.4%；固定资产投资额由20.5亿元增至2606亿元, 约占全国总投资额的1.8%。

本文研究区域定义为丝绸之路经济带核心区新疆下辖的14地、州、市, 约占全疆域总面积的99.5%, 不包括生产建设兵团“师市合一”城市性质的8个县级市。另天山北坡城市群是国家主体功能区划的重点开发区, 也是中国—中亚—西亚国际经济走廊的国内段起点, 喀什都市圈是新疆“十三五”规划纲要的重点建设区, 也是中巴经济走廊的国内段起点, 因此将核心区这两个重点区域作为案例区, 对比分析了二者用地扩展的区域特征、主控要素及其影响机理的共性与分异性。

1.2 数据来源与影响要素选择

研究区空间行政边界矢量数据, 取自全国地理信息资源目录服务系统(http://www.webmap.cn)1:100万基础地理信息数据库。土地利用栅格数据, 采用中国科学院资源环境数据云共享平台(http://www.resdc.cn)1980—2015年7期1 km遥感监测数据, 各期数据重建的用地类型综合判别精度均符合制图要求, 解译分类系统一致利于长时间序列用地信息获取。DEM数据来源于国家地理信息中心。各年份有关社会经济数据, 来源于《中国县(市)社会经济统计年鉴》、《新疆五十年》、《新疆统计年鉴》等文献资料。

城镇建设用地扩展是社会经济发展和环境本底条件约束共同作用结果。丝绸之路经济带核心区新疆三大山系包围两大盆地的特殊地貌轮廓是城市空间扩展的刚性环境约束基底, 海拔高程、坡度及地形位等自然要素的地域分异特征明显。自改革开放后, 研究区经济体量和城市体量不断壮大, 产业结构优化明显, 现代化建设成效显著, 非均衡性城市发展差异仍客观存在。基于已有城镇建设用地空间扩展因素研究的理论与实证分析结果^[26-33], 考虑近35年研究区人文因子的数据可得性与自然因子的可定量性质和对城镇建设用地时空演变的影响程度, 本文影响要素指标的选择如下：选取X1(城镇化率指数)、X2(人口密度指数)、X3(GDP指数)、X4(二三产业结构指数)、X5(固定资产投资指数), 作为表征社会经济发展的人文驱动要素；选取X6(地形起伏度)、X7(坡度)和X8(地形位), 作为表征区域本底环境的自然驱动要素。所选指标均采用县域统计数据, 市单元不细分至辖区。

1.3 研究方法

(1) 城镇建设用地时空演变信息提取与统计

文中定义的城镇建设用地包含城镇用地与其它建设用地2类, 现丝绸之路经济带核心区新疆城镇建设用地总面积约为6700 km², 约占全疆域总面积0.4%, 面积数值变化难以揭示时空扩展形态演进, 遥感监测产品和GIS技术结合能直观刻画用地演变过程。对1980—2015年7期土地利用栅格数据, 经重分类将城镇建设用地像元类型定义为1, 非城镇建设用地定义为0。针对城镇建设用地的动态演变特征, 经栅格运算和制图融合识别出稳定型、扩展型和转化为其它用地三种类型^[34], 提取各年份城镇建设用地时空信息。根据丝绸之路经济带核心区新疆地域面积大, 城镇建设用地零星分散特点, 经比对不同尺寸的网格单元大小对县域尺度用地信息的空间表达效果, 最终确定, 以创建10 km×10 km渔网(Fishnet)为宜, 在研究区域形成15704个完整的百公里方格单元, 与各年份用地数据作叠置分析, 记录每个网格内城镇建设用地面积及其占网格面积百分比值；运用属性选择功能和字段计算器工具, 对叠置后的数据图层属性表结果融合处理, 统计丝绸之路经济带核心区新疆城镇建设用地时空演变信息。用城市扩展强度指数UEII(Urban Expansion Intensity Index)表征起止时相的年均用地规模数量在时间维度的相对变化幅度, 表达式为^[35]：

(1)

式中, UA_n+i与UA_i分别表示起止时相的城镇建设用地面积, n为起止时相间隔年份。

(2) 地形要素指数提取

地形起伏度、坡度和地形位指数是地形因子最常用的指标^[36]。其中, 地形起伏度, 是特定地域内最高海拔与最低海拔高程的差值, 宏观表征地表形态和地貌类型的空间分异程度。地形位指数, 是高程和坡度的复合地理因子, 综合表征城镇建设用地空间格局与地形梯度的关系, 表达式为^[37]：

(2)

式中, T为地形位指数；E和S分别为某点的海拔高程值(m)和坡度(°), E₀和S₀分别为某点所在区域的平均海拔高程值(m)和平均坡度(°), 本文定义该点为各行政单元的政府所在地。海拔高程越高、坡度越大的区域地形位指数愈大, 反之则愈小。对研究区DEM数据利用表面分析工具提取研究区地形起伏度、坡度和地形位指数, 然后对各地形要素进行等级划分。

(3) 地理探测器模型

基于空间分层异质原理开发的地理探测器, 是定量表达不同区域之间地理现象差异性因子解释效应的有力工具, 且可揭示任意两个自变量对因变量的交互作用类型, 已被广泛运用在自然环境科学和社会经济等领域^[38-42]。该模型可将用地扩展的多类影响要素在同一空间尺度下离散归一化处理, 能有效解析城市扩展空间分异的驱动机理^[26-28]。其优点在于影响因素的地理探测力q值没有线性假设, 其明确的物理含义能客观地探测出自变量解释100×q%的因变量^[43-44], 地理探测力q值的表达公式为：

(3)

式中, q_{D, U}为城镇建设用地变化影响因素的解释力指标, n为整个区域样本数, m为次级区域个数, n_{D, i}为次一级区域样本数, σ²为全域建设用地变化量的方差, σ_{UD, i}²为次一级区域的方差。假设σ_{UD, i}²≠0, 模型成立, q_{D, U}取值区间在[0, 1], q_{D, U}=0时, 表示城镇建设用地不受影响因素驱动而呈随机分布, q_{D, U}越趋近于1, 说明驱动因素对城镇建设用地的解释效力越大。本文应用地理探测器的因子(Factor)探测和交互(Interaction)探测两个子模块分析人文要素和自然要素对丝绸之路经济带核心区新疆城镇建设用地扩展的影响程度及其交互作用。

2 丝绸之路经济带核心区新疆城镇建设用地扩展的时空演变特征

基于各期遥感解译数据集提取的用地类型信息, 从总量趋势、分地州市级区域、分二级用地类型及特殊性方面把握丝绸之路经济带核心区新疆城镇建设用地扩展的时空格局与演变过程。研究区全域城镇建设用地总面积由1980年的881 km²扩展至2015年的3416 km², 35年间共扩展约2.9倍, 时间序列下的空间演变特征如下图 2。

2.1 扩展总量不断增加, 但强度波动增长且以分散小斑块为主

据时间序列的动态视角, 丝绸之路经济带核心区新疆城镇建设用地面积总量不断增加, 但扩展强度呈波动增长态势。具体表现为：1980—1990年和2000—2005年为较快扩展, 年均扩展强度为4.9%、5.4%；1990—1995年为近乎停滞状态, 年均仅扩展0.2%；1995—2000年和2005—2010年为缓慢扩展, 年均扩展强度1.2%、1.6%；2010年后城镇建设用地扩展为近35年来最强, 年均扩展强度15.4%。据空间过程的动态视角, 整体呈现扩展态势, 但城镇建设用地斑块以小面积单元为主, 且在核心区全域空间极为分散。如表 1, 丝绸之路经济带核心区新疆城镇建设用地斑块网格数量在1990—1995年扩展停滞阶段略减少, 其他年份均增加；百公里方格网内分布的城镇建设用地斑块面积以1—10 km²为主导, 占总网格数57%—67%；百公里方格网内分布的城镇建设用地在10—50 km²斑块数量逐年增加, 50 km²以上斑块数量有扩增潜力。说明城镇建设用地总体仍呈现扩展态势, 用地斑块以小面积单元为主, 且在核心区全域空间极为分散。

2.2 分地州市城镇建设用地扩展的时空差异性显著

丝绸之路经济带核心区新疆各地州市城镇建设用地扩展的时空差异性显著。在面积规模上, 核心区扩展面积最大的地州市乌鲁木齐市城镇建设用地面积值476 km², 城镇建设用地超过10 km²方格网数量占全市的1/2, 超过50 km²方格网数为4个, 占了核心区全域的1/2；博尔塔拉蒙古自治州、克孜勒苏柯尔克孜自治州、和田地区三地州的城镇建设用地面积都不足60 km²。在扩展程度上, 有“华夏第一州”之称的巴音郭楞蒙古自治州为最高, 昌吉回族自治州和乌鲁木齐市次之；南疆克孜勒苏柯尔克孜自治州、和田地区1980—2015年仅扩展了8 km²、19 km², 是全疆城镇建设用地扩展最慢的地州。在时序变化上, 如图 3, 1990—1995年和1995—2000年, 少数地州市单元出现不同幅度城镇建设用地向其它用地类型转化情况；2005—2010年, 除巴音郭楞蒙古自治州、喀什地区增势较为明显外, 其他地州市均无明显扩展态势；进入“十二五”后, 各地州市城市建设步伐加快, 除克拉玛依市、塔城地区外, 其余12地州市的城镇建设用地扩展强度均为近35年来最高。

2.3 次一级城镇建设用地类型扩展的动态变化明显

城镇建设用地分次一级类型来看, 丝绸之路经济带核心区新疆城镇用地与其它建设用地两类扩展过程有明显动态变化特征。如下图 4所示, 1980年城镇用地和其它建设用地面积占比约为2:1, 城镇用地分布为主导；1990—2000年城镇用地面积占比逐年增加, 2000年占比69.8%达到最大值；2005—2015年其它建设用地扩展明显增强, 其面积占比逐年增加至2015年超过城镇用地, 成为次一级用地的主导类型。20世纪改革开放初期, 丝绸之路经济带核心区新疆其它建设用地主要零星分布在首府乌鲁木齐市、克拉玛市、吐鲁番市、哈密市、昌吉回族自治州等少数资源型城市和工业化较为发达的地州市。大型工业区、矿业区、采石场、油田等工业场地开发建设及交通道路铺设与基础设施新建活动, 促进其它建设用地大面积扩展。尤其自“十二五”以来, 其它建设用地在在巴音郭楞蒙古自治州、阿勒泰地区、阿克苏地区、伊犁州直属县(市)、塔城地区扩展强度也均增幅明显。

2.4 城镇建设用地扩展受自然本底胁迫和经济社会要素影响的特殊性

丝绸之路经济带核心区新疆受自然本底胁迫条件和经济社会发展进程的多重因素, 城镇建设用地扩展有其特殊性。全域高山与盆地相排列, 盆地被高山所环抱, 海拔高程、坡度及地形位等自然要素的地域分异特征明显, 决定了适宜城镇建设用地比重较低的特殊性。孤立且分散的绿洲多分布在塔里木盆地和准噶尔盆地的南缘、西缘、北缘的近山地区, 是人类生产生活空间的主要承载地域。研究区2015年城镇化水平低于全国同期近9个百分点, 二三产业比重不足60%的县域占全域的15%, 平均每10000km²土地面积仅分布14人, 社会基础设施建设和运营投资成本较高, 决定了城镇建设用地布局在空间零星散布的不均衡性和扩展缓慢的特殊性。城镇建设用地扩展区域主要在首府乌鲁木齐市、工业化水平较高的克拉玛依市、条带状县域单元密集的昌吉回族自治州及吐鲁番市和哈密市；相较之, 自然环境约束程度更高和经济社会发展水平更低的和田地区、克孜勒苏柯尔克孜自治州等县域城镇建设用地面积较少、扩展较慢。

3 丝绸之路经济带核心区新疆城镇建设用地扩展的影响要素与机理分析 3.1 城镇建设用地扩展程度与各影响要素的空间耦合分析

利用城镇建设用地扩展程度与各影响要素的分级在县域尺度上进行空间耦合匹配分析, 其中乌鲁木齐市和克拉玛依市2地级市, 因缺少早期市辖区统计资料按现各市级行政区为统计单元；哈密市和吐鲁番市2地级市, 在2015—2016年更设地级市前分别为哈密地区和吐鲁番地区各下辖1市2县, 现以各县域行政区为统计单元。

与各类要素进行空间匹配结果如图 5所示, 城镇建设用地扩展较低的县域单元集中在和田地区、喀什地区、克孜勒苏柯尔克孜自治州、阿克苏地区、博尔塔拉蒙古自治州、伊犁州直属县(市)、阿勒泰地区、巴音郭楞蒙古自治州大部分及哈密地区北部。中城镇建设用地扩展县域主要集中在昌吉回族自治州东部、塔城地区部分县域、吐鲁番市东西两翼及巴音郭楞蒙古自治州库尔勒市。城镇建设用地扩展程度较高的地区集中在乌鲁木齐市和克拉玛依市2个传统的地州级城市, 及巴音郭楞蒙古自治州东部的若羌县和哈密市伊州区(2016年撤地设市前为原哈密地区哈密市)。

3.2 丝绸之路经济带核心区新疆城镇建设用地扩展的主控影响要素

根据地理探测器模型因子分析模块, 8项要素对城镇建设用地扩展的影响力q_{D, U}值的探测结果(p value = 0.00)分别为X1城镇化率(0.41)、X2人口密度(0.20)、X3 GDP(0.33)、X4二三产业比重(0.37)、X5固定资产投资(0.31)、X6地形起伏度(0.40)、X7坡度(0.37)、X8地形位(0.57)。由此可知, 地形位、地形起伏度、城镇化发展水平、二三产业比重、坡度的要素影响效力较大；GDP要素、固定资产投资和人口密度要素影响较小。

交互作用探测结果显示各要素对丝绸之路经济带核心区新疆城镇建设用地扩展的影响起到正向强化效应, 各类要素两两交互作用后的综合因子解释力显著大于各因子独立作用之和(表 2)。在28对交互因子中, 共计有17对交互因子综合作用的解释力超过0.990, GDP与二三产业比重、地形起伏度、地形位指数, 人口密度与地形位指数、城镇化率与固定资产投资5对交互因子的影响力q_{D, U}=1, 说明该种组合综合影响力与城镇建设用地扩展变化在空间上有显著一致性, 极强地解释了县域建设用地扩展的分异程度。各交互要素组合效果显示为非线性协同与双要素协同两类, 非线性协同效果要更显著。其中, 城镇化率∩地形位指数的交互作用为双要素协同类型, 其他所有要素交互组合为非线性协同类型。

综合地理探测器因子分析模块和交互作用模块的探测结果, 判定地形位指数、城镇化率、地形起伏度、二三产业比重为影响丝绸之路经济带核心区新疆城镇建设用地扩展的主控要素。

3.3 丝绸之路经济带核心区新疆城镇建设用地扩展的主控要素作用机理

进一步分析研究区影响城镇建设用地扩展主控要素的作用机理, 为城镇因地制宜、科学规划和调控未来发展提供参考依据。

(1) 高地形位指数限制城镇建设用地场所的选择。该指数是地形与坡度的综合表征, 地形位特征为用地格局提供基础, 地形梯度也影响着用地结构与变化方向。该项单因子作用的解释力为0.57, 是影响丝绸之路经济带核心区新疆县域城镇建设用地分布与扩展规模的首要主控要素。受地形位指数限制较大的县域主要分布在塔里木盆地北缘、天山南麓的阿克苏地区、喀什地区、巴音郭楞蒙古自治州北部, 及阿尔泰山一带阿勒泰地区。地形位指数低的区域通常交通较为便利, 经济水平较为发达, 信息化条件也较好, 人类居住生活和从事生产活动的适宜度与便利度更高, 是城镇建设用地新建的倾向性选择；而高地形梯度区限制城镇建设用地场所的选择与扩建, 且高地形位区用地规模扩展对所需的水资源供给、交通道路等配套设施及建筑安全系数皆增加了难度。

(2) 城镇化率提升刺激城镇建设用地扩展。从县域尺度看, 中高城镇化水平增长刺激城镇建设用地扩展集中在天山北坡城市群县域及东疆若羌县和伊州区。近35年以来, 丝绸之路经济带核心区新疆建制市数量明显增长, 城市体系以小城镇为主向大中城市规模转变, 城市行政地域规模及其常住人口总量直接决定其城镇化水平的高低, 进而影响其实体地域范围扩展。1980年仅有乌鲁木齐1个地级市和克拉玛依等7个县级市, 2015年末发展至乌鲁木齐、克拉玛依、吐鲁番3个地级市和22个县级市。同期, 全域城镇化水平提升18.2%, 城镇建设用地面积扩展2535 km², 即城镇化水平每提升1%, 所新增的城镇建设用地面积为139 km²。自治区“撤县建市”、“撤乡建镇”、小城镇试点改革、新型城镇化规划及行动计划等工作推进及全疆农业转移人口落户城镇的措施与机制制度的补全健全, 农村转移和迁移人口落户城镇后的居住、从业、教育、娱乐等功能区的土地需求量对城镇建设用地扩展起正向协同促进效应。

(3) 高地形起伏度胁迫城镇建设用地扩展。平坦开阔之地, 用地选址的自由度较大, 且规模扩展的可操作性较强。核心区全域以中海拔梯度为主要地形地貌特征, 其占比超过1%的类型有：低海拔(< 1 km)占22.4%、中海拔(1—2 km)占42.3%、中高海拔(2—4 km)占14.1%、高海拔(4—6 km)占15.78%、极高海拔(>6 km)占2.8%及现代冰川占1.4%。高地形起伏地势对城镇建设用地扩的开发和扩展具有相当大程度的胁迫性, 连片的低城镇建设用地扩展高、中地形起伏度匹配区域集中在昆仑山和天山山脉西侧的边境城镇带及北部阿尔泰山麓的部分县域, 高地形起伏度区城镇建设适宜用地占比也较低。

(4) 二三产业比重增高拉动城镇建设用地扩展。区域产业结构水平, 是体现生产社会化程度与市场经济水平高低的重要标志。全域三次产业结构比由1980年40.4:40.3:19.3变动为2015年16.7:38.6:44.7, 实现自“一二三”以农业为主向“三二一”以服务业为主的产业结构升级；同期就业结构也由70.0:14.8:15.2变动为44.1:15.2:40.8。在县域尺度单元, 富蕴县(60.55%)、柯坪县(60.2%)、阿克陶县(57.0%)的二三产业比重提升幅度均实现50%的跨越式发展。一方面, 工业、建筑业、各类服务业等二三产业投入产出效率相对较高, 生产要素集聚和规模发展均需以土地作为载体, 并对与之配套的交通铺设、商业及公共服务空间等设施提出更高的需求, 势必致使城镇建设用地扩展的刚性需求；同时, “西部大开发”战略及内地省市援疆计划引进建设项目和资金, 大规模工业园区、商业园区、基础设施及民生工程设施等新建在很大程度上拉动城镇建设用地定向扩展。

3.4 丝绸之路经济带核心区“一群一圈”城镇建设用地扩展分异的影响要素及机理

天山北坡城市群, 区域总面积19.28万km², 城镇建设用地总面积由1980年的472 km²扩展至2015年的1598 km², 35年间共扩展约2.4倍。喀什都市圈, 区域总面积18.46万km², 城镇建设用地总面积则由1980年的66 km²扩展至2015年的253 km², 35年间共扩展约2.8倍。

根据地理探测器模型分析结果, 各要素对“一群一圈”城镇建设用地扩展的影响力如下图 6。经对比可知, 固定资产投资、坡度、地形位指数均是两重点发展区域共同的主控影响要素(q_{D, U}均高于0.95), GDP指数的影响程度都较低；有明显分异性的是, 天山北坡城市群二三产业比重要素影响力明显(q_{D, U}=0.86), 而喀什都市圈人口密度要素影响力明显(q_{D, U}=0.94), 说明产业结构升级与人口集聚分别是二者城镇建设用地扩展分异的显著驱动因素, 这与其各自发展阶段有关。天山北坡城市群历来是全疆人口和产业最密集的区域, 也是核心区内唯一的城市群建设单元, 产业结构升级的利益驱动主导着城镇建设用地扩展, 即处于城市群发展的较高级阶段；而喀什都市圈尚处于区域发展的初级阶段, 以人口数量的增加导致城镇建设用地的自然扩展。此外, 城镇化率驱动作用天山北坡城市群强于喀什都市圈, 而喀什都市圈的地形起伏度限制程度则要高于天山北坡城市群。

4 结论与讨论

本文将丝绸之路经济带核心区新疆近35年以来的城镇建设用地扩展作为一个时空变化系统进行分析, 提取1980—2015年七个年份城镇建设用地演变信息, 借助10 km×10 km网格单元对时间过程和空间扩展程度加以统计, 从总量趋势、分地州市级区域、分用地类型及特殊性全面探析扩展的时空演变特征, 采用地理探测器模型, 在县域尺度定量诊断城镇建设用地扩展的人文要素和自然要素的影响程度及要素之间的交互作用, 进而识别主控要素并解释其作用机理, 探析天山北坡城市群和喀什都市圈两个重点发展区域的分异性。主要结论如下：

(1) 丝绸之路经济带核心区新疆城镇建设用地总面积由1980年的881 km²扩展至2015年的3416 km², 35年间共扩展约2.9倍。在时间序列上, 面积总量不断增加, 但扩展并非逐年增强而呈现波动增长态势, 2010—2015年全域扩展强度最大；在空间过程上, 总体呈现扩展态势, 斑块以小面积单元为主且极为分散。各地州市城镇建设用地的面积规模和扩展程度时空差异性显著。城镇用地与其它建设用地两类次一级用地扩展具有动态变化特征。受自然本底胁迫条件和经济社会发展过程的多重影响, 丝绸之路经济带核心区新疆城镇建设用地扩展有其特殊性。

(2) 对城镇建设用地扩展程度与各影响要素在县域尺度进行空间耦合匹配分析, 全域以城镇建设用地低水平扩展与各级影响要素水平的匹配关系为主, 中高程度城镇建设用地扩展集中在北疆天山北坡城市群大部分县市及东疆部分县市和巴音郭楞蒙古自治州东部的若羌县。该结论与前人在全疆域地州尺度的研究结果基本吻合, 且在速率程度、结构变化及空间分异的大体态势也趋近^[31]。

(3) 综合地理探测器因子识别和交互作用探测结果, 判定地形位指数、城镇化率、地形起伏度、二三产业比重为影响丝绸之路经济带核心区新疆城镇建设用地扩展的主控要素。综合来看, 以自然环境要素为主导, 经济发展要素次之, 人口集聚要素对城镇建设用地扩展的影响最小。进一步分析各主控要素的作用机理, 为丝绸之路经济带核心区新疆城镇因地制宜、科学规划和调控未来发展提供参考依据。

(4) 重点发展区域天山北坡城市群和喀什都市圈城镇建设用地扩展的各类要素影响力有共性也有分异性。两区域共同的主控影响要素为固定资产投资、坡度和地形位指数, 而二三产业比重要素在天山北坡城市群影响力明显, 人口密度要素在喀什都市圈影响力明显；且城镇化率驱动和地形起伏度胁迫程度在两个区域内有分异性。该结论也验证了已有研究发现因区域发展阶段不同而使城市扩展的主控驱动要素有别的客观差异, 通常经济发达区域的城市扩展以经济增长为主要驱动力, 而欠发达区域则以人口涌入等自然扩展驱动力为主^[45], 较为发达的天山北坡城市群与尚处在初级阶段的喀什都市圈的对比结果恰验证了这一客观规律。

1980年以来, 处于西部省区的丝绸之路经济带核心区新疆与我国城市扩展集中发生地的东中部地区相比较, 其城镇建设用地时空演变特征共性体现在扩展总量不断增加, 且主要集中在城市群地区, 得益于改革开发政策后中国城镇化与城市群建设的辉煌成就^[46]。有明显差异的是, 2010年后, 国家主体功能区规划对国土空间开发格局的有效管控抑制了“摊大饼”传统扩展模式, 东、中部地区建设用地扩展面积分别呈现大幅减少、略有减少态势, 扩展速度放缓^[18]；而该时段丝绸之路经济带核心区新疆扩展幅度为近35年来最强。东中部地区城市扩展的县域差异日趋缩小, 中小城市发展加快, 两级分化现象缓解, 且新增城镇建设用地的次一级类型始终以城镇用地为主^[47-48]；而新疆城市扩展以建设用地基数大的城市为主, 尤其南疆及部分边境中小城市扩展仍相对缓慢, 次一级类型独立工矿等其它建设用地自2005年后扩展明显增强, 至今面积占比量已超过城镇用地。东中部湿润、半湿润地区城镇建设用地扩展主要在地形起伏度较小区域^{[37, 47]}, 而干旱区新疆自然本底基础决定了城镇发育以分散绿洲为主体依托, 绿洲城镇沿三大山脉麓前盘踞, 环南北疆两大盆地边缘呈串珠状展布。究其根本, 丝绸之路经济带核心区新疆与东中部地区的城镇建设用地扩展差异由地域分异自然基础与社会经济发展历史进程和现状阶段差异所致。

科学地监测和认知是可持续发展系统解决方案分解步骤的第一步, 本文以分析1980—2015年长时间序列下丝绸之路经济带核心区新疆城镇建设用地扩展时空演变特征——识别城镇建设用地扩展与各影响要素的耦合匹配关系——判定影响城镇建设用地扩展的主控要素——探究其作用机理为框架思路展开研究。将陆地表层的自然要素与人类社会经济要素结合的综合研究, 是当代地理学科发展的重要方向。地理探测器模型的引入, 不仅定量诊断核心区县域尺度的不同类人文经济要素和自然本底要素对城镇建设用地扩展的影响及交互作用进而识别其主控驱动要素, 而且对比分析了核心区两个重点区域天山北坡城市群和喀什都市圈用地扩展的共性与差异, 进一步揭示用地扩展分异主控要素及其影响机理, 兼顾了城镇建设用地扩展的区域异质性与时空非平稳性的考虑。但研究仍存在一定的局限性：1)本文仅以县域尺度时间序列完整、可定量的经济社会因素指标和自然因素指标作以实证剖析, 国家及自治区层面的土地差别化宏观政策及具体指令性用地配额等行政因素^[49]未计入探测范围, 一方面因其在县域尺度难以定量化, 另一方面已采用的固定资产投资额指标、二三产业结构优化是政策实施结果的重要表现方面。2)局限于对近35年以来现有数据的探究, 后续分析可运用系统耦合模型和计量经济学模型对耦合关系与状态、预测及预警方面加以深入研究。在丝绸之路经济带核心区新疆经济社会稳定发展的前提下, 充分考虑自然环境本底胁迫性、经济社会的发展程度及区域特殊性, 把握分区域分城市的主控要素, 通过一定的政策手段引导城镇建设用地的理性扩展和总体布局, 提高用地效率和协同发展水平, 因地因城因类差别化地制定区域调控措施, 进一步加快丝绸之路经济带国内段边疆支点和走廊门户建设。