The first complete volcano-structural map of El Hierro (Canary Islands, Spain) has been developed in order to provide a tool for volcano-tectonic analyses and volcanic hazard evaluation on the island. This map is a synthesis of collated and interpreted field data and bathymetric maps. We have integrated information obtained from: (1) high-resolution digital elevation models, (2) bathymetric information, (3) topographic, geologic and geomorphological maps, (4) aerial photographs and orthoimages, (5) previous reports and scientific publications, and (6) new detailed field surveys. The 1:100,000-scale map includes the main volcano-structural elements such as vents, eruptive fissures, dykes, faults, and landslides scars. This information has been used for analysing the volcano-tectonic evolution of El Hierro and for estimating the probability of new vent opening on the island (i.e. volcanic susceptibility). We expect that this map will underpin future geological studies and future volcanic risk assessment.

This study used spatiotemporal land gravity data to investigate the 2021 eruption that occurred in the Cumbre Vieja volcano (La Palma, Canary Islands). First, we produced a density model by inverting the local gravity field using data collected in July 2005 and July 2021. This model revealed a low-density body beneath the western flank of the volcano that explains a highly fractured and altered structure related to the active hydrothermal system. Then, we retrieved changes in gravity and GNSS vertical displacements from repeated measurements made in a local network before (July 2021) and after (January 2022) the eruption. After correcting the vertical surface displacements, the gravity changes produced by mass variation during the eruptive process were used to build a forward model of the magmatic feeding system consisting of dikes and sills based on an initial model defined by the paths of the earthquake hypocenters preceding the eruption. Our study provides a final model of the magma plumbing system, which establishes a spatiotemporal framework tracing the path of magma ascent from the crust–mantle boundary to the surface from 11–19 September 2021, where the shallowest magma path was strongly influenced by the low-density body identified in the inversion process.

Risk management stakeholders in high-populated volcanic islands should be provided with the latest high-quality volcanic information. We present here the first volcanic susceptibility map of Lanzarote and Chinijo Islands and their submarine flanks based on updated chronostratigraphical and volcano structural data, as well as on the geomorphological analysis of the bathymetric data of the submarine flanks. The role of the structural elements in the volcanic susceptibility analysis has been reviewed: vents have been considered since they indicate where previous eruptions took place; eruptive fissures provide information about the stress field as they are the superficial expression of the dyke conduit; eroded dykes have been discarded since they are single non-feeder dykes intruded in deep parts of Miocene-Pliocene volcanic edifices; main faults have been taken into account only in those cases where they could modified the superficial movement of magma. The application of kernel density estimation via a linear diffusion process for the volcanic susceptibility assessment has been applied successfully to Lanzarote and could be applied to other fissure volcanic fields worldwide since the results provide information about the probable area where an eruption could take place but also about the main direction of the probable volcanic fissures.