Skip to article content

La Palma Seismicity 2021

An analysis of earthquake swarms in relation to the 2021 eruption

Abstract

In September 2021, a significant jump in seismic activity on the island of La Palma (Canary Islands, Spain) signaled the start of a volcanic crisis that still continues at the time of writing. Earthquake data is continually collected and published by the Instituto Geográphico Nacional (IGN). We have created an accessible dataset from this and completed preliminary data analysis which shows seismicity originating at two distinct depths, consistent with the model of a two reservoir system feeding the currently very active volcano.

Key Points

You may specify 1 to 3 keypoints for this PDF template

These keypoints are complete sentences and less than or equal to 140 characters

They are specific to this PDF template, so they will not appear in other exports plain

Plain Language Summary

In September 2021, heightened seismic activity on La Palma signaled an ongoing volcanic crisis; analysis of Instituto Geográfico Nacional data reveals seismicity at two depths, supporting a dual-reservoir system for the active volcano.

Keywords:La PalmaEarthquakes

Introduction

The content of your article is written in MyST markdown and supports standard markdown typography and many directives and roles for figures, tables, equations, etc.

La Palma is one of the west most islands in the Volcanic Archipelago of the Canary Islands, a Spanish territory situated is the Atlantic Ocean where at their closest point are 100km from the African coast Figure 1 The island is one of the youngest, remains active and is still in the island forming stage.

Figures may be added to your article using the figure directive. They may refer to images saved in your images/ folder, images from the web, or notebook cell outputs referenced by label. The :name: is used to reference the figure in your text; a reference to the following figure is found in the paragraph above. The figure caption is given as the body of this directive.

Map of La Palma in the Canary Islands. Image credit NordNordWest

Figure 1:Map of La Palma in the Canary Islands. Image credit NordNordWest

La Palma has been constructed by various phases of volcanism, the most recent and currently active being the Cumbre Vieja volcano, a north-south volcanic ridge that constitutes the southern half of the island.

Eruption History

A number of eruptions were recorded since the colonization of the islands by Europeans in the late 1400s, these are summarized in Table 1.

Simple tables may be created using the list-table directive. Similar to figures, tables may be referenced in the text by their name. The caption for this table is the first line of the directive.

Table 1:Recent historic eruptions on La Palma

NameYear
Current2021
Teneguía1971
Nambroque1949
El Charco1712
Volcán San Antonio1677
Volcán San Martin1646
Tajuya near El Paso1585
Montaña Quemada1492

This equates to an eruption on average every 79 years up until the 1971 event. The probability of a future eruption can be modeled by a Poisson distribution (1).

Numbered equations may be defined using the math directive or in line. Equations defined with the math directive may be reference in the text by label.

p(x)=eλλxx!p(x)=\frac{e^{-\lambda} \lambda^{x}}{x !}

Where λ is the number of eruptions per year, λ=179\lambda=\frac{1}{79} in this case. The probability of a future eruption in the next tt years can be calculated by:

pe=1etλp_e = 1-\mathrm{e}^{-t \lambda}

So following the 1971 eruption the probability of an eruption in the following 50 years — the period ending this year — was 0.469. After the event, the number of eruptions per year moves to λ=175\lambda=\frac{1}{75} and the probability of a further eruption within the next 50 years (2022-2071) rises to 0.487 and in the next 100 years, this rises again to 0.736.

Magma Reservoirs

You may add citations two ways. First, you may simply insert a markdown link to a DOI like so: Thompson et al. (1994). No additional bibliographic information is required for this approach; the reference will be looked up by DOI and added implicitly to the references. Alternatively, you may provide the bibliography directly as references.bib BibTeX file, then embed the citation by BibTeX key in your text using the @cite2023 or [@cite2023; @cite2023b] for narrative or parenthetical citations, respectively. The following paragraph provides an example of this. A single paper may combine both DOI and BibTeX citations.

Studies of the magma systems feeding the volcano, such as Marrero et al. (2019) has proposed that there are two main magma reservoirs feeding the Cumbre Vieja volcano; one in the mantle (30-40km depth) which charges and in turn feeds a shallower crustal reservoir (10-20km depth).

Proposed model from .

Figure 2:Proposed model from Marrero et al. (2019).

In this paper, we look at recent seismicity data to see if we can see evidence of such a system action, see Figure 2.

Dataset

All data used in the notebook should be present in the data/ folder so notebooks may be executed in place with no additional input.

The earthquake dataset used in our analysis was generated from the IGN web portal this is public data released under a permissive license. Data recorded using the network of Seismic Monitoring Stations on the island. A web scraping script was developed to pull data into a machine-readable form for analysis. That code tool is available on GitHub along with a copy of recently updated data.

Results

The dataset was loaded into a Jupyter notebook visualization and filtered down to La Palma events only. This results in 5465 data points which we then visualized to understand their distributions spatially, by depth, by magnitude and in time.

<Figure size 7200x3600 with 3 Axes>

Figure 3:Earthquake data over time (n=5465) to understand their distributions spatially, by depth, by magnitude and in time. This figure uses cell output from the visualization notebook. The first line of the cell is #| label: eq-timeline. Referencing that label pulls in the output of the cell as a figure.

From our analysis in Figure 3, we can see 3 different systems in play.

Firstly, the shallow earthquake swarm leading up to the eruption on 19th September, related to significant surface deformation and shallow magma intrusion.

After the eruption, continuous shallow seismicity started at 10-15km corresponding to magma movement in the crustal reservoir.

Subsequently, high magnitude events begin occurring at 30-40km depths corresponding to changes in the mantle reservoir. These are also continuous but occur with a lower frequency than in the crustal reservoir.

Conclusions

From the analysis of the earthquake data collected and published by IGN for the period of 11 September through to 9 November 2021. Visualization of the earthquake events at different depths appears to confirm the presence of both mantle and crustal reservoirs as proposed by Marrero et al. (2019).

Acknowledgments

The authors would like to thank the Notebooks Now working groups for their help on creating these templates and advising on metadata.

Data Availability

A web scraping script was developed to pull data into a machine-readable form for analysis. That code tool is available on GitHub along with a copy of recently updated data.

References
  1. Thompson, J. D., Higgins, D. G., & Gibson, T. J. (1994). CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Research, 22(22), 4673–4680. 10.1093/nar/22.22.4673
  2. Marrero, J., García, A., Berrocoso, M., Llinares, Á., Rodríguez-Losada, A., & Ortiz, R. (2019). Strategies for the development of volcanic hazard maps in monogenetic volcanic fields: the example of La Palma (Canary Islands). Journal of Applied Volcanology, 8. 10.1186/s13617-019-0085-5
La Palma Seismicity 2021
Seismic Monitoring Stations