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

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

Uncertainty assessment applied to marine subsurface datasets

A recently released voxel model quantifying aggregate resources of the Belgian part of the North Sea includes lithological properties of all Quaternary sediments and modelling-related uncertainty. As the underlying borehole data come from various sources and cover a long time-span, data-related uncertainties should be accounted for as well. Applying a tiered data-uncertainty assessment to a composite lithology dataset with uniform, standardized lithological descriptions and rigorously completed metadata fields, uncertainties were qualified and quantified for positioning, sampling and vintage. The uncertainty on horizontal positioning combines navigational errors, on-board and off-deck offsets and underwater drift. Sampling-gear uncertainty evaluates the suitability of each instrument in terms of its efficiency of sediment yield per lithological class. Vintage uncertainty provides a likelihood of temporal change since the moment of sampling, using the mobility of fine-scale bedforms as an indicator. For each uncertainty component, quality flags from 1 (very uncertain) to 5 (very certain) were defined and converted into corresponding uncertainty percentages meeting the input requirements of the voxel model. Obviously, an uncertainty-based data selection procedure, aimed at improving the confidence of data products, reduces data density. Whether or not this density reduction is detrimental to the spatial coverage of data products, will depend on their intended use. At the very least, demonstrable reductions in spatial coverage will help to highlight the need for future data acquisition and to optimize survey plans. By opening up our subsurface model with associated data uncertainties in a public decision support application, policy makers and other end users are better able to visualize overall confidence and identify areas with insufficient coverage meeting their needs. Having to work with a borehole dataset that is increasingly limited with depth below the seabed, engineering geologists and geospatial analysts in particular will profit from a better visualization of data-related uncertainty

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

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

Corrigendum : Differential subsidence in the urbanised coastal-deltaic plain of the Netherlands (Geologie en Mijnbouw/Netherlands Journal of Geosciences (2018) DOI: 10.1017/njg.2018.11)

An error has occurred in Table 1 of this article (Koster, K; Stafleu, J; Stouthamer, E). The correct Table 1 is given

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

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

Distribution of Holocene Marine Mud and Its Relation to Damage from the 1923 Earthquake Disaster in the Tokyo Metropolitan Area, Japan

Tokyo, which is located near the boundary between the North American and Philippine Sea plates, has been frequently struck by large earthquakes throughout the Holocene. The 1923 Taisho Kanto Earthquake is a rare historical earthquake that can be reconstructed in detail because abundant datasets were collected by investigations performed just after the earthquake. We examined 13,000 borehole logs from the Tokyo and Nakagawa lowlands to clarify the distribution and thickness of incised-valley fills and soft marine mud that had accumulated since the Last Glacial Maximum (LGM) on a grid with a resolution of 150 m × 150 m. We compared these datasets with the distribution of wooden house damage ratios caused by the Taisho Kanto Earthquake. Our results showed that the thickness of the soft mud, but not that of the incised-valley fills, was strongly correlated with the wooden house damage ratio. The mud content was >60%, water content was >30%, and S-wave velocity was ca. 100 m/s in the soft Holocene marine mud. The wooden house damage ratio was highest where the soft mud thickness was 20 m, because in those areas, both the soft mud and the wooden houses resonated with a natural period of ca. 1 s

R. Soekmono : Gurah. The link between the Central and the East-Javanese arts M. Soekarto Kartoatmodjo : The discovery of three new inscriptions in the district of Klaten (South central Java)

Lombard Denys. R. Soekmono : Gurah. The link between the Central and the East-Javanese arts M. Soekarto Kartoatmodjo : The discovery of three new inscriptions in the district of Klaten (South central Java). In: Bulletin de l'Ecole française d'Extrême-Orient. Tome 60, 1973. pp. 423-424

Models created as 3-D cellular voxel arrays

This chapter describes the different aspects of voxel modeling based on the GeoTOP program developed and run by TNO - Geological Survey of the Netherlands. The GeoTOP program produces a national voxel model that describes the architecture and properties of the shallow subsurface down to a maximum depth of 50 m below ordnance datum. The GeoTOP modeling procedure is followed by discussions on how to evaluate model uncertainty. The chapter describes the way in which voxel models can be populated with properties such as hydraulic conductivity, geotechnical parameters, organic matter or even absolute geological age. Property attributes convert the lithostratigraphic-unit voxel models into powerful instruments for a wide range of applications including groundwater management, ground risk assessment, the planning of infrastructure works, and aggregate assessments. A basic geological voxel has lithostratigraphic and lithologic attributes. Voxels can subsequently be populated with any parameter that correlates with combinations of these basic attributes, and in some cases depth