一个用于从空间监测、量化和分类火山活动的协同热框架

《Remote Sensing of Environment》：A synergistic thermal framework for monitoring, quantifying, and classifying volcanic activity from space

　　

结合中红外（MIR）和热红外（TIR）观测捕获了273–1500 K范围内的总辐射热通量。仅基于MIR的清单低估了全球火山热预算2–20倍。VRPRatio可区分热结构，从而实现火山活动分类。该框架在萨班卡亚火山喷发前8个月检测到了热前兆。

本研究针对火山活动卫星监测中单波段热红外遥感方法存在的温度局限性问题，提出了一种结合中红外（MIR）与热红外（TIR）观测的协同热框架。研究背景指出，全球约1400座潜在活动火山中近半数缺乏地面监测，而现有MIR-based火山辐射功率（VRP_MIR）方法仅对>600 K的高温源有效，TIR-based方法（VRP_TIR）适用于<600 K的低温源，但两者独立应用导致全球火山热预算被系统性低估2–20倍，且无法区分热结构变化。研究人员通过合成模拟（10⁶种亚像素温度分布）和VIIRS卫星数据对15座不同活动类型火山（包括Isluga、Nevado del Ruiz、Fuego、La Palma、Nyiragongo等）的2012–2025年时间序列分析，构建了总火山辐射功率（VRP_Tot）和热结构指示因子（VRP_Ratio）两个互补指标。研究得出以下结论：VRP_Tot通过比例因子k≈0.82将MIR和TIR估计合并，在273–1500 K范围内实现温度鲁棒性量化，误差从±30%降至±20%；VRP_Ratio（VRP_MIR/VRP_TIR）作为无量纲有效温度代理，可区分火山活动类型（水热系统<0.1、熔岩穹丘0.2–0.8、开放系统~1、熔岩流1–3、熔岩湖>3）；在Sabancaya火山2012–2025年扰动-喷发序列中，VRP_Ratio在喷发前8个月出现下降趋势（最低值~0.07），揭示了岩浆-水热相互作用导致的热结构前兆信号，而单独VRP_MIR或VRP_TIR序列无法识别。该框架发表于《Remote Sensing of Environment》，对修订全球火山能量预算、实现活动状态定量分类以及早期预警具有重要价值。

关键技术方法：研究人员主要使用了MIROVA算法（基于MIR通道，VIIRS 375m分辨率）获取VRP_MIR，以及TIRVolcH算法（基于TIR通道，VIIRS夜间数据）获取VRP_TIR。通过合成模拟生成10⁶种亚像素温度分布（包括273–1500 K的多种高斯成分），验证了VRP_Tot和VRP_Ratio的物理基础。采用时间序列谐波化处理（线性插值填补缺失值、10点移动平均平滑、剔除天顶角>50°的数据）对齐MIR和TIR数据。样本队列包括来自全球15座火山（如Isluga、Nevado del Ruiz、Fuego、La Palma、Nyiragongo等）的2012–2025年观测数据，覆盖五种活动类型。

研究结果：
3.1 温度鲁棒的总火山辐射功率量化（VRP_Tot）
通过合成模拟，研究人员发现VRP_MIR和VRP_TIR的偏差在结合后相互补偿，VRP_True与（VRP_MIR+VRP_TIR）的比值在273–1500 K范围内保持恒定（k≈0.82±20%），因此定义VRP_Tot = k·(VRP_MIR+VRP_TIR)，消除了单波段温度约束，误差降至±20%。

3.2 VRP_Ratio揭示火山热特征（VTF）的热结构
理论计算和合成模拟表明，VRP_Ratio随有效温度（T_Eff）单调递增，且与VTF面积无关，可作为热结构代理。在VRP_MIR-VRP_TIR参数空间中，沿等值线移动反映面积变化，垂直移动反映温度变化。

3.3 不同活动类型的热体制与能量分配
对15座火山（代表性案例：Isluga水热系统、Nevado del Ruiz熔岩穹丘、Fuego开放系统、La Palma熔岩流、Nyiragongo熔岩湖）的VRP_MIR和VRP_TIR时间序列分析显示：水热系统VRP_Ratio<0.1，TIR占能量95%；熔岩穹丘VRP_Ratio≈0.4，MIR占比30%；开放系统VRP_Ratio≈1，MIR和TIR各半；熔岩流喷发期VRP_Ratio>1，冷却期降至0.1；熔岩湖VRP_Ratio>3，MIR占80%。

3.4 协同热框架用于火山活动分类
在VRP_MIR vs. VRP_TIR参数空间密度分布中，不同活动类型占据不同VRP_Ratio域：水热系统（<0.1）、熔岩穹丘（0.2–0.8）、开放系统（~1）、熔岩流（1–3）、熔岩湖（>3）。冷却熔岩流沿垂直方向迁移，VRP_Ratio下降而VRP_TIR保持高值。

3.5 喷发前水热场热结构的前兆变化：Sabancaya案例
对Sabancaya火山2012–2025年扰动-喷发序列的分析发现，2016年3月VRP_Ratio出现下降（最低~0.07），而VRP_MIR和VRP_TIR同时上升。这一信号由岩浆侵入触发水热系统沸腾导致，增强了低温和扩散性热辐射，在喷发前8个月被检测到，而单独VRP_MIR或VRP_TIR无法识别。喷发开始后VRP_Ratio升至0.3–1，与熔岩穹丘活动一致。

讨论部分总结了框架对全球火山热预算的修订意义：现有MIR-only清单低估2–20倍，尤其水热系统、冷却熔岩等中低温源贡献被忽略，需将全球火山热预算上调2–5倍（中温）至10–20倍（低温）。VRP_Ratio作为热结构诊断指标，可区分不同火山源或过程，如低温柔汽扩散与高温熔岩流的热信号。在Sabancaya案例中，VRP_Ratio下降提供了前兆检测能力，其优势包括定量分类、提前数月预警以及无需火山特定校准的可迁移性。此外，该框架可扩展至行星火山学，如针对木卫一Io的MIR-only低估问题。

研究结论翻译：本研究提出了一个整合MIR-TIR的协同框架，用于获取总火山辐射功率（VRP_Tot）并通过VRP_Ratio无量化指标表征火山热结构。该框架克服了单波段方法的长期温度局限性，表明当前全球火山热预算被显著低估2–20倍，并将VRP_Ratio确立为火山活动及状态转换的一阶诊断指标，已通过15座火山（涵盖全部活动类型）和Sabancaya火山喷发前兆检测得到验证。未来，VRP_Ratio与地震、形变、气体排放等参数的集成，结合机器学习和人工智能方法，将实现多参数系统的自动热结构异常检测和改进的概率性灾害评估。扩展至地球静止平台可提供捕捉快速热转换的时间分辨率，持续应用于行星火山学将有助于完善跨太阳系比较热分析的标准指标。