南水北调中线工程作为国家的重大基础工程，对于缓解我国北方地区的水资源紧张的问题具有重要的现实意义。丹江流域(河南段)作为其主要的水源地，其生态质量的好坏引起人们的高度关注，而近20年来随着城市化进程的加快，丹江流域(河南段)部分区域的生态环境明显遭到破坏，如何对流域的生态质量进行大范围、长时间周期性监测与评估已成为首要的考虑问题。遥感因其可以对地面进行实时快速的周期性观测的优点在生态环境领域得到广泛应用，已成为监测区域生态质量的重要手段^[1-4]。

对于遥感数据，通常采用计算各类遥感指数来监测区域的生态质量，例如：孟丹等^[5]利用NDVI数据并结合降水、气温等资料对京津冀等地区NDVI的变化特征和空间分布进行研究，并分析影响植被覆盖度变化的气候因子驱动力；罗敏等^[6]以NDVI为指示因子，分析塔里木河2000—2013年间NDVI的时空变化特征，探讨区分降水、气温气候因子以及人类活动对植被覆被变化的影响；丁小松等^[7]利用Landsat系列影像提取水体指数并计算水陆分割阈值，对大连市獐子岛的海岸线进行自动提取并分析其时空变化原因；王小平等^[8]通过提取各类水体指数值并结合实测数据建模分析，对艾比湖流域的地表水水质进行监测；王乐^[9]对西安市不同时期的地表温度进行反演，并结合其他研究数据，分析研究西安市的城市热岛效应。但是，单一的遥感指数并不能完全反映区域的生态质量状况，学者开始考虑综合多个遥感指数信息监测区域的生态质量。由徐涵秋提出的遥感生态指数(RSEI)目前已经得到了较为广泛的应用^[10-13]，例如刘智才等^[12]利用遥感生态指数对杭州市的生态变化进行监测，并分析引起杭州市生态质量降低的原因。该指数集成了由遥感影像中提取的某些能反映区域生态质量的指标，对其进行长时间周期性的监测，具有较高的参考价值；但是目前遥感生态指数主要以城市为研究对象，监测分析其生态质量变化，对于流域尺度的生态质量变化的研究较少，特别是对于丹江流域(河南段)的生态质量研究尚未见报道。

因此，为进一步研究南水北调中线工程水源地的生态质量变化情况，本文以丹江流域(河南段)为研究区，利用1998、2008和2017年的Landsat系列卫星影像，集成耦合得到研究区各年份的遥感生态指数，对研究区历年的生态质量状况进行监测，并结合研究区的经济发展状况分析其变化原因，以便为水源地的保护提供合理的科学依据和数据支撑。本研究为保护水源地的水质安全，以及监测区域的植被覆盖和水土流失状况提供一种新的技术手段。

丹江流域(河南段)位于河南省西南部，秦岭余脉伏牛山南麓，与湖北省、陕西省相临，包含三门峡市卢氏县部分乡镇、洛阳市栾川县部分乡镇，南阳市的西峡县、淅川县，以及内乡县、邓州市部分乡镇，总面积为8 218 km², 其中南阳市6 668 km²，三门峡市1 238 km²，洛阳市312 km²，地理坐标为E 110°52′~112°，N 33°30′~34°之间。研究区地处北亚热带向暖温带过渡地带，林地面积约为4 280 km²，主要分布于西部和北部山区，耕地面积约为1 470 km², 主要分布于中部和东南部的平原及丘陵地带，其中大部分耕地为旱地，水浇地仅占耕地面积的5%左右。

研究选取3景Landsat系列卫星遥感影像作为数据源，日期分别为1998.09.06 (TM)、2008.09.01(TM)和2017.08.09(OLI和TIRS)(图 1)。所选的3期影像季节相近，避免因植被生长状态不同造成的影响。

图 1(Fig. 1)

图 1 丹江流域(河南段)遥感影像 Fig. 1 Danjiang River Basin(Henan section) remote sensing images

为减少地形、光照和大气等因素对不同时相影像的影像，保证影像在空间叠加分析中的准确性，需要对三期影像进行预处理：1)对影像进行辐射定标，采用Chander和Chavez的模型和参数将影像的灰度值转换为传感器处的反射率^[14-15]；2)为消除大气和光照等因素对地物反射率的影响，需要对辐射定标后的影像进行Flaash大气校正；3)采用二次多项式和最邻近像元法对不同时相的3期影像进行几何校正，使其均方根误差满足小于0.5个像元的精度要求，避免后期进行差值计算时出现误差。

作者采用遥感生态指数(RSEI)监测研究区的生态环境状态及变化情况。该指数由徐涵秋^[10]提出，选取绿度、湿度、热度和干度4个与人类日常生活密切相关的指标并通过主成分分析进行集成，克服了各指标之间的不相关性。遥感生态指数(RSEI)与2006年国家环保部颁发的《生态环境状况评价技术规范》中的生态环境状况指数(EI)^[16-18]相比具有明显的优势，避免了人为的主观因素造成的影响，可以对研究区的生态质量进行可视化表达和建模分析。其中绿度指标的计算方法见文献[19], 湿度指标的计算方法见文献[20-21]，热度指标的计算方法见文献[12, 22-23], 干度指标的计算方法见文献[24]。利用上述方法计算得到的各指标能够反映研究区不同方面的生态环境状况，但是其量纲各不统一，因此在构建遥感生态指数前需要进行归一化处理，避免因权重失衡造成的影响^[25-26]。通用的归一化公式为

式中：NI_i为上述4个指标归一化后的数值；I_ni为该指标在像元n处的数值；I_maxi为该指标的最大值；I_mini为该指标的最小值。归一化后各指标的数值范围应在[0, 1]之间。

笔者采用主成分分析集成归一化处理后的指标数据，主成分分析法(PCA)是一种将多个变量通过正交线性变换来选出少数重要变量的多维数据压缩技术，得到的第1个主分量是对原始多变量数据集贡献最大的新变量。该方法最大的优点就是根据数据本身的性质和各个指标对各主分量的贡献度来自动客观地确定集成各指标的权重，较大地避免了由于主观因素造成的误差，保证结果的可靠性。

由于丹江口水库的水域面积占研究区总面积的4%以上，在进行主成分分析前需要利用改进的归一化差异水体指数(MNDWI)对水体进行掩膜处理^[27]，避免影响PCA的荷载分布。再将归一化后的4个指标合成一副影像并进行主成分分析。

利用ENVI中的主成分分析模块可以得到的第1主成分PC1，进而对PC1进行归一化处理得到遥感生态指数，其计算公式为

式中：PC1_max为第一主成分中的最大值；PC1_min为第一主成分中的最小值。RSEI的值越高，代表研究区的生态状况越好，反之说明生态状况越差。

通过对丹江流域(河南段)的3期Landsat系列影像数据处理，得到其绿度指标、湿度指标、热度指标和干度指标，利用ENVI软件中的主成分分析模块得到的主成分分析结果并统计，利用SPSS软件对RSEI与4个指标之间的相关性进行统计分析。

主成分分析后的统计结果如表 1所示。从表 1中可以看出，研究区3个年份主成分变换后的第1主成分(PC1)贡献率较高，分别为81.83%、76.7%和83.83%，说明变换后的第1主成分已经包含原影像的大部分属性特征，并且绿度和湿度对PC1的贡献值均为正数，说明绿度指标和湿度指标对生态质量有正面的影响，热度和干度对PC1的贡献值均为负数，说明温度指标和干度指标对生态质量存在负面影响。

表 1(Tab. 1)

表 1 主成分分析结果 Tab. 1 Principal component analysis results 指标 Index 1998 2008 2017 PC1 PC2 PC3 PC4 PC1 PC2 PC3 PC4 PC1 PC2 PC3 PC4 NDVI 0.467 0.527 0.452 0.546 0.527 0.417 -0.580 0.460 0.617 0.231 -0.707 0.257 WET 0.365 -0.057 -0.826 0.426 0.248 0.195 0.787 0.531 0.190 0.092 0.500 0.840 LST -0.671 0.676 -0.268 0.146 -0.698 0.706 -0.034 0.117 -0.372 0.928 -0.023 -0.004 NDSI -0.444 -0.513 0.204 0.706 -0.417 -0.538 -0.209 0.702 -0.667 -0.278 -0.499 0.479 特征值Eigenvalues 0.027 9 0.003 4 0.002 4 0.000 4 0.015 9 0.003 4 0.001 2 0.000 3 0.016 3 0.002 5 0.000 6 0.000 012 贡献率Contribution rate/% 81.83 10.09 6.89 1.19 76.70 16.17 5.87 1.26 83.83 12.91 3.20 0.06 注：NDVI为绿度指标，WET为湿度指标，LST为温度指标，NDSI为干度指标。下同。Notes: NDVI is the greenness index, WET is the humidity index, LST is the temperature index, and NDSI is the dryness index. The same below.

表 1 主成分分析结果 Tab. 1 Principal component analysis results

由表 2可见，研究区近20年的遥感生态指数平均值逐步升高，说明研究区总体的生态环境质量不断变好；热度指标和干度指标的荷载值之和的绝对值大于绿度指标和湿度指标的荷载值之和的绝对值，表明热度和干度对研究区生态环境的抑制作用大于绿度和湿度对研究区生态环境的促进作用；通过对比绿度和湿度对PC1的荷载值以及结合遥感生态指数的变化情况，可以得出绿度对研究区生态环境状况的提高作用更为明显。

表 2(Tab. 2)

表 2 遥感生态指数与各指标统计值 Tab. 2 Statistical value of remote sensing ecological index (RSEI) and each index 指标 Index 1998 2008 2017 均值 Mean 方差 Variance PC1荷载 PC1 Load 均值 Mean 方差 Variance PC1荷载 PC1 Load 均值 Mean 方差 Variance PC1荷载 PC1 Load NDVI 0.658 3 0.073 0.467 0.800 6 0.070 0.527 0.900 7 0.084 0.617 WET 0.599 5 0.069 0.365 0.683 5 0.044 0.248 0.861 4 0.014 0.190 LST 0.451 8 0.120 -0.671 0.443 9 0.097 -0.698 0.483 5 0.068 -0.372 NDSI 0.447 0 0.060 -0.444 0.553 9 0.063 -0.417 0.290 0 0.087 -0.667 RSEI 0.632 0 0.127 0.636 6 0.096 0.694 4 0.102

表 2 遥感生态指数与各指标统计值 Tab. 2 Statistical value of remote sensing ecological index (RSEI) and each index

将各年份的遥感生态指数结果以0.2为间隔分为5个级别，分别代表研究区生态差[0~0.2)、较差[0.2~0.4)、中等[0.4~0.6)、良好[0.6~0.8)和优[0.8~1.0]的5个等级，生成遥感生态指数等级图(图 2)并统计各级别面积，结果见表 3。

图 2(Fig. 2)

图 2 遥感生态指数等级空间分布图 Fig. 2 Spatial distribution of RSEI level

表 3(Tab. 3)

表 3 生态质量等级面积和比例 Tab. 3 Area and percentage of ecological quality level 级别 Level 1998 2008 2017 面积Area/km2 比例Percentage/% 面积Area/km2 比例Percentage/% 面积Area/km2 比例Percentage/% 差Poor[0.0~0.2) 2.268 0.028 2.291 0.028 0.653 0.008 较差Inferior[0.2~0.4) 278.038 3.383 305.098 3.713 140.631 1.711 中等Medium[0.4~0.6) 3 144.181 38.260 3 777.028 45.960 1 319.867 16.061 良好Good[0.6~0.8) 4 157.957 50.596 4 036.072 49.113 5 958.461 72.505 优Excellent[0.8~1.0) 635.556 7.734 97.511 1.187 798.388 9.715

表 3 生态质量等级面积和比例 Tab. 3 Area and percentage of ecological quality level

由表 3可见，研究区的生态等级位于良好级别的面积占研究区总面积的比例较大，1998和2008年占的比例在50%左右，2017年达到74.5%，占研究区总面积的近3/4；对其他的几个生态等级进行分析发现：研究区处于中等等级的面积比例呈现先上升后下降的趋势，从1998年的38.26%上升到2008年的45.96%，到2017年迅速下降为16.06%；位于优级别的面积占研究区总面积较少，其比例处于先下降后上升的趋势，2017年最高，达到9.72%；处于差和较差等级的面积所占比例更少，且总体呈现先上升后下降的趋势，最高为2008年的3.66%，最低为2017年的1.72%。综上所述，丹江流域(河南段)总体的生态质量较好，遥感生态指数处于良好级别的较多，2017年研究区的生态质量明显提高，研究区的生态环境已经得到改善。

由图 2可见，生态状况较差的地区主要分布于人口聚集区、耕地以及裸土区域，说明人为的生产活动和植被覆盖状况对生态质量的好坏有明显影响，西部和北部山区的总体生态质量明显高于其他地方。

为了更好地了解丹江流域(河南段)生态环境的变化状况，以遥感生态指数等级图为基础对各年份的遥感生态指数进行差值计算，将得到的结果按区间分为7个等级进行变化检测，包括急剧变差(-1, -0.4]、明显变差(-0.4, -0.2]、略微变差(-0.2, -0.05]、无明显变化(-0.05, 0.05)、略微变好[0.05, 0.2)、明显变差[0.2, 0.4)、急剧变差[0.4, 1)，并统计各变化级别的面积，统计结果见表 4。

表 4(Tab. 4)

表 4 生态质量等级和面积变化统计 Tab. 4 Ecological quality level and statistics of area change 变化级别 Change level 1998—2008 2008—2017 1998—2017 面积 Area/ km2 比例 Percentage/ % 面积 Area/ km2 比例 Percentage/ % 面积 Area/ km2 比例 Percentage/ % 急剧变差Rapidly worse 0.371 0.005 1.335 0.016 0.029 0.001 生态变差Ecology worse 明显变差 Obviously worse 16.665 0.203 24.057 0.293 8.390 0.101 略微变差 Slightly worse 1 632.749 19.868 814.145 9.907 122.215 1.487 生态不变Unchanged 无明显变化No changes 6 125.482 74.537 5 005.831 60.913 4 536.136 55.198 略微变好 Slightly better 431.908 5.256 2 257.548 27.471 3 474.587 42.280 生态变好Ecology better 明显变好 Obviously better 10.782 0.131 112.933 1.374 76.103 0.926 急剧变好 Rapidly better 0.043 0.001 2.151 0.026 0.540 0.007

表 4 生态质量等级和面积变化统计 Tab. 4 Ecological quality level and statistics of area change

由表 4可见，近20年来研究区的生态质量总体处于稳定状态，大部分区域的生态质量无明显变化。1998年到2008年间研究区无明显变化区域的比例占到研究区总面积的74.54%，2008年到2017年间研究区无明显变化区域的面积比例为60.91%，而近20年研究区无明显变化区域的面积比例为55.20%。此外，1998年到2008年间研究区略微变差的面积达到了1 632.749 km²，占研究区总面积的19.87%，略微变好的面积仅为431.908 km²，其他4个级别变化的面积极少，总计面积为27.861 km²，比例仅为0.34%。2008年到2017年间研究区略微变差的区域面积为814.145 km²，而略微变好的面积区域达到2 257.548 km², 占到研究区总面积的27.47%，而明显变好的区域与前10年相比明显上升，达到了112.933 km²，其他变化面积仍然极少。对近20年的变化情况进行分析，发现略微变好的区域面积达到3 474.587 km²，略微变差的面积仅为122.215 km²，其他4类变化面积的比例仅占1.035%。综上所述，从1998—2017研究区生态质量提高的区域所占的面积大于生态质量下降的面积，且上升比例远远高于下降比例，表明研究区的生态质量得到了明显提升。

由图 3(a)可见，1998—2008年间研究区生态质量略微变差的区域较多，且分布于整个研究区，结合国家气象局提供的研究区2年的降雨量数据分析发现，2008年的降雨量明显低于1998年，说明RSEI与降雨量成正相关，与实际情况相符；而略微变好的区域主要在东南部的平原区域，该区域的主要土地利用类型为耕地，说明受到人类活动的影响明显减少。图 3(b)为2008—2017年间研究区生态质量变化图，其中变差的区域主要分布于西峡县、淅川县城区及沿河村镇区域，说明城市扩张是当前区域生态质量变差的主要原因。图 3(c)为1998—2017年的研究区生态质量变化图，生态质量变好的区域集中于东南部及中部的耕地区域，这与退耕还林等国家政策的实施有着密不可分的关系。而变差的区域主要分布于北部的林地和中部的西峡、淅川等县的城市扩张地带，表明城市化进程的加快对区域生态质量变化具有明显的影响，此外北部林地生态质量的降低与人类的活动密切相关，因此需要加强对北部林地的保护，避免水土流失等情况的发生，保护生态系统的稳定性。

图 3(Fig. 3)

图 3 遥感生态指数变化监测图 Fig. 3 RSEI change monitoring map

对丹江流域(河南段)1998—2017年的遥感生态指数分析发现，遥感生态指数均值由1998年的0.632 0上升至2017年的0.694 4，上升幅度达到9.87%，表明研究区的生态质量呈现提升的趋势。其中生态质量变好的面积达到3 551.23 km²，占到研究区总面积43.21%，生态质量变差的面积为130.634 km², 仅占研究区面积的1.59%。在空间格局的变化上，丹江流域(河南段)生态质量变好的区域主要分布于东南部以及中部的耕地区域，生态质量退化的区域主要分布于北部林地以及中部的城区扩张地带。在反映研究区生态质量的4个指标中，绿度指标和干度指标对生态质量的变化影响最为明显，城市的扩张活动会导致区域生态质量显著降低，植被覆盖度的增加对研究区生态质量的提升有明显的促进作用。

利用Landsat系列遥感数据构建的遥感生态指数主要以自然状态指标为主，各指标计算简便且不需要人为设置权重，为大范围、周期性的区域生态状况监测提供了一种新的思路和方法。但是目前仍存在一些不足之处，例如遥感影像的数据质量问题、未结合研究区的经济发展情况等原因；因此在后续的研究中可以考虑以遥感生态指数为基础，结合研究区的经济发展等方面的数据对区域环境状况的协调发展进行监测分析，以期为区域的生态保护提供更加合理的科学依据。