Lava flow hazard map of Piton de la Fournaise volcano

Piton de la Fournaise, situated on La Réunion island (France), is one of the most active hot spot basaltic shield volcanoes worldwide, experiencing at least two eruptions per year since the establishment of the volcanological observatory in 1979. Eruptions are typically fissure-fed and form extensive lava flow fields. About 95 % of some ∼ 250 historical events (since the first confidently dated eruption in 1708) have occurred inside an uninhabited horseshoe-shaped caldera (hereafter referred to as the Enclos), which is open to the ocean on its eastern side. Rarely (12 times since the 18th century), fissures have opened outside of the Enclos, where housing units, population centers, and infrastructure are at risk. In such a situation, lava flow hazard maps are a useful way of visualizing lava flow inundation probabilities over large areas. Here, we present the up-to-date lava flow hazard map for Piton de la Fournaise based on (i) vent distribution, (ii) lava flow recurrence times, (iii) statistics of lava flow lengths, and (iv) simulations of lava flow paths using the DOWNFLOW stochastic numerical model. The map of the entire volcano highlights the spatial distribution probability of future lava flow invasion for the medium to long term (years to decades). It shows that the most probable location for future lava flow is within the Enclos (where there are areas with up to 12 % probability), a location visited by more than 100 000 visitors every year. Outside of the Enclos, probabilities reach 0.5 % along the active rift zones. Although lava flow hazard occurrence in inhabited areas is deemed to be very low (< 0.1 %), it may be underestimated as our study is only based on post-18th century records and neglects older events. We also provide a series of lava flow hazard maps inside the Enclos, computed on a multi-temporal (i.e., regularly updated) topography. Although hazard distribution remains broadly the same over time, some changes are noticed throughout the analyzed periods due to improved digital elevation model (DEM) resolution, the high frequency of eruptions that constantly modifies the topography, and the lava flow dimensional characteristics and paths. The lava flow hazard map for Piton de la Fournaise presented here is reliable and trustworthy for long-term hazard assessment and land use planning and management. Specific hazard maps for short-term hazard assessment (e.g., for responding to volcanic crises) or considering the cycles of activity at the volcano and different event scenarios (i.e., events fed by different combinations of temporally evolving superficial and deep sources) are required for further assessment of affected areas in the future – especially by atypical but potentially extremely hazardous large-volume eruptions. At such an active site, our method supports the need for regular updates of DEMs and associated lava flow hazard maps if we are to be effective in keeping up to date with mitigation of the associated risks.

Lava flow hazard map of Piton de la Fournaise volcano

&lt;p&gt;Piton de la Fournaise, situated on La Re&amp;#769;union Island (France), is one of the most active hot spot basaltic shield volcanoes worldwide, experiencing at least two eruptions per year since the establishment of the observatory in 1979. Eruptions are typically fissure-fed and form extensive lava flow fields. About 95 % of some ~250 historical events (since the first confidently dated eruption in 1708) have occurred inside an uninhabited horse-shoe shaped caldera (hereafter referred to as the Enclos) which is open to the ocean on its eastern side. Rarely (12 times since the 18th century), fissures have opened outside of the Enclos where housing units, population centers and infrastructure are at risk. In such a situation, lava flow hazard maps are a useful way of visualizing lava flow inundation probabilities over large areas. Here, we present a lava flow hazard map for Piton de la Fournaise volcano based on: i) vent distribution, ii) statistics of lava flow lengths, iii) lava flow recurrence times, and iv) simulations of lava flow paths across multi-temporal (i.e., regularly updated) topography using the DOWNFLOW stochastic numerical model. A map of the entire volcano highlights that the most probable (up to 12 %) location for future lava flow inundation is within the Enclos, where about 100,000 visitors are present each year. Hazard distribution changes throughout the analysis period due to the high frequency of eruptions that constantly modifies the vent opening distribution as well as the topography and the lava flow dimensional characteristics. Outside of the Enclos, probabilities reach 0.5 % along the well-defined rift zones and, although hazard occurrence in inhabited areas is deemed to be very low (&lt;0.1 %), it may be underestimated here, as our study is only based on post-18th century records and neglects cycles of activity at the volcano. Specific hazard maps considering different event scenarios (i.e., events fed by different combinations of temporally evolving superficial and deep sources) are required to better assess affected areas in the future &amp;#8211; especially by atypical, but potentially extremely hazardous, large volume eruptions. At such an active site, our method supports the need for regular updates of DEMs and associated lava flow hazard maps if we are to be effective in mitigating the associated risks.&lt;/p&gt;

Lava flow hazard map of Piton de la Fournaise volcano

Piton de la Fournaise, situated on La Réunion Island (France), is one of the most active hot spot basaltic shield volcanoes worldwide, experiencing at least two eruptions per year since the establishment of the observatory in 1979. Eruptions are typically fissure-fed and form extensive lava flow fields. About 95 % of some ~250 historical events (since the first confidently dated eruption in 1708) have occurred inside an uninhabited horse-shoe shaped caldera (hereafter referred to as the Enclos) which is open to the ocean on its eastern side. Rarely (12 times since the 18th century), fissures have opened outside of the Enclos where housing units, population centers and infrastructure are at risk. In such a situation, lava flow hazard maps are a useful way of visualizing lava flow inundation probabilities over large areas. Here, we present a lava flow hazard map for Piton de la Fournaise volcano based on: i) vent distribution, ii) statistics of lava flow lengths, iii) lava flow recurrence times, and iv) simulations of lava flow paths across multi-temporal (i.e., regularly updated) topography using the DOWNFLOW stochastic numerical model. A map of the entire volcano highlights that the most probable (up to 12 %) location for future lava flow inundation is within the Enclos, where about 100,000 visitors are present each year. Hazard distribution changes throughout the analysis period due to the high frequency of eruptions that constantly modifies the vent opening distribution as well as the topography and the lava flow dimensional characteristics. Outside of the Enclos, probabilities reach 0.5 % along the well-defined rift zones and, although hazard occurrence in inhabited areas is deemed to be very low (

Optimizing barrier placement for lava flow hazard and risk mitigation

&lt;p&gt;Mitigating hazards when lava flows threaten infrastructure is one of the most challenging fields of volcanology, and has an immediate and practical impact on society. Lava flow hazard is determined by the probability of inundation, and essentially controlled by the topography of the area of interest. The most common actions of intervention for lava flow hazard mitigation are therefore the construction of artificial barriers and ditches that can control the flow direction and advancement speed. Estimating the effect a barrier or ditch can have on lava flow paths is non-trivial, but numerical modelling can provide a powerful tool by simulating the eruptive scenario and thus assess the effectiveness of the mitigation action. We present a numerical method for the design of optimal artificial barriers, in terms of location and geometric features, aimed at minimizing the impact of lava flows based on the spatial distribution of exposed elements. First, an exposure analysis collects information about elements at risk from different datasets: population per municipality, distribution of buildings, infrastructure, routes, gas and electricity networks, and land use; numerical simulations are used to compute the probability for these elements to be inundated by lava flows from a number of possible eruptive scenarios&amp;#160; (hazard assessment) and computing the associated economic loss and potential destruction of key facilities (risk assessment). We then generate several intervention scenarios, defined by the location, orientation and geometry (width, length, thickness and even shape) of multiple barriers, and compute the corresponding variation in economic loss. Optimality of the barrier placement is thus considered as a minimization problem for the economic loss, controlled by the barrier placement and constrained by the associated costs. We demonstrate the operation of this system by using a retrospective analysis of some recent effusive eruptions at Mount Etna, Sicily.&lt;/p&gt;

Fast short-term lava flow hazard assessment with graph theory

&lt;p&gt;We explore the use of graph theory to assess short-term hazard of lava flow inundation, with Mt Etna as a case study. In the preparation stage, we convert into a graph the long-term hazard map produced using about 30,000 possible eruptive scenarios calculated by simulating lava flow paths with the physics-based MAGFLOW model. Cells in the original DEM-based representation are merged into graph vertices if reached by the same scenarios, and for each pair of vertices, a directed edge is defined, with an associated lava conductance (probability of lava flowing from one vertex to the other) computed from the number of scenarios that reach both the start and end vertex. In the application stage, the graph representation can be used to extract short-term lava flow hazard maps in case of unrest. When a potential vent opening area is identified e.g. from monitoring data, the corresponding vertices in the graph are activated, and the information about lava inundation probability is iteratively propagated to neighboring vertices through the edges, weighted according to the associated lava conductance. This allows quick identification of potentially inundated areas with little computational time. A comparison with the deterministic approach of subsetting and recomputing the weights in the long-term hazard map is also presented to illustrate benefits and downsides of the graph-based approach.&lt;/p&gt;

Lava flow hazard at Fogo Volcano, Cabo Verde, before and after the 2014–2015 eruption

Lava flow simulations help to better understand volcanic hazards and may assist emergency preparedness at active volcanoes. We demonstrate that at Fogo Volcano, Cabo Verde, such simulations can explain the 2014–2015 lava flow crisis and therefore provide a valuable base to better prepare for the next inevitable eruption. We conducted topographic mapping in the field and a satellite-based remote sensing analysis. We produced the first topographic model of the 2014–2015 lava flow from combined terrestrial laser scanner (TLS) and photogrammetric data. This high-resolution topographic information facilitates lava flow volume estimates of 43.7 ± 5.2 × 106 m3 from the vertical difference between pre- and posteruptive topographies. Both the pre-eruptive and updated digital elevation models (DEMs) serve as the fundamental input data for lava flow simulations using the well-established DOWNFLOW algorithm. Based on thousands of simulations, we assess the lava flow hazard before and after the 2014–2015 eruption. We find that, although the lava flow hazard has changed significantly, it remains high at the locations of two villages that were destroyed during this eruption. This result is of particular importance as villagers have already started to rebuild the settlements. We also analysed satellite radar imagery acquired by the German TerraSAR-X (TSX) satellite to map lava flow emplacement over time. We obtain the lava flow boundaries every 6 to 11 days during the eruption, which assists the interpretation and evaluation of the lava flow model performance. Our results highlight the fact that lava flow hazards change as a result of modifications of the local topography due to lava flow emplacement. This implies the need for up-to-date topographic information in order to assess lava flow hazards. We also emphasize that areas that were once overrun by lava flows are not necessarily safer, even if local lava flow thicknesses exceed the average lava flow thickness. Our observations will be important for the next eruption of Fogo Volcano and have implications for future lava flow crises and disaster response efforts at basaltic volcanoes elsewhere in the world.