Options
Taroni, Matteo
Loading...
Preferred name
Taroni, Matteo
Email
matteo.taroni@ingv.it
Staff
staff
ORCID
49 results
Now showing 1 - 10 of 49
- PublicationRestrictedWhen time and faults matter: towards a time-dependent probabilistic SHA in Calabria, Italy(2017)
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; In this study, we attempt to improve the standards in Probabilistic Seismic Hazard Assessment (PSHA) towards a time-dependent hazard assessment by using the most advanced methods and new databases for the Calabria region, Italy. In this perspective we improve the knowledge of the seismotectonic framework of the Calabrian region using geologic, tectonic, paleoseismological, and macroseismic information available in the literature. We built up a PSHA model based on the long-term recurrence behavior of seismogenic faults, together with the spatial distribution of historical earthquakes. We derive the characteristic earthquake model for those sources capable of rupturing the entire fault segment (full-rupture) independently with a single event of maximum magnitude. We apply the floating rupture model to those earthquakes whose location is not known sufficiently constrained. We thus associate these events with longer fault systems, assuming that any such earthquake can rupture anywhere within the particular fault system (floating partial-rupture) with uniform probability. We use a Brownian Passage Time (BPT) model characterized by mean recurrence, aperiodicity, or uncertainty in the recurrence distribution and elapsed time since the last characteristic earthquake. The purpose of this BPT model is to express the time-dependence of the seismic processes to predict the future ground motions in the region. Besides, we consider the influence on the probability of earthquake occurrence controlled by the change in static Coulomb stress (ΔCFF) due to fault interaction; to pursue this, we adopt a model built on the fusion of BPT model (BPT + ΔCFF). We present our results for both time-dependent (renewal) and time-independent (Poisson) models in terms of Peak Ground Acceleration (PGA) maps for 10% probability of exceedance in 50 years. The hazard may increase by more than 20% or decrease by as much as 50% depending on the different occurrence model. Seismic hazard in terms of PGA decreases about 20% in the Messina Strait, where a recent major earthquake took place, with respect to traditional time-independent estimates. PGA near the city of Cosenza reaches ~ 0.36 g for the time-independent model and 0.40 g for the case of the time-dependent one (i.e. a 15% increase). Both the time-dependent and time-independent models for the period of 2015–2065 demonstrate that the city of Cosenza and surrounding areas bear the highest seismic hazard in Calabria.1588 18 - PublicationRestrictedEarthquake forecasting during the complex Amatrice-Norcia seismic sequence(2017-09)
; ; ; ; ; Earthquake forecasting is the ultimate challenge for seismologists, because it condenses the scientific knowledge about the earthquake occurrence process, and it is an essential component of any sound risk mitigation planning. It is commonly assumed that, in the short term, trustworthy earthquake forecasts are possible only for typical aftershock sequences, where the largest shock is followed by many smaller earthquakes that decay with time according to the Omori power law. We show that the current Italian operational earthquake forecasting system issued statistically reliable and skillful space-time-magnitude forecasts of the largest earthquakes during the complex 2016-2017 Amatrice-Norcia sequence, which is characterized by several bursts of seismicity and a significant deviation from the Omori law. This capability to deliver statistically reliable forecasts is an essential component of any program to assist public decision-makers and citizens in the challenging risk management of complex seismic sequences.210 5 - PublicationOpen AccessHow many strong earthquakes will there be tomorrow?In this note, we study the distribution of earthquake numbers in both worldwide and regional catalogs: in the Global Centroid Moment Tensor catalog, from 1980 to 2019 for magnitudes Mw 5. 5+ and 6.5+ in the first case, and in the Italian instrumental catalog from 1960 to 2021 for magnitudes Mw 4.0+ and 5.5+ in the second case. A subset of the global catalog is also used to study the Japanese region. We will focus our attention on short-term time windows of 1, 7, and 30 days, which have been poorly explored in previous studies. We model the earthquake numbers using two discrete probability distributions, i.e., Poisson and Negative Binomial. Using the classical chi-squared statistical test, we found that the Poisson distribution, widely used in seismological studies, is always rejected when tested against observations, while the Negative Binomial distribution cannot be disproved for magnitudes Mw 6.5+ in all time windows of the global catalog. However, if we consider the Japanese or the Italian regions, it cannot be proven that the Negative Binomial distribution performs better than the Poisson distribution using the chi-squared test. When instead we compared the performances of the two distributions using the Akaike Information Criterion, we found that the Negative Binomial distribution always performs better than the Poisson one. The results of this study suggest that the Negative Binomial distribution, largely ignored in seismological studies, should replace the Poisson distribution in modeling the number of earthquakes.
63 47 - PublicationOpen AccessWhat is the impact of the August 24, 2016 Amatrice earthquake on the seismic hazard assessment in central Italy?(2016)
; ; ; ; ; ; ; The recent Amatrice strong event (Mw6.0) occurred on August 24, 2016 in Central Apennines (Italy) in a seismic gap zone, motivated us to study and provide better understanding of the seismic hazard assessment in the macro area defined as “Central Italy”. The area affected by the sequence is placed between the Mw6.0 1997 Colfiorito sequence to the north (Umbria-Marche region) the Campotosto area hit by the 2009 L’Aquila sequence Mw6.3 (Abruzzo region) to the south. The Amatrice earthquake occurred while there was an ongoing effort to update the 2004 seismic hazard map (MPS04) for the Italian territory, requested in 2015 by the Italian Civil Protection Agency to the Center for Seismic Hazard (CPS) of the Istituto Nazionale di Geofisica e Vulcanologia INGV. Therefore, in this study we brought to our attention new earthquake source data and recently developed ground-motion prediction equations (GMPEs). Our aim was to validate whether the seismic hazard assessment in this area has changed with respect to 2004, year in which the MPS04 map was released. In order to understand the impact of the recent earthquakes on the seismic hazard assessment in central Italy we compared the annual seismic rates calculated using a smoothed seismicity approach over two different periods; the Parametric Catalog of the Historical Italian earthquakes (CPTI15) from 1871 to 2003 and the historical and instrumental catalogs from 1871 up to 31 August 2016. Results are presented also in terms of peak ground acceleration (PGA), using the recent ground-motion prediction equations (GMPEs) at Amatrice, interested by the 2016 sequence.814 96 - PublicationOpen AccessApplication of an ensemble earthquake rate model in Italy, considering seismic catalogs and fault moment release(2020)
; ; ; ; ; ; ; We develop an ensemble earthquake rate model that provides spatially variable time-independent (Poisson) long-term annual occurrence rates of seismic events throughout Italy, for magnitude bin of 0.1 units from Mw ≥4.5 in spatial cells of 0.1° x 0.1°. We weighed seismic activity rates of smoothed seismicity and fault-based inputs to build our earthquake rupture forecast model, merging it into a single ensemble model. Both inputs adopt a tapered Gutenberg-Richter relation with a single b-value and a single corner magnitude estimated by earthquakes catalog. The spatial smoothed seismicity was obtained using the classical kernel smoothing method with the inclusion of magnitude dependent completeness periods applied to the Historical (CPTI15) and Instrumental seismic catalogs. For each seismogenic source provided by the Database of the Individual Seismogenic Sources (DISS), we computed the annual rate of the events above Mw4.5, assuming that the seismic moments of the earthquakes generated by each fault are distributed according to the tapered Gutenberg-Richter relation with the same parameters of the smoothed seismicity models. Comparing seismic annual rates of the catalogs with those of the seismogenic sources, we realized that there is a good agreement between these rates in Central Apennines zones, whereas the seismogenic rates are higher than those of the catalogs in the north east and south of Italy. We also tested our model against the strong Italian earthquakes (Mw5.5+), in order to check if the total number (N-test) and the spatial distribution (S-test) of these events was compatible with our model, obtaining good results, i.e. high p-values in the test. The final model will be a branch of the new Italian seismic hazard map.259 112 - PublicationOpen AccessWhat Is the Effect of Seismic Swarms on Short-Term Seismic Hazard and Gutenberg-Richter b-Value Temporal Variation? Examples from Central Italy, October–November 2023A seismic hazard can be quantified by using probabilities. Modern seismic forecasting models (e.g., Operational Earthquake Forecasting systems) allow us to quantify the short-term variations in such probabilities. Indeed these probabilities change with time and space, in particular after strong seismic events. However, the short-term seismic hazard could also change during seismic swarms, i.e., a sequence with several small-/medium-sized events. The goal of this work is to quantify these changes, using the Italian Operational Earthquake Forecasting system, and also estimate the variations in the Gutenberg–Richter b-value. We focus our attention on three seismic swarms that occurred in Central Italy in October–November 2023. Our results indicate that short-term variations in seismic hazard are limited, less than an order of magnitude, and also that b-value variations are not significant. Placing our findings in a more general context, we can state that according to currently available models and catalogs, the occurrence of seismic swarms does not significantly affect the short-term seismic hazard.
53 7 - PublicationOpen AccessEditorial: Physical and statistical approaches to earthquake modeling and forecasting(2023)
; ; ; ; ; ; ; earthquake forecasting and testing, ETAS model, operational earthquake forecasting, statistical seismology, short-term seismic hazard assessmen22 2 - PublicationOpen AccessGood practices in PSHA: declustering, b-value estimation, foreshocks and aftershocks inclusion; a case study in ItalyThe classical procedure of the Probabilistic Seismic Hazard Analysis (PSHA) requires a Poissonian distribution of earthquakes. Seismic catalogs follow a Poisson distribution just after the application of a declustering algorithm that leaves only one earthquake for each seismic sequence (usually the stronger, i.e. the main shock). Removing earthquakes from the seismic catalogs leads to underestimation of the annual rates of the events and consequently associate with low seismic hazard as indicated by several studies. In this study, we aim investigating the performance of two declustering methods on the Italian instrumental catalog and the impact of declustering on estimation of the b-value and on the seismic hazard analysis. To this end, first the spatial variation in the seismicity rate was estimated from the declustered catalogs using the adaptive smoothed seismicity approach, considering small earthquakes (Mw≥3.0). We then corrected the seismicity rates using new approach that allows for counting all events in the complete seismic catalog by simply changing the magnitude frequency distribution. The impact of declustering on seismic hazard analysis is illustrated using PSHA maps in terms of peak ground acceleration (PGA) and spectral acceleration (SA) in 2 s, with 10% and 2% probability of exceedance in 50 years, for Italy. We observed that the hazard calculated from the declustered catalogs was always lower than the hazard computed using the complete catalog. These results are in agreement with previous results obtained in different parts of the world.
570 63 - PublicationOpen AccessTemporal Variations of Seismicity Rates and Gutenberg–Richterb-Values for a Stochastic Declustered Catalog: An Example in Central Italy(2023)
; ; ; ; ; ; ; ;; ; ; ; ; One important aspect of the seismicity is the spatiotemporal clustering; hence, the distinction between independent and triggered events is a critical part of the analysis of seismic catalogs. Stochastic declustering of seismicity allows a probabilistic distinction between these two kinds of events. Such an approach, usually performed with the epidemic‐type aftershock sequence (ETAS) model, avoids the bias in the estimation of the frequency–magnitude distribution parameters if we consider a subset of the catalog, that is, only the independent or the triggered events. In this article, we present a framework to properly include the probabilities of any event to be independent (or triggered) both in the temporal variation of the seismic rates and in the estimation of the b‐value of the Gutenberg–Richter law. This framework is then applied to a high‐definition seismic catalog in the central part of Italy covering the period from April 2010 to December 2015. The results of our analysis show that the seismic activity from the beginning of the catalog to March 2013 is characterized by a low degree of clustering and a relatively high b‐value, whereas the following period exhibits a higher degree of clustering and a smaller b‐value.52 73 - PublicationRestrictedEnsemble Smoothed Seismicity Models for the New Italian Probabilistic Seismic Hazard MapWe develop a long-term (a few decades or longer) earthquake rate forecast for Italy based on smoothed seismicity for incor- poration in the 2017–2018 Italian Probabilistic Seismic Haz- ard Maps (IPSHM). Because the earthquake rate models from previous IPSHM were computed using source zones that were drawn around seismicity and tectonic provinces, the present model will be the first introduction of the smoothed seismic- ity method into the IPSHM. Smoothed seismicity models are constructed from both historical CPTI15 (Catalogo Parame- trico dei Terremoti Italiani, 1000–2014) and instrumental (1981–2016) earthquake catalogs and use both fixed and adaptive smoothing methods. We compute spatial likelihood values comparing the spatial distribution of observed earth- quakes with a suite of trial earthquake rate models to optimize smoothing parameters and catalogs. Then we produce an en- semble model using two different smoothing models (adap- tive and fixed) and two earthquake catalogs (historical and instrumental), which are weighted equally through a logic- tree approach to improve the forecast capability. We also compare our optimized smoothed seismicity models with the best two models of the Italian Collaboratory for the Study of Earthquake Predictability (CSEP) experiment and retro- spectively test them with the CSEP methodology. We ob- served that the ensemble model performs slightly better than the optimized fixed and the adaptive smoothing seismic- ity models obtained in this study and the best time-indepen- dent model of the CSEP Italian experiment. The preferred ensemble model forecasts an annual rate of 1.47 M ≥ 5:0 earthquakes, with higher rates mainly concentrating along the Apennines chain, eastern Alps, Calabria, and northeast Sicily. Finally, six ensemble models are created from the differ- ent smoothing methods using different weights through a logic-tree approach to estimate the uncertainty associated with the model.
921 6