生态与环境

基于改进遥感生态指数的准东地区生态环境质量时空变化分析

  • 邓文彬 ,
  • 宋森 ,
  • 易红梅
展开
  • 新疆大学建筑工程学院, 新疆 乌鲁木齐 830046
宋森(2002-),男,硕士研究生,主要从事生态遥感等方面的研究. E-mail:

邓文彬(1977-),男,博士,教授,主要从事变形监测与灾害预警、摄影测量与遥感等方面的教学与研究. E-mail:

收稿日期: 2025-04-28

  修回日期: 2025-07-22

  网络出版日期: 2026-03-11

基金资助

国家自然科学基金项目(52404188)

新疆维吾尔自治区自然科学基金项目(2024D01C244)

新疆重点研发项目(2022B03033-1)

新疆维吾尔自治区天池英才项目资助

Spatiotemporal dynamics of ecological quality in the Zhundong region based on an arid modified remote sensing ecological index model

  • Wenbin DENG ,
  • Sen SONG ,
  • Hongmei YI
Expand
  • School of Architectural Engineering, Xinjiang University, Urumqi 830046, Xinjiang, China

Received date: 2025-04-28

  Revised date: 2025-07-22

  Online published: 2026-03-11

摘要

干旱区资源开发与生态保护的矛盾日益突出,科学量化新疆准东煤炭基地生态损伤机制,对破解干旱区资源开发与保护矛盾、保障国家能源安全与区域可持续发展具有关键意义。在此背景下以煤炭资源富集的新疆准东地区为对象,针对其煤炭开发引发的土地沙化、盐渍化及PM2.5污染等生态问题,提出改进型遥感生态指数(Arid modified remote sensing ecological index,ARSEI),通过谷歌地球引擎,使用Landsat时序数据生成了ARSEI,并利用最优参数地理探测器等工具,揭示了2000—2023年准东地区生态环境质量的时空演化特征。结果表明:(1) 2000—2023年准东地区ARSEI整体呈“先降后升”趋势:2000年为峰值0.368,2018年下降至谷底0.225,2023年回升至0.289;空间分布上,生态环境质量等级以“差”与“较差”为主,整体表现为南优北劣。(2) 生态环境质量普遍下降,其中显著下降区域占83.5%,与中北部地区、中度及以上生态脆弱区高度重合;而生态环境质量上升区仅占4.6%,且集中于南部生态非脆弱区的农田与城镇周边。(3) 驱动力分析表明,人类活动强度、蒸散发量与气温为主导单因子,解释力均大于0.35;交互作用后,气温∩人类活动强度的解释力最强,解释力均大于0.7,是生态环境变化的关键驱动因素。

本文引用格式

邓文彬 , 宋森 , 易红梅 . 基于改进遥感生态指数的准东地区生态环境质量时空变化分析[J]. 干旱区地理, 2025 , 48(12) : 2073 -2086 . DOI: 10.12118/j.issn.1000-6060.2025.236

Abstract

The Zhundong region in Xinjiang, a typical arid and semi-arid area, is a vital national energy base with abundant coal resources. However, intensive coal mining and related industrial activities have caused severe ecological issues, including land desertification, salinization, and PM2.5 pollution, intensifying the conflict between resource exploitation and ecological preservation. To accurately evaluate the ecological environment quality (EEQ) of this region, an arid modified remote sensing ecological index (ARSEI) is proposed herein. The ARSEI improves upon the traditional remote sensing ecological index by incorporating indicators for salinity, desertification, and air pollution, besides greenness, wetness, and heat. Using the Google Earth Engine platform, Landsat 5/8 images from 2000 to 2023 were processed to generate the ARSEI. The first principal component derived via principal component analysis was used to construct the ARSEI, with the direction of the wetness component applied to resolve the eigenvector ambiguity. Temporal trends were analyzed using Sen’s slope estimator and Mann-Kendall trend test, and ecological vulnerability was assessed based on the sensitivity and adaptability of the ARSEI. The optimal parameter-based geographical detector was employed to identify key driving factors and their interactions. Results indicated that the ARSEI of the Zhundong region showed a trend of initial decline followed by a partial recovery from 2000 to 2023, with its mean value decreasing from 0.368 in 2000 to 0.225 in 2018, before increasing to 0.289 in 2023. Spatially, the EEQ was generally poor, with over 60% of the area classified as “poor” or “relatively poor”, exhibiting a clear north-south gradient where the southern part exhibited better ecological conditions. Trend analysis revealed that 83.5% of the region experienced significant ecological degradation, whereas only 4.6% showed improvement, mainly in the southern agricultural and urban areas with low ecological vulnerability. Single-factor detection identified human activity intensity (q>0.6) and evapotranspiration (q>0.35) as the primary drivers. The interaction between temperature and human activity intensity had the strongest explanatory power (q>0.7), indicating nonlinear enhancement effects. In summary, large-scale resource development has exerted considerable pressure on the already fragile ecosystem. Although ecological restoration measures since 2019 have led to a partial recovery, the northern mining areas remain highly vulnerable due to water scarcity and soil degradation. The ARSEI model demonstrates enhanced applicability in arid resource-based regions and offers a scientific basis for ecological protection and sustainable development strategies.

干旱半干旱地区占中国国土面积的47%,生态系统结构简单、抗干扰能力弱,水资源短缺与植被覆盖低下导致其生态环境具有先天脆弱性。随着干旱区资源开发强度不断加大,开发与生态系统承载力之间的矛盾日益尖锐,土地荒漠化、盐渍化、草场退化等生态环境问题的不断加剧,使区域可持续发展面临严峻挑战[1]。因此,准确、定量地刻画生态环境质量时空演化并识别关键驱动因素,已成为协调资源开发与生态保护的前提[2]
对生态系统准确定量评价是一项复杂而困难的任务。传统生态评价多采用单一指标,难以综合反映复杂生态系统状态[3-4];2006年我国国家环保局发布的《生态环境状况评价技术规范》中提出了生态环境状况指数[5],其虽整合多个要素,但存在权重主观、指标冗余以及无法可视化等局限。2013年徐涵秋提出了遥感生态指数[6](Remote sensing ecological index,RSEI),通过主成分分析(Principal component analysis,PCA)集成了绿度(Normalized differential vegetation index,NDVI)、湿度(Wetness,WET)、热度(Land surface temperature,LST)和干度(Normalized difference built-up and soil index,NDBSI)4个生态指标,能够更全面地反映区域生态环境质量的变化。并且RSEI基于自然因子,无需人为赋权,评价结果更加客观合理,且能够实现空间可视化,被广泛用于我国生态环境质量评价领域[7-8],然而,RSEI在城市生态环境评价中表现优异,但极端环境下却面临指标适用性不足、无法提取目标生态特征等局限。为了更好地适用于特定区域,国内外学者从指标重构与算法优化进行探索。在指标重构方面,王杰等[9]针对干旱区生态环境,引入盐度与土地退化度构建干旱遥感生态指数,显著提升了沙漠地区的识别精度;马鹏飞等[10]考虑矿区的大气污染问题,引入PM2.5差值指数,并对矿区进行长时序监测,并证实其可以较好反映矿区生态。在算法优化方面,汤从沧等[11]指出PCA在塔里木河流域存在信息损失,所以采用熵值法加权构建改进遥感生态指数,对塔里木河流域的生态环境进行了准确的评价;Zhang等[12]通过湿度分量稳定性较好的特点,解决PCA过程中特征向量方向歧义问题,避免结果倒置。尽管上述改进提升了RSEI的适用性与精度,但已有研究多聚焦指数构建本身,缺乏对RSEI结果的挖掘。同时,针对干旱资源开发区与其周边生态环境变化的评价较少。
准东地区位于新疆维吾尔自治区昌吉回族自治州,是典型的干旱区生态,生态系统先天脆弱。但该区域拥有丰富的煤炭资源,是国家级煤炭基地,预计存量可达3.9×1012 t,是新疆重要的能源核心区和“西气东输”“疆电外送”的重要基地[13]。然而,大规模露天煤炭开采及煤化工产业发展加剧了土地沙化、盐碱化及大气污染等问题[14],形成了开发与保护的典型矛盾体,其生态环境质量制约准东地区可持续发展。尽管近年来生态修复力度加强,但生态系统的脆弱性与资源开发的矛盾仍未根本缓解。基于此,为精准刻画准东地区生态环境演化特征,本文构建改进型干旱遥感生态指数(Arid modified remote sensing ecological index,ARSEI)。在RSEI基础上,以盐度(Modified salinity index,MSI)、沙度(Desertification modified index,DMI)替代NDBSI,表征干旱区土地退化,PM2.5差值指数(Difference index,DI),反应矿区工业污染,并与NDVI、WET、LST耦合形成多维指标,实现对准东地区的综合表征。并依托谷歌地球引擎(Google Earth Engine,GEE)平台处理2000—2023年的Landsat时序数据实现ARSEI,并将ARSEI生成的生态脆弱性引入分析过程并结合Sen+Mann-Kendall趋势分析,对准东地区生态变化进行更加深入的讨论。最后,借助最优地理探测器分析准东地区生态质量空间变化的驱动因素,解释生态质量变化与区域发展的内在联系,以期助力建设资源型区域实现可持续发展。

1 数据与方法

1.1 研究区概况

研究区位于新疆维吾尔自治区昌吉回族自治州(图1),地理位置介于88°30′~90°30′E,44°00′~45°30′N之间,总面积为50148 km2,横跨吉木萨尔、奇台、木垒3县。该区域属典型的温带大陆性干旱气候,年均气温7.6 ℃,年降水量低于150 mm,而年蒸发量则大于2000 mm。地势呈南北高、中间低的特征,植被覆盖度低,植被类型以耐旱灌木与低覆盖草地为主。但研究区煤炭资源富集,为国家级经济技术开发区,2023年煤炭产能超2×108 t,贡献昌吉回族自治州35%的GDP,且为“西电东送”的重要基地,已建火电装机25 GW。在研究区快速发展的同时,其资源开发、经济发展与生态环境保护之间的矛盾日益凸显[15]
图1 研究区概况

Fig. 1 Overview of the study area

1.2 数据来源与预处理

数据来源方面,遥感影像基于GEE平台获取Landsat 5、Landsat 8卫星影像(2000—2023年),选择每年6—9月云量低于10%的影像,经去云掩膜等处理后提取地表参数,其他数据来源见表1。此外,为确保高精度与时空一致性,生成ARSEI时,使用30 m的Landsat 5/8数据,而在进行分析时,将所有数据统一为WGS 1984坐标系,空间分辨率重采样至1 km。
表 1 数据信息与来源

Tab. 1 Data information and source

数据类型 数据名称 分辨率 来源
遥感影像 Landsat 5 30 m GEE(LANDSAT/LT05/C02/T1_L2)
Landsat 8 30 m GEE(LANDSAT/LC08/C02/T1_L2)
自然数据 高程 90 m 地理空间数据云(https://www.gscloud.cn/
坡度 90 m 地理空间数据云(https://www.gscloud.cn/
坡向 90 m 地理空间数据云(https://www.gscloud.cn/
气温 1 km NOAA(https://www.ncei.noaa.gov
降水 1 km 国家青藏高原科学数据中心(https://data.tpdc.ac.cn/home
蒸散发量 1 km 国家青藏高原科学数据中心(https://data.tpdc.ac.cn/home
经济数据 人口密度 1 km LandScan Global数据集(https://landscan.ornl.gov/
夜间灯光 500 m 国家地球系统科学数据中心(http://www.geodata.cn
土地利用类型 1 km 中国科学院地理科学与资源研究所(http://www.resdc.cn

1.3 研究方法

1.3.1 改进遥感生态指数

ARSEI是在RSEI的基础上,考虑到研究区地广人稀特点去除NDBSI[9],并因研究区土地沙化、盐碱化的生态特点加入DMI和MSI指标,再将DI引入以反映准东经济开发区产生的污染对生态的影响[16],进而使用NDVI、WET、LST、MSI、DMI与DI这6个指标,从绿度、湿度、热度、土地盐碱化、土地沙化与空气污染6个方面对生态环境质量进行全方面因地制宜的评价。其相关指标计算方法如表2所示。
表 2 ARSEI模型生态指标计算公式

Tab. 2 Calculation formula for ARSEI model ecological indicators

遥感参数 计算公式 参数解释
NDVI[17] $\mathrm{N}\mathrm{D}\mathrm{V}\mathrm{I}=\left({\rho }_{\mathrm{N}\mathrm{I}\mathrm{R}}-{\rho }_{\mathrm{R}\mathrm{e}\mathrm{d}}\right)/\left({\rho }_{\mathrm{N}\mathrm{I}\mathrm{R}}+{\rho }_{\mathrm{R}\mathrm{e}\mathrm{d}}\right)$ NDVI为归一化植被指数;${\rho }_{\mathrm{N}\mathrm{I}\mathrm{R}}$${\rho }_{\mathrm{R}\mathrm{e}\mathrm{d}}$分别为近红外、红波段的反射率
WET[18] $\begin{array}{l}\mathrm{W}\mathrm{E}{\mathrm{T}}_{5}=0.0315{\rho }_{\mathrm{B}\mathrm{l}\mathrm{u}\mathrm{e}}+0.2021{\rho }_{\mathrm{G}\mathrm{r}\mathrm{e}\mathrm{e}\mathrm{n}}+0.3102{\rho }_{\mathrm{R}\mathrm{e}\mathrm{d}}+0.1594{\rho }_{\mathrm{N}\mathrm{I}\mathrm{R}}-\\ 0.6806{\rho }_{\mathrm{S}\mathrm{W}\mathrm{I}\mathrm{R}1}-0.6109{\rho }_{\mathrm{S}\mathrm{W}\mathrm{I}\mathrm{R}2}\end{array}$
$\begin{array}{l}\mathrm{W}\mathrm{E}{\mathrm{T}}_{8}=0.1511{\rho }_{\mathrm{B}\mathrm{l}\mathrm{u}\mathrm{e}}+0.1972{\rho }_{\mathrm{G}\mathrm{r}\mathrm{e}\mathrm{e}\mathrm{n}}+0.3283{\rho }_{\mathrm{R}\mathrm{e}\mathrm{d}}+0.3407{\rho }_{\mathrm{N}\mathrm{I}\mathrm{R}}-\\ 0.7117{\rho }_{\mathrm{S}\mathrm{W}\mathrm{I}\mathrm{R}1}-0.4559{\rho }_{\mathrm{S}\mathrm{W}\mathrm{I}\mathrm{R}2}\end{array}$
WET为湿度指数;$\mathrm{W}\mathrm{E}{\mathrm{T}}_{5}$$\mathrm{W}\mathrm{E}{\mathrm{T}}_{8}$分别为Landsat 5与Landsat 8影像下的湿度计算方法;${\rho }_{\mathrm{B}\mathrm{l}\mathrm{u}\mathrm{e}}$${\rho }_{\mathrm{G}\mathrm{r}\mathrm{e}\mathrm{e}\mathrm{n}}$${\rho }_{\mathrm{S}\mathrm{W}\mathrm{I}\mathrm{R}1}$${\rho }_{\mathrm{S}\mathrm{W}\mathrm{I}\mathrm{R}2}$分别为蓝、绿、短波红外1、短波红外2波段的反射率
LST[19] $\mathrm{L}\mathrm{S}\mathrm{T}={K}_{2}/\left(\frac{{K}_{1}}{T}+1\right)$
${L}_{\lambda }=\left[\epsilon T+\left(1-\epsilon \right){L}_{\mathrm{d}\mathrm{o}\mathrm{w}\mathrm{n}}\right]+{L}_{\mathrm{u}\mathrm{p}}$
$T=\left[{L}_{\lambda }-{L}_{\mathrm{u}\mathrm{p}}-\tau \left(1-\epsilon \right){L}_{\mathrm{d}\mathrm{o}\mathrm{w}\mathrm{n}}\right]/\tau \epsilon $
LST为地表温度;${L}_{\lambda }$T分别为热红外、黑体热辐射亮度值;$\epsilon $为地表比辐射率;$\tau $为大气透过率;Lup${L}_{\mathrm{d}\mathrm{o}\mathrm{w}\mathrm{n}}$为大气上行、下行辐射亮度;K1K2为传感器定标参数
MSI[20] $\mathrm{S}\mathrm{I}=\left({\rho }_{\mathrm{B}\mathrm{l}\mathrm{u}\mathrm{e}}\times {\rho }_{\mathrm{R}\mathrm{e}\mathrm{d}}\right)/2$
$\mathrm{M}\mathrm{S}\mathrm{A}\mathrm{V}\mathrm{I}=\frac{(2{\rho }_{\mathrm{N}\mathrm{I}\mathrm{R}}+1)-\sqrt{(2{\rho }_{\mathrm{N}\mathrm{I}\mathrm{R}}{+1)}^{2}-8({\rho }_{\mathrm{N}\mathrm{I}\mathrm{R}}-{\rho }_{\mathrm{R}\mathrm{e}\mathrm{d}})}}{2}$
$\mathrm{M}\mathrm{S}\mathrm{I}=\sqrt{{\left({M}_{\mathrm{N}\mathrm{I}}-{M}_{\mathrm{m}\mathrm{a}\mathrm{x}}\right)}^{2}+\mathrm{S}{\mathrm{I}}_{\mathrm{N}\mathrm{I}}^{2}}$
MSI为土壤盐渍化指数;MSAVI为修改型土壤调整植被指数;SI为盐分指数;MNI${M}_{\mathrm{m}\mathrm{a}\mathrm{x}}$分别为MSAVI的归一化与最大值;$\mathrm{S}{\mathrm{I}}_{\mathrm{N}\mathrm{I}}$为SI的归一化
DMI[11] $\mathrm{a}\mathrm{l}\mathrm{b}\mathrm{e}\mathrm{d}\mathrm{o}=0.356\times {\rho }_{\mathrm{B}\mathrm{l}\mathrm{u}\mathrm{e}}+0.130{\rho }_{\mathrm{R}\mathrm{e}\mathrm{d}}+0.085{\rho }_{\mathrm{S}\mathrm{W}\mathrm{I}\mathrm{R}1}+0.072{\rho }_{\mathrm{S}\mathrm{W}\mathrm{I}\mathrm{R}2}-0.0018$
$\mathrm{S}\mathrm{F}\mathrm{I}=\left({\rho }_{\mathrm{S}\mathrm{W}\mathrm{I}\mathrm{R}1}-{\rho }_{\mathrm{B}\mathrm{l}\mathrm{u}\mathrm{e}}\right)/\left(200-{\rho }_{\mathrm{S}\mathrm{W}\mathrm{I}\mathrm{R}1}-{\rho }_{\mathrm{S}\mathrm{W}\mathrm{I}\mathrm{R}2}\right)$
$\mathrm{D}\mathrm{M}\mathrm{I}=\sqrt{{\left({M}_{\mathrm{N}\mathrm{I}}-{M}_{\mathrm{m}\mathrm{a}\mathrm{x}}\right)}^{2}+({a}_{\mathrm{N}\mathrm{I}}-{a}_{\mathrm{m}\mathrm{i}\mathrm{n}}{)}^{2}+({S}_{\mathrm{N}\mathrm{I}}-{S}_{\mathrm{m}\mathrm{i}\mathrm{n}}{)}^{2}}$
albedo为地表反射率;SFI为土壤肥力指数;DMI为土壤沙化指数;SNI${S}_{\mathrm{m}\mathrm{i}\mathrm{n}}$分别为SI的归一化与最小值;aNI${a}_{\mathrm{m}\mathrm{i}\mathrm{n}}$分别为albedo的归一化与最小值
DI[21] $\mathrm{D}\mathrm{I}={\rho }_{\mathrm{R}\mathrm{e}\mathrm{d}}-{\rho }_{\mathrm{N}\mathrm{I}\mathrm{R}}$ DI为PM2.5差值指数

1.3.2 构建ARSEI

构建ARSEI需要结合NDVI、WET、LST、MSI、DMI和DI,由于这6个指标量纲不同,所以在进行结合之前需要对各项指标进行归一化,避免各项指标失衡。归一化公式为:
$\mathrm{N}{\mathrm{I}}_{n}=({I}_{n}-{I}_{\mathrm{m}\mathrm{i}\mathrm{n}})/({I}_{\mathrm{m}\mathrm{a}\mathrm{x}}-{I}_{\mathrm{m}\mathrm{i}\mathrm{n}})$
式中:NIn为目标归一化后的结果;In为原始数据值;${I}_{\mathrm{m}\mathrm{i}\mathrm{n}}$${I}_{\mathrm{m}\mathrm{a}\mathrm{x}}$分别为In的最小值和最大值。
本文采用通过主成分分析法得到的第一主成分(PC1)为ARSEI的初始值,不仅避免主观确定权重还最大程度上表示出原始变量的信息。但由于主成分分析方法特征向量的非唯一性,可能会导致相反的结果[22],参考Zhang等[12]的研究,选择使用受季节性变化影响较小的WET分量来判断主成分PC1的特征向量方向,使得改进后的模型无论输入指标如何变化,最终的结果都是合理的。计算公式如下:
$\begin{array}{l}\mathrm{A}\mathrm{R}\mathrm{S}\mathrm{E}{\mathrm{I}}_{0}=\\ \left\{\begin{array}{l}\mathrm{P}\mathrm{C}1\left[f\left(\begin{array}{l}\mathrm{N}\mathrm{D}\mathrm{V}{\mathrm{I}}_{\mathrm{N}\mathrm{I}},\mathrm{W}\mathrm{E}{\mathrm{T}}_{\mathrm{N}\mathrm{I}},\mathrm{L}\mathrm{S}{\mathrm{T}}_{\mathrm{N}\mathrm{I}},\\ \mathrm{M}\mathrm{S}{\mathrm{I}}_{\mathrm{N}\mathrm{I}},\mathrm{D}\mathrm{M}{\mathrm{I}}_{\mathrm{N}\mathrm{I}},\mathrm{D}{\mathrm{I}}_{\mathrm{N}\mathrm{I}}\end{array}\right)\right],{\lambda }_{\mathrm{W}\mathrm{E}\mathrm{T}}\ge 0\\ 1-\mathrm{P}\mathrm{C}1\left[f\left(\begin{array}{l}\mathrm{N}\mathrm{D}\mathrm{V}{\mathrm{I}}_{\mathrm{N}\mathrm{I}},\mathrm{W}\mathrm{E}{\mathrm{T}}_{\mathrm{N}\mathrm{I}},\mathrm{L}\mathrm{S}{\mathrm{T}}_{\mathrm{N}\mathrm{I}},\\ \mathrm{M}\mathrm{S}{\mathrm{I}}_{\mathrm{N}\mathrm{I}},\mathrm{D}\mathrm{M}{\mathrm{I}}_{\mathrm{N}\mathrm{I}},\mathrm{D}{\mathrm{I}}_{\mathrm{N}\mathrm{I}}\end{array}\right)\right],{\lambda }_{\mathrm{W}\mathrm{E}\mathrm{T}}<0\end{array}\right.\end{array}$
式中:$\mathrm{A}\mathrm{R}\mathrm{S}\mathrm{E}{\mathrm{I}}_{0}$为归一化前的改进型干旱遥感生态指数;$\mathrm{N}\mathrm{D}\mathrm{V}{\mathrm{I}}_{\mathrm{N}\mathrm{I}}$$\mathrm{W}\mathrm{E}{\mathrm{T}}_{\mathrm{N}\mathrm{I}}$$\mathrm{L}\mathrm{S}{\mathrm{T}}_{\mathrm{N}\mathrm{I}}$$\mathrm{M}\mathrm{S}{\mathrm{I}}_{\mathrm{N}\mathrm{I}}$$\mathrm{D}\mathrm{M}{\mathrm{I}}_{\mathrm{N}\mathrm{I}}$$\mathrm{D}{\mathrm{I}}_{\mathrm{N}\mathrm{I}}$分别为归一化后的NDVI、WET、LST、MSI、DMI、DI;${\lambda }_{\mathrm{W}\mathrm{E}\mathrm{T}}$为WET对PC1的特征向量方向。
最后,将$\mathrm{A}\mathrm{R}\mathrm{S}\mathrm{E}{\mathrm{I}}_{0}$进行归一化,得到ARSEI,并将其按照[0.0,0.2)、[0.2,0.4)、[0.4,0.6)、[0.6,0.8)、[0.8,1.0]对其进行划分为差、较差、一般、良、优。

1.3.3 Sen+Mann-Kendall趋势分析

本文采用Sen方法[23](Sen’s slope estimator)用于ARSEI趋势分析。
${Q}_{\mathrm{A}\mathrm{R}\mathrm{S}\mathrm{E}\mathrm{I}}=\mathrm{M}\mathrm{e}\mathrm{d}\mathrm{i}\mathrm{a}\mathrm{n}\left(\frac{\mathrm{A}\mathrm{R}\mathrm{S}\mathrm{E}{\mathrm{I}}_{j}-\mathrm{A}\mathrm{R}\mathrm{S}\mathrm{E}{\mathrm{I}}_{i}}{j-i}\right),{\forall }_{j}>i$
式中:Median()为取中值;ij为时间序列的序号;${Q}_{\mathrm{A}\mathrm{R}\mathrm{S}\mathrm{E}\mathrm{I}}$为ARSEI的变化趋势,当${Q}_{\mathrm{A}\mathrm{R}\mathrm{S}\mathrm{E}\mathrm{I}}$大于0,则表示ARSEI为增长趋势;反之,则表示ARSEI为下降趋势。
由于Sen方法无法对数据的显著性进行检验,故使用Mann-Kendall趋势检验[24]对Sen方法进行显著性检。其计算公式如下:
$\mathrm{s}\mathrm{i}\mathrm{g}\mathrm{n}\left({x}_{j}-{x}_{i}\right)=\left\{\begin{array}{l}+1,\mathrm{ }\mathrm{ }\mathrm{A}\mathrm{R}\mathrm{S}\mathrm{E}{\mathrm{I}}_{j}-\mathrm{A}\mathrm{R}\mathrm{S}\mathrm{E}{\mathrm{I}}_{i}>0\\ 0,\mathrm{A}\mathrm{R}\mathrm{S}\mathrm{E}{\mathrm{I}}_{j}-\mathrm{A}\mathrm{R}\mathrm{S}\mathrm{E}{\mathrm{I}}_{i}=0\\ -1,\mathrm{ }\mathrm{ }\mathrm{A}\mathrm{R}\mathrm{S}\mathrm{E}{\mathrm{I}}_{j}-\mathrm{A}\mathrm{R}\mathrm{S}\mathrm{E}{\mathrm{I}}_{i}<0\end{array}\right.$
$S=\stackrel{n-1}{\displaystyle\sum _{i=1}} \stackrel{n}{\displaystyle\sum _{j=i+1}} \mathrm{s}\mathrm{i}\mathrm{g}\mathrm{n}\left(\mathrm{A}\mathrm{R}\mathrm{S}\mathrm{E}{\mathrm{I}}_{j}-\mathrm{A}\mathrm{R}\mathrm{S}\mathrm{E}{\mathrm{I}}_{i}\right)$
$Z=\left\{\begin{array}{cc}\frac{S}{\sqrt{\mathrm{V}\mathrm{a}\mathrm{r}\left(S\right)}},& S>0\\ 0,& \mathrm{ }S=0\\ \frac{S+1}{\sqrt{\mathrm{V}\mathrm{a}\mathrm{r}\left(S\right)}},& S<0\end{array}\right.$
$\mathrm{V}\mathrm{a}\mathrm{r}\left(S\right)=\frac{n\left(n-1\right)\left(2n+5\right)-\stackrel{m}{\displaystyle\sum _{k=1}}{t}_{k}({t}_{k}-1)(2{t}_{k}+5)}{18}$
式中:sign()为2个年份ARSEI值的变化方向;xixj分别为第i和第j个时间点对应的值;S为所有年份对之间变化方向的累计和;n为时间序列的样本量;Z为标准化后的检验统计量;Var为方差,可使用$\mathrm{V}\mathrm{a}\mathrm{r}\left(S\right)$式估算;tk为第k组中相同值的数量。
本文采用显著性水平α=0.05进行显著性检验,若|Z|大于临界值1.96则表示存在显著趋势,其判断方法如表3所示。
表3 Sen+Mann-Kendall 趋势等级划分

Tab. 3 Sen+Mann-Kendall trend classification criteria

${Q}_{\mathrm{A}\mathrm{R}\mathrm{S}\mathrm{E}\mathrm{I}}$ Z 趋势特征
${Q}_{\mathrm{A}\mathrm{R}\mathrm{S}\mathrm{E}\mathrm{I}}>0$ 1.96<Z 显著上升
Z<1.96 不显著上升
${Q}_{\mathrm{A}\mathrm{R}\mathrm{S}\mathrm{E}\mathrm{I}}=0$ Z 稳定不变
${Q}_{\mathrm{A}\mathrm{R}\mathrm{S}\mathrm{E}\mathrm{I}}<0$ Z<1.96 不显著下降
1.96<Z 显著下降

注:${Q}_{\mathrm{A}\mathrm{R}\mathrm{S}\mathrm{E}\mathrm{I}}$为ARSEI的变化趋势;Z为ARSEI变化的显著性水平。

1.3.4 基于ARSEI的生态脆弱性

生态系统的脆弱性指生态系统在一定时空间尺度下受到外界干扰后所表现出的敏感反应和自我恢复能力,开展生态脆弱性评估对于理解生态系统的结构与功能、制定生态环境整治政策具有重要意义。本文参考生态脆弱性定义[25]与已有研究[26],以能够综合反映生态系统健康状态并对外界环境变化做出响应的ARSEI为指标,选取ARSEI的年际波动情况来表示生态系统敏感性,ARSEI年际变率线性拟合趋势线的斜率来代表生态系统的适应性,并将其标准化后进行脆弱性的计算。最后,使用自然断点法将其分为极度、重度、中度、轻度脆弱与不脆弱,计算公式为:
敏感性[27]
$S=\frac{\stackrel{n}{\displaystyle\sum _{i=1}}\left|{F}_{i}-\stackrel{-}{F}\right|}{\stackrel{-}{F}}$
适应性[28]
y=Ax+B
$A=\frac{n\stackrel{n}{\displaystyle\sum _{i=1}}{x}_{i}{y}_{i}-\left(\stackrel{n}{\displaystyle\sum _{i=1}}{x}_{i}\right)\left(\stackrel{n}{\displaystyle\sum _{i=1}}{y}_{i}\right)}{n\stackrel{n}{\displaystyle\sum _{i=1}}{x}_{i}^{2}-{\left(\stackrel{n}{\displaystyle\sum _{i=1}}{x}_{i}\right)}^{2}}$
标准化与脆弱性[29]
$\mathrm{S}\mathrm{I}=\frac{{Z}_{i}-\mu }{\sigma }$
V=S'-A'
式中:i为第i年(n=23);Fi为第i年参数;$\stackrel{-}{F}$为参数平均值;S为生态系统敏感性;y为参数年际变率,是由每年的参数减去2000年到2023年参数的平均值;x为系统时间序列,为自然数1、2、3等,对应2000年到2023年;A为生态系统适应性;B为线性拟合后的常数项;Zi为原始栅格值;$\mu $$\sigma $分别为原始栅格的平均值和标准差;SI为标准化后数值;S'A'分别为对SA进行标准化;V为生态脆弱性。

1.3.5 参数最优地理探测器

地理探测器是由王劲峰等[30]提出的空间统计模型,广泛用于进行驱动力分析和因子分析。由于地理探测器的数据离散化过程易受经验和离散方法影响[31],因此使用多种离散方法,将每个参数分为3~7类,并对参数进行排序,从而实现自变量数据的最优离散化方法选取与执行,进而得到最优的解释力以及交互作用强度。其中交互作用类型如表4所示,地理探测器原理如下:
$q=1-\frac{\stackrel{L}{\displaystyle\sum _{h}}{N}_{h}{\sigma }_{h}^{2}}{{\sigma }^{2}{N}^{2}}$
式中:q为影响因子对生态环境质量的解释力;L为指标的分类数或分区;NhN分别为层数h和全局的单元数;${\sigma }_{h}^{2}$${\sigma }^{2}$分别为层数h和全局的方差。
表4 地理探测器交互作用探测类型

Tab. 4 Types of Geodetector interactions

判断依据 交互作用类型
$q\left({X}_{1}\bigcap {X}_{2}\right)<\mathrm{m}\mathrm{i}\mathrm{n}\left[q\left({X}_{1}\right),q\left({X}_{2}\right)\right]$ 非线性减弱
$\mathrm{m}\mathrm{i}\mathrm{n}\left[q\left({X}_{1}\right),q\left({X}_{2}\right)\right]<q\left({X}_{1}\bigcap {X}_{2}\right)<\mathrm{m}\mathrm{a}\mathrm{x}\left[q\left({X}_{1}\right),q\left({X}_{2}\right)\right]$ 单因子非线性减弱
$q\left({X}_{1}\bigcap {X}_{2}\right)>\mathrm{m}\mathrm{a}\mathrm{x}\left[q\left({X}_{1}\right),q\left({X}_{2}\right)\right]$ 双因子增强
$q\left({X}_{1}\bigcap {X}_{2}\right)=q\left({X}_{1}\right)+q\left({X}_{2}\right)$ 独立
$q\left({X}_{1}\bigcap {X}_{2}\right)>q\left({X}_{1}\right)+q\left({X}_{2}\right)$ 非线性增强

注:q为影响因子对生态环境质量的解释力;X1X2分别为2个进行交互作用探测的影响因子。

2 结果与分析

2.1 主成分分析结果

主成分分析PC1的有效荷载能够反映各指标分量对生态环境质量的影响效果,通过分析第一主成分的有效荷载,可以直接评估出ARSEI构建的合理性[32]。由表5可得,ARSEI的PC1贡献率虽因数据波动出现变化,但仍均高于80%,可见PC1可以集中6个指标的大部分特征。并在PC1中各个成分载荷分布稳定,且大小相对稳定,未出现载荷值剧烈波动而引发的特征信息丢失问题,表明ARSEI模型在时序分析中具有稳健性,能够可靠反映准东地区生态环境质量的长期变化。
表5 主成分分析结果

Tab. 5 Principal component analysis of ARSEI

年份 成分 各指标载荷值 特征值 贡献率/%
NDVI WET LST MSI DMI DI
2000 PC1 -0.574 -0.127 0.526 0.262 0.476 0.286 0.0951 87.84
PC2 -0.374 0.061 -0.824 0.224 0.163 0.316 0.0078 7.20
2005 PC1 -0.592 -0.135 0.473 0.301 0.462 0.330 0.0992 85.18
PC2 -0.304 -0.054 -0.847 0.245 0.093 0.344 0.0106 9.10
2010 PC1 -0.623 -0.164 0.475 0.316 0.379 0.339 0.1110 85.37
PC2 -0.336 0.115 -0.821 0.252 -0.014 0.370 0.0102 7.88
2015 PC1 -0.554 -0.098 0.485 0.316 0.502 0.309 0.0871 84.64
PC2 -0.335 -0.026 -0.833 0.301 0.044 0.317 0.0098 9.50
2018 PC1 -0.588 -0.104 0.526 0.373 0.261 0.399 0.0743 81.83
PC2 -0.322 0.092 -0.789 0.321 -0.136 0.378 0.0112 12.38
2023 PC1 -0.575 -0.159 0.493 0.355 0.374 0.367 0.0891 80.74
PC2 -0.333 0.070 -0.694 0.259 -0.299 0.496 0.0120 10.79

注:PC1、PC2分别为主成分分析后得到的第一、第二主成分。

通过图2表5进一步分析发现,在构建PC1的过程中,NDVI与WET的特征向量值为负值,且与ARSEI呈显著正相关,表明其对生态环境质量具有正向作用;而LST、MSI、DMI与DI在构建PC1过程中的特征向量值为正值,且与ARSEI呈显著负相关,表明其对生态环境质量具有负向影响。该结果与实际情况相符,符合生态学原理。此外,由PC1构建的ARSEI与各分项指数间的相关性系数均高于0.55,进一步验证了其合理性。相比之下,其他主成分中参数的特征向量值难以从生态学角度进行合理解释。综上所述,利用PC1构建的ARSEI表征区域生态环境质量是科学合理的。
图2 相关性分析结果

注:ARSEI为改进型干旱遥感生态指数;NDVI为归一化植被指数;WET为湿度指数;LST为地表温度指数;MSI为土壤盐渍化指数;DMI为土壤沙化指数;DI为PM2.5差值指数。下同。

Fig. 2 Correlation analysis results

2.2 ARSEI提取信息能力分析

为验证ARSEI在准东地区的适用性,研究选取工矿用地、裸地与耕地3个典型区域,进行RSEI与ARSEI细节提取的对比。从图3的信息中可得,ARSEI与RSEI在总体趋势上保持一致,但是ARSEI在可以反映出更多细节,可以更好的区分地物且更加符合事实。例如,在影像A区中,ARSEI相较于RSEI显著增强了露天矿区的识别能力,RSEI没有将矿区与附近的地物较好的区分;在影像B区中,ARSEI可以较好的识别有稀疏的植被覆盖与无植被覆盖的区域,进而更好的显示了盐碱地与沙地的轮廓特征与分布形态;在影像C区中,ARSEI对农田与居民区的反馈相较于RSEI在整体趋势相同的前提下,还突出了图像纹理细节。综上,本文对ARSEI修改是合理且具有实际价值,且ARSEI的评价结构相比较于RSEI更加贴近真实地表情况。
图3 2023年准东地区局部细节对比

注:RSEI为遥感生态指数。

Fig. 3 Comparison of local-scale details in the Zhundong region in 2023

2.3 准东地区生态环境分布特征以及变化趋势

研究区内生态环境质量差异显著,大致可分为3个区域:北部戈壁工矿密集区、中部生态过渡带和南部农业生态区。北部区域以戈壁与工矿用地为主,生态环境质量为“差”和“较差”,且由于矿产开采和工业建设导致该区域沙漠化、土壤侵蚀[33]以及土壤重金属累积[34]问题愈发严重,对本就脆弱的生态环境形成较大胁迫;中部区域作为荒漠-草地交错带,生态环境质量为“较差”和“一般”,起着一定的生态缓冲作用,同时也是对气候变化与人为干扰响应最敏感的区域;南部区域为草地、耕地和城镇所在区域,得益于天山冰川融水,该地区生态环境质量为“一般”以上,但该地区也面临过度开发和生态退化的威胁。
图4图5可得,2000—2023年准东地区不同等级生态环境质量的占比以及ARSEI均值发生了较为显著的变化。ARSEI均值在0.225到0.368之间波动,总体处于较低水平,其中,2000年为最高值0.368,2018年为最低值0.225。整体趋势呈现先稳定后下降最后上升的态势:2000—2010年ARSEI虽有下降但仍较为稳定;2010—2018年ARSEI持续下降;2018—2023年ARSEI缓慢恢复。具体来看,2000—2018年生态环境质量“差”的区域占比从17.7%上升至67.3%,主要由北部戈壁工矿密集区与中部生态过渡带退化而来;而生态环境质量“一般”以上的区域则从29.3%下降至18.8%,主要与南部农业生态区生态环境质量退化有关。2018—2023年经生态修复,生态环境质量“差”的区域占比从67.3%降至53.38%,这与中部生态过渡带区域生态环境质量改善有关;生态环境质量“一般”以上的区域从18.8%升至23.67%,主要来自南部农业生态区的提升。尽管如此,2023年ARSEI均值仍低于2000年,优质生态区域占比偏低。表明2018年以来的生态修复虽遏制恶化的趋势但恢复幅度有限,区域生态环境质量仍具较大提升空间[35]
图4 2000—2023年改进型干旱遥感生态指数空间分布

Fig. 4 Spatial distributions of arid modified remote sensing ecological index from 2000 to 2023

图5 生态环境质量等级构成

Fig. 5 Composition of ecological environment quality grades

由Sen+Mann-Kendall趋势分析与生态脆弱性分析可得图6,其与现有研究相符合[36-37]。2000—2023年准东地区间生态环境变化趋势以“显著下降”为主,占总面积83.5%,“不显著下降”“不显著上升”和“显著上升”区域分别占11.9%、3.12%和1.48%。具体而言,上升的区域,主要位于南部农业生态区的农田与城镇附近区域,而在其他区域均出现不同程度的下降。由生态脆弱性分析可得,“极度脆弱”与“重度脆弱”区域主要分布在北部戈壁工矿密集区,占总面积11.8%与36.3%;“中度”脆弱区域主要分布于中部生态过渡带,占总面积的22.7%;“轻度脆弱”与“不脆弱”区域主要分布于南部农业生态区,占总面积23.1%与6.12%。结合Sen+Mann-Kendall趋势分析与生态脆弱性分析可得,在生态环境质量中呈现下降的区域与生态“轻度脆弱”以上的区域高度重合,仅有部分生态不脆弱的地区生态环境质量上升。综上所述,准东地区的人类活动对原本脆弱生态环境产生了一定的负面影响;尽管存在生态保护与修复措施,但目前这些措施仅对人类活动密集且生态较稳定区域产生影响,对人类活动稀少且生态脆弱的区域积极影响尚未完全显现[38-39]。总体而言,2020—2023年准东地区生态环境仍以下降趋势为主,局部地区改善但未能扭转整体退化格局。
图6 生态状况变化趋势与生态脆弱性分布

Fig. 6 Trend of ecological changes and distribution of ecological fragility

2.4 准东地区最优地理探测器探测结果

本文采用参数最优的地理探测器,探究各个因子对准东地区生态环境质量变化的影响,共选取9个参数(X1~X9)涵盖人为与自然因素,其中不同土地利用类型下人类活动影响强度参数如表6所示。
表6 不同土地利用类型下人类活动影响强度参数

Tab. 6 Parameters of human activity impact intensity under different land use types

土地利用类型 人类活动影响强度参数
耕地 0.200
林地 0.082
草地 0.067
水体 0.143
建设用地 1.000
未利用土地 0.000
经自然断点法等5种分类方法分类选取最优方法后,多数参数最优分类数目为7种,少数为6种。图7结果显示,人类活动强度[40]、蒸散发量与气温为主要生态驱动力,对生态环境起着关键性作用。在自然因素中,蒸散发量与气温起关键作用,其平均解释力(q值)均高于0.35;在人为因素中,人类活动强度起关键作用,其q值均高于0.6。交互探测结果显示,对生态环境质量影响最大的为气温与人类活动强度的交互,其每年的交互作用q值均高于0.7,进一步凸显两者对准东地区生态环境变化的关键作用。此外,蒸散发量、降水、人口密度与其他因素结合对生态环境的影响显著高于其他因素。这些因素交互作用类型为非线性增强与双因子增强,它们会通过相互作用和协同影响形成特定的影响机制,进而调节其他驱动因子,最终对生态环境产生影响。这种复杂的相互作用和协同影响,揭示了准东地区生态环境变化的多因子驱动特性,强调了在生态保护和管理中综合考虑多种因素及其交互作用的重要性[41]
图7 参数最优地理探测器交互探测结果

注:X1为蒸散发量;X2为坡向;X3为坡度;X4为气温;X5为高程;X6为降水;X7为人口密度;X8为人类活动强度;X9为夜间灯光。

Fig. 7 Interactive detection results of optimal parameters-based geographical detector

2.5 准东地区生态现状成因分析

2.5.1 南优北劣格局的成因

准东地区南优北劣的生态格局并非单一因素所致,而是其自然底色与人类的活动强度共同作用的结果。准东地区位于欧亚大陆腹地,气候干旱少雨,北部深居戈壁,使得植被以稀疏耐旱灌木为主,导致生态脆弱;而南部地区则凭天山冰川融水,水分条件相对优越,形成较为良好的生态环境,这奠定南优北劣的自然底色[42]。在此基础上,人类活动进一步放大差异,北部工业基地年产大量原煤且煤化工产业集聚,导致植被退化、沙漠化、土壤侵蚀等问题不断加剧;而南部以农田和城镇为主,人类活动强度较低,生态扰动有限,退化压力显著较小。生态质量与产业强度错位,共同塑造了南优北劣的鲜明分异。

2.5.2 生态环境质量时空演变特征

自2005年自治区党委决定加快准东煤炭资源开发以来就拉开了准东大开发、大建设序幕,但在初期实际开发强度较低,直到2012年准东地区原煤产量突破3.1372×107 t。在2012年准东经济技术开发区成立后,准东地区经济发展迎来新的阶段,煤炭产能急剧扩张,2019年准东地区原煤产量达1.022×108 t。同时,由于人口持续增长和人类活动强度加大,南部草地大规模垦殖为农田,北部煤矿开采规模不断扩大[43],不节制的开发导致准东地区植被减少、水资源缺乏、土壤污染等生态问题开始激化。所以,2010年原国家环境保护部出具了文件对准东地区煤炭资源开发提出一定环保要求;2017年,将《准东环境保护法》列入地方性法规立法计划,并于2019年施行[44]。同时,准东地区也通过油田退出、采坑修复、荒漠绿化、水资源循环、绿色矿山和生物多样性保护6大核心工程,形成了生产退出、地貌恢复、生态重建与物种保护的全链条修复模式。结合本文内容,可知准东地区生态变化的全貌。2000—2010年准东地区进行了一定程度的开发,但因规模较小没有对生态环境质量造成较大的影响;2010—2018年准东地区开发快速,但由于环保措施滞后,生态环境质量持续下降;2018—2023年环境保护法落地,准东地区强化空间管控,进行一系列生态环境治理行动,准东地区生态环境质量有所回升。但修复效果呈现南快北慢的空间异质性,南部农田与城镇地区生态修复效果显著,生态恢复迅速,而北部戈壁区因水资源限制和土壤基质破坏[45],生态环境治理效果短时间内无法深入,生态恢复缓慢[46]。总体而言,区域已遏制生态环境恶化趋势,但整体质量仍显著低于2000年水平,生态恢复任务依然艰巨。

2.6 准东地区可持续发展的建议

本研究通过ARSEI与地理探测器揭示了准东地区生态环境演化规律以及驱动机制,为区域可持续发展提供3点对策建议。首先,分区施策,实施差异化生态修复。北部工矿区应加强绿色矿山建设,修复历史遗留采坑,遏制地表水土流失;中部生态过渡带需构建以耐旱灌木为主体的缓冲带,增强地表稳定性;南部农业区应发挥水资源优势建设生态环境,同时也应警惕过度开发导致退化。其次,应强化政策引导与法规执行,准东地区生态稳定离不开合理的政策引导。自《准东环境保护法》实行以来,生态质量总体向好,但形势依然严峻,应持续强化政策引导,健全生态质量动态监测与评估制度,建立以南部生态优势区为重点的生态补偿机制,激励区域生态保护行为。最后,围绕气温、蒸散发量与人类活动强度等主导因子,构建点面结合的调控路径,通过关键节点治理带动区域整体生态恢复,实现生态系统的精准保护与可持续利用。

3 结论

(1) 2000—2023年准东地区生态环境质量呈现先下降后上升的趋势,ARSEI均值变化范围为0.225到0.368,2018年为关键拐点达到最低点0.225,此后ARSEI均值呈现显著回升趋势,2023年回升至0.289;空间分布上,生态环境质量等级以“差”与“较差”为主,长期占比超过60%,且呈现明显的南优北劣格局,南部地区质量等级普遍优于北部。
(2) 生态环境退化范围广泛,2000—2023年准东地区83.5%的区域生态环境质量“显著下降”趋势,生态“轻度脆弱”以上的大多数区域生态环境质量均存在下降。生态环境质量上升的区域仅占4.6%,且集中分布于南部生态“不脆弱”的区域。这表明大规模资源开发活动对当地本就脆弱的生态环境产生了负面影响,而现有的生态保护措施的效果,目前主要体现于人类活动密集且生态脆弱性较低的南部周围区域。
(3) 单因子探测中人类活动强度、蒸散发量、气温是影响准东地区生态环境变化的关键因子,其q值平均大于0.35;因子交互作用后,q值普遍增强,以气温∩人类活动强度对生态环境影响最为显著,其q值达到0.7以上,是驱动该区域生态环境质量时空演变的关键因素。其他因子交互后也表现出增强效应,凸显了生态环境变化的多因子协同驱动特性。
本研究改进了ARSEI并评估了准东生态质量,但存在不足。首先,以ARSEI替代常用植被生产力作为脆弱性核心指标,却缺少可行性论证;其次,未对既有修复措施做长期动态评估与预测,建议的长期有效性待验。后续将聚焦这两方面深化。
[1]
Wu G, Tian J, Feng X, et al. Spatiotemporal variation and driving factors of ecological vulnerability in arid and semi-arid regions: A case study of Ningxia, China[J]. Catena, 2025, 259: 109378, doi: 10.1016/j.catena.2025.109378.

[2]
屈楠, 卢杰, 王子明, 等. 毛乌素沙地生态环境质量时空演变及驱动因素[J]. 水土保持通报, 2025, 45(5): 336-346.

[Qu Nan, Lu Jie, Wang Ziming, et al. Spatio-temporal evolution and driving factors of ecological environment quality in Mu Us sandy land[J]. Bulletin of Soil and Water Conservation, 2025, 45(5): 336-346.]

[3]
徐勇, 戴强玉, 黄雯婷, 等. 2000—2020年西南地区植被NDVI时空变化及驱动机制探究[J]. 环境科学, 2023, 44(1): 323-335.

[Xu Yong, Dai Qiangyu, Huang Wenting, et al. Spatio-temporal variation in vegetation cover and its driving mechanism exploration in southwest China from 2000 to 2020[J]. Environmental Science, 2023, 44(1): 323-335.]

DOI

[4]
刘博, 韩庆功, 张杰琳, 等. 变化环境下黄河中游土壤湿度时空格局及其影响因素[J]. 水土保持学报, 2025, 39(1): 73-83.

[Liu Bo, Han Qinggong, Zhang Jielin, et al. Spatial-temporal patterns and factors of soil moisture in the middle reaches of the Yellow River under changing environments[J]. Journal of Soil and Water Conservation, 2025, 39(1): 73-83.]

DOI

[5]
国家环保总局. 生态环境状况评价技术规范(试行)[EB/OL]. [2006-03-09]. https://www.mee.gov.cn/gkml/zj/gg/200910/t20091021_171618.htm.

[State Environmental Protection Agency. Te- chnical criterion for eco-environment status evaluation[EB/OL]. [2006-03-09]. https://www.mee.gov.cn/gkml/zj/gg/200910/t20091021_171618.htm.]

[6]
徐涵秋. 城市遥感生态指数的创建及其应用[J]. 生态学报, 2013, 33(24): 7853-7862.

[Xu Hanqiu. A remote sensing urban ecological index and its application[J]. Acta Ecologica Sinica, 2013, 33(24): 7853-7862.]

[7]
Liu Y, Xu W H, Hong Z H, et al. Integrating three-dimensional greenness into RSEI improved the scientificity of ecological environment quality assessment for forest[J]. Ecological Indicators, 2023, 156: 111092, doi: 10.1016/j.ecolind.2023.111092.

[8]
么廷禹, 王东升, 邰姗姗. 基于遥感生态指数的鞍山市生态环境质量变化研究[J]. 环境科学与技术, 2024, 47(12): 199-210.

[Yao Tingyu, Wang Dongsheng, Tai Shanshan. Evaluation of ecological environment quality and analysis of influencing factors in Anshan City based on remote sensing ecological index[J]. Environmental Science & Technology, 2024, 47(12): 199-210.]

[9]
王杰, 马佳丽, 解斐斐, 等. 干旱地区遥感生态指数的改进——以乌兰布和沙漠为例[J]. 应用生态学报, 2020, 31(11): 3795-3804.

DOI

[Wang Jie, Ma Jiali, Xie Feifei, et al. Improvement of remote sensing ecological index in arid regions: Taking Ulan Buh Desert as an example[J]. Chinese Journal of Applied Ecology, 2020, 31(11): 3795-3804.]

DOI

[10]
马鹏飞, 张安兵, 王贺封, 等. 基于改进遥感生态指数的矿区生态环境监测及驱动因素分析[J]. 生态与农村环境学报, 2023, 39(12): 1580-1590.

[Ma Pengfei, Zhang Anbing, Wang Hefeng, et al. Ecological environment monitoring and driving factors analysis of mining area based on improved remote sensing ecological index[J]. Journal of Ecology and Rural Environment, 2023, 39(12): 1580-1590.]

[11]
汤从沧, 李巧, 陶洪飞, 等. 基于改进遥感生态指数模型的塔里木河流域生态环境质量评价[J]. 环境科学, 2025, 46(7): 4485-4498.

[Tang Congcang, Li Qiao, Tao Hongfei, et al. Evaluating the ecological environment quality in the Tarim River Basin based on the modified remote sensing ecological index model[J]. Environmental Science, 2025, 46(7): 4485-4498.]

[12]
Zhang Y Q, Fang J. Developing a remote sensing-based ecological index based on improved biophysical features[J]. Journal of Applied Remote Sensing, 2022, 16(1): 12008, doi: 10.117/1.JRS.16.012008.

[13]
Yu Z H, Jin J, Hou F X, et al. Understanding effect of phosphorus-based additive on ash deposition characteristics during high-sodium and high-calcium Zhundong coal combustion in drop-tube furnace[J]. Fuel, 2021, 287: 119462, doi: 10.1016/j.fuel.2020.119462.

[14]
曾强, 李根生, 董敬宣, 等. 新疆煤炭资源开采典型生态环境问题及对策[J]. 矿业安全与环保, 2017, 44(1): 106-110.

[Zeng Qiang, Li Gensheng, Dong Jingxuan, et al. Typical ecological and environmental issues and countermeasures in coal mining in Xinjiang region[J]. Mining Safety and Environmental Protection, 2017, 44(1): 106-110.]

[15]
张东升, 刘洪林, 范钢伟, 等. 新疆大型煤炭基地科学采矿的内涵与展望[J]. 采矿与安全工程学报, 2015, 32(1): 1-6.

[Zhang Dongsheng, Liu Honglin, Fan Gangwei, et al. Connotation and prospection on scientific mining of large Xinjiang coal base[J]. Journal of Mining and Safety Engineering, 2015, 32(1): 1-6.]

[16]
薛一波, 张小啸, 雷加强, 等. 新疆大气颗粒物污染时空演变及沙尘组分研究[J]. 中国环境科学, 2024, 44(6): 3012-3020.

[Xue Yibo, Zhang Xiaoxiao, Lei Jiaqiang, et al. Study on the spatio-temporal variation of atmospheric particulate pollution and its dust components in Xinjiang, northwest China[J]. China Environmental Science, 2024, 44(6): 3012-3020.]

[17]
赵旭, 李浩, 朱益虎, 等. 融合高分影像与时序NDVI的农作物语义分割模型[J]. 农业工程学报, 2025, 41(14): 216-227.

[Zhao Xu, Li Hao, Zhu Yihu, et al. Crop semantic segmentation model integrating high-resolution imagery and time-series NDVI[J]. Transactions of the Chinese Society of Agricultural Engineering, 2025, 41(14): 216-227.]

[18]
Peng X, Li X, Shen Y Y B, et al. Applicability of different indices for delineating the boundary of arid region: A case study in northwestern China[J]. Journal of Earth Science, 2025, 36(1): 212-222.

DOI

[19]
Wen X L, Ming Y L, Gao Y G, et al. Dynamic monitoring and analysis of ecological quality of pingtan comprehensive experimental zone, a new type of sea island city, based on RSEI[J]. Sustainability, 2019, 12(1): 21, doi: 10.3390/su12010021.

[20]
李世娇, 张珂珂, 谢宝妮, 等. 基于遥感的库尔勒地区生态环境质量评价及成因分析[J]. 干旱区地理, 2024, 47(12): 2064-2074.

DOI

[Li Shijiao, Zhang Keke, Xie Baoni, et al. Ecological environment quality evaluation and driving factors of Korla region based on remote sensing[J]. Arid Land Geography, 2024, 47(12): 2064-2074.]

DOI

[21]
冯海英, 冯仲科, 冯海霞. 一种基于差值指数的颗粒物PM2.5浓度反演新方法[J]. 光谱学与光谱分析, 2018, 38(10): 3012-3016.

[Feng Haiying, Feng Zhongke, Feng Haixia. One new method of PM2.5 concentration inversion based on difference index[J]. Spectroscopy and Spectral Analysis, 2018, 38(10): 3012-3016.]

[22]
Li N, Wang J Y, Qin F. The improvement of ecological environment index model RSEI[J]. Arabian Journal of Geosciences, 2020, 13(11): 403, doi: 10.1007/s12517-020-05414-7.

[23]
Şen Z, Şişman E. Risk attachment Sen’s slope calculation in hydrometeorological trend analysis[J]. Natural Hazards, 2024, 120(4): 3239-3252.

DOI

[24]
王南, 刘泽轩, 郑江华, 等. 天山冰湖分布时空特征及驱动力分析[J]. 干旱区地理, 2024, 47(9): 1472-1481.

DOI

[Wang Nan, Liu Zexuan, Zheng Jianghua, et al. Spatiotemporal characteristics and driving forces of glacial lakes in Tianshan Mountains[J]. Arid Land Geography, 2024, 47(9): 1472-1481.]

DOI

[25]
IPCC Climate change 2014: Impacts, adaptation, and vulnerability. Part A:Global and sectoral aspects. Contribution of working group II to the fifth assessment report of the intergovernmental panel on climate change[M]. Cambridge, UK: Cambridge University Press, 2014.

[26]
王鹤松, 何敏, 闫薇, 等. 基于遥感总初级生产力的天山-塔里木绿洲地区生态系统脆弱性研究[J]. 生态学报, 2021, 41(24): 9729-9737.

[Wang Hesong, He Min, Yan Wei, et al. Ecosystem vulnerability in the Tianshan Mountains and Tarim Oasis based on remote sensed gross primary productivity[J]. Acta Ecologica Sinica, 2021, 41(24): 9729-9737.]

[27]
尤南山, 蒙吉军. 基于生态敏感性和生态系统服务的黑河中游生态功能区划与生态系统管理[J]. 中国沙漠, 2017, 37(1): 186-197.

DOI

[You Nanshan, Meng Jijun. Ecological functions regionalization and ecosystem management based on the ecological sensitivity and ecosystem service in the middle reaches of the Heihe River[J]. Journal of Desert Research, 2017, 37(1): 186-197.]

DOI

[28]
Sun J, Qin X J. Precipitation and temperature regulate the seasonal changes of NDVI across the Tibetan Plateau[J]. Environmental Earth Sciences, 2016, 75: 1-9.

DOI

[29]
於琍, 曹明奎, 陶波, 等. 基于潜在植被的中国陆地生态系统对气候变化的脆弱性定量评价[J]. 植物生态学报, 2008, 32(3): 521-530.

DOI

[Yu Li, Cao Mingkui, Tao Bo, et al. Quantitative assessment of the vulnerability of terrestrial ecosystems of China to climate change based on potential vegetation[J]. Chinese Journal of Plant Ecology, 2008, 32(3): 521-530.]

DOI

[30]
王劲峰, 徐成东. 地理探测器: 原理与展望[J]. 地理学报, 2017, 72(1): 116-134.

DOI

[Wang Jinfeng, Xu Chengdong. Geodetector: Principle and prospective[J]. Acta Geographica Sinica, 2017, 72(1): 116-134.]

DOI

[31]
魏伟, 尹力. 东北地区“三区空间”格局演化特征及驱动机制[J]. 地理科学, 2023, 43(2): 324-336.

DOI

[Wei Wei, Yin Li. Evolutionary characteristics and driving mechanisms of the “three-zone space”pattern in northeast China[J]. Scientia Geographica Sinica, 2023, 43(2): 324-336.]

DOI

[32]
冯东梅, 夏郡, 张步林. 多尺度视角下基于MRSEI的东北地区生态环境质量时空演变及分区管控[J/OL]. 环境科学. [2025-07-22]. https://doi.org/10.13227/j.hjkx.202504177.

[Feng Dongmei, Xia Jun, Zhang Bulin. Spatial-temporal evolution and zoning management of ecological environment quality in northeast China based on MRSEI from a multi-scale perspective[J/OL]. Environm-ental Science. [2025-07-22]. https://doi.org/10.13227/j.hjkx.202504177.]

[33]
卞正富, 于昊辰, 侯竟, 等. 西部重点煤矿区土地退化的影响因素及其评估[J]. 煤炭学报, 2020, 45(1): 338-350.

[Bian Zhengfu, Yu Haochen, Hou Jing, et al. Influencing factors and evaluation of land degradation of 12 coal mine areas in western China[J]. Journal of China Coal Society, 2020, 45(1): 338-350.]

[34]
袁婷婷, 王志强, 汪溪远, 等. 准东红沙泉矿区重金属生态风险缓冲区分析[J]. 土壤通报, 2020, 51(1): 227-233.

[Yuan Ting-ting, Wang Zhiqiang, Wang Xiyuan, et al. Buffer analysis of heavy metal ecological risk in the Hongshaquan mining area of east Junggar Basin[J]. Soil Bulletin, 2020, 51(1): 227-233.]

[35]
刘英, 许萍萍, 毕银丽, 等. 新疆戈壁煤矿露天开采对生态环境扰动定量分析[J]. 煤炭学报, 2023, 48(2): 959-974.

[Liu Ying, Xu Pingping, Bi Yinli, et al. Quantitative analysis of coal mining disturbance on environment in Xinjiang Gobi open-pit mining area[J]. Journal of China Coal Society, 2023, 48(2): 959-974.]

[36]
马丽莎, 郑江华, 彭建, 等. 新疆生态脆弱性特征及其驱动力[J]. 生态学报, 2024, 44(20): 9053-9066.

[Ma Lisha, Zheng Jianghua, Peng Jian, et al. Research on the characteristics and driving forces of ecological vulnerability in Xinjiang[J]. Acta Ecologica Sinica, 2024, 44(20): 9053-9066.]

[37]
于昊辰, 陈浮, 尹登玉, 等. 采矿与气候变化对戈壁矿区土地生态系统的影响——以准东煤炭基地为例[J]. 煤炭学报, 2021, 46(8): 2650-2663.

[Yu Haochen, Chen Fu, Yin Dengyu, et al. Effects of mining activities and climate change on land ecosystem in Gobi mining area: A case study of Zhundong coal base[J]. Journal of China Coal Society, 2021, 46(8): 2650-2663.]

[38]
Jiang J, Abulizi A, Abliz A, et al. Construction of landscape ecological security pattern in the Zhundong region, Xinjiang, NW China[J]. International Journal of Environmental Research and Public Health, 2022, 19(10): 6301, doi: 10.3390/ijerph19106301.

[39]
Shen L, Zeng Q. Multiscenario simulation of land use and land cover in the Zhundong mining area, Xinjiang, China[J]. Ecological Indicators, 2022, 145: 109608, doi: 10.1016/j.ecolind.2022.109608.

[40]
徐勇, 孙晓一, 汤青. 陆地表层人类活动强度: 概念、方法及应用[J]. 地理学报, 2015, 70(7): 1068-1079.

DOI

[Xu Yong, Sun Xiaoyi, Tang Qing. Human activity intensity of land surface: Concept, method and application in China[J]. Acta Geographica Sinica, 2015, 70(7): 1068-1079.]

DOI

[41]
李怡乐, 张乐艺, 李霞, 等. 基于改进遥感生态指数的粤港澳大湾区生态环境质量监测及驱动力分析[J/OL]. 环境科学. [2025-07-26]. https://doi.org/10.13227/j.hjkx.202504062.

[Li Yi-le, Zhang Leyi, Li Xia, et al. Ecological environment monitoring and driving factors analysis of the greater bay area based on improved remote sensing ecological index[J/OL]. Environmental Science. [2025-07-26]. https://doi.org/10.13227/j.hjkx.202504062.]

[42]
李梦怡, 邓铭江, 凌红波, 等. 新疆绿洲格局变化与生态耗水结构分析[J]. 水科学进展, 2023, 34(5): 719-730.

[Li Mengyi, Deng Mingjiang, Ling Hongbo, et al. Analysis on the change of oasis pattern and ecological water consumption structure in Xinjiang[J]. Advances in Water Science, 2023, 34(5): 719-730.]

[43]
刘芳, 仇韪, 夏楠, 等. 准东南煤炭资源地区土地利用转型及其对生态系统服务价值的影响[J]. 中国煤炭, 2024, 50(9): 144-153.

[Liu Fang, Qiu Wei, Xia Nan, et al. Land use transition of coal-based areas in Southeastern Junggar Basin and its impact on ecosystem service value[J]. China Coal, 2024, 50(9): 144-153.]

[44]
新疆维吾尔自治区发展和改革委员会. 新疆首个开发区生态环境保护条例施行[EB/OL]. [2025-07-29]. https://xjdrc.xinjiang.gov.cn/xjfgw/c108390/201911/ecef101606d946f7b4ce74d7ebc5155a.shtml.

[Development and Reform Commission of Xinjiang Uygur Autonomous Region. The first environmental protection regulation for development zones in Xinjiang is implemented[EB/OL]. [2025-07-29]. https://xjdrc.xinjiang.gov.cn/xjfgw/c108390/201911/ecef101606d946f7b4ce74d7ebc5155a.shtml.]

[45]
张俊杰, 白中科, 杨博宇. 新疆荒漠露天矿区生态受损及砾幕层重构方法研究[J]. 地学前缘, 2021, 28(4): 142-152.

DOI

[Zhang Junjie, Bai Zhongke, Yang Boyu. Gravel curtain layer in the desert open pit mining area of Xinjiang: Ecological damage and reconstruction method[J]. Earth Science Frontiers, 2021, 28(4): 142-152.]

DOI

[46]
Ou B, Abulizi A, Zayiti A, et al. Ecological risk zoning control in Zhundong economic development zone based on landscape pattern changes[J]. Sustainability, 2023, 15(22): 15972, doi: 10.3390/su152215972.

文章导航

/