Search

Fault zones are major sources of hazard for many populated regions around the world. Earthquakes still occur unanticipated, and research has started to observe fault properties with increasing spatial and temporal resolution, having the goal of detecting signs of stress accumulation and strength weakening that may anticipate the rupture. The common practice is monitoring source parameters retrieved from measurements; however, model dependence and strong uncertainty propagation hamper their usage for small and microearthquakes. Here, we decipher the ground motion (i.e., ground shaking) variability associated with microseismicity detected by dense seismic networks at a near-fault observatory in Irpinia, Southern Italy, and obtain an unprecedentedly sharp picture of the fault properties evolution both in time and space. We discuss the link between the ground-motion intensity and the source parameters of the considered microseismicity, showing a coherent spatial distribution of the ground-motion intensity with that of corner frequency, stress drop, and radiation efficiency. Our analysis reveals that the ground-motion intensity presents an annual cycle in agreement with independent geodetic displacement observations from two Global Navigation Satellite System stations in the area. The temporal and spatial analyses also reveal a heterogeneous behavior of adjacent fault segments in a high seismic risk Italian area. Concerning the temporal evolution of fault properties, we highlight that the fault segment where the 1980 Ms 6.9 Irpinia earthquake nucleated shows changes in the event-specific signature of ground-motion signals since 2013, suggesting changes in their frictional properties. This evidence, combined with complementary information on the earthquake frequency–magnitude distribution, reveals differences in fault segment response to tectonic loading, suggesting rupture scenarios of future moderate and large earthquakes for seismic hazard assessment.

We estimate the source parameters of small-magnitude earthquakes that occurred during 2008–2020 in the Irpinia faults area (southern Italy). We apply a spectral decomposition approach to isolate the source contribution from propagation and site effects for ∼3000 earthquakes in the local magnitude range between ML 0 and 4.2. We develop our analyses in three steps. First, we fit the Brune (1970) model to the nonparametric source spectra to estimate corner frequency and seismic moment, and we map the spatial distribution of stress drop across the Irpinia area. We found stress drops in the range 0.4–8.1 MPa, with earthquakes deeper than 7 km characterized by higher average stress drop (i.e., 3.2 MPa). Second, assuming a simple stress-release model (Kanamori and Heaton, 2000), we derive fracture energy and critical slip-weakening distance. The spatial variability of stress drop and fracture energy allows us to image the present stress conditions of fault segments activated during the 23 November 1980 Ms 6.9 earthquake. The variability of the source parameters shows clear patterns of the fault mechanical properties, suggesting that the Irpinia fault system can be divided into three main sectors, with the northern and southern ones showing different properties from the central one. Our results agree with previous studies indicating the presence of fluids with different composition in the different sectors of the Irpinia fault system. In the third step, we compare the time evolution of source parameters with a time series of geodetic displacement recorded near the fault system. Temporal trends in the correlation between geodetic displacement and different source parameters indicate that the poroelastic deformation perturbation generated by the karst aquifer recharge is modulating not only the occurrence rate of microseismicity (D’Agostino et al., 2018) but may lead to rupture asperities with different sizes and characteristics.

Estimating local site effects and modifying the uncertainty in ground motion predictions are two indispensable parts of partially non-ergodic site-specific PSHA. Local site effects can be estimated using site response simulations or recorded ground motions at the site. When such predictions are available, the aleatory variability of ground motions used in PSHA can be changed to the single station sigma value. However, in these cases, the epistemic uncertainty in predicting site effects must be incorporated into the hazard analyses. This research focuses on the challenges specific to conducting partially non-ergodic site-specific PSHA using recorded ground motions or site response analysis. The main challenge in estimating local site effects using recorded data is whether ground motions collected in a relatively short time can be used to estimate site effects for long return period events. We first develop a database for recorded ground motions at the KiK-net array to investigate this question and use this database to develop a predictive model for the Fourier Amplitude Spectra of ground motions. The ground motion model (GMM) residuals are used to investigate the stability of site terms across different tectonic regimes. We observe that empirical site terms are stable across different tectonic regimes. This observation allows the use of ground motions from any tectonic regime (whether they belong to the tectonic regime that controls the hazard or not) to estimate local site effects. Moreover, in Fourier amplitude, site effects are not dependent on event magnitude and source to site distance; therefore, estimates of site effects from low magnitude events can be easily extrapolated to larger events. The Fourier amplitude GMM developed in this study adds to the library of Fourier amplitude models to be used in future partially non-ergodic site-specific PSHAs. In practice, one of the most common tools for simulating wave propagation is 1-D site response analysis. Two central assumptions in 1-D site response analysis are that the soil profile is comprised of horizontal soil layers of infinite extent and that the vertically propagating SH-waves control the horizontal component of ground motion. SH-waves tend to propagate vertically near the surface because as earthquake waves hit softer layers traveling from the source to the site, they refract until the path becomes steeply inclined. The validity of both assumptions in 1-D site response depends on the geological setting at the site and the geology between the earthquake source and the site, raising the question of which sites are suitable for 1-D site response analysis and what the model error in 1-D site response analysis is. We use the GMM developed for FAS to estimate observed and empirical site terms. The empirical site effects are then compared with the theoretical site effects to determine whether sites are amenable to 1-D site response analyses, and to quantify the model error in the analyses. Doctor of Philosophy It is impossible to predict future earthquake-induced ground motions due to randomness in the process and a lack of knowledge. In fact, there are significant uncertainties not only in predicting the location, time, and magnitude of a future earthquake but also in predicting the intensity of ground motion induced by a given future earthquake. Therefore, assessing the safety of the human environment against earthquake hazards requires a method that considers all sources of uncertainties. To this end, Earthquake Engineers have developed Probabilistic Seismic Hazard Analysis(PSHA) framework. Structural engineers use the results of PSHA to design a new structure or assess the safety of an existing building. The accuracy of PSHA estimations leads to designs that are both safe and cost-efficient. The distribution of possible ground motions induced by a given earthquake scenario significantly controls the result of PSHA. This distribution should consider the effect of source, source to site path, and local site effects. This research focuses on improving PSHA results by estimating local site effects using recorded ground motions or simulating wave propagation in the site. In estimating local site effects using recorded data, the local site effect observed in ground motions collected in a relatively short time window is used to estimate hazards from all scenarios. However, the collected ground motions usually belong to frequent low magnitude events that are different from large magnitude events that control the hazard. This difference requires either using a measure of local site effect that is independent of the magnitude and distance of the earthquake or considering the effect of magnitude and distance on the local site effect estimate. Moreover, since frequent events sample different sources and paths than large events, we need to make sure the local site effect is consistent across different sources and paths. This research develops Ground Motion Models(GMMs) for Fourier amplitude, a linear function of ground motion times series, using Japanese ground motions. The ratio of Fourier amplitude at the surface over bedrock is a measure of local site effect that is not dependant on magnitude and distance. The model is then used to see if the trade-off between source and site effect and path and site effect is significant or not. In practice, one of the most common tools for simulating wave propagation is 1-D site response analysis. Two central assumptions in 1-D site response analysis are that the soil profile comprises horizontal soil layers of infinite extent and that the vertically propagating horizontal shear waves (SH-waves) control the horizontal component of ground motion. SH-waves tend to propagate vertically near the surface because as earthquake waves hit softer layers traveling from the source to the site, they refract until the path becomes vertically inclined. The validity of both assumptions in 1-D site response depends on the geological setting at the site and the geology between the earthquake source and the site, raising the question of which sites are suitable for 1-D site response analysis and what the model error in 1-D site response analysis is. We use the GMM developed for FAS to estimate empirical local site effects. The empirical site effects are then compared with the theoretical site effects to determine whether sites are amenable to 1-D site response analyses and quantify the model error in the analyses.

Ground‐motion prediction equations (GMPEs) are calibrated to predict the intensity of ground shaking at any given location, based on earthquake magnitude, source‐to‐site distance, local soil amplifications, and other parameters. GMPEs are generally assumed to be independent of time; however, evidence is increasing that large earthquakes modify the shallow soil conditions and those of the fault zone for months or years. These changes may affect the intensity of shaking and result in time‐dependent effects that can potentially be resolved by analyzing between‐event residuals (residuals between observed and predicted ground motion for individual earthquakes averaged over all stations). Here, we analyze a data set of about 65,000 recordings for about 1400 earthquakes in the moment magnitude range 2.5–6.5 that occurred in central Italy from 2008 to 2017 to capture the temporal variability of the ground shaking at high frequency. We first compute between‐event residuals for each earthquake in the Fourier domain with respect to a GMPE developed ad hoc for the analyzed data set. The between‐events show large changes after the occurrence of mainshocks such as the 2009 Mw 6.3 L'Aquila, the 2016 Mw 6.2 Amatrice, and Mw 6.5 Norcia earthquakes. Within the time span of a few months after the mainshocks, the between‐event contribution to the ground shaking varies by a factor 7. In particular, we find a large drop in the between‐events in the aftermath of the L'Aquila earthquake, followed by a slow positive trend that leads to a recovery interrupted by a new drop at the beginning of 2014. We also quantify the frequency‐dependent correlation between the Brune stress drop Δσ and the between‐events. We find that the temporal changes of Δσ resemble those of the between‐event residuals; in particular, during the period when the between‐events show the positive trend, the average logarithm of Δσ increases with an annual rate of 0.19 (i.e., the amplification factor for Δσ is 1.56 per year). Breakpoint analysis located a change in the linear trend coefficients of Δσ versus time in February 2014, although no large earthquakes occurred at that time. Finally, the temporal variability of Δσ mirrors the relative seismic‐velocity variations observed in previous studies for the same area and period, suggesting that both crack healing along the main fault system and healing of microcracks distributed at shallow depths throughout the surrounding region might be necessary to explain the wider observations of postearthquake recovery.

This report describes the main conclusions of the STREST project (and associated guidelines) to evaluate hazard of Low-Probability High-Consequences (LP-HC) events used to define stress tests for non-nuclear Critical Infrastructures (CIs). Several new approaches have been developed to assess these extreme hazard scenarios and to evaluate the associated uncertainties. This report presents a summary of the developments, results and products issued from Work Package 3 (WP3) of STREST. It is given as a set of “recommendations” for potential users responsible of the estimation of hazard for a particular non-nuclear CI in the European Union and other countries. The methods and guidelines are dedicated to two different target-users: project managers and hazard experts. It poses the main differences with a traditional Probabilistic Hazard Assessment analysis, the benefits and extra challenges, and the particular information requirements for the three selected infrastructure classes covered in STREST. In a simple and understandable manner, it summarizes the principal available tools, the main references and the application examples issued from the project in order to help the users in the realization of theirs studies., JRC.E.4-Safety and Security of Buildings