Temporal models for the occurrence of Etna eruptions and implications for hazard assessment

Sandri, Laura; Garcia, Alexander; Scollo, Simona; Mereu, Luigi; Prestifilippo, Michele

doi:10.5194/egusphere-2025-5875

Preprints

https://doi.org/10.5194/egusphere-2025-5875

Preprints

04 Dec 2025

| 04 Dec 2025

Status: this preprint is open for discussion and under review for Natural Hazards and Earth System Sciences (NHESS).

Temporal models for the occurrence of Etna eruptions and implications for hazard assessment

Laura Sandri, Alexander Garcia, Simona Scollo, Luigi Mereu, and Michele Prestifilippo

Abstract. Mt Etna volcanic activity is broadly divided into flank eruptions and summit paroxysms. Here, building on previously-available literature and data on the start time of these events, we collate two separate catalogs of the two activity types. Then we separately model their temporal occurrence. The catalog of flank eruptions, spanning the last 400 years, has been modelled by means of the most widely used renewal models, among which the best one (through Akaike Information Criterion) is the Brownian Passage Time. The catalog of summit paroxysms, covering the period 1986–2022, according to our cluster analysis is best characterized by 12 clusters of paroxysms. We separately analyze the inter-event times between onset times of successive clusters of paroxysms (inter-cluster inter-event times) and the inter-event times between successive paroxysms within clusters (intra-cluster inter-event times). Again, the Brownian Passage Time is the best-fitting model, obviously with very different parameters in the two cases. We test the best-fitting models by checking their ability to reproduce features of the real catalogs. Finally, we provide an example of how to use in practice such temporal models in the context of probabilistic hazard assessment, showing a possible use in the case of tephra fallout hazard from summit paroxysms.

Received: 27 Nov 2025 – Discussion started: 04 Dec 2025

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.

Download & links

Preprint (PDF, 1876 KB)

Supplement (2163 KB)

Download & links

Laura Sandri, Alexander Garcia, Simona Scollo, Luigi Mereu, and Michele Prestifilippo

Status: open (until 27 Apr 2026)

Post a comment Subscribe to comment alert

RC1:
'Comment on egusphere-2025-5875', Ian Main, 24 Feb 2026 reply
Temporal models for the occurrence of Etna eruptions and implications for hazard assessment
by Laura Sandri, Alexander Garcia, Simona Scollo, Luigi Mereu, and Michele Prestifilippo
This manuscript develops a new database for eruptions at Mount Etna, separates it into flank and summit eruptions, and analyses the data to determine the best fitting model in the two cases. There are many examples of good practice in the analysis, including the testing of multiple hypotheses, the use of an appropriate information criterion to choose the best one, and the simulation of model variability and averaging by running an ensemble of models, not just using the best fit. The results show a clear difference between the two types of eruption, with summit eruptions showing clear evidence of short-term clustering, even in the primary data, and the best fit model for flank eruptions consistent with long-term anti-clustering indicating a renewal process with some memory. The results could in principle be used in managing the future risk of eruptions of the different types in operational forecasting mode for a time varying hazard rate, and some examples of this are given to illustrate this.
The paper is of a suitable standard for publication, and I recommend publication after addressing the following mostly minor points.
The abstract could be more specific about the example application to hazard forecasting.

Section 2.1. Can you say something about the uncertainty in dating events in such a long catalogue, and how this varies in time?

Line 50. Can you also demonstrate and justify the statement of completeness and how this was determined? (In any case you change the completeness starting date later and explain why, so this is an initial guess, not the final version). In fact, the completeness is determined effectively by trial and error with stationarity as a criterion.

Line 58. How was the catalogue merging done? Was this simply by combining these databases, or did it involve some recalibration and/or decisions on which to use in the case of overlapping events? Would you say this resulted in a homogeneous catalogue after the merger, and why?

Line 70. ‘does not show evident clustering’. Where is the evidence for this statement? This assertion (in fact all unsupported assertions) should be justified by a logical narrative and cross-reference to the evidence.

Line 82. ‘due to their tendency to cluster’ - see previous comment. This assertion should be justified.

Line 93. Can you explain how the compound model works? Is this a sum of two different processes or is this a temporary change where one process switches off and another switches on?

Line 121. Can you compare these annual rates with what you would get from a standard Poisson model? What is the probability gain on the forecasts compered to this null hypothesis?

Lines 169-171. The evidence and chain of logic behind this inference on the hazard function is not presented, but it should be. Do you mean a steady increase, or a temporary increase and transient decay to the background level after the cluster stops? I think you should cross-refer to the relevant figure for the hazard functions here for both types of eruption from the two variants of the BPT model, and compare and contrast these with each other rather than simply assert the inference. It would also be useful to add a time-independent rate to the relevant figure and discussion for comparison.

Figure 1. It would be good to see the incremental probability function here as well as the cumulative version.

Figure 3. There is clear evidence of clustering here, but you don’t refer to this in the main text as evidence when you introduce the notion of clustering.

Figure 4. I think it would be useful to also show the variability or uncertainty in the best fit curves (red) from the simulations, similar to Figure 1.

Figure 8. Can you explain why the lower percentile curve is shorter than the upper one?

Ian Main, University of Edinburgh, 24 February 2026

Reply
Citation: https://doi.org/10.5194/egusphere-2025-5875-RC1

Laura Sandri, Alexander Garcia, Simona Scollo, Luigi Mereu, and Michele Prestifilippo

Supplement

https://doi.org/10.5194/egusphere-2025-5875-supplement

Data sets

OTETNA Laura Sandri et al. https://doi.org/10.13127/etna/otetna

Laura Sandri, Alexander Garcia, Simona Scollo, Luigi Mereu, and Michele Prestifilippo

Viewed

Total article views: 309 (including HTML, PDF, and XML)

HTML	PDF	XML	Total	Supplement	BibTeX	EndNote
176	111	22	309	48	24	23

HTML: 176
PDF: 111
XML: 22
Total: 309
Supplement: 48
BibTeX: 24
EndNote: 23

Views and downloads (calculated since 04 Dec 2025)

Month	HTML	PDF	XML	Total
Dec 2025	77	37	12	126
Jan 2026	38	28	4	70
Feb 2026	31	24	4	59
Mar 2026	22	20	2	44
Apr 2026	8	2	0	10

Cumulative views and downloads (calculated since 04 Dec 2025)

Month	HTML	PDF	XML	Total
Dec 2025	77	37	12	126
Jan 2026	38	28	4	70
Feb 2026	31	24	4	59
Mar 2026	22	20	2	44
Apr 2026	8	2	0	10

Viewed (geographical distribution)

Total article views: 300 (including HTML, PDF, and XML) Thereof 300 with geography defined and 0 with unknown origin.

Country	#	Views	%

Latest update: 09 Apr 2026

Short summary

We model how eruptive events occur in time at Mt Etna. We distinguish between eruptions from flank fissures and paroxysmal eruptions from the summit craters. Our model allows to compute the expected annual frequency of flank events and the expected number of summit paroxysms in a future period of time, with uncertainty estimates. We show how this model can improve the assessment of probabilistic hazard from these eruptions.


Total:	0
HTML:	0
PDF:	0
XML:	0