Analysis of Land Use Change and Habitat Quality Evolution Based on InVEST and PLUS Models: Example from Hanzhong Basin
-
摘要:
汉中盆地是中国南水北调工程中的重要源头区，分析区内土地利用及生境质量，能提高对该地区演变规律的认识，进一步保护下游研究区的自然资源和生态环境安全。笔者基于2000~2020年20年内共3期土地利用数据，建立了研究区多类型土地利用数据库和生境质量评估体系，从土地利用转移矩阵、景观格局指数、生境质量等多个方面对研究区内的土地利用的时空变化进行研究，并利用PLUS模型进行土地利用扩张和预测分析。结果表明：2000~2020年，土地利用类型变化较大，变化转移主要发生在水田和农村居民点用地、城镇用地之间；蔓延度指数逐年减少，景观破碎化程度较高；生境质量结果呈逐年下降趋势，高值面积减少并向南部扩展，低值面积增多并以城镇为中心向外辐射扩散，高等级的生态质量景观由水域和林草为主，主要沿汉江流域分布。
Abstract:Hanzhong Basin is an important source area in China's South-to-North Water Diversion Project. Analyzing land use and habitat quality within the region can enhance our understanding of its evolutionary patterns and further protect the natural resources and ecological environment security of the downstream working areas. This paper establishes a multi-type land use database and a habitat quality assessment system for the working area based on three periods of land use data from 2000 to 2020. It studies the spatiotemporal changes of land use in the research area from multiple aspects including land use transition matrix, landscape pattern indices, and habitat quality, and uses the PLUS model for land use expansion and prediction analysis. The results show that from 2000 to 2020, there were significant changes in land use types, mainly transfers between paddy fields, rural residential land, and urban land. The sprawl index decreased year by year, indicating a high degree of landscape fragmentation. Habitat quality results showed a downward trend year by year, with the area of high values decreasing and expanding southward, and the area of low values increasing, radiating outward from urban centers. High-grade ecological quality landscapes, dominated by water bodies and forests and grasslands, are mainly distributed along the Han river basin.
-
Keywords:
- InVEST model /
- PLUS model /
- land use /
- transfer matrix /
- landscape pattern /
- simulated prediction
-
中国生态地质调查起步于20世纪80年代，近年来国内学者在相关研究方面成果丰硕（王尧等，2019；张茂省等，2019；王京彬等，2020；聂洪峰等，2021；李文明等，2022）。南水北调中线工程是中国重要的跨流域调水工程，陕西省商洛市丹江源地区是南水北调中线工程重要水源涵养区（图1a），因地处秦岭保护区，又是全国生态环境建设重点试点示范区，具有重要的生态功能地位。同时，土地利用/覆被变化是导致生态格局演变和生态系统变化的重要因素（Veldkamp et al.，2004；岳德鹏等，2007；陈利顶等，2008，2013），景观指数可以体现景观结构组成和空间配置等特征，能高度浓缩景观格局信息，定量获取要素的空间分布特征，为生态格局动态变化研究提供基础信息（Robert，1997；邬建国，2000；Moser et al.，2002；陈文波等，2002；王天山等，2016）。张雁等（2016，2018）对商洛市水源地土地利用类型结构及生态环境效应进行分析，但多侧重于水生态评价研究。李文明等（2022）梳理了西北地区的生态地质调查进展并分析了生态安全屏障区存在的问题，尚未对研究区生态地质格局进行过深入调查。因此，摸清丹江源地区的生态地质本底，分析生态格局动态演变，以及科学合理利用配置土地资源意义重大。笔者依托丹江源地区生态地质调查项目，采用坡面调查方法，梳理不同地质建造上的典型生态关键带生态地质特征，利用ENVI 5.3、ArcGIS 10.2及Fragstats4.2软件，解译分析区内2000～2020年3期土地利用时空演变和景观格局演变特征，掌握全区地下—地表生态地质格局动态演变，为丹江源地区的生态地质环境治理、土地利用结构优化、生态系统保护及生态质量评价提供重要的数据支撑和科学依据，对充分发挥其生态服务功能具有重要现实意义。
![]()
图 1 研究区地理位置（a）及丹江源地区水系分布图（b）Figure 1. (a) Geographical location of the study area and (b) drainage distribution map of Danjiangyuan Area1. 生态地质背景
按照流域划分原则，在ArcGIS 10.2软件支持下，界定商洛市商州区内丹江及支流的流域面积为1 757 km2（图1a）。研究区地形总体呈西北高、东南低，海拔为699～1 684 m，主要有中山区、低山丘陵区、河谷川塬区3种地貌单元，具有暖温带南缘过渡带季风性、半湿润性山地气候特点。区内多年平均降水量为699.44 mm，年际间降水量分布不均，受地形地貌影响，区域降水呈西多东少，南多北少的地域特征，总体表现为山地多，河谷阶地少。山区年均降雨量800 mm以上，年均气温为10～11 ℃，河谷川塬区年均降雨量在730 mm以下，年均气温在13 ℃以上（国家气象科学数据中心，1951～2020）。降水量随着高度增加而增多，表现为降水量的垂直差异性。区内水系发育，河网密度0.69～1.28 km/km2（图1b），多年平均年径流深为230 mm，年径流量为4.48×108 m3，平均流量为8 m3/s，径流年内分配不均匀。
研究区大地构造位置属华北板块，以商丹断裂为界，划分为华北陆块和商丹地壳对接带2个二级构造单元（图2）。基于时代+物质组成的划分方案，将区内地质建造划分为13类（图2a），前寒武纪砂泥质碎屑岩建造、侏罗纪—古近纪陆相砂泥质碎屑岩建造、前寒武纪酸性岩浆岩建造、寒武纪—奥陶纪海相碳酸盐岩建造等为主（图2b）。研究区内岩浆岩较发育，火山岩岩性为变玄武岩、变安山岩、凝灰岩等。侵入岩主要有新元古代混合片麻岩、志留纪辉长岩体等，另有少量三叠纪、侏罗纪花岗岩体。
![]()
图 2 丹江源地区地质建造图（a）及不同地质建造面积分布图（b）1.第四纪陆相松散堆积建造；2.寒武纪海相碳酸盐岩建造；3.古生代中酸性岩浆岩建造；4.寒武纪基性火山岩建造；5.侏罗纪-古近纪陆相砂泥质碎屑岩建造；6.奥陶纪泥质碎屑岩建造；7.前寒武纪碳酸盐岩建造；8.三叠纪陆相砂泥质碎屑岩建造；9.古生代基性火山岩建造；10.前寒武纪砂泥质碎屑岩建造；11.三叠纪中酸性岩浆岩建造；12.古生代基性-超基性岩浆岩建造；13.前寒武纪酸性岩浆岩建造；14.实测剖面Figure 2. (a) Geological construction map and (b) distribution map of different geological construction areas in Danjiangyuan area区内土壤类型有黄棕壤、紫色土、褐土、潮土、粗骨土和新积土（潘贤章等，2015）。丹江源地区植被空间分布明显，丹江源地区用材林和主要水源涵养林区主要分布在中山区（陕西省商洛地区地理志，1981）；土壤较瘠簿的低山区主要分布有常绿阔叶树的落叶林带、栓皮栎林带、经济林、用材林和草场草坡；低山丘陵区森林覆盖率较低，以天然草地、人工林为主。区内产业结构以农业、林业为主导产业，经济基础较为薄弱。
2. 研究方法
2.1 生态地质剖面与生态关键带解析
在初步查明生态地质背景的基础上，结合微地貌的影响因素，采用野外路线调查结合坡面调查的方法（王京彬等，2020；聂洪峰等，2021），选取代表性的地段开展生态地质剖面测量解析，控制重要的生态地质类型，结合浅钻及垂向剖面等形式予以揭露，观察、测量生态地质现象，记录不同地质、地形地貌、生态、土壤等生态地质信息，绘制生态地质剖面，梳理区域生态地质条件现状及制约因素，垂向上查明不同地质建造上的典型生态关键带相关生态要素间的相互作用过程。
2.2 遥感解译
文中遥感影像数据选择2000年Landsat TM5，2010年Landsat TM5，2020年Landsat OLI影像，时相以9～10月为主，分辨率为30 m，云量均3%以下，使用的数字高程模型下载自地理空间数据云官网（中国科学院计算机网络信息中心科学数据中心，http://www.gscloud.cn）。运用ENVI 5.3对3期遥感航片进行解译。结合商州区土地利用现状，对研究区土地利用类型进行划分，并与GLOBELLAND30数据（Zhu et al.，2010；Chen et al.，2011，2014，2016）进行比对，准确率达到84.6%以上，解译结果符合精度要求，在ArcGIS 10.2软件支持下，绘制3期土地利用分类图，分析区域土地利用时空演变。
2.3 景观格局分析及指数选取
结合研究实际需求，笔者参考前人研究成果的基础上（李煜东等，2020；张晓宁等，2020；苏明伟等，2021；张建等，2021；张林等，2022），选取相关景观格局指数进行计算（表1），主要从斑块类型和景观水平层次开展研究。利用 ArcGIS 10.2的空间分析模块将土地利用类型图转换为Grid格式后，采用目前广泛使用的景观格局计算分析软件 Frastats 4.2进行丹江源地区2000～2020年3期各景观组分的相关景观指数计算。
表 1 景观格局指数选取及指示意义Table 1. Selection and indicative significance of landscape pattern index景观单元特征指数 指标选取 指示意义 景观面积度量指标 斑块类型面积（CA） CA值的大小影响着斑块类型聚集地中的物种数量及丰度 斑块类型比（PLAND） PLAND指某一斑块类型占整个景观面积的相对比例 最大斑块指数（LPI） LPI主要表示某一景观类型最大斑块占整个景观面积的比例，决定了景观优劣斑块，反应景观变化受人类活动干扰程度 景观面积（TA） TA定义景观幅度，是监测生态系统是否稳定的重要指标 景观形状指标 景观形状指数（LSI） LSI反映景观和斑块形状的分散和规则程度，值越大说明景观形状越复杂 景观邻近度指标 香农多样性指数（SHDI） SHDI表示景观类型的复杂程度，值越小斑块类型越少，值增大说明斑块类型增加或各斑块类型在景观中呈均衡化趋势分布 景观聚集与分散度
测量指标蔓延度指数（CONTAG） CONTAG反映景观类型的聚集程度和延展程度，高蔓延度值表明某种优势斑块类型具有良好的连通性 散布与并列指数（IJI） IJI反映斑块类型的隔离分布情况 聚集度指数（AI） AI是基于栅格数量来测度景观或者某种斑块类型的聚集程度 注：各景观格局指数计算公式和详细意义可参阅相关文献（何鹏等，2009；孙天成等，2019）。 3. 结果与分析
3.1 典型生态地质剖面与生态关键带解剖
基于生态关键带基岩–成土母质–土壤–植被垂向调查思路（图3），笔者绘制区内基岩–成土母质–土壤–植被变化图（图4）和典型生态地质剖面（图5）。
![]()
图 3 丹江源地区不同生态关键带基岩–成土母质–土壤–植被变化图Figure 3. Changes of bedrock, parent material, soil and vegetation in different ecological key zones of Danjiangyuan area![]()
图 4 丹江源地区金陵寺镇–柿园街生态地质剖面Figure 4. Eco–geological profile of Jinlingsi Town–Shizyuan Street in Danjiangyuan area![]()
图 5 丹江源地区杨峪河镇麻街岭–任家村生态地质剖面Figure 5. Eco–geological profile of Majeiling–Renjia village in Yangyuhe town, Danjiangyuan area金陵寺镇–柿园街生态地质剖面显示（图2a、图4），研究区南部的三叠纪中酸性侵入岩建造、古生代酸性侵入岩建造、前寒武纪酸性侵入岩建造的基岩节理发育，成土母质为中酸性岩类风化物，中酸性花岗岩的成土母质层和淀积层相对较薄（图3），在中山区形成的土壤类型为棕壤，质地为沙质土，结构疏松，透水性好，含水性、保水性差，上覆植被类型为针阔叶混交林，因土壤pH值呈酸性，针叶油松较发育（图3、图4）。在低山丘陵区主要形成粗骨土，砂砾含量高，植被郁闭度稍低，植被类型为以阔叶为主的针阔混交林。南部的古生代基性火山岩、基性-超基性侵入岩建造和前寒武纪基性火山岩建造的基岩易风化解体，形成较厚的风化层，质地一般为壤质，保水性能强，呈中性或弱碱性，在中山区形成棕壤，在低山丘陵区形成淋溶褐土。
杨峪河镇麻街岭–任家村生态地质剖面显示（图2a、图5），区内中部的白垩纪—古近纪陆相砂泥质碎屑岩建造，其地层产状平缓，节理不发育或仅发育风化节理等，成土母质为白垩纪–古近纪砂泥岩类风化物，例如含砾砂岩风化形成的淀积层较厚（图3），在低山丘陵区和中山区均形成紫色土，土体浅薄，质地为沙质土，有机质含量低，保水性差，水土流失易发，植被覆盖类型主要为草本植物（图3、图5）。研究区中部的前寒武纪碎屑岩建造的片理化和节理发育，泥质含量高，易风化，质地一般为壤质，在中山区形成棕壤，在低山丘陵区形成淋溶褐土。上覆植被郁闭度高，林型为以栓皮栎等阔叶林为主的针阔混交林（图3、图5）。前寒武纪碳酸盐岩建造的岩石质地坚硬，不易风化，例如大理岩的成土母质层较厚（图3），在中山区和低山丘陵区形成淋溶褐土，在河谷川塬区形成石质土，多处出现裸露基岩，植被覆盖以草本为主，零星分布少量低矮灌木和油松（图5）。
不同的地质建造，因为岩性不同、形成时代不同，受控于不同的构造地质背景和环境背景，导致形成不同的地形地貌、成土母质、土壤类型、水文类型、局部小气候类型，在以上立地条件的综合影响下，最终形成不同的植被覆盖类型和景观格局，以期为地方政府实施国土空间分区管控提供地学依据。
3.2 土地利用时空演变分析
3.2.1 土地利用动态变化
借助ENVI 5.3提取研究区的土地利用类型、分布、规模等信息。依据《土地利用现状分类》标准，将研究区土地类型分为林地、草地、耕地、水域、建设用地及未利用土地6类。通过2000、2010、2020年3期土地利用类型的遥感解译，结合研究区年际土地变更成果数据进行土地利用类型监测，研究区土地利用现状整体表现为不同土地利用类型面积随时间变化存在差异。
从土地利用类型数量结构（图6）来看，丹江源地区位于秦巴山地，气候温暖湿润，降水相对丰沛，是中国南水北调中线工程重要水源地，也是重要的生物多样性保护和水源涵养生态功能区，加之多山少田的地理特征，土地利用结构中林地比例占绝对优势；其次是耕地，林地和耕地占比约85%以上，其他用地类型比例偏低。研究区内的自然景观和人造景观之间转入转出频繁，自然景观转出面积大小排序为：草地＞林地，林地基数大，变化不显著。人造景观耕地转出明显，建设用地、水域大幅转入，转入面积大小排序为：建设用地＞水域，因研究区未利用土地多为利用度不高的裸地，其转入转出幅度变化不大。
![]()
图 6 丹江源地区2000、2010、2020年不同土地利用类型数量结构Figure 6. Quantitative structure of different land use types of Danjiangyuan area in 2000, 2010 and 2020土地利用动态度差异较明显（表2），建筑用地和水域变化最明显。建设用地动态表现为持续扩张趋势，2000～2020年建设用地动态增长度为11.70%，其中2000～2010年动态增长度为3.98%；2010～2020年增长速度较快，动态增长度为13.89%。这表明建设用地扩张主要发生在2010～2020年。2000～2020年水域动态增长度为8.51%，水域呈持续扩张趋势，其中2000～2010年的动态增长度为3.54%；2010～2020年的动态增长度为3.66%，在后期变化较大。林地基数大，2000～2020年动态增长度为0.01%，其中，2000～2010年增长度仅为0.1%，2010～2020年递减率为−0.09%。草地的动态变化不大，呈先增长后递减趋势，耕地呈持续递减趋势，未利用土地动态变化不明显。
表 2 丹江源地区2000～2020年土地动态度变化Table 2. Changes of land dynamic attitude from 2000 to 2020 in Danjiangyuan area
土地利用类型土地动态度 2000～2010年 2010～2020年 2000～2020年 耕地
林地
草地
水域
建设用地
未利用土地−0.32%
0.10%
−0.67%
3.54%
3.98%
−0.08%−1.05%
−0.09%
0.20%
3.66%
13.89%
0.07%−0.67%
0.01%
−0.24%
8.51%
11.70%
−0.05%3.2.2 土地利用空间变化
丹江源地区受地形条件，土地利用格局具有特殊性（图7）。林地占绝对优势，以天然林区为主，主要分布在杨斜镇、麻街镇、板桥镇、沙河子镇，草地主要分布在腰市镇、城关街道，耕地和建设用地分布规律明显，耕地主要分布在河谷川塬两岸及东北部的腰市镇、大荆镇，建设用地均分布于水系两岸，城区建设用地分布于丹江主河道两岸。不同土地利用类型时空转移特征明显，其中，耕地与林地之间的时空转移特征变化明显，且耕地向林地转移较明显，变化区域随时间变化先由中部向两侧扩散；耕地与建设用地之间的变化区域主要集中在商州城区、腰市镇和大荆镇，呈带状分布（图7），这与城市发展和经济发展有关。研究区未来的土地利用格局还以林地为主，但林地有退化趋势，且耕地后备资源不足，说明土地利用结构尚需进一步优化。
![]()
图 7 丹江源地区2000、2010、2020年土地利用类型遥感解译变化图Figure 7. Remote sensing interpretation changes of land use types of Danjiangyuan area in 2000, 2010 and 20202000～2020年土地利用时空演变结果总体显示丹江源的生态趋势向好有所改善。早期，当地居民伐林开垦，大规模修建坡耕地。2006年，商洛市积极调整土地利用结构及方式，以生态清洁小流域为生态建设单元，经过多年的综合治理和发展，当地生态环境得到了极大的改善。2011年底，圆满完成陕西省丹江口库区及上游水土保持一期工程的建设。2012年，商洛市启动陕西省丹江口库区及上游水土保持二期工程，2014年已顺利实现“清水进京”的宏伟目标，自通水至今，丹江的监测断面水质全部达到功能区标准。在生态建设保护的政策引导下，丹江源地区的土地利用变化对生态环境扰动较少。林地和耕地将长期是水源地主要土地利用类型，直接影响着其生态环境效应。
由于南水北调中线工程建设的需要，一方面，当地政府及民众对水源地和林地的保护意识不断增强，对公共资源的使用日趋合理，特别是对“山水林田湖草”生态文明建设的贯彻执行，进一步保障了研究区林地基数和水域面积的稳定增长，对可持续发展政策的实施有很大的正面效应。另一方面，由于调水工程及城市化进程建设以及保护水源地的移民搬迁，区内建设用地大幅度增加，存在草地退化和耕地减少现象。作为南水北调中线工程典型的山区水源地，今后应合理配置土地资源，优化土地利用结构，加强林地和耕地保护，将有利于丹江源地区社会经济和资源环境的可持续发展。受自然环境条件和社会经济水平限制，该区也是经济限制发展区。水源保护工作取得了一定成效，但还存在生态保护与经济发展矛盾突出、治污项目资金缺口大等困难和问题。
3.3 景观格局动态变化分析
3.3.1 斑块类型水平下景观格局指数
结合土地利用变化的解译，利用Frastats 4.2软件，计算丹江源地区2000～2020年3期景观格局指数。在类型水平上（表3），各土地利用类型的景观指数呈现出不同的变化趋势。林地景观具有很强的优势性，以其为代表的自然景观斑块呈现破碎状态。从2000～2020年的斑块类型比PLAND数据来看，未发生太大变化，呈微弱的减少趋势。人造景观耕地则呈现先增加后减少的变化趋势，草地变化趋势相反，先减后增，表明区内景观生态过程较活跃。
表 3 丹江源地区2010～2020年3期斑块类型水平的景观格局指数Table 3. Landscape pattern indices of patch types in the Danjiangyuan area from 2010 to 2020斑块类型 年份 CA（km2） PLAND（%） LPI（%） LSI IJI AI 林地 2000 1240.4300 76.0300 72.6400 17.4570 39.6831 98.0450 2010 1253.2700 75.8700 72.4600 17.6510 34.3240 98.0204 2020
20001242.2100
321.120075.1900
20.810071.4900 16.8271 37.6794 98.1100
88.0625耕地 20.2100 53.5204 53.9354 2010 310.8700 21.0200 20.4700 53.1127 47.0259 88.2095 2020 278.1700 20.8100 17.1600 57.3163 48.3957 87.2000 草地 2000 171.9900 2.7400 0.1900 84.8531 28.3333 47.1916 2010 160.3900 2.5300 0.0500 85.8279 23.7656 44.2992 2020 163.6400 2.7500 0.0500 95.8660 19.9750 40.3615 水域 2000 2.6800 0.0400 0.0300 8.1220 71.1135 62.7075 2010 3.6300 0.0600 0.0400 6.5652 62.6626 73.7705 2020 4.9600 0.0800 0.0700 6.7037 72.2046 77.8098 建设用地 2000 20.1800 0.3400 0.1000 12.7699 47.8730 78.5033 2010 28.2200 0.4700 0.1400 13.2256 42.9460 81.0820 2020 67.4100 1.1300 0.8000 14.8398 49.5618 86.4070 未利用土地 2000 0.4700 0.0400 0.0800 2.1938 31.5743 38.8427 2010 0.5000 0.0500 0.0900 2.7472 32.2674 39.3754 2020 0.4900 0.0400 0.0800 2.3584 31.6832 38.7120 注：由于各土地利用类型斑块类型面积统计存在四舍五入，面积总和与研究区总面积有略微出入。 从景观形状指数LSI变化趋势可知，耕地、草地的景观形状指数较大且总体呈现上升趋势，反映了区内耕地、草地景观和斑块形状分散且不规则，说明二者的景观形状比较复杂。林地和水域的景观形状指数有所下降，建设用地的景观形状指数呈上升趋势。以建设用地为代表的人造景观斑块面积呈现增大趋势， 2000～2010年增长不明显；2010～2020年增长迅速，斑块呈现片状分布。合理的土地利用规划让建设用地聚集程度提高，边界规则化，特别是移民搬迁政策的执行，使得建设用地聚集程度高。水域斑块类型面积呈现增加的趋势，斑块类型比和最大斑块指数总体呈上升趋势，不断增加并趋于平稳，水域占研究区面积的比例持续上升，散布与并列指数减小，景观形状指数和聚集度指数与水域相邻的景观要素变少，连通性增强。
3.3.2 景观水平下景观格局指数
景观水平的景观格局指数显示（表4），2000～2020年间丹江源地区的景观结构发生了变化，景观形状指数LSI呈现先有微弱减小后又增加趋势，总体呈增加，各类景观趋于离散，且形状趋于规则；蔓延度指数CONTAG在2000～2020年间先增加后减小，总体呈增大趋势，反映景观类型的聚集程度和延展程度，各景观之间连通性较好。散布与并列指数IJI呈减小趋势，反映了区内斑块类型的隔离分布情况，丰富度变低。香农多样性指数SHDI在2000～2020年间增加了0.031，说明区内的各斑块类型呈均衡趋势分布，土地利用类型越来越丰富，景观邻近度变好。聚集度指数AI在2000～2010年有上升，在2010～2020年有下降，先增加后减小。
表 4 丹江源地区2000~2020年景观水平的景观格局指数Table 4. Landscape pattern index of the Danjiangyuan area from 2000 to 2020年份 TA（km2） LSI CONTAG IJI SHDI AI 2000 1757.0000 27.8229 73.2957 45.5026 0.6557 94.4950 2010 1757.0000 27.7720 75.9058 39.2190 0.6596 94.5070 2020 1757.0000 29.7378 74.7040 40.4132 0.6963 94.1028 注：TA. 景观面积；LSl. 景观形状指数；CONTAG. 蔓延度指数；IJI. 散布与并列指数；SHDI. 香农多样性指数；Al. 聚集度指数。 2000～2010年，水域面积、林地面积及建设用地增加，耕地和草地面积减少，南水北调中线工程的建设，及从保护和改善生态环境出发退耕还林政策的实施，保障了丹江源地区水域面积和林地面积的增长；另一方面是移民搬迁政策，使原本分散的聚居地转为城镇、移民新村等，整体上建设用地面积增加，聚集度提升。2010～2020年，水域、草地及建设用地增加，林地、耕地面积减少，当地政府加大生态文明建设，对水源地及特殊土地类型湿地的保护政策不断推进，保障了水域面积的稳定增长，尤其是湿地保护成效显著，但随着近十年城市发展建设的需求，建设用地增加幅度较大，人口增长和社会需求也加速了土地利用程度的变化。
4. 结论
（1）各地质建造因构造地质背景差异导致上覆地形地貌、成土母质、土壤类型、水文类型、局部小气候类型，在以上立地条件的综合影响下形成了不同的生态地质格局。
（2）丹江源地区林地比例占绝对优势，林地和耕地占比85%以上，其他土地类型比例偏低。2000～2020年土地利用类型变化显著，水域、耕地和建设用地面积变化较为剧烈。区内自然景观和人造景观之间转入转出频繁，但整体土地利用变化对生态环境扰动较少。林地和耕地将长期是水源地主要的土地利用类型，直接影响着其生态环境效应。
（3）丹江源地区2000～2020年景观空间格局变化明显，尤其在2010～2020年景观生态过程较活跃。区内各斑块类型趋于规则呈均衡趋势分布，斑块类型间形成了良好的连接性，景观聚集程度逐渐提升，空间分布趋向集中。
（4）南水北调中线工程建设进一步保障了研究区林地基数和水域面积的稳定增长，丹江源地区的生态趋势向好；调水工程及城市化进程建设及保护水源地的移民搬迁导致区内建设用地大幅度增加，存在草地退化和耕地减少现象。建议合理配置土地资源，优化土地利用结构，加强林地和耕地保护及生态补偿，将有利于丹江源地区社会经济和生态环境的可持续发展。
致谢：在本文撰写过程中得到国际竹藤中心漆良华教授、博士生张建的帮助，商洛市自然资源局、商洛市林业局、商州区林业局等兄弟单位在资料收集方面给予了大力支持，在此一并致谢。
-
![]()
图 2 2020～2020年汉中盆地土地利用类型分布图
a. 2000年土地利用类型分布图；b. 2010年土地利用类型分布图；c. 2020年土地利用类型分布图
Figure 2. Distribution of land use types in Hanzhong basin in 2020-2022
![]()
图 3 土地转移矩阵弦图
a. 2000～2010年；b. 2010～2020年；c. 2000～2020年
Figure 3. Chord diagram of land stochastic matrix
![]()
图 4 生境质量变化转移矩阵弦图
a. 2000～2010年；b. 2010～2020年；c. 2000～2020年
Figure 4. Chord diagram of land stochastic matrix
![]()
图 5 生境质量时空变化图
a. 2000年；b. 2010年；c. 2020年
Figure 5. Spatial and temporal variation map of habitat quality
![]()
图 6 土地利用扩张时空变化图
a. 2000～2010年；b. 2010～2020年
Figure 6. Temporal and spatial changes in land use expansion
![]()
图 7 建设类用地扩张驱动因素排名（a）和建设类用地扩张因子贡献分析图（b）
Figure 7. (a) Ranking of driving factors for expansion of construction land and（b) Analysis of the Contribution of Construction Land Expansion Factors
![]()
图 8 2030年土地利用预测分析图
a.勉县区域；b.汉中区域；c.固城区域；d.洋县区域
Figure 8. Land use forecast analysis chart in 2030
表 1 驱动因子信息表
Table 1 Driver information table
类型 数据 数据来源 社会经济数据 人口 https://www.resdc.cn/Default.aspx GDP 到一级道路距离 Openstreetmap（https://github.com/openstreetmap） 到二级道路距离 到三级道路距离 到铁路距离 到政府距离 气候与环境数据 土壤类型 HWSDv1.2（https://www.fao.org/soils-portal/en/） 年平均温度 WorldClim v2.0（https://www.fao.org/soils-portal/en/） 年平均降水量 高程 NASA SRTM1 v3.0 坡度 到水域距离 （https://github.com/openstreetmap） 
下载: 导出CSV
表 2 威胁因子信息表
Table 2 Threat factor information table
威胁因子 最大影响距离（km） 权重 衰减类型 耕地 1 0.2 线性衰退 水域 1 0.2 线性衰退 建设类用地 5 0.5 指数衰退 农村居民点 6 0.7 指数衰退 城镇用地 10 1 指数衰退
下载: 导出CSV
表 3 威胁因子敏感度
Table 3 Threat factor sensitivity
土地利用类型 生境适宜度 威胁因子 耕地 水域 建设类用地 农村居民点 城镇用地 耕地 0.3 0.6 0.5 0.6 0.7 0.5 林地 0.9 0.5 0.8 0.7 0.3 0.8 草地 0.8 0.8 0.8 0.4 0.7 0.4 水域 0.7 0.7 0.2 0.1 0.3 0.3 城乡工矿居民 0 0 0 0.2 0 0 未利用土地 0.3 0.2 0.2 0.1 0.2 0.3
下载: 导出CSV
表 4 汉中盆地景观类型水平
Table 4 Landscape type level of Hanzhong basin
类型 2000年 2010年 2020年 PLAND LPI TE ED PLAND LPI TE ED PLAND LPI TE ED 水田 71.28 39.4 5320878.1 31.79 68.07 32.80 5576238.6 33.31 66.03 24.3 5963113.8 35.63 旱地 5.05 1.36 610927.42 3.65 5.18 1.39 646091.46 3.86 5.16 1.35 675959.51 4.04 有林地 0.14 0.02 48228.74 0.29 0.13 0.02 45581.76 0.27 0.13 0.02 45701.71 0.27 灌木林 0.05 0.03 12057.19 0.07 0.05 0.02 13194.72 0.08 0.04 0.02 11995.20 0.07 疏林地 0.76 0.03 271406.66 1.62 0.76 0.03 276039.54 1.65 0.76 0.03 275679.68 1.65 其他林地 0.89 0.04 289552.42 1.73 1.02 0.06 312325.02 1.87 1.01 0.06 317812.82 1.90 高覆盖度草地 0.06 0.03 19345.49 0.12 0.07 0.03 22790.88 0.14 0.08 0.03 26629.34 0.16 中覆盖度草地 9.12 0.17 3140447.0 18.76 8.53 0.17 3031306.9 18.1 8.89 0.17 3139263.7 18.76 低覆盖度草地 0.00 0.00 1439.66 0.01 0.00 0.00 1379.45 0.01 0.00 0.00 1319.47 0.01 河渠 1.22 0.99 333072.27 1.99 0.86 0.14 243802.44 1.46 1.27 1.01 294332.22 1.76 湖泊 0.69 0.08 254730.55 1.52 0.65 0.03 244462.18 1.46 0.65 0.03 243622.51 1.46 水库坑塘 0.27 0.05 68773.95 0.41 0.59 0.25 156447.40 0.93 0.61 0.26 159386.22 0.95 滩地 3.04 0.61 535225.09 3.20 3.18 0.71 549859.97 3.29 2.78 0.64 509915.95 3.05 城镇用地 0.78 0.30 55217.11 0.33 1.90 0.75 126849.24 0.76 2.29 0.76 167183.10 1.00 农村居民点 6.51 0.06 1590768.7 9.50 8.69 0.28 1959835.7 11.71 8.73 0.20 1945261.5 11.62 其他建设用地 0.10 0.06 14756.56 0.09 0.28 0.09 51969.20 0.31 1.51 0.10 520801.60 3.11 沼泽 0.03 0.02 5398.74 0.03 0.02 0.02 5457.82 0.03 0.02 0.02 5457.82 0.03 裸土地 − − − − 0.01 0.00 2938.82 0.02 0.04 0.01 18112.75 0.11
下载: 导出CSV
表 5 汉中盆地景观水平
Table 5 Landscape type level of Hanzhong basin
年份（a） CA LPI TE ED CONTAG SHDI SHEI 2000 167379.642 39.486 6286112.898 37.556 75.089 1.143 0.403 2010 167379.866 32.804 6683185.656 39.928 73.938 1.247 0.432 2020 167379.866 24.304 7210674.576 43.080 72.106 1.330 0.460
下载: 导出CSV
-
陈颐, 林毅伟, 林丽丽, 等. 基于Markov和LoGIStic模型的莆田市土地利用变化及林地转出空间模拟[J]. 中国农业大学学报, 2017, 22（2）: 87−97. doi: 10.11841/j.issn.1007-4333.2017.02.011 CHEN Yi, LIN Yiwei, LIN Lili, et al. Spatial simulation of land use change and forest land transfer in Putian City based on Markov and LoGIS models[J]. Journal of China Agricultural University,2017,22(2):87−97. doi: 10.11841/j.issn.1007-4333.2017.02.011
曹雪, 罗平, 李满春, 等. 基于扩展CA模型的土地利用变化时空模 拟研究——以深圳市为例[J]. 资源科学, 2011, 33（1）: 127−133. CAO Xue, LUO Ping, LI Manchun, et al. Study on spatiotemporal simulation of land use change based on extended CA model: A case study of Shenzhen City[J]. Resource Science,2011,33(1):127−133.
陈竹安, 刘子强, 张立亭, 等. 南昌市LUCC多情景模拟 和生境质量时空演变与预测[J]. 农业机械学报, 2023, 54（5）: 170−180. CHEN Zhu’an, LIU Ziqiang, ZHANG Liting, et al. Multi scenario simulation and spatiotemporal evolution and prediction of habitat quality in Nanchang LUCC[J]. Journal of Agricultural Machinery,2023,54(5):170−180.
代克志. 基于3S技术的花溪河流域景观格局变化及影响研究[D]. 贵州: 贵州师范大学, 2015. DAI Kezhi. Research on Landscape Pattern Changes and Impact in the Huaxi River Basin Based on 3S Technology [D]. Guizhou: Guizhou Normal University, 2015.
冯雨林, 杨佳佳, 吴梦红. 基于景观转移矩阵的黑龙江双河自然护区土地覆被转移研究[J]. 地质与资源, 2016, 25（5）: 1−5. FENG Yulin, YANG Jiajia, WU Menghong. Research on land cover transfer in the Shuanghe Nature Reserve of Heilongjiang Province based on landscape transfer matrix[J]. Geology and Resources,2016,25(5):1−5.
范树平, 程从坤, 刘友兆, 等. 中国土地利用/土地覆盖研究综述与展望[J]. 地域研究与开发, 2017, 36（2）: 94−101. doi: 10.3969/j.issn.1003-2363.2017.02.018 FAN Shuping, CHENG Congkun, LIU Youzhao, et al. A Review and Prospects of Land Use/Land Cover Research in China[J]. Regional Research and Development,2017,36(2):94−101. doi: 10.3969/j.issn.1003-2363.2017.02.018
高文龙, 张景华, 刘洪, 等. 基于转移矩阵的大凉山区土地利用变化研究[J/OL]. 西北地质, 2023: 1−13. doi: 10.12401/j.nwg.2023167. GAO Wenlong,ZHANG Jinghua,LIU Hong,et al. Research on Land Use Change in Daliangshan District Based on Transfer Matrix [J/OL]. Northwest Geology,2023: 1−13. Doi: 10.12401/j.nwg.2023167.
黄钰清, 李骁尧. 1995-2020年黄河流域土地利用变化及驱动力分析[J]. 西北林学院学报, 2020, 37（6）: 113−121. HUANG Yuqing, LI Xiaoyao. Analysis of Land Use Changes and Driving Forces in the Yellow River Basin from 1995 to 2020[J]. Journal of Northwest Forestry University,2020,37(6):113−121.
韩朝辉, 王郅睿, 田辉, 等. 汉中盆地地下水水化学特征及其成因研究[J]. 西北地质, 2023, 56（4）: 263−273. HAN Chaohui, WANG Zhirui, TIAN Hui, et al. Hydrochemical Characteristics and Genesis of Groundwater in the Hanzhong Basin[J]. Northwestern Geology,2023,56(4):263−273.
黄金亭, 曹艳萍, 秦奋. 基于土地利用/覆盖变化模拟的黄河流域生态环境质量分析[J]. 河南大学学报: 自然科学版, 2020, 50（2）: 127−138. HUANG Jinting, CAO Yanping, QIN Fen. Analysis of Ecological Environment Quality in the Yellow River Basin Based on Land Use/Cover Change Simulation[J]. Journal of Henan University: Natural Science Edition,2020,50(2):127−138.
胡丰, 张艳, 郭宇, 等. 基于PLUS和InVEST模型的渭河流域土地利用与生境质量时空变化及预测[J]. 干旱区地理, 2011, 45（4）: 1125−1136. HU Feng, ZHAGN Yan, GUO Yu, et al. Spatial and temporal changes and prediction of land use and habitat quality in the Wei River Basin based on PLUS and InVEST models[J]. Arid Region Geography,2011,45(4):1125−1136.
吕桂军. 长时间序列LandsatTM影像应用于地理国情监测中的 技术方法探讨[D]. 焦作: 河南理工大学, 2014. LV Guijun. Exploration of Technical Methods for Applying Long Time Series LandsatTM Images to Geographical and National Monitoring [D]. Jiaozuo: Henan University of Technology, 2014.
梁甜, 黄茜, 杨霏. 基于InVEST-PLUS模型的三峡库区(重庆段)生境质量演变及预测[J]. 长江流域资源与环境, 2023, 32（10）: 2184−2195. LIANG Tian, HUANG Qian, YANG Fei. Habitat quality evolution and prediction in the Three Gorges Reservoir Area (Chongqing section) based on the InVEST-PLUS model[J]. Yangtze River Basin Resources and Environment,2023,32(10):2184−2195.
刘纪远, 宁佳, 匡文慧, 等. 2010-2015年中国土地利用变化的时空格局与新特征[J]. 地理学报, 2020, 73（5）: 789−802. LIU Jiyuan, NING Jia, KUANG Wenhui, et al. The spatiotemporal pattern and new characteristics of land use change in China from 2010 to 2015[J]. Journal of Geography,2020,73(5):789−802.
李灿, 曾和平. 基于面积转移矩阵的龙川江流域土地利用变化研究[J]. 人民长江, 2018, 49（17）: 06−13. LI Can, ZENG Heping. Research on Land Use Change in the Longchuan River Basin Based on Area Transfer Matrix[J]. Yangtze River,2018,49(17):06−13.
陶蕴之, 张甜, 梁琦珍, 等. 基于转移矩阵的环渤海土地利用研究[J]. 绿色科技, 2016, 37（2）: 4−7. doi: 10.3969/j.issn.1674-9944.2016.02.002 TAO Yunzhi, ZHANG Tian, LIANG Qizhen, et al. Research on Land Use of Bohai Rim Based on Transfer Matrix[J]. Green Science and Technology,2016,37(2):4−7. doi: 10.3969/j.issn.1674-9944.2016.02.002
卫建军, 李新平, 杨勤科. 基于遥感的土地利用与退耕还林时空变化研究—以延河流域下游地区为例[J]. 水土保持通报, 2007（2）: 87−90. doi: 10.3969/j.issn.1000-288X.2007.02.020 WEI Jianjun, LI Xinping, YANG Qinke. Research on spatiotemporal changes of land use and returning farmland to forests based on remote sensing: A case study of the lower reaches of the Yanhe River Basin[J]. Soil and Water Conservation Bulletin,2007(2):87−90. doi: 10.3969/j.issn.1000-288X.2007.02.020
王颖维, 张亚峰, 钱信禹, 等. 南水北调中线工程丹江源地区生态地质格局动态演变[J]. 西北地质, 2023, 56（3）: 129−140. WANG Yingwei, ZHANG Yafeng, QIAN Xinyu, et al. Dynamic Evolution of Eco−geological Pattern: Taking Danjiangyuan Area of Central Line Project of South−to−North Water Diversion[J]. Northwestern Geology,2023,56(3):129−140.
杨朔, 苏昊, 赵国平. 基于PLUS模型的城市生态系统服务价值多情景模拟: 以汉中市为例[J]. 干旱区资源与环境, 2022, 36（10）: 86−95. YANG Shuo, SU Hao, ZHAO Guoping. Multi scenario simulation of urban ecosystem service value based on PLUS model: A case study of Hanzhong City[J]. Resources and Environment in Arid Areas,2022,36(10):86−95.
阳柏苏. 景区土地利用格局及生态系统服务功能研究[D]. 长沙: 中南林业科技大学, 2005. YANG Baisu. Research on Land Use Patterns and Ecosystem Services in Scenic Spots [D]. Changsha: Central South University of Forestry and Technology, 2005.
岳东霞, 杜军, 刘俊艳, 等. 基于RS和转移矩阵的泾河流域生态承载力时空动态评[J]. 生态学报, 2011, 31（9）: 2550−2558. YUE Dongxia, DU Jun, LIU Junyan, et al. Spatial and temporal dynamic evaluation of ecological carrying capacity in the Jinghe River Basin based on RS and transfer matrix[J]. Journal of Ecology,2011,31(9):2550−2558.
张冉, 王义民, 畅建霞, 等. 基于水资源分区的黄河流域土地利用变化对人类活动的响应[J]. 自然资源学报, 2019, 34（2）: 274−287. doi: 10.31497/zrzyxb.20190205 ZHANG Ran, WANG Yimin, CHANG Jianxia, et al. Response of land use change in the Yellow River Basin to human activities based on water resource zoning[J]. Journal of Natural Resources,2019,34(2):274−287. doi: 10.31497/zrzyxb.20190205
张立新, 段文凯. 经济转型背景下城市土地城镇化格局 及影响因素—以长江经济带城市为例[J]. 国农业大学学报, 2021, 26（9）: 206−215. ZHANG Lixin, DUAN Wenkai. Urban Land Urbanization Pattern and Influencing Factors under the Background of Economic Transformation: A Case Study of Cities in the Yangtze River Economic Belt[J]. Journal of China Agricultural University,2021,26(9):206−215.
张恒义, 刘卫东, 王世忠, 等. “省公顷”生态足迹模型中均衡因子及产量因子的计算: 以浙江省为例[J]. 自然资源学报, 2009, 24（1）: 82−92. doi: 10.3321/j.issn:1000-3037.2009.01.010 ZHANG Hengyi, LIU Weidong, WANG Shizhong, et al. Calculation of Equilibrium Factors and Yield Factors in the Ecological Footprint Model of "Provincial Hectares": Taking Zhejiang Province as an Example[J]. Journal of Natural Resources,2009,24(1):82−92. doi: 10.3321/j.issn:1000-3037.2009.01.010
Faichia C, Tong Z, Zhang J, et al. Using RS Data-Based CA–Markov Model for Dynamic Simulation of Historical and Future LUCC in Vientiane, Laos[J]. Sustainability,2020,12(20):8410−8419. doi: 10.3390/su12208410
Lambin EF, Turbner B, Geist HJ, et al. The cauyses of land-use and land-over change: moviong beyong the myths[J]. Global enviroment change,2001,11(4):261−269. doi: 10.1016/S0959-3780(01)00007-3
Newbold T, Hudson L N, Phillips HR, et al. A global model of the response of tropical and sub-tropical forest biodiversity to anthropo genic pressures[J]. Proceedings Biological Sciences,2014,281(1792):1−10.
Verburg P H, De konnig G, Kok K, et al. A spatial explicit allocation procedure for modeling the pattern of land use change based upon actual land use[J]. Ecoligical modelling,1999,116(1):45−61. doi: 10.1016/S0304-3800(98)00156-2
-
期刊类型引用(4)
1. 付柯锦. 生态环境保护背景下南水北调中线工程水源区生态旅游发展的路径选择. 地域研究与开发. 2025(01): 137-142+150 . 
百度学术
2. 袁江龙，刘晓煌，李洪宇，邢莉圆，雒新萍，王然，王超，赵宏慧. 1990—2050年黄河中游伊洛河流域不同土地利用类型碳储量时空分异特征. 现代地质. 2024(03): 559-573 . 百度学术
3. 刘小玉，李士杰，何海洋，秦昊洋，王思琪，孙旭. 基于InVEST和PLUS模型下的土地利用变化及生境质量演变分析：以汉中盆地为例. 西北地质. 2024(04): 271-284 . 本站查看
4. 欧阳渊，刘洪，张景华，唐发伟，张腾蛟，黄勇，黄瀚霄，李富，陈敏华，宋雯洁. 西南山区生态地质调查技术方法研究. 西北地质. 2023(04): 218-242 . 本站查看
其他类型引用(2)
计量
- 文章访问数: 119
- HTML全文浏览量: 7
- PDF下载量: 36
- 被引次数: 6
