基于土地利用变化的黄土丘陵区生境质量时空演变特征

ＡＣＴＡＥＣＯＬＯＧＩＣＡＳＩＮＩＣＡ

生态学报

Ｖｏｌ．３８，Ｎｏ．２０Ｏｃｔ．，２０１８

ＤＯＩ：１０．５８４６／ｓｔｘｂ２０１７１０１０１８００

刘春芳，王川．基于土地利用变化的黄土丘陵区生境质量时空演变特征———以榆中县为例．生态学报，２０１８，３８（２０）：７３００⁃７３１１．ｃｏｕｎｔｙ．ＡｃｔａＥｃｏｌｏｇｉｃａＳｉｎｉｃａ，２０１８，３８（２０）：７３００⁃７３１１．

ＬｉｕＣＦ，ＷａｎｇＣ．Ｓｐａｔｉｏ⁃ｔｅｍｐｏｒａｌｅｖｏｌｕｔｉｏｎｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｈａｂｉｔａｔｑｕａｌｉｔｙｉｎｔｈｅＬｏｅｓｓＨｉｌｌｙＲｅｇｉｏｎｂａｓｅｄｏｎｌａｎｄｕｓｅｃｈａｎｇｅ：ａｃａｓｅｓｔｕｄｙｉｎＹｕｚｈｏｎｇ

———以榆中县为例

基于土地利用变化的黄土丘陵区生境质量时空演变特征

刘春芳１，２，∗，王　川１，２

１西北师范大学地理与环境科学学院，兰州　７３００７０２甘肃省土地利用与综合整治工程研究中心，兰州　７３００７０

摘要：生境质量是生态系统服务水平的重要表征，是区域生态安全保障和人类福祉提升的关键环节，评估土地利用变化引起的生境质量时空演变可为区域生态保护和土地管理提供科学依据。以典型黄土丘陵区县域———榆中县为例，基于ＡｒｃＧＩＳ平台和ＩｎＶＥＳＴ模型探讨了该县１９９５、２００５和２０１５年生境质量空间格局及其演变特征。结论如下：①研究区土地类型和景观格局变化明显，主要土地利用变化是耕地转向建设用地和林地，以及草地和耕地之间的相互转变，这些变化导致景观异质性增强和破碎度上升。②生境质量空间格局一直表现为“南部高、中部居中、北部低”的态势。生境质量上升区域主要在北部，原因是“退耕还林”及“生态北山”工程的实施使得林地面积增加，下降区域主要在中部，主要是由建设用地扩张占用耕地造成的。③生境质量热点以及“Ｈ⁃Ｈ”型聚集区主要位于南部山区，而冷点和“Ｌ⁃Ｌ”聚集区主要位于北部黄土丘陵区，２０１５年生境质量空间集聚程度较１９９５年有所增强，“Ｈ⁃Ｈ”聚集区增加０．６８％，“Ｌ⁃Ｌ”聚集区增加０．６７％。对该类型区域生境质量优化提出了相应的政策建议。

关键词：土地利用变化；ＩｎＶＥＳＴ模型；生境质量；时空格局；榆中县

Ｓｐａｔｉｏ⁃ｔｅｍｐｏｒａｌｅｖｏｌｕｔｉｏｎｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｈａｂｉｔａｔｑｕａｌｉｔｙｉｎｔｈｅＬｏｅｓｓＨｉｌｌｙＲｅｇｉｏｎｂａｓｅｄｏｎｌａｎｄｕｓｅｃｈａｎｇｅ：ａｃａｓｅｓｔｕｄｙｉｎＹｕｚｈｏｎｇＣｏｕｎｔｙ

１ＧｅｏｇｒａｐｈｙａｎｄＥｎｖｉｒｏｎｍｅｎｔＣｏｌｌｅｇｅｏｆＮｏｒｔｈｗｅｓｔＮｏｒｍａｌＵｎｉｖｅｒｓｉｔｙ，Ｌａｎｚｈｏｕ７３００７０，Ｃｈｉｎａ

ＬＩＵＣｈｕｎｆａｎｇ１，２，∗，ＷＡＮＧＣｈｕａｎ１，２

２ＧａｎｓｕＥｎｇｉｎｅｅｒｉｎｇＲｅｓｅａｒｃｈＣｅｎｔｅｒｏｆＬａｎｄＵｔｉｌｉｚａｔｉｏｎａｎｄＣｏｍｐｒｅｈｅｎｓｉｏｎＣｏｎｓｏｌｉｄａｔｉｏｎ，Ｌａｎｚｈｏｕ７３００７０，Ｃｈｉｎａ

Ａｂｓｔｒａｃｔ：Ｈａｂｉｔａｔｑｕａｌｉｔｙｉｓａｎｉｍｐｏｒｔａｎｔｒｅｐｒｅｓｅｎｔａｔｉｏｎａｎｄｒｅｆｌｅｃｔｉｏｎｏｆｔｈｅｌｅｖｅｌｏｆｅｃｏｓｙｓｔｅｍｓｅｒｖｉｃｅｓ，ａｓｗｅｌｌａｓｔｈｅｋｅｙｌｉｎｋｉｎｇｕａｒａｎｔｅｅｉｎｇｒｅｇｉｏｎａｌｅｃｏｌｏｇｉｃａｌｓｅｃｕｒｉｔｙａｎｄｉｍｐｒｏｖｉｎｇｈｕｍａｎｗｅｌｆａｒｅ．Ａｓｓｅｓｓｍｅｎｔｏｆｔｈｅｔｅｍｐｏｒａｌａｎｄｓｐａｔｉａｌｌａｎｄｍａｎａｇｅｍｅｎｔ，ｅｓｐｅｃｉａｌｌｙｉｎｅｃｏｌｏｇｉｃａｌｌｙｆｒａｇｉｌｅａｒｅａｓ．Ｆｏｒａｌｏｎｇｔｉｍｅ，ｔｈｅＬｏｅｓｓＨｉｌｌｙＲｅｇｉｏｎｈａｓｂｅｅｎａｔｙｐｉｃａｌｃｈａｎｇｅｓｉｎｈａｂｉｔａｔｑｕａｌｉｔｙｃａｕｓｅｄｂｙｌａｎｄｕｓｅｃｈａｎｇｅｃａｎｐｒｏｖｉｄｅａｓｃｉｅｎｔｉｆｉｃｂａｓｉｓｆｏｒａｒｅｇｉｏｎ′ｓｅｃｏｌｏｇｉｃａｌｐｒｏｔｅｃｔｉｏｎａｎｄｅｃｏｌｏｇｉｃａｌｌｙｖｕｌｎｅｒａｂｌｅａｒｅａｉｎＣｈｉｎａ，ａｎｄｓｉｎｃｅｔｈｅ１９９０ｓｔｈｅｓｏｃｉａｌａｎｄｅｃｏｎｏｍｉｃｃｏｎｄｉｔｉｏｎｓｏｆＣｈｉｎｅｓｅｃｏｕｎｔｉｅｓｈａｖｅｕｎｄｅｒｇｏｎｅｒａｐｉｄｄｅｖｅｌｏｐｍｅｎｔ．Ｃｏｎｓｉｄｅｒｉｎｇｔｈｅｓｅｆａｃｔｏｒｓ，ｗｅｕｓｅｄａｔｙｐｉｃａｌＬｏｅｓｓＨｉｌｌｙＲｅｇｉｏｎｃｏｕｎｔｙ，ＹｕｚｈｏｎｇＣｏｕｎｔｙｉｎＧａｎｓｕｐｒｏｖｉｎｃｅ，ａｓａｎｅｘａｍｐｌｅ，ａｎｄｂａｓｅｄｏｎＡｒｃＧＩＳｐｌａｔｆｏｒｍａｎｄｔｈｅＩｎＶＥＳＴｍｏｄｅｌ，ｕｓｉｎｇｔｈｅｌａｎｄｕｓｅｄａｔａｏｆＹｕｚｈｏｎｇＣｏｕｎｔｙｉｎ１９９５，２００５，ａｎｄ２０１５，ｓｔｕｄｉｅｄｆｏｒｔｈｅｆｉｒｓｔｔｉｍｅｔｈｅｃｈａｎｇｅｓｉｎｌａｎｄｕｓｅ，ｌａｎｄｓｃａｐｅｐａｔｔｅｒｎｓ，ａｎｄｈａｂｉｔａｔｑｕａｌｉｔｙｆｒｏｍ１９９５ｔｏ２０１５；ｗｅｔｈｅｎａｎａｌｙｚｅｄｔｈｅｒｅｓｐｏｎｓｅｏｆｈａｂｉｔａｔｑｕａｌｉｔｙｔｏｌａｎｄｕｓｅｃｈａｎｇｅａｎｄｌａｎｄｓｃａｐｅｐａｔｔｅｒｎｃｈａｎｇｅ，ｕｓｉｎｇ

基金项目：国家自然科学基金项目（４１７７１１３０）

收稿日期：２０１７⁃１０⁃１０；　　网络出版日期：２０１８⁃０７⁃０９

∗通讯作者Ｃｏｒｒｅｓｐｏｎｄｉｎｇａｕｔｈｏｒ．Ｅ⁃ｍａｉｌ：ｌｉｕｃｈｕｎｆａｎｇ＠ｎｗｎｕ．ｅｄｕ．ｃｎ

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

　２０期　　　刘春芳　等：基于土地利用变化的黄土丘陵区生境质量时空演变特征———以榆中县为例　７３０１

ｈｏｔｓｐｏｔａｎａｌｙｓｉｓａｎｄＬＩＳＡｃｌｕｓｔｅｒｉｎｇａｎａｌｙｓｉｓｔｏｓｔｕｄｙａｎｄｉｄｅｎｔｉｆｙｔｈｅｓｐａｔｉａｌｄｉｓｔｒｉｂｕｔｉｏｎ，ｃｌｕｓｔｅｒｉｎｇｔｙｐｅ，ａｎｄｅｖｏｌｕｔｉｏｎｆｏｌｌｏｗｉｎｇ：１）Ｄｕｒｉｎｇｔｈｅｓｔｕｄｙｐｅｒｉｏｄ，ｂｏｔｈｔｈｅｌａｎｄｕｓｅｔｙｐｅａｎｄｌａｎｄｓｃａｐｅｐａｔｔｅｒｎｃｈａｎｇｅｄｏｂｖｉｏｕｓｌｙ，ａｎｄｔｈｅｍａｉｎｌａｎｄｕｓｅｃｈａｎｇｅｓｉｎｔｈｅｓｔｕｄｙａｒｅａｗｅｒｅｔｈａｔｃｕｌｔｉｖａｔｅｄｌａｎｄｗａｓｃｏｎｖｅｒｔｅｄｔｏｃｏｎｓｔｒｕｃｔｉｏｎｌａｎｄａｎｄｆｏｒｅｓｔ，ａｎｄｔｈｅｒｅｗａｓｍｕｔｕａｌ

ｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｈａｂｉｔａｔｑｕａｌｉｔｙｏｆｔｈｅＬｏｅｓｓＨｉｌｌｙＲｅｇｉｏｎｉｎ１９９５，２００５，ａｎｄ２０１５．Ｔｈｅｍａｉｎｃｏｎｃｌｕｓｉｏｎｓａｒｅｔｈｅ

ｃｏｎｖｅｒｓｉｏｎｂｅｔｗｅｅｎｇｒａｓｓｌａｎｄａｎｄｃｕｌｔｉｖａｔｅｄｌａｎｄ；ｔｈｅｓｅｃｈａｎｇｅｓｌｅｄｔｏａｎｉｍｐｒｏｖｅｍｅｎｔｉｎｌａｎｄｓｃａｐｅｈｅｔｅｒｏｇｅｎｅｉｔｙａｎｄｔｈｅｄｅｇｒｅｅｏｆｆｒａｇｍｅｎｔａｔｉｏｎ；２）Ｔｈｅｓｐａｔｉａｌｄｉｓｔｒｉｂｕｔｉｏｎｏｆｈａｂｉｔａｔｑｕａｌｉｔｙｗａｓ：ｈｉｇｈｉｎｔｈｅｓｏｕｔｈｅｒｎａｒｅａｓ，ｍｅｄｉｕｍｉｎｔｈｅｍｉｄｄｌｅａｒｅａｓ，ａｎｄｌｏｗｉｎｔｈｅｎｏｒｔｈｅｒｎａｒｅａｓ．Ｄｕｒｉｎｇｔｈｅｓｔｕｄｙｐｅｒｉｏｄ，ｔｈｅａｒｅａｉｎｗｈｉｃｈｈａｂｉｔａｔｑｕａｌｉｔｙｉｍｐｒｏｖｅｄｗａｓｍａｉｎｌｙｉｎｔｈｅｎｏｒｔｈｏｆｔｈｅｃｏｕｎｔｙｂｅｃａｕｓｅｏｆｔｈｅｉｍｐｌｅｍｅｎｔａｔｉｏｎｏｆｔｈｅ“ＧｒａｉｎｆｏｒＧｒｅｅｎＰｒｏｊｅｃｔ”ａｎｄｔｈｅ“ＥｃｏｌｏｇｉｃａｌＢｅｉｓｈａｎ”ｐｒｏｊｅｃｔ；ｔｈｅａｒｅａｏｆｄｅｔｅｒｉｏｒａｔｉｏｎｗａｓｐｒｉｍａｒｉｌｙｌｏｃａｔｅｄｉｎｔｈｅｃｅｎｔｅｒｏｆｔｈｅｃｏｕｎｔｙ，ｍａｉｎｌｙｄｕｅｔｏｔｈｅｅｘｐａｎｓｉｏｎｏｆ

ｃｏｎｓｔｒｕｃｔｉｏｎｌａｎｄ，ｗｈｉｃｈｏｃｃｕｐｉｅｄｔｈｅｆａｒｍｌａｎｄａｎｄｌｅｄｔｏａｎｉｎｃｒｅａｓｅｄｄｅｇｒｅｅｏｆｆｒａｇｍｅｎｔａｔｉｏｎａｎｄａｄｅｃｌｉｎｅｉｎｓｏｕｔｈｅｒｎｍｏｕｎｔａｉｎａｒｅａ，ｗｈｉｌｅｔｈｅｃｏｌｄｓｐｏｔａｒｅａｓａｎｄａｒｅａｓｏｆ“Ｌｏｗ⁃Ｌｏｗ”ｈａｂｉｔａｔｑｕａｌｉｔｙｔｙｐｅｃｌｕｓｔｅｒｓｗｅｒｅｍａｉｎｌｙｌｏｃａｔｅｄ

ｃｏｎｎｅｃｔｉｖｉｔｙ；３）Ｔｈｅｈｏｔｓｐｏｔａｒｅａｓａｎｄａｒｅａｓｏｆ“Ｈｉｇｈ⁃Ｈｉｇｈ”ｈａｂｉｔａｔｑｕａｌｉｔｙｔｙｐｅｃｌｕｓｔｅｒｓｗｅｒｅｍａｉｎｌｙｌｏｃａｔｅｄｉｎｔｈｅｉｎｔｈｅＬｏｅｓｓＨｉｌｌｙＲｅｇｉｏｎｉｎｔｈｅｎｏｒｔｈ．Ｃｏｍｐａｒｅｄｗｉｔｈ１９９５，ｔｈｅｄｅｇｒｅｅｏｆｈａｂｉｔａｔｑｕａｌｉｔｙｃｌｕｓｔｅｒｉｎｇｉｎ２０１５ｗａｓｅｎｈａｎｃｅｄ：Ｈｉｇｈ”ｔｙｐｅｃｌｕｓｔｅｒｉｎｇａｒｅａｉｎｃｒｅａｓｉｎｇｂｙ０．６８％ａｎｄｔｈｅ“Ｌｏｗ⁃Ｌｏｗ”ｔｙｐｅｃｌｕｓｔｅｒｉｎｇａｒｅａｉｎｃｒｅａｓｉｎｇｂｙ０．６７％．Ｔｈｅｐａｐｅｒａｌｓｏｐｕｔｓｆｏｒｗａｒｄｐｏｌｉｃｙｓｕｇｇｅｓｔｉｏｎｓｆｏｒｔｈｅｏｐｔｉｍｉｚａｔｉｏｎｏｆｈａｂｉｔａｔｑｕａｌｉｔｙｉｎｓｉｍｉｌａｒａｒｅａｓ．

ｔｈｅａｒｅａｓｏｆｂｏｔｈｔｈｅ“Ｈｉｇｈ⁃Ｈｉｇｈ”ｔｙｐｅｃｌｕｓｔｅｒｉｎｇａｎｄｔｈｅ“Ｌｏｗ⁃Ｌｏｗ”ｔｙｐｅｃｌｕｓｔｅｒｉｎｇｉｎｃｒｅａｓｅｄｆｒｏｍ１９９５；ｔｈｅ“Ｈｉｇｈ⁃

ＫｅｙＷｏｒｄｓ：ｌａｎｄｕｓｅｃｈａｎｇｅ；ＩｎＶＥＳＴｍｏｄｅｌ；ｈａｂｉｔａｔｑｕａｌｉｔｙ；ｓｐａｔｉａｌａｎｄｔｅｍｐｏｒａｌｐａｔｔｅｒｎ；ＹｕｚｈｏｎｇＣｏｕｎｔｙ

工业化、城市化的快速增长与资源配置流动性的不断增强对区域土地利用空间结构产生了深远影响，同时也对区域生境产生了较大冲击。土地利用变化会改变区域生境的结构与组成要素，随之对生境斑块之间的物质流、能量流循环过程产生影响，进而改变了区域生境的生产能力和服务功能［１］。生境质量是指在一定的时间和空间中生态系统提供适宜个体与种群持续发展与生存条件的能力，是对人类生存繁衍、生产生活的适宜程度，是一切生态系统功能和服务的前提和基础，更是人类福祉提升的关键环节与地理学领域研究的热点议题［４⁃７］。

［２⁃３］

境质量的状态、趋势及其对人类福祉的影响，是当前国际生态系统服务研究领域的重要任务，也是国际生态学

早期的研究侧重于对野生动物生境质量的评价，以及人类活动对生境质量的影响。如Ｇｏｅｒｔｚ等［８］通过

。评估、模拟和预测生

草原棉鼠种群数量来评价生境质量的高低，Ｔｈｉｅｌ［９］探讨了乡村道路密度与狼的生境适宜性的关系，Ｅｎｇｅｌ等

［１０］

研究了人类活动对鱼类生境结构和功能的影响。随着国际地圈生物圈计划（ＩＧＢＰ）、国际全球环境变化人

文因素计划（ＩＨＤＰ）、世界气候研究计划（ＷＣＲＰ）等国际项目将土地利用／覆被变化（ＬＵＣＣ）作为全球变化的核心研究问题，国内外学者开始广泛关注土地利用变化对区域生境质量的影响。学者们探讨了土地利用、气候变化与生物多样性之间关系，指出了土地转化和气候变化已经对生物多样性和相关生态系统服务产生重大影响。农用地面积的增长会使区域生境受损，也会引起生境破碎化并增加对周围自然或半自然生境的农业化学品的投入［１１］。许多研究也关注了土地利用变化对野生动物生境的影响，如Ｊｅｔｚ等［１２］分析了在未来土地利用变化下鸟类可能遭受因气候变化和土地转换而带来的危害。Ｏｔｔｏ等［１３］的实证研究表明土地利用的变化降低了蜜蜂对生境的适应性。Ｓｉｅｂｅｒ等［１４］则研究了苏联解体后全国土地利用变化对大型哺乳动物的栖息地生境的影响。２０００年以来，国内学者针对土地利用变化产生的自然生境影响进行了实证研究［１５⁃１６］，近几年开始一是基于模型的评价方法，如ＨＩＳ模型［１９］，ＩｎＶＥＳＴ模型［２０］、ＳｏｌＶＥＳ模型［２１］等；二是基于指标体系的方法［２２⁃２３］。总体来看，早期的生境质量研究以单一时点的野外动植被调查数据为基础，建立包括生物丰度指数、植被覆盖度等生境质量评价指标系统和评分标准测算生境质量指数，对区域或群落生境质量水平进行静态研究。由于数据积累及研究方法等的限制，对生境质量动态变化及其空间集聚状况关注较少。近些年来，

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

关注区域、流域等尺度土地利用变化所带来的生境质量效应［１７⁃１８］。目前对生境质量的研究方法主要分两类：

７３０２　生　态　学　报　　　３８卷　

数理模型方法以及３Ｓ技术的应用，学者们开始关注生境质量的时空差异，尤其是ＩｎＶＥＳＴ模型的使用，对国家、省、流域等较大尺度生境质量的研究逐渐增多。

黄土丘陵区作为我国典型的生态脆弱区，集中了大量发展需求迫切的贫困人口，社会经济发展与生态保护之间的矛盾非常突出［２４］，１９９０年以来，中国县域经济开始得到较快发展，其生境质量变化也最为显著和典型，而已有研究对县域尺度的关注较少。基于此，本文以榆中县为例，通过对该县１９９５—２０１５年生境质量空间格局及演变特征进行分析，探索了生境质量的变化规律及其对黄土丘陵区土地利用变化的响应。这一方面能丰富黄土丘陵区生境质量变化的理论探索，拓展县域尺度的生境质量研究；另一方面也能为黄土丘陵区生态安全格局构建和空间优化提供重要的理论基础与实践指导。１　研究区概况

榆中县地处黄土高原西部、甘肃省中部、兰州市东郊（１０３°４９′１５″—１０４°３４′４０″Ｅ，３５°３４′２０″—３６°２６′３０″４４．４５万（２０１５年）。该县由地势地貌迥异的３个区域组成（图１）：南部为高寒石质山区，气候湿润，最高海拔保护区，区内植被茂密，是全县的生态屏障；北部为干旱少雨、植被稀疏、沟壑纵横的黄土丘陵区，最高海拔２４９５ｍ；中部是川塬河谷区，地势平坦、交通便利、光照海拔介于１３６３—２０００ｍ。该县属于典型的温带干旱、半干旱大陆性气候，年均气温６．７℃，年均降雨量区域差异明显，从南部的５００ｍｍ左右向北逐渐减少到２００ｍｍ左右，降雨量季节分配不均，７０％集中在７—９月，蒸发量１３５０ｍｍ左右。

２　数据获取与研究方法２．１　数据获取与预处理

选择１９９５、２００５和２０１５年３个时段，对研究区的生境质量空间格局及演变特征进行分析。土地利用数据源于美国地质调查局（ＵＳＧＳ）对地观测资源与科学中心提供的空间分辨率为３０ｍ、秋季无云覆盖的ＬａｎｄｓａｔＴＭ遥感数据（行列号１３０／３５）（ｈｔｔｐ：／／ｇｌｏｖｉｓ．ｕｓｇｓ．ｇｏｖ／）。为提高解译精度，基于２０１５年榆中县土地利用变更１９８４坐标系，通过野外实地考察和谷歌地球高清影像验证，三期分类数据判读精度优于９３％。参考土地利用现状分类标准（ＧＢ／Ｔ２１０１０—２００７）及研究区土地利用特点建立二级分类体系（表１）。数据处理和图形生成通过ＡｒｃＭａｐ１０．２和ＧｅｏＤａ１．１０软件共同实现。２．２　研究方法

２．２．１　土地利用变化及景观格局变化分析

首先基于ＡｒｃＧＩＳ平台分析１９９５—２０１５年土地转移变化，进一步利用Ｆｒａｇｓｔａｔｓ软件计算景观指数研究榆中县景观格局特征和演变趋势。通常，景观中的斑块大小、连通性及景观多样性等对生物多样性及各种生态过程都有重要影响［２５］。结合类似研究［１７⁃１８，２０］，考虑黄土丘陵区土地利用斑块数量大、破碎度较高等特点，本文选取斑块数量、破碎度等七个相关指数来度量该区域的景观变化（表２）。利用Ｆｒａｇｓｔａｔｓ４．２软件，可以计算

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

调查数据，首先解译２０１５年数据，再参照２０１５年解译结果对１９９５和２００５年影像进行解译，统一采用ＷＧＳ＿

图１　研究区位置及高程图

Ｆｉｇ．１　Ｌｏｃａｔｉｏｎａｎｄｅｌｅｖａｔｉｏｎｏｆｓｔｕｄｙａｒｅａ

Ｎ），总面积３３０８．６４ｋｍ２，海拔１３６３—３６７０ｍ，总人口

３６７０ｍ，位于西南部的兴隆山自然保护区是国家级自然

充足、水源丰富，是县内主要农业区，县城也位于该区，

　２０期　　　刘春芳　等：基于土地利用变化的黄土丘陵区生境质量时空演变特征———以榆中县为例　７３０３

出除破碎度以外的所有指标。破碎度则是在已有斑块数目的基础上，利用公式（１）计算得到［２５］：

ＦＮ＝（ＮＰ－１）／ＮＣ

式中，ＦＮ为研究区景观破碎度；ＮＰ为景观中各类斑块总数；ＮＣ为景观总面积。

表１　土地利用分类体系Ｔａｂｌｅ１　Ｌａｎｄｕｓｅｃｌａｓｓｉｆｉｃａｔｉｏｎ一级地类Ｃｌａｓｓ１耕地Ｃｕｌｔｉｖａｔｅｄｌａｎｄ二级地类Ｃｌａｓｓ２水浇地旱地有林地灌木林地其他林地天然草地其他草地果园一级地类Ｃｌａｓｓ１水域Ｗａｔｅｒ二级地类Ｃｌａｓｓ２水库坑塘内陆滩涂河流城镇　农村居民点公路　铁路工矿用地裸地（１）

林地Ｆｏｒｅｓｔ建设用地Ｃｏｎｓｔｒｕｃｔｉｏｎｌａｎｄ未利用地Ｕｎｕｓｅｄｌａｎｄ草地Ｇｒａｓｓｌａｎｄ园地Ｇａｒｄｅｎ表２　景观指数及其生态学含义

Ｔａｂｌｅ２　Ｔｈｅｌａｎｄｓｃａｐｅｉｎｄｉｃｅｓａｎｄｔｈｅｉｒｅｃｏｌｏｇｉｃａｌｍｅａｎｉｎｇ

景观指数Ｌａｎｄｓｃａｐｅｉｎｄｉｃｅｓ斑块数量ＮＰ

平均斑块面积ＡＲＥＡ＿ＭＮ形状指数ＳＨＡＰＥ内聚力指数ＣＯＨＥＳＩＯＮ蔓延度指数ＣＯＮＴＡＧ香农多样性指数ＳＨＤＩ破碎度ＦＮ

生态学含义Ｅｃｏｌｏｇｉｃａｌｍｅａｎｉｎｇ

景观中斑块总数，值越大表明景观异质性越高，同时景观破碎度也越高景观总面积与斑块数量的比值，值越大景观破碎度越低

用来反映各斑块形状的复杂程度，其值越接近１表示斑块形状越规则、集聚程度越高用来测度同一类型斑块的自然连通度，其值越高表示斑块集聚度越高表征景观类型的蔓延程度，值越高表明景观中存在某种连通性较好的斑块

反映景观的丰富度，描述不同景观面积大小的均匀程度，值越大表明各类斑块在景观中分布越均衡

表征景观被分割的程度，反映景观空间结构的复杂性，可在一定程度体现人类对景观的干扰度

　　ＮＰ：ＮｕｍｂｅｒｏｆＰａｔｃｈｅｓ；ＡＲＥＡ＿ＭＮ：ＭｅａｎｏｆＰａｔｃｈＡｒｅａ；ＳＨＡＰＥ：ＳｈａｐｅＩｎｄｅｘ；ＣＯＨＥＳＩＯＮ：ＣｏｈｅｓｉｏｎＩｎｄｅｘ；ＣＯＮＴＡＧ：ＣｏｎｔａｇｉｏｎＩｎｄｅｘ；ＳＨＤＩ：Ｓｈａｎｎｏｎ′ｓＤｉｖｅｒｓｉｔｙＩｎｄｅｘ；ＦＮ：ＦｒａｇｍｅｎｔａｔｉｏｎＩｎｄｅｘ

２．２．２　基于ＩｎＶＥＳＴ模型的生境质量评估

威胁源和生境类型，在威胁源设定时考虑到建设用地是所有土地利用类型中人类活动较为集中的地类，它会对区域生物多样性造成威胁，因此将所有建设用地设为威胁源，而耕地作为半人工半自然且经常受到人为干扰的地类，也被设定为威胁源，其余不同的地类代表了不同的生境类型。模型结合生境的敏感度和外界威胁强度来计算生境质量。首先计算生境退化度：

Ｄｘｊ＝

ｗｒ／∑ｗｒ)ｒｙｉｒｘｙβｘＳｊｒ(∑∑ｒ＝１ｙ＝１ｒ＝１

Ｒ

Ｙｒ

Ｒ

ＩｎＶＥＳＴ模型中的生境是指被物种占有并可为其提供资源和生存及繁育的空间［２６］。模型运行需要设定

æｄｘｙö

ｉｒｘｙ＝１－ç÷（线性衰减）

èｄｒｍａｘøæ－２．９９ｄｘｙö

ｉｒｘｙ＝ｅｘｐç÷（指数衰减）

ｄèøｒｍａｘ

（２）（３）（４）

式中，Ｄｘｊ为生境类型ｊ中ｘ栅格的生境退化度；Ｒ为威胁源个数；Ｗｒ为威胁源ｒ的权重；Ｙｒ为威胁源的栅格数；ｒｙ为栅格ｙ的胁迫值；ｉｒｘｙ为栅格ｙ的胁迫值ｒｙ对栅格ｘ的胁迫水平；βｘ为威胁源对栅格ｘ的可达性（受法律保护的区域为０，其余区域为１）；Ｓｊｒ为生境类型ｊ对威胁源ｒ的敏感度；ｄｘｙ为栅格ｘ与栅格ｙ的直线距离；ｄｒｍａｘ为

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

７３０４　生　态　学　报　　　３８卷　

威胁源ｒ的最大胁迫距离。

在此基础上计算生境质量：

Ｄｘｊéùú－Ｑｘｊ＝Ｈｊê（５）１（）ｚｚúê＋Ｄｘｊｋûë

式中，Ｑｘｊ为生境类型ｊ中ｘ栅格的生境质量指数；Ｈｊ为生境类型ｊ的生境适宜度，值域为［０，１］；ｋ为半饱和常

ｚ

数，取最大生境退化度（由模型运算一次得到）的一半；ｚ为归一化常量，通常设为２．５［２６］。

由于模型中相关参数需要根据研究区具体情况进行调整，因此在考虑模型使用说明基础之上，向熟悉该

区域的中国科学院寒区旱区环境与工程研究所、兰州大学、甘肃农业大学、榆中县国土资源局、榆中县环保局等单位的２５位相关领域专家发放调查问卷，根据研究区的实际情况对威胁源的胁迫距离和权重、各土地利用类型的生境适宜度以及对各类威胁源的敏感度进行打分，采用德尔菲法（Ｄｅｌｐｈｉ）综合确定各项参数最终值（表３、表４）。

表３　威胁源及其最大威胁距离、权重及衰减类型

Ｔａｂｌｅ３　Ｔｈｒｅａｔａｎｄｉｔｓｍａｘｉｍｕｍｅｆｆｅｃｔｄｉｓｔａｎｃｅ，ｗｅｉｇｈｔａｎｄｓｐａｔｉａｌｄｅｃａｙｔｙｐｅ

威胁源Ｔｈｒｅａｔｓ

工矿用地Ｉｎｄｕｓｔｒｉａｌａｎｄｍｉｎｉｎｇｌａｎｄ城镇Ｔｏｗｎ村庄Ｖｉｌｌａｇｅ耕地Ｃｕｌｔｉｖａｔｅｄｌａｎｄ铁路Ｒａｉｌｗａｙ公路Ｒｏａｄ

最大威胁距离／ｋｍ

Ｍａｘｉｍｕｍｄｉｓｔａｎｃｅｏｆｉｎｆｌｕｅｎｃｅ

４５３１２３

权重Ｗｅｉｇｈｔ０．５１１０．１５０．４０．６

空间衰减类型Ｓｐａｔｉａｌｄｅｃａｙｔｙｐｅ

指数指数指数线性线性线性

表４　生境适宜度及其对不同威胁源的相对敏感程度

Ｔａｂｌｅ４　Ｈａｂｉｔａｔｓｕｉｔａｂｉｌｉｔｙｄｅｇｒｅｅａｎｄｉｔｓｒｅｌａｔｉｖｅｓｅｎｓｉｔｉｖｉｔｙｔｏｅａｃｈｔｈｒｅａｔ

生境类型Ｈａｂｉｔａｔｔｙｐｅ

有林地Ａｒｂｏｒｆｏｒｅｓｔｌａｎｄ灌木林地Ｓｈｒｕｂｌａｎｄ其他林地Ｏｔｈｅｒｆｏｒｅｓｔｌａｎｄ天然草地Ｎａｔｕｒｅｇｒａｓｓｌａｎｄ其他草地Ｏｔｈｅｒｇｒａｓｓｌａｎｄ河流Ｒｉｖｅｒ

水库坑塘Ｒｅｓｅｒｖｏｉｒａｎｄｐｏｎｄ内陆滩涂Ｉｎｌａｎｄｂｅａｃｈ旱地Ｄｒｙｌａｎｄ水浇地Ｉｒｒｉｇａｂｌｅｌａｎｄ果园Ｏｒｃｈａｒｄ裸地Ｂａｒｅｌａｎｄ

生境适宜度Ｈａｂｉｔａｔｓｕｉｔａｂｉｌｉｔｙ

１０．７０．６０．７０．４０．８０．７０．６０．２０．３０．４０

村庄Ｖｉｌｌａｇｅ０．８０．７０．６０．７０．７０．６０．５０．７０．７０．８０．７—

工矿用地Ｉｎｄｕｓｔｒｉａｌａｎｄｍｉｎｉｎｇｌａｎｄ

０．８０．８０．５０．６０．６０．８０．７０．７０．５０．６０．５—

城镇Ｔｏｗｎ０．７０．８０．７０．７０．７０．８０．８０．８０．７０．８０．７—

铁路Ｒａｉｌｗａｙ０．５０．６０．７０．４０．４０．４０．４０．５０．４０．６０．５—

公路Ｒｏａｄ０．６０．５０．８０．４０．４０．６０．５０．６０．５０．６０．５—

耕地Ｃｕｌｔｉｖａｔｅｄｌａｎｄ０．３０．３０．２０．３０．５０．３０．５０．６０．３０．２０．３—

２．２．３　生境质量热点分析及空间自相关分析

著的高值和低值。基于ＡｒｃＧＩＳ平台，本文使用Ｇｅｔｉｓ⁃ＯｒｄＧ∗指数来刻画榆中县生境质量冷热点［２７］。该统计量通过计算一个栅格的生境质量值及其周围栅格的生境质量值的总和与所有栅格生境质量值的总和进行比较，用来分析生境质量在局部空间上的集聚情况。当Ｇ∗值显著为正时，生境质量呈高值集聚，为热点区域，当Ｇ∗值显著为负时，生境质量呈低值集聚，为冷点区域［２８］，其中９９％置信水平上的Ｇ∗值对应的区域视为为

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

研究局部区域聚类分布特征的方法称为热点分析，可用来表示生境质量在空间分布上是否存在统计上显

　２０期　　　刘春芳　等：基于土地利用变化的黄土丘陵区生境质量时空演变特征———以榆中县为例　７３０５

热点和冷点，９５％置信水平上的Ｇ∗值对应的区域为次热点和次冷点。计算公式为［２７］：

＝Ｇ∗ｉ

Ｓ

－（∑ｗｉｊ）］［ｎ∑ｊ＝１ｗ２ｉｊｊ＝１

ｎ

ｎ

∑ｊ＝１

ｎｗｉｊｘｊ－ｘ􀭰∑ｊ＝１ｗｉｊ

ｎ

２

（６）

ｎ－１

式中，ｘｊ为栅格ｊ的生境质量；ｗｉｊ为栅格ｉ与栅格ｊ的空间权重矩阵（通过ｑｕｅｅｎｓ准则计算空间权重矩阵），若ｉ与ｊ相邻，其空间权重为１，否则为０；ｘ为生境质量平均值；Ｓ为生境质量标准差；为ｎ为栅格总数。􀭰

空间自相关是指地理事物的某一属性在不同空间位置上的相关程度，分为全局自相关和局部自相关。生境质量在空间分布上会表现出一定的规律性，全局自相关可用于描述生境质量在整个区域上有无集聚效应，本文采用全局Ｍｏｒａｎ′ｓＩ指数来估计，其值域为［－１，１］，＞０为正相关且越接近１集聚程度越高，＜０为负相关，等于０表示随机分布［２９］。局部空间自相关分析可以用来衡量观测单元生境质量与周边单元生境质量的差异程度及其显著性。本文采用Ａｎｓｅｌｉｎ［３０］提出的空间关系局域指标ＬＩＳＡ（Ｌｏｃａｌｉｎｄｉｃａｔｏｒｓｏｆｓｐａｔｉａｌａｓｓｏｃｉａｔｉｏｎｓ）来识别榆中县每个栅格的生境质量与其周围栅格生境质量的局部关联和差异，揭示生境质量在空间上的集聚和离散特征，进而为土地利用优化与生境质量提升提供决策依据。通过ＬＩＳＡ集聚分析可知道生境质量空间集聚或离散的具体位置，它将生境质量的空间格局分为５种类型，即Ｈ⁃Ｈ（高高集聚，表示一个生境质量高值被高值包围）、Ｈ⁃Ｌ（高低离散，生境质量高值被低值包围）、Ｌ⁃Ｈ（低高离散，生境质量低值被高值包围）、Ｌ⁃Ｌ（低低集聚，生境质量低值被低值包围）和ＮＳ（不显著，不存在显著的空间集聚现象），其中，Ｈ⁃Ｈ和Ｌ⁃Ｌ为正相关类型，Ｈ⁃Ｌ和Ｌ⁃Ｈ为负相关类型，详细计算原理请参见文献［３０］。３　结果与分析３．１　土地利用变化分析

１９９５—２０１５年间，草地、耕地和林地一直是榆中县的主要地类，合计占土地总面积的８０％以上，是该县景

观基质的重要构成。２０年来，变化较大的地类是耕地、林地和建设用地。１９９５—２００５年，耕地减少４３５６ｈｍ２，减幅为３．３２％；林地增加４０４１ｈｍ２，增幅为１３．０２％；建设用地增加１３６１ｈｍ２，增幅为９．９８％。２００５—２０１５年，耕１８７２ｈｍ２，减幅为４．４７％和１．３７％；建设用地和林地分别增加３０５５、４５０７ｈｍ２，增幅为２９．７４％和１４．５％；园地减用变化的主要特点。

地减少１１６９ｈｍ２，建设用地增加１６９４ｈｍ２，林地增加４６６ｈｍ２。总体来看，２０年间耕地和草地分别减少５５２５、少２０１ｈｍ２，裸地和水域变化不大。可见，建设用地急剧扩张，林地逐渐增加和耕地持续减少是榆中县土地利

构建土地利用转移矩阵可以进一步解释２０年来榆中县土地利用转移情况（表５）。１９９５—２０１５年，耕地向建设用地和林地分别转移了３１５０ｈｍ２和３９４８ｈｍ２。这说明耕地是榆中县建设用地增加的重要来源，且“退

表５　１９９５—２０１５年榆中县土地利用类型转移矩阵／ｈｍ２

Ｔａｂｌｅ５　ＴｒａｎｓｉｔｉｏｎｍａｔｒｉｘｏｆｌａｎｄｕｓｅｏｆＹｕｚｈｏｎｇＣｏｕｎｔｙｆｒｏｍ１９９５ｔｏ２０１５

年份

Ｙｅａｒ

土地类型Ｌａｎｄｔｙｐｅ

耕地Ｃｕｌｔｉｖａｔｅｄｌａｎｄ园地Ｇａｒｄｅｎ林地Ｆｏｒｅｓｔ

１９９５

草地Ｇｒａｓｓｌａｎｄ水域Ｗａｔｅｒ裸地Ｂａｒｅｌａｎｄ

建设用地Ｃｏｎｓｔｒｕｃｔｉｏｎｌａｎｄ

２０１５

建设用地Ｃｏｎｓｔｒｕｃｔｉｏｎ

ｌａｎｄ３１５０２７０７７１２５６８７８１５４３８

耕地

Ｃｕｌｔｉｖａｔｅｄｌａｎｄ

１２４２２６

５６２３５１００４４７６２９０

园地Ｇａｒｄｅｎ１１１５３９０１４０００

林地Ｆｏｒｅｓｔ３９４８０３０７１５７０５１６６００

草地Ｇｒａｓｓｌａｎｄ１１７００１３４２９７２５８８２３４

水域Ｗａｔｅｒ０００２０３４２８１３０

裸地Ｂａｒｅｌａｎｄ

０００３２２２５０１２５６０

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

７３０６　生　态　学　报　　　３８卷　

耕还林”工程也取得了较好的效果。研究期内有约８７％的林地较为稳定，仅与耕地、建设用地和草地发生了较少数量的转移。草地与其他地类转换较为频繁，其中主要向耕地转移了１００４ｈｍ２，向林地转移了７０５ｈｍ２，是因为２００５年以来，榆中县实施了一系列的土地开发项目以及“生态北山”项目。水域和裸地的变动均较小。总体来看，１９９５—２０１５年榆中县各类土地均发生一定程度的变化，主要表现为耕地转向林地和建设用地，以及草地转向耕地和林地。３．２　景观格局变化分析

在类型水平上分析表明：研究期内草地和耕地的ＮＰ、ＡＲＥＡ＿ＭＮ和ＣＯＨＥＳＩＯＮ均较大，说明草地和耕地是榆中县的优势景观。草地的ＮＰ和ＳＨＡＰＥ逐渐增加而ＡＲＥＡ＿ＭＮ减小，表明破碎化程度增加；林地的ＮＰ、ＡＲＥＡ＿ＭＮ和ＳＨＡＰＥ增加而ＣＯＨＥＳＩＯＮ减小，表明林地面积增加但形状趋于复杂化，集聚程度有所下降，这是由于２０００年以来各乡镇陆续实施“退耕还林”工程的结果。建设用地的ＮＰ、ＡＲＥＡ＿ＭＮ和ＣＯＨＥＳＩＯＮ均逐渐增大，反映了２０年来建设用地呈集聚性扩张趋势。水域和园地的ＮＰ、ＡＲＥＡ＿ＭＮ和ＣＯＨＥＳＩＯＮ均较小，说明榆中县水域和园地数量较少且分散。３个时期草地、耕地和林地的ＣＯＨＥＳＩＯＮ均较高，说明这３种地类分布较为集中。

景观水平上分析结果（表６）表明：１９９５—２０１５年，榆中县ＮＰ持续增加，ＡＲＥＡ＿ＭＮ持续下降，表明景观趋于破碎化，ＦＮ指数的增加也进一步验证了这一特征；ＣＯＮＴＡＧ在研究期内的降低表明景观中斑块数目增加、集聚度下降，同时也说明这一时期榆中县主要地类的连通度有所下降，这与ＮＰ增加和ＡＲＥＡ＿ＭＮ减小是一致的；１９９５—２０１５年，ＳＨＤＩ从１．２７增加到１．３１，表明景观中各组分面积差异减小，景观异质性增强，这主要是由于耕地、草地面积减少的同时林地和建设用地面积增加，从而使得这些主要地类之间的面积差异减小导致的。

表６　景观水平上景观格局指数

Ｔａｂｌｅ６　Ｌａｎｄｓｃａｐｅｉｎｄｅｘｅｓｏｆｌａｎｄｓｃａｐｅｌｅｖｅｌ

年份Ｙｅａｒ１９９５２００５２０１５

斑块数量ＮＰ／个３４８１３５９８３６２８

平均斑块面积ＡＲＥＡ＿ＭＮ／ｈｍ２

９４．６６９２．８７９１．５８

蔓延度ＣＯＮＴＡＧ／％４１．７５４０．５４４０．１８

香农多样性指数

ＳＨＤＩ

１．２７１．２９１．３１

破碎度ＦＮ

０．０１０５６０．０１０７７０．０１０９２

３．３　生境质量时空格局分析３．３．１　生境质量时空格局变化分析

运行ＩｎＶＥＳＴ模型，得到榆中县１９９５—２０１５年生境质量空间分布，为便于比较生境质量的变化情况，在ＡｒｃＧＩＳ中使用自然断点法将其分类为高（Ｉ）、较高（ＩＩ）、中（ＩＩＩ）、较低（ＩＶ）和低（Ｖ）５个等级（图２），统计３个时期各等级的栅格比例及生境质量平均值（表７），并结合土地利用变化及景观指数变化来分析生境质量对ＬＵＣＣ的响应。

表７　各等级生境比例及生境平均值

Ｔａｂｌｅ７　Ｐｅｒｃｅｎｔａｇｅｏｆｅａｃｈｈａｂｉｔａｔｌｅｖｅｌａｎｄａｖｅｒａｇｅｖａｌｕｅｏｆｈａｂｉｔａｔｑｕａｌｉｔｙ

等级ＬｅｖｅｌＩ

分值区间Ｖａｌｕｅｉｎｔｅｒｖａｌ０．５—０．８０．３—０．５０．１—０．３０—０．１０．８—１

１９９５年

生境质量平均值Ａｖｅｒａｇｅｖａｌｕｅｏｆｈａｂｉｔａｔｑｕａｌｉｔｙ

２００５年

生境质量平均值Ａｖｅｒａｇｅｖａｌｕｅｏｆｈａｂｉｔａｔｑｕａｌｉｔｙ

２０１５年

比例Ｐｒｏｐｏｒｔｉｏｎ／％

４．６３３８．５７３９．８３８．０２８．９５

比例Ｐｒｏｐｏｒｔｉｏｎ／％１０．１７３８．２６３８．６８．３６４．６１

比例Ｐｒｏｐｏｒｔｉｏｎ／％１０．３２３７．８１３８．０４９．２８４．５６

生境质量平均值Ａｖｅｒａｇｅｖａｌｕｅｏｆｈａｂｉｔａｔｑｕａｌｉｔｙ

ＩＩ

ＩＩＩＶ

ＩＶ

０．３５２３０．３５５５０．３５３０

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

　２０期　　　刘春芳　等：基于土地利用变化的黄土丘陵区生境质量时空演变特征———以榆中县为例　７３０７

图２　１９９５—２０１５年榆中县生境质量空间分布

Ｆｉｇ．２　ＳｐａｔｉａｌｄｉｓｔｒｉｂｕｔｉｏｎｏｆｈａｂｉｔａｔｑｕａｌｉｔｙｏｆＹｕｚｈｏｎｇＣｏｕｎｔｙｆｒｏｍ１９９５ｔｏ２０１５

从空间格局来看（图２），全县生境质量整体表现为南部高、中部居中、北部低的态势。生境质量分布与土地利用类型分布大体一致，林地和天然草地分布区域的生境质量较高，裸地和旱地生境质量较低。全县以ＩＩＩ等级和ＩＶ等级生境为主，其中ＩＩＩ等级分布在中部和北部，ＩＶ等级主要分布在中部，土地类型以耕地和荒草地为主。Ｉ等级和ＩＩ等级生境主要分布在南部山区，该区域有大量的林地和天然草地，并且是兴隆山自然保护区所在地，受到人为干扰较小，生境质量普遍较高。Ｖ等级生境在北部有较大面积的分布，主要土地类型是裸地和农村居民点；在中部也有分布，主要的土地类型为居住地及工矿用地。

从时间尺度上看（表７），１９９５—２００５年，Ｉ等级变化较小，ＩＩ等级斑块数量有明显的上升，其比例从８．９５％上升到１０．１７％。主要原因是榆中县自２０００年以来大力推进“退耕还林”政策，尤其是“生态北山”等工程项目的具体实施，促使耕地、荒草地等较低等级生境向林地这一高等级生境转化。ＩＩＩ等级和ＩＶ等级斑块比例则分别下降０．３１％和１．２３％，Ｖ等级的斑块增加０．３４％，究其原因：一方面可从土地利用变化来解释，１９９５—２００５年建设用地增加１３６１ｈｍ２，其来源主要是耕地和园地。建设用地占用耕地、园地不仅会导致生境受损还会胁迫周围的生境，此外，有少部分荒草地退化为裸地也导致了生境质量下降；另一方面，从景观指数来看，这一时期草地和耕地斑块数目分别增加８个和３２个，平均斑块面积分别减小２．６ｈｍ２和１．２ｈｍ２，景观蔓延度指生境的胁迫水平，从而导致低等级生境斑块增加。２００５—２０１５年，生境变化与上一阶段相似，表现为ＩＩ、Ｖ等级的斑块增加，ＩＩＩ、ＩＶ等级生境减少，其中ＩＩ等级增加０．１５％，增幅较上一阶段有所下降。这主要是因为该时期林地面积增幅下降且平均斑块面积有轻微下降，Ｖ等级斑块增加０．９２％。这一时期建设用地的增幅较上一时期更大，从而对生境的胁迫作用也更加强烈。

总体来看，全县１９９５、２００５和２０１５年平均生境质量分别为０．３５２３、０．３５５５和０．３５３０，表明该县生境质量状况略有提升且呈先升后降的变化特征，由于组成该县景观基质的草地和耕地的比例较为稳定，全县整体生１８．０１％，增幅较大，且斑块数和平均斑块面积分别增加２３个和３．３ｈｍ２，而林地面积增幅较小，仅为１．３３％。１９９５—２０１５年，榆中县Ｉ等级斑块比例变化较小，“退耕还林”、“生态北山”等工程的实施促使林地平均斑块面积增大、集聚度增加，进而使得ＩＩ等级生境斑块增加１．３７％。ＩＩＩ、ＩＶ等级生境分别减少０．７６％和１．７９％，Ｖ

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

境变化较小。１９９５—２００５年，在“退耕还林”的推动下林地面积大幅增加，而２００５—２０１５年建设用地增加数下降１．２１％，生境斑块之间分离度增大。同时，建设用地的斑块数目增加在一定程度上也加剧了威胁源对

７３０８　生　态　学　报　　　３８卷　

等级生境增加１．２６％，究其原因，该时期内有３１５０ｈｍ２耕地转为建设用地，建设用地斑块数的增加，势必对生境斑块造成胁迫并使生境趋于破碎化。

为更加深入了解生境质量变化的时空特征，使用ＡｒｃＧＩＳ得到生境质量变化图（图３）。从图３可看出，１９９５—２００５年，榆中县生境质量上升区域主要位于北部，其中又以贡井和园子岔分布最为广泛，而生境质量下降的区域位于中部，以和平、夏官营和城关等乡镇为主。２００５—２０１５年，生境质量呈现普遍小幅下降，尤其是中部的和平、定远和金崖有较大面积的下降区域出现，而上升区域仅在北部和南部有少量分布。总体来看，１９９５—２０１５年，大部分地区生境质量无明显变化，生境质量上升区域主要分布于北部，主要原因是“退耕还林”以及“生态北山”等工程的实施。生境质量下降区域主要位于中部，特别是定远、和平、连塔这几个乡镇分布较多，原因在于兰州市近年来实施城市东扩，居住地、工矿用地以及交通运输用地等建设用地快速增加，进而胁迫周围生境，并导致生境破碎度加剧、连通性变差，最终使得这些区域的生境质量下降。

图３　榆中县生境质量变化空间分布

Ｆｉｇ．３　ＨａｂｉｔａｔｑｕａｌｉｔｙｃｈａｎｇｅｍａｐｏｆＹｕｚｈｏｎｇＣｏｕｎｔｙ

３．３．２　生境质量空间热点分析

从图４看出，榆中县生境质量空间冷热点分布呈现出“南热北冷”的格局。热点区域主要位于南部，包括马坡、新营和小康营３个乡镇，这一区域雨水充沛，植被覆盖度高，且是兴隆山保护区所在地，生境质量普遍较高（平均生境质量０．６以上）；冷点主要分布于北部黄土丘陵区，以园子岔、上花岔等乡镇为主，这一区域常年干旱少雨，植被稀疏，地形破碎，主要土地类型为荒草地和旱地，因而生境质量普遍较低（平均生境质量０．２５左右）；中部有次冷点分布，这些区域主要分布有大量的居民点、道路等建设用地，由于胁迫效应导致区内生境质量较低。从时间变化来看，１９９５—２００５年，中部次冷点范围有所下降，而北部次冷点范围增加，热点区域变化不大；２００５—２０１５年，中部和北部部分次冷点转变为不显著，南部山区与平原区交界处有次冷点出现，西南部有小部分热点转为次热点。３．３．３　生境质量空间自相关分析

由全局空间自相关分析可知，榆中县１９９５、２００５和２０１５年生境质量的全局Ｍｏｒａｎ′ｓＩ值分别为０．１６０、０．１５６和０．１７２（Ｐ＜０．０１），表明该县生境质量分布具有一定空间正相关性，在空间上存在集聚现象，且２０１５年Ｍｏｒａｎ′ｓＩ的增大表明生境质量的空间集聚性有所增强。

根据生境质量ＬＩＳＡ集聚分析得到图５，从中可看出在９５％的置信水平上，榆中县生境质量局部空间关联

类型特征为：大部分区域表现为ＮＳ（不显著），“Ｈ⁃Ｈ”集聚区域主要分布于南部的兴隆山和马啣山，表明该区

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

　２０期　　　刘春芳　等：基于土地利用变化的黄土丘陵区生境质量时空演变特征———以榆中县为例　７３０９

图４　榆中县生境质量热点分布

Ｆｉｇ．４　ＨｏｔｓｐｏｔｓｄｉｓｔｒｉｂｕｔｉｏｎｏｆｈａｂｉｔａｔｑｕａｌｉｔｙｏｆＹｕｚｈｏｎｇＣｏｕｎｔｙ

图５　榆中县生境质量ＬＩＳＡ集聚图

Ｆｉｇ．５　ＬＩＳＡｃｌｕｓｔｅｒｉｎｇｍａｐｏｆｈａｂｉｔａｔｑｕａｌｉｔｙｏｆＹｕｚｈｏｎｇＣｏｕｎｔｙ

域总体生境质量较好，可定义为生态涵养区。１９９５—２０１５年，“Ｈ⁃Ｈ”集聚区增加０．６８％，并向西扩展，体现了近年来生态保护的效果。“Ｈ⁃Ｌ”和“Ｌ⁃Ｈ”离散区未出现明显的集中分布，主要位于中部川区靠近兴隆山周边的区域，且在１９９５—２０１５年间呈逐渐增加趋势。这一方面体现了城镇扩张对生境的威胁，另一方面也说明未来需要对此类区域采取更为严格的生态治理与保护。“Ｌ⁃Ｌ”集聚区主要位于北部园子岔、上花岔和哈岘等乡镇，１９９５—２０１５年间上升了０．６７％，表明这一区域生境质量有所下降，未来应加强这一区域的退耕还林还草及生态保育。

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

７３１０　生　态　学　报　　　３８卷　

４　结论与讨论

本文以典型黄土丘陵区———甘肃省榆中县为例，采用ＩｎＶＥＳＴ模型分析了该县生境质量的时空分异特征，并探讨了生境质量对土地利用变化及景观格局变化的响应。论文将生境质量的研究拓展到了县域尺度，丰富了黄土丘陵区生境质量时空变化的研究，对该区域生态安全构建与格局优化也有着积极的政策价值。

分析表明：１９９５—２０１５年，榆中县土地利用变化呈现出建设用地急剧增加，耕地快速减少，林地逐渐增加的特点。主要的土地转变类型是耕地转变为建设用地和林地，草地转变为耕地和林地。上述土地利用变化导致景观斑块数量增大，平均斑块面积降低，景观破碎度上升。总体来看，全县生境质量稳中略有提升，但县域内部生境质量由南向北逐渐递减，呈现出“南热北冷”的空间格局。研究期内，生境质量升高区域主要位于北部，原因是近年来实施的“退耕还林”工程和“生态北山”工程；生境质量下降区域主要位于中部，原因是城镇化的推进导致区内耕地、园地及水域面积的下降。２０年来，全县生境质量的空间集聚性有所增强，表现为南部“Ｈ⁃Ｈ”集聚区有所增加，北部“Ｌ⁃Ｌ”集聚区有所增加，而随着城镇化的推进，中部“Ｌ⁃Ｈ”离散区有所增加。论文对县域尺度生境质量变化的时空特征进行了分析，但对生境威胁源之间的相互作用关系以及它们对生境胁迫的内在机制方面还有待进一步的研究。

区域生境质量直接影响居民的生计与福祉。随着社会经济的发展和人口的增长，人地矛盾不断增大，人为干扰对生境的影响越来越大。尽管部分乡镇由于“退耕还林”和“生态北山”等工程的实施，区内生境质量有明显改善。但快速城镇化使得该区域大量耕地被建设用地占用，导致生境的丧失和破碎化，进而使得该区域低生境质量区域有所增加，并且与２００５年相比，２０１５年山区与平原交界处生境存在降低的风险。因此，需要针对区域内部不同地域的生境特征与差异性制定相应的调整与优化策略。

区域生态建设要加强对山区的连片植树造林及已有林地的保育，同时还应重点关注山区和平原交界地区，可通过限制浅山区开发及建立生态缓冲区来保护这一区域的生境。就榆中县生态建设与保护而言，未来要加强对北部地区的生态整治，南部兴隆山可划定生态保护红线，以减小城镇扩张对该区域生境质量的威胁。由于不同土地类型提供物种持续发展的能力有差异，土地的空间布局也会影响区域生境质量，今后的土地规划和生态保护工作要统筹安排与合理布局各类土地，重点保护森林、草地等生态源地，积极推进退耕还林还草，降低建设用地扩张对生态源地和生态廊道的破坏，提高生境的斑块面积及斑块间的连通度，进而提高黄土丘陵区区域整体生境质量。

参考文献（Ｒｅｆｅｒｅｎｃｅｓ）：

［１］　ＭｃＫｉｎｎｅｙＭＬ．Ｕｒｂａｎｉｚａｔｉｏｎ，Ｂｉｏｄｉｖｅｒｓｉｔｙ，ａｎｄＣｏｎｓｅｒｖａｔｉｏｎ：ｔｈｅｉｍｐａｃｔｓｏｆｕｒｂａｎｉｚａｔｉｏｎｏｎｎａｔｉｖｅｓｐｅｃｉｅｓａｒｅｐｏｏｒｌｙｓｔｕｄｉｅｄ，ｂｕｔｅｄｕｃａｔｉｎｇａ

８８３⁃８９０．

ｈｉｇｈｌｙｕｒｂａｎｉｚｅｄｈｕｍａｎｐｏｐｕｌａｔｉｏｎａｂｏｕｔｔｈｅｓｅｉｍｐａｃｔｓｃａｎｇｒｅａｔｌｙｉｍｐｒｏｖｅｓｐｅｃｉｅｓｃｏｎｓｅｒｖａｔｉｏｎｉｎａｌｌｅｃｏｓｙｓｔｅｍｓ．ＢｉｏＳｃｉｅｎｃｅ，２００２，５２（１０）：

［２］　ＦｏｒｍａｎＲＴＴ，ＳｐｅｒｌｉｎｇＤ，ＢｉｓｓｏｎｅｔｔｅＪＡ，ＣｌｅｖｅｎｇｅｒＡＰ，ＣｕｔｓｈａｌｌＣＤ，ＤａｌｅＶＨ，ＦａｈｒｉｇＬ，ＦｒａｎｃｅＲＬ，ＧｏｌｄｍａｎＣＲ，ＨｅａｎｕｅＫ，ＪｏｎｅｓＪ，［３］　王建华，田景汉，吕宪国．挠力河流域河流生境质量评价．生态学报，２０１０，３０（２）：４８１⁃４８６．

ＳｗａｎｓｏｎＦ，ＴｕｒｒｅｎｔｉｎｅＴ，ＷｉｎｔｅｒＴＣ．ＲｏａｄＥｃｏｌｏｇｙ：ＳｃｉｅｎｃｅａｎｄＳｏｌｕｔｉｏｎｓ．Ｗａｓｈｉｎｇｔｏｎ，ＤＣ：ＩｓｌａｎｄＰｒｅｓｓ，２００３．

［４］　严岩，朱捷缘，吴钢，詹云军．生态系统服务需求、供给和消费研究进展．生态学报，２０１７，３７（８）：２４８９⁃２４９６．［５］　于丹丹，吕楠，傅伯杰．生物多样性与生态系统服务评估指标与方法．生态学报，２０１７，３７（２）：３４９⁃３５７．

［６］　ＳａｌａＯＥ，ＣｈａｐｉｎＩＩＩＦＳ，ＡｒｍｅｓｔｏＪＪ，ＢｅｒｌｏｗＥ，ＢｌｏｏｍｆｉｅｌｄＪ，ＤｉｒｚｏＲ，Ｈｕｂｅｒ－ＳａｎｗａｌｄＥ，ＨｕｅｎｎｅｋｅＬＦ，ＪａｃｋｓｏｎＲＢ，ＫｉｎｚｉｇＡ，Ｌｅｅｍａｎｓ

２１００．Ｓｃｉｅｎｃｅ，２０００，２８７（５４５９）：１７７０⁃１７７４．

Ｒ，ＬｏｄｇｅＤＭ，ＭｏｏｎｅｙＨＡ，ＯｅｓｔｅｒｈｅｌｄＭ，ＰｏｆｆＮＬ，ＳｙｋｅｓＭＴ，ＷａｌｋｅｒＢＨ，ＷａｌｋｅｒＭ，ＷａｌｌＤＨ．Ｇｌｏｂａｌｂｉｏｄｉｖｅｒｓｉｔｙｓｃｅｎａｒｉｏｓｆｏｒｔｈｅｙｅａｒ

［７］　王大尚，郑华，欧阳志云．生态系统服务供给、消费与人类福祉的关系．应用生态学报，２０１３，２４（６）：１７４７⁃１７５３．

［８］　ＧｏｅｒｔｚＪＷ．Ｔｈｅｉｎｆｌｕｅｎｃｅｏｆｈａｂｉｔａｔｑｕａｌｉｔｙｕｐｏｎｄｅｎｓｉｔｙｏｆｃｏｔｔｏｎｒａｔｐｏｐｕｌａｔｉｏｎｓ．ＥｃｏｌｏｇｉｃａｌＭｏｎｏｇｒａｐｈｓ，１９６４，３４（４）：３５９⁃３８１．

［９］　ＴｈｉｅｌＲＰ．ＲｅｌａｔｉｏｎｓｈｉｐｂｅｔｗｅｅｎｒｏａｄｄｅｎｓｉｔｉｅｓａｎｄｗｏｌｆｈａｂｉｔａｔｓｕｉｔａｂｉｌｉｔｙｉｎＷｉｓｃｏｎｓｉｎ．ＡｍｅｒｉｃａｎＭｉｄｌａｎｄＮａｔｕｒａｌｉｓｔ，１９８５，１１３（２）：４０４⁃４０７．

Ｐｏｌｉｃｙ，１９９９，２（６）：４６５⁃４７５．

［１０］　ＥｎｇｅｌＤＷ，ＴｈａｙｅｒＧＷ，ＥｖａｎｓＤＷ．Ｌｉｎｋａｇｅｓｂｅｔｗｅｅｎｆｉｓｈｅｒｙｈａｂｉｔａｔｑｕａｌｉｔｙ，ｓｔｒｅｓｓｏｒｓ，ａｎｄｆｉｓｈｅｒｙｐｏｐｕｌａｔｉｏｎｓ．ＥｎｖｉｒｏｎｍｅｎｔａｌＳｃｉｅｎｃｅ＆

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

　２０期　　　刘春芳　等：基于土地利用变化的黄土丘陵区生境质量时空演变特征———以榆中县为例　７３１１

［１１］　ＤｏｎａｌｄＰＦ，ＧｒｅｅｎＲＥ，ＨｅａｔｈＭＦ．ＡｇｒｉｃｕｌｔｕｒａｌｉｎｔｅｎｓｉｆｉｃａｔｉｏｎａｎｄｔｈｅｃｏｌｌａｐｓｅｏｆＥｕｒｏｐｅ′ｓｆａｒｍｌａｎｄｂｉｒｄｐｏｐｕｌａｔｉｏｎｓ．ＰｒｏｃｅｅｄｉｎｇｓｏｆｔｈｅＲｏｙａｌ［１２］　ＪｅｔｚＷ，ＷｉｌｃｏｖｅＤＳ，ＤｏｂｓｏｎＡＰ．Ｐｒｏｊｅｃｔｅｄｉｍｐａｃｔｓｏｆｃｌｉｍａｔｅａｎｄｌａｎｄ－ｕｓｅｃｈａｎｇｅｏｎｔｈｅｇｌｏｂａｌｄｉｖｅｒｓｉｔｙｏｆｂｉｒｄｓ．ＰＬｏＳＢｉｏｌｏｇｙ，２００７，５（６）：［１３］　ＯｔｔｏＣＲＶ，ＲｏｔｈＣＬ，ＣａｒｌｓｏｎＢＬ，ＳｍａｒｔＭＤ．Ｌａｎｄ－ｕｓｅｃｈａｎｇｅｒｅｄｕｃｅｓｈａｂｉｔａｔｓｕｉｔａｂｉｌｉｔｙｆｏｒｓｕｐｐｏｒｔｉｎｇｍａｎａｇｅｄｈｏｎｅｙｂｅｅｃｏｌｏｎｉｅｓｉｎｔｈｅ［１４］　ＳｉｅｂｅｒＡ，ＵｖａｒｏｖＮＶ，ＢａｓｋｉｎＬＭ，ＲａｄｅｌｏｆｆＶＣ，ＢａｔｅｍａｎＢＬ，ＰａｎｋｏｖＡＢ，ＫｕｅｍｍｅｒｌｅＴ．Ｐｏｓｔ－Ｓｏｖｉｅｔｌａｎｄ－ｕｓｅｃｈａｎｇｅｅｆｆｅｃｔｓｏｎｌａｒｇｅ［１５］　盛学斌，刘云霞，孙建中．近５０年冀北高原土地利用变化的土壤生境效应．应用生态学报，２００４，１５（４）：５８９⁃５９２．

ｍａｍｍａｌｓ′ｈａｂｉｔａｔｉｎＥｕｒｏｐｅａｎＲｕｓｓｉａ．ＢｉｏｌｏｇｉｃａｌＣｏｎｓｅｒｖａｔｉｏｎ，２０１５，１９１：５６７⁃５７６．

ｎｏｒｔｈｅｒｎＧｒｅａｔＰｌａｉｎｓ．ＰｒｏｃｅｅｄｉｎｇｓｏｆｔｈｅＮａｔｉｏｎａｌＡｃａｄｅｍｙｏｆＳｃｉｅｎｃｅｓｏｆｔｈｅＵｎｉｔｅｄＳｔａｔｅｓｏｆＡｍｅｒｉｃａ，２０１６，１１３（３７）：１０４３０⁃１０４３５．ｅ１５７．

ＳｏｃｉｅｔｙＢＢｉｏｌｏｇｉｃａｌＳｃｉｅｎｃｅｓ，２００１，２６８（１４６２）：２５⁃２９．

［１６］　吴未，张敏，许丽萍，欧名豪．土地利用变化对生境网络的影响———以苏锡常地区白鹭为例．生态学报，２０１５，３５（１４）：４８９７⁃４９０６．［１８］　刘春芳，王川，刘立程．三大自然区过渡带生境质量时空差异及形成机制———以榆中县为例．地理研究，２０１８，３７（２）：４１９⁃４３２．［１９］　刘红玉，李兆富，白云芳．挠力河流域东方白鹳生境质量变化景观模拟．生态学报，２００６，２６（１２）：４００７⁃４０１３．

２０１５，３２（３）：６２２⁃６２９．

［１７］　吴健生，曹祺文，石淑芹，黄秀兰，卢志强．基于土地利用变化的京津冀生境质量时空演变．应用生态学报，２０１５，２６（１１）：３４５７⁃３４６６．

［２０］　包玉斌，刘康，李婷，胡胜．基于ＩｎＶＥＳＴ模型的土地利用变化对生境的影响———以陕西省黄河湿地自然保护区为例．干旱区研究，［２１］　杨芝歌，周彬，余新晓，樊登星，张艺，王奇．北京山区生物多样性分析与碳储量评估．水土保持通报，２０１２，３２（３）：４２⁃４６．

１０５⁃１１４．

［２２］　杨建强，朱永贵，宋文鹏，张娟，张龙军，罗先香．基于生境质量和生态响应的莱州湾生态环境质量评价．生态学报，２０１４，３４（１）：［２３］　刘世梁，尹艺洁，杨珏婕，安南南，王聪，董世魁．漫湾库区景观破碎化对区域生境质量的影响．生态学报，２０１７，３７（２）：６１９⁃６２７．［２４］　刘春芳，谢国林，李锋．新常态下黄土丘陵区土地整治规划探索与实践．北京：科学出版社，２０１７．［２５］　郑新奇，付梅臣．景观格局空间分析技术及其应用．北京：科学出版社，２０１０．

［２６］　ＴａｌｌｉｓＨ，ＲｉｃｋｅｔｔｓＴ，ＧｕｅｒｒｙＡ，ＷｏｏｄＳ，ＳｈａｒｐＲ，Ｃｈａｐｌｉｎ－ＫｒａｍｅｒＲ．ＩｎＶＥＳＴ３．３．０Ｕｓｅｒ′ｓＧｕｉｄｅ：ＩｎｔｅｇｒａｔｅｄＶａｌｕａｔｉｏｎｏｆＥｎｖｉｒｏｎｍｅｎｔａｌ［２７］　ＧｅｔｉｓＡ，ＯｒｄＪＫ．Ｔｈｅａｎａｌｙｓｉｓｏｆｓｐａｔｉａｌａｓｓｏｃｉａｔｉｏｎｂｙｕｓｅｏｆｄｉｓｔａｎｃｅｓｔａｔｉｓｔｉｃｓ．ＧｅｏｇｒａｐｈｉｃａｌＡｎａｌｙｓｉｓ，１９９２，２４（３）：１８９⁃２０６．［２９］　ＭｏｒａｎＰＡＰ．Ｎｏｔｅｓｏｎｃｏｎｔｉｎｕｏｕｓｓｔｏｃｈａｓｔｉｃｐｈｅｎｏｍｅｎａ．Ｂｉｏｍｅｔｒｉｋａ，１９５０，３７（１／２）：１７⁃２３．

ＳｅｒｖｉｃｅｓａｎｄＴｒａｄｅｏｆｆｓ．２０１５．［２０１７⁃１０⁃１０］．ｈｔｔｐ：／／ｄａｔａ．ｎａｔｕｒａｌｃａｐｉｔａｌｐｒｏｊｅｃｔ．ｏｒｇ／ｉｎｖｅｓｔ－ｒｅｌｅａｓｅｓ／ｄｏｃｕｍｅｎｔａｔｉｏｎ／３＿０＿０／．

［２８］　ＯｒｄＪＫ，ＧｅｔｉｓＡ．Ｌｏｃａｌｓｐａｔｉａｌａｕｔｏｃｏｒｒｅｌａｔｉｏｎｓｔａｔｉｓｔｉｃｓ：ｄｉｓｔｒｉｂｕｔｉｏｎａｌｉｓｓｕｅｓａｎｄａｎａｐｐｌｉｃａｔｉｏｎ．ＧｅｏｇｒａｐｈｉｃａｌＡｎａｌｙｓｉｓ，１９９５，２７（４）：２８６⁃３０６．［３０］　ＡｎｓｅｌｉｎＬ．Ｌｏｃａｌｉｎｄｉｃａｔｏｒｓｏｆｓｐａｔｉａｌａｓｓｏｃｉａｔｉｏｎ—ＬＩＳＡ．ＧｅｏｇｒａｐｈｉｃａｌＡｎａｌｙｓｉｓ，１９９５，２７（２）：９３⁃１１５．

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ