昌 亭，周生路，戴 靓，陈 龙

（南京大学 地理与海洋科学学院，南京210023）

近年来，我国在工业化、城镇化快速发展的背景下，土地利用方式更加多样化，土地利用更加强调集约节约利用，区域土地生态问题凸显并逐渐成为研究的热点［1－2］。目前国内外关于土地生态系统做了较多的研究，主要集中在土壤污染［3］、土地利用景观格局［4］、土地生态评价（承载力评价［5］、安全评价［6－7］、风险评价［8］、适宜性评价［9］）、土地生态预警［10］等方面。城市化是驱动土地生态系统变化的重要因素，为了探究城镇化是否会影响土地生态的质量，本文提出区域土地生态质量分布具有城乡梯度规律这一假说，分析土地生态质量的空间分布特征，从而探究这一假说是否成立。

本文基于国土资源部科研专项经费项目“长三角经济发达地区土地生态状况调查与评估”项目提出的土地生态指数指标体系，计算得到金坛市行政村尺度的各准则层和综合层分值，首先通过土地生态指数热点分析验证城镇化对土地生态质量影响的规律；其次通过环形梯度分析计算土地生态指数距离衰减量；再通过条带梯度分析，对从城镇中心区辐射出的各个方向进行重采样，分析土地生态各准则层与综合层的分值状况，进一步探讨城镇化对土地生态的影响机理。

1 研究区概况

研究区金坛市（31°33′42″—31°53′22″N，119°17′45″—119°44′59″E）位于江苏省南部，是南京、上海和杭州三角地带的中心，属长三角经济发达地区，金坛市是我国农村综合实力百强市（县）。金坛全市总面积975.49km2，金坛东西相距约50km，南北相距约40km。金坛市市级行政区域界线主要涉及7个镇和1个开发区的区域单位，共165个行政村单元。2011年末，金坛市行政区域内人口54.9万人，近十年来人口年平均增长率为0.17‰，城镇化率达到49.60%，城镇和农村居民人均收入逐年上涨。金坛市经济快速发展过程中，土地生态问题凸显。

2 数据来源与方法

2.1 数据来源

构建土地生态指数的指标体系，从基础、胁迫、结构与效益角度构建4个准则层和15个指标（表1），权重确定方法选择主观赋权较为典型的层次分析法和客观赋权较为典型的熵权法综合得到各元指标层和准则层权重值［11－12］。研究时段静态指标基期年选取2011年，动态指标基期年选取2009—2011年。各指标计算公式及解释、数据来源见表1。根据土地生态指数指标体系，计算得到金坛市行政村单元各准则层和土地生态指数的值，进而在此基础上探讨金坛市土地生态质量分布的城乡梯度规律。

表1 土地生态指数指标体系及数据来源

2.2 研究方法

2.2.1 冷热点分析 空间热／冷点布局的分析旨在探究各指标类型属性值的“中心”（高值集聚，hotspots）和“洼地”（低值集聚，coldspots），主要基于空间计量地理学的空间自相关思想，其数理依据为Anselin在1995年所提出的空间局部自相关统计量（Local Indicators of Spatial Association，LISA），通过计算局部Getis－Ord G值定量判断区域极高或极低值要素的聚集情况。Getis－Ord表达式为

2.2.2 梯度分析模型 常用的城乡梯度分析，主要包括条带梯度分析（transect gradient analysis）和环形梯度分析（circle gradient analysis）两种［14］；为研究梯度分析中的各向异性问题，也有学者开展了四向与八向环形梯度分析［15］。条带梯度分析是以区域城镇中心为基准，取东西向或南北向或其他方向上特定长宽的连续单位空间，进而分析特定方向上研究对象的梯度变化；环形梯度分析也是以区域城镇中心为基础，以特定长度为增量半径向外做圆而等距环，进而分析由城镇中心向远郊乡村研究对象的梯度变化。条带梯度分析有利于揭示研究对象变化的细节信息，而环形梯度分异有利于揭示研究对象变化的一般趋势。因此，本文将首先运用环形梯度分析确定金坛市土地生态指数空间分异的城乡梯度变化规律和梯度力，再结合条带梯度分析城镇化对土地生态影响的细节。

环形梯度分析与条带梯度分析的步骤相似，包括城镇中心的确定，增量半径或样方大小的选取，以及各梯度研究对象信息统计对比三个主要步骤。其具体操作较为简单，关键在于城镇中心的确定和增量半径或样方大小的选取两步骤。已有相关研究中，对城镇中心的确定或回避［16］或以地区政府经纬度为代表，缺乏科学性与普适性。本文认为，选取区域中心商业区（Central Business District，CBD）范围，运用重心法计算城镇中心经纬度可以较好地满足相关研究及推广需求。

重心模型的基本原理为：假设某区域CBD由n个相对集聚的不连续平面空间单元构成，其中第i个单元的地理中心坐标为（Xi，Yi），Zi为该平面单元的占地面积，则CBD重心坐标（¯X，¯Y）为

重心模型实质上对离散要素空间坐标的面积加权，该过程不涉及主观判断，主要是相对的客观计算，可以基于ArcGIS软件平台等完成操作。此外，增量半径或样方大小问题在已有研究中也通常被忽视或回避［17－18］，选取的主观性较强。本文认为，尺度选择主要应依据信息解译原则，即所选择的研究尺度应足够小以致充分揭示研究对象的梯度变化规律，但尺度不宜比基础数据的“最粗”尺度小。据此，本文提出基于平均单元面积法确定所选尺度。具体公式为

式中：Si——基础数据的统计单元的面积；n——其数量。就本文而言，Si为金坛市165个行政村的相应面积，n即等于165。

3 结果与分析

3.1 土地生态指数冷热点分析

在ArcGIS 9.3中的Hot Spot Analysis下，进行土地生态指数的冷热点计算，可以得到其空间分异的冷热点图（图1）。

图1 土地生态指数空间冷热点分布

显而易见，金坛市县域范围内的土地生态指数呈现出的“核心—边缘”结构，即由城镇中心向远郊土地生态指数下降，主城区附近片区土地生态指数显著低值集聚，而外围远郊各村的区域土地生态指数则呈显著高值集聚。这一空间分异现象与空间经济学或新经济地理学所证实，由于规模报酬递增与交易成本所导致的区域发展过程中自生演化的城乡二元“核心—边缘”结构［19－20］相契合。

3.2 土地生态城乡梯度分析

土地生态指数热点分析表明金坛市土地生态质量具有城乡分布规律，为了进一步探讨金坛市细节上的土地生态质量城乡梯度规律，运用城乡梯度法对金坛市作土地环形梯度和条带梯度的分析。城乡梯度分析法在景观生态学及相关领域中得到广泛运用，并被证明为揭示事物空间分布模式及其人文影响因素的有效研究方法［20］。其基本思想是，通过空间换时间的策略，使得特定事物在不同发展阶段下的不同状况在同一地理范围内得到空间显化，从而可以通过进一步的对比研究发现事物的发展规律与驱动机制。

3.2.1 土地生态环形梯度分析 首先，在2011年金坛市土地利用数据库中，提取所有城镇用地，通过Fragstats定位出其最大斑块。在ArcGIS 9.3中，利用mean center工具定位出最大城镇斑块的重心，即得到CBD重心。其次，根据165个行政村面积统计，得到平均面积为5.42km2，则每个评价单元平均外切步长为2 328.83m。同时，基础数据中栅格数据基本重采样到30m空间分辨率。所以，环形梯度增量半径的选取满足为外切步长和基础数据空间分辨率的最小公倍数，即2 400m。最后，以CBD重心为中心，每2 400m为增量半径，通过multiple ring buffer向外作12个圈层以进行环形梯度分析（图2）。

图2 金坛市环形梯度分析示意图

对这12个从CBD中心向外拓展的各圈层所包含的行政区域的土地生态基础分值、土地生态结构分值、土地生态胁迫分值、土地生态效益分值和土地生态指数5个值进行面积加权，以得到各圈层可比的5个指标值，如图3所示。

图3 各圈层土地生态状况分值对比

研究发现：沿着金坛市“城镇—城乡结合部—近郊—远郊”梯度变化，从CBD中心向外，总体上区域土地生态指数不断升高；各圈层土地生态基础相差不大，由内向外呈逐步升高趋势，即就区域生态本底而言，外围地区略优于核心地区；金坛市外圈层比内圈层土地生态结构优良，城镇远郊地带较近郊地区明显要好；各区域由于直接人为驱动的生态胁迫程度不具备明显的梯度特征，而具有显著的局域特性；中远郊地带土地利用的生态效益优于近郊地带，但差别不大。

进一步地对土地生态指数与各圈层半径做回归分析，结果如图4所示，方程拟合显著，城乡空间梯度变化能解解释92.3%的区域土地生态指数差异。结果表明，沿地区CBD中心向外围地带的梯度来看，区域土地生态质量的改善非常显著，每2 400m向外总土地生态指数提高1.6个百分点，金坛市最边缘地带土地生态指数比其最核心圈层要高17.6个点。

图4 土地生态指数的城乡梯度拟合分析

3.2.2 土地生态条带梯度分析 为更好地研究金坛市土地生态指数城乡梯度细节，本节尝试建立如图5所示的8方向4条带分析梯度，并以金坛市CBD中心为零点，对每个样方编号。其中，东—西样带有17个样方，西11个东5个；南—北样带有14个样方，南8个北5个，但由于南半样带跨非评价区的钱资荡和长荡湖地区，导致样方南—4数据失真、样方南—5数据缺失；西北—东南样带有11个样方，西北6个东南4个；东北—西南样带有10个样方，东北2个西南7个。

同样利用面积加权求和法，计算每个样方的土地生态基础分值、土地生态结构分值、土地生态胁迫分值、土地生态效益分值和土地生态指数5个指数值进行比较，各样带统计曲线对比如图6所示。

观察4条土地综合分值曲线是否大致呈倒U型曲线，以及横坐标0值处曲线位于最低点，可见无论从哪个方向取条带梯度样方，从边缘地区向核心地区过渡，局地土地生态分值将显著下降，其原因在于生态结构的差异。图5表明东—西向、南—北向两条样带的综合分值曲线呈明显倒U型，西北—东南向和东北—西南向两样带的综合分值曲线也近0值范围内呈倒U型。和生态结构曲线的对比观察发现，综合分值曲线的整体走势主要取决于生态结构曲线。

图5 金坛市条带梯度分析示意图

图6 各样带土地生态状况分值对比

另一方面，观察四条土地综合分值曲线的反常值所属样方，对应图5找到其空间位置，根据其是否包含城镇建成区以及包含面积或与建成区临近程度，可以发现包含建制镇的样方区域，其局地土地生态质量将会出现异常波动，比正常状态要差。图中有东—7样方、西—3样方、南—2样方、南—4样方、北—3样方、北4样方、西南—3样方、西南—4样方、西北—2样方、西北—4样方和东南—3样方等计11个样方出现逆城乡梯度规律现象。结合图6分析可知，北4样方受边缘效应影响、南—2样方和南—4样方由于包含钱资荡水域或长荡湖局域而造成指数信息失真，其他逆城乡梯度样方除西北—4样方和西南—4样方外都主要受到城镇建成区的影响，这两个样方的反常可能是由于区域土壤酸化较严重导致的，需要后续研究进一步探究。

为定量揭示区域土地生态指数与各准则层分值的相关关系，以4条带51个样方为样本通过逐步回归，建立了如下作用方程。

式中：y——各样方土地生态指数；x1——生态结构分值；x2——生态胁迫分值。方程通过显著性检验，拟合精度超过90%，表明生态结构与生态胁迫是影响区域土地生态指数的主要因素，且由回归系数可知区域土地生态质量对土地生态结构的敏感性比对土地生态胁迫的敏感性更强。

结合以上分析，可以得出结论：从区域经济社会发展的核心地区向边缘地区过渡，局地土地生态质量具有显著改变的城乡梯度响应规律，其改变主要源于边缘地区的土地生态结构优于核心地区，而不同地区土地利用生态结构分值的不同则是导致局地土地生态质量波动变化的主要原因。

4 结论

（1）金坛市的土地生态指数冷热点分析显示，土地生态质量空间分异特征符合城乡梯度规律，土地生态指数呈现“核心—边缘”结构，城镇化是导致土地生态质量分布城乡梯度规律的重要原因。

（2）金坛市的环形梯度分析表明，土地生态指数沿着金坛市“城镇—城乡结合部—近郊—远郊”梯度变化，总体上随该梯度区域土地生态指数不断升高，沿城镇中心向外每增加2 400m土地生态指数提高1.6%。

（3）金坛市的条带梯度表明土地生态质量分布的城乡梯度规律：土地综合分值曲线大致呈倒U型曲线，最大城镇斑块处曲线位于最低点；包含建制镇的样方区域土地生态质量处于低值，主要原因是受到城镇建成区的影响，由此说明不仅城镇中心土地生态质量状况处于最低值，而且城镇所在地土地生态质量都会比正常状态差。区域土地生态指数与各准则层分值之间的回归分析表明区域土地生态质量主要影响途径是通过对区域土地生态结构的影响和经济社会发展过程中对土地利用的人文胁迫作用，且主要取决于结构分值的影响。

（4）本文对土地生态指数空间分异规律作了深入探究，而对土地生态指数空间分异的驱动力分析和主控因子识别将在进一步研究中进行探索。

［1］ 戴靓，姚新春，周生路，等.长三角经济发达区金坛市土地生态状况评价［J］.农业工程学报，2013，29（8）：249－257.

［2］ 刘焱序，李春越，任志远，等.基于LUCC的生态型城市土地生态敏感性评价［J］.水土保持研究，2012，19（4）：125－130.

［3］ 钟晓兰，周生路，李江涛，等.长江三角洲地区土壤重金属污染的空间变异特征：以江苏省太仓市为例［J］.土壤学报，2007，44（1）：33－40.

［4］ 李保杰，顾和和，纪亚洲.矿区土地复垦景观格局变化和生态效应［J］.农业工程学报，2012，28（3）：251－256.

［5］ 许月卿.基于生态足迹的北京市土地生态承载力评价［J］.资源科学，2007，29（5）：37－42.

［6］ 黄海，刘长城，陈春，等.基于生态足迹的土地生态安全评价研究［J］.水土保持研究，2013，20（1）：193－196.

［7］ 许国平.中国土地资源安全评价研究进展及展望［J］.水土保持研究，2012，19（2）：276－279.

［8］ 孙洪波，杨桂山，苏伟忠，等.沿江地区土地利用生态风险评价：以长江三角洲南京地区为例［J］.生态学报，2010，30（20）：2253－2260.

［9］ 刘孝富，舒俭民，张林波.最小累积阻力模型在城市土地生态适宜性评价中的应用：以厦门为例［J］.生态学报，2010，30（2）：421－428.

［10］ 徐美，朱翔，刘春腊.基于RBF的湖南省土地生态安全动态预警［J］.地理学报，2012，67（10）：1411－1422.

［11］ 李晓倩，刘学录.土地评价指标权重赋值方法的比较分析［J］.甘肃农业大学学报，2012，27（5）：129－133.

［12］ Chen－jing FAN，Shi－guang S，Si－hui W，et al.Research on urban land ecological suitability evaluation based on Gravity－Resistance Model：A case of deyang city in China［J］.Procedia Engineering，2011，21：676－685.

［13］ 杨宇，刘毅，金凤君，等.塔里木河流域绿洲城镇发展与水土资源效益分析［J］.地理学报，2012，67（2）：157－168.

［14］ Luck M，Wu J.A gradient analysis of urban landscape pattern：a case study from the Phoenix metropolitan region，Arizona，USA［J］.Landscape Ecology，2002，17（4）：327－339.

［15］ Aithal B H，Sanna D D.Insights to urban dynamics through landscape spatial pattern analysis［J］.International Journal of Applied Earth Observation and Geoinformation，2012，18：329－343.

［16］ 张利权，吴健平，甄彧，等.基于GIS的上海市景观格局梯度分析［J］.植物生态学报，2004，28（1）：78－85.

［17］ Hahs A K，McDonnell M J.Selecting independent measures to quantify Melbourne′s urban－rural gradient［J］.Landscape and Urban Planning，2006，78（4）：435－448.

［18］ 尹海伟，孔繁花.济南市城市绿地时空梯度分析［J］.生态学报，2005，25（11）：218－226.

［19］ 宋成舜，周惠萍，刘成武，等.鄂东南丘陵地区优化城乡用地结构和布局研究：以湖北省崇阳县为例［J］.水土保持研究，2010，17（1）：72－76.

［20］ 李俊祥，王玉洁，沈晓虹，等.上海市城乡梯度景观格局分析［J］.生态学报，2004，24（9）：1973－1980.