scholarly journals Development of a seismic site-response zonation map for the Netherlands

2022 ◽  
Vol 22 (1) ◽  
pp. 41-63
Author(s):  
Janneke van Ginkel ◽  
Elmer Ruigrok ◽  
Jan Stafleu ◽  
Rien Herber

Abstract. Earthquake site response is an essential part of seismic hazard assessment, especially in densely populated areas. The shallow geology of the Netherlands consists of a very heterogeneous soft sediment cover, which has a strong effect on the amplitude of ground shaking. Even though the Netherlands is a low- to moderate-seismicity area, the seismic risk cannot be neglected, in particular, because shallow induced earthquakes occur. The aim of this study is to establish a nationwide site-response zonation by combining 3D lithostratigraphic models and earthquake and ambient vibration recordings. As a first step, we constrain the parameters (velocity contrast and shear-wave velocity) that are indicative of ground motion amplification in the Groningen area. For this, we compare ambient vibration and earthquake recordings using the horizontal-to-vertical spectral ratio (HVSR) method, borehole empirical transfer functions (ETFs), and amplification factors (AFs). This enables us to define an empirical relationship between the amplification measured from earthquakes by using the ETF and AF and the amplification estimated from ambient vibrations by using the HVSR. With this, we show that the HVSR can be used as a first proxy for site response. Subsequently, HVSR curves throughout the Netherlands are estimated. The HVSR amplitude characteristics largely coincide with the in situ lithostratigraphic sequences and the presence of a strong velocity contrast in the near surface. Next, sediment profiles representing the Dutch shallow subsurface are categorised into five classes, where each class represents a level of expected amplification. The mean amplification for each class, and its variability, is quantified using 66 sites with measured earthquake amplification (ETF and AF) and 115 sites with HVSR curves. The site-response (amplification) zonation map for the Netherlands is designed by transforming geological 3D grid cell models into the five classes, and an AF is assigned to most of the classes. This site-response assessment, presented on a nationwide scale, is important for a first identification of regions with increased seismic hazard potential, for example at locations with mining or geothermal energy activities.

2021 ◽  
Author(s):  
Janneke van Ginkel ◽  
Elmer Ruigrok ◽  
Jan Stafleu ◽  
Rien Herber

Abstract. Earthquake site-response is an essential part of seismic hazard assessment, especially in densely populated areas. The shallow geology of the Netherlands consists of a very heterogeneous soft sediment cover, which has a strong effect on seismic wave propagation and in particular on the amplitude of ground shaking, resulting in significant damage on structures despite the fact that the events are of small magnitude. Even though it is a low-to-moderate seismicity area, the seismic risk cannot be neglected, in particular, because shallow induced earthquakes occur. The aim of this study is to establish a nationwide site-response zonation by using the lithostratigraphy, earthquake- and ambient vibration recordings. In the first step, we constrain the parameters (velocity contrast and shear-wave velocity) that are indicative of ground-motion amplification in the Groningen area. For this, we combine ambient vibration and earthquake recordings using resp. the horizontal-to-vertical spectral ratio method (HVSR), borehole empirical transfer functions (ETFs) and amplification factors (AFs). This enables us to define an empirical relationship between measured earthquake amplification from the ETF and AF, and amplification estimated with the HVSR derived from the ambient seismic field. Therewith, we show that the HVSR can be used as a first proxy for amplification. Subsequently, HVSR curves throughout the Netherlands are estimated. The resulting peak amplitudes largely coincide with the in-situ lithostratigraphic sequences and the presence of a strong velocity contrast in the near-surface. Next, sediment profiles representing the Dutch shallow subsurface are categorized into five classes, where each class is representing a level of expected amplification. The mean amplification for each class, and its variability, is quantified using 66 sites with measured earthquake amplification (ETF and AF) and 115 sites with HVSR curves. The site-response (amplification) zonation map for the Netherlands is designed by transforming published geological 3D grid cell models into the five classes and an AF is assigned to most of the classes. This presented site-response assessment on a national scale is important for a first identification of regions with increased seismic hazard potential, for example at locations with mining or geothermal energy activities.


2021 ◽  
Author(s):  
Janneke van Ginkel ◽  
Elmer Ruigrok ◽  
Rien Herber

<p>Local site conditions can strongly influence the level of amplification of ground-motion at the surface during an earthquake. Especially near-surface low velocity sediments overlying stiffer seismic bedrock modify earthquake ground motions in terms of amplitudes and frequency content, the so-called site response. Earthquake ground-motion site response is of great concern because it can lead to amplified surface shaking resulting in significant damage on structures despite small magnitude events. The Netherlands has tectonically related seismic activity in the southern region with magnitudes up to 5.8 measured so far. In addition, gas extraction in the Groningen field in the northern part of the Netherlands, is regularly causing shallow (3 km), low magnitude (Mw max= 3.6), induced earthquakes. The shallow geology of the Netherlands consists of a very heterogeneous soft sediment cover, which has a strong effect on seismic wave propagation and in particular on the amplitude of ground shaking.</p><p> </p><p>The ambient seismic field and local earthquakes recorded over 69 borehole stations in Groningen are used to define relationships between the subsurface lithological composition, measured shear-wave velocity profiles, horizontal-to-vertical spectral ratios (HVSR) and empirical transfer functions (ETF). For the Groningen region we show that the HVSR matches the ETF well and conclude that the HVSR can be used as a first proxy for earthquake site-response. In addition, based on the ETFs we observe that most of the seismic wave amplification occurs in the top 50 m of the much thicker sediment layer. Here, a velocity contrast is present between the very soft Holocene clays and peat on top of the stiffer Pleistocene sands.</p><p> </p><p>Based on the learnings from Groningen we first constructed sediment type classes for the Dutch subsurface, each class representing a level of expected amplification. Secondly, the HVSR curves are estimated for all surface seismometers in the Netherlands seismic network and a sediment class is assigned to each location. Highest HVSR peak amplitudes are measured at sites with the highest level of amplification of the sediment classification. Based on this correlation and the presence of a detailed shallow geological model at most sites in the Netherlands, a simplistic approach is presented to predict amplification at any location with sufficient lithologic information. With this approach based on the shallow sediment composition, we can obtain constraints on the seismic hazard in areas that have limited data availability but have potential risk of seismicity, for example due to geothermal energy extraction.</p>


2006 ◽  
Vol 87 (1-2) ◽  
pp. 105-121 ◽  
Author(s):  
Torild van Eck ◽  
Femke Goutbeek ◽  
Hein Haak ◽  
Bernard Dost

2015 ◽  
Vol 42 (7) ◽  
pp. 463-476 ◽  
Author(s):  
Hadi Ghofrani ◽  
Gail M. Atkinson ◽  
Luc Chouinard ◽  
Philippe Rosset ◽  
Kristy F. Tiampo

Montreal has significant seismic risk due to the combination of moderate seismicity, high population density, and vulnerable infrastructure. An important tool in damage and risk assessment is a scenario shakemap, which shows the expected ground shaking intensity distribution patterns. In this study, we use regional ground motion and site response evaluations to generate scenario shakemaps for Montreal. The impact of event location on expected ground motions and intensities was tested by considering the occurrence of a scenario (a given magnitude event) at various locations, where the scenarios are defined based on an analysis of the most likely future event locations. Variability in near surface geology plays an important role in earthquake ground shaking; we use microzonation information from Montreal to assess the expected site amplification effects. The results of this study may be used as input to seismic risk studies for Montreal.


Author(s):  
Carlos Mendoza ◽  
Stephen Hartzell

ABSTRACT We invert the shear‐wave displacement spectra obtained from 30 three‐component, broadband waveforms recorded within 300 km of the 6 November 2011 Mw 5.7 Prague, Oklahoma, earthquake to recover the site‐response contribution using an inversion method that simultaneously inverts for source, path, and site effects. Site‐response functions identify resonant frequencies within a range of 0.1–10 Hz that generally coincide with spectral peaks in horizontal‐to‐vertical ratio curves derived from the recorded waveforms. S‐wave velocity profiles available for several sites were also used to calculate theoretical SH transfer functions that predict the site amplification due to the near‐surface soil structure down to depths of 30–50 m. The transfer functions do not provide resonance information below about 5–8 Hz, indicating that the spectral peaks in the site response obtained from the waveform analysis result from deeper velocity variations. A 0.3 Hz spectral peak observed at several stations, for example, coincides with the strong, surface‐wave amplitudes observed at 3 s periods for induced M≥3 earthquakes in Oklahoma and Kansas, suggesting that this resonant peak may be due to surface waves trapped in the upper ∼2  km sedimentary layer of the crust. Both shallow and deep contributions to the site response are important for the characterization of ground motion from central and eastern North America (CENA) earthquakes. We obtain a corner frequency of 0.229, consistent with independent observations of the size of the event. A frequency‐dependent attenuation relation of Q(f)=1107f0.398 consistent with prior CENA path measurements is also derived.


Solid Earth ◽  
2019 ◽  
Vol 10 (2) ◽  
pp. 379-390 ◽  
Author(s):  
Yaniv Darvasi ◽  
Amotz Agnon

Abstract. Instrumental strong motion data are not common around the Dead Sea region. Therefore, calibrating a new attenuation equation is a considerable challenge. However, the Holy Land has a remarkable historical archive, attesting to numerous regional and local earthquakes. Combining the historical record with new seismic measurements will improve the regional equation. On 11 July 1927, a rupture, in the crust in proximity to the northern Dead Sea, generated a moderate 6.2 ML earthquake. Up to 500 people were killed, and extensive destruction was recorded, even as far as 150 km from the focus. We consider local near-surface properties, in particular, the shear-wave velocity, as an amplification factor. Where the shear-wave velocity is low, the seismic intensity far from the focus would likely be greater than expected from a standard attenuation curve. In this work, we used the multichannel analysis of surface waves (MASW) method to estimate seismic wave velocity at anomalous sites in Israel in order to calibrate a new attenuation equation for the Dead Sea region. Our new attenuation equation contains a term which quantifies only lithological effects, while factors such as building quality, foundation depth, topography, earthquake directivity, type of fault, etc. remain out of our scope. Nonetheless, about 60 % of the measured anomalous sites fit expectations; therefore, this new ground-motion prediction equation (GMPE) is statistically better than the old ones. From our local point of view, this is the first time that integration of the 1927 historical data and modern shear-wave velocity profile measurements improved the attenuation equation (sometimes referred to as the attenuation relation) for the Dead Sea region. In the wider context, regions of low-to-moderate seismicity should use macroseismic earthquake data, together with modern measurements, in order to better estimate the peak ground acceleration or the seismic intensities to be caused by future earthquakes. This integration will conceivably lead to a better mitigation of damage from future earthquakes and should improve maps of seismic hazard.


2021 ◽  
pp. 875529302110194
Author(s):  
Daniel Verret ◽  
Denis LeBœuf ◽  
Éric Péloquin

Eastern North America (ENA) is part of a region with low-to-moderate seismicity; nonetheless, some significant seismic events have occurred in the last few decades. Recent events have reemphasized the need to review ENA seismicity and ground motion models, along with continually reevaluating and updating procedures related to the seismic safety assessment of hydroelectric infrastructures, particularly large dams in Québec. Furthermore, recent researchers have shown that site-specific characteristics, topography, and valley shapes may significantly aggravate the severity of ground motions. To the best of our knowledge, very few instrumental data from actual earthquakes have been published for examining the site effects of hydroelectric dam structures located in eastern Canada. This article presents an analysis of three small earthquakes that occurred in 1999 and 2002 at the Denis-Perron (SM-3) dam. This dam, the highest in Québec, is a rockfill embankment structure with a height of 171 m and a length of 378 m; it is located in a narrow valley. The ground motion datasets of these earthquakes include the bedrock and dam crest three-component accelerometer recordings. Ground motions are analyzed both in the time and frequency domains. The spectral ratios and transfer functions obtained from these small earthquakes provide new insights into the directionality of resonant frequencies, vibration modes, and site effects for the Denis-Perron dam. The crest amplifications observed for this dam are also compared with previously published data for large dams. New statistical relationships are proposed to establish dam crest amplification on the basis of the peak ground acceleration (PGA) at the foundation.


Geophysics ◽  
2012 ◽  
Vol 77 (4) ◽  
pp. E301-E315 ◽  
Author(s):  
Thomas Kalscheuer ◽  
Juliane Hübert ◽  
Alexey Kuvshinov ◽  
Tobias Lochbühler ◽  
Laust B. Pedersen

Magnetotelluric (MT), radiomagnetotelluric (RMT), and, in particular, controlled-source audiomagnetotelluric (CSAMT) data are often heavily distorted by near-surface inhomogeneities. We developed a novel scheme to invert MT, RMT, and CSAMT data in the form of scalar or tensorial impedances and vertical magnetic transfer functions simultaneously for layer resistivities and electric and magnetic galvanic distortion parameters. The inversion scheme uses smoothness constraints to regularize layer resistivities and either Marquardt-Levenberg damping or the minimum-solution length criterion to regularize distortion parameters. A depth of investigation range is estimated by comparing layered model sections derived from first- and second-order smoothness constraints. Synthetic examples demonstrate that earth models are reconstructed properly for distorted and undistorted tensorial CSAMT data. In the inversion of scalar CSAMT data, such as the determinant impedance or individual tensor elements, the reduced number of transfer functions inevitably leads to increased ambiguity for distortion parameters. As a consequence of this ambiguity for scalar data, distortion parameters often grow over the iterations to unrealistic absolute values when regularized with the Marquardt-Levenberg scheme. Essentially, compensating relationships between terms containing electric and/or magnetic distortion are used in this growth. In a regularization with the minimum solution length criterion, the distortion parameters converge into a stable configuration after several iterations and attain reasonable values. The inversion algorithm was applied to a CSAMT field data set collected along a profile over a tunnel construction site at Hallandsåsen, Sweden. To avoid erroneous inverse models from strong anthropogenic effects on the data, two scalar transfer functions (one scalar impedance and one scalar vertical magnetic transfer function) were selected for inversion. Compared with a regularization of distortion parameters with the Marquardt-Levenberg method, the minimum-solution length criterion yielded smaller absolute values of distortion parameters and a horizontally more homogeneous distribution of electrical conductivity.


2017 ◽  
Vol 33 (2) ◽  
pp. 481-498 ◽  
Author(s):  
Julian J. Bommer ◽  
Peter J. Stafford ◽  
Benjamin Edwards ◽  
Bernard Dost ◽  
Ewoud van Dedem ◽  
...  

The potential for building damage and personal injury due to induced earthquakes in the Groningen gas field is being modeled in order to inform risk management decisions. To facilitate the quantitative estimation of the induced seismic hazard and risk, a ground motion prediction model has been developed for response spectral accelerations and duration due to these earthquakes that originate within the reservoir at 3 km depth. The model is consistent with the motions recorded from small-magnitude events and captures the epistemic uncertainty associated with extrapolation to larger magnitudes. In order to reflect the conditions in the field, the model first predicts accelerations at a rock horizon some 800 m below the surface and then convolves these motions with frequency-dependent nonlinear amplification factors assigned to zones across the study area. The variability of the ground motions is modeled in all of its constituent parts at the rock and surface levels.


Sign in / Sign up

Export Citation Format

Share Document