近30年贵阳市土地利用冲突时空变化分析.pdf

第3 0卷第6期2 0 2 3年1 2月水土保持研究R e s e a r c ho fS o i l a n dW a t e rC o n s e r v a t i o nV o l.3 0,N o.6D e c.,2 0 2 3 收稿日期:2 0 2 2-0 9-1 2 修回日期:2 0 2 2-1 0-1 2 资助项目:国家自然科学基金委员会贵州省人民政府喀斯特科学研究中心项目(U 1 8 1 2 4 0 1);贵州省科学技术项目(2 0 1 71 1 3;贵州省科技支撑项目(2 0 1 72 8 5 5)第一作者:陈大蓉(1 9 9 6),女(布依族),贵州贵阳人,硕士研究生,研究方向为地理信息与遥感应用。E-m a i l:1 4 0 3 5 6 5 5 6 3q q.c o m 通信作者:周旭(1 9 8 1),男,四川古蔺人,博士,副教授,研究方向为遥感水文与流域管理。E-m a i l:z x z y 8 1 7 81 6 3.c o mh t t p:s t b c y j.p a p e r o n c e.o r gD O I:1 0.1 3 8 6 9/j.c n k i.r s w c.2 0 2 3.0 6.0 4 0.陈大蓉,周旭,胡锋,等.近3 0年贵阳市土地利用冲突时空变化分析J.水土保持研究,2 0 2 3,3 0(6):3 3 7-3 4 4,3 5 6.C h e nD a r o n g,Z h o uX u,H uF e n g,e t a l.A n a l y s i so f S p a t i a l a n dT e m p o a lC h a n g e s i nL a n dU s eC o n f l i c t s i nG u i y a n gC i t y i n t h eL a s t 3 0Y e a r sJ.R e s e a r c ho fS o i l a n dW a t e rC o n s e r v a t i o n,2 0 2 3,3 0(6):3 3 7-3 4 4,3 5 6.近3 0年贵阳市土地利用冲突时空变化分析陈大蓉1,周 旭1,胡 锋2,裴 宇1,胡玉雪1,罗为维1(1.贵州师范大学 地理与环境科学学院,贵阳5 5 0 0 2 5;2.贵州省生态气象和卫星遥感中心,贵阳5 5 0 0 0 2)摘 要:目的 科学识别和掌握土地利用冲突的时空演变特征是区域寻求可持续发展的基础。方法 基于相对生态风险评价概念模型和景观生态学理论构建土地利用空间冲突测度模型,定量分析了贵阳市1 9 9 02 0 2 0年的空间冲突水平及时空演变特征。结果1 9 9 02 0 2 0年贵阳市土地利用冲突整体处于可控级别且变化明显,稳定可控占比由2 4.7 6%上升到6 1.1 0%,严重失控占比由2.2 4%波动上升到2.4 3%,有增长趋势。贵阳市土地利用冲突呈“南高北低”的空间格局,稳定可控和基本可控向息烽县、开阳县、乌当区等山地丘陵聚集;基本失控和严重失控趋于向城市市区及其周边区域集中。贵阳市土地利用冲突具有明显的空间聚集效应,整体呈“南高北低”相对的分布特征;热点区由清镇市、修文县、花溪区向城市中心迁移,冷点区主要聚集在乌当区、开阳县。结论 贵阳市近3 0年土地利用冲突时空变化显著,基本可控向稳定可控转变,失控地区向城市中心聚集,未来要重视城市中心生态建设,协调生态与经济的平衡发展。关键词:土地利用冲突;冲突测度;时空变化;贵阳市中图分类号:X 3 7;T P 7 9 文献标识码:A 文章编号:1 0 0 5-3 4 0 9(2 0 2 3)0 6-0 3 3 7-0 8A n a l y s i so fS p a t i a l a n dT e m p o a lC h a n g e s i nL a n dU s eC o n f l i c t s i nG u i y a n gC i t y i nt h eL a s t 3 0Y e a r sC h e nD a r o n g1,Z h o uX u1,H uF e n g2,P e iY u1,H uY u x u e1,L u oW e i w e i1(1.S c h o o l o fG e o g r a p h ya n dE n v i r o n m e n t a lS c i e n c e s,G u i z h o uN o r m a lU n i v e r s i t y,G u i y a n g5 5 0 0 2 5,C h i n a;2.E c o l o g i c a lM e t e o r o l o g ya n dS a t e l l i t eR e m o t eS e n s i n gC e n t e ro fG u i z h o uP r o v i n c e,G u i y a n g5 5 0 0 0 2,C h i n a)A b s t r a c t:O b j e c t i v eT h es c i e n t i f i c i d e n t i f i c a t i o na n dm a s t e r yo f t h es p a t i a l a n dt e m p o r a l e v o l u t i o n a r yc h a r-a c t e r i s t i c so f l a n du s e c o n f l i c t s a r e t h eb a s i s f o r t h e r e g i o n t os e e ks u s t a i n a b l ed e v e l o p m e n t.M e t h o d sB a s e do nt h ec o n c e p tm o d e l o f r e l a t i v ee c o l o g i c a l r i s ka s s e s s m e n ta n dt h et h e o r yo f l a n d s c a p ee c o l o g y,am e a s u r e-m e n tm o d e l o f l a n du s es p a t i a l c o n f l i c tw a s c o n s t r u c t e dt oq u a n t i t a t i v e l ya n a l y z e t h es p a t i a l c o n f l i c t l e v e l a n ds p a t i o t e m p o r a l e v o l u t i o nc h a r a c t e r i s t i c so fG u i y a n gC i t yf r o m1 9 9 0t o2 0 2 0.R e s u l t sF r o m1 9 9 0t o2 0 2 0,t h eo v e r a l l l a n du s e c o n f l i c t i nG u i y a n gC i t yw a s a t a c o n t r o l l a b l e l e v e l a n dc h a n g e ds i g n i f i c a n t l y.T h ep r o p o r-t i o no f s t a b l ea n dc o n t r o l l a b l er o s ef r o m2 4.7 6%t o6 1.1 0%.a n dt h ep r o p o r t i o no fs e r i o u s l yo u to fc o n t r o lr o s e f r o m2.2 4%t o2.4 3%,s h o w i n gag r o w t ht r e n d.T h e l a n du s ec o n f l i c t i nG u i y a n gC i t yp r e s e n t e das p a-t i a lp a t t e r no fh i g h l e v e l i nt h e s o u t ha n d l o wl e v e l i nt h en o r t h,w h i c hw a s s t a b l e,c o n t r o l l a b l ea n db a s i c a l l yc o n t r o l l a b l ea n dc o n v e r g e s i nX i f e n gC o u n t y,K a i y a n gC o u n t y,Wu d a n gD i s t r i c t,a n do t h e rm o u n t a i n sh i l l.B a s i c a l l yo u to f c o n t r o l a n ds e r i o u s l yo u to f c o n t r o l t e n d e dt oc o n c e n t r a t e i nu r b a na r e a sa n dt h e i r s u r r o u n d-i n ga r e a s.T h e l a n du s ec o n f l i c t i nG u i y a n gC i t yh a do b v i o u s s p a t i a l a g g r e g a t i o ne f f e c t,w h i c hw a s c h a r a c t e r-i z e db yar e l a t i v e l ys o u t h-n o r t hd i s t r i b u t i o n.T h eh o ts p o t sm o v e df r o m Q i n g z h e nC i t y,X i u w e nC o u n t y,a n dH u a x iD i s t r i c tt ot h ec i t yc e n t e r.T h ec o l ds p o t sm a i n l yc o n c e n t r a t e di n Wu d a n gD i s t r i c ta n dK a i y a n gC o u n t y.C o n c l u s i o nT h es p a t i o t e m p o r a l c h a n g e so f l a n du s ec o n f l i c t sw e r es i g n i f i c a n t i nG u i y a n gC i t yo v e rt h ep a s t 3 0y e a r s,b a s i c l l yc o n t r o l l a b l ea r e ac h a n g e dt os t a b l ea n dc o n t r o l l a b l ea r e a,a n dt h eo u t-o f-c o n t r o la r e a sg a t h e r e d i nu r b a nc e n t e r.G u i y a n gC i t ys h o u l da t t a c hi m p o r t a n c et ot h ee c o l o g i c a lc o n s t r u c t i o no f t h ec i t yc e n t e ra n dr e c o n c i l e t h eb a l a n c e dd e v e l o p m e n to f e c o l o g ya n de c o n o m y i nt h e f u t u r e.K e y w o r d s:l a n du s ec o n f l i c t;c o n f l i c tm e a s u r e m e n t;s p a t i a l t e m p o r a l c h a r a c t e r i s t i c s;G u i y a n gC i t y 土地是人类活动的基本场所,是社会经济发展和自然生态环境的重要资源1。土地资源的多宜性和土地供给的有限性往往使其在社会经济和生态环境之间存在博弈,衍生出许多土地利用之间的矛盾和冲突2。随着城市化和工业化进程加快,土地资源的多样化需求日益增加,不同利益相关者基于目标需求不同,造成土地利用冲突,使人与地之间的关系紧张,对区域稳定性和可持续发展构成了严重威胁。土地利用冲突的实质是不同利益相关者在利用土地资源的过程中因需求存在差异,在如何利用土地、所需土地数量等方面不一致、不和谐,以及不同土地利用方式与环境之间存在的矛盾状态2。因此,科学识别土地利用冲突类型,揭示冲突发生的原因和内在机制,探讨土地资源可持续管理方式是解决土地利用冲突问题的关键,是实现人与自然和谐共生的必要条件。2 0世纪7 0年代以来,土地利用冲突成为全球土地利用 及 相 关 学 科 研 究 关 注 的 热 点 问 题 之 一3。2 0 0 1年“自然资源管理和利用中的冲突管理办法”会议的召开,让土地利用冲突问题逐渐受到国内学者的关注。至今,土地利用冲突经历了从社会学、经济学、管理学到地理学的发展历程4-5,其研究也由冲突的概念、识别方法向时空演变、驱动因素和优化调控逐步扩展5-7。定性和定量识别土地利用冲突的方法很多,主要包括参与式方法8、博弈论9、压力状态响应概念模型1 0、多目标评价1 1、空间叠加法1 2和景观生态风险评价1 3等。不同的研究方法所适应的研究场景不同,参与式方法能深入了解冲突产生的原因,但易受被调查者主观影响,且耗时耗力。博弈论能探索土地利用冲突的内在机理,但不能识别冲突类型和空间分布。随着遥感技术的引入,实现了土地利用冲突的快速精准监测1 4。景观生态风险评价基于遥感监测的土地利用数据,借鉴生态风险评价概念模型和景观格局指数构建土地利用冲突模型。这种方法可以准确识别土地利用冲突的位置,揭示土地利用冲突的空间异质性,追踪区域时空变化特征。目前,从景观生态风险评价开展土地利用冲突的研究多从快速城市化出发,包括中国1 5、京津冀地区1 6、长株潭城市群1 7-1 8、黑河中游等1 3,1 9中宏观尺度,以及江苏省里下河2 0、重庆市江津区2 1等局地区域,对特殊地貌的山地城市关注较少。山地城市土地资源有限,人口密度增加和社会经济的不发达,造成土地资源的不可持续发展,导致生态失衡,景观连通性下降,土地利用冲突问题频频发生,影响区域生态安全。如何实现社会经济和生态环境和谐共存,筑牢生态安全屏障,是现阶段山地城市面临问题与挑战。经济和人口的快速增长加快了贵阳市的城市化步伐,土地利用空间冲突显著,亟需研究冲突的原因和空间分异特征。因此,本文以贵阳市为研究对象,采用景观生态学中的生态风险评价模型识别快速城市化和生态建设双重背景下的土地利用冲突及其时空变化特征,旨在协调山地城市的生产发展、生活宜居和生态安全格局,为山地城市绿色可持续发展提供科学依据。1 研究区概况贵阳市地处云贵高原黔中山原丘陵中部,地理坐标为东经1 0 6 0 7 1 0 7 1 7,北纬2 6 1 1 2 6 5 5,是贵州省政治、经济、文化、科教、交通中心(图1)。贵阳市国土总面积约为80 4 3k m2,下辖六区(云岩区、南明区、观山湖区、花溪区、白云区、乌当区)三县(修文县、息烽县、开阳县)和一个县级市(清镇市)。贵阳市位于长江流域和珠江流域的分水岭,属亚热带季风性湿润气候,气候温和,降水丰沛,雨热同季,年平均气温约1 5.3,多年平均降水量约10 9 6mm。地形以山地、丘陵为主,山多坡陡、景观异质性强,山地面积占5 2.3 0%2 2-2 3。区内喀斯特面积占8 5.0 2%,生态环境脆弱,破坏后难以恢复2 4。近3 0年来贵阳市发展迅速,据 贵阳市统计年鉴 和 贵州省统计年鉴统计,1 9 9 02 0 2 0年,人口由2 8 5.1 5万人增加到5 9 8.9 8万人,G D P由6 0.2 2亿元增长至43 1 1.6 5亿元,翻了近7 0倍;城市化率由4 3.9 1%增加到8 0.0 7%,增长率达8 2.3 5%。随着城市空间发展需要,建成区面积也由8 0.8 5k m2扩张到3 6 9k m2。同时,贵阳生态环境持续改善,森林覆盖率由1 4.5%增长到5 5%。833 水 土 保 持 研 究 第3 0卷图1 研究区示意图(审图号:G S(2 0 1 9)1 8 2 2号)F i g.1 S k e t c hm a po f t h e r e s e a r c ha r e a2 研究数据与方法2.1 数据获取与处理研究使用1 9 9 0年、2 0 0 0年、2 0 1 0年、2 0 2 0年四期土地利用数据是基于L a n d s a t卫星在G E E平台衍生的中国年度土地覆盖产品C L U D s(h t t p s:d o i.o r g/1 0.5 2 8 1/z e n o d o.4 4 1 7 8 1 0)2 5。以G o o g l eE a r t h遥感影像为底图,对其目视修正,经检验K a p p a系数均在8 0%以上。依据C L U D s分类体系、土地利用现状分类(G B/T 2 0 1 0-2 0 1 7)将土地利用类型分为林地、草地、耕地、水域、建设用地和未利用地6类。其中,采用A r c G I S将上述数据重采样为3 0m3 0m的栅格单元;投影坐标系均为WG S_1 9 8 4_U TM_Z o n e_4 8N。2.2 土地利用冲突指数构建基于生态安全视角,土地利用空间冲突指数可借鉴“风险源风险受体风险效应”的相对生态风险评价概念模型和景观生态学理论进行测度1 3。风险源是指可能对区域生态系统造成不利作用的风险来源2 6,选用土地利用强度相关的外部压力因子作为指标;风险受体指生态系统中受到风险源不利作用的载体2 7,选用空间脆弱性作为指标;风险效应指不同风险源作用于不同风险受体后的反应2 7,选用空间稳定性指数作为指标。空间单元的外部压力越大,风险暴露程度越高,内部稳定性越低,生态风险的可能性越大,也就意味着生态系统会遭受更大的干扰,加深各类土地利用间的矛盾,土地利用冲突问题亦更加突出。同理,空间格局变化对区域生态系统功能的影响越小,带来的空间生态风险越小,则土地利用冲突程度越低2 8。参考以往研究1 8,土地利用冲突指数L U C I(L a n dU s eC o n f l i c t I n d e x)可以表示为:L U C I=P I+-S I(1)式中:L U C I为土地利用冲突指数;P I为外部压力指数;为脆弱性指数;S I为稳定性指数。(1)空间外部压力指数(P I)。城市化不断加快及人类活动的过度干预使得区域内的土地利用变得更加复杂和破碎,降低土地利用效率的同时加剧了土地利用冲突。面积加权平均斑块分维数(A WM P F D)基于周长面积的关系来度量景观斑块的复杂性,即表示邻域景观对当前土地利用的干扰影响程度2 9。因此,借鉴A WM P F D度量空间外部压力指数,分维数越大则表示景观斑块的形状越复杂,不易受人类活动影响;反之,分维数越小则容易受人类活动影响1 8。其计算公式为:AWMP F D=mi=1nj=12 l n(0.2 5Pi jl n(ai j)(ai jA)(2)式中:Pi j为斑块周长;ai j为斑块面积;A为景观总面积;m为空间单元总数;n为土地利用类型。(2)空间脆弱性指数()。空间脆弱性指数是用来表征不同风险受体(不同土地利用类型)抵抗外界干扰能力的大小,其值越高,生态风险越大2 9。不同土地利用类型对外部压力的抗干扰能力不同,借鉴景观脆弱性指数(F I)来度量空间脆弱性,参考已有文献3 0,确定建设用地、其他用地、水域、耕地、草地、林地的脆弱度指数分别为5,4,3,2,1,1。其计算公式为:=ni=1FiaiS(3)式中:Fi为景观i的脆弱性指数;ai为单元内景观i的面积;S为空间单元总面积。(3)空间稳定性指数(S I)。景观破碎化是区域空间形态中土地利用冲突的主要后果之一,景观斑块越破碎,稳定性越差、空间冲突越强烈,同时越容易诱发生物多样性丧失,导致生态系统失衡3 1。因此,借鉴景观生态学中的斑块密度(P D)从负面衡量区域景观的稳定性7。S I=1-P D=1-niA(4)式中:ni为单元内景观i的斑块数目;A为空间单元面积。2.3 土地利用冲突分析尺度景观 格 局 指 数 对 景 观 尺 度 具 有 明 显 的 依 赖性1 3,贵阳市南北最大纵距和东西最大横距分别为1 3 0k m,1 1 5k m,参考彭佳捷等1 8选择空间单元尺度的方法,文中选择24 0 0m格网尺度为评价单元。其中,研究区边界没有覆盖整个单元面积的斑块按一个完整单元进行计算。为方便计算,将公式(2)(4)的外部压力指数、脆弱性指数和稳定性指数分别标准933第6期 陈大蓉等:近3 0年贵阳市土地利用冲突时空变化分析化处理至(0,1),代入公式(1)中计算各空间单元的土地利用冲突指数并标准化至(0,1)。通过测算4个时期各空间单元网格的土地利用冲突强度,参考吴蒙等2 8提出的空间冲突分级方法将其划分为4个等级,其中0,0.5 为稳定可控、(0.5,0.7 为基本可控、(0.7,0.8 为基本失控、(0.8,1 为严重失控。3 结果与分析3.1 土地利用冲突的时间演变特征从土地利用冲突指数结果(表1)可以看出:贵阳市1 9 9 0年、2 0 0 0年、2 0 1 0年、2 0 2 0年土地利用冲突的平均值分别为0.5 7 3,0.5 6 7,0.4 7 9,0.4 8 5,即贵阳市在研究期间以稳定可控和基本可控为主,土地利用冲突呈现出“先降低后增加”的变化趋势。其中,稳定可控的 空 间 单 元 数 量 占 比 在1 9 9 02 0 0 0年 稳 定 在2 4%,到2 0 1 0年和2 0 2 0年空间单元数量占比增加了近1.5倍,占贵阳市空间单元数量的3/5左右;基本可控级别空间单元数量所占比例在研究期间内呈先增后减趋势,在2 0 0 0年达到最大值,但增加的比例远远小于减少的比例,对调控空间冲突生态风险及预防区域生态安全问题发生发挥积极作用。与之相反,基本失控和严重失控的空间单元数量所占比例呈先减少后增加趋势。1 9 9 02 0 1 0年基本失控的空间单元占比大幅度减少,由1 1.0 7%降到2.1 1%,2 0 2 0年又上升至3.7 1%,基本失控是土地利用冲突的初始表现,对区域生态安全存在潜在威胁,调控不当容易加剧冲突;严重失控空间单元数量占比变化趋势与基本失控一致,2 0 1 02 0 2 0年从0.7 7%增加到2.4 3,增幅1.6 6%,有逐渐加重的发展趋势,是今后要重点调控土地利用方式的空间单元,以降低生态安全风险。表1 贵阳市1 9 9 0-2 0 2 0年空间冲突测算统计结果T a b l e1 T h e s t a t i s t i c a l r e s u l t so f s p a t i a l c o n f l i c tm e a s u r e m e n t i nG u i y a n gC i t y f r o m1 9 9 0t o2 0 2 0冲突级别冲突值空间单元个数/个1 9 9 0年2 0 0 0年2 0 1 0年2 0 2 0年空间单元百分比/%1 9 9 0年2 0 0 0年2 0 1 0年2 0 2 0年稳定可控0.00.5 03 8 73 7 49 1 99 5 52 4.7 62 3.9 35 8.8 06 1.1 0基本可控0.5 10.7 09 4 01 0 2 15 9 95 1 26 0.1 46 5.3 23 8.3 23 2.7 6基本失控0.7 10.8 01 7 31 3 33 35 81 1.0 78.5 12.1 13.7 1严重失控0.8 11.06 33 51 23 84.0 32.2 40.7 72.4 3总计1 5 6 31 5 6 31 5 6 31 5 6 31 0 01 0 01 0 01 0 0冲突平均值0.5 7 30.5 6 70.4 7 90.4 8 53.2 土地利用冲突的空间演变特征从贵阳市1 9 9 02 0 2 0年3 0年间土地利用冲突空间分布(图2)可以看出:空间上具有显著的南北差异,即“南高北低”。土地利用冲突的高值区主要分布于城市中心及周边区域,这些区域地势相对较低,人口和资本相对集中,城市建成区不断向外延伸,侵占周边的林地、耕地等,造成土地利用冲突水平提高;土地利用冲突低值区主要分布于息烽县、开阳县、乌当区、修文县东部和花溪区中部等山地丘陵地区,这些区域林地资源丰富且连片、植被覆盖度高,受人类活动扰动小,因此土地利用冲突水平较低。结合土地利用冲突空间变化趋势和面积转移特征(图3)来看,1 9 9 02 0 0 0年贵阳市土地利用冲突失控级别区域主要集中在农村,各冲突级别转移情况较少,以维持自身冲突状态为主。此时人口快速增加,且聚集于农村,人们为满足生活需要盲目开荒毁林,耕地大面积扩张的同时加剧土地的破碎化,导致土地冲突水平较高。2 0 0 02 0 1 0年贵阳市土地利用冲突问题得到明显的控制,处于研究期间的最低水平,失控级别区域主要集中在中心城区及周边区域。各冲突级别变化显著,基本可控 转 向 稳 定 可 控 的 面 积 有31 5 0.7 2k m2,占基本可控单元的5 3.5 7%,转向基本失控和严重失控的较少,分别占0.7 8%和0.3 9%;基本失控分别转向稳定可控和基本可控5 7.6k m2,6 3 9.3 6k m2,二者占比达到9 0.9 8%。贵阳市作为西部重要的经济发展区,在西部大开发战略的实施中,对土地利用需求日益增多,建成区向白云区、观山湖区和花溪区方向蔓延,但土地利用冲突并未呈现增加态势,这主要归功于贵阳市在发展经济的时候注重生态文明建设,同步实施了石漠化治理、退耕还林等生态工程,有目标的提升区域生态环境质量和生态韧性。2 0 1 02 0 2 0年贵阳市土地利用冲突开始出现“反弹”,中心城区土地利用冲突问题上升,与其他各区县冲突问题对比显著。失控级别变化比较明显,严重失控有5.7 6k m2的面积转向基本失控,但远不及基本失控向严重失控转移面积9 7.9 2k m2,且可控级别也有少量空间单元转为失控级别。2 0 1 0年以后,一系列生态工程虽有效缓解土地利用冲突,但贵阳市城市化水平仍不断提升,建设用地持续扩张,侵占周边土地,以及大量农村自建房占用耕地,进一步加剧人地矛盾,进而加重土地利用冲突。043 水 土 保 持 研 究 第3 0卷图2 贵阳市1 9 9 0-2 0 2 0年空间冲突指数变化图F i g.2 T h e s p a t i a l c o n f l i c t i n d e xc h a n g e s i nG u i y a n gC i t y f r o m1 9 9 0t o2 0 2 0图3 贵阳市1 9 9 0-2 0 2 0年不同空间冲突面积转移弦图F i g.3 T h ea r e a t r a n s f e rc h o r dd i a g r a mo fd i f f e r e n t s p a t i a l c o n f l i c t s i nG u i y a n gC i t y f r o m1 9 9 0t o2 0 2 03.3 空间自相关分析采用A r c G I S1 0.2中空间统计的M o r a n sI对贵阳市1 9 9 02 0 2 0年土地利用冲突综合指数进行全局空间自相关分析,结果显示:各期的全局自相关M o r a n sI指数分别为0.5 1 61,0.4 9 36,0.5 6 04,0.6 5 24,均为正值,Z值分别为2 8.3 1,2 7.0 8,3 0.7 5,3 5.7 9,p值均小于0.0 0 1,分析表明贵阳市不同空间单元的空间冲突水平在空间格局呈正相关,存在明显的聚集趋势,并非随机分布。随着时间的推移,M o r a n sI指数先降后升,Z值亦是如此,呈现“V”型变化趋势,在2 0 0 0年达到最低值,表明2 0 0 0年空间冲突聚集性下降。2 0 0 0年以后宏观政策对土地利用变化起导向作用,使冲突的空间分布具有明显的聚集趋势。为进一步探究不同冲突级别的集聚和离散规律,采用A n s e l i nL o c a lM o r a n sI工具得到1 9 9 02 0 2 0年局部空间自相关L I S A聚集图(图4)。高高类型在1 9 9 02 0 2 0年呈现零散逐渐集中的变化趋势,1 9 9 0年高高类型聚集状态较分散,主要分布在清镇市东部和西部、花溪区南部和北部,云岩区、南明区和修文县零星分布;2 0 0 0年高高类型小幅度减少,主要分布在清镇市西部,花溪区南部由高高类型转为不显示类型;2 0 1 0年的高高类型在中部出现明显的聚集现象,且清镇市西部分布范围再次减少,修文县亦减少;2 0 2 0年高高类型在中部范围向花溪区、白云区和观山湖区延伸,小部分零星分布,说明贵阳市生态保护政策对土地利用冲突的正向影响显著。低低类型主要分布在贵阳市西北方向,即乌当区和开阳县,3 0年变化不大。这些区域受交通条件和地形影响,城市扩张受到限制,同时离城市中心较远,土地利用冲突发生可能性较低。3.4 冲突冷热点分析采用A r c G I S1 0.2中空间统计的关联指数G e t i s-O r dGi*,选取置信度分别为9 0%,9 5%和9 9%的冷热点来分析贵阳市近3 0年空间冲突的冷热点空间分布可知(图5):贵阳市空间冲突冷热点空间分布与局部空间自相关L I S A聚集分布特征较为一致,冷热点沿西北东南轴相对分布,轴线以上为冷点分布区,反之为热点分布区。1 9 9 02 0 2 0年3 0年的热点区域差异明显,呈现由四周零散分布向中心城区集中分布的变化趋势,冷点区域变化不明显。具体而言,1 9 9 0年贵阳市热点主要分143第6期 陈大蓉等:近3 0年贵阳市土地利用冲突时空变化分析布于中部和西部清镇市、西北部修文县、南部花溪区边缘地区零星分布;2 0 0 0年与1 9 9 0年较为接近,差异之处在于修文县西部出现明显热点区;2 0 1 0年清镇县和修文县热点区明显减少,中部显著增加,除息烽县外,各区县零星分布;2 0 2 0年热点区集中分布于中部平缓地,该区域根据地形特征将建成区布局沿中曹司盆地南北向拓展,以及设立观山湖新区以解决人口增长对住宅和城市设施的需求。冷点区域主要分布于乌当区、开阳县东部和息烽县、开阳县交界处,且逐渐向乌当区和开阳县方向偏移。图4 贵阳市1 9 9 0-2 0 2 0年空间冲突局部自相关L I S A聚集图F i g.4 T h e s p a t i a l c o n f l i c t l o c a l a u t o c o r r e l a t i o nL I S Aa g g r e g a t i o nm a pf r o mi nG u i y a n gC i t y1 9 9 0t o2 0 2 0图5 贵阳市1 9 9 0-2 0 2 0年空间冲突冷热点分布F i g.5 T h ed i s t r i b u t i o no f s p a t i a l c o n f l i c t c o l dh o t s p o t s i nG u i y a n gC i t yf r o m1 9 9 0t o2 0 2 04 讨论与结论4.1 讨 论本文分析了贵阳市土地利用冲突的时空演变特征,研究表明,贵阳市基本失控和严重失控区域向基本可控和稳定可控转变,冲突区域高度集中在中心城区及周边区域,土地利用冲突先减弱后增加,有逐渐恶化趋势。土地利用冲突变化情况与曾昉3 2、李雪冬3 3等研究的生态冲突及生态安全变化情况基本一致。随着喀斯特山地城市经济的快速发展,生产生活空间持续扩张,加剧生态空间的“破碎化”,人地矛盾日益突出,严重影响城市的生态安全。从2 0 0 1年“环境立市”发展战略的确立,到2 0 1 4年提出建设全国生态文明示范城市,贵阳市一直奋力走出一条绿色可持续发展道路,以实现经济发展与生态保护的有利局面。在这些目标的驱动下,退耕还林、环城林带、石漠化治理和天然林林业生态建设,采石场、矿山等环境治理和复垦,城市绿化工程等多项措施的实施增加森林总量,森林覆盖率由1 9 9 0年的1 4.5%提升至2 0 2 0年的5 5%,建成区绿地率由2 9.4%增长至3 9.6%,生态环境进一步优化,土地利用冲突逐渐减弱,城市生态承载力增强。但2 0 2 0年中心城区及周边区域,基本失控和严重失控区域空间单元占比增加,空间分布特征与贵阳市“东扩、西连、南延、北拓”的城市扩展路径一致,土地利用冲突开始有恶化趋势。因此,未来的城市发展要秉承“两山”理论的发展理念推进城市生态文明建设,集约化工业用地和住宅以提高土地利用效率,坚持不以牺牲环境为代价发展经济,对已开发区域开展生态修复工程,改善土地利用冲突问题,提高区域的生态系统服务效益和生态环境承载力,筑牢城市生态安全。本文借鉴生态风险评价概念模型和景观生态学243 水 土 保 持 研 究 第3 0卷理论揭示研究区土地利用冲突的时空变化规律。然而,不同利益相关者开发利用土地资源的原因并不唯一,还涉及了社会、经济、生态等各方面。本文从土地利用类型角度进行空间冲突的测度,结果获取快速、直观、省时省力,但没有考虑社会、经济等因素,一定程度上限制了其适用性,其计算指标需进一步改进。4.2 结 论(1)1 9 9 02 0 2 0年贵阳市土地利用空间冲突均值由0.5 7 3波动下降到0.4 8 5,冲突整体处于可控水平且变化明显。稳定可控2 0 1 0年的占比是2 0 0 0年的1.5倍;严重失控空间单元由1 9 9 0年的6 3个下降到2 0 1 0年的1 2个,又上升到2 0 2 0年的3 8个,冲突有加剧趋势。(2)1 9 9 02 0 2 0年贵阳市土地利用冲突空间差异显著,呈“南高北低”的空间格局,稳定可控和基本可控向息烽县、开阳县、乌当区等山地丘陵聚集特征;基本失控和严重失控集中于城市市区及其向外蔓延区域。(3)1 9 9 02 0 2 0年贵阳市土地利用冲突的全局M o r a n sI指数均大于0,表明空间冲突具有显著的空间正相关性和聚集效应;不同阶段其空间聚集分布特征不同,但整体上与空间冲突“南高北低”相对的空间分布特征一致;热点区由清镇市、修文县、花溪区向城市中心迁移,冷点区主要聚集在乌当区、开阳县,且逐渐向城市边缘偏移。参考文献:1 刘纪远.中国资源环境遥感宏观调查与动态研究M.北京:中国科学技术出版社,1 9 9 6.L i uJY.S t u d yo nr e s o u r c e s&e n v i r o n m e n ts u r v e ya n dd y n a m i cu s i n gr e m o t es e n s i n gi nC h i n aM.B e i j i n g:C h i n aS c i e n c ea n dT e c h n o l