The 2021 volcanic eruption at La Palma, Canary Islands, was the island’s most voluminous historical eruption. Little is known about this volcano’s feeding system. During the eruption, seismicity was distributed in two clusters at ~10-14 km and ~33-39 km depth, separated by an aseismic zone. This gap coincides with the location of weak seismic swarms in 2017-2021 and where petrological data have implied pre-eruptive magma storage. Here we use seismological methods to understand the seismic response to magma transfer, with 8,488 hypocentral relocations resolving small-scale seismogenic structures, and 156 moment tensors identifying stress heterogeneities and principal axes flips. Results suggest a long-lasting preparatory stage with the progressive destabilisation of an intermediate, mushy reservoir, and a co-eruptive stage with seismicity controlled by the drainage and interplay of two localised reservoirs. Our study provides new insights into the plumbing system that will improve the monitoring of future eruptions in the island.

Summary Seismicity along subduction interfaces is usually dominated by large mainshock-aftershocks sequences indicative of a continuum distribution of highly coupled large asperities. In the past decades, however, the increased resolution of seismic catalogues at some subduction zone seems to indicate instead a more complex rheological segmentation of the interface. Large and megathrust earthquake ruptures seem interspersed among regions of low seismic coupling and less stress build-up. In this weaker zones, the strain is primarily released via a combination of moderate size swarm-like seismicity and aseismic slip. Along the Chilean subduction zone, the densification of the seismic network allowed for the identification of localized seismic clusters, some of them appearing in the form of swarms before megathrust earthquakes. The origin and driving processes of this seismic activity have not yet been identified. In this study, we follow a systematic approach to characterize the seismicity at two persistent clusters in Central Chile, one located offshore Navidad and one inland, at ∼40 km depth beneath Vichuquén, which occurred throughout ∼20 years. We investigated these clusters, by deriving high resolution hypocentral locations and moment tensors and performing a detailed analysis of spatio-temporal patterns, magnitude, and inter-event time distributions of the clustered earthquakes. Both clusters are characterized by weak to moderate seismicity (below MW 6) and stand out as clear seismicity rate and Benioff strain anomalies. At the Navidad cluster, seismicity occurs in the form of swarms, with a characteristic duration of 2–7 days and location and thrust mechanisms compatible with activity on the slab interface. Conversely, we find at Vichuquén activity dominated by thrust earthquakes occurring as repeaters on the slab interface, with a slip rate of approximately ∼5.0 cm/yr. We attribute these clusters to local features of the subducting plate: the Navidad swarms are likely driven by repeated high pore pressure transients along a pre-fractured patch of the slab, while the seismicity at the Vichuquén cluster is interpreted as the result of a subducting seamount. Both clusters have been active before and after the MW 8.8 Maule earthquake and persisted afterwards with the seismicity decay following the Omori law. These interactions are especially evident for the Vichuquén cluster, where the seismicity rate increased considerably after the Maule earthquake and continues to be an area of clearly elevated seismicity rate compared to its surroundings.

No abstract