Tracking Lava Flow Cooling from Space: Implications for Erupted Volume Estimation and Cooling Mechanisms

Accurate estimation of erupted lava volumes is essential for understanding volcanic processes, interpreting eruptive cycles, and assessing volcanic hazards. Traditional methods based on Mid-Infrared (MIR) satellite imagery require clear-sky conditions during eruptions and are prone to sensor saturation, limiting data availability. Here, we present an alternative approach based on the post-eruptive Thermal InfraRed (TIR) signal, using the recently proposed VRPTIR method to quantify radiative energy loss during lava flow cooling. We identify thermally anomalous pixels in VIIRS I5 scenes (11.45 µm, 375 m resolution) using the TIRVolcH algorithm, this allowing the detection of subtle thermal anomalies throughout the cooling phase, and retrieve lava flow area by fitting theoretical cooling curves to observed VRPTIR time series. Collating a dataset of 191 mafic eruptions that occurred between 2010 and 2025 at (i) Etna and Stromboli (Italy); (ii) Piton de la Fournaise (France); (iii) Bárðarbunga, Fagradalsfjall, and Sundhnúkagígar (Iceland); (iv) Kīlauea and Mauna Loa (United States); (v) Wolf, Fernandina, and Sierra Negra (Ecuador); (vi) Nyamuragira and Nyiragongo (DRC); (vii) Fogo (Cape Verde); and (viii) La Palma (Spain), we derive a new power-law equation describing mafic lava flow thickening as a function of time across five orders of magnitude (from 0.02 Mm3 to 5.5 km3). Finally, from knowledge of areas and episode durations, we estimate erupted volumes. The method is validated against 68 eruptions with known volumes, yielding high agreement (R2 = 0.947; ρ = 0.96; MAPE = 28.60%), a negligible bias (MPE = −0.85%), and uncertainties within ±50%. Application to the February-March 2025 Etna eruption further corroborates the robustness of our workflow, from which we estimate a bulk erupted volume of 4.23 ± 2.12 × 106 m3, in close agreement with preliminary estimates from independent data. Beyond volume estimation, we show that VRPTIR cooling curves follow a consistent decay pattern that aligns with established theoretical thermal models, indicating a stable conductive regime during the cooling stage. This scale-invariant pattern suggests that crustal insulation and heat transfer across a solidifying boundary govern the thermal evolution of cooling basaltic flows.