Refine
Has Fulltext
- no (69) (remove)
Year of publication
Language
- English (69) (remove)
Is part of the Bibliography
- yes (69)
Keywords
- site effects (4)
- KiK-net (3)
- earthquake (3)
- Aleatory variability (2)
- Earthquake hazards (2)
- Epistemic uncertainty (2)
- Europe (2)
- Ground-motion prediction equations (2)
- HVSR (2)
- Site effects (2)
In this study, we analyzed 10 yrs of seismicity in central Italy from 2008 to 2017, a period witnessing more than 1400 earthquakes in the magnitude range 2.5≤Mw≤6.5. The data set includes the main sequences that have occurred in the area, including those associated with the 2009 Mw 6.3 L'Aquila earthquake and the 2016–2017 sequence (Mw 6.2 Amatrice, Mw 6.1 Visso, and Mw 6.5 Norcia earthquakes). We calibrated a local magnitude scale, investigating the impact of changing the reference distance at which the nonparametric attenuation is tied to the zero‐magnitude attenuation function for southern California. We also developed an attenuation model to compute the radiated seismic energy (Es) from the time integral of the squared ground‐motion velocity. Seismic moment (M0) and stress drop (Δσ) were estimated for each earthquake by fitting a ω‐square model to the source spectra obtained by applying a nonparametric spectral inversion. The Δσ‐values vary over three orders of magnitude from about 0.1 to 10 MPa, the larger values associated with the mainshocks. The Δσ‐values describe a lognormal distribution with mean and standard deviation equal to log(Δσ)=(−0.25±0.45) (i.e., the mean Δσ is 0.57 MPa, with a 95% confidence interval from 0.08 to 4.79 MPa). The Δσ variability introduces a spread in the distribution of seismic energy versus moment, with differences in energy up two orders of magnitudes for earthquakes with the same moment. The variability in the high‐frequency spectral levels is captured by the local magnitude (ML), which scales with radiated energy as ML=(−1.59+0.52logEs) for logEs≤10.26 and ML=(−1.38+0.50logEs) otherwise. As the peak ground velocity increases with increasing Δσ, local and energy magnitudes perform better than moment magnitude as predictors for the shaking potential. The availability of different magnitude scales and source parameters for a large earthquake population will help characterize the between‐event ground‐motion variability in central Italy.
We present the results of a consistency check performed over the flatfile extracted from the engineering strong motion (ESM) database. The flatfile includes 23,014 recordings from 2179 earthquakes in the magnitude range from 3.5 to 7.8 that occurred since the 1970s in Europe and Middle East, as presented in the companion article by Lanzano et al. (Bull Earthq Eng, 2018a). The consistency check is developed by analyzing different residual distributions obtained from ad-hoc ground motion prediction equations for the absolute spectral acceleration (SA), displacement and Fourier amplitude spectra (FAS). Only recordings from earthquakes shallower than 40 km are considered in the analysis. The between-event, between-station and event-and-station corrected residuals are computed by applying a mixed-effect regression. We identified those earthquakes, stations, and recordings showing the largest deviations from the GMPE median predictions, and also evaluated the statistical uncertainty on the median model to get insights on the applicable magnitude–distance ranges and the usable period (or frequency) range. We observed that robust median predictions are obtained up to 8 s for SA and up to 20 Hz for FAS, although median predictions for Mw ≥ 7 show significantly larger uncertainties with ‘bumps’ starting above 5 s for SA and below 0.3 Hz for FAS. The between-station variance dominates over the other residual variances, and the dependence of the between-station residuals on logarithm of Vs30 is well-described by a piece-wise linear function with period-dependent slopes and hinge velocity around 580 m/s. Finally, we compared the between-event residuals obtained by considering two different sources of moment magnitude. The results show that, at long periods, the between-event terms from the two regressions have a weak correlation and the overall between-event variability is dissimilar, highlighting the importance of magnitude source in the regression results.
Applying conservation of energy to estimate earthquake frequencies from strain rates and stresses
(2020)
Estimating earthquake occurrence rates from the accumulation rate of seismic moment is an established tool of seismic hazard analysis. We propose an alternative, fault-agnostic approach based on the conservation of energy: the Energy-Conserving Seismicity Framework (ENCOS). Working in energy space has the advantage that the radiated energy is a better predictor of the damage potential of earthquake waves than the seismic moment release. In a region, ENCOS balances the stationary power available to cause earthquakes with the long-term seismic energy release represented by the energy-frequency distribution's first moment. Accumulation and release are connected through the average seismic efficiency, by which we mean the fraction of released energy that is converted into seismic waves. Besides measuring earthquakes in energy, ENCOS differs from moment balance essentially in that the energy accumulation rate depends on the total stress in addition to the strain rate tensor. To validate ENCOS, we exemplarily model the energy-frequency distribution around Southern California. We estimate the energy accumulation rate due to tectonic loading assuming poroelasticity and hydrostasis. Using data from the World Stress Map and assuming the frictional limit to estimate the stress tensor, we obtain a power of 0.8 GW. The uncertainty range, 0.3-2.0GW, originates mainly from the thickness of the seismogenic crust, the friction coefficient on preexisting faults, and models of Global Positioning System (GPS) derived strain rates. Based on a Gutenberg-Richter magnitude-frequency distribution, this power can be distributed over a range of energies consistent with historical earthquake rates and reasonable bounds on the seismic efficiency.
The simulation of broad-band (0.1 to 10 + Hz) ground-shaking over deep and spatially extended sedimentary basins at regional scales is challenging. We evaluate the ground-shaking of a potential M 6.5 earthquake in the southern Lower Rhine Embayment, one of the most important areas of earthquake recurrence north of the Alps, close to the city of Cologne in Germany. In a first step, information from geological investigations, seismic experiments and boreholes is combined for deriving a harmonized 3D velocity and attenuation model of the sedimentary layers. Three alternative approaches are then applied and compared to evaluate the impact of the sedimentary cover on ground-motion amplification. The first approach builds on existing response spectra ground-motion models whose amplification factors empirically take into account the influence of the sedimentary layers through a standard parameterization. In the second approach, site-specific 1D amplification functions are computed from the 3D basin model. Using a random vibration theory approach, we adjust the empirical response spectra predicted for soft rock conditions by local site amplification factors: amplifications and associated ground-motions are predicted both in the Fourier and in the response spectra domain. In the third approach, hybrid physics-based ground-motion simulations are used to predict time histories for soft rock conditions which are subsequently modified using the 1D site-specific amplification functions computed in method 2. For large distances and at short periods, the differences between the three approaches become less notable due to the significant attenuation of the sedimentary layers. At intermediate and long periods, generic empirical ground-motion models provide lower levels of amplification from sedimentary soils compared to methods taking into account site-specific 1D amplification functions. In the near-source region, hybrid physics-based ground-motions models illustrate the potentially large variability of ground-motion due to finite source effects.
We analyze the spatiotemporal evolution of seismicity during a sequence of moderate (an M-w 4.7 foreshock and M-w 5.8 mainshock) earthquakes occurring in September 2019 at the transition between a creeping and a locked segment of the North Anatolian fault in the central Sea of Marmara, northwest Turkey. To investigate in detail the seismicity evolution, we apply a matched-filter technique to continuous waveforms, thus reducing the magnitude threshold for detection. Sequences of foreshocks preceding the two largest events are clearly seen, exhibiting two different behaviors: a long-term activation of the seismicity along the entire fault segment and a short-term concentration around the epicenters of the large events. We suggest a two-scale preparation phase, with aseismic slip preparing the mainshock final rupture a few days before, and a cascade mechanism leading to the nucleation of the mainshock. Thus, our study shows a combination of seismic and aseismic slip during the foreshock sequence changing the strength of the fault, bringing it closer to failure.
1-D site response analysis dominates earthquake engineering practice, while local 2-D/3-D models are often required at sites where the site response is complex. For such sites, the 1-D representation of the soil column can account neither for topographic effects or dipping layers nor for locally generated horizontally propagating surface waves. It then remains a crucial task to identify whether the site response can be modelled sufficiently precisely by 1-D analysis. In this study we develop a method to classify sites according to their 1-D or 2-D/3-D nature. This classification scheme is based on the analysis of surface earthquake recordings and the evaluation of the variability and similarity of the horizontal Fourier spectra. The taxonomy is focused on capturing significant directional dependencies and interevent variabilities indicating a more probable 2-D/3-D structure around the site causing the ground motion to be more variable. While no significant correlation of the 1-D/3-D site index with environmental parameters and site proxies seems to exist, a reduction in the within-site (single-station) variability is found. The reduction is largest (up to 20 per cent) for purely 1-D sites. Although the taxonomy system is developed using surface stations of the KiK-net network in Japan as considerable additional information is available, it can also be applied to any (non-downhole array) site.
In this article, we address the question of how observed ground-motion data can most effectively be modeled for engineering seismological purposes. Toward this goal, we use a data-driven method, based on a deep-learning autoencoder with a variable number of nodes in the bottleneck layer, to determine how many parameters are needed to reconstruct synthetic and observed ground-motion data in terms of their median values and scatter. The reconstruction error as a function of the number of nodes in the bottleneck is used as an indicator of the underlying dimensionality of ground-motion data, that is, the minimum number of predictor variables needed in a ground-motion model. Two synthetic and one observed datasets are studied to prove the performance of the proposed method. We find that mapping ground-motion data to a 2D manifold primarily captures magnitude and distance information and is suited for an approximate data reconstruction. The data reconstruction improves with an increasing number of bottleneck nodes of up to three and four, but it saturates if more nodes are added to the bottleneck.
The main Marmara fault (MMF) extends for 150 km through the Sea of Marmara and forms the only portion of the North Anatolian fault zone that has not ruptured in a large event (Mw >7) for the last 250 yr. Accordingly, this portion is potentially a major source contributing to the seismic hazard of the Istanbul region. On 26 September 2019, a sequence of moderate-sized events started along the MMF only 20 km south of Istanbul and were widely felt by the population. The largest three events, 26 September Mw 5.8 (10:59 UTC), 26 September 2019 Mw 4.1 (11:26 UTC), and 20 January 2020 Mw 4.7 were recorded by numerous strong-motion seismic stations and the resulting ground motions were compared to the predicted means resulting from a set of the most recent ground-motion prediction equations (GMPEs). The estimated residuals were used to investigate the spatial variation of ground motion across the Marmara region. Our results show a strong azimuthal trend in ground-motion residuals, which might indicate systematically repeating directivity effects toward the eastern Marmara region.
The basic seismic load parameters for the upcoming national design regulation for DIN EN 1998-1/NA result from the reassessment of the seismic hazard supported by the German Institution for Civil Engineering (DIBt). This 2016 version of the national seismic hazard assessment for Germany is based on a comprehensive involvement of all accessible uncertainties in models and parameters and includes the provision of a rational framework for integrating ranges of epistemic uncertainties and aleatory variabilities in a comprehensive and transparent way. The developed seismic hazard model incorporates significant improvements over previous versions. It is based on updated and extended databases, it includes robust methods to evolve sets of models representing epistemic uncertainties, and a selection of the latest generation of ground motion prediction equations. The new earthquake model is presented here, which consists of a logic tree with 4040 end branches and essential innovations employed for a realistic approach. The output specifications were designed according to the user oriented needs as suggested by two review teams supervising the entire project. Seismic load parameters, for rock conditions of nu(S30) = 800 m/s, are calculated for three hazard levels (10, 5 and 2% probability of occurrence or exceedance within 50 years) and delivered in the form of uniform hazard spectra, within the spectral period range 0.02-3 s, and seismic hazard maps for peak ground acceleration, spectral response accelerations and for macroseismic intensities. Results are supplied as the mean, the median and the 84th percentile. A broad analysis of resulting uncertainties of calculated seismic load parameters is included. The stability of the hazard maps with respect to previous versions and the cross-border comparison is emphasized.
One paragraph of the manuscript of the paper has been inadvertently omitted in the very final stage of its compilation due to a technical mistake. Since this paragraph discusses the declustering of the used earthquake catalogue and is therefore necessary for the understanding of the seismicity data preprocessing, the authors decided to provide this paragraph in form of a correction. The respective paragraph belongs to chapter 2 of the paper, where it was placed originally, and should be inserted into the published paper before the second to the last paragraph. The omitted text reads as follows:
The Seismic Hazard Inferred from Tectonics based on the Global Strain Rate Map (SHIFT_GSRM) earthquake forecast was designed to provide high-resolution estimates of global shallow seismicity to be used in seismic hazard assessment. This model combines geodetic strain rates with global earthquake parameters to characterize long-term rates of seismic moment and earthquake activity. Although SHIFT_GSRM properly computes seismicity rates in seismically active continental regions, it underestimates earthquake rates in subduction zones by an average factor of approximately 3. We present a complementary method to SHIFT_GSRM to more accurately forecast earthquake rates in 37 subduction segments, based on the conservation of moment principle and the use of regional interface seismicity parameters, such as subduction dip angles, corner magnitudes, and coupled seismogenic thicknesses. In seven progressive steps, we find that SHIFT_GSRM earthquake-rate underpredictions are mainly due to the utilization of a global probability function of seismic moment release that poorly captures the great variability among subduction megathrust interfaces. Retrospective test results show that the forecast is consistent with the observations during the 1 January 1977 to 31 December 2014 period. Moreover, successful pseudoprospective evaluations for the 1 January 2015 to 31 December 2018 period demonstrate the power of the regionalized earthquake model to properly estimate subduction-zone seismicity.
A proper assessment of seismic reference site conditions has important applications as they represent the basis on which ground motions and amplifications are generally computed. Besides accounting for the average S-wave velocity over the uppermost 30 m (V-S30), the parameterization of high-frequency ground motions beyond source-corner frequency received significant attention. kappa, an empirical parameter introduced by Anderson and Hough (1984), is often used to represent the spectral decay of the acceleration spectrum at high frequencies. The lack of hard-rock records and the poor understanding of the physics of kappa introduced significant epistemic uncertainty in the final seismic hazard of recent projects. Thus, determining precise and accurate regional hard-rock kappa(0) values is critical. We propose an alternative procedure for capturing the reference kappa(0) on regional scales by linking thewell-known high-frequency attenuation parameter kappa and the properties of multiple-scattered coda waves. Using geological and geophysical data around more than 1300 stations for separating reference and soft soil sites and based on more than 10,000 crustal earthquake recordings, we observe that kappa(0) from multiple-scattered coda waves seems to be independent of the soil type but correlated with the hard-rock kappa(0), showing significant regional variations across Europe. The values range between 0.004 s for northern Europe and 0.020 s for the southern and southeastern parts. On the other hand, measuring kappa (and correspondingly kappa(0)) on the S-wave window (as classically proposed), the results are strongly affected by transmitted (reflected, refracted, and scattered) waves included in the analyzed window biasing the proper assessment of kappa(0). This effect is more pronounced for soft soil sites. In this way, kappa(coda)(0) can serve as a proxy for the regional hard-rock kappa(0) at the reference sites.
Ground-motion prediction equations (GMPE) are essential in probabilistic seismic hazard studies for estimating the ground motions generated by the seismic sources. In low-seismicity regions, only weak motions are available during the lifetime of accelerometric networks, and the equations selected for the probabilistic studies are usually models established from foreign data. Although most GMPEs have been developed for magnitudes 5 and above, the minimum magnitude often used in probabilistic studies in low-seismicity regions is smaller. Disaggregations have shown that, at return periods of engineering interest, magnitudes less than 5 may be contributing to the hazard. This paper presents the testing of several GMPEs selected in current international and national probabilistic projects against weak motions recorded in France (191 recordings with source-site distances up to 300 km, 3:8 <= M-w <= 4:5). The method is based on the log-likelihood value proposed by Scherbaum et al. (2009). The best-fitting models (approximately 2:5 <= LLH <= 3:5) over the whole frequency range are the Cauzzi and Faccioli (2008), Akkar and Bommer (2010), and Abrahamson and Silva (2008) models. No significant regional variation of ground motions is highlighted, and the magnitude scaling could be the predominant factor in the control of ground-motion amplitudes. Furthermore, we take advantage of a rich Japanese dataset to run tests on randomly selected low-magnitude subsets, and confirm that a dataset of similar to 190 observations, the same size as the French dataset, is large enough to obtain stable LLH estimates. Additionally we perform the tests against larger magnitudes (5-7) from the Japanese dataset. The ranking of models is partially modified, indicating a magnitude scaling effect for some of the models, and showing that extrapolating testing results obtained from low-magnitude ranges to higher magnitude ranges is not straightforward.
This study aims to identify the best-performing site characterization proxy alternative and complementary to the conventional 30 m average shear-wave velocity V-S30, as well as the optimal combination of proxies in characterizing linear site response. Investigated proxies include T-0 (site fundamental period obtained from earthquake horizontal-to-vertical spectral ratios), V-Sz (measured average shear-wave velocities to depth z, z = 5, 10, 20 and 30 m), Z(0.8) and Z(1.0) (measured site depths to layers having shear-wave velocity 0.8 and 1.0 km/s, respectively), as well as Z(x-infer) (inferred site depths from a regional velocity model, x = 0.8 and 1.0, 1.5 and 2.5 km/s). To evaluate the performance of a site proxy or a combination, a total of 1840 surface-borehole recordings is selected from KiK-net database. Site amplifications are derived using surface-to-borehole response-, Fourier- and cross-spectral ratio techniques and then are compared across approaches. Next, the efficacies of 7 single-proxies and 11 proxy-pairs are quantified based on the site-to-site standard deviation of amplification residuals of observation about prediction using the proxy or the pair. Our results show that T-0 is the best-performing single-proxy among T-0, Z(0.8), Z(1.0) and V-Sz. Meanwhile, T-0 is also the best-performing proxy among T-0, Z(0.8), Z(1.0) and Z(x-infer) complementary to V-S30 in accounting for the residual amplification after V-S30-correction. Besides, T-0 alone can capture most of the site effects and should be utilized as the primary site indicator. Though (T-0, V-S30) is the best-performing proxy pair among (V-S30, T-0), (V-S30, Z(0.8)), (V-S30, Z(1.0)), (V-S30, Z(x-infer)) and (T-0, V-Sz), it is only slightly better than (T-0, V-S20). Considering both efficacy and engineering utility, the combination of T-0 (primary) and V-S20 (secondary) is recommended. Further study is needed to test the performances of various proxies on sites in deep sedimentary basins.
Evaluation of a novel application of earthquake HVSR in site-specific amplification estimation
(2020)
Ground response analyses (GRA) model the vertical propagations of SH waves through flat-layered media (1DSH) and are widely carried out to evaluate local site effects in practice. Horizontal-to-vertical spectral ratio (HVSR) technique is a cost-effective approach to extract certain site-specific information, e.g., site fundamental frequency (f(0)), but HVSR values cannot be directly used to approximate the levels of S-wave amplifications. Motivated by the work of Kawase et al. (2019), we propose a procedure to correct earthquake HVSR amplitudes for direct amplification estimations. The empirical correction compensates HVSR by generic vertical amplification spectra categorized by the vertical fundamental frequency (f(0v)) via kappa-means clustering. In this investigation, we evaluate the effectiveness of the corrected HVSR in approximating observed linear amplifications in comparison with 1DSH modellings. We select a total of 90 KiK-net (Kiban Kyoshin network) surface-downhole sites which are found to have no velocity contrasts below their boreholes and thus of which surface-to-borehole spectral ratios (SBSRs) can be taken as their empirical transfer functions (ETFs). 1DSH-based theoretical transfer functions (TTFs) are computed in the linear domain considering uncertainties in Vs profiles through randomizations. Five goodness-of-fit metrics are adopted to gauge the closeness between observed (ETF) and predicted (i.e., TTF and corrected HVSR) amplifications in both amplitude and spectral shape over frequencies from f(0) to 25 Hz. We find that the empirical correction to HVSR is highly effective and achieves a "good match" in both spectral shape and amplitude at the majority of the 90 KiK-net sites, as opposed to less than one-third for the 1DSH modelling. In addition, the empirical correction does not require a velocity model, which GRAs require, and thus has great potentials in seismic hazard assessments.
The aim of this paper is to investigate the ability of various site-condition proxies (SCPs) to reduce ground-motion aleatory variability and evaluate how SCPs capture nonlinearity site effects. The SCPs used here are time-averaged shear-wave velocity in the top 30 m (V-S30), the topographical slope (slope), the fundamental resonance frequency (f(0)) and the depth beyond which V-s exceeds 800 m/s (H800). We considered first the performance of each SCP taken alone and then the combined performance of the 6 SCP pairs [V-S30-f(0)], [V-S30-H-800], [f(0)-slope], [H-800-slope], [V-S30-slope] and [f(0)-H-800]. This analysis is performed using a neural network approach including a random effect applied on a KiK-net subset for derivation of ground-motion prediction equations setting the relationship between various ground-motion parameters such as peak ground acceleration, peak ground velocity and pseudo-spectral acceleration PSA (T), and Mw, RJB, focal depth and SCPs. While the choice of SCP is found to have almost no impact on the median groundmotion prediction, it does impact the level of aleatory uncertainty. VS30 is found to perform the best of single proxies at short periods (T < 0.6 s), while f(0) and H-800 perform better at longer periods; considering SCP pairs leads to significant improvements, with particular emphasis on [V-S30-H-800] and [f(0)-slope] pairs. The results also indicate significant nonlinearity on the site terms for soft sites and that the most relevant loading parameter for characterising nonlinear site response is the "stiff" spectral ordinate at the considered period.
The creation of building exposure models for seismic risk assessment is frequently challenging due to the lack of availability of detailed information on building structures. Different strategies have been developed in recent years to overcome this, including the use of census data, remote sensing imagery and volunteered graphic information (VGI). This paper presents the development of a building-by-building exposure model based exclusively on openly available datasets, including both VGI and census statistics, which are defined at different levels of spatial resolution and for different moments in time. The initial model stemming purely from building-level data is enriched with statistics aggregated at the neighbourhood and city level by means of a Monte Carlo simulation that enables the generation of full realisations of damage estimates when using the exposure model in the context of an earthquake scenario calculation. Though applicable to any other region of interest where analogous datasets are available, the workflow and approach followed are explained by focusing on the case of the German city of Cologne, for which a scenario earthquake is defined and the potential damage is calculated. The resulting exposure model and damage estimates are presented, and it is shown that the latter are broadly consistent with damage data from the 1978 Albstadt earthquake, notwithstanding the differences in the scenario. Through this real-world application we demonstrate the potential of VGI and open data to be used for exposure modelling for natural risk assessment, when combined with suitable knowledge on building fragility and accounting for the inherent uncertainties.
On April 29, 2017 at 0:56 UTC (2:56 local time), an M (W) = 2.8 earthquake struck the metropolitan area between Leipzig and Halle, Germany, near the small town of Markranstadt. The earthquake was felt within 50 km from the epicenter and reached a local intensity of I (0) = IV. Already in 2015 and only 15 km northwest of the epicenter, a M (W) = 3.2 earthquake struck the area with a similar large felt radius and I (0) = IV. More than 1.1 million people live in the region, and the unusual occurrence of the two earthquakes led to public attention, because the tectonic activity is unclear and induced earthquakes have occurred in neighboring regions. Historical earthquakes south of Leipzig had estimated magnitudes up to M (W) ae 5 and coincide with NW-SE striking crustal basement faults. We use different seismological methods to analyze the two recent earthquakes and discuss them in the context of the known tectonic structures and historical seismicity. Novel stochastic full waveform simulation and inversion approaches are adapted for the application to weak, local earthquakes, to analyze mechanisms and ground motions and their relation to observed intensities. We find NW-SE striking normal faulting mechanisms for both earthquakes and centroid depths of 26 and 29 km. The earthquakes are located where faults with large vertical offsets of several hundred meters and Hercynian strike have developed since the Mesozoic. We use a stochastic full waveform simulation to explain the local peak ground velocities and calibrate the method to simulate intensities. Since the area is densely populated and has sensitive infrastructure, we simulate scenarios assuming that a 12-km long fault segment between the two recent earthquakes is ruptured and study the impact of rupture parameters on ground motions and expected damage.