《Water》：Screening Potential Control-Ring Interference in Remote-Sensing Assessments of Cascade Reservoir Microclimates: A Cautionary Case from Yunnan, China
编辑推荐：
遥感评估水库微气候的研究常将水库周边区域与周围对照区进行比较，而对照区应既未受处理影响又与环境条件相当。本研究在一组密集的梯级水库系统中，通过几何筛查情景和示例性结果分析评估了上述要求。针对中国云南的62座水库，96.8%的主要20–50 km对照环与至少一座
遥感评估水库微气候的研究常将水库周边区域与周围对照区进行比较，而对照区应既未受处理影响又与环境条件相当。本研究在一组密集的梯级水库系统中，通过几何筛查情景和示例性结果分析评估了上述要求。针对中国云南的62座水库，96.8%的主要20–50 km对照环与至少一座邻近水库的4倍水面面积（SWA）缩放代理缓冲区相交；在2倍–5倍半径情景下，该比例为88.7%–98.4%，而两个更宽的对照环则达到100%。由于这些缓冲区属于未经验证的几何体，重叠仅被解读为警示信号，而非物理溢出的证据。在一项包含34座水库的案例研究（2000–2024年）中，原始对照环的最高温度（TMAX）原始估计值为?0.048 °C（95%置信区间[?0.073, ?0.023]，p = 0.00042）。在排除代理重叠像元后，该估计值变化不大（?0.044 °C），但在进行高程–坡度平衡后减弱（?0.025 °C，95%置信区间[?0.053, 0.002]，p = 0.065），在联合排除与平衡后进一步减弱（?0.023 °C，95%置信区间[?0.049, 0.004]，p = 0.092）。去趋势估计在全部四种对照条件下均不显著，其余六项结果亦不具备规格稳定性。日期精度情景保持了负向原始方向，但未能解决趋势敏感性。因此，该案例展示的是对照设计的脆弱性，而非确定性的水库冷却效应。在进行因果解读之前，几何暴露筛查应配合地形诊断、替代对照以及明确的误差报告。
本研究聚焦于梯级水库遥感微气候评估中对照区选择的方法论问题，系统揭示了密集梯级水库系统中传统等距对照环设计所面临的双重威胁：其一，某一水库的对照环可能与其他水库的潜在暴露区相交，导致近/远对比所比较的位置受到不同程度的处理暴露；其二，等距对照环可能跨越山脊和高程带，使其气候条件与焦点水库的近区环境不具可比性。上述问题源于稳定单元处理价值假设（SUTVA）的潜在违反，在沿同一河流密集布设的水坝系统中尤为突出。中国西南地区，特别是云南省的金沙江、澜沧江和怒江流域，梯级水电开发迅速扩张，金沙江中下游干流即建有11座大型水坝，这为检验对照环设计可靠性提供了理想的研究场景。既有研究虽已报道水库蓄水后局部温度、降水、湿度及风速的变化，但普遍将距离焦点水库区域之外的对照区视为适用，未系统评估邻近水库处理暴露和地形可比性问题。
研究人员利用中国大型水库属性与流域边界数据集（CLEAR v1.2）确定云南省62座候选水库样本，基于MODIS NDWI时序数据判定蓄水时间，采用ChinaMet 0.01°网格化气象强迫数据获取温度、降水等气象变量，并借助Copernicus DEM GLO-30数据进行地形诊断。研究构建了三种近区/对照环设计方案（0–20 km对20–50 km、0–30 km对30–80 km、0–50 km对50–100 km），并针对全国913座水库生成水面面积（SWA）缩放代理缓冲区，在2倍至5倍半径乘数情景下评估对照环与邻近水库代理缓冲区的几何重叠。在34座具有可用蓄水日期的水库中，研究人员开展2000–2024年微气候案例研究，对最高温度（TMAX）、最低温度（TMIN）、潜在蒸散（PET）、相对湿度（RHU）、降水（PREC）、植被净初级生产力（NPP）和NDWI七项结果变量进行前后对比分析，并设计了四种预设统计规格及多种控制定义。
控制环代理缓冲区重叠广泛存在。在主要对照环设计下，62座候选水库中有60座（96.8%）的对照环与至少一座邻近水库的SWA缩放代理缓冲区相交，并集重叠面积占对照环面积的均值为4.71%，中位数为3.27%。在0–30 km对30–80 km和0–50 km对50–100 km设计下，计数层面重叠率均达到100%，并集面积占比均值分别为5.05%和5.78%。坝点敏感性检验显示主要对照环的重叠率为82.3%，但总体格局一致。乘数敏感性分析表明，主要对照环的重叠率从2倍因子下的88.7%递增至5倍因子下的98.4%。该模式并非局限于金沙江梯级最密集河段，而是出现在全部五个流域，澜沧江、元江和南盘江流域的重叠率均达到100%，金沙江流域为95.8%，瑞丽江/南宛河（伊洛瓦底江）流域为83.3%，表明这是梯级系统普遍存在的结构性问题而非单一河段的偶然现象。
警示性案例研究显示TMAX信号对对照设计和趋势处理敏感。在主要对照环的原始全部规格下，TMAX平均水库级效应为?0.048 °C（95%置信区间[?0.073, ?0.023]，p = 0.00042），呈现名义上显著的冷却信号。在clean75和strict50预趋势筛选子集中，信号虽保持负向但幅度减小，分别为?0.029 °C和?0.032 °C。然而，在去除各水库蓄水前自身线性趋势后，估计效应减弱至?0.014 °C且不再显著（p = 0.5878）。其余六项变量未呈现一致的相同符号显著模式，经Benjamini–Hochberg多重性调整后仅三个孤立单元低于q = 0.05，去趋势单元无一通过调整。替代对照定义显著改变了原始TMAX结果：排除代理重叠像元后估计值为?0.044 °C，与原始结果接近；联合高程–坡度平衡将估计值减弱至?0.025 °C（p = 0.065），结合排除与平衡后为?0.023 °C（p = 0.092），两者置信区间均包含零。去趋势估计在全部四种对照定义下均不显著（p = 0.59–0.80）。等河流加权汇总呈现相同解读，原始对照环均值为?0.040 °C，地形平衡后减弱至?0.018 °C，联合控制下为?0.014 °C，全部去趋势等河流区间均包含零。日期洗牌安慰剂检验中，500次置换中仅1次原始效应达到观测值的极端程度（经验p = 0.004），而去趋势效应则完全落在安慰剂分布之内（经验p = 0.679）。日期精度敏感性分析显示，全局移位?2至+2个月保持原始均值负向（?0.048至?0.023 °C），但去趋势估计在?0.030至+0.018 °C范围内均不显著；限制于16个日精度记录时原始均值缩小至?0.032 °C且不显著（p = 0.093）。
讨论部分指出，几何重叠筛查与协变量平衡评估需分别进行，两者均不足以单独确立平行趋势假设，但共同使用可使近/远比较的假设更易评估。代理排除单独未能解释表观冷却效应，而地形平衡产生了更大影响，表明控制区环境可比性是更关键的限制因素。研究局限性包括：SWA缩放代理半径属于假设性筛查距离而非经校准的影响边界；地形加权未考虑坡向、土地覆盖、山谷朝向等因素；34座水库分布于18条河流且14条为单水库河流，无法充分估计河内相关性；蓄水日期精度不一且可能指向不同的运行里程碑；存档月度NDWI合成产品未重新处理像元级质量保证层。研究建议采用分阶段对照选择方法：先以透明半径范围报告代理重叠筛查诊断，再检查地形及其他基线差异，最后在可行的替代对照下重新估计结果。该关切同样适用于聚集性可再生能源设施、保护区网络、矿业开发和城市扩展廊道等研究场景，距离本身不能确立可信的比较，对照必须既合理未暴露又环境可比，否则近/远及前后估计应作为描述性证据而非因果效应呈现。
结论部分翻译：本研究在密集梯级背景下，使用显式筛查情景评估了传统对照环。在62座水库的云南样本中，96.8%的主要对照环与至少一座邻近水库的4倍SWA缩放代理缓冲区相交，更宽的对照环达到100%；在2倍–5倍乘数下该比例仍保持高位。这些数值表明几何筛查风险，而非实测物理干扰。示例性34座水库案例展示了对照构建为何重要。原始对照环TMAX估计为负向，但在高程–坡度平衡或联合平衡与代理排除后减弱并失去显著性。去趋势估计在全部对照定义下均不显著，其余六项结果不具备规格稳定性。河流结构、日期精度及未重新处理的MOD09A1像元质量保证进一步限制了因果解读。因此，遥感DID、BACI及相关设计应记录对照区为何合理未暴露且可比，而非仅与焦点干预有一定距离。代理缓冲区筛查可识别候选风险，但因果解读需要地形与协变量诊断、替代对照估计及明确的误差报告。本研究结果宜解读为对照设计的警示案例，而非固定水库微气候效应的证据。
