Evaluation of eco-environmental quality in typical rocky desertification areas in the upper reaches of the Yangtze River

云南省曲靖市会泽县(25°48′~27°04′N，103°03′~103°55′E)位于云南省东北部、金沙江东岸、曲靖市西北部，是石漠化综合治理试点县和石漠化重点监测县。全县地势西高东低，由西向东阶梯状递减，平均海拔为2 220 m。立体气候分布明显，同时存在温带高原季风气候和南亚热带气候至寒温带气候，年平均气温为12.7 ℃，年降水量为800 mm。境内植被以亚高山灌丛、草甸、阔叶类混交林、亚热带稀树草原旱生植被为主。县境内碳酸盐岩分布广泛，多为灰岩和白云岩，为典型的喀斯特山区。

遥感数据为2002年2月25日Landsat-7 ETM+、2010年2月7日Landsat-5 TM、2018年3月1日Landsat-8 OLI_TIRS等卫星影像数据。数字高程模型数据来自于地理空间数据云(http://www.gscloud.cn)。云南省地质图来自于中国科学院地球化学研究所(http://www.gyig.cas.cn/)的喀斯特数据中心。对原始影像进行辐射校正、大气校正、几何校正等预处理，结合行政区划矢量数据拼接影像并裁剪出研究区。

基于植被覆盖度(V_F)、岩石裸露率(R_B)、坡度(s)等3个指标，通过综合分析法构建石漠化遥感监测指标体系^[16-18]。计算归一化植被指数、归一化岩石指数，并采用像元二分模型分别计算植被覆盖度和岩石裸露率。按照表1分别对植被覆盖度、岩石裸露率、坡度赋值，依据公式计算石漠化综合指标(I_R)：I_R=0.4V_F+0.4R_B+0.2s。采用决策树分类法划分石漠化等级：0≤I_R≤1为无石漠化、1＜I_R≤2为潜在石漠化、2＜I_R≤4为轻度石漠化、4＜I_R≤6为中度石漠化、6＜I_R≤8为强度石漠化、8＜I_R≤10为极强度石漠化^[18]。依据云南省地质图剔除非喀斯特区，得到石漠化等级分布图。

采用遥感生态指数(I_RSE)进行生态环境质量评价^[7]，该指数耦合了绿度、湿度、干度、热度4个指标，分别用归一化植被指数、湿度分量、建筑裸土指数、地表温度表示，通过主成分分析评价当地生态环境。

植被状况是生态环境质量评价的重要方面。归一化植被指数(I_NDV)作为应用最广的植被指数，能够很好地反映土地上的植被覆盖状况^[19]。计算公式为：I_NDV=(ρ_NIR−ρ_R)/(ρ_NIR+ρ_R)。其中：ρ_NIR和ρ_R分别为卫星影像数据的近红外和红外波段的反射率。

缨帽变换得到的第三分量(湿度分量)与植被湿度、土壤湿度紧密相关^[20]。利用湿度分量(W)代表湿度指标，计算公式为：W=C₁ρ_B+C₂ρ_G+C₃ρ_R+C₄ρ_NIR+C₅ρ_SWIR1+C₆ρ_SWIR2。其中：ρ_B、ρ_G、ρ_R、ρ_NIR、ρ_SWIR1、ρ_SWIR2分别为卫星影像数据的蓝波段、绿波段、红波段、近红外波段、短波红外1波段、短波红外2波段的反射率。C₁~C₆为湿度分量的计算系数，对于TM^[21]，C₁=0.031 5，C₂=0.202 1，C₃=0.310 2，C₄=0.159 4，C₅ =−0.680 6，C₆=−0.610 9；对于ETM+^[22]，C₁=0.262 6，C₂=0.214 1，C₃=0.092 6，C₄=0.065 6，C₅=−0.762 9，C₆=−0.538 8；对于OLI^[23]，C₁ =0.151 1，C₂=0.197 3，C₃=0.328 3，C₄=0.340 7，C₅ =−0.711 7，C₆ =−0.455 9。

研究区存在大量石漠化裸土地区，因此干度指标以裸土指数(I_S)与建筑指数(I_IB)的均值表示，记为建筑裸土指数(I_NDBS)^[24-25]。干度指标对生态环境质量起负面作用。计算公式为：

地表温度与植被的生长发育、城市热岛效应等密切相关，因此本研究采用地表温度表示热度指标。热度指标对生态环境质量起负面作用。地表温度的反演采用大气校正法^[26]。把卫星传感器接收的辐射亮度L_λ 分为大气向上辐射亮度L^↑、地面的真实辐射亮度经过大气层之后到达卫星传感器的能量、大气向下辐射到达地面后反射的能量L_↓，得到公式：

其中：ε为地表比辐射率；T_S卫星传感器处的辐射亮度值；B(T_S)为T_S下的黑体辐射亮度值；τ为大气在热红外波段的透过率。黑体辐射亮度B(T_S)的计算公式为：

地表比辐射率(ε)的计算公式根据植被覆盖度(F_v)分为水体、城镇、自然3部分^[27]，ε_water为水体地表比辐射率；ε_building为城镇地表比辐射率；ε_natural为自然地表比辐射率。

地表真实温度(T_LS)通过普朗克公式求得：

其中：K₁、K₂分别为定标参数。对于TM数据，K₁=607.76 W·m⁻²·sr⁻¹·μm⁻¹，K₂=1 260.56 K；对于ETM+数据，K₁=666.09 W·m⁻²·sr⁻¹·μm⁻¹、K₂ = 1 282.71 K；对于TIRS band10数据，K₁=774.89 W·m⁻²·sr⁻¹·μm⁻¹、K₂=1 321.08 K。

分别将绿度、湿度、干度、热度4个生态指标正规化，减少因数值大小不同带来的影响，使其值位于[0，1]。采用主成分变换分析集合4个指标的新影像，计算得到第一主成分载荷值(P_C1)及相关统计结果。为了便于结果的对比分析，用1−P_C1，获得初始遥感生态指数I_RSE0，并将其正规化得到遥感生态指数I_RSE。I_RSE值越大，生态环境质量越好，反之越差。即：

以0.2为间隔将I_RSE划分为5个区间^[7]，分别代表生态环境质量差(0~0.2)、较差(0.2~0.4)、中等(0.4~0.6)、良好(0.6~0.8)、优(0.8~1.0)。

使用ArcGIS对会泽县进行均匀采样(1 km×1 km)，去除非喀斯特地区后，共得到 2 366个采样点，获得对应点位的I_R与I_RSE。为使石漠化综合指标大的值代表石漠化程度较低地区，对I_R进行归一化，用1减去归一化后的I_R，得到转换后的石漠化综合指标(I_R0)，对I_R0与I_RSE进行相关性分析。

由表2可知：2002—2018年，研究区已石漠化地区(极强度、强度、中度、轻度石漠化)面积减少了583.33 km²，占全域面积的24.57%；无石漠化、潜在石漠化面积增加，轻度石漠化面积减少，中度、强度、极强度石漠化面积先增加后减少，研究区内石漠化程度整体得到改善。具体来看，2002—2010年，已石漠化地区面积总体减少了253.22 km²，主要为轻度石漠化面积减少(432.73 km²) ，年均降率为8.58%。2010—2018年，无石漠化面积显著增加，达284.48 km²，年均增率为10.13%；极强度石漠化面积减少了28.38 km²，年均降率为8.23%。

2002—2018年研究区I_RSE均值变化幅度不大，整体呈上升趋势，生态环境质量总体呈中等水平(表3)。16 a间生态等级为优和差的地区均有增加。与2012年相比，2018年研究区I_RSE略有所下降，其中绿度指标、湿度指标均为正值，提示此两者对生态环境起到了积极的作用；干度指标、热度指标均为负值，提示此两者对环境起到了消极作用。同时，2002、2010年绿度、湿度的P_C1之和大于干度、热度P_C1之和的绝对值，表明绿度、湿度对生态环境的改善作用强于干度、热度的负面作用。2018年绿度、湿度的P_C1之和小于干度、热度P_C1之和的绝对值，表明干度、热度的负面作用强于绿度、湿度的改善作用。16 a间，负面作用的加强主要来自于干度指标，其P_C1由−0.029变为−0.622。

对不同年份不同生态环境等级地区进行面积统计(表4)可知：16 a间研究区生态环境恶化与优化并存，中等、较差和良好等级地区占比之和均超过了90%，但区域等级异动较大，优和差等级比例不断提高。其中，中等生态等级面积占比始终最高，但逐年来呈减少趋势，16 a间降率达14.53%。较差及以下等级面积先减少后增加，但整体变化幅度不大。良好及优等级面积增加，增幅达12.92%。由表5可知：2018年，生态等级为差的区域主要分布在迤车镇、马路乡、火红乡、矿山镇、者海镇、大井镇、娜姑镇、大海乡，优等级主要分布在老厂乡、五星乡、金钟镇、新街回族乡、待补镇、鲁纳乡。

基于遥感生态指数的5个等级，对会泽县2002和2018年遥感生态指数进行差值变化检测，按照级差符号的正、0、负依次归为生态变好、不变、变差。由表6可知： 2002—2018年，研究区生态环境质量等级变好面积为1 605.37 km²，占全域面积的27.42%；生态变差地区面积为883.58 km²，占全域面积的15.09%，可见研究区生态环境总体呈变好趋势。由表7可知：生态环境质量等级变好地区主要分布在老厂乡、五星乡、金钟镇、新街回族乡、待补镇、纳鲁乡、上村乡、架车乡、田坝乡等地。这些乡镇海拔较高，人类活动相对较少，全县范围内大力实施生态环境保护，如天然林保护、公益林补植补造、石漠化治理等，使得生态环境得到了极大改善^[28]。生态环境质量等级变差地区主要分布在研究区东北部及西南部等地势较平坦的乡镇，包括纸厂乡、迤车镇、马路乡、火红乡、乐业镇、矿山镇、者海镇、大井镇、娜姑镇、大海乡。这些乡镇大力发展三大产业，国内生产总值猛增，城镇化发展迅速，建设用地激增，一定程度上影响了当地生态状况^[28]。

由表8可知：极强度石漠化地区，生态环境质量总体较差，强度、中度石漠化地区遥感生态指数等级以较差、中等为主，轻度、潜在石漠化地区等级主要表现为中等、良好，无石漠化地区，生态环境质量主要为良好，而非喀斯特地区，较差、中等、良好等级的占比均较大。可见石漠化与生态环境之间存在一定相关性。2002—2018年，无、潜在、轻度石漠化地区，生态环境均明显好转。其中，无石漠化和潜在石漠化地区，生态环境主要由中等、良好转为良好、优；轻度石漠化地区，生态环境主要由较差、中等、良好转为中等、良好。反之，强度、极强度石漠化地区生态环境质量逐步变差。其中，极强度石漠化地区，生态环境质量为差等级的面积增率达28.11%。这可能与干旱有关。16 a间研究区干度指标 的P_C1由−0.029变为−0.622，而湿度指标P_C1由0.608下降为0.361 (表3)，均对I_RSE产生负面作用，从而表现为以强度、极强度石漠化为代表的干旱地区生态环境质量恶化。因此，在生产建设、生态环境恢复中，应进一步注意石漠化严重地区的保护，减少人为破坏干扰。

分别计算2002、2010、2018年转换后I_R0与I_RSE的相关性系数可知：2002年两者相关性系数为0.688，2010年为0.750，2018年为0.873，3期均在0.01水平上显著相关。可见，生态环境质量与石漠化之间存在正相关，随着石漠化改善，生态环境质量逐渐向好，即通过RSEI法进行石漠化地区生态环境质量评价具有很好的效果。

本研究基于2002—2018年Landsat 5、 Landsat 7、Landsat 8卫星影像数据，划分会泽县石漠化等级，并利用遥感生态指数法定量评价会泽县生态环境质量。结果表明：①会泽县石漠化状况整体得到改善，已石漠化地区总面积减少了583.33 km²，其中，轻度石漠化地区面积减少432.73 km²。②干度指标的载荷值绝对值不断增大，P_C1载荷值由−0.029变为−0.622，对I_RSE的贡献度不断增加，主要为地势平坦地区的城镇化建设用地增加，成为制约会泽县生态环境质量优化的重要因素。③ 2002、2010、2018年遥感生态指数均值分别为0.458、0.490、0.488，生态环境质量整体表现为中等，生态环境呈优化趋势。④会泽县内生态环境质量优化面积占总面积的27.42%，主要分布在老厂乡、五星乡、金钟镇、新街回族乡、待补镇、纳鲁乡、上村乡、架车乡、田坝乡；生态环境质量变差地区面积占总面积的15.09%，分布在纸厂乡、迤车镇、马路乡、火红乡、乐业镇、矿山镇、者海镇、大井镇、娜姑镇、大海乡。⑤石漠化与生态环境质量呈显著正相关，相关性指数为0.688~0.873，表明通过遥感生态指数法评价长江上游典型石漠化地区的生态环境质量效果较好。