重庆市幅员面积8.24×10⁴ km²，辖38个区(县)、707个镇街、252个乡. 2014年常住人口3016.55万人，人口城镇化率60.94%.重庆市城市布局集中，农村地区地域广大，发展不均衡性明显(图 1).以重庆市的主体功能区、区(县)、乡(镇)3级单元为研究样本.采用重庆市区(县)2009年、2014年土地利用变更调查数据，利用ARCGIS软件在1:10 000土地利用现状数据库中，提取各区(县)、乡(镇)行政边界、2009年及2014年城镇建设用地、农村居民点用地面积.行政区界线有所调整的区域，以2014年行政边界为基础进行修正.为凸显土地城镇化变化的政策意义，按照重庆市主体功能区域划分(图 2)，兼顾行政界限的完整性，将重庆市划分为主城区、渝西地区、渝东南和渝东北四大功能区.

1) 土地城镇化水平测度

土地城镇化是土地利用方式、景观形态向城镇的转变，也是土地承载的经济社会活动从乡到城的转换.对土地城镇化的测度主要体现为对土地利用结构改变的速度和规模的量化.现有研究采用区域城镇建设用地的占比或者增速^[19]、城乡建设用地中城镇工矿用地占比^[20]等单一指标，或景观变化、投入产出强度等综合指标^[21]衡量土地城镇化水平.参考现有方法，城镇建设用地在城镇建设用地及农村居民点用地总和中的占比这一算法将城乡建设用地统筹考虑，体现了城镇化进程中土地利用结构的改变^[22]，也与人口城镇化的计算方法相对应.因此，本文选用该方法计算土地城镇化水平.公式如下：

式中：G_it为i区域在t年度的土地城镇化水平；U_it，R_it分别为i单元在t年度的城镇建设用地、农村居民点用地规模.

2) 土地城镇化变化

采用土地城镇化变化指数计算土地城镇化变化幅度.公式如下：

式中：S_i为i区域土地镇化变化指数；G_it1，G_it2分别为i区域在t₁，t₂年份土地城镇化水平.

1) 全局空间自相关

通过计算Global Moran's I反映土地城镇化空间集聚总体特征.公式如下：

式中：I为全局Moran指数；n为参与分析的镇域单元数；X_i和X_j分别为空间单元i和j的土地城镇化水平；$\bar X $为全市各乡镇土地城镇化率的均值；W_ij为镇域i和j之间的空间邻接矩阵. Global Moran's I的正负分别表示计算单元与周边单元的土地城镇化水平正相关有统计学意义或有较大的差异.

2) 局部空间自相关

通过计算局部Local Moran's I_i反映土地城镇化与其周边单元的相关程度，以此分析相邻镇域土地城镇化水平是否在空间上存在高值集聚和低值集聚的区域.公式如下：

式中：Z_i，Z_j为镇i和j的土地城镇化的标准化值，其余参数同式(3).按照I和Z_j的取值，将土地城镇化发展程度划分为HH，LL，HL，LH 4个象限.

尺度方差及其分解是尺度效应分析测度的主要方法之一，用以说明差别化的尺度变化过程对研究整体的变异性的贡献水平^[23].模型如下^[24]：

根据以上模型，得出主体功能区、区(县)、乡(镇)3级方差分解：

式6、式7中：SSα，SSβ分别为区县级、乡镇级方差分解；I是功能区个数；J_i是第i个功能区所包含的区(县)单元数；K_ij是第i个区(县)所包含的乡(镇)单元数.

各层级平均方差：

2009-2014年，全市土地城镇化水平由27.67%上升至32.18%(图 3)，受大量在城市常住而未退出农村土地的“半市民化”人口“两栖占地”的影响，全市土地城镇化水平总体较人口城镇化水平低.主体功能区尺度上表现为功能区之间差距稳定，功能区内部增长平缓的格局态势. 2009年以来，四大主体功能区土地城镇化水平从大到小总体表现为主城区、渝西地区、渝东北地区、渝东南地区的格局特性.

区(县)土地城镇化水平最低的仅有7.98%，最高的达到100%(图 3)，呈以主城区为中心向外围区域逐步递减的“圈层—外围”结构.并表现为“S”型特性变化规律，即处于城镇化水平初期阶段和后期阶段的地区发展相对平缓，在中期阶段的地区发展加剧.高土地城镇化区域集中分布在主城区，且土地城镇化水平趋于稳定，是高水平低变化幅度区域.以主城区为核心，向周边地区扩展，土地城镇化水平逐次降低.土地城镇化水平30%以上的区域集中分布在渝西地区，并呈现快速增长的趋势.在重庆市经济发展水平较低、人口城镇化率较低的渝东北、渝东南地区，土地城镇化水平大多在20%以下，并表现低速增长的趋势.万州、黔江等区域中心城市，在两翼地区具有较大的经济集聚能力，土地城镇化水平相对较高.

乡(镇)级尺度的土地城镇化水平，与区(县)尺度表现出的“圈层-外围”特征基本一致，但高土地城镇化水平圈层的集聚特征更加显著，向外围递减趋势呈空间断层，并围绕主要的交通干线呈跳跃性带状、点状分布特点. 2009年和2014年两个时间节点，土地城镇化水平较高的乡镇集中分布在主城区(图 4)，沿主城区向外扩展，在渝西地区形成较高土地城镇化乡镇的梯级分布空间.其余地区高土地城镇化水平的乡镇分布较少.在主城区内部，高土地城镇化水平的乡镇以长江、嘉陵江交汇点(渝中半岛)为核心，呈更加显著的集聚圈层.在主城区外围，高土地城镇化水平的乡镇集中趋势减缓，分布分散.沿长江和主要交通干线，形成带状、点状分布的高土地城镇化水平乡镇.在土地城镇化水平整体不高的区县，围绕中心城区，形成高土地城镇化水平的点状分布乡镇.

由于功能区级的样本数量少，且土地城镇化率的变化相对稳定，因此着重对比区(县)、乡(镇)土地城镇化空间格局特征(图 5).

1) 垂直方向，3级尺度单元所表现的土地城镇化空间格局“圈层-外围”特性一致. 2009年、2014年全市区(县)尺度土地城镇化空间自相关全局指数为0.623 4，0.660 6，在5%水平上有统计学意义.表明区(县)级尺度土地城镇化水平空间集聚分布特征显著且总体稳定. 2009年、2014年全市乡(镇)尺度土地城镇化空间自相关全局指数为0.526 4，0.469 7，表明乡(镇)级尺度上土地城镇化分布呈现高土地城镇化水平值相邻，低土地城镇化水平值相邻的空间集聚特征.区(县)尺度空间自相关全局指数大于乡(镇)，表明乡(镇)尺度总体集聚的态势稍有减弱，但总体集聚趋势一致.

2) 水平方向，同一尺度单元内，小尺度单元土地城镇化集聚类型变化显著，大尺度单元土地城镇化集聚类型相对稳定.区(县)尺度，高高集聚(HH)的区域主要为主城区.与2009年相比，2014年区(县)尺度高高集聚区域保持不变，低低集聚(LL)的区域也基本稳定.乡(镇)尺度高高集中的区域主要在主城区，在外围区(县)呈少量分布.低低集聚的区域主要分布在两翼地区，以及在渝西地区的边缘地带.与2009年相比，2014年乡(镇)尺度低低值区域出现了明显的增加或者减少区域.

尺度方差分解结果中，区县级最大(2009年为0.064 7，2014年为0.072 4)、功能区级次之(2009年为0.048 1，2014年为0.054 8)、乡级最小(2009年为0.038 4，2014年为0.042 4)，表明功能区、区(县)级、乡(镇)级尺度方差均呈扩大趋势.对土地城镇化的差异贡献来说，区县尺度土地城镇化差异贡献最大，而乡镇级的差异贡献最小.尺度方差分解表明，区(县)土地城镇化格局是形成全市土地城镇化空间分异的核心尺度，区(县)土地城镇化空间集聚特性奠定了全市土地城镇化的基础格局.在此基础上，乡(镇)尺度土地城镇化体现出变化幅度大、变化类型多的特点，丰富了土地城镇化的空间格局形态.

不同于单体城市城镇空间扩展的同向性机理，区域土地城镇化的演变是在区域经济发展中的空间分工新格局，同时受到自上而下的制度变迁和自下而上的城镇化双向驱动^[25].土地作为经济社会发展的载体，土地城镇化自然受政策环境、经济增长、产业结构、人口集聚等因素的综合影响^[26]，参考已有研究，并考虑数据的可获取性，选取其中受到较多关注的人口城镇化率、人均地区生产总值、全社会固定投资总额3个因子，分别分析其对土地城镇化空间格局的影响.结合以上3个因素的空间分布，可以明显发现土地城镇化水平的高值区和人口城镇化水平高值区、人均地区生产总值高值区、全社会固定资产投资总额的高值区有一定的空间耦合关系.相关关系分析表明，整体而言，区县尺度土地城镇化水平与以上3个因素相关关系均有统计学意义，相关系数分别为0.94，0.85，0.49.经直辖以来的发展，重庆市的经济、产业、人口等经济社会水平明显形成了明显的功能区集聚的格局.这种社会经济集聚发展的格局也同样映射到土地上，形成与区域的经济发展总体格局基本一致土地城镇化格局.表明注重设施建设同步发展等新型城镇化的内涵特性已经在重庆有所体现^[27].区(县)级尺度上土地城镇化水平也形成了以主城区为高值区，外围区域水平相对较低的格局结构.

受自然条件和经济基础的双重约束，山区、库区和少数民族地区离主城区和区域中心城市较远，地形复杂，投资带动力度小，城市扩展力度有限，县域土地城镇化率普遍较低.受城市集聚区域的辐射带动，主城区、六大区域性中心城市的县域土地城镇化水平相对较高.在主体功能区战略驱动下，渝西地区受主城区退二进三产业支撑，基础设施投资力度加大，人口集聚能力加强，土地城镇化水平快速增长.

在同一时段，区(县)级行政区域的制度政策、市场化程度等城市化政策、经济驱动力具有较大的同质性，县域范围内不同乡(镇)土地城镇化受地形条件、区位特征等自然要素的影响更大.地貌等自然条件是导致土地城镇化宏观分布格局的基本因素.一般说来，随着海拔的上升，地形起伏度的加大，开发成本增大，土地城镇化增速会随着海拔的提高而降低.分析发现，重庆市镇域土地城镇化水平与海拔并不存在明显的线性拟合关系.但随着海拔的升高，乡镇的土地城镇化水平逐渐降低.在平均海拔超过800 m的488个乡镇中，仅有6.9%的乡镇土地城镇化水平达到40%以上.而在海拔低于800 m的乡镇中，土地城镇化水平超过40%的乡镇个数达到海拔800 m以下乡镇总数的20%.将各个乡镇的平均海拔与土地城镇化水平绘制散点图，可以明显看出高土地城镇化水平的区域在800 m海拔以下的区域集中(图 6).绝大部分土地城镇化水平大于90%的乡镇位于海拔500 m以下.高海拔区域的较高土地城镇化水平的乡镇大多数位于风景名胜区，具有优良的资源开发优势，对当地的城镇化有显著带动作用.长江是重要的交通轴线，也是重庆市主要的经济发展轴带.长江沿线分布的99个乡镇中，土地城镇化水平高于40%的乡镇达到43个，其中全市土地城镇化水平达到100%的乡镇中的50%分布于长江沿线.区(县)中心城区是各个区县的经济、文化中心，各区(县)中心城区所在乡镇土地城镇化均远远超过区(县)土地城镇化的平均水平.

1) 土地城镇化空间格局演变具有明显的尺度依赖性. 2009-2014年，重庆市主体功能区、区(县)、乡(镇)3个尺度单元呈现出垂直方向格局特性从上至下逐层制约；水平方向上呈现尺度越小，异质性越显著，集聚趋势变化越剧烈的特性.区(县)尺度土地城镇化水平不均衡明显，空间分布呈现“圈层-外围”的格局特性，并具有明显的空间集聚特征.乡(镇)尺度土地城镇化水平不均衡程度加剧，在“圈层-外围”总体格局的基础上，表现出明显的点状、带状分布特性.

2) 土地城镇化空间格局的成因也具有明显的尺度差异.区(县)尺度土地城镇化水平受到经济发展、产业集聚、人口吸纳等多因素的综合影响，土地城镇化的空间分异与区域经济发展格局一致.在中低海拔集中分布、沿长江沿线带状分布、各区(县)中心城区点状分布的高土地城镇化乡(镇)格局说明，乡(镇)尺度土地城镇化水平受地形、区位等因素影响更大.

3) 尺度方差分解表明，尽管乡(镇)级小尺度单元的土地城镇化格局变化更加剧烈，但区(县)行政单元对城镇化差异程度的贡献度最大，乡(镇)级最小.表明我国的城镇化战略以县级行政单元实施效率最高，但乡镇级行政单元的城镇化发展能有效优化局部地区的城镇化格局.

1) 对乡(镇)尺度土地城镇化的格局特征和演变规律还有待深入探讨.乡镇是社会经济活动的基础行政单元，土地城镇化格局和演变还应受到乡(镇)产业发展、人口集聚等多因素的综合影响.但限于数据收集难度大，特别是限于对乡镇土地利用的长时序数据收集困难，乡(镇)行政单元的土地城镇化格局演变机理还有待深入研究.

2) 对乡(镇)尺度土地城镇化和人口、经济城镇化的耦合关系还有待深入研究.土地是社会经济发展的载体，对于土地城镇化的研究，出发点是解决我国土地城镇化冒进、土地城镇化和人口城镇化失调等问题.因此，相关研究在分析土地城镇化其特征和形成机理的基础上，更注重深入探讨其与人口城镇化的协调度、剖析土地城镇化和人口城镇化失调原因.之后的研究中，有必要在乡镇尺度上，深入探寻土地城镇化与人口城镇化的耦合度，为寻求土地、人口、经济协调发展的城镇化提供依据.