苏南沿江城镇扩展的空间模式及其测度_马荣华
第62卷第10期2007年10月
地理学报
ACTAGEOGRAPHICASINICA
Vol.62,No.10
Oct.,2007
苏南沿江城镇扩展的空间模式及其测度
马荣华1,顾朝林2,蒲英霞3,马晓冬4,朱传耿4
(1.中国科学院南京地理与湖泊研究所,南京210008;2.清华大学城市规划系,北京100084;3.南京大学地理与海洋科学学院,南京210093;4.徐州师范大学城市与环境学院,徐州221009)
摘要:定义了城镇边界的划分规则,通过卫星遥感影像获取了苏南沿江地区1984年、1991年、2000年和2005年的城镇边界,利用分形维数(半径维数、网格维数、相关维数、边界维数)、紧凑度指数以及空间自相关指数(GlobalMoranI、LocalMoranI),定量揭示了区域城镇扩展的空间模式。结果表明:分形度量可以定量揭示区域城镇的空间分布特征以及空间组织结构,空间自相关可以定量揭示城镇的空间扩展过程、聚簇模式,二者相互结合,可以测度区域内城镇扩展的空间模式,揭示城市群/城市连绵区的形成过程;苏南沿江地区城镇扩展从改革开放初期自发的以大城市为中心的点状扩展模式,经历了由政策驱动下以融合填充为主导的点-轴延伸模式,发展到以经济辐射联系为策动力基于区域产业空间组织的城镇连绵集聚模式,一个大的城市群正在形成。特别自2000年以来,城市群的空间结构特征日益显现。区域内城市的多中心现象始终存在,但空间极化作用逐渐削弱。关键词:城市化;空间扩展模式;分形;空间自相关;江苏
1引言
诺贝尔经济学奖获得者JosephStiglitz曾说过:“中国的城市化和美国的高新技术发展是21世纪深刻影响人类发展的两大关键因素…,中国的城市化是推动区域经济发展的发动机”城市化过程”理论,城市形态作为城市化的结果,空间模式的。根据Knox的“
演变明确地反映了城市化过程[1]。中国的城市化过程既不与发达国家相吻合,也不同于第三世界国家(工业化明显滞后于城市化),还没有完全适合中国城市化模式和机理研究的理论可供借鉴[2]。在不同国家的不同发展阶段,城镇扩展的空间模式不同,一般遵循多中心—中心集聚—中心扩散—相对平衡的演化过程[3]。研究城镇扩展的空间模式,可以帮助我们正确把握城市化进程所处的阶段,预测城市的未来发展趋势。
城市地理学对城镇扩展空间模式的研究大多以扩展动力分析为基础,从资本、产业、信息等角度入手,结合社会经济统计资料,以定性和半定量分析为主[4-8]。随着遥感和GIS等技术手段的广泛应用,土地利用/土地覆被变化为城市空间扩展研究提供了新的视角。借鉴土地利用/土地覆被变化的定量化测度方法或指标[9],城市扩展的空间模式研究逐渐走向定量化。目前普遍使用扩展速度、扩展强度以及边界维数等指标,对比分析不同时期的扩展程度和城市形态特征,进而揭示各个阶段的扩展模式。这些分析大多集中在单一城市,如北京[10,11]、上海[12]。也有学者使用上述方法研究城镇群体的扩展过程,如北京地区[13]、常熟地区[14]、南通地区[15]、长江三角洲地区[16],但对城镇群体扩展的空间组织演化与模式发现还较少涉及,而空间数据挖掘方法为这一研究提供了新的技术支撑[17,18]。城镇形态作为城镇空间扩展的表征,是定量测度城镇扩展模式的重要载体。城市形态具有
收稿日期:2007-01-15;修订日期:2007-05-28
基金项目:国家自然科学基金项目(40435013;40601074;40671053)[Foundation:NationalNaturalScienceFoundation
ofChina,No.40435013;No.40601074;No.40671053]
作者简介:马荣华(1972-),男,山东费县人,博士,副研究员。E-mail:rhma@niglas.ac.cn
1011-1022页
1012
地理学报62卷
明显的分形特征[19],城市的空间扩展模式可以通过分形的幂定律来度量[20]。区域城镇的空间自组织也具有无标度性,显示出分形特征[21,22]。同时,城镇实体地域的演化在区域范围还具有空间自相关的特征,并呈现出不同尺度、不同方位上的空间趋同和空间异质的结构,这种特性可以揭示区域内城市聚簇(城市群)组织演化的空间模式。本研究以中国经济发达地区—苏南沿江地带内的所有城镇为研究对象,通过不同时期的卫星遥感影像获取城镇实体地域的形态,通过分形和空间自相关分析等方法,度量区域内不同时期的城市空间分布状况及其空间扩展模式,为经济发达地区城市化的空间格局与过程提供定量化研究方法,为世界上不同区域的城市空间扩展提供可比较的典型研究案例。
2数据与方法
2.1数据来源
研究区为江苏苏南濒临长江的市县(区),行政范围包括常州市区、无锡市区、苏州市区、江阴市、张家港市、常熟市、太仓市、昆山市(图1),面积约9200km2,是长江三角洲地区的核心区。该区域是中国城市化水平最高的地区之一,也是20世纪80年代以来,中国快速城市化的典型地区之一。目前,该地区有3座特大城市(常州、无锡和苏州)、5座中等城市(江阴、张家港、常熟、昆山和太仓),1万人以上的建制镇或大的居民点205个。
改革开放以来,中国经历了三个重要的发展阶段。第一阶段开始于20世纪80年代初期,源于邓小平提出的“改革开放”政策;第二阶段始于90年代初期,源于社会主义市场经济体制的确立,特别是邓小平的“南巡讲话”有力的推动了中国的对外开放和市场化改革;第三阶段开始于2000年,经济全球化使长三角、珠三角等我国东部沿海地区成为FDI的密集区,特别是上海、广州等核心城市成为跨国公司全球营销体系的重要结点,上海城市战略定位的调整使长三角区域发展从中心城市单向扩散和单一动力主体模式,朝着多层次、多维向、多中心、多动力联动整合模式转变[23],上海对它周围的地区显示出越来越强的影响力。根据这三个发展阶段,选择1984年、1991年、2000年和2005年作为分析的时间节点。另外,为了避免行政境界变迁对分析结果的影响,使用1991年城镇行政边界作为基本的分析和计算单元。统计数据引自江苏省公安厅汇编、1985年4月出版的江苏省人口资料(1984)以及中国统计出版社分别于1992年、2002年、2006年出版的江苏统计年鉴(1992)、江苏统计年鉴(2001)、江苏统计年鉴(2005)。
在城市空间扩展的研究过程中,大都把城市看作一个具有明确边界、内部均质的空间实体,城市边界成为城市空间扩展的重要载体,是城市扩展的一个重要指标。因此,准确定义城市空间范围、确定城市边界就显得尤为重要,但比较困难,并且没有统一的标准[24]。我们根据提取城镇建设用地的遥感影像的空间分辨率,同时考虑到该区域水网密集和高强度工业开发的显著特点,把下述空间或土地利用类型看作城镇建设用地,并入本研究中城镇实体地域的范畴:①完全被其它城市土地利用类型(如居住地、工业用地、商业用地)所包围的开放空间(如水体、绿地);②紧靠城镇的工业开发区或经济技术开发区(包含这些开发区的建成区以及由建成区包围的非建成区);③与城镇紧密相连的村镇。另外,如果两个或几个城镇因扩展而连为一体,城市扩展后面积按照行政边界切割计算。
使用空间分辨率分别为30m和23m的LandsatMSS/TM/ETM(1984年8月4日、1991年7月23日、2000年5月4日)和IRS-P6(2005年5月3日)卫星影像数据,利用1:50000地形图,通过最近邻算法,进行几何校正,误差在1个象元之内;之后在ERDAS支持下进行图像拉伸、增强,突出纹理特征,通过目视判读获取城市建成区数据(图2)。
10期马荣华等:苏南沿江城镇扩展的空间模式及其测度1013
通过ArcGIS编辑处理,并转换为bmp文件(4724×4724个像元,250m×250m的像元尺寸),便于分形分析。2.2研究方法
不2.2.1分形维数
同类型的分形维数可以定量度量区域内城镇的空间分布特征,进而揭高速公路
铁路示城镇的空间组织结
长江岛屿构。基于分形维数测度
河流湖泊
方法的城市空间扩展可
省界km
以为世界上不同区域的研究区城市分类提供典型案
图1研究区位置
例,分形维数已成为城
Fig.1Locationofthestudyarea
市形态对比的一个实用
指标[19],包括半径维数、网格维数、相关维数和边界维数。前3种维数归属计数法,遵循逐步叠代原理,主要涉及两个变量:计数窗口内计算对象(黑色像元)的数量(N)以及计数窗口的大小(! ),有
N=" D或N=#-D(1)
式中:D是分维值。一般而言,按下式估测分维数:
N=a$D+c(2)
式中:a被称作形状因子(pre-factoroftheshape),是局部偏离的综合指标,理论上a取值为1,一些情况下等于0.5或3;如果其值超过10或小于0.1,则不具有分形性质[19]。通常情况下,可以使用非线性回归获取最佳拟合曲线,以获取分形维数D,拟合优度通过相关系数来判断,如果相关系数小于0.9999,拟合曲线与理论曲线不能很好吻合,表
明研究对象不具有分形性质或者具有分段分形的性质。
半径维数(Dr)表明整个区域内城镇从中心(计数中心)到外围的空间分布衰减特征。如果Dr小于2,表明区域城镇的空间分布在密度上从中心向外围衰减;如果等于2,表明在半径方向上均一;如果大于2,表明从中心到外围增加。网格维数(Dg)表示一定区域范围内城市分布的均匀度,理论上位于0到2之间。如果Dg等于0,表明所有城镇集中于一点,
1984年城镇
1984-1991年扩展1991-2000年扩展2000-2005年扩展
km
图21984-2005年研究区内城镇的空间扩展分布
Fig.2Spatialdistributionofcitiesandtownsinthestudyareain1984-2005
1014
地理学报62卷
即仅有1个城市;如果等于1,表明城镇沿线(如铁路、河流等)集聚分布;如果等于2,表明城镇在整个区域内均匀分布;一般地,Dg位于1到2之间,随着Dg值的不断增加,城镇的空间分布均质水平愈显著。与Dg、Dc一样,相关维数(Dc)也表明一定区域内城镇分布的均匀度,但Dc可以揭示更多的细节。一般而言,Dc位于0到2之间;如果接近2,表明城镇分布较为均一;如果接近但大于0,表明区域内存在首位城市。边界维数又称面积周长维数;如果城市是一个简单的几何形状,其边界维数(线)为1,表面维数(面)为2,面与线的比率约为1.05[19],这与欧几里德几何学相矛盾,但与分形几何学相一致。对每一个表示城镇的多边形而言,周长(P)和面积(A)之间存在如式(3)的基本分形关系[25]:
P=kAD/2
式中:D是分形维数(Da),k是比例常数。上式可以转换为:
(3)
lnA=(2/Da)×lnP+c(4)
式中:c是线性回归截距。一般而言,Da位于1到2之间。
本研究使用上述4种分维方法,从不同的视角揭示经济发达地区不同阶段的城镇空间分布和组织体系。另外,使用标度曲线(ScalingBehaviorCurve,SBC)! (" )帮助识别城
镇聚簇形态变化的空间阈值,以揭示城市集聚在不同规模下的标度行为。
借用景观生态学中的紧凑度,定量测度区域内城镇的整体空间集聚2.2.2紧凑度指数
形态。一般使用修正的紧凑度指数(revisedcompactnessindex,CI')[26]:
CI'=
! 2" i
i
i
N
2
(6)
式中:Si和pi是第i个城市(包括城镇)的面积和周长,pi是第i个城市内接圆的周长,N是城市的总个数。
扩展强度指数(SprawlIntensityIndex,SII)定义为[14]:2.2.3扩展强度
SIIi=
Asi
×100
ti×(7)
式中:Ati是第i个城镇的镇域总面积,Asi是第i个城镇在时间内镇区扩展的面积。
区域城镇群体在其生长发育的过程中所表现出的空间格局或组织是2.2.4空间自相关
由城镇社会经济发展的空间关联特质引起的,从而造成不同层级城镇之间的空间相互作用,进而形成区域城市化进程的某种空间模式,单一要素之间的模式可以通过空间自相关校验的方法来定量度量。GlobalMoranI(GMI)和LocalMoranI(LMI)[27,28]可以定量测度城镇的聚簇模式,识别城镇扩展“热点区”的分布,进而探测城镇扩展的空间模式[15,18]。GMI定义为:
I=n
0
! ! w(x
i
j≠i
ij
nn
i
-xj-)-)2
! (x
i
(8)
i
式中:n是样本总数,xi和xj分别是位置i和j处的观测值(本研究中指不同城镇的扩展强度),是观测值xi在所有位置处的平均,wij是空间权重矩阵(n×n),S0是空间权重矩阵wij中所有元素之和。利用MoranI可以测度空间自相关性,发现观测值在空间分布的差异性和相关性。当位于一定距离d内的观测值相近时,MoranI显著而且为正,不相近时为负,当观察值随机排列时为零。GMI是对观测值空间模式的整体定量描述。另外,在全局空间随机分布的样本中,也可能存在局部空间自相关的观测值。因此在全局分析
10期马荣华等:苏南沿江城镇扩展的空间模式及其测度1015
的基础上,采用LMI进行局部自相关的测度。LMI定义为:
Ii=
! w'
ij
ZiZj
(11)
式中:Zi和Zj是观测值的标准化值,w'ij是wij的行标准化。在显著性水平下(如p值小于0.05),当Ii和Zi均为正时,表明位置i处的观测值和它周围的观测值均为高值区,即高高集聚(HH);当Ii为正、Zi为负时,表明位置i处的观测值和它周围的观测值均为低值区,即低低集聚(LL);当Ii为负、Zi为正时,表明位置i处的观测值大于它周围的观测值,即高低集聚(HL);当Ii和Zi均为负时,表明位置i处的观测值小于它周围的观测值,即低高集聚(LH)。
3结果分析
苏南沿江地区2005年的城镇总面积分别比1984年、1991年、2000年增长了11.42、4.34、2.41倍(图3a)。该区城镇总面积(TUA)和城镇总人口(TUP)之间的线性关系非常显著(R2=0.9992)(图3b);从每万人拥有的城镇面积来看,属于异速增长;城镇建设用地的增长速度快于城镇人口的增长,表现出与人口增长并不协调的城镇圈地式增长模式。3.1异质性、紧凑度和中心性
以区域质心[310335]、常州[113229]、无锡[253322]和苏州[388441]作为计数中心,分析半径维数、网格维数和相关维数,其中相关维数的计算选择方形计数窗口。
以区域质心3.1.1城镇空间分布整体上逐渐趋于均匀,局部范围内异质特征较为显著
为中心的全局分形半径维数(GFRD)的分析表明(图4a):整体上GFRD值逐渐增大,区域内城镇的空间分布逐渐趋于均匀;仅1991年GFRD出现异常(小于1),表明1991年及其前后城镇扩展的无序化程度加剧,区域内城镇的空间组织形态类似谢尔宾斯基地毯的形状[20]。1984年、1991年和2000年的SBC曲线变化趋势相似(图5a),沿半径方向城镇的空间分布异质性较大;但2005年变化明显,区域内城镇在向着区域质心的方向扩展,空间分布的同质性增加,尤其在半径400像元(100km)内,具有更好的紧凑度。分别以苏州、无锡和常州为中心的GFRD分析表明(图5b,5c),从1984年到2005年,稀释半径范围不断向外推移,表明三个城市的中心范围在不断扩大,促使区域内城镇的空间分布整体上向着均匀、同质的方向发展。1984年的全局相关维数(GFCD)约为1(图6a),局部相关维数(LFCD)随计数窗口的像元尺寸在0.55到1.43之间的变化(图6b),特别当计数窗口尺寸范围在23到86、118到130以及130到154个像元内时,LFCD均小于1,
T U A (k m )
(k m )
(a)
TUP (×10)
图3不同年份的城镇总面积及其与城镇总人口的关系
Fig.3Relationshipbetweenthetotaltownareaandthetotaltownpopulationindifferentyears
1016
地理学报
62卷
图4全局半径维数及其对应的相关系数:(a)以区域质心[310335]为中心;(b)以苏州[388441]
为中心;(c)以无锡[253322]为中心;(d)以常州[113229]为中心
Fig.4Globalfractalradiusdimensionsandthecorrespondingcorrelationcoefficientscenteredat:(a)regioncentroid[310335],(b)Suzhou[388441],(c)Wuxi[253322]and(d)Changzhou[113229]
a
a
图51984-2005年半径分维的尺度变化曲线SBC:(a)以区域质心[310335]为中心;(b)以苏州[388
441]为中心;(c)以无锡[253322]为中心;(d)以常州[113229]为中心
Fig.5Scalingbehaviorcurves(SBC)offractalradiusdimensionin1984,1991,2000and2005,centeredat:(a)regioncentroid[310335],(b)Suzhou[388441],(c)Wuxi[253322]and(d)Changzhou[113229]
10期马荣华等:苏南沿江城镇扩展的空间模式及其测度1017
1.6
1.000 0.999 0.998 0.997
2.0
1.4 1.3 1.2 1.1 1.0
1984
1991
2000
2005
a
1.5
1.81.61.41.21.00.845161 [***********]91
( )
图61984-2005年区域城镇的相关维数及其尺度变化:(a)相关维数;(b)尺度曲线SBC
Fig.6Correlationdimensionsandtheircorrespondingscalingbehaviorcurves(SBC)in1984,1991,2000and2005
区域内城镇空间组织的异质性特征更为明显,这与城 镇空间上的相互分离相对应。1991年和2000年的GFCD大于1但小于1.3,SBC的变化趋势基本相同,
稍有差异的是,1991年计数窗口尺寸! 在28到90以
及112到126个像元范围内时LFCD小于1,而2000
年均大于1,表明2000年城镇空间组织在上述标度范
围内优于1991年。2005年的相关维数最高,表明
2005年的城镇空间分布更为均质。紧凑度指数的分析
表明(图7),从1984到2005年城镇空间分布越来越
图7修正的紧凑度指数
紧凑、连接越来越紧密,验证了分形分析的相关分析
Fig.7Revisedcompactnessindex
结果;与其不同的是,分形边界维数的分析表明(图8),城镇形态的复杂性在增加,但波动较大,1991年
前后最高,2000年、2005年有所下降。1991年较多
的乡镇新增建设用地出现,到2000年以后城镇快速扩展导致城镇空间粘合过程显著,形态的稳定性又有所增强。格网分维分析同样揭示出:区域内城镇空间
分布在1984-2005年间整体上逐渐均一化,局部异质
性特征也较为明显,围绕不同的中心城市,城镇的空
间组织结构存在差异。篇幅关系,不再对格网维数进
行详细分析。
图8边界分形维数
3.1.2区域内城市三中心特征始终存在,但中心极化Fig.8Fractalboundarydimension
以苏州为中心作用逐渐削弱,有走向融合的趋势
的SBC曲线表明(图5b):①4个年份具有相似的空间组织结构,存在明显的中心和外围差别,中心城市的特征显著;②1984-2000年,苏州的城市中心范围相近,具有相同的稀释半径范围(12.5 ̄47.5km),最大有效稀释半径(47.5km)远离无锡和昆山;③2005年苏州的城市中心范围明显扩大,稀释半径范围(35 ̄60km)向外推移,接近无锡和昆山的城市边界。以无锡为中心的SBC曲线表明(图5c):①1984-2000年,无锡的城市中心范围相近,在50到160个像元(12.5到40km)的范围内存在较为明显的稀释带,稀释带的最大有效半径(40km)均远离苏州和常州,稀释带内LFRD都小于1,表明稀释带内的城镇空间分布存在明显的异质特征;②2005年SBC相对平缓,城市中心范围扩大,稀释带仍存在,但范围变窄(约31到50km)并向外推移,最大有效半径(50km)接近苏州
CC T '?
1018
地理学报62卷
和常州边界,稀释带内LFRD为1.489,表明稀释带内城镇空间分布较为均质;③以无锡为中心的城镇空间组织结构存在明显的中心和外围差别,中心城市的特征较为明显。以常州为中心的SBC曲线表明(图5d):①城市中心范围不断向外围扩展,稀释半径逐渐外移,已接近江阴和无锡;②以常州为中心的城镇空间组织结构存在中心和外围差别,中心城市的特征较为明显。3.2空间扩展模式
2000年以来,苏南沿江区城镇扩展的强度迅速增加,约是20世纪80年代中末期的6.8倍。三个时段内(1984-1991年、1991-2000年、2000-2005年)城镇扩展强度的平均值分别为0.54、0.34和3.65。当以拓扑邻接关系构建空间权重矩阵时,上述三个时段内扩展强度SII的全局MoranI值(GMI)分别为0.427、0.176和0.294;当以5km(最佳邻域距离)和10km为相关距离构建空间权重矩阵时,三个时段的计算结果如表1。
表1利用不同的邻域距离构建的空间邻接矩阵计算所得的空间扩展强度的全局MoranI值(GMI)Tab.1GlobalMoranI(GMI)ofspatialsprawlintensitycalculatedbydifferentspatialadjacencymatrix
basedonadjacencydistanceof5kmand10km,respectively
I(d) E(d) Z Score
1984-1991 1991-2000 2000-2005 1984-1991 1991-2000 2000-2005
1.620 -0.005 15.317
0.200 -0.005 1.878
0.456 -0.005 4.165
0.707 -0.005 18.517
0.167 -0.005 4.363
0.222 -0.005 5.666
三种计算结果均表明:①1984-1991年城镇扩展强度的聚集度最高,其次是2000-2005年,1991-2000年最低;②1984-1991年城镇的空间扩展主要集中表现在少数城市;③1991-2000年扩展强度的聚集度明显下降,城镇扩展表现出一定的离散性;④与1991-2000年相比,2000-2005年扩展强度的聚集度增加,空间分布的异质性有所增强,空间上有较大扩展的城镇的数量增加。扩展强度的局部MoranI值(LMI)的计算结果表明,区域内空间扩展的局部聚集模式存在较大差异(图9a-9c):①三个时段内,呈现出明显的HH集聚区逐渐扩展、增生,LL区域不断缩小的过程,HH区集聚规模、分布的位置和集聚形态的时相演化特征也十分明显;②1984-1991年,HH区为鲜明的中心地结构,集中在常州、无锡、苏州市区周围,其余大部分城镇位于LL集聚区,表明这个时段内城镇发展主要为核心集聚型;③1991-2000年,HH区开始增多,多个HH集聚核开始出现,特别是无锡的江阴地区、苏州的昆山地区成为新的HH集聚区,而且随着无锡及江阴的HH区扩展,研究区的西北部即无锡、常州地区呈现HH连接区;④2000-2005年,在苏州地区临近上海的太仓、昆山附近的城镇形成连片,形成大的HH集聚区,并向西与苏州市区、无锡市区粘连为非常连续的HH区集聚带,常州地区的HH区逐渐缩小;⑤整体上从城镇扩展的HH区演化上看,研究区的城镇发展经历了一个中心集聚(1984-1991年),即城镇的快速扩展区集中在少数几个中心城市,苏州、无锡、常州三个地级中心城市主导着研究区的城市化进程。然后多点集片(1991-2000年),即多个城镇快速扩展区出现,并与原快速扩展区的外围集结成片,特别是无锡各市县的城镇发展表现出强势的特征,并辐射影响到常州,有与常州连接成片的趋势。最后连片成带(2000-2005年),即城镇快速扩展区不断扩张,并与新增的城镇快速扩展区相互粘合,逐渐演变为较连续的带状集聚区,在此过程中,苏州各市县的发展步入快车道,替代无锡成为研究区城市化进程的引领者。3.3模式驱动机制
城市(城镇)扩展是一个复杂的过程,最终表现出不同的空间形态。历史形态学派、区位经济学派、决策行为学派和政治经济学派对城市扩展/城市增长模式都有不同的理
10期马荣华等:苏南沿江城镇扩展的空间模式及其测度1019
解[9]。从一个城市整体的角度,城市的空间扩展呈现轴向增长、同心圆增长、扇形扩展以及多核增长等多种形态,历史学派认为圆形城市是城市形态和增长模式的理想类型[9];考虑城市形态的形成过程,城市增长一般表现为填充式、扩张式、独立式、线状式、簇状式等五种方式[29]。从区域的角度,一般表现为城市群、都市圈、城市连绵区、大都市伸展区等四种方式,区域内的重要交通轴线在城市扩展过程中起着重要作用,往往成为城市扩展的牵引方向。中国的城市化属于双轨城市化模式[30,31],即政府主导城市化(自上而下)和自发城市化(自下而上);前者以非农人口的增加为特征,后者以小城镇城市化和流动人口迁移为特征,二者共同决定了城市化的空间模式。改革开放以来,包括资本和技术的流动、生产的全球重构与转移等要素在内的经济全球化为城市化带来了新的动力[4-6,8,32],促使江苏城市化动力机制由自上而下和自下而上转变为多元城市机制,外资成为江苏城市化新的动力机制[33]。苏南一直是我国城镇比较密集的地区,特别进入1990年代以来,苏南地区设立了大量的吸纳工业集聚的经济开发区,多数开发区的规模超过了原城镇建成区,这些开发区或在城镇的边缘布局,或依托城镇对外交通通道形成组团或新城,经济开发区往往成为城镇的新增长空间和扩张牵引方向,城镇实体地域的快速城市化极大地推动了区域城镇群体形态的演化,尤为明显地表现在沿交通线的各级经济开发区在快速扩张中相互粘合。进入新世纪以后,以上海为核心的长三角地区的城市化进程进入功能提升阶段,一方面许多开发区由单纯的工业区向综合城区演化,进入开发区的功能城市化阶段;另一方面,全球化浪潮使得世界城市以其高新技术及其创新升级产业结构、促进城市集群、推动区域一体化的内生性而产生的辐射力、聚合力,在配置区域资源方面的功能将愈益强化。苏南沿江区,尤其
(a)1984-1991年
(b)1991-2000年
(c)2000-2005年
图91984-2005年LMI散点地图Fig.9ScattermapoflocalMoranI(LMI)ofurbanspatialsprawlintensityin1984-2005
1020
地理学报62卷
是各市县中心城市,再一次抓住了上海为打造世界城市而进行战略调整的机遇,进一步提高了与上海都市区经济的产业耦合与集聚化程度,并通过技术互补、产业配套、要素配置,加强了城市间的经济联系,在不断促进沿重要经济联系走廊上工业集聚带或城镇连绵带生长的基础上,进一步推动了区域城镇集群和大都市伸展区的发育。
4结论
(1)1984年以来,苏南沿江区城镇的总面积扩大了近12倍;2000年以来,扩展强度迅速增加,增加了近7倍。区域内城市的多中心现象始终存在,但空间极化作用逐渐削
弱。
(2)总体上,苏南沿江区的城镇空间扩展从改革开放初期自发的以大城市为中心的点状扩展模式,经历了由政策驱动下以融合填充为主导的点-轴延伸模式,到以经济辐射
联系为策动力基于区域产业空间组织的城镇连绵集聚模式,一个大的城市群正在形成。特别自2000年以来,城市群的空间特征日益显现,区域内城镇的空间分布逐渐均质化、紧凑化。
(3)改革开放以来,苏南沿江区的发展热点不断转换,常州始终处于发展的次热点区。最初,常州、无锡和苏州三个地区的发展较为均衡;1990年代后,无锡地区率先成为苏南沿江区的发展热点,常州也取得了较快发展;2000年代以后,苏州成为苏南沿江区的新热点。
(4)区域内城镇的空间分布具有分形特征,分形度量可以定量揭示区域城镇的空间分布特征,辅助空间组织结构分析;空间自相关可以定量揭示城镇的空间扩展过程、聚簇模式和热点转换;二者相互结合,可以测度区域内城镇的空间扩展模式,揭示城市群/城市连绵区的空间形成过程。
(5)分形计算中的标度曲线分析非常重要,是揭示围绕不同中心城市的不同城镇空间组织结构的关键;扩展强度是进行城镇空间扩展分析的一个实用标准化指数,是利用空间自相关的测度方法探测城市空间过程、聚簇模式和热点转换的基础,其中的关键指数是LocalMoranI和LocalGetis-OrdG,关键的可视化表现手段是Moran散点地图。参考文献(References)
[1]KmoxPL.Urbanization:AnIntroductiontoUrbanGeography.EnglewoodCliffsNJ,1994.
[2]GuChaolin.Urbanizationstudies:Aninternationalapproach.CityPlanningReview,2003,27(6):19-24.[顾朝林.城市
化的国际研究.城市规划,2003,27(6):19-24.]
[3]GuChaolin,ZhenFeng,ZhangJingxiang.AgglomerationandDiffusion:NewTheoryonUrbanSpaceStructure.
Nanjing:SoutheastUniversityPress,2000.[顾朝林,甄峰,张京祥.集聚与扩散:城市空间结构新论.南京:东南大学出版社,2000.]
[4]NingYuemin.Newprocessofurbanization:DynamicsandfeaturesofurbanizationinChinasince1990.Acta
GeographicaSinica,1998,53(5):410-477.[宁越敏.新城市化进程:90年代中国城市化动力机制和特点探讨.地理学报,1998,53(5):470-477.]
[5]ZhouYixing,CaoGuangzhong.CityPlanningReview,1999,23(12):8-13.[周一星,曹广忠.改革开放20年来的中国
城市化进程.城市规划,1999,23(12):8-13.]
[6]GuChaolin,ZhangQin,CaiJianming.EconomicGlobalizationandUrbanDevelopmentinChina:ResearchonUrban
DevelopmentStrategiesintheNewCentury.Beijing:TheCommercialPress,1999.[顾朝林,张勤,蔡建明.经济全球化与中国城市发展:跨世纪城市发展战略研究.北京:商务印书馆,1999.]
[7]XiaoYizhuo.Spatialgrowthpatternofinternationalcities:CasestudiesofTokyoandParis.UrbanProblems,2003,(3):
30-33.[肖亦卓.国际城市空间扩展模式.城市问题,2003,(3):30-33.]
[8]GuChaolin.NewtrendsofChina'surbandevelopment.CityPlanningReview,2006,30(3):26-31.[顾朝林.中国城市
发展的新趋势.城市规划,2006,30(3):26-31.]
10期马荣华等:苏南沿江城镇扩展的空间模式及其测度1021[9]LiuShenghe,HeShujin.Aspatialanalysismodelformeasuringtherateoflandusechange.JournalofNaturalResources,2002,17(5):533-540.[刘盛和,何书金.土地利用动态变化的空间分析测算模型.自然资源学报,2002,17(5):533-540.]
[10]GuChaolin.Studyonthelanduse/coverchangemechanisminBeijing.JournalofNaturalResources,1999,14(4):
307-312.[顾朝林.北京土地利用/覆盖变化机制研究.自然资源学报,1999,14(4):307-312.]
[11]LiuShenghe,WuChuanjun,ShenHongquan.AGISbasedmodelofurbanlandusegrowthinBeijing.Acta
GeographicaSinica,2000,55(4):407-416.[刘盛和,吴传钧,沈洪泉.基于GIS的北京城市土地利用扩展模式.地理学报,2000,55(4):407-416.]
[12]LiXiaowen,FangJingyun,PuShilong.TheintensityandmodesofurbanlandusegrowthinShanghai.Journalof
NaturalResources,2003,18(4):412-422.[李晓文,方精云,朴世龙.上海城市用地扩展强度、模式及其空间分异特征.自然资源学报,2003,18(4):412-422.]
[13]HeChunyang,ShiPeijun,ChenJinetal.ProcessandmechanismofurbanizationinBeijingarea.ActaGeographica
Sinica,2002,57(3):363-371.[何春阳,史培军,陈晋等.北京地区城市化过程与机制研究.地理学报,2002,57(3):363-371.]
[14]MaRonghua,ChenWen,ChenXiaohuietal.TownandconstructionlandgrowthofChangshucity.ActaGeographica
Sinica,2004,59(3):418-426.[马荣华,陈雯,陈小卉等.常熟市城镇用地扩展分析.地理学报,2004,59(3):418-426.]
[15]ChuJinlong,MaXiaodong,GaoShuetal.Sptiao-temporalcharacteristicsofurbanizationareagrowthinNantongCity,
JiangsuProvince.JournalofNaturalResources,2006,21(1):55-63.[储金龙,马晓冬,高抒等.南通地区城镇用地扩展时空特征分析.自然资源学报,2006,21(1):55-63.]
[16]HeJianfeng,ZhuangDafang.GeographicalResearch,2006,25(3):388-396.[何剑锋,庄大方.长江三角洲地区城镇时
空动态格局及其环境效应.地理研究,2006,25(3):388-396.]
[17]MaXiaodong,MaRonghua,XuJiangang.SpatialstructureofcitiesandtownswithESDA-GISframework.Acta
GeographicaSinica,2004,59(6):1048-1057.[马晓冬,马荣华,徐建刚.基于ESDA-GIS的城镇群体空间结构.地理学报,2004,59(6):1048-1057.]
[18]MaRonghua,PuYingxia,MaXiaodong.MiningSpatialAssociationPatternsfromGISDatabase.Beijing:Science
Press,2006.[马荣华,蒲英霞,马晓冬.GIS空间关联模式发现.北京:科学出版社,2007.]
[19]TannierC,PumainD.Fractalsinurbangeography:Ageneraloutlineandanempiricalexample.Cybergeo,2005,307:
22.
[20]FrankhauserP.ComparingthemorphologyofurbanpatternsinEurope:Afractalapproach.In:EuropeanCities
InsightsonOutskirts,ReportCOSTAction10UrbanCivilEngineering,Vol.2,Structures,BorsdorfA,ZembriP(eds.),Brussels,2004:79-105
[21]ChenYanguang,LiuJisheng.Fractaldimensionsofspatialstructureofanurbansystemsandthemethodsoftheir
determination.GeographicalResearch,1999,18(2):171-177.[陈彦光,刘继生.城镇体系空间结构的分形维数及其测算方法.地理研究,1999,18(2):171-177.]
[22]LiuJisheng,ChenYanguang.Fractaldimensionsofscalestructureofanurbansystemandthemethodsoftheir
determination.GeographicalResearch,1998,17(1):82-89.[刘继生,陈彦光.城镇体系空间结构的分形维数及其测算方法.地理研究,1998,17(1):82-89.]
[23]ZhangLiang.Theregionaldevelopmentpattern:Transitionfromcenterspreadingpatterntolinkingconformitypattern.
华东理工大学学报(社会科学版),2005,(1):78-84.[张良.区域发展模式:从中心扩散型到联动整合型的转变.华东理工大学学报(社会科学版),2005,(1):78-84.]
[24]FishmanR.America'sNewCity:MegalopolisUnbound.WilsonQuarterly.1990,14:25-45.
[25]JohnsonGD,TempelmanA,PatilGP.Fractalbasedmethodsinecology:Areviewforanalysisatmultiplespatial
scales.Coenoses,1995,10:123-131
[26]LiX,YehAG.AnalyzingspatialrestructuringoflandusepatternsinafastgrowingregionusingRemoteSensingand
GIS.LandscapeandUrbanPlanning,2004,69:335-354.
[27]CliffAD,OrdJK.SpatialProcesses,ModelsandApplications.London:Pion,1981.
[28]AnselinL.Localindicatorsofspatialassociation:LISA.GeographicalAnalysis,1995,27:93-115.
[29]DaiQin,ChenXue,MaJianwenetal.Definitionsofnatural-languagespatialrelationsinGIS.Geomaticsand
InformationScienceofWuhanUniversity,2005,30(6):525-538.[戴芹,陈雪,马建文等.基于遥感数据SOFM网络分类的五种城市增长方式鉴别方法应用研究.武汉大学学报(信息科学版),2005,30(6):525-538.]
[30]ShenJ.ThepoliticaleconomyofdualurbanizationinChina.The2ndEastAsianRegionalConferenceinAlternative
Geography.HongKong,2000.
[31]WongK,ShenJ,FengZetal.Ananalysisofdual-trackurbanizationinthePearlRiverDeltasince1980.Tijdschrift
1022地理学报62卷
VoorEconomischeenSocialeGeografie,2003,94:205-218.
[32]LiXiaojian.ForeigndirectinvestmentanditsimpactoneconomicdevelopmentofcoastalChina.ActaGeographica
Sinica,1999,54(5):420-430.[李小建.外商直接投资对中国沿海地区经济发展的影响.地理学报,1999,54(5):420-430.]
[33]WuLiya,GuChaolin.Globalization,foreigninvestmentandurbanizationindevelopingcountries:Acasestudyof
Jiangsu.CityPlanningReview,2005,29(7):28-33.[吴莉娅,顾朝林.全球化、外资与发展中国家城市化:江苏个案研究.城市规划,2005,29(7):28-33.]
UrbanSpatialSprawlPatternandMetricsinSouthofJiangsu
ProvincealongtheYangtzeRiver
MARonghua1,GUChaolin2,PUYingxia3,MAXiaodong4,ZHUChuangeng4
(1.NanjingInstituteofGeographyandLimnology,CAS,Nanjing210008,China;
2.DepartmentofCityPlanning,TsinghuaUniversity,Beijing100084,China;
3.SchoolofGeographicandOceanographicSciences,NanjingUniversity,Nanjing210093,China;
4.CollegeofUrbanandEnvironmentalScience,XuzhouNormalUniversity,Xuzhou221009,China)
Abstract:Satelliteimagesacquiredin1984,1991,2000and2005,respectively,wereusedtointerprettheboundariesofcitiesandtownsinthesouthofJiangsuprovincealongtheYangtzeRiveronthebasisofourboundarydefinition.Andthenregionalurbanspatialsprawlpatternwasuncoveredquantificationallybysomeindexes,suchasfractaldimensionsincludingradiusdimension,griddimension,correlationdimension,boundarydimension,compactnessindex,andspatialautocorrelationindexincludingglobalMoranIandlocalMoranI.Theresultshowsthat:(1)fractaldimensionmetricsandspatialautocorrelationindex,respectively,canquantificationallyrevealregionalurbanspatialdistributionfeaturesandspatialsprawlprocessandclusterpattern,andtheirfusecanmeasureregionalurbanspatialsprawlpatternanduncovertheprocessofhowtoformtheurbancluster;and(2)attheinitialstageofreformandopen-uppolicy,theurbansprawlinthesouthofJiangsuprovincealongtheYangtzeRiverwasspontaneouslypoint-patterncenteredaroundbigcities,andthenpoint-axis-patterndominatedbypolicy-drivenurbanspatialfuseandfill,atpresentcluster-stretchedpatternbasedonregionalindustryspatialorganizationdrivenbyeconomicalradiationandcontact,andabigurbangroupiscomingintobeing.Especiallysince2000,spatialstructurefeaturesofanurbangrouparemoreandmoreevidentwithtimeelapsing.Additionally,themainfeatureofurbanizationinthestudyareaisthatindustrializationpushesurbanizationbasedonthebackgroundofglobaleconomy.Inthestudyareathreecentriccitieshaveexistedsincethe1980s,however,theirspatialpolarizationshavealessandlesseffectonothers.
Keywords:urbanization;spatialsprawlpattern;fractaldimension;spatialautocorrelation;Jiangsuprovince