1980-2020年松嫩平原生态脆弱性时空分异.pdf

第4 3卷第4期2 0 2 3年8月水土保持通报B u l l e t i no fS o i l a n dW a t e rC o n s e r v a t i o nV o l.4 3,N o.4A u g.,2 0 2 3 收稿日期:2 0 2 2-1 0-0 8 修回日期:2 0 2 2-1 2-2 6 资助项目:吉林省自然科学基金项目“松嫩沙地灌丛沙丘发育的地貌生态过程”(2 0 2 1 0 1 0 1 3 9 8 J C)第一作者:宋春威(1 9 9 8),男(满族),吉林省伊通县人,硕士研究生,研究方向为生态环境遥感。E m a i l:s c h w 1 8 1 71 6 3.c o m。通讯作者:杜会石(1 9 8 3),男(汉族),吉林省伊通县人,博士,教授,硕士生导师,主要从事干旱区地貌方面的研究。E m a i l:d u h s j l n u.e d u.c n。1 9 8 02 0 2 0年松嫩平原生态脆弱性时空分异宋春威1,2,相恒星3,陈美宏1,2,杜会石1,2(1.吉林师范大学 旅游与地理科学学院,吉林 四平1 3 6 0 0 0;2.吉林省生态环境遥感大数据重点实验室,吉林 四平1 3 6 0 0 0;3.中国科学院 东北地理与农业生态研究所,吉林 长春1 3 0 1 0 2)摘 要:目的探究1 9 8 02 0 2 0年松嫩平原生态脆弱性时空变化特征和驱动机制,为该区生态系统管理与可持续发展提供科学依据。方法从自然和人文两个角度出发,采用“压力状态响应”(P S R)模型构建评价指标体系,结合主成分分析法和M o r a nI指数开展研究。结果松嫩平原生态脆弱性以极度、重度和中度脆弱性为主,共占8 1.7 2%。2 0 0 0年之前生态脆弱性整体水平呈现上升趋势,2 0 0 0年之后生态脆弱性整体水平不断下降。松嫩平原生态脆弱性以轻度和中度脆弱性变化为主,变幅分别为6 2.1 9%和9 1.3 0%。生态脆弱性变化是由于草地、水域、未利用地、建筑用地和耕地间的转化。生态脆弱性存在空间自相关性,低低集聚主要分布在研究区西部林地区域。结论近4 0a松嫩平原生态脆弱性总体呈现波动降低趋势,空间上生态脆弱性自东向西逐渐增加,气温、降水等自然因素和土地利用等人文因素对生态脆弱性影响较大。关键词:生态脆弱性;空间自相关;空间主成分分析;生态环境;松嫩平原文献标识码:A 文章编号:1 0 0 0-2 8 8 X(2 0 2 3)0 4-0 3 6 6-0 9 中图分类号:X 8 2 6,X 2 2文献参数:宋春威,相恒星,陈美宏,等.1 9 8 02 0 2 0年松嫩平原生态脆弱性时空分异J.水土保持通报,2 0 2 3,4 3(4):3 6 6-3 7 4.D O I:1 0.1 3 9 6 1/j.c n k i.s t b c t b.2 0 2 3.0 4.0 4 2;S o n gC h u n w e i,X i a n g H e n g x i n g,C h e nM e i h o n g,e ta l.S p a t i o t e m p o r a lv a r i a b i l i t yi ne c o l o g i c a lv u l n e r a b i l i t yo fS o n g n e nP l a i nf r o m1 9 8 0t o2 0 2 0J.B u l l e t i no 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,4 3(4):3 6 6-3 7 4.S p a t i o t e m p o r a lV a r i a b i l i t y i nE c o l o g i c a lV u l n e r a b i l i t yo fS o n g n e nP l a i nf r o m1 9 8 0t o2 0 2 0S o n gC h u n w e i1,2,X i a n gH e n g x i n g3,C h e nM e i h o n g1,2,D uH u i s h i1,2(1.C o l l e g eo fT o u r i s ma n dG e o g r a p h i c a lS c i e n c e,J i l i nN o r m a lU n i v e r s i t y,S i p i n g,J i l i n1 3 6 0 0 0,C h i n a;2.K e yL a b o r a t o r yo fR e m o t eS e n s i n go fE c o l o g i c a lE n v i r o n m e n t i nJ i l i nP r o v i n c e,J i l i nN o r m a lU n i v e r s i t y,S i p i n g,J i l i n1 3 6 0 0 0,C h i n a;3.N o r t h e a s tI n s t i t u t eo fG e o g r a p h ya n dA g r o e c o l o g y,C h i n e s eA c a d e m yo fS c i e n c e s,C h a n g c h u n,J i l i n1 3 0 1 0 2,C h i n a)A b s t r a c t:O b j e t i v eT h ec h a r a c t e r i s t i c sa n dd r i v i n g m e c h a n i s m so fs p a t i a la n dt e m p o r a lc h a n g e si nt h ee c o l o g i c a l v u l n e r a b i l i t yo f t h eS o n g n e nP l a i nf r o m1 9 8 0t o2 0 2 0w e r es t u d i e di no r d e rt op r o v i d eas c i e n t i f i cb a s i s f o r e c o s y s t e m m a n a g e m e n t a n ds u s t a i n a b l ed e v e l o p m e n t i n t h e r e g i o n.M e t h o d sF r o mt h ep e r s p e c t i v eo fn a t u r a lf a c t o r sa n dh u m a nf a c t o r s,t h e“p r e s s u r e-s t a t e-r e s p o n s e”(P S R)m o d e lw a su s e dt ob u i l da ne v a l u a t i o n i n d e xs y s t e m.P r i n c i p a lc o m p o n e n ta n a l y s i sw a sc o m b i n e dw i t ht h e M o r a nIi n d e x.R e s u l t s T h ee c o l o g i c a lv u l n e r a b i l i t yo ft h eS o n g n e nP l a i n w a sd o m i n a t e db ye x t r e m e,s e v e r e,a n d m o d e r a t ev u l n e r a b i l i t i e s,a c c o u n t i n gf o r8 1.7 2%o f t h e t o t a lv u l n e r a b i l i t y.T h eo v e r a l l l e v e lo fe c o l o g i c a lv u l n e r a b i l i t ys h o w e da nu p w a r dt r e n db e f o r e2 0 0 0,a n dt h eo v e r a l l l e v e lo fe c o l o g i c a lv u l n e r a b i l i t yd e c l i n e da f t e r2 0 0 0.T h ee c o l o g i c a lv u l n e r a b i l i t yo ft h eS o n g n e n P l a i n w a sd o m i n a t e db yl i g h t(6 2.1 9%)a n d m o d e r a t e(9 1.3 0%)v u l n e r a b i l i t yc h a n g e s.C h a n g e si ne c o l o g i c a lv u l n e r a b i l i t yw e r ed u et oc o n v e r s i o n sb e t w e e ng r a s s l a n d s,w a t e r s,u n u s e d l a n d,c o n s t r u c t i o n l a n d,a n dc u l t i v a t e d l a n d.S p a t i a l a u t o c o r r e l a t i o ne x i s t e d i ne c o l o g i c a l v u l n e r a b i l i t y,a n dl o w-l o wa g g r e g a t i o nw a sm a i n l yl o c a t e di nt h ef o r e s t l a n da r e ai nt h ew e s t e r np a r to f t h es t u d ya r e a.C o n c l u s i o nT h ee c o l o g i c a lv u l n e r a b i l i t yo fS o n g n e nP l a i nd u r i n gt h ep a s t4 0y e a r sg e n e r a l l ys h o w e da f l u c t u a t i n ga n dd e c r e a s i n gt r e n d.T h es p a t i a le c o l o g i c a lv u l n e r a b i l i t yg r a d u a l l yi n c r e a s e df r o me a s t t ow e s t.N a t u r a l f a c t o r s s u c ha s t e m p e r a t u r ea n dp r e c i p i t a t i o n,a n dh u m a n f a c t o r s s u c ha s l a n du s eh a dag r e a t i m p a c to ne c o l o g i c a l v u l n e r a b i l i t y.K e y w o r d s:e c o l o g i c a lv u l n e r a b i l i t y;s p a t i a la u t o c o r r e l a t i o n;s p a t i a lp r i n c i p a lc o m p o n e n ta n a l y s i s;e c o l o g i c a le n v i r o n m e n t;S o n g n e nP l a i n 生态脆弱性理论起源于2 0世纪提出的生态过渡带理论,是生态系统自身固有的属性,是指在特定的时空尺度上,由于外界环境变化生态系统所表现出的抗干扰能力以及自我恢复的能力1。近几十年来,随着全球环境变化和人类活动加剧,生态系统的自我调节能力逐渐减弱破坏了原本平衡的生态系统,导致生态脆弱性问题日益凸显,因此生态脆弱性评估工作成为全球研究的热点2。生态脆弱性能够在一定程度上掌握生态系统变化的原因和特点,可为生态可持续发展和生态保护提供很大程度的借鉴作用和科学依据,并为生态恢复工程提供决策支持3。国内外学者针对生态脆弱性从不同区域和不同尺度开展相关的研究工作,运用不同的评价模型和方法,取得了大量研究成果4-5。I p p o l i t o等6运用脆弱性指数方法,对意大利北部的两个河流系统在面对多种压力源时的生态脆弱性进行评估,揭示生态系统对多种压力源的潜在响应。张金茜等7选取8个定量指标构建生态脆弱性评价体系,基于S R P模型结合景观格局指数和生态系统服务类型,从3种不同的尺度对2 0 0 2和2 0 1 4年甘肃白龙江流域的生态脆弱性的空间关系进行深入分析。上述学者从评价方法和驱动因素等多方面开展研究,加深了对区域的生态脆弱性现状以及分布特征的理解。但也存在一些局限性,在构建指标体系时,忽略各指标间的相关性,将选取的指标全部加入评价模型中,影响最终评价的精度。压力状态响应”模型(P S R)在人地关系及可持续发展中应用广泛,从自然环境本身和人文干扰两个角度出发,反映人类干扰改变了自然环境本身状态,进而通过决策等改善生态 环境,从而促 进生 态 系 统 的 良 性 循环8。因此,亟需采用P S R模型,掌握研究区生态脆弱性变化趋势,分析空间异质性特征,对区域生态脆弱性状况进行综合评估。此外,当前的研究多集中于干旱区9、高寒区1 0、煤炭矿区等1 1自然条件较差或人为干扰剧烈的区域,但是对松嫩平原这类典型黑土农业区的生态脆弱性研究较少。松嫩平原是松花江和嫩江冲击形成的平原,是中国重要的粮食生产基地,对于国家及区域粮食安全和生态安全至关重要。松嫩平原生态环境十分脆弱,它是世界三大盐碱土分布区之一,且盐碱程度不断加剧造成农田减产、草地退化、湿地锐减,生态环境日趋恶化影响松嫩平原地区社会经济的发展1 2。近年来,为满足当地社会发展的需求和国家粮食生产要求,松嫩平原受到不合理的生产方式影响,该区的土地利用变化显著,使得该区域黑土层退化严重,生态脆弱性表现出较强的典型性和复杂性,受到社会各界的广泛关注。基于此,本文以 松 嫩 平 原 作 为 研 究 区,基 于P S R模 型 选 取1 9 8 02 0 2 0年气温、降水、坡度等1 3个指标因子构建评价指标体系,运用空间主成分分析法和空间自相关法对松嫩平原生态脆弱性进行评价,综合分析松嫩平原的生态脆弱性时空演化格局、空间异质性特征以及驱动机制,旨在为松嫩平原的生态安全、生态环境保护和修复提供理论依据。1 研究区概况松嫩平原(4 2 4 9 4 9 1 3 N,1 2 1 3 8 1 2 8 1 9 E)与大兴安岭和小兴安岭接壤,位于吉林省西部和黑龙江省西南部,总面积约为2.2 41 05k m2。研究区共包括9个市级行政区,分别是吉林省的长春市、四平市、松原市和白城市以及黑龙江省的哈尔滨市、大庆市、绥化市、黑河市和齐齐哈尔市1 3-1 4。该区域属于温带大陆性半湿润半干旱季风气候,四季分明,降水主要集中在夏季,年均降水量为3 5 0 6 5 0mm,年均气温为-2 7,年均蒸发量为12 0 016 0 0mm。研究区地势平缓,土壤以黑土和黑钙土为主1 5。2 数据来源及研究方法2.1 数据来源本研究选取1 9 8 0,1 9 9 0,2 0 0 0,2 0 1 0和2 0 2 0年5个时期的数据(表1)。不同的数据来源和精度导致选取的评价指标存在空间差异,为方便后续计算,本研究对数据进行初步整理并以栅格作为基本评价单元,将分辨率不同的所有数据重新采样为1 0 0 m1 0 0m,并将空间数据设置为K r a s o v s k y坐标A l b e r s投影,保证各指标具有良好的空间重合性1 6。763第4期 宋春威等:1 9 8 02 0 2 0年松嫩平原生态脆弱性时空分异表1 生态脆弱性评价指标来源及处理T a b l e1 S o u r c ea n dp r o c e s s i n go f e v a l u a t i n g i n d i c a t o r s f o r e c o l o g i c a l v u l n e r a b i l i t y指 标 计算方法及说明 分辨率/比例尺数据来源 地形起伏度数字高程模型提取9 0m地理空间数据云平台河网密度数字高程模型提取9 0m地理空间数据云平台年平均降水量克里金插值法5 0 0m中国气象数据网年均温度克里金插值法5 0 0m中国气象数据网植被覆盖度像元二分模型5 0 0m美国航天局网站植被净初级生产力C A S A模型5 0 0m美国航天局网站土地利用类型重分类为耕地、林地、草地、水域、建筑用地和未利用地6类3 0m中国科学院资源环境科学数据中心土壤类型分布土壤数据库11 0 0万中国科学院资源环境科学数据中心土壤侵蚀强度土壤类型矢量数据11 0 0万中国科学院资源环境科学数据中心人口密度人口数量/土地面积 中国统计年鉴人均G D PG D P/人口数量 中国统计年鉴人均耕地面积耕地面积/人口数量 中国统计年鉴人均林地面积林地面积/人口数量 中国统计年鉴2.2 生态脆弱性指标体系构建P S R模型即压力状态响应模型,通过人类和自然环境间相互作用相互影响,分析生态环境内部因果关系1 7-1 8。压力指标反映人类活动对环境所造成的影响,本文通过选取人口密度和土地利用类型来反映社会活动对松嫩平原环境带来的压力。状态指标是指生态系统内部各种生态因子长期作用的结果1 9。本文通过选取平均降水量和平均温度反映气候脆弱因子,地形起伏度和河网密度反映地貌脆弱因子,土壤侵蚀强度数据反映水土流失状况,土壤类型分布数据反映松嫩平原土壤脆弱性因子,植被覆盖度和植被净初级生产力反映植被生长状况和松嫩平原生态系统的质量。响应指标是指人类为了提高生态环境质量,采取一系列措施减缓或阻止给生态环境带来的负面影响。本文通过选取人均G D P值、人均耕地面积和人均林地面积来反映生态脆弱性评价的响应因子。2.3 指标数据标准化在生态脆弱性评价过程中,由于选取的各评价指标的性质和固有属性不同,所以无法对生态系统脆弱性直接评价2 0。因此,本研究对不同评价指标数据统一进行标准化处理,通过选取的不同评价指标对松嫩平原生态环境的影响,将各指标划分为正向指标和负向指标两类2 1。计算公式为:RT=X-Xm i nXm a x-Xm i n(1)RG=Xm a x-XXm a x-Xm i n(2)式中:RT为正向指标标准化值;RG为负向指标标准化值;Xm a x,Xm i n分别为指标X的最大值和最小值。研究选取的1 3个指标通过公式得到正向指标8个和负向指标5个。正向指标包括:人口密度、地形起伏度、土地利用类型、土壤侵蚀强度、土壤类型分布、人均G D P值、人均耕地面积、人均林地面积。负向指标包括:植被初级净生产力、植被覆盖度、河网密度、平均温度、平均降水量。2.4 生态脆弱性指数空间主成分分析(S P C A)是基于数理统计手段,旋转原始特征光谱空间的坐标轴,将相关的多个空间指标数据转换为互不关联的几个综合性指标数据,达到减少数据冗余的目的2 2。通过S P C A法计算P S R模型选取的各指标权重,可在一定程度上减少人为赋予指标权重,同时对松嫩平原1 9 8 0,1 9 9 0,2 0 0 0,2 0 1 0和2 0 2 0年5a生态脆弱性进行主成分分析,根据主成分大于8 5%的前5个主成分,计算松嫩平原生态脆弱性指数,计算公式为2 3:An=R1C1n+R2C2n+RnCk n(3)式中:An是生态脆弱性指数;R1,R2Rn是对应指标权重;C1n,C2nCk n表示累计贡献率大于8 5%的主成分;n为年份。生态脆弱性指数越小,表明研究区生态环境越好;反之,越差。为了便于比较1 9 8 0,1 9 9 0,2 0 0 0,2 0 1 0和2 0 2 0年5期数据的生态脆弱性,本文将5期数据的结果进行标准化处理,计算公式为2 4:SE=E-Em i nEm a x-Em i n1 0(4)式中:SE为生态脆弱性指数标准化值;E为研究区生态脆弱性指数;Em a x,Em i n分别为研究区范围内生态脆弱性指数的最大值和最小值。参考松嫩平原生态环境特征及北方生态区生态863 水土保持通报 第4 3卷脆弱性直方图分布和标准差,综合考虑松嫩平原特殊的生态环境特征,基于自然断点法对松嫩平原生态脆弱性指数进行分级,将生态脆弱性进行分级划分为5级(表2)。表2 生态脆弱性分级标准T a b l e2 C l a s s i f i c a t i o ns t a n d a r d s f o r e c o l o g i c a l v u l n e r a b i l i t y脆弱性等级脆弱性指数生态特征 微度脆弱性0,1.8生态环境状况好,具有很强的抵抗外界干扰能力轻度脆弱性(1.8,3.1生态环境状况良好,具有较强的抵抗外界干扰能力中度脆弱性(3.1,5.7生态环境状况一般,具有一般的抵抗外界干扰能力重度脆弱性(5.7,7.5生态环境状况较差,具有较差的抵抗外界干扰能力极度脆弱性(7.51 0生态环境状况差,具有很差的抵抗外界干扰能力2.5 空间自相关分析空间自相关分析是一系列空间数据分析方法和技术的集合,是监测空间属性及其变化的相关性程度的重 要 方 法2 5。本 研 究 针 对 单 一 要 素 选 取 局 部M o r a nI指数来表达空间集散结构及分布规律。空间关联局部指标聚类图(L I S A)是在局部M o r a nI指数基础上计算同类对象的相似度和不同对象间的差异的结果,主要包含5种集聚模式:高高集聚、高低离散、低高离散、低低集聚和不显著2 6。局部M o r a nI指数的计算公式为:I=(Xi-X)S2jWi j(Xj-X)(5)式中:I为M o r a nI指数;Xi,Xj分别为第i个、第j个评价单元内的脆弱性指数的均值;X为全部评价单元的脆弱性均值;Wi j为空间权重矩阵;S为空间权重矩阵各元素之和。3 结果与分析3.1 松嫩平原生态脆弱性空间分异从空间上来看,1 9 8 02 0 2 0年松嫩平原生态脆弱性总体呈现自西向东呈现逐渐递减的趋势(图1)。这可能与松嫩平原所处的地理位置和自然环境的本底条件息息相关。研究区自西向东呈现多年蒸散量逐渐增加的趋势,土壤类型也由研究区西部盐碱土向研究区中部黑钙土和研究区东部黑土过渡,植被类型也呈现从草原向林地依次过渡。2 0 0 0年松嫩平原的极度和重度脆弱性占研究区总面积的6 8.3 1%,而中度脆弱性研究区面积的1 5.1 4%。1 9 8 0年,松嫩平原生态环境脆弱性等级的面积大小由高到低依次为:中度脆弱性轻度脆弱性极度脆弱性重度脆弱性微度脆弱性,其对应面积比例分别为:2 6.4 5%2 1.6 9%2 0.0 4%1 9.4 9%1 2.2 3%。而 到2 0 2 0年,松嫩平原生态环境脆弱性等级的面积大小由高到低依次为:中度脆弱性重度脆弱性极度脆弱性微度脆弱性轻度脆弱性,其对应面积比例分别 为:5 0.6 0%1 6.9 4%1 4.1 8%1 0.0 8%8.2 0%。研究区西部地区以极度脆弱性为主,包括白城市、大庆市等地区,其主要原因是该区土壤主要是由风沙土组成,持水能力较差,土壤水分极易流失。由于松嫩平原的研究区西部地处内陆,离海洋较远,具有很强的大陆气候,年均降水量少蒸散量大,加剧了研究区西部的土壤缺水程度。轻度和微度脆弱性主要分布在研究区东部地区,包括四平市和长春市中东部、哈尔滨市东部、黑河市以及绥化市东部,这主要是因为该区与研究区西部相比,气候较为湿润,植被覆盖度较高,所以土壤涵养水源的能力较强。重度和中度脆弱性分布在研究区中部地区,包括松原市、齐齐哈尔市、绥化市中西部以及哈尔滨市中西部,主要是因为该区处于研究区东部和研究区西部的过渡区域,土壤处于中度缺水状态。松嫩平原是典型的农业区,面积广阔,地势平坦,研究区自西向东降水量逐渐增加,土壤水分含量也与降水量相关。3.2 松嫩平原生态脆弱性时间序列变化特征1 9 8 02 0 2 0年,松嫩平原中度脆弱性面积呈波动增加趋势,变化幅度为9 1.3 0%。而微度、轻度、重度和极度脆弱性面积均呈现波动减少的趋势,变化幅度分别为1 7.5 8%,6 2.1 9%,1 3.5 3%和2 9.2 4%(图2)。以2 0 0 0年为分割点,1 9 8 02 0 0 0年,松嫩平原地区微度、轻度和中度脆弱性面积呈降低趋势,总体上分别下降6 0.8 3%,3 7.4 8%和4 2.7 6%,而重度和极度脆弱性面积呈上升趋势,总体上分别升高3 7.7 9%和7 5.5 5%。齐齐哈尔市、大庆市、绥化市、松原市、长春市、四平市、黑河市南部和哈尔滨市中西部由原来的微度、轻度、中度和重度脆弱性转化为轻度、中度、重度和极度脆弱性。这是因为中国社会经济高速发展,人类对生态系统的干扰作用增强,如农业生产为主的区域,主要是对森林和草原植被区的开垦,导致土壤的有机质含量和土地利用类型的变化,从而引起生态963第4期 宋春威等:1 9 8 02 0 2 0年松嫩平原生态脆弱性时空分异环境和生态脆弱性的变化。2 0 0 02 0 2 0年,松嫩平原极度脆弱性、重度脆弱性面积呈下降趋势,降幅分别为5 9.6 9%和4 5.9 3%,而中度脆弱性、轻度脆弱性和 微 度 脆 弱 性 面 积 呈 上 升 趋 势,增 幅 分 别 为2 3 4.2 1%,3 9.5 3%和5 2.4 8%。其中,齐齐哈尔市、绥化市、松原市、大庆市、黑河市,以及四平市和哈尔滨市西部极度脆弱性和重度脆弱性大幅降低为中度脆弱性,这是因为针对区域的生态问题开展了退耕还林还草政策的实施,使原本被破坏环境得到了弥补,从而降低了生态系统的脆弱性。图1 松嫩平原生态脆弱性空间分布F i g.1 S p a t i a l d i s t r i b u t i o no f e c o l o g i c a l v u l n e r a b i l i t y i nS o n g n e nP l a i n图2 松嫩平原生态脆弱性不同等级时间变化特征F i g.2 T i m ec h a n g ec h a r a c t e r i s t i c so fd i f f e r e n t l e v e l so fe c o l o g i c a l v u l n e r a b i l i t y i nS o n g n e nP l a i n3.3 松嫩平原不同土地利用结构的生态脆弱性时空变异不同的土地利用类型间相互转换致使研究区生态脆弱性发生不同程度的变化。2 0 2 0年,依照不同的土地利用面积由低到高的依次为:建筑用地水域草地未利用地林地农田(图34)。从总体上分析,1 9 8 02 0 2 0年建筑用地面积呈现增加的趋势,增幅1 9 5.5 6%。1 9 8 02 0 2 0年水域面积呈减小的趋势,降幅为3 7.8 1%。1 9 8 02 0 2 0年松嫩平原的林地和草地面积呈现先下降后上升的趋势,2 0 0 0年都达最小值,林地和草地面积分别占研究区总面积的1 6.6 8%和4.9 0%。1 9 8 02 0 2 0年,松嫩平原耕地面073 水土保持通报 第4 3卷积呈现先上升后下降的趋势,2 0 0 0年时耕地面积达到最大值占研究区面积的6 4.9 9%。以2 0 0 0年作为分割点,1 9 8 02 0 0 0年 林 地 和 草 地 面 积 分 别 减 少60 9 2.8 0和70 5 6.0 0k m2,降 幅 分 别 为1 4.0 2%和3 9.1 3%,而 耕 地 面 积 增 加1 92 1 9.2 0 k m2,增 幅1 5.2 1%。2 0 0 02 0 2 0年,由于退耕还林还草政策的实施,耕地面积减少73 4 7.2 0k m2,降幅为5.0 5%,林地和草地面积分别增加52 8 6.4 0和28 6 7.2 0k m2,增幅分别为1 4.1 5%和2 6.1 2%,主要原因是人类不断的开垦,导致林地草地转化耕地。研究区东部生态区微度脆弱性、轻度脆弱性和中度脆弱性主要分布在林地、水域和耕地,这主要是因为该区的生态本底条件较好,所处地带区域气候湿润,土壤条件较好,气候因子与土壤因子的空间差异性不大。研究区西部以极度脆弱区为主,大部分集中在未利用地和草地,其原因在于该生态区的生态系统状况受自然因素的气温、降水等因素影响较大,降水稀少蒸发量大,植被覆盖量和物种多样性较低,生态系统较单一。人类活动导致土壤盐渍化和沙漠化,使本就恶劣的生态环境变得更加脆弱,进而改变了该区生态脆弱性。图3 松嫩平原土地利用动态监测结果F i g.3 D y n a m i cm o n i t o r i n gr e s u l t so f l a n du s e i nS o n g n e nP l a i n3.4 松嫩平原生态脆弱性空间集聚特征本研究以栅格为空间研究单元,通过M o r a ns指数分析1 9 8 02 0 2 0年中国松嫩平原生态脆弱性指数的空间相关性,生成集聚图(图5)。生态脆弱性在空间上呈现正相关关系且相关关系较强,表明研究区生态脆弱性呈现出空间集聚性特征,且整体趋势大致相同。近年来,水土开发利用增强及城市化建设的加快,人工建设也进入大规模生态外部性的空间发展阶段,空间集聚性特征也存在不同程度上向外的扩张和迁移。低低集聚特征主要位于研究区东部,土地利用类型以林地为主,该区以微度、轻度和中度脆弱性为173第4期 宋春威等:1 9 8 02 0 2 0年松嫩平原生态脆弱性时空分异主,1 9 8 02 0 2 0年低低集聚区空间基本相似,但不同年份间也存在集聚性的扩张和迁移。高高集聚特征主要位于研究区西部,土地利用类型以未利用地和草地为主,该区以极度和重度脆弱性为主。1 9 8 02 0 2 0年,以2 0 0 0年为分割点,高高集聚特征呈现先增加后下降的趋势,由于人类活动的影响和土地利用类型的变化,高高集聚特征也发生了不同程度的迁移,如大庆市发生的变化特征较大,主要是因为石油的开采,导致该地区人口密度,城镇化和经济发展较快,说明人类活动的集聚程度与脆弱性的空间集聚程度也具有一致性。图4 松嫩平原不同年份土地利用面积比例F i g.4 L a n du s ea r e ar a t i o i nd i f f e r e n t y e a r s i nS o n g n e nP l a i n图5 松嫩平原局部空间自相关L I S A集聚图F i g.5 L o c a l s p a c e s e l f-r e l a t e dL I S Ac l u s t e rd i a g r a mo fS o n g n e nP l a i n4 讨 论4.1 生态脆弱性自然因子驱动力分析松嫩平原作为典型的脆弱区,受多种自然因子和人文因子的综合作用的结果2 7。自然因子是影响区域生态脆弱性的成因之一,在一定程度上决定了空间分布格局2 8。本研究发现,生态脆弱性空间格局分布主要与地形因子、气象因子和土壤因子相关,该结论与前人的研究结论一致,如郭兵等2 9利用综合指数法对中国北方半干旱荒漠草原的生态脆弱性进行分析后得出地形、气温、降水等自然因素等变化强度与脆弱性变化显著相关。本研究通过地形、气温和降273 水土保持通报 第4 3卷水等自然因子相关数据进行分析,松嫩平原的不同地区的气温降水等自然因子不同,这导致该地区的水热平衡也存在较大的差异,这对研究区的生态环境变化具有深远的影响,进而改变了该地区生态系统过程,该结论也印证了前人得出松嫩平原地区的生态脆弱区的结论。如张良侠等3 0对黄土高原区生态系统的脆弱性进行研究后,认为植被覆盖度和降水等自然因子是导致黄土高原区生态脆弱性时空变化的主要驱动因子,且不同地区和不同土地利用类型脆弱性差异较大。因此本研究结果表明自然因素对松嫩平原的生态脆弱性的影响具有一定的说服力3 1。4.2 生态脆弱性人文因子驱动力分析人类活动也是影响生态脆弱性的因子3 2。研究表明,由于人类活动的影响和社会经济的发展导致土地的利用类型发生剧烈的变化,进一步影响了松嫩平原的生态脆弱性,该结论与前人的研究结论一致3 3。如田盛等3 4利用R S与G I S技术,对环杭州湾土地利用变化进行生态脆弱性评价,研究结果表明环杭州湾耕地面积减少建筑用地面积大幅增加,导致城市化进程加快。本研究表明,松嫩平原生态脆弱性的变化,可能与松嫩平原耕地面积不断扩张有关。耕地和建筑用地的变化是人类活动对生态环境影响后的映射,由于人类活动的加剧,耕地和建设用地的面积变化较大。松嫩平原特殊的地理位置以及近年来该地区人口数量的不断增长,大面积的林地、草地、水域转化为耕地和建设用地,这表明社会经济与城乡建设不断发展,城市化进程加快,人类活动影响越来越强烈,这与前人的研究结果一致3 5。如周岩等3 6利用层次分析法对中国辉南县生态脆弱性时空变化进行评估,证明土地利用变化空间分布差异对生态环境脆弱性具有重要影响,林地面积减小转变为耕地和建筑用地,导致生态环境脆弱性发生改变。近年来,中国对“山水田林湖草”等进行系统性保护和恢复措施,通过协调松嫩平原农业生产与生态脆性间的平衡,采取正确的适当的措施抑制生态脆弱性上升的趋势,这符合国家政策对于生态建设和环境保护的支持3 7。4.3 松嫩平原生态环境优化策略松嫩平原是复杂的自然经济社会复合体,在该区应继续维持其生态保护政策,本着“绿水青山就是金山银山”的发展理念,加强宣传教育,增强群众的保护意识,创建全民保护环境的良好氛围3 8;严格控制人类过度开垦林地、草地、水域等行为造成土壤退化、植被覆盖度降低和地表径流减少等现象;继续实施设立生态育林,把握农林牧业结构;在耕作的基础上也应该保护土壤环境质量,保持生态系统的稳定性,加强土地利用和区域生态治理的综合性规划,促进区域综合发展,实施退耕还林还草还牧还湖等生态工程3 9,以维护生态环境稳定性;同时,注意充分利用太阳能、地热能、风能等各种自然资源,促进区域生态系统的功能整合和社会生态经济的协调发展。5 结 论松嫩平原生态环境脆弱性整体以极度、重度和中度脆 弱性为主,三 者的面积均 超过研究区 面 积 的6 6.0 8%。1 9 8 02 0 2 0年,松嫩平原生态脆弱性总体呈现 先 增 加 后 降 低 的 趋 势,以2 0 0 0年 为 分 水 岭,1 9 8 02 0 0 0年生态脆弱性呈现增加趋势。而2 0 0 02 0 2 0年生态脆弱性呈现降低趋势,这表明松嫩平原生态系统稳定性以及环境状况明显好转。松嫩平原生态脆弱性呈现出自西向东逐渐减弱的趋势。研究期内,松嫩平原生态脆弱性轻度和中度脆弱性面积变化较大,变幅分别为6 2.1 9%和9 1.3 0%。生态环境脆弱性存在着显著的空间自相关性,且为显著的正相关,空间集聚性特征也有所迁移和扩张。松嫩平原生态脆弱性时空变异,主要受该地区自然因素和人类活动共同影响。从数理统计分析结果看,生态脆弱性的变化与气温、降水的变化显著相关。同时,应该指出,社会经济因素对生态脆弱性起到促进作用,其影响有逐渐加强的趋势。松嫩平原生态脆弱性评价是一个复杂的过程