编辑推荐：
研究人员研究了六条智利中部的冰川物质平衡序列的极端行为，这些冰川因气候变率和变化而面临显著退缩与冰体损耗，导致旱季融水减少并加剧下游水资源竞争。理解冰川物质平衡对预测未来高需水情景下的水资源可用性至关重要。研究人员采用极值理论(Extreme Value Th
研究人员研究了六条智利中部的冰川物质平衡序列的极端行为，这些冰川因气候变率和变化而面临显著退缩与冰体损耗，导致旱季融水减少并加剧下游水资源竞争。理解冰川物质平衡对预测未来高需水情景下的水资源可用性至关重要。研究人员采用极值理论(Extreme Value Theory, EVT)工具分析数据并识别极端事件，主要目标是对智利冰川物质亏损进行统计分析，所用物质平衡数据来源于国内外公开资料。结果表明各冰川极端行为高度异质：部分冰川呈近似指数尾（Guanaco冰川），部分呈稳定且弱尾特征（Echaurren Norte与Mocho冰川），一条冰川（Amarillo冰川）具高度不稳定结构，其余（Juncal Norte与Juncal Sur冰川）呈弱且可能有限尾。结果证实深入研究安第斯山脉冰川对理解其气候变化响应及未来冰川融水资源的流域管理具有重要意义。
论文解读：《On the Mean Excess Plot Measures of Chilean Glacier Mass Balance Data》
研究背景与立项依据
智利拥有约24,000 km2冰川面积，占南美洲冰川总面积的76%，是干旱及半干旱地区关键的水资源来源。近几十年来冰川持续退缩与分裂，尤其在中央山谷（Valparaíso、Metropolitana、O'Higgins及Maule大区），如Juncal Sur冰川在1955–2018年间损失大量面积。现有研究多依赖线性趋势分析，无法刻画极端气候事件（如2010年以来的超级干旱megadrought）引发的非线性加速物质流失，且存在观测数据稀缺、空间变异性大等局限。极值理论(Extreme Value Theory, EVT)可专门建模罕见高影响事件，弥补常规趋势分析的不足。本研究旨在运用EVT中的均值超出图(Mean Excess Plot/Mean Excess Function, MEF)及峰值超过阈值(Peaks Over Threshold, POT)方法，定量表征智利六条典型冰川（Echaurren Norte、Guanaco、Amarillo、Mocho、Juncal Norte、Juncal Sur）年度物质平衡的负极端（极端物质亏损）分布形态，为水资源风险管理提供统计依据。本文发表于《Water》(MDPI)。
主要技术方法概述
研究人员收集世界冰川监测局(World Glacier Monitoring Service, WGMS) 2025版数据库及智利水资源总局(Dirección General de Aguas, DGA)公布的Juncal Norte冰川资料，以及文献中间接来源的Juncal Norte与Juncal Sur历史物质平衡记录，单位为米水当量每年(m w.e./year)。对直接数据筛选2000–2024年共有时段，剔除缺失值与异常值后构建年度物质平衡时间序列。为分析左尾（极端亏损），将负物质平衡取绝对值转化为正值。采用经验均值超出函数(empirical Mean Excess Function, MEF)以中位数初设阈值绘制均值超出图，判断尾部分布类型（正斜率→重尾/heavy tail，ξ>0；零斜率→指数尾/exponential tail，ξ≈0；负斜率→短尾/short tail或bounded tail，ξ<0）。辅以指数QQ图、超值变异系数(Coefficient of Variation, CV)平稳性检验（Castillo & Padilla法，检验广义帕累托分布Generalized Pareto Distribution, GPD拟合适用性），以及Hill估计量、t-Hill估计量和Zipf估计量分析尾指数(tail index, 1/α或ξ)。POT模型选用广义帕累托分布(GPD)拟合超出量，通过MEF线性段选取阈值u，GPD形状参数(shape parameter, ξ)表征尾重。使用R 4.1.1与Python 3.12.13进行数据分析与绘图。
研究结果（按原文Results小节）
3.1. Direct Data（直接数据/WGMS来源）
Echaurren Norte冰川有49条完整年度记录（2000–2024年分析），年均物质亏损?1.27 m w.e.，标准差1.00；Guanaco冰川年均亏损?0.62 m w.e.，标准差0.40（最稳定）；Amarillo冰川年均亏损?0.29 m w.e.但标准差达1.66（高年际波动）；Mocho、Toro 1/2、Esperanza均值介于?1.10至?1.24 m w.e.。Juncal Norte与Tapado无有效直接数据。直接数据显示各冰川均为净物质亏损趋势。
3.2. Indirect Data（间接数据/Juncal Norte与Juncal Sur）
Juncal Norte（Aconcagua河流域）与Juncal Sur（Maipo河流域）采用文献历史记录补充，二者均呈负物质平衡趋势，Juncal Norte亏损幅度略大于Juncal Sur，与历史报道一致。
3.3. Coefficient of Variance Analysis（变异系数分析）
以50%~95%分位数初筛阈值，进阶检验95%~99%分位数。Guanaco与Amarillo的CV统计量接近0（表明中阈度下CV较稳），Echaurren Norte与Mocho稍大但仍兼容POT。各冰川在0.95~0.99最高分位段CV稳定性下降（样本少致估参不确定增加）。间接数据的Juncal Norte与Juncal Sur在中段阈值CV≈稳定，高阈段同样失稳。
3.4. Heavy Tail Analysis（重尾分析——Hill/t-Hill/Zipf估计量）
对负物质平衡取绝对正值后应用三种尾指数估计量。Guanaco在k≈n/4内Hill与Zipf曲线呈稳定平台，t-Hill后段略升，指示近似指数尾（ξ≈0）。Echaurren Norte有小段稳定区但t-Hill随k增大明显上翘，反映中度重尾（ξ略>0但偏轻）。Mocho类似Echaurren Norte但t-Hill发散更早。Amarillo三估计量均无明确稳定平台，t-Hill在较小k后剧增，表明尾结构不稳定、难归属经典Pareto或指数族。
讨论与结论翻译总结
讨论部分指出：MEF线性拟合均方误差(Mean Squared Error, MSE) Guanaco最低(≈0.0018，高度线性→GPD适配佳，近似指数尾ξ≈0)，Mocho次之(≈0.0033)，Echaurren Norte中等(≈0.0128，弱重/轻尾ξ≤0)，Amarillo最高(≈0.0451，MEF非线→极不稳定)。GPD形状参数ξ归纳——Guanaco: ξ≈0（指数尾）；Echaurren Norte与Mocho: ξ≈0或略负（轻尾/light tail）；Juncal Norte与Juncal Sur: ξ<0（短尾/potentially bounded tail）；Amarillo: 无可靠ξ估计（结构不稳定）。安第斯冰川极端响应高度异质，需逐冰川单独做EVT诊断；Guanaco最适POT-GPD建模，Amarillo需谨慎对待阈值选取及极端风险推估。间接数据（Juncal系列）尾更短、阈值敏感性强。
结论部分原文翻译：
本研究利用极值理论(EVT)工具，在峰值超过阈值(POT)方法与广义帕累托分布(GPD)拟合框架下，分析了智利Guanaco、Echaurren Norte、Mocho、Amarillo、Juncal Norte及Juncal Sur冰川的物质平衡极端行为。通过均值超出函数(MEF)、指数QQ图、变异系数(CV)图及Hill、t-Hill与Zipf尾指数估计量综合评价了模型适用性。MEF分析及线性拟合MSE揭示各冰川极端行为差异明显：Guanaco冰川MSE≈0.0018显示高度稳定的尾结构与GPD优良适配（近似指数尾）；Mocho冰川MSE≈0.0033表现亦稳健；Echaurren Norte冰川MSE≈0.0128反映可接受但变异性较大的稳定尾；Amarillo冰川MSE≈0.0451表明平均超出函数非线性及极端结构高度不稳定。Guanaco冰川表现为持续规整的物质亏损、平稳的CV图及与理论指数分布吻合良好的QQ图，对应形状参数ξ≈0。Echaurren Norte与Mocho冰川MEF递减、CV低于1，属轻尾（ξ≈0或略负）。Amarillo冰川MEF振荡剧烈、CV图无平稳段、QQ图显著偏离，无法在POT/GPD框架内可靠赋予ξ值。Juncal Norte与Juncal Sur（间接数据）CV单调递减且远低于1，暗示短尾（ξ<0），POT适用但阈值选取敏感。研究证实智利（尤其中安第斯）冰川极端物质平衡行为具高度异质性，须逐个冰川开展EVT诊断，成果为安第斯冰川极端事件建模、水文气候预估及变化环境下水资源适应规划提供了统计学基础。
