We describe the space-time pattern of seismicity occurring on Campi Flegrei Caldera (CFC), Italy, where ground deformations and seismicity represent the drivers of its current bradyseismic crisis, well known and extensively studied at an international level. In detail we consider the seismicity in the time interval starting on 1.1.2000 and ending on 31.12.2023. We revise the statistics of the earthquake occurrence, focusing at possible precursory time changes of the b-parameter of the Gutenberg and Richter (G&R b−value) distribution and at the time distribution of the total seismic moment inside any swarm. To estimate the G&R b−value we use a Monte Carlo method instead of the ordinary Least Squares or Maximum likelihood methods, to easily measure the uncertainty on the b−value taking into account uncertainties on the magnitude estimates. Results show that G&R a−value and b−value calculated for cumulative and discrete distributions of Mw, the moment-magnitude, and Md, the so-called duration-magnitude, are the same inside the uncertainties; a−value and b−value for Md are significantly different from the same parameters estimated for Mw, being b−value for Mw close to the value of 1.0 and b−value for Md close to 0.8. The “bounded” G&R distribution fits the data yielding a−value and b−value close to those for the unbounded distribution. The mean annual rate of exceedance, calculated for the entire catalogue, results to be 0.033±0.015 (years−1) corresponding to a return period of 30±14 years for Mw=4.5. The time dependence of G&R b-parameter show a b−value time pattern characterized by variations slightly outside 1-σ uncertainty bar, tending to the value of 1 approaching present. As evidenced by several past studies, earthquakes in CFC occur in space-time clusters, with a mean time duration of 1 day. We selected the swarms with a selection algorithm, based on the joint estimation of inter-arrival times and inter-event distance for the consecutive event couples. The plot of the event number in each cluster vs the time of occurrence, clearly show that in CFC the number of cluster occurrences and the event number in each cluster accelerates during time starting from 2010. The time-pattern of the total seismic moment in each cluster shows that, contrarily to the event number, there is no evident striking increase of the total swarm moment as a function of time.We also consider distribution of the Energy Space Density of the CFC earthquakes, ESD. This quantity shows a clearly visible enlargement of the fractured rock volumes in the last 15 months, toward West, at a depth around 3000 m below sea level. The most fractured zone coincides with the greatest contrasts in seismic attenuation.The present study confirms that the current unrest phase is still ongoing, with an enlargement of the rock fracturing zone which extends southward and westward as compared with that measured 22 months before.
Read full abstract