基于土地利用变化的贵安新区景观生态风险评价

基于土地利用变化的贵安新区景观生态风险评价【摘要】以贵州省贵安新区2000、2005、2010、2015年土地利用数据为基础，划分2kmX2km的格网作为评价单元，借助Fragstats4.2软件计算景观指数，并构建景观生态风险评价模型，采用ArcGIS空间分析法，分析贵安新区生态风险时空分异及其对土地利用变化的响应。结果表明：①土地利用数量和结构变化明显：耕地和林地面积减少，草地和水域面积增加，建设用地面积显著增加；耕地主要转变为林地和建设用地；除水域外，其他类型分布集中程度增强。②生态风险集聚趋势增强，热点分布在贵安新区的西部和东部边缘地区，冷点分布在红枫湖水域；不同生态风险等级面积变化显著,低、较低和中生态风险区面积增加，较高和高生态风险区面积减少，总体生态环境质量改善。③低、较低生态风险区重心向东部的贵阳方向移动，中、较高和高生态风险区重心向西部的安顺方向移动，呈相反的转移方向。④除蔡官镇、城关镇和羊昌乡外，其他乡镇生态风险对土地利用响应弹性降低，影响生态风险的因素呈现多样化的特征。关键词：土地利用；生态风险；重心转移；贵安新区LandscapeecologicalriskassessmentofGui'anNewAreabasedonlandusechange.Abstract:BasedonthelandusedataofGui'anNewAreainGuizhouProvincein2000,2005,2010and2015,theecologicalriskassessmentunitsweredividedona2km乂2kmgridscale.WiththesupportofFragstats4.2,alandscapeecologicalriskassessmentindexwasconstructedtoanalyzethetemporal-spatialdiversityofecologicalriskanditsresponsetothelandusechangeinGui'anNewAreabyArcGISspatialanalysismodel.Theresultsshowedthat:(1)Thelandusequantityandstructurechangedsignificantly:thereductionofcultivatedlandandforestland,theincreaseofgrasslandandwaterarea,andtheobviousincreaseofconstructionland.Thecultivatedlandwasmainlytransferredtoforestlandandconstructionland.Exceptforthewaterarea,theconcentrationofothertypesincreased.(2)Theecologicalriskagglomerationtrendwasenhanced.ThehotspotsweredistributedinthewestandeastedgeofGui'anNewArea,andthecoldspotsweredistributedinthewatersofHongfengLake.Theareaofdifferentecologicalriskgradeschangedsignificantly,andtheareaoflow,lowmediumandmediumecologicalriskareasincreased,theareaofmediumhighandhighecologicalriskareasdecreased,andtheoverallecologicalenvironmentqualityimproved.(3)ThecentersofgravityoflowandlowmediumecologicalriskareasmovedtowardGuiyangintheeast,whilethatofmedium,mediumhighandhighecologicalriskareasmovedtowardAnshuninthewest,showingtheoppositedirection.(4)InadditiontoCaiguanTown,ChengguanTownandYangchangTown,theecologicalriskofothertownswaslessresilienttolanduse,andthefactorsaffectingecologicalriskshoweddiversifiedcharacteristics.Keywords:landuse;ecologicalrisk;centerofgravitytransfer;Gui'anNewArea图1研究区概况Figure1OverviewofthestudyareaNboOOENnooozzNbo。图1研究区概况Figure1OverviewofthestudyareaNboOOENnooozzNbo。芯匚二I贵安新区乡镇级行政区贵安新区DEMiW-1647低：1018 ]106°40'EN0Z09ZN0Z09Z0前言区域生态风险评价是指在区域尺度上，生态系统受到人类活动和自然环境变化等外界因素干扰时,对生态环境造成负面效应可能性的评价［1］o与传统的生态风险评价不同，区域生态风险评价更重视景观空间结构对生态风险过程的影响［2］o国内外关于区域生态风险的评价，研究对象主要针对生态风险终点评价［3〕、生态风险受体〔4］和整个区域的综合生态风险评价［5］o研究区域方面，以流域为评价范围的研究较多［6'7］,随着城市化的推进，研究区域也扩展至行政区范围［8'9］o评价方法可归纳为两类：一是基于风险源汇的评价方法，从风险受体和风险源两个方面构建评价的指标模型［10］o二是基于景观格局的评价方法，采用土地利用数据和景观格局指标构建生态风险评价指标［11］o20世纪90年代提出的土地利用/土地覆被变化已成为目前国际上研究的热点之一［12］。随着经济的发展和人口的增长，城市化进程加快，土地利用强度增加，生态安全受到影响［13］o研究发现土地利用与生态环境问题之间存在相关性，对景观格局［14］、生态服务价值［15］和生态风险的时空演变产生重要影响［16］o土地利用的强度体现了人类对自然环境影响的大小，呈现出区域性和累积性，直接反映在生态系统的结构和组成上［17］o所以在缺乏生态资料时，采用相对容易获取的土地利用数据来研究区域生态风险具有重要的意义。贵安新区位于长江流域和珠江流域交界的黔中地区，境内碳酸盐岩发育，生态环境脆弱，石漠化潜在威胁大。自2000年国家实施"西部大开发”战略以来，黔中经济区得以迅速发展，尤其是2014年国务院将黔中核心地带的贵安新区批复为第八个国家级新区以后，贵安新区始终肩负着“西部重要经济增长极、内陆开放型经济新高地、生态文明示范区”三大建设使命，其土地利用格局发生了较大改变，随之生态环境质量受到影响。因此，利用2000年、2005年、2010年和2015年的土地利用数据，从土地利用变化速度、转移方向和空间分异特征等方面，系统分析贵安新区2000-2015年土地利用变化规律，从景观格局角度构建生态风险评价指标体系,采用ArcGIS空间分析法，分析生态风险的时空分异及其对土地利用的响应，以期为土地资源的优化配置、提高土地的生态安全、评估该区域生态文明建设成效提供科学依据。1材料与方法1.1研究区概况贵安新区地处贵州高原中部，位于贵州省贵阳、安顺两市之间，分布范围26°11'7”—26°36'56"N,105°56'20”—106°39'10”E,包括清镇市、花溪区、平坝区、西秀区的20个乡镇，总面积约1895km2（图1）o地处高原型亚热带季风湿润气候区，年均降雨量1113—1367mm,平均相对湿度79%,平均气

温12.8—16.2°C,土壤种类繁多，组合复杂。贵安新区属于喀斯特低丘缓坡地貌，地势西高东低，总体较为平缓，海拔在1018—1647m之间。境内喀斯特地貌广泛发育，地形条件复杂，高原山地、丘陵盆地交错分布，且境内丘陵、山地大部分属于石漠化敏感地区。地表河流有马场河、麻线河、羊昌河、乐平河、青岩河等，流域总面积占新区总面积的80%,森林覆盖率达42%。1.2数据来源贵安新区土地利用数据包含4期，其中2000、2005、2010年土地利用数据来源于贵州省环保厅的“贵州省生态环境十年（2000—2010年变化遥感调查与评估）''。2015年土地利用数据是以2015年精度为30m的Landsat8遥感影像为数据源，对图像进行融合、矫正、裁剪等预处理，通过人机交互式解译得来，并对其进行野外验证，解译精度在80%以上。根据《土地利用现状分类标准》（GB/T21010—2017）,借鉴相关学者的研究成果［18-20］,以及研究区的实际,将土地利用类型合并为5个一级地类，最终建立4期土地利用数据库（图2）。贵州省级、市级、县乡级行政区划的矢量数据根据行政区划图提取。DEM数据（分辨率为30m）来源于地理空间数据云。1.3研究方法1.3.1土地利用动态变化土地利用变化量和变化率是表征土地变化状况的动态指标。土地利用变化量表示研究时段内某种土地利用类型面积变化数量，土地利用变化率反映土地利用类型变化的剧烈程度，计算公式参考文献［21］。土地利用转移矩阵是基于多个时期土地利用数据，对各类土地类型间数量转移和方向定量识别的方法,表达式参考文献a。1.3.2生态风险小区划分将研究区划分为若干个格网，对研究区进行空间重采样［23］。已有研究显示，渔网一般采用平均斑块面积的2—5倍［24］。结合研究区实际，将研究区划分为2kmx2km的网格，得到609个生态风险小区，计算每个网格的生态风险指数，以此作为生态风险小区中心点的生态风险指数。1.3.3景观生态风险指数构建景观格局反映了人类活动对生态系统的影响程度，不同景观抗外界的干扰能力也不同，利用景观结构指数、景观脆弱度指数和景观组分的面积比重,从景观格局的角度构建生态风险计算模型［25-26］,计算公式如下：nAERiii=1Ak式中:ERIi为第，个样本单元内的景观生态风险指数;Ak,为第k个样本单元内景观类型i的面积;Ak为第k个样本单元的面积;Ri为景观损失度指数。（1）景观结构指数（月）：建设用地■林地■水域图22000—2015年贵安新区土地利用空间分布Figure2LandusespatialdistributionofGui'anNewAreafrom2000to2015表示不同景观所代表的生态系统，受外界干扰

的程度I27四，计算公式如下：Si=aCi+bNj+cDi式中：C为景观破碎度指数；N,为景观分离度指数;Di为景观优势度指数；a,b,c为对应的权重,a+b+c=1,根据已有研究［28-29］和研究区实际情况，a、b、建设用地■林地■水域图22000—2015年贵安新区土地利用空间分布Figure2LandusespatialdistributionofGui'anNewAreafrom2000to2015景观脆弱度指数(F,):表示不同景观对外界的抗干扰能力［27］，根据已有研究成果［24'26］，并结合研究区特点，将5种景观格局按照脆弱程度由高到低赋值：水域5、耕地4、草地3、林地2、建设用地1,通过归一化得出景观脆弱度指数。景观损失度指数(Ri):表示不同景观所代表的生态系统，受外界影响时自身的损失程度㈣，公式为：Ri=JSi*Fi式中：Si表示第i景观类型的景观结构指数;F,表示第i景观类型的景观脆弱度指数。1.3.4空间分析法为了分析生态风险的时空演变，以构建的生态风险指数为空间变量，采用全局空间自相关(Moran,sI指数)和热点分析(Getis一OrdG：')探讨生态风险的自相关和集聚程度，公式参考文献网。基于地统计学中的半变异函数模拟分析［24］,再利用ArcGIS地统计模块中的克里金插值法，生成生态风险图［31］。重心转移模型分析不同等级生态风险空间重心的转移情况，公式参考文献［32］。1.3.5响应弹性分析引用弹性系数反映生态风险对土地利用的响应程度，以贵安新区乡镇为研究单元，基于土地利用和生态风险指数构建弹性指数,公式为：ERCERIb-ERIa/ERIaLC式中：ERC为研究时段内生态风险对土地利用变化的响应弹性系数;ERIb、ERIa分别为研究初期和末期的生态风险指数；LC为研究时段内的综合土地利用动态度，公式参考文献［33］。2结果与分析土地利用时空演变土地利用面积变化2000—2015年贵安新区土地结构变化较为明显,4个时期的耕地面积呈持续减少趋势，林地面积先增加后减少，草地面积持续增加，水域面积先不变后增加，建设用地变化最明显，呈持续增加趋势。从变化率上看，2000—2015年，建设用地的变化率最大，林地、草地和水域变化较大，耕地变化不明显,土地类型后期变化幅度比前期大。土地利用转移分析通过叠加分析得出贵安新区近15年的土地利用转移矩阵，结果显示:2000—2005年，以耕地转为草地和林地为主，其他类型转移较少。2005—2010年，主要以耕地转为建设用地为主，其次是草地转林地，其他类型转移不明显。2010—2015年，由于新区的成立转移面积明显增加，耕地转为林地、建设用地和草地,林地转为耕地和建设用地,草地主要转为耕地，建设用地主要转为耕地，水域转为耕地。土地利用变化的空间分异特征年份 -Moran5sI系数草地建设用地耕地林地水域2000年0.0058-0.0042-0.0019-0.00030.00712005年0.0035-0.0010-0.0029-0.00060.00712010年0.00590.0016-0.0033-0.00030.00712015年0.32750.00590.03860.02990.0004表12000—2015年贵安新区各土地利用类型Moran's/系数Table1Moran'sIcoefficientoflandusetypesinGui'anNewAreafrom2000to20152000、2005、2010和2015年各土地利用类型的Moran,sI系数如表1所示。总体来看，草地、耕地和林地的Moran,sI指数变化幅度较大，其他土地类型的Moran,sI变化不显著。草地的Moran,sI指数呈波动上升趋势，年份 -Moran5sI系数草地建设用地耕地林地水域2000年0.0058-0.0042-0.0019-0.00030.00712005年0.0035-0.0010-0.0029-0.00060.00712010年0.00590.0016-0.0033-0.00030.00712015年0.32750.00590.03860.02990.0004表12000—2015年贵安新区各土地利用类型Moran's/系数Table1Moran'sIcoefficientoflandusetypesinGui'anNewAreafrom2000to2015点主要分布在贵安新区北部，特别是红枫湖水域冷点集聚最明显，是因为政府重视红枫湖水域的生态保护，通过《贵州省红枫湖百花湖水资源保护条例》,使得水域面积上升，人类活动对其影响较小，从而2.2景观生态风险评价形成生态风险低值集聚区。2.2.1生态风险空间聚集分析利用ArcGIS空间自相关(Moran'sI)的方法，分析生态风险的集聚特征。经过计算,2000—2015年4期生态风险指数全局自相关指数，分别为0.2245,0.2867,0.2736,0.7213,且P值均小于0.05。Moran,sI为正值，表明生态风险与土地利用变化在空间上呈正相关的特征，表现出明显的集聚区，即高生态风险区分布在高生态风险区周围，低生态风险区分布在低生态风险区周围。2000—2015年全局莫兰系数呈上升趋势，表明贵安新区整体生态风险集聚程度增强。局部空间自相关可以表达出空间中的高值和低值聚集，采用热点分析(Getis—OrdG*)的方法，识别贵安新区2000—2015年生态风险的冷热点。由图3可看出，2000—2005年热点减少，冷点增加；2005—2010年热点和冷点均有增加，且趋于集中；2010—2015年热点和冷点均有减少。2000—2015年，热点主要分布在贵安新区西部和东部的边缘地区，这是由于城市快速发展，经济发展迅速，交通便利，土2.2.2生态风险时空变化由于2010年贵安新区的成立，土地利用强度地利用强度大，从而出现生态风险高值集聚区；冷图3贵安新区生态风险冷热点空间分布Figure3SpatialdistributionofhotandcoldspotsforecologicalriskinGui'anNewArea生态风险指数均值由2000年的0.1014下降到2015年的0.0983,表明2000—2015年生态环境转好。对生态风险指数做正态分布检验，表明空间上存在自相关性，采用半方差函数对2000、2005、2010和2015年生态风险指数进行拟合，通过对比，发现球状模型拟合效果最好，再利用ArcGIS克里金插值生成生态风险分布图，基于自然断点法并结合实际,将插值后的结果分为5个等级：低生态风险区(<0.090),较低生态风险区(0.090—0.095),中生态风险区(0.095—0.100),较高生态风险区(0.100—0.105),高生态风险区(20.105)。由表2和图4可知，2000—2015年贵安新区生态风险空间分异与演变明显。总体来看，15年间，低生态风险区面积先增加后减少,2000—2010年面积持续增加，主要分布在贵安新区东北部的红枫湖水域，以及北部乐平乡和十字乡的边缘，受人类活动影响小，生态环境好;2010年东部的石板镇和党武镇也有低生态风险区存在；2015年低生态风险区集中于红枫湖水域和党武镇的东部,

增大，造成低生态风险区面积减少。较低生态风险区面积持续上升，分布在低生态风险区的周围，较低生态风险区主要是由中生态风险区和低生态风险区转化而来的，这说明在政府环境保护政策的影响下，生态环境在转好。中生态风险区面积持续增加,从2010年开始占主导，各时期面积占比分别为27.05%、33.09%、37.57%、39.41%。高生态风险区面积总体呈下降趋势，各时期面积占比分别为17.13%、8.87%、4.54%、6.86%;2000年和2005年高风险区分布在大西桥镇北部，高峰镇北部，红枫湖镇东北部，马场镇东南部，党武镇西南部以及刘官乡附近；2010年和2015年仅分布在大西桥镇，马场镇东南部，党武镇西南以及刘官乡周围，由分散转为集中。较高生态风险区面积总体呈减少趋势,分布在高生态风险区周围，由分散到集中，2015年主要分布在贵安新区的中西部地区。总体来看,2000—2015年，低、较低生态风险区，从北部边缘和红枫湖水域，扩展至东部边缘以及中部地区；较高、高生态风险区，从东中西部都有分布，到仅分布在西部和东南边缘。2015年与2000年相比，低、较低和中生态风险区面积均有增加，分别增加1.89%、14.23%、由于2010年贵安新区的成立，土地利用强度地利用强度大，从而出现生态风险高值集聚区；冷图3贵安新区生态风险冷热点空间分布Figure3SpatialdistributionofhotandcoldspotsforecologicalriskinGui'anNewArea生态风险等级-2000年2005年2010年2015年面积/km2占比/%面积/km2占比/%面积/km2占比/%面积/km2占比/%低68.343.61127.886.75295.9615.62104.165.50较低131.226.92186.099.82349.2918.43400.8621.15中512.6727.05627.1433.09711.9237.57746.8639.41较高858.6545.31786.1541.49452.1523.86513.4927.10高324.5417.13168.168.8786.104.54130.056.86表22000—2015年不同生态风险等级面积和占比Table2Areasandproportionsofdifferentecologicalrisklevelsfrom2000to2015■较低生心风险区 ■高生态风险区图4贵安新区2000—2015年生态风险空间分布Figure4SpatialdistributionofecologicalrisksinGui'anNewAreafrom2000to2015从生态风险重心的转移角度和距离方面定量分析（图5）,低生态风险重心，2000—2005年，向西南转移5.44km,2005—2010年，转移速度明显加快,向东南转移32.9km,2010—2015年转移速度略微下降，向西北转移23.01km,2000—2015年，重心由十

字乡的中部转移到了东南部，由此可见低生态风险区总体有向东南转移的趋势。较低生态风险重心,2000—2005年，向东南转移14.99km,2005—2010年，向西南转移21.79km,2010—2015年，转移速度减缓，向东北转移13.09km,2000—2015年，较低生态风险重心由夏云镇的西北部转到东南部，总体向东南方向转移。中生态风险重心,2000—2005生态风险等级-2000年2005年2010年2015年面积/km2占比/%面积/km2占比/%面积/km2占比/%面积/km2占比/%低68.343.61127.886.75295.9615.62104.165.50较低131.226.92186.099.82349.2918.43400.8621.15中512.6727.05627.1433.09711.9237.57746.8639.41较高858.6545.31786.1541.49452.1523.86513.4927.10高324.5417.13168.168.8786.104.54130.056.86表22000—2015年不同生态风险等级面积和占比Table2Areasandproportionsofdifferentecologicalrisklevelsfrom2000to2015■较低生心风险区 ■高生态风险区图4贵安新区2000—2015年生态风险空间分布Figure4SpatialdistributionofecologicalrisksinGui'anNewAreafrom2000to20152.3生态风险对土地利用的响应分析由表3可知，2000—2015年贵安新区各乡镇生态风险对土地利用变化的响应弹性呈现出空间分异特征。2000—2005年，各乡镇的弹性系数都为正，是正响应，即研究区内各乡镇土地利用变化会引起生态风险的增加；其中大西桥镇和夏云镇的弹性系数较大，说明较小的土地利用变化则会导致生态风险的显著增加，生态环境脆弱；七眼桥镇、蔡官镇、十字乡和旧州镇弹性系数较低，表明土地利用变化导致生态风险小幅增加，生态环境相对稳定。2005—2010年，各乡镇的弹性指数仍为正，是正响应，即各乡镇土地利用变化会导致生态风险的增加；其中城关镇和大西桥镇弹性系数较大，七眼桥镇、蔡官镇、十字乡和天龙镇弹性系数较小。2010—2015年,麦坪镇和石板镇弹性系数为负，是负响应，即土地利用变化有利于生态风险的降低；其他乡镇弹性系数为正，马场镇和湖潮乡弹性较大，蔡官镇、旧州镇和白云镇弹性系数较小。2000—2015年，蔡官镇、城关镇和羊昌镇弹性系数先升高后降低，表明前期土地利用是生态风险增加的主要驱动力，后期土地利用对生态风险的影响减弱；其他乡镇弹性指数呈现下降趋势，虽然建设用地扩张会提高生态风险,但生态风险对土地利用的响应弹性降低，表明影响生态风险的因素呈现多样化的特征。3讨论低生态风险重心较低生态风险重心中生态风险重心较高生态风险重心高生态风险重心转移路径贵安新区乡镇界限图例O研究表明，土地利用变化会引起生态环境的变低生态风险重心较低生态风险重心中生态风险重心较高生态风险重心高生态风险重心转移路径贵安新区乡镇界限图例O表32000—2015年各乡镇生态风险对土地利用变化的响应弹性Table3Elasticityofecologicalriskofvarioustownsinresponsetolandusechangefrom2000to2015区域-响应弹性系数2000—2005年2005—2010年2010—2015年红枫湖镇12.128.201.40七眼桥镇9.714.081.32蔡官镇4.474.810.31高峰镇17.1214.490.91马场镇46.9510.22.22十字乡8.160.000.97旧州镇7.948.050.42夏云镇40.548.621.24麦坪镇14.712.99-0.55城关镇16.5851.281.43天龙镇20.830.001.04湖潮乡28.1714.292.10白云镇27.7817.610.46大西桥镇78.9533.330.82