环境科学  2026, Vol. 47 Issue (9): 6264-6278   PDF    
基于FLUS-InVEST模型的黑河流域三生空间生境质量演变及多情景模拟预测
杨鸿魁1,2,3, 李业1, 马玉欣1, 祁万强1,2, 杨明新1,2, 欧立鹏1     
1. 中国地质调查局西宁自然资源综合调查中心,西宁 810021;
2. 祁连山中部亚高山生态系统野外科学观测研究站,张掖 734000;
3. 青海师范大学地理科学学院,西宁 810008
摘要: 研究三生空间生境质量演变时空格局,对于掌握区域生态现状和可持续发展具有实际意义. 以黑河流域1990~2023年8期土地利用数据为基础对三生空间类型进行划分. 基于FLUS-InVEST模型和GIS平台,通过耦合11个驱动要素来模拟预测黑河流域2035年生产空间优先、生活空间优先和生态空间优先发展等3种情景下的三生空间及生境质量时空分布格局. 结果表明:①1990~2023年黑河流域农业生产空间面积净增1 184.23 km2,主要来自未利用地和草地的转入;水域生态用地面积增幅达20.96%;林、草地生态用地面积分别减少167.04 km2和606.31 km2;②1990~2023年黑河流域生境质量年际变化率为-0.013%·a-1. 较高质量和高质量区主要向中等质量和较低质量区发生转移,面积减少786.15 km2,其重心分别向北西和北东迁移17.61 km和21.66 km;低质量区面积增加7 178.28 km2,低质量区重心向中上游绿洲区迁移52.77 km;③较2020年,2035年生产空间优先情景中林地和草地面积分别减少了271.23 km2和800.71 km2;生活空间优先发展下,城镇生活空间显著增加,林地和草地生态用地面积持续减少;生态空间优先发展下林地面积增加17.76 km2. 该情景下草地面积减少幅度最低;④2035年生产、生活和生态空间优先发展情景的HQ均值分别为0.376 9、0.377 1和0.381 5. 其中,生产和生活两种情景的HQ较2020年降幅分别为1.13%和1.08%. 主因在于较高质量和高质量区面积的减少. 而生态优先情景下,高质量区面积增加234.70 km2,HQ较2020年和其他两种情景增幅分别为0.08%、1.22%和1.17%. 研究结果可为黑河流域生态治理和未来三生空间的协调发展提供多视角参考.
关键词: 三生空间      生境质量(HQ)      FLUS-InVEST模型      时空格局      黑河流域     
Evolution and Multi-scenario Simulation Prediction of Habitat Quality of Production Living Ecological Spaces in the Heihe River Basin Based on the FLUS-InVEST Model
YANG Hong-kui1,2,3 , LI Ye1 , MA Yu-xin1 , QI Wan-qiang1,2 , YANG Ming-xin1,2 , OU Li-peng1     
1. Xining Natural Resources Comprehensive Survey Center, China Geological, Xining 810021, China;
2. Observation Station of Subalpine Ecology Systems in the Middle Qilian Mountains, Zhangye 734000, China;
3. School of Geography Science, Qinghai Normal University, Xining 810008, China
Abstract: Studying the spatio-temporal pattern of the evolution of the quality of the Production Living Ecological Spaces in the Heihe River Basin is of practical significance for understanding the current ecological situation and promoting sustainable development in the region. This study classifies the Production Living Ecological Spaces based on the land use data of the Heihe River Basin from 1990 to 2023. Using the FLUS-InVEST model and GIS platform and coupling 11 driving factors, it simulates and predicts the spatio-temporal distribution patterns of the Production Living Ecological Spaces and their habitat quality under the three scenarios of production space priority, living space priority, and ecological space priority in 2035. The results show that: ① From 1990 to 2023, the area of agricultural production space in the Heihe River Basin increased by 1 184.23 km2, mainly due to the conversion of unused land and grassland; the area of water ecological land increased by 20.96%; and the areas of forest and grassland ecological land decreased by 167.04 km2 and 606.31 km2, respectively. ② From 1990 to 2023, the annual change rate of habitat quality in the Heihe River Basin was -0.013%·a-1. High-quality and very high-quality areas mainly transferred to medium-quality and low-quality areas, with a reduction of 786.15 km2 in area. Their centers of gravity moved 17.61 km to the northwest and 21.66 km to the northeast, respectively; the area of low-quality areas increased by 7 178.28 km2, and the center of gravity of low-quality areas moved 52.77 km towards the oasis area in the middle and upper reaches. ③ Compared with those in 2020, in the production space priority scenario in 2035, the areas of forest and grassland decreased by 271.23 km2 and 800.71 km2, respectively, under the living space priority development scenario, and the area of urban living space increased significantly, while the areas of forest and grassland ecological land continued to decrease. Under the ecological space priority development scenario, the area of forest land increased by 17.76 km2. The decrease in grassland area was the smallest in this scenario. ④ In 2035, the average HQ values of the Production Living Ecological Spaces priority development scenarios were 0.376 9, 0.377 1, and 0.381 5, respectively. Among them, the HQ values of the Production Living Ecological Spaces priority scenarios decreased by 1.13% and 1.08% compared with those in 2020, mainly due to the reduction in the area of high-quality and very high-quality areas. However, under the ecological space priority scenario, the area of high-quality areas increased by 234.70 km2, and the HQ value increased by 0.08%, 1.22%, and 1.17% compared with that in 2020 and the other two scenarios, respectively. The research results can provide multi-perspective references for ecological governance and the coordinated development of the Production Living Ecological Spaces in the Heihe River Basin in the future.
Key words: Production Living Ecological Spaces      habitat quality(HQ)      FLUS-InVEST model      spatial and temporal pattern      Heihe River Basin     

生境质量是指在一定空间范围内,生态系统可为生物群体提供生存繁衍所需的能量水平[1],是地理科学与环境科学等相关领域研究的核心内容,与区域土地利用结构的空间分布密切相关[2]. 三生空间(即生产空间、生活空间和生态空间)是国土空间规划中具有空间功能属性定位的核心类型[3],是承载生态文明建设的功能单元,是实现人与自然和谐共生的“顶层设计”[4],而土地利用类型则是三生空间的具体落地形式[5]. 不科学、不规划的土地利用必然会造成用地结构破坏,加剧水土流失和土地碎斑化等功能性衰退局面的产生[6],进而影响区域整体生态状况.

国内外学者对于生境质量的研究主要分为野外调查[7]和模型评估[8]两类. 野外调查是比较传统的工作手段,数据采集质量好、精度高,但在大尺度研究上需要消耗更多的人力物力财力. 模型评估方面较常见的有InVEST模型[9]、生境适宜性指数HSI模型[10]、最大熵MAXENT模型[11]及生态位模型[12]等. 其中,InVEST模型兼具综合性、多尺度适用性,便于可视化且数据需求简洁,运行机制流畅等优势,在区域尺度上更适合整体生境质量评价,特别是在地类转移频繁和人类干扰剧烈的区域[13,14]. 土地利用模拟预测模型的发展从初期的统计模型发展为空间显性模型和多模型耦合,如Markov模型、元胞自动机CA模型、CLUE-S模型、FLUS模型以及PLUS模型等,各具优势. FLUS模型则是通过整合人工神经网络ANN与自适应惯性机制,在处理非线性驱动机制和降低模拟噪声方面占优,可有效提升模拟精度,在行政区划[15]、自然保护区[16]、重要经济区[17]以及流域[18]等不同尺度的相关研究中得到应用. 以往工作更多集中于历史时期上游、中游或者下游某个子流域尺度土地利用和生态系统服务方面的研究[19,20],且研究方向多集中在土地利用类型这个实际载体,研究时段多关注历史时期,而针对黑河流域三生空间生境质量模拟预测研究相对较少,以动态视角来研究生产、生活和生态空间的功能属性、探索未来情景下的研究区域生态状况相对匮乏. 因此,在长时间序列方面进行生境质量的模拟预测具可操作性.

目前,针对黑河流域的相关研究成果丰硕,主要包括生境质量评估[21]、碳储量模拟预测[22]、生态格局构建[23]、生态风险评估[24]、气候变化背景下水文水资源特征与“三生空间”转换[25,26]以及国土空间生态修复识别[27]等多个方面,研究成果强有力地支撑了黑河流域的可持续发展,但未来情景下的生境质量演变预测还需进一步研究,将对气候变化与人类活动双重影响下的黑河流域生态响应具有指导意义. 本文利用FLUS-InVEST模型,通过耦合自然地理和社会经济等11个驱动要素,基于“三生空间”视角来探究1990~2023年历史背景下的生境质量状况,并设置生产优先、生活优先以及生态优先发展等3种未来情景模式,分析2035年黑河流域“三生空间”与生境质量时空演变规律,旨在为黑河流域的生态治理、国土规划和可持续发展提供基础数据支撑.

1 研究区概况与数据来源 1.1 研究区概况

黑河流域为我国第二大内陆河流域,行政区划横跨青海、甘肃及内蒙古三省区,地理位置介于98°~101°30′E、38°~42°N之间,流域面积约14.29万km2[22](图 1). 因深居内陆,中下游一带气候干燥,降水稀少,蒸发强烈,昼夜温差大. 自上游而下,呈现出不同的生态系统景观类型,上游是重要的水源涵养区,山地垂直地带特征明显,分布有寒漠、草甸及森林等垂直带植被型;中游人口密集,经济相对发达,是重要的农业灌溉区,农业用水需求高;下游则是以戈壁、荒漠和沙漠等未利用地为主,额济纳旗绿洲区是重要的防风固沙生态屏障区,生长有胡杨、柽柳、骆驼刺和芦苇等耐旱植被. 近年来,受中游农业用水量增加的制约,下游居延海湿地的生态需水面临极大的风险挑战.

图 1 研究区概况 Fig. 1 Overview of the study area

1.2 数据来源与预处理

本文的土地利用数据来源于资源环境科学数据平台,CNLUCC数据分类主要是基于Landsat遥感影像,通过人工目视解译得到,其总体分类精度在89%以上,具有较高的分类精度. 该数据集采用二级分类系统,涵盖6个一级地类和25个二级地类,空间分辨率为30 m. 结合研究区实际并参考相关研究论文[25,28],对区域“三生空间”的地类划分进行界定(表 1);DEM数据来自地理空间数据云平台的ASTER GDEM,空间分辨率为30 m,坡向和坡度通过GIS工具计算获取;年均气温、年平均降雨量、植被覆盖度来源于地球资源数据云平台(http://www.gis5g.com/),空间分辨率为1 000 m;河流水系矢量数据来自黑河计划数据管理中心发布的黑河流域河流分布数据集[29],时间为2009年;道路矢量数据来自OpenStreetMapp(https://www.openstr eetmap.org);政府驻地矢量数据来自全国地理信息资源目录服务系统的1∶100万全国基础地理数据库;人口密度数据来自Landscan数据集(https://landscan.ornl.gov/),空间分辨率为1 000 m.

表 1 黑河流域“三生空间”土地利用分类体系 Table 1 Classification system of land use for Production Living Ecological Spaces in the Heihe River Basin

因模型运行需求,所有数据均需栅格化处理,距河流、距道路、距铁路以及距政府驻地的欧氏距离均由ArcGIS Pro3.0的欧氏距离工具处理得到. 行列号为14 250×18 361,空间分辨率30 m×30 m,投影坐标系统一设置为WGS_1984_UTM_46N. 驱动因子见表 2.

表 2 主要数据类型及来源说明 Table 2 Explanation of main data types and sources

2 材料与方法 2.1 FLUS模型

FLUS模型常用于自然条件与人为影响下的土地利用类型变化及未来情景下的地类模拟预测,在多种尺度下均有较好的模拟预测能力[30]. 该模型基于一期土地利用数据和多源驱动因子,采用神经网络算法提取到每种类型的适应性发展概率. 而后基于轮盘赌选择的自适应惯性竞争机制,通过调整邻域因子、惯性系数以及转换成本等,得到每个地类总转移概率[30],并结合Markov模型预测获取研究时段的土地利用类型需求数量. 通过设置类型预测需求总数、转换原则和领域权重,将FLUS与Markov模型融合,进行黑河流域三生空间的模拟预测,保证其模拟精度满足要求. 结合其他学者研究和试验对比分析,本文中FLUS模型采取随机采样法,采样参数为20,隐藏层数默认12层. 模拟参数保持默认:最大迭代次数300,领域为3×3,线程数结合计算机性能设置为8,提高运行速度.

TPp,kt=Pp,k×Ωp,kt×Ikt×SCc→k (1)

式中,TPp,kt表示单元p经过t时刻后转变为k类型的总概率;Pp,k表示单元p转变为k类型的发展概率;Ωp,kt表示邻域影响因子;Ikt表示自适应惯性系数;SCc→k表示c类型转换为k类型的转换成本.

2.1.1 领域权重

领域权重是一种表征三生空间类型转移或扩张能力强弱的指标[31],取值为0~1,值越接近1,表示该类型发生转移的概率越高,将优先分配进行类型间的转移. 参考其他学者在领域权重设置的方法[13,17],本文结合2010~2020年三生空间转移矩阵,通过两期数据面积变化差值确定领域权重值. 在符合模拟精度前提下,未来情景的领域权重依据初期(2020年)与末期(2035年)三生空间面积变化来设置(表 3). 领域权重的计算公式如下[22]:

Pi=Si-SminSmax-Smin (2)
表 3 2010~2020年黑河流域三生空间领域权重 Table 3 Weight of the Production Living Ecological Spaces domain in the Heihe River Basin from 2010 to 2020

式中,Pi为i类三生空间的领域权重;Si为i类三生空间的面积;Smax、Smin分别为i类三生空间面积的最大值与最小值.

2.1.2 模型精度验证

本文通过Kappa系数来检验模拟预测的准确性. Kappa系数基本原理是通过混淆矩阵来判断预测结果中类型正确的数量与实际结果中正确分类数量的比值,比值越高,分类精度越高,其模拟结果越准确. 当Kappa系数高于0.70时,说明模型精度较好[30]. 首先基于2010年三生空间数据和11个驱动因子模拟预测2020年三生空间类型分布,而后将模拟结果与2020年实际三生空间进行比较,判读模型预测精度(图 2). 经对比,Kappa系数为0.83,总体精度达0.92,符合精度要求,可进行研究区2035年多情景三生空间模拟预测.

图 2 黑河流域2020年三生空间预测分布与实际分布空间对比 Fig. 2 Spatial comparison between the predicted distribution and the actual distribution of the Production Living Ecological Spaces in the Heihe River Basin in 2020

2.1.3 三生空间多情景模拟设置

三生空间的科学合理规划旨在优化生产、生活及生态空间之间的协调关系,达到可持续发展的目的. 黑河流域上游、中游和下游具有着不同的生态角色,分别为水源涵养区、农灌区和防风固沙区. 2001年《黑河流域近期治理规划》中提出通过实施统一调度用水,减轻遏制黑河下游生态区退化与恶化等问题[32];阿拉善盟《黑河流域水生态环境保护“十四五”规划》明确落实黑河生态补水目标,通过水利工程完善和信息化技术提升,达到下游额济纳绿洲区面积的持续增加,并优化中游农业用水的灌溉效率,保障下游生态需水[33]. 基于此,本文共设置3种情景:生产空间优先情景、生活空间优先情景和生态空间优先情景,并结合实际将各情景转换原则进行设定(表 4).

表 4 2035年黑河流域三生空间多情景转换原则1) Table 4 Multi-scenario transformation principles for the Production Living Ecological Spaces in the Heihe River Basin in 2035

生产空间优先发展情景. 该情景在遵循2010~2020年三生空间地类转移规律的基础上,通过其他用地空间的转移来保护农业生产与工矿企业生产空间用地. 由于黑河流域的绝大部分耕地分布在中游一带,为确保足够的粮食供给安全,提取耕地作为该情景的限制区. 将农业生产空间和水域生态空间向林地和草地生态空间的转移概率降低20%;林地和草地生态空间向农业生产空间的概率增加10%,并将其他生态空间向生产空间的转移概率增加20%,以此来实现生产空间优先发展的目的.

生活空间优先发展情景. 仍以2010~2020年三生空间地类转移规律为基础,确保未来时期内拥有足够的生活空间. 该情景下不设置限制区,通过适当调整其他空间用地向生活空间用地的转移概率,实现优先发展. 将城镇和农村生活空间向农业生产空间、草地生态空间以及其他生态空间的转移概率降低20%,并将草地生态和其他生态空间向城镇生活与农业生产空间的转移概率增加20%.

生态空间优先发展情景. 生态优先是区域可持续发展的重要抓手,是提升区域生境质量的可行性举措[34]. 黑河流域生态治理修复是今后一定时期内常抓不懈的系统工程. 在该情景下,优先将自然保护区设置为限制区,降低生态地类向其他空间用地的转移,并适度调整其他空间用地向林草等生态用地的转移概率. 将农业生产空间向城镇和农村生活空间的转移降低10%,林地、草地和水域等生态空间向农业生产空间和其他生态空间转移降低40%,工矿企业空间向草地生态空间转移增加20%,未利用地等其他生态空间向林地和草地生态空间的转移概率增加40%,通过提升低质量生产用地与生活空间用地使用效率来实现生态优先发展的目的.

2.2 基于InVEST模型的生境质量评估

本文利用InVEST模型的生境质量模块(habitat quality,HQ)对黑河流域生境质量进行系统评估,包括1990~2023年历史时期的生境质量状况和2035年多情景下的生境质量时空分布特征,通过生境质量指数和生境退化度指数来衡量. 二者的取值范围均为[0, 1],值越接近1,表明其生境质量或者生境退化程度均较高. 该模块主要基于土地利用类型数据进行赋值计算,参考Zhang等[35]和Yan等[36]在黑河流域的相关研究工作以及InVEST模型指导手册,结合研究区三生空间分布特征,将农业生产空间、工矿企业生产空间、城镇生活空间、农村生活空间以及未利用地等其他生态空间设置为威胁因子,并对最大威胁距离、权重、生境适宜度和生境敏感性等指标进行赋值(表 5和表 6).

表 5 三生空间威胁因子与对应参数 Table 5 Threatening factors of the Production Living Ecological Spaces and their corresponding parameters

表 6 三生空间对不同威胁因子的生境适宜度与敏感性 Table 6 Habitat suitability and sensitivity of the Production Living Ecological Spaces to different threatening factors

2.3 三生空间转移矩阵

转移矩阵可以用来量化不同时期不同类型之间相互转化的空间关系,能够清晰揭示不同类型的数量变化与转移方向,在土地利用变化[37]和生态系统服务[38]等方面应用广泛,是一种含有丰富统计信息的定量描述方法. 利用ArcGIS 10.7的叠加分析工具,计算得到1990~2023年黑河流域三生空间转移矩阵. 其表达式如下[39]:

Sij=S11⋯S1n⋮⋮Sn1⋯Snn (3)

式中,S 表示三生空间面积;i和j分别表示研究时段内初期和末期的三生空间类型;n表示前后时段三生空间类型数.

2.4 三生空间重心迁移

重心迁移常用来揭示人口、经济以及国土规划等在空间范围上的运行轨迹[40]. 通过其迁移轨迹和迁移距离来判断三生空间用地在空间上的变化方向,可以更好优化用地结构和国土规划治理. 本文利用ArcGIS 10.7平均中心与椭圆标准差工具来完成黑河流域生境质量等级的重心计算. 其重心坐标计算公式如下[40]:

Xt=∑i=1n(Cti×Xti)∑i=1nCti; Yt=∑i=1n(Cti×Yti)∑i=1nCti (4)

式中,Xt和Yt为t时段某一类生境质量等级的经纬度坐标;n为生境质量分级数;Cti为t时段某生境质量等级第i个斑块的面积;Xti和Yti为t时段某生境质量等级第i个斑块的几何中心坐标.

3 结果与分析 3.1 黑河流域1990~2023年三生空间分布特征

黑河流域三生空间分布和面积在1990~2023年间发生了较显著的变化(表 7和图 3). 生产、生活和生态空间均呈波动变化,其中农业生产空间呈激增略减趋势,整体净增1 184.23 km2,增长率为23.76%;工矿企业用地、城镇生活空间和农村生活空间用地显著增加,增长率分别为464.69%、226.55%和12.89%,用地面积的增加与区域经济发展存在显著正相关性,主要集中在中游一带的人口集中区;生态用地中除水域生态用地面积增长外(20.96%),林地、草地和其他生态用地面积整体呈减少趋势,分别减少167.04、606.31和1 186.14 km2. 2010年和2015年林地、草地与水域生态用地面积显著减少,而未利用地等其他生态用地面积陡增(表 7),黑河下游额济纳旗绿洲的东部大面积的低覆盖度草地向未利用地发生了转移,是其他生态用地面积增加的主因之一.

表 7 黑河流域1990~2023年三生空间面积1) Table 7 Area of the Production Living Ecological Spaces in the Heihe River Basin from 1990 to 2023

图 3 黑河流域1990~2023年三生空间分布 Fig. 3 Distribution of Production Living Ecological Spaces in the Heihe River Basin from 1990 to 2023

从空间转移来看(图 4),农业生产用地增加主要来自未利用地和草地的转入,转入面积分别为785.20 km2和564.91 km2,主要集中在黑河中游张掖与金塔县等一带;工矿企业生产用地主要来自未利用地和草地等生态用地的转入,转入面积分别为248.47 km2和38.92 km2,分布在张掖市、酒泉市以及额济纳旗策克口岸一带;城镇农村生活用地面积主要由耕地和未利用地等转入而增加,在嘉峪关、酒泉市以及张掖市等主城区范围扩张明显;林地生态用地主要由草地和未利用地转入,转入面积分布为120.50 km2和110.78 km2,在上游的祁连县境内较为集中. 此外,受退耕还林等工程影响,由耕地转为林地的面积约12.83 km2. 草地生态用地是黑河流域面积占比排第二的类型,分别由未利用地、林地、耕地以及水域等用地类型转入,在上中下游均有不同程度的分布,主要呈面状和片状展布;水域生态用地的转入主要集中在上游一带,面状或条带状分布特征明显;未利用地等其他生态用地是区内分布最广的用地类型,草地、水域以及林地等生态用地的转入是主要来源,三者迁移面积分别为780.73、385.90和155.63 km2. 整体来说,黑河流域三生空间之间的转移较为频繁且复杂.

数值单位为km2 图 4 黑河流域1990~2023年三生空间转移桑基图 Fig. 4 Sankey diagram of Production Living Ecological Space transitions in the Heihe River Basin from 1990 to 2023

3.2 黑河流域2035年三生空间多情景模拟

在满足模拟精度的前提下,本文基于FLUS模型分别设置生产空间优先、生活空间优先和生态空间优先发展等3种情景,对研究区2035年三生空间时空分布进行预测. 预测结果显示,三生空间分布格局与1990~2023年具相似性(图 5). 同时,为保证各情景下不同用地类型发生转移的概率准确性,其领域权重值通过2020~2035年三生空间面积变化进行赋值(表 3).

图 5 黑河流域2035年三生空间多情景模拟 Fig. 5 Multi-scenario simulation of Production Living Ecological Spaces in the Heihe River Basin in 2035

生产空间优先发展情景下,2035年三生空间用地时空分布特征如下:较2020年,生产空间用地面积呈减少趋势,农业与工矿企业生产空间面积分别减少了15.37 km2和15.34 km2,主要向草地和其他生态用地发生了转移,在中游一带较为集中. 城镇生活空间的面积增加了3.00 km2,而农村生活用地减少18.34 km2,其中生产用地向城镇生活流入,导致其面积增加;而农村生活空间主要流向耕地和未利用地等. 生态空间中林地和草地面积损失较大,分别减少了271.23 km2和800.71 km2,二者之间互有转移,其次是向水域和其他生态用地发生了转移. 水域和其他生态空间用地面积均呈增加趋势,分别增加68.89 km2和1 050.53 km2,主要来自耕地、林地和草地的大量转入. 从区域来看,上游高海拔区域的草地退化较为严重,可能受寒漠化碎石流覆盖的影响[41].

生活空间优先发展情景下,随着经济的快速发展,相较于2020年,城镇生活空间呈显著增加趋势,面积净增85.44 km2,主要来自生产空间和草地生态空间的转入,集中在黑河中游的张掖市、嘉峪关和酒泉市主城区一带. 生产空间合计减少211.70 km2,除向生活空间转入外,还有部分流向草地和其他类生态空间. 该情景下,由于生活空间与其他生态空间面积的扩张,林地和草地生态用地面积在持续减少,分别损失245.31 km2和783.34 km2. 其他生态空间的面积增加了1 041.17 km2,主要来自草地、水域及林地等生态空间和农业生产空间的转入. 特别是,在上游肃南县、中游金塔盆地以及下游额济纳绿洲区均有分布.

生态空间优先发展情景下,2035年黑河流域三生空间分布格局为:限制了生态用地向生产空间的转入,农业生产与工矿企业生产用地的面积分别减少168.64 km2和88.33 km2,转换区域主要集中在中游绿洲盆地内. 生活空间用地面积也呈减少的趋势,其中农村生活用地向农业生产用地转移面积达146.17 km2,呈碎斑状分布在中游农灌区. 受益于对生态用地向其他用地转移的制约,最大程度上降低了草地面积的损失,该情景下林地生态用地面积增加了17.76 km2. 相较于其他情景,草地面积减少216.43 km2,损失最小,主要来自其他生态用地的转入,其次为耕地等生产用地的转入,可能与退耕还林还草以及“三北”防护林等生态工程的持续推进相关.

3.3 黑河流域1990~2023年三生空间生境质量特征

基于InVEST模型的生境质量模块得到黑河流域1990~2023年生境质量指数和生境退化度指数,通过GIS平台实现可视化(图 6). 生境质量指数取值范围为[0, 1],为进一步直观展示生境质量指数在空间上的分布特征,结合徐云等[42]研究成果并结合区域实况,将其划分为[0,0.2]、(0.2,0.4]、(0.4,0.6]、(0.6,0.8]和(0.8,1],并对应为低质量、较低质量、中等质量、较高质量和高质量5个等级.

图 6 1990~2023年生境质量时空分布 Fig. 6 Spatial and temporal distribution of habitat quality from 1990 to 2023

从时间尺度来看,1990~2023年黑河流域HQ均值分别由0.381 5下降至0.378 6,多年HQ均值为0.370 4,年际变化率为-0.013%·a-1,呈明显波动变化,生境质量状况有衰退迹象. 通过生境质量等级面积转移矩阵来看(表 8),低质量区域面积累计转入8 530.57 km2,主要来自较低质量等级的转入;中等质量区域有所增加,主要由较高质量(46.54%)和低质量(39.16%)转入. 较高质量区域主要向中等质量和较低质量区发生转移,整体面积损失44.20 km2;而高质量等级总体呈减少趋势,合计减少741.95 km2.

表 8 黑河流域1990~2023年生境质量等级面积转移矩阵1)/km2 Table 8 Transfer matrix of area of habitat quality grades in the Heihe River Basin from 1990 to 2023/km2

从各等级空间分布来看,高质量等级主要集中分布在上游祁连山保护区内,是重要的水源涵养区,生态本底良好,人为干扰程度较低. 较高等级区主要围绕高质量区外围分布,在上游、中游绿洲灌区以及额济纳旗绿洲河谷和平原较为集中分布. 中等级分布在黑河中游山前冲洪积平原和绿洲区. 较低质量区主要以其他生态空间用地为自然本底,自金塔盆地至黑河下游呈面状展布,是分布范围最广的一类. 而低质量区主要集中在人类活动频繁,城市建设和工业发展密集区. 图 6(i)为1990~2023年生境质量等级变化情况,中游农业区人为活动频繁,对环境干扰强烈,是生境变化较为剧烈的区域;上游在祁连山国家公园的持续保护下,生境质量显著增加区域扩张明显.

1990~2023年黑河流域三生空间生境退化指数见图 7. 生境退化度均值由0.021 9增长至0.031 8,年平均增长率达0.03%. 轻度退化区域在黑河流域上游和下游为主要分布区,从三生空间类型来看,主要以其他生态空间和林草等生态空间用地为主. 较轻度退化区主要集中在轻度退化外围边缘,在过渡带分布较为集中. 而研究区中度退化、较高度退化和高度退化主要分布在中游主要绿洲农业区、经济发达区,包括下游额济纳旗达来呼布镇、东风镇一带. 从退化等级面积变化来看,1990~2023年间除轻度退化区面积减少了8 713.42 km2,其余等级均呈增加趋势,中度至高度退化面积共增加了7 078.05 km2,表现为中游绿洲区和山前冲洪积平原生境退化程度增加,特别是随着经济高速发展和人口数量增长,对于绿洲农业区和其他生态用地的改造加剧,对其生境状况产生直接影响.

图 7 1990~2023年生境退化度时空分布 Fig. 7 Spatial and temporal distribution of habitat degradation degree from 1990 to 2023

基于GIS工具对黑河流域1990~2023年不同生境质量等级的重心空间分布和迁移特征进行分析(图 8). 分析认为不同等级重心均发生了较为明显的空间移动,生境低质量重心向南移动52.77 km,表现为低质量区往中上游绿洲区的转移. 生境较低质量区从黑河下游向中上游空间聚集加剧,其重心向东部迁移了14.47 km. 中等质量区保持相对稳定,其迁移距离最小,仅向北西向迁移4.53 km. 生境质量较高、高等级区向中下游扩张,分别向北西向、北东向迁移17.61 km和21.66 km,也反映出下游生境质量有向好发展的趋势.

图 8 不同生境质量等级重心空间分布 Fig. 8 Spatial distribution of the centers of gravity of different habitat quality grades

3.4 黑河流域2035年生境质量多情景预测

基于FLUS-InVEST模型对黑河流域2035年3种情景下的三生空间生境质量进行预测估算. 图 9(a)~9(c)是2035年各情景下生境质量空间分布特征;图 9(d)~9(f)是以2020年为基准,到2035年各情景下生境质量空间变化分布情况;图 10为2035年各情景下不同质量等级面积. 整体来看,2035年3种情景生境质量的空间分布与1990~2023年具相似性,整体呈现“东南高,西北低”的格局.

图 9 2035年多情景下生境质量分布与时空演变 Fig. 9 Distribution and spatiotemporal evolution of habitat quality under multiple scenarios in 2035

图 10 2035年各情景下不同生境质量等级面积 Fig. 10 Area of different habitat quality grades under various scenarios in 2035

生产空间优先发展情景下[图 9(a)和图 10],相比于2020年,不同生境质量等级在空间位置和面积上均有所变化. 生境质量指数均值为0.376 9,相对下降了1.13%;其生境退化程度有所减缓,下降了4.50%. 主要表现在:低质量区面积减少了1 048.21 km2,主要向较低质量发生转移;较低质量区增加1 951.35 km2,大部分来自较高质量和低质量等级的转入,集中在黑河中上游生态用地的转出;中等级增加区域在黑河中游较为集中,增加82.80 km2;较高质量和高质量区面积均有损失,分别减少了729.72 km2和256.24 km2. 主要以二者的相互转移和较低质量区域的转入为主,在上游的肃南县、祁连县以及中游绿洲区外围表现较为突出.

生活空间优先发展情景下[图 9(b)和图 10],相比于2020年,三生空间的生境质量指数仍呈下降趋势,其均值由0.381 2下降至0.377 1,下降幅度1.08%. 生境退化风险较2020年和生产空间优先情景均有所下降. 在空间转移和面积变化来看:该情景下低质量区面积仍表现为减少,其面积减少2 404.33 km2,主要向较低质量等级发生转移,呈环状展布,可能与局部区域的生态修复相关;生境较低质量区主要受低质量区的大面积转入,整体面积增加3 338.67 km2,在金塔盆地和下游绿洲区一带较集中;中等质量区整体保持稳定;生境较高质量和高质量区面积分别减少了588.21 km2和359.60 km2. 特别是在上游高海拔冰川区、中游绿洲区平原以及下游额济纳旗人类活动区等减少较明显. 在不同地类交界处,各生境等级相互转移也较频繁.

生态空间优先发展情景下[图 9(c)和图 10],2035年生境质量指数均值为0.381 5,较2020年和其他两种情景增幅分别为0.08%、1.22%和1.17%,生境改善较明显. 其生境退化度为0.028,风险指数继续呈下降趋势. 在生态优先前提下,低质量区继续保持减少,主要向较低质量区发生转移,转换区域主要集中在中游一带的山前冲洪积平原,在未利用地与林草生态用地、城郊交互带以及河道两侧转换频繁. 中等质量区分布面积由增转减,较其他情景面积减少明显,面积减少了57.66 km2;生境较高质量区面积虽有减少,但较其他情景减少幅度明显下降,减少面积为211.08 km2;高质量区由减转增,主要由较高质量区转入,也向中等、低质量区方向扩张,其面积增加了234.70 km2. 特别是下游的居延海湿地一带高质量区有所增加,生态改善成效较好.

2020~2035年各情景下生境质量变化等级的空间分布格局具相似性[图 9(d)~9(f)]. 上游和中游是生境质量变化集中区,下游以未利用地等生态空间分布为主,保持相对稳定. 同时,在生活优先发展下生境质量显著减少区明显增加,特别是在城市建设区和农灌区. 而生态优先发展下未显现冲突区域,在上游一带以生境质量明显增加等级为主.

4 讨论 4.1 黑河流域子流域生境质量差异分析

黑河流域上游、中游及下游三生空间类型存在显著差异(图 3). 1990~2023年黑河流域HQ均值下降了0.003,其原因在于生境低质量区扩张,而高质量和较高质量区缩减. 上游林地、草地和水域等生态用地分布集中,祁连山国家公园的建立保障了生态本底相对完整性[43]. 其HQ均值明显高于中下游地区,这与Zhang等[35]的研究结论相一致. 但已有研究[41]认为高海拔冻融地质灾害形势严峻,在一定程度上制约了上游区域生境质量的提升.

中游一带经济发展迅速,人口在绿洲城市不断聚集,生产用地和生活用地扩张明显,地区干扰程度不断加剧[44]. 图 7显示,黑河中游是生境退化较集中的区域,显然与生产空间和生活空间高度开发密切相关[45]. 林草地等生态用地更多向未利用地等转移导致其生境质量变差,这与禹丝思等[46]的研究结论相一致. 下游额济纳绿洲承担着阻隔巴丹吉林沙漠与腾格里沙漠汇合的生态重任. Wang等[47]认为未来黑河下游的生物丰度与植被盖度形势严峻,土地退化势必会阻碍生态环境的恢复并加剧生态风险. 此外,地下水水位埋深变化对河岸胡杨林、柽柳等植被生长具有明显的调控作用[48],水量调度政策出台实施对下游生态状况改善极具现实意义[49].

2035年多情景预测中,生态空间优先产生的结果是生境质量得到好转,对比于生产空间和生活空间优先这两种情景,区域内生境低质量、较低质量区面积呈明显减少或者增幅变缓,生境等级向高质量等级扩张,主要分布在额济纳旗绿洲东部,靠近巴丹吉林沙漠一带,形成防风固沙生态屏障. 这与其他学者认为生态保护情景能够增强区域生态系统服务功能的结论相一致[50]. 同时,生态优先发展也是践行习近平生态文明思想的有力举措[51].

4.2 不足与展望

在研究过程中也存在一定不足. 一是模型驱动因子在选取时尽可能考虑其全面性、合理性和可获取性,以保证模型模拟精度符合要求. 但研究区地理区位特殊,地形地貌差异大,生态功能不同,具点多线杂面广的特点,驱动因子还需进一步补充完善. 二是FLUS模型精度与参数设置高度相关[52],如成本转移原则、领域权重值、未来情景的像元数等,但随着社会经济高速发展,受国家政策调控[53]、区域产业结构优化调整[54]等影响较大,三生空间用地转移频繁,功能结构更加复杂多样,因此,模型参数优化方面还有进步空间,以此来实现区域科学治理的指导目标.

此外,黑河流域三生空间体系规划要因地制宜. 上游是重要的水源涵养区,要以祁连山国家公园为主体,核心区限牧禁牧,避免过多人类干预,确保生态本底质量. 中游是河西走廊重要的农业区和交通枢纽,用水矛盾突出,人类活动频繁,更是遏制巴丹吉林沙漠南侵的生态屏障区. 应当考虑集约高效的水利设施建设,提高用地用水效率,为黑河下游生态用水做好保障. 下游生态本底有限,绿洲区人地矛盾尖锐,生态用地形势严峻. 要建立健全生态风险阈值监测预警机制,充分做好生态保护修复评价工作,为黑河流域生态治理提供支撑.

5 结论

(1)1990~2023年黑河流域三生空间转移频繁且复杂. 其他生态用地是区内分布最广的类型,主要来源是草地、水域以及林地等生态用地的转入. 农业生产空间主要来自未利用地和草地的转入,呈先激增后略减的趋势;工矿企业与城镇农村生活用地的增加与区域经济水平显著正相关,主要来自未利用地和草地的转入;林、草地生态用地面积分别减少了167.04 km2和606.31 km2.

(2)1990~2023年黑河流域生境质量整体呈下降趋势,较高质量和高质量区向中等和较低质量区转移是主因. 两者重心分别向北西向、北东向迁移17.61 km和21.66 km;低质量区重心向南移动52.77 km,有向中上游绿洲区转移的趋势,绿洲区的生境退化程度明显增强.

(3)2035年生产空间优先发展时,林地和草地面积分别减少了271.23 km2和800.71 km2;生活空间优先发展下,城镇生活空间显著增加,林地和草地生态用地面积持续减少,其他生态空间面积则增加了1 041.17 km2;生态空间优先发展时,生产生活用地面积均呈减少趋势,林地面积增加17.76 km2. 相较于其他情景,草地面积减少幅度最低,主要来自其他生态用地和耕地等生产用地的转入.

(4)2035年生产空间优先和生活空间优先情景下,生境质量指数分别下降1.13%和1.08%. 不同地类交界处,生境等级转换频繁,城建区和农灌区的生境质量冲突显著. 生态优先情景下,其生境质量指数较2020年和其他两种情景增幅分别为0.08%、1.22%和1.17%,生境改善明显. 低质量区分布持续减少,高质量区面积由减转增,增加234.70 km2,在上游一带以生境质量明显增加等级为主.

参考文献
[1] 王梓洋, 石培基, 李雪红, 等. 河西走廊地区生境质量对土地利用变化的响应机制及提升路径[J]. 环境科学, 2024, 45(12): 6910-6921.
Wang Z Y, Shi P J, Li X H, et al. Response mechanism and promotion path of habitat quality to land use change in Hexi Corridor Area[J]. Environmental Science, 2024, 45(12): 6910-6921. DOI:10.13227/j.hjkx.202312250
[2] 翟玉鑫, 张飞云, 马丽娜. 基于三生空间的博斯腾湖流域生境质量时空演变及预估[J]. 干旱区地理, 2023, 46(11): 1792-1802.
Zhai Y X, Zhang F Y, Ma L N. Evolution and prediction of habitat quality in the Bosten Lake Basin based on production-living-ecological space[J]. Arid Land Geography, 2023, 46(11): 1792-1802.
[3] 王超越, 郭先华, 陆晨东, 等. 基于FLUS的成渝城市群三生空间演变及模拟预测[J]. 生态科学, 2024, 43(5): 157-168.
Wang C Y, Guo X H, Lu C D, et al. Evolution and simulation of production-living-ecological space in Chengdu-Chongqing urban agglomeration based on FLUS[J]. Ecological Science, 2024, 43(5): 157-168.
[4] 张振龙, 侯琰珍, 孙鸿鹄, 等. 乡镇街道尺度国土"三生"空间功能评价及协调关系[J]. 自然资源学报, 2022, 37(11): 2898-2914.
Zhang Z L, Hou Y Z, Sun H H, et al. Study on the evaluation of the spatial function and coordination relationship of the territorial "production-living-ecological" spaces at the township-street scale[J]. Journal of Natural Resources, 2022, 37(11): 2898-2914.
[5] 武子豪, 祖健, 史云扬, 等. 城市职能视角下"三生"空间的识别与评价——以京津冀城市群为例[J]. 资源科学, 2022, 44(11): 2247-2259.
Wu Z H, Zu J, Shi Y Y, et al. Identification and evaluation of production-living-ecological space from the perspective of urban function: taking the Beijing-Tianjin-Hebei urban agglomeration as an example[J]. Resources Science, 2022, 44(11): 2247-2259.
[6] 潘婷, 王懿祥, 刘宪钊, 等. 雄安新区土地利用变化及其对生态质量的影响[J]. 浙江农林大学学报, 2023, 40(5): 1102-1110.
Pan T, Wang Y X, Liu X Z, et al. Land use change and its impact on ecological quality in Xiong'an New Area[J]. Journal of Zhejiang A&F University, 2023, 40(5): 1102-1110.
[7] Wanjiru C, Nagelkerken I, Rueckert S, et al. Where to fish in the forest? Tree characteristics and contiguous seagrass features predict mangrove forest quality for fishes and crustaceans[J]. Journal of Applied Ecology, 2023, 60(7): 1340-1351. DOI:10.1111/1365-2664.14421
[8] Wang H, Liu W H, Dai G H, et al. The effectiveness of natural reserves from the perspective of habitat quality in the southern section of the Hengduan Mountains, Southwestern China[J]. Journal of Environmental Management, 2024, 363. DOI:10.1016/J.JENVMAN.2024.121380
[9] Zhao L S, Yu W Y, Meng P, et al. InVEST model analysis of the impacts of land use change on landscape pattern and habitat quality in the Xiaolangdi Reservoir area of the Yellow River Basin, China[J]. Land Degradation & Development, 2022, 33(15): 2870-2884.
[10] 孟庆林, 李明玉, 任春颖, 等. 基于HSI模型的吉林省东部地区生境质量动态评价[J]. 国土资源遥感, 2019, 31(3): 140-147.
Meng Q L, Li M Y, Ren C Y, et al. Dynamic assessment of habitat quality in eastern Jilin Province based on HSI model[J]. Remote Sensing for Land & Resources, 2019, 31(3): 140-147.
[11] Yang H X, Zhang H Y, Wang Y R, et al. Urban bird diversity conservation plan based on the MaxEnt model and InVEST model: a case study of Jinan, China[J]. Ecological Indicators, 2025, 174. DOI:10.1016/J.ECOLIND.2025.113463
[12] Araya-Donoso R, Biddy A, Munguía-Vega A, et al. Habitat quality or quantity? Niche marginality across 21 plants and animals suggests differential responses between highland and lowland species to past climatic changes[J]. Ecography, 2024, 2024(9). DOI:10.1111/ECOG.07391
[13] 张廷, 胡玉柱, 胡海辉, 等. 基于PLUS-InVEST模型的哈尔滨市土地利用及生境质量预测[J]. 环境科学, 2024, 45(8): 4709-4721.
Zhang T, Hu Y Z, Hu H H, et al. Prediction of land use and habitat quality in Harbin City based on the PLUS-InVEST model[J]. Environmental Science, 2024, 45(8): 4709-4721. DOI:10.13227/j.hjkx.202308212
[14] 张能能, 刘兆刚. 基于InVEST模型的小兴安岭生态系统服务功能权衡/协同效应[J]. 环境科学, 2025, 46(7): 4628-4640.
Zhang N N, Liu Z G. Trade-offs/synergies of ecosystem services in the Xiaoxing'an Mountains based on the InVEST model[J]. Environmental Science, 2025, 46(7): 4628-4640. DOI:10.13227/j.hjkx.202405314
[15] 左翔, 刘晖. 基于FLUS-InVEST模型的西安市水源涵养功能动态分析与多情景预测[J]. 西北林学院学报, 2025, 40(4): 256-266.
Zuo X, Liu H. Dynamic analysis and multi-scenario prediction of water conservation function of Xi'an City based on FLUS-InVEST model[J]. Journal of Northwest Forestry University, 2025, 40(4): 256-266.
[16] Xu Q, Guo P, Jin M T, et al. Multi-scenario landscape ecological risk assessment based on Markov-FLUS composite model[J]. Geomatics, Natural Hazards and Risk, 2021, 12(1): 1449-1466. DOI:10.1080/19475705.2021.1931478
[17] 袁雪松, 周俊, 胡蓓蓓, 等. 基于FLUS模型粤港澳大湾区"三生空间"多情景模拟预测[J]. 地理科学, 2023, 43(3): 564-574.
Yuan X S, Zhou J, Hu B B, et al. Multi-scenario simulation and prediction of ecological-productive-living spaces in Guangdong-Hong Kong-Macao Greater Bay Area based on FLUS model[J]. Scientia Geographica Sinica, 2023, 43(3): 564-574.
[18] Yang J, Xie B P, Zhang D G, et al. Habitat quality assessment and multi-scenario prediction of the Gansu-Qinghai section of the Yellow River Basin based on the FLUS-InVEST model[J]. Frontiers in Ecology and Evolution, 2023, 11. DOI:10.3389/FEVO.2023.1228558
[19] 王生棠, 盖迎春, 王少昆, 等. 黑河流域植物物种多样性对土地利用变化的响应[J]. 中国沙漠, 2022, 42(6): 221-232.
Wang S T, Gai Y C, Wang S K, et al. Responses of plant species diversity to land use change: a case study in the Heihe River Basin[J]. Journal of Desert Research, 2022, 42(6): 221-232.
[20] 白羽萍, 王佳斌, 刘恋, 等. 黑河流域生态系统服务价值预测及权衡分析——以张掖市为例[J]. 生态学报, 2022, 42(13): 5187-5199.
Bai Y P, Wang J B, Liu L, et al. Prediction and trade-off analysis of ecosystem service value in Heihe River Basin: a case study in Zhangye City[J]. Acta Ecologica Sinica, 2022, 42(13): 5187-5199.
[21] 王嘉伟, 董国涛, 蒋晓辉, 等. 生态输水背景下黑河流域生态系统服务的时空演变特征及驱动力研究[J]. 干旱区地理, 2025, 48(8): 1406-1420.
Wang J W, Dong G T, Jiang X H, et al. Spatiotemporal evolution characteristics and driving forces of ecosystem services in the Heihe River Basin under the context of ecological water conveyance[J]. Arid Land Geography, 2025, 48(8): 1406-1420.
[22] 杨鸿魁, 寇明珠, 祁万强, 等. 基于LUCC变化的黑河流域碳储量时空模拟预测及驱动力分析[J]. 环境科学, 2025, 46(11): 7070-7082.
Yang H K, Kou M Z, Qi W Q, et al. Spatiotemporal simulation and driving force analysis of carbon storage in Heihe River Basin based on LUCC changes[J]. Environmental Science, 2025, 46(11): 7070-7082. DOI:10.13227/j.hjkx.202411234
[23] 许澳康, 吉天琪, 石晶. 甘肃黑河流域生态网络构建与优化[J]. 中国沙漠, 2025, 45(2): 61-70.
Xu A K, Ji T Q, Shi J. Construction and optimization of ecological network of the Heihe River Basin (Gansu section)[J]. Journal of Desert Research, 2025, 45(2): 61-70.
[24] Yu J S, Zhou J, Zhao J, et al. Agroecological risk assessment based on coupling of water and land resources—a case of Heihe River Basin[J]. Land, 2023, 12(4). DOI:10.3390/LAND12040794
[25] 吴绮琦, 蒙吉军. 近40年来黑河流域三生空间格局演变及其驱动因素分析[J]. 北京大学学报(自然科学版), 2023, 59(6): 970-980.
Wu Q Q, Meng J J. Analysis of the evolution and driving factors of production-living-ecological space pattern in the Heihe River Basin from 1980 to 2020, China[J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2023, 59(6): 970-980.
[26] 邹易, 赵一阳, 蒙吉军. 气候变化背景下黑河流域"三生"空间转型对生态环境质量的影响[J]. 北京大学学报(自然科学版), 2024, 60(3): 528-538.
Zou Y, Zhao Y Y, Meng J J. Impacts of production-living-ecological space transformation in the Heihe River Basin on ecological environment quality in the context of climate change[J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2024, 60(3): 528-538.
[27] 窦寒梅, 赵锐锋, 陈喜东, 等. 西北干旱区国土空间生态修复优先区识别——以黑河流域张掖市为例[J]. 应用生态学报, 2024, 35(2): 469-479.
Dou H M, Zhao R F, Chen X D, et al. Identification of priority areas for territorial space ecological restoration in arid area of Northwest China: a case study of Zhangye City in Heihe River Basin[J]. Chinese Journal of Applied Ecology, 2024, 35(2): 469-479.
[28] 刘继来, 刘彦随, 李裕瑞. 中国"三生空间"分类评价与时空格局分析[J]. 地理学报, 2017, 72(7): 1290-1304.
Liu J L, Liu Y S, Li Y R. Classification evaluation and spatial-temporal analysis of "production-living-ecological" spaces in China[J]. Acta Geographica Sinica, 2017, 72(7): 1290-1304.
[29] 年雁云. 黑河流域河流分布数据集[DB/OL]. 国家冰川冻土沙漠科学数据中心. https://cstr.cn/CSTR:11738.11.ncdc.Westdc.2020.465, 2020-12-18.
[30] 李雪, 李文, 高宇. 基于FLUS-InVEST模型的大庆市土地覆盖模拟与碳储量评估[J]. 环境科学, 2024, 45(10): 5983-5993.
Li X, Li W, Gao Y. Land Cover simulation and carbon storage assessment in Daqing City based on FLUS-InVEST model[J]. Environmental Science, 2024, 45(10): 5983-5993. DOI:10.13227/j.hjkx.202312007
[31] 李小军, 车良革, 胡宝清. FLUS-InVEST模型的北海市生态系统碳储量时空差异分析[J]. 测绘通报, 2023(6): 117-123, 183.
Li X J, Che L G, Hu B Q. Spatio-temporal difference analysis of carbon storage in Beihai secosystem based on FLUS-InVEST models[J]. Bulletin of Surveying and Mapping, 2023(6): 117-123, 183.
[32] 国务院. 国务院关于黑河流域近期治理规划的批复[EB/OL]. https://www.gov.cn/gongbao/content/2001/content_60988.htm, 2001-08-03.
[33] 阿拉善盟生态环境局. 阿拉善盟黑河流域水生态环境保护"十四五"规划[R]. 阿拉善盟: 阿拉善盟生态环境局, 2022.
[34] Zhou S H, Qu Y L, Wang Y X, et al. Ecosystem service bundles under SSP-RCP and local scenarios: a pathway to comprehensive spatial planning for sustainability[J]. Resources, Environment and Sustainability, 2025, 20. DOI:10.1016/J.RESENV.2025.100211
[35] Zhang L G, Hu N K. Spatial variation and terrain gradient effect of ecosystem services in Heihe River Basin over the past 20 years[J]. Sustainability, 2021, 13(20). DOI:10.3390/su132011271
[36] Yan Y X, Liu X H, Lin T, et al. Analysis of the Interrelationships and drivers of ecosystem services in the Heihe River Basin[J]. Sustainability, 2025, 17(5). DOI:10.3390/su17051942
[37] 魏婷婷, 张尧, 马其浩, 等. 基于FLUS-Markov模型的宝鸡峡灌区土地利用变化预测[J]. 排灌机械工程学报, 2025, 43(6): 604-611.
Wei T T, Zhang Y, Ma Q H, et al. Prediction of land use change in Baojixia irrigation district based on FLUS-Markov model[J]. Journal of Drainage and Irrigation Machinery Engineering, 2025, 43(6): 604-611.
[38] 韦绍音, 陆汝成, 林晓楠, 等. 基于土地利用转型的广西边境地区生态服务功能交叉敏感性分析[J]. 水土保持研究, 2022, 29(3): 308-316.
Wei S Y, Lu R C, Lin X N, et al. Cross-sensitivity analysis of ecological service function in Guangxi border region based on land use transition[J]. Research of Soil and Water Conservation, 2022, 29(3): 308-316.
[39] 杜怀玉, 俞金凤, 张媛, 等. 石羊河流域多情景土地利用优化及碳储量评估[J]. 环境科学, 2024, 45(7): 4164-4176.
Du H Y, Yu J F, Zhang Y, et al. Land use optimization and carbon reserve assessment in the Shiyang River Basin[J]. Environmental Science, 2024, 45(7): 4164-4176. DOI:10.13227/j.hjkx.202308015
[40] 张艳, 江海云, 杨维新. 渭河流域"三生"空间格局演变特征及其驱动因素[J]. 地球科学与环境学报, 2025, 47(4): 718-732.
Zhang Y, Jiang H Y, Yang W X. Evolution characteristics and driving factors of production-living-ecological space pattern in Weihe River Basin, China[J]. Journal of Earth Sciences and Environment, 2025, 47(4): 718-732.
[41] 杨鸿魁, 庞文龙, 李有三, 等. 祁连山黑河源区央隆草原一带冻融地质作用及其生态环境问题的探讨[J]. 地质灾害与环境保护, 2023, 34(1): 95-100.
Yang H K, Pang W L, Li Y S, et al. Discussion on the geological process of freeze-thaw and its eco-environmental problems of Yanglong grassland in Heihe source region of Qilian mountains[J]. Journal of Geological Hazards and Environment Preservation, 2023, 34(1): 95-100.
[42] 徐云, 高明坤, 张增信. 基于InVEST模型的长江中下游地区生境质量变化研究[J]. 水土保持研究, 2024, 31(5): 355-364.
Xu Y, Gao M K, Zhang Z X. Land use change and its impact on habitat quality in the middle and lower reaches of the Yangtze River based on InVEST model[J]. Research of Soil and Water Conservation, 2024, 31(5): 355-364.
[43] 田新苗, 冯起. 祁连山国家公园区域生态保护与经济发展路径探究[J]. 中国科学院院刊, 2025, 40(6): 1103-1112.
Tian X M, Feng Q. Exploration on path of ecological protection and economic development in Qilian Mountain National Park region[J]. Bulletin of Chinese Academy of Sciences, 2025, 40(6): 1103-1112.
[44] 姚礼堂, 张学斌, 周亮, 等. "山地-绿洲-荒漠"复合系统土地利用变化的生态系统服务权衡与协同效应——以张掖市为例[J]. 生态学报, 2022, 42(20): 8138-8151.
Yao L T, Zhang X B, Zhou L, et al. Ecosystem service tradeoffs and synergies effects of land use change in mountain-oasis-desert complex system: a case study of Zhangye City[J]. Acta Ecologica Sinica, 2022, 42(20): 8138-8151.
[45] 迟宏旭, 陈宇, 杨旭, 等. 东北三省土地利用及生境质量的时空演变效应[J]. 环境科学与技术, 2023, 46(6): 187-196.
Chi H X, Chen Y, Yang X, et al. The spatio-temporal evolution effect of land use function change and habitat quality in the three northeastern provinces of China[J]. Environmental Science & Technology, 2023, 46(6): 187-196.
[46] 禹丝思, 蔡尚书, 张生瑞, 等. 松嫩鹤乡国家公园创建潜力区土地生境质量演化及多情景预测[J]. 林业科学, 2025, 61(7): 208-219.
Yu S S, Cai S S, Zhang S R, et al. Evolution and multi-scenario prediction of land habitat quality in the potential area of the proposed Songnen Plain Crane Homeland National Park[J]. Scientia Silvae Sinicae, 2025, 61(7): 208-219.
[47] Wang C L, Jiang Q O, Shao Y Q, et al. Ecological environment assessment based on land use simulation: a case study in the Heihe River Basin[J]. Science of the Total Environment, 2019, 697. DOI:10.1016/j.scitotenv.2019.133928
[48] 张晓龙, 周继华, 来利明, 等. 黑河下游胡杨群落多样性沿河岸距离的变化特征[J]. 生态环境学报, 2021, 30(10): 1952-1960.
Zhang X L, Zhou J H, Lai L M, et al. Variation characteristics of Populus euphratica community diversity along the downstream riverbanks of the Heihe River[J]. Ecology and Environmental Sciences, 2021, 30(10): 1952-1960.
[49] 郭瑞琪. 科学配置水资源、推进东居延海生态治理——黑河下游生态保护和高质量发展方案研究[J]. 灌溉排水学报, 2020, 39(S2): 136-139.
Guo R Q. Scientific allocation of water resources and promotion of ecological management of Dongjuyanhai: study on ecological protection and high quality development scheme in the Lower reaches of Heihe River[J]. Journal of Irrigation and Drainage, 2020, 39(S2): 136-139.
[50] Li W, Chen X, Zheng J H, et al. A multi-scenario simulation and dynamic assessment of the ecosystem service values in key ecological functional areas: a case study of the Sichuan province, China[J]. Land, 2024, 13(4). DOI:10.3390/land13040468
[51] 张厚美, 刘娟, 蒙俊霖. 坚定不移走生态优先绿色发展之路——以习近平新时代中国特色社会主义思想引领青川生态文明建设[J]. 环境保护, 2023, 51(7): 48-50.
Zhang H M, Liu J, Meng J L. Unswervingly take the road of ecological priority green development—the Qingchuan practice of the construction of ecological civilization leading by Xi Jinping Thought on Socialism with Chinese Characteristics for a New Era[J]. Environmental Protection, 2023, 51(7): 48-50.
[52] Shi G, Chen C, Cao Q C, et al. Spatiotemporal dynamics and prediction of habitat quality based on land use and cover change in Jiangsu, China[J]. Remote Sensing, 2024, 16(22). DOI:10.3390/RS16224158
[53] 王乐, 熊昌盛, 田宇, 等. 中国"三生"空间功能异质性评价、时空演变及分区调控[J]. 农业工程学报, 2024, 40(10): 265-275.
Wang L, Xiong C S, Tian Y, et al. Spatial heterogeneity in functional evaluation, spatio-temporal evolution and zoning schemes of production-living-ecological space in China[J]. Transactions of the Chinese Society of Agricultural Engineering, 2024, 40(10): 265-275.
[54] 王成新, 于雷, 吕红迪. 城市环境空间管控体系探索与思考[J]. 环境科学与管理, 2019, 44(6): 184-189.
Wang C X, Yu L, Lv H D. Research on urban environmental space control system[J]. Environmental Science and Management, 2019, 44(6): 184-189.