Search

The potential of the island of La Palma (Canary Islands) to host geothermal resources is very high, mainly due to its high volcanic activity. The primary goal of this study is to get a tridimensional image of the seismic intrinsic attenuation using ambient seismic noise and to identify anomalies that may be linked to active geothermal reservoirs on La Palma island. For this purpose, we developed a new Ambient Noise Attenuation Tomography (ANAT) technique, which uses seismic ambient noise for imaging intrinsic attenuation in 3-D at a local scale down to 5 km depth. Our research identifies two areas with high attenuation in the island's southern region. One area could be associated with hydrothermal alteration zones beneath the Cumbre Vieja volcanic complex. Another high-attenuation zone was observed in the island's southern part, which could be associated with extensively fractured rocks that might facilitate the circulation of heated fluids. Furthermore, we discuss the geothermal relevance of such anomalies, making a comparison with previous resistivity, S-wave velocity and density models. This study confirms that intrinsic attenuation retrieved from the coda of Rayleigh waves is more sensitive to the presence of fluids than velocity. Fluids being a key component of active geothermal reservoirs, it is reasonable to expect intrinsic attenuation anomalies in these systems. Therefore, we conclude that ANAT can be useful for geothermal exploration., (© 2024. The Author(s).)

La Palma island is located in the NW of the Canary Islands and is one of the most volcanically active of the archipelago, therefore the existence of geothermal resources on the island is highly probable. The main objective of this work is to detect velocity anomalies potentially related to active geothermal reservoirs on La Palma island, by achieving a high-resolution seismic velocity model of the first few kilometres of the crust using Ambient Noise Tomography (ANT). The obtained ANT model is merged with a recent local earthquake tomography model. Our findings reveal two high-velocity zones in the island's northern and southern parts, that could be related to a plutonic intrusion and old oceanic crust materials. Conversely, four low-velocity zones are imaged in the southern part of the island. Two of them can be related to hydrothermal alteration zones located beneath the Cumbre Vieja volcanic complex. This hypothesis is reinforced by comparing the S-wave velocity model with the seismicity recorded during the pre-eruptive phase of the 2021 Tajogaite eruption, which revealed an aseismic volume coinciding with these low-velocity zones. Another low-velocity zone is observed in the southern part of the island, which we interpret as highly fractured rocks which could favour the ascent of hot fluids. A last low-velocity zone is observed in the central part of the island and associated with loose deposits generated by the Aridane valley mega landslide., (© 2023. Springer Nature Limited.)

On Sept. 19th, 2021, a volcanic eruption began on the island of La Palma (Canary Islands, Spain). The pre-eruptive episode was characterized by seismicity and ground deformation that started only 9.5 days before the eruption. In this study, we applied seismic interferometry to the data recorded by six broadband seismic stations, allowing us to estimate velocity variations during the weeks preceding the eruption. About 9.5 days before the eruption, we observed a reduction in the seismic velocities is registered next to the eruptive centers that opened later. Furthermore, this zone overlaps with the epicenters of a cluster of volcano-tectonic earthquakes located at shallow depth (< 4 km) and detached from the main cluster of deeper seismicity. We interpret the decrease in seismic velocities and the occurrence of such a shallow earthquake cluster as the effect of hydrothermal fluid released by the ascending magma batch and reaching the surface faster than the magma itself., (© 2023. The Author(s).)

The contrasting dynamics between Mauna Loa and Kīlauea have been studied over the last 100 years from multiple viewpoints. The fact that dynamic changes of one volcano trigger a dynamic response of the other volcano indicates a connection may exist. Petrological works show a direct relationship between the magmatic systems of these two volcanoes is not possible. We analysed DInSAR data and GPS measurements of ground deformation patterns associated with the activity of Mauna Loa and Kīlauea volcanoes. The DInSAR SBAS dataset spans the interval between 2003 and 2010, and was acquired along ascending and descending orbits of the ENVISAT (ESA) satellite under different look angles. Of the 10 tracks that cover the Big Island (Hawai'i), 4 cover both volcanic edifices. Using GPS measurements, we computed the areal strain on 15 triplets of stations for Kīlauea volcano and 11 for Mauna Loa volcano. DInSAR data was analysed by applying Independent Component Analysis (ICA) to decompose the time-varying ground deformation pattern of both volcanoes. The results revealed anticorrelated ground deformation behaviour of the main calderas of Mauna Loa and Kīlauea, meaning that the opposite response is seen in the ground deformation of one volcano with respect to the other. At the same time, Kīlauea exhibits a more complex pattern, with an additional component, which appears not to be correlated with the dynamics of Mauna Loa. The GPS areal strain time series support these findings. To corroborate and help interpret the results, we performed inverse modelling of the observed ground deformation pattern using analytical source models. The results indicate that the ground deformation of Mauna Loa is associated with a dike-shaped source located at 6.2 km depth. In comparison, the anticorrelated ground deformation of Kīlauea is associated with a volumetric source at 1.2 km depth. This excludes a hydraulic connection as a possible mechanism to explain the anticorrelated behaviour; instead, we postulate a stress-transfer mechanism. To support this hypothesis, we performed a 3D numerical modelling of stress and strain fields in the study area, determining the elastic interaction of each source over the others. The most relevant finding is that the Mauna Loa shallow plumbing system can affect the shallowest magmatic reservoir of Kīlauea, while the opposite scenario is unlikely. Conversely, the second independent component observed at Kīlauea is associated to a sill-shaped source located at a depth of 3.5 km, which is less affected by this interaction process., (© 2022. The Author(s).)

For the first time, we obtained high-resolution images of Earth's interior of the La Palma volcanic eruption that occurred in 2021 derived during the eruptive process. We present evidence of a rapid magmatic rise from the base of the oceanic crust under the island to produce an eruption that was active for 85 days. This eruption is interpreted as a very accelerated and energetic process. We used data from 11,349 earthquakes to perform travel-time seismic tomography. We present high-precision earthquake relocations and 3D distributions of P and S-wave velocities highlighting the geometry of magma sources. We identified three distinct structures: (1) a shallow localised region (< 3 km) of hydrothermal alteration; (2) spatially extensive, consolidated, oceanic crust extending to 10 km depth and; (3) a large sub-crustal magma-filled rock volume intrusion extending from 7 to 25 km depth. Our results suggest that this large magma reservoir feeds the La Palma eruption continuously. Prior to eruption onset, magma ascended from 10 km depth to the surface in less than 7 days. In the upper 3 km, melt migration is along the western contact between consolidated oceanic crust and altered hydrothermal material., (© 2022. The Author(s).)

The study of geothermal systems is nowadays a topic of great importance because of the huge amount of energy that could be converted in electricity for human consumption from such sources. Among the various geophysical methods employed to study geothermal reservoirs, the magnetotelluric (MT) method is capable to reveal the internal structures of the subsurface and interpret the geological structures from the electrical resistivity. We present the first 3D resistivity model of La Palma (Canary archipelago, Spain) obtained from a dataset of 44 broadband magnetotelluric soundings distributed around the island. Our results highlight the presence of resistivity anomalies, spatially coinciding with density anomalies present in literature. In the north of the island, a high resistivity anomaly can be interpreted as the signature of an old intrusive body beneath the Taburiente caldera. In the south, a complex resistivity structure around the Cumbre Vieja volcanic ridge could be indicative of presence of an active geothermal system. In particular, low-resistivity anomalies, located in a high-fractured zone, have values compatible with clay alteration caps (illite and illite-smectite). Such a result suggests the presence of hot rocks, or a dike system, heating fluids in the interior of Cumbre Vieja volcanic system.