基于空间统计的珠三角地区土地利用生态风险时空特征.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-1 0-2 5 修回日期:2 0 2 2-1 1-1 0 资助项目:国家自然科学基金“珠三角城市群的城市空间集聚时空分析与预测及空间优化决策”(4 1 0 0 1 0 7 8)第一作者:王彬(1 9 9 8),女(壮族),广西南宁人,硕士生,研究方向为土地利用生态风险评价。E-m a i l:w a n g b i n b i n_7 1 11 6 3.c o m 通信作者:梅志雄(1 9 7 6),男,湖北黄梅人,博士,教授,主要从事空间分析与建模、土地利用模拟与优化研究。E-m a i l:z h i x i o n g m e i 7 6 1 2 6.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 2 9.王 彬,梅志雄,魏俊超,等.基于空间统计的珠三角地区土地利用生态风险时空特征J.水土保持研究,2 0 2 3,3 0(6):4 2 3-4 3 1,4 4 0.W a n gB i n,M e i Z h i x i o n g,W e i J u n c h a o,e t a l.S p a t i o t e m p o r a lC h a r a c t e r i s t i c s o fL a n d-U s eE c o l o g i c a lR i s k i n t h eP e a r lR i v e rD e l t aR e g i o nB a s e do nS p a t i a lS t a t i s t i c 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):4 2 3-4 3 1,4 4 0.基于空间统计的珠三角地区土地利用生态风险时空特征王 彬,梅志雄,魏俊超,马 杰(华南师范大学 地理科学学院,广州5 1 0 6 3 1)摘 要:目的 揭示珠三角地区土地利用生态风险的空间结构模式和时空变化规律,为促进土地资源可持续利用和保障区域生态安全提供参考依据。方法 基于景观指数构建土地利用生态风险评价模型,测算2 0 0 0年、2 0 1 0年、2 0 2 0年珠三角地区多种尺度下的土地利用生态风险,利用半变异函数确定最适宜的研究尺度,运用空间自相关和G I S空间分析方法分析了适宜尺度下土地利用生态风险的时空特征。结果(1)3k m样方尺度是最佳研究尺度,小流域尺度也较合适;两个尺度下土地利用生态风险空间变异均以结构性因素影响为主,即存在较强的空间自相关性。(2)研究区土地利用生态风险总体上具较强的正空间集聚特征,但集聚程度趋于减弱;大多单元局部空间集聚特征相对稳定,以高高集聚和低低集聚为主,低高和高低集聚区零星分布且变化不大。(3)区域土地利用生态风险总体呈略微下降趋势,以低风险和较低风险等级区为主且面积逐渐增加,中风险等级以上面积减小。结论 珠三角地区土地利用生态风险整体上呈现正空间集聚特征,且土地利用生态风险水平逐渐降低,研究结果可为区域国土空间优化利用和生态风险管控提供决策参考。关键词:土地利用生态风险;半变异函数;空间自相关;珠三角中图分类号:X 8 2 6 文献标识码:A 文章编号:1 0 0 5-3 4 0 9(2 0 2 3)0 6-0 4 2 3-0 9S p a t i o t e m p o r a lC h a r a c t e r i s t i c so fL a n d-U s eE c o l o g i c a lR i s k i nt h eP e a r lR i v e rD e l t aR e g i o nB a s e do nS p a t i a l S t a t i s t i c sW a n gB i n,M e iZ h i x i o n g,W e i J u n c h a o,M aJ i e(S c h o o l o fG e o g r a p h y,S o u t hC h i n aN o r m a lU n i v e r s i t y,G u a n g z h o u5 1 0 6 3 1,C h i n a)A b s t r a c t:O b j e c t i v eT h ea i mo ft h i ss t u d yi st oe x p l o r et h es p a t i a ls t r u c t u r ep a t t e r na n ds p a t i a l-t e m p o r a lc h a n g er u l eo f l a n du s ee c o l o g i c a l r i s k i nt h eP e a r lR i v e rD e l t a(P R D)r e g i o n,a n dt h e nt op r o v i d es c i e n t i f i cs u p p o r t f o r t h e s u s t a i n a b l eu t i l i z a t i o no f l a n d r e s o u r c e s a n de c o l o g i c a l s e c u r i t y.M e t h o d sAl a n d-u s e e c o l o g-i c a l r i s ka s s e s s m e n tm o d e lw a sc o n s t r u c t e dw i t hl a n d s c a p em e t r i c s,a n dl a n d-u s ee c o l o g i c a lr i s ki n d e x e sa tm u l t i p l es c a l e s i nt h eP e a r lR i v e rD e l t ai n2 0 0 0,2 0 1 0a n d2 0 2 0w e r em e a s u r e db yt h i sm o d e l.T h e n,t h es e m i-v a r i o g r a m,s p a t i a l a u t o c o r r e l a t i o na n dG I Ss p a t i a l a n a l y s i sm e t h o d sw e r eu s e d t od e t e r m i n ea p p r o p r i a t er e s e a r c hs c a l e sa n de x p l o r et h es p a t i o t e m p o r a lc h a r a c t e r i s t i c so fl a n d-u s ee c o l o g i c a lr i s ka tt h e s es c a l e s.R e s u l t s(1)T h e3k mq u a d r a t i cs c a l ew a s t h eb e s t r e s e a r c hs c a l e,a n dt h es m a l lw a t e r s h e ds c a l ew a sa l s oa p p r o p r i a t e.A t b o t hs c a l e s,t h e s p a t i a l v a r i a b i l i t yo f l a n d-u s e e c o l o g i c a l r i s k f r o m2 0 0 0 t o2 0 1 0w a sm a i n l y i n f l u-e n c e db ys t r u c t u r a l f a c t o r s,i.e.,t h e r ew a sas t r o n gs p a t i a l a u t o c o r r e l a t i o n.(2)T h e l a n d-u s ee c o l o g i c a l r i s ki nt h es t u d ya r e ag e n e r a l l yh a ds t r o n gp o s i t i v es p a t i a la g g l o m e r a t i o nc h a r a c t e r i s t i c s,b u tt h ed e g r e eo fa g g l o m e r a t i o nt r e n d e dt ow e a k e n.T h el o c a ls p a t i a la g g l o m e r a t i o nc h a r a c t e r i s t i c so fm o s tu n i t sw e r er e l a t i v e l ys t a b l e.T h ed i s t r i b u t i o no fh i g h-h i g ha g g l o m e r a t i o na n dl o w-l o wa g g l o m e r a t i o na r e a sw e r eo b v i o u s,a n dt h ed i s t r i b u t i o no f l o w-h i g ha g g l o m e r a t i o na n dh i g h-l o wa g g l o m e r a t i o na r e a sw e r es p o r a d i ca n dd i dn o tc h a n g em u c h.(3)T h ed e g r e eo f l a n d-u s ee c o l o g i c a l r i s kh a dg r a d u a l l yd e c r e a s e d.L o w-r i s kl e v e l a n dl o w e r-r i s kl e v e lw e r ed o m i n a n t r i s k s a n d t h e d i s t r i b u t i o na r e a sw e r e i n c r e a s i n g,w h i l e t h e a r e a s a b o v em e d i u m-r i s k l e v e lw e r e d e c r e a-s i n g.C o n c l u s i o nT h eo v e r a l l l a n du s ee c o l o g i c a lr i s ki nt h eP R Dr e g i o ns h o w e dp o s i t i v es p a t i a lc l u s t e r i n gc h a r a c t e r i s t i c s,a n d t h e l e v e l o f l a n du s ee c o l o g i c a l r i s kg r a d u a l l yd e c r e a s e d,a n d t h e r e s u l t so f t h i s s t u d yc a np r o v i d ed e c i s i o n-m a k i n gr e f e r e n c e f o r t h eo p t i m a lu s eo f r e g i o n a l l a n ds p a c ea n de c o l o g i c a l r i s kc o n t r o l.K e y w o r d s:l a n d-u s ee c o l o g i c a l r i s k;s e m i-v a r i o g r a m;s p a t i a l a u t o c o r r e l a t i o n;P e a r lR i v e rD e l t a 当前,我国的城市化进程加快推进,加速了对土地的索取,土地资源利用程度持续增加,对自然和生态环境造成了极大压力,人地矛盾日渐突出,其中以土地利用结构不合理、土地生态风险加剧等问题最为明显1-3。研究表明:土地利用是人类活动与自然环境交互过程中的关键部分,与很多生态系统问题密切相关,甚至对区域生态安全起到首要作用4-5。不同土地利用方式和过程对区域内生态系统产生累积性的影响,其影响直观地反映在生态系统的结构和组成变化当中,并涉及到许多生态过程,由此产生了生态风险问题6-7。在此背景下,为探究区域土地利用变化对生态环境引起的不利影响,寻找降低其影响的方法和途径,开展区域土地利用生态风险分析并提出相应的管控措施,有利于促进土地资源可持续利用和保障生态安全。土地利用生态风险是指土地利用及其变化造成的生态系统结构破坏及功能降低的可能性8,近年来受到诸多学者的关注,且日益注重表现其时空分异特性9-1 0。土地利用生态风险在不同尺度上表现出的空间格局差异较大,具有尺度效应1 1。但现有相关研究通常基于一种空间尺度,如以格网作为评价单元(样方尺度)1 2-1 3或以行政区(省、区县、乡镇)为评价单元1 4,或小流域单元尺度等1 5。这些研究一定程度上推动了土地利用生态风险研究的进展,但仍存在如下缺陷:仅限于单一尺度的土地利用生态风险时空变化分析,普遍缺乏多种空间尺度下的土地利用生态风险规律的对比分析而未真正揭示其尺度效应,且主要依赖经验及主观因素选取研究尺度,缺乏客观判断依据1 6。目前土地利用生态风险评价方法主要有基于风险源汇和基于景观格局两类4,8。基于风险源汇的方法需要识别风险受体和风险源,涉及变量多且变量间关系复杂,数据获取困难,且对生态风险的空间表达能力不足、未考虑空间异质性;基于景观格局的土地利用生态风险评价方法是以土地利用变化为诱因,直接从景观格局角度定量识别和评估生态本底安全风险4,1 7-1 8。这类方法所需数据易获取,便于通过空间采样将景观空间结构转化为空间化的生态风险变量,与G I S结合可深入表达生态风险的空间特征、空间异质性及其时序变化4,7-8,1 9-2 0,逐渐成为近年土地利用生态风险评价的主流方法。目前国内外相关研究热点区域主要集中在城市1 8、流域7,2 1、沿海地带2 2以及重点自然保护区2 3等,而针对珠三角地区开展的相关研究较少。仅见叶长盛等2 4评价了1 9 9 02 0 0 6年珠三角地区土地利用生态风险,但缺乏2 0 0 6年后的分析,且只基于5k m5k m格网尺度。为此,本文再次从土地利用视角,选取城镇化快速发展、土地利用结构剧烈变化导致生态风险较为突出的珠三角地区为研究区,以2 0 0 0年、2 0 1 0年、2 0 2 0年的土地利用数据为基础构建生态风险指数,通过半变异函数工具对比分析样方尺度、行政单元尺度和小流域尺度下土地利用生态风险的空间异质性,并借助空间自相关和G I S空间分析方法揭示其时空变化特征。研究将力图揭示不同尺度下区域土地利用生态风险的空间结构模式、时空规律及其差异表现,探测其尺度效应,弥补现有相关研究的不足,并识别和探寻出需要关注的重点区域和潜在生态风险区,进而根据生态风险等级制定不同土地规划管理措施,为区域土地资源合理利用和规划管理及促进区域生态安全提供决策参考。1 研究区及数据珠三角(2 1 4 0 2 4 3 0 N,1 1 1 2 0 1 1 5 3 0 E)位于广东省中南部,毗邻南海,与东南亚隔海相望,总面积约5.4 8万k m2,包括广州、深圳、佛山、惠州、东莞、珠海、中山、肇庆和江门9个地级市。改革开放尤其是新世纪以来,伴随着经济和城镇化的快速发展,该区大量农业用地和生态用地转为建设用地2 5、土地利用结构变化剧烈,严重扰动了生态环境,区域内生态风险问题备受关注2 6。因此,从土地利用视角分析新世纪以来该区土地利用变化所引起的生态风险时空特征,对生态风险防控具有重要理论及实践指导意义。本文所用的数据主要包括数字高程模型(D E M)和土地利用数据。从地理空间数据云平台(h t t p:w w w.g s c l o u d.c n/)获得2 0 0 9年3 0m分辨率的D E M数据,用424 水 土 保 持 研 究 第3 0卷于划分分水岭小流域单元。通过全球地理信息公共产品平台(h t t p:w w w.g l o b a l l a n d c o v e r.c o m/)获取2 0 0 0年、2 0 1 0年、2 0 2 0年3 0m分辨率的地表覆盖栅格数据,根据研究区实际和本文需要,将该地表覆盖数据重新归类得到1 0 0m分辨率下耕地、林地、草地、建设用地、水体和未利用地6种地类。2 研究方法2.1 评价单元划分选取格网单元、行政区划单元和小流域单元,分别对应样方、行政单元和小流域3种空间尺度,来实现多尺度土地利用生态风险值的空间表达。具体划分情况为:(1)根据景观生态学原理,通常将研究区内斑块平均面积的25倍来确定格网单元的划分尺度1 9。根据研究区土地景观斑块大小实际,选取0.5k m为步长,从0.5k m0.5k m到5k m5k m创建1 0种格网单元,分别得到2 1 60 6 0,5 39 4 0,2 39 9 8,1 34 8 8,86 2 9,60 2 5,45 4 2,33 8 3,30 6 1,21 6 0个格网评价单元;(2)以研究区下辖5 6 9个乡镇为行政区评价单元;(3)使用A r c G I S的水文分析模块,从D EM提取出34 1 5个小分水岭作为小流域评价单元。2.2 土地利用生态风险指数为综合表达不同景观类型在外界影响下产生的潜在危害性,参考前人研究6,8,2 7,选择土地利用类型面积比例、景观干扰度指数和景观脆弱度指数来构建土地利用生态风险指数,计算公式为:E R Ik=ni=1Ak iAk(EiFi)(1)式中:E R Ik为单元k的土地利用生态风险指数;Ak和Ak i分别为单元k的面积(k m2)和单元k内地类i的面积(k m2);Ei和Fi分别为i地类的景观干扰度指数和景观脆弱度指数。2.2.1 景观干扰度指数 景观干扰度指数Ei表示景观要素受到外部环境干扰的程度2 8-2 9,可表达为:Ei=a Ci+b Si+c Ui(2)式中:Ci,Si,Ui分别为i地类的景观破碎度指数、分离度指数和优势度指数;a,b,c为权重,参考相关研究2 8并结合专家咨询分别赋值为0.5,0.3,0.2。(1)景观破碎度指数(Ci)指某一地类在外界干扰下从均匀的整体状态转变为异质的破碎状态,其值越大,单元内部越不稳定,可由i地类的斑块数与单元面积的比值表示。(2)景观分离度指数(Si)指某地类的斑块单元在空间上的分离程度,其值越高,表示该地类的空间散布程度越大,容易受到外界干扰。其计算公式为:Si=DiPI(3)Di=12niA(4)Pi=AiA(5)式中:Si为i地类的景观分离度指数;Di为i地类的距离指标;Pi为i地类的面积比值;A和Ai分别为单元面积(k m2)和单元内i地类的面积(k m2);ni为单元内i地类的斑块数。(3)景观优势度指数(Ui)指某一地类对景观整体的影响程度,其值越大,表明该地类对景观格局的影响程度越大,可表达为:Ui=(1+nin+Pi)/3(6)式中:Ui为i地类的景观优势度;n为单元内总斑块数;ni,Pi含义同上。2.2.2 景观脆弱度指数 景观脆弱度指数Fi表示景观要素结构受到外部因素扰动时的易损性,一般通过专家咨询法并结合实际对不同地类赋值并归一化获得。参考相关研究2 0,2 8和专家意见,并进行归一化得到6种地类的脆弱度指数如下:耕地0.1 9 05,林地0.0 9 52,草地0.1 4 2 9,建设用地0.0 4 7 2,水体0.2 3 8 1,未利用地0.2 8 57。2.3 半变异函数半变异函数是地统计学中研究区域化变量空间异质性的重要方法,能够表述变量的结构性和随机性变化特征3 0。其计算公式为:(h)=12N(h)N(h)i=1z(xi)-z(xi+h)2(7)式中:(h)为半变异函数;h为样本点间隔距离即步长;N(h)为步长为h时的样本点对数量;z(xi)和z(xi+h)分别为位置xi和(xi+h)处的土地利用生态风险值。半变异函数有3个重要拟合参数:(1)块金值C0,表示区域化变量的不确定性;(2)基台值C+C0(包括块金值C0和结构方差C),为半变异函数保持平稳时对应的数值,描述变量整体变异情况;(3)变程A,表示变量保持空间自相关性的距离。常用的半变异函数模型有指数模型、球状模型和高斯模型等,可通过模型拟合后的决定系数R2和残差R S S大小对比来确定最优拟合模型3 0。2.4 空间自相关方法空间自相关用来探测空间单元与邻域单元在相同属性值上的关联程度,并揭示其空间异质和空间集聚,以检验区域化变量的空间结构关系,包括全局和524第6期 王彬等:基于空间统计的珠三角地区土地利用生态风险时空特征局部空间自相关,常用M o r a n sI指数表示3 1。全局M o r a n sI反映区域单元观测值在整个研究区的空间自相关特征3 2,其表达式为:I=NNi=1Nj=1Wi jNi=1Nj=1Wi j(Xi-X)(Xj-X)Ni=1(Xi-X)2(8)式中:Xi为i单元的观测值;Wi j为行标准化的空间权重矩阵;I值范围落在-1,1,在定义的显著性水平下,当I显著为正时,单元观测值呈空间集聚分布;当I显著为负时,单元观测值呈空间散布状态;当I值为0时,表明单元观测值为独立随机分布。全局空间自相关虽可探测区域中具有同一属性单元的总体集聚程度,但无法表示局域单元的空间集聚特征。而局部空间自相关不仅可探测空间异质,也可准确反映出集聚区的具体位置和范围3 1。局部M o r a n sI表达式如下:Ii=Xi-XS2Nj=1,jiWi j(Xj-X)(9)式中:S为各单元观测值的标准差。局部M o r a n sI计算结果可划分为高高集聚、低低集聚、高低集 聚和低高 集 聚4种 类 型。高高(低低)型即某单元与邻域单元间表现为正空间自相关,高(低)值在局域空间呈现集聚现象;高低(低高)型即某单元与邻域单元间表现为负空间自相关,高(低)值单元被邻域的低(高)值包围,是空间异质区域3 1。3 结果与分析3.1 多尺度比较及适宜评价尺度的确定根据公式(1)(6),借助F r a g s t a t s软件计算得到选择的不同空间尺度下2 0 0 0年、2 0 1 0年、2 0 2 0年各评价单元的土地利用生态风险值,将其作为样本点数据进行后续研究。以2 0 0 0年土地利用生态风险指数计算结果为例,使用G S+软件计算得到其半变异函数拟合结果(表1)并对比不同尺度的参数变化特征,以选取最适宜研究尺度。表1 2 0 0 0年不同尺度下土地利用生态风险指数半变异函数模型拟合结果T a b l e1 L a n du s ee c o l o g i c a l r i s k i n d e xs e m i-v a r i a n c e f u n c t i o nm o d e l f i t t i n gr e s u l t sa td i f f e r e n t s c a l e s i n2 0 0 0尺度最优拟合模型C0C+C0C/(C+C0)A/mR2R S S0.5k m指数模型3.1 71 0-61.9 01 0-50.8 3 32 0 7 0 00.4 3 11.9 61 0-81k m指数模型6.0 01 0-63.4 01 0-50.8 2 42 1 5 0 00.5 7 55.4 91 0-1 11.5k m指数模型1.5 01 0-59.4 01 0-50.8 4 02 2 5 0 00.4 8 15.2 31 0-1 02k m指数模型4.7 01 0-52.8 61 0-40.8 3 62 5 5 0 00.4 9 26.1 01 0-9样方尺度 2.5k m指数模型8.8 01 0-55.2 61 0-40.8 3 32 9 4 0 00.5 3 52.2 31 0-83k m指数模型1.9 21 0-41.1 21 0-30.8 2 93 2 4 0 00.8 2 99.7 81 0-83.5k m球状模型7.4 51 0-41.6 31 0-30.5 4 34 3 6 0 00.6 0 02.3 91 0-74k m球状模型1.3 51 0-32.9 01 0-30.5 3 44 2 6 0 00.5 7 98.0 41 0-74.5k m球状模型1.8 21 0-33.9 71 0-30.5 4 04 3 7 0 00.5 6 21.7 71 0-65k m球状模型2.4 91 0-35.7 01 0-30.5 6 34 7 4 0 00.6 2 33.5 01 0-6小流域尺度 指数模型2.2 81 0-31.4 21 0-20.8 3 92 4 9 0 00.6 1 71.2 61 0-7行政单元尺度高斯模型1.9 91 0-28.6 71 0-20.7 7 01 7 6 6 60.4 5 31.8 61 0-7 由表1可知,0.5k m样方尺度下,土地利用生态风险指数的C0和C+C0均最小,说明该尺度下土地利用生态风险的随机性较小;随着样方尺度评价单元从0.5k m到5k m变化,C0和C+C0逐渐增加,说明土地利用生态风险的随机性逐渐增大,人类活动或自然灾害等外部干扰因素即偶然因素造成的空间异质性逐渐增强。行政单元尺度下的C0和C+C0均最大,说明该尺度下土地利用生态风险的随机性特征比其他尺度都明显。0.53k m样方尺度下,其结构方差/基台值即C/(C+C0)呈小幅波动,但均维持在0.8 2以上。根据地统计学原理,这说明结构性因素引起的空间异质性在样本间总变异中起主导作用,也反映了土地利用生态风险的空间自相关性较强1 6。而3k m5k m样方尺度下,C/(C+C0)总体呈快速下降态势,说明结构性因素引起的空间异质性逐渐减弱,随机因素的影响程度加深,土地利用生态风险的随机性变异增强。小流域尺度的C/(C+C0)达到0.8 3 9,行政单元尺度的C/(C+C0)也较高。从不同尺度拟合结果的决定系数R2来看,0.52.5k m样方尺度下的R2总体上比3.55k m样方尺度下的R2小,而3k m样方尺度下的R2为0.8 2 9,为所有尺度中的最大值;小流域尺度下的R2为0.6 1 7,拟合效果尚可;行政单元尺度下的R2为0.4 5 3,拟合效果较差。因此,3k m样方尺度是研究区2 0 0 0年土地利用生态风险分析的最佳尺度。以同样的方法分别对2 0 1 0年、2 0 2 0年的土地利用生态风险进行适宜评价尺度分析,均得到3k m样方尺度为最佳评价尺度。但考虑到小流域尺度的半624 水 土 保 持 研 究 第3 0卷方差函数拟合效果也较好以及为了保持研究区地理表面特征,最终选择3k m样方尺度和小流域尺度作为本文的分析尺度。3.2 适宜尺度下土地利用生态风险的空间变异分析由表2可知,3k m样方尺度和小流域尺度下2 0 0 0年、2 0 1 0年、2 0 2 0年半变异函数最优拟合模型均为指数模型。这两个尺度下各年拟合结果的决定系数R2均较高,且残差R S S均较小,表明各年拟合效果均较好。观察各年最佳拟合参数中的块金值C0,不论是3k m样方尺度还是小流域尺度,均为2 0 1 0年C0值最小、2 0 2 0年C0值最大,表明2 0 2 0年土地利用生态风险的空间异质性受到随机因素的影响程度较强。两个尺度下基台值在2 0 0 02 0 2 0年均逐渐降低,表明土地利用生态风险的空间变异程度有所减弱,C+C0变化不大说明结构趋于稳定。两个尺度下结构方差和基台值的比值C/(C+C0)均为2 0 1 0年最大,且2 0 0 0年和2 0 1 0年的C/(C+C0)均达8 0%以上,表明研究区土地利用生态风险的空间异质性受随机部分因素影响较小,以结构性因素影响为主,也反映该区域土地利用生态风险的空间自相关性较强;2 0 1 02 0 2 0年,两个尺度的C/(C+C0)均明显下降,说明随机因素对土地利用生态风险的影响逐渐增强。两个尺度下2 0 0 0年变程A均最大,表明该年研究区土地利用生态风险的空间相关性的作用范围较大。表2 2 0 0 0-2 0 2 0年土地利用生态风险指数半变异函数模型拟合结果T a b l e2L a n du s e e c o l o g i c a l r i s k i n d e xs e m i-v a r i a n c e f u n c t i o nm o d e l f i t t i n gr e s u l t s f o r2 0 0 0-2 0 2 0尺度年份最优拟合模型C0C+C0C/(C+C0)A/mR2R S S3k m样方尺度2 0 0 0指数模型1.9 21 0-41.1 21 0-30.8 2 93 2 4 0 00.8 2 99.7 81 0-82 0 1 0指数模型1.6 81 0-41.0 21 0-30.8 3 52 8 5 0 00.7 0 63.5 11 0-82 0 2 0指数模型2.8 41 0-41.1 01 0-30.7 4 22 7 5 0 00.9 3 97.7 51 0-9小流域尺度2 0 0 0指数模型2.2 81 0-31.4 21 0-20.8 3 92 4 9 0 00.6 1 71.2 61 0-72 0 1 0指数模型2.0 81 0-31.3 91 0-20.8 5 02 0 1 0 00.8 7 65.8 41 0-82 0 2 0指数模型2.7 11 0-31.1 61 0-20.7 6 72 0 4 0 00.8 5 33.6 61 0-83.3 土地利用生态风险空间自相关分析3.3.1 全局空间自相关分析 根据公式(8)计算各年两个尺度下土地利用生态风险的全局M o r a n sI值和相关指标(表3)。由表3可知,两个尺度下各年的全局M o r a n sI值均显著为正,说明2 0 0 02 0 2 0年珠三角土地利用生态风险整体呈现正空间集聚模式,即高高值集聚、低低值集聚。3k m样方尺度下全局空间正相关性整体比小流域尺度下更强;相较于2 0 0 0年、2 0 1 0年两个尺度下的全局M o r a n sI值小幅上升,说明2 0 0 02 0 1 0年两个尺度下的这种正空间集聚程度均略微增强。但2 0 2 0年两个尺度下的全局M o r a n sI值又略有下降,说明2 0 1 02 0 2 0年两个尺度下的这种正空间集聚程度趋于减弱。表3 3k m样方尺度和小流域尺度下土地利用生态风险的全局M o r a n sI和相关指标T a b l e3 G l o b a lM o r a n sIa n dr e l a t e d i n d i c a t o r so f l a n du s e e c o l o g i c a l r i s ka t t h e3k mq u a d r a t i c s c a l ea n ds m a l lw a t e r s h e ds c a l e尺度年份M o r a n sIZ得分p值3k m样方尺度2 0 0 00.5 7 55 9.4 5 20.0 0 12 0 1 00.5 8 86 0.0 3 70.0 0 12 0 2 00.5 7 18 3.1 4 40.0 0 1小流域尺度2 0 0 00.3 7 53 9.2 1 70.0 0 12 0 1 00.3 8 83 6.9 0 80.0 0 12 0 2 00.3 6 73 4.9 4 80.0 0 13.3.2 局部空间自相关分析 利用公式(9)计算各年两个尺度下各单元土地利用生态风险值的局部M o r a n sI值,然后借助A r c G I S软件空间化得到局部空间集聚格局演化图(图1)。研究期内,两个尺度下研究区大部分评价单元土地利用生态风险的局部空间结构变化较小,高高集聚和低低集聚单元相对较多且表现相对突出。高高集聚区主要集中在较高风险和高风险区域,其主要地类为耕地和水域等,但其面积随着时间推移有减少趋势,说明土地利用生态风险有一定程度的减弱趋向;低低集聚区主要集中分布在研究区中部,其土地利用类型较单一、多为建设用地,同时该集聚区有向外扩散趋向;低高集聚区零星分布于高高集聚区周围,主要地类为耕地和草地等;高低集聚区零星分布于低低集聚区周围,主要地类为建设用地和耕地,建设用地的景观脆弱度水平较低、内部结构稳定,对减轻土地利用生态风险值有一定作用。低高和高低集聚区零星分布且无多大变化,说明研究期内区域土地利用生态风险的空间异质区域并未发生剧烈的变动,区域内部土地利用导致的生态风险不均衡性虽然总体表现不强烈,但依然存在。3.4 土地利用生态风险的时空变化分析前文已证明3k m样方尺度和小流域尺度下3个年份的土地利用生态风险指数均具有显著的空间自相关性,因此可利用克里金法进行空间化插值。由724第6期 王彬等:基于空间统计的珠三角地区土地利用生态风险时空特征于两个尺度下划分的评价单元大小不同,生态风险指数值计算结果相差较大,但两个尺度的计算结果是相对的而不影响生态风险等级对比分析。为方便对比各年土地利用生态风险空间分布状况,采用自然断裂点法将两个尺度下的插值结果均划分为5个等级,其中3k m样方尺度下插值结果划分为低风险区(E R I 0.0 6 5)、较低风险区(0.0 6 5E R I 0.0 8 6)、中风险区(0.0 8 6E R I 0.1 1 2)、较高风险区(0.1 1 2E R I 0.1 5 0)、高风险区(0.1 5 0E R I 0.2 3 0),小流域尺度下的等级划分为低风险区(E R I 0.0 0 62)、较低风险区(0.0 0 62E R I 0.0 0 83)、中风险区(0.0 0 83E R I 0.0 1 06)、较高风险区(0.0 1 06E R I 0.0 1 4)、高风险区(0.0 1 4E R I 0.0 2 4)。最终得到3个年份土地利用生态风险等级空间格局图(图2)。图1 3k m样方尺度和小流域尺度下土地利用生态风险的局部空间集聚格局演化F i g.1 E v o l u t i o no f l o c a l i z e ds p a t i a l a g g l o m e r a t i o np a t t e r n so f l a n du s ee c o l o g i c a l r i s k sa t t h e3k mq u a d r a t i c s c a l ea n ds m a l lw a t e r s h e ds c a l e 从图2可以看出,3k m样方尺度和小流域尺度下的土地利用生态风险值的空间分布状况比较相似,都随着时间的推移,总体风险程度有所降低。研究区大多都以低风险和较低风险等级为主,