《Heliyon》：Deciphering the fate of Gahar Lake: Integrating remote sensing and climate models for future projections
编辑推荐：
全球淡水湖泊日益受到气候变化与人为压力叠加效应的威胁，凸显了建立综合监测与预测框架以保护这些重要生态系统的迫切性。Gahar Lake是伊朗一个关键淡水生态系统，近几十年来经历了严重退化，正是这些挑战的典型例证。本研究对2001年至2023年期间该湖泊的水面面
全球淡水湖泊日益受到气候变化与人为压力叠加效应的威胁，凸显了建立综合监测与预测框架以保护这些重要生态系统的迫切性。Gahar Lake是伊朗一个关键淡水生态系统，近几十年来经历了严重退化，正是这些挑战的典型例证。本研究对2001年至2023年期间该湖泊的水面面积动态进行了综合评估，并预测了未来的气候驱动变化因素。研究人员利用谷歌地球引擎（Google Earth Engine, GEE）平台中的多传感器卫星影像（Landsat 7、Landsat 8和Sentinel-2），应用多种水文指数（归一化差异水体指数（Normalized Difference Water Index, NDWI）、改进归一化差异水体指数（Modified Normalized Difference Water Index, MNDWI）和归一化差异湿度指数（Normalized Difference Moisture Index, NDMI））追踪地表水变化。分析结果显示湖泊面积急剧缩减，MNDWI显示其从约190公顷下降至约64公顷。Landsat 7时间序列的趋势分析证实MNDWI指数存在统计学显著的负斜率（p < 0.1），而多传感器合并数据集呈不显著下降趋势。为预测未来状况，研究人员在高温室气体排放情景——代表性浓度路径8.5（Representative Concentration Pathway 8.5, RCP8.5）下，使用Long Ashton Research Station天气发生器（LARS-WG）模型对2021-2040年期间进行预测。预测结果表明威胁将持续存在，预计降水量将减少15%，其中5月份最为显著，同时冬季和春季气温将上升4%-13%。这些气候变迁预计将加剧湖泊退化，威胁其生物多样性和生态系统服务。研究结果证明了将遥感与气候模型相结合以监测和预测地表水动态的有效性。本研究为制定Gahar Lake的主动保护策略提供了关键数据，并可作为评估全球脆弱水生生态系统的方法框架。
一、研究背景与问题
淡水湖泊是全球水资源系统的重要组成部分，但正日益遭受气候变化与人为活动的双重胁迫。Gahar Lake位于伊朗洛雷斯坦省奥什托兰库赫（Oshtorankuh）保护区内、扎格罗斯山脉腹地，海拔2350米，是该国独特的高山湖泊之一，由大地震及滑坡活动形成，分为大Gohar（水深4-28米，长1800米）和小Gohar（水深1-3米，长500米）两部分。湖泊依赖地下泉水与河流补给，入流量280升/秒，出流量570升/秒，区域年均降水量约933毫米，降水多以积雪形式储存于高山，春季融雪形成季节性溪流。然而近几十年来，该湖出现严重退化。周边Borujerd-Dorud平原约12500公顷的灌溉农业在旱季大量取水，季节性旅游也加剧了水资源压力。此前研究存在两大方法缺陷：一是缺乏多传感器卫星数据的综合运用——Landsat与Sentinel-2具有互补的空间（10-30米）和时间（5-16天）分辨率，尤其适用于云覆盖多、地形阴影复杂的山区；二是未将历史地表水动态与未来气候预测进行定量关联。为此，研究人员开展了这项集成遥感与气候模型的综合评估研究，旨在为这一脆弱生态系统的适应性管理提供科学依据，论文发表在《Heliyon》。
二、主要技术方法
研究基于谷歌地球引擎（Google Earth Engine, GEE）云平台，整合了NASA/USGS的Landsat 7（2001-2012）、Landsat 8（2013-2014）和ESA的Sentinel-2（2015-2023）多传感器卫星影像，统一进行云掩膜与地表反射率校正。应用三种光谱指数——归一化差异水体指数（NDWI）、改进归一化差异水体指数（MNDWI）和归一化差异湿度指数（NDMI）提取水体范围，并采用Mann-Kendall检验与Sen斜率估计器进行趋势显著性分析。气候预测方面，使用Long Ashton Research Station天气发生器（LARS-WG）v6模型，对EC-EARTH全球气候模式（GCM）输出进行降尺度，在RCP8.5高排放情景下生成2021-2040年的逐日气象数据。历史气象数据来自伊朗气象组织（IRIMO）邻近站点，经Kriging插值估算湖区气象要素。
三、研究结果
3.1 基于NDWI、MNDWI和NDMI指数的Gahar Lake水面面积变化
NDWI结果显示，湖泊面积在2001年约180公顷、2002年约188公顷，2012年异常增至约385公顷（可能由雪与水在可见光和近红外波段的光谱相似性造成），此后波动下降至2023年约59公顷，线性斜率为-0.7。MNDWI结果显示，湖泊面积从2001年约190公顷降至2023年约64公顷，拟合方程斜率为-0.19，证实了明确的下降趋势。NDMI结果则呈现较大波动，从2016年9公顷增至2018年86.5公顷，2021年达约138公顷，随后降至2023年约67公顷，反映出该指数对土壤湿度的高度敏感性。
3.1.1 光谱指数误差与不确定性分析
误差分析表明，NDWI平均误差最大（22.4公顷），因其对植被和土壤湿度敏感；MNDWI表现稳定（平均误差14.6公顷），利用短波红外（SWIR）波段有效降低了非水体要素干扰；NDMI虽平均误差最小（9.2公顷），但在干旱年份接近零值、湿润年份高估，不适宜用于开放水体监测。从传感器角度，Landsat 7时期误差最大（NDWI平均29.8公顷），Sentinel-2时期误差显著降低（NDWI平均13.2公顷）。综合比较后，研究推荐MNDWI结合Sentinel-2数据作为Gahar Lake长期监测的最佳方案。
3.2 Gahar Lake空间变化
基于MNDWI的秋季影像显示，2001年湖泊面积约190公顷；2016年缩减至约78公顷；2017年因有利的降雪和降雨条件回升至约123公顷；2023年再次收缩至约64公顷。湖泊边界萎缩在影像上清晰可见。
3.3 线性趋势分析
NDWI（y = ?0.7047x + 85.997）和MNDWI（y = ?0.1988x + 72.78）的线性方程斜率均为负，表明湖泊面积呈下降趋势。NDMI因Sentinel-2时期数值整体高于Landsat时期而呈现正斜率（y = 2.6713x + 13.885），属于传感器间系统性差异所致。
3.4 基于Mann-Kendall和Sen检验的Gahar Lake水面面积趋势分析
Mann-Kendall检验显示，MNDWI基于Landsat 7的数据序列在90%置信水平上呈显著下降趋势（Z = ?1.95），Sen斜率为?1.57公顷/年；多传感器合并数据呈非显著下降趋势（Z = ?0.24，Sen斜率?0.28公顷/年），差异可能源于传感器分辨率和定标差异。NDWI在Landsat 7与8合并数据中呈99%置信水平的显著下降（Z = ?2.65）；NDMI在合并数据中呈95%置信水平的显著上升（Z = 2.53），反映其对周边土壤水分的响应而非水体面积变化。
3.5 2021-2040年气候变量变化
RCP8.5情景下，EC-EARTH模式预测显示：5月降水较基准期减少15%，4月减少2%，其余月份变化较小；最高气温在1月和2月增幅最显著，最低气温在3月增幅最大（达52%），冬季和春季整体升温幅度为4%-13%。空间分布上，研究区东南部最高温呈下降趋势，西北部最低温上升，这些变化可能与大气环流、植被覆盖和温室气体浓度升高有关。
3.6 人为与水文背景
湖流域旱季（5-9月）农业取水、无序旅游（废弃物排放、物种干扰）对湖泊构成额外压力。有限的径流数据（2012-2019）显示下降趋势，与湖泊面积萎缩一致，但受数据长度限制无法定量分离气候与人为因素的相对贡献。
四、讨论与结论
讨论部分指出，Gahar Lake过去二十年损失66%的湖面面积，退化速率高于伊朗Maharloo湖（64%）、Mighan湿地（56%）、Lake Urmia（50%）以及土耳其Burdur湖（41.77%）、Beysehir湖（4.68%）和Egirdir湖（5.44%），与扎格罗斯山脉雪盖年减少0.25%-0.68%的趋势相吻合。气候预测结果与印度（降水减少12%-18%）、中国沿海（8%-15%）及其他干旱半干旱地区基于LARS-WG的研究具有数值一致性。研究存在多传感器集成不确定性、单一GCM与单一情景、缺乏实地验证数据及云覆盖导致的时间缺口等局限。
研究结论可翻译概括为：本研究综合运用多传感器卫星数据（Landsat 7、Landsat 8和Sentinel-2）与多种遥感指数（NDWI、MNDWI、NDMI），在谷歌地球引擎平台上评估了Gahar Lake 2001-2023年的长期水面面积动态，并结合LARS-WG模型完成RCP8.5情景下的未来气候预测（2021-2040）。核心发现包括：（1）MNDWI显示湖泊面积从约190公顷减至约64公顷，损失近66%；（2）Mann-Kendall检验确认90%置信水平上的显著下降趋势（Z = ?1.95），年均损失1.57公顷；（3）NDWI和MNDWI均呈负趋势，而NDMI呈上升趋势（Sen斜率2.23公顷/年），反映其对土壤水分的敏感性；（4）RCP8.5情景下5月降水减少15%、冬春气温升高4%-13%，将加速湖泊蒸发与水位下降；（5）湖面萎缩趋势与伊朗其他水体及扎格罗斯山脉雪盖减少的研究结论具有区域一致性。建议将Gahar Lake列为高优先级保护区，建立卫星与实地水位计相结合的监测系统，规范旱季河流取水，并制定适应气候变化（雪盖减少与
