Skip to main content
Ronald  Harting
  • Princetonlaan 6
    3584 CB Utrecht
  • +31 88 866 48 24
The accelerations experienced at the surface as a result of earthquakes induced by the production of gas from the Groningen gasfield depend on the local shallow geological and soil conditions. This is called the 'site response effect'. In... more
The accelerations experienced at the surface as a result of earthquakes induced by the production of gas from the Groningen gasfield depend on the local shallow geological and soil conditions. This is called the 'site response effect'. In order to improve our knowledge of this effect, the NAM invited Deltares to build a detailed model of the shallow subsurface below Groningen. This report prepared by Deltares describes the quaternary geology of the Groningen area.
Research Interests:
Download (.pdf)
Het beschikbaar maken van kennis en informatie van de Nederlandse ondergrond is een kerntaak van de Geologische Dienst Nederland - TNO (GDN). Het is daarbij van groot belang dat deze informatie is gebaseerd op een consistente dataset,... more
Het beschikbaar maken van kennis en informatie van de Nederlandse ondergrond is een kerntaak van de Geologische Dienst Nederland - TNO (GDN). Het is daarbij van groot belang dat deze informatie is gebaseerd op een consistente dataset, waarin parameters over de ondergrond worden beschreven en geïnterpreteerd. In Nederland is voor veel toepassingen geen homogene, landsdekkend bruikbare collectie parameters beschikbaar. Een centrale opslag ontbreekt en de wel beschikbare gegevens zijn met verschillende (vaak niet precies bekende) methoden verzameld. Om een solide verzameling basisgegevens samen te stellen is de GDN gestart met het systematisch verzamelen van deze data - in heel Nederland – door het uitvoeren van steek-boringen, deze te bemonsteren, te analyseren en te rapporteren. Daarbij is de aandacht vooral gericht op de lithologie (grondsoort met alle daarin aanwezige kenmerken), korrelgrootte, doorlatendheid, geochemische en geotechnische eigenschappen.
In dit rapport worden lithologische gegevens beschreven en geanalyseerd die door de GDN op een systematische wijze zijn verzameld in Noord-Nederland binnen het onderzoeksprogramma TopIntegraal. Doel is om voor de geologische eenheden, die voorkomen in het bereik tot 30 meter beneden maaiveld, gegevens over de lithologische samenstelling te geven en deze geologisch te interpreteren. Zowel de “ruwe” gegevens als de interpretatie beogen de gebruikers een middel te geven om maatschappelijke vraagstukken die met de ondergrond te maken hebben, beter te kunnen onderzoeken en op te lossen. In Noord-Nederland zijn 75 boringen geplaatst, die lithostratigrafisch zijn ingedeeld. Uit de boringen zijn in totaal meer dan 2500 monsters genomen en geanalyseerd op o.a. korrelgrootte verdeling, organische stof en een aantal andere parameters (o.a. doorlatendheid en geochemische analyses). Deze laatste worden separaat gerapporteerd. De resultaten hebben geleid tot een consistente dataset, die een uniek beeld geeft van de lithologische kenmerken van de verschillende geologische eenheden in Noord-Nederland.
Research Interests:
Download (.pdf)
20150224_DEF_TNO_2014-R10680_BIJLAGEN.pdf
20150224_DEF_TNO_2014-R10680.pdf
Kruiver, P.P., Wiersma, A., Kloosterman, F.H., De Lange, G., Korff, M., Stafleu, J., Busschers, F.S., Harting, R., Gunnink, J.L., Green, R.A., Van Elk, J. and Doornhof, D., 2017. Characterisation of the Groningen subsurface for seismic... more
Kruiver, P.P., Wiersma, A., Kloosterman, F.H., De Lange, G., Korff, M., Stafleu, J., Busschers, F.S., Harting, R., Gunnink, J.L., Green, R.A., Van Elk, J. and Doornhof, D., 2017. Characterisation of the Groningen subsurface for seismic hazard and risk modelling. Netherlands Journal of Geosciences,  96: 215-233, doi:10.1017/njg.2017.11

The shallow subsurface of Groningen, the Netherlands, is heterogeneous due to its formation in a Holocene tidal coastal setting on a periglacially and glacially inherited landscape with strong lateral variation in subsurface architecture. Soft sediments with low, small-strain shear wave velocities (V S30 around 200 m s −1) are known to amplify earthquake motions. Knowledge of the architecture and properties of the subsurface and the combined effect on the propagation of earthquake waves is imperative for the prediction of geohazards of ground shaking and liquefaction at the surface. In order to provide information for the seismic hazard and risk analysis, two geological models were constructed. The first is the 'Geological model for Site response in Groningen' (GSG model) and is based on the detailed 3D GeoTOP voxel model containing lithostratigraphy and lithoclass attributes. The GeoTOP model was combined with information from boreholes, cone penetration tests, regional digital geological and geohydrological models to cover the full range from the surface down to the base of the North Sea Supergroup (base Paleogene) at ∼800 m depth. The GSG model consists of a microzonation based on geology and a stack of soil stratigraphy for each of the 140,000 grid cells (100 m × 100 m) to which properties (V S and parameters relevant for nonlinear soil behaviour) were assigned. The GSG model serves as input to the site response calculations that feed into the Ground Motion Model. The second model is the 'Geological model for Liquefaction sensitivity in Groningen' (GLG). Generally, loosely packed sands might be susceptible to liquefaction upon earthquake shaking. In order to delineate zones of loosely packed sand in the first 40 m below the surface, GeoTOP was combined with relative densities inferred from a large cone penetration test database. The marine Naaldwijk and Eem Formations have the highest proportion of loosely packed sand (31% and 38%, respectively) and thus are considered to be the most vulnerable to liquefaction; other units contain 5–17% loosely packed sand. The GLG model serves as one of the inputs for further research on the liquefaction potential in Groningen, such as the development of region-specific magnitude scaling factors (MSF) and depth–stress reduction relationships (r d).
Download (.pdf)
Amsterdam is situated on the coastal-deltaic plain of the western Netherlands. Its geographical position brought the city prosperity, but also created huge challenges associated with heterogeneous and often adverse ground conditions. This... more
Amsterdam is situated on the coastal-deltaic plain of the western Netherlands. Its geographical position brought the city prosperity, but also created huge challenges associated with heterogeneous and often adverse ground conditions. This paper explores the geology of Amsterdam to a depth of c. 100 m, based on the output of the 3D geological subsurface models DGM and GeoTOP. The model results are used to create a geological map of the area, to determine the extent and depth of the foundation levels that are in use for buildings in the city centre and to detect the source of filling sand on which part of the more recent expansion of the city was founded. It is shown that subsurface conditions have had a profound effect on both landscape development and historical city growth. Geomodels like DGM and GeoTOP provide an easily accessible way to obtain a better understanding of the shallow subsurface.
Research Interests:
Download (.pdf)
Download (.pdf)
A key sedimentary unit from the northern Netherlands is the Saalian glacial till (in Dutch: “keileem”). Its shallow occurrence, lithological heterogeneity and relatively high hydraulic resistance make it an important unit to be considered... more
A key sedimentary unit from the northern Netherlands is the Saalian glacial till (in Dutch: “keileem”). Its shallow occurrence, lithological heterogeneity and relatively high hydraulic resistance make it an important unit to be considered in regional and local groundwater models. In 2013, TNO – Geological Survey of the Netherlands completed a new high-resolution 3D model of the glacial till, encompassing geometry, lithological composition and hydraulic resistance (Vernes et al., 2013). The high data density generally permits the till model to be used on a sub-regional to near-local scale.
Besides the direct application in the regional MIPWA groundwater model, the construction of the till model has resulted in several geologically relevant observations. This presentation discusses new insights regarding the geometry and heterogeneity of the Saalian glacial till. We show the exceptional variation in till depth and thickness over very short distances that was mapped in 3D for the first time at a regional scale, including the directionality present in the till geometry. Also, lithological variation of the till can now be studied over a range of local to regional scales. Furthermore, model construction has led to the spatial mapping of a second, older Saalian till in the south-east of Drenthe.
Research Interests:
Download (.pdf)
The Dutch shallow subsurface is extensively being used for natural aggregate resources, groundwater extraction and construction works. This makes that insights into the composition and heterogeneity of the subsurface and its physical and... more
The Dutch shallow subsurface is extensively being used for natural aggregate resources, groundwater extraction and construction works. This makes that insights into the composition and heterogeneity of the subsurface and its physical and chemical properties are of vital importance for Dutch society (Vernes et al. 2010, Van de Meulen et al. 2013). Although a large amount of physical and chemical data of the subsurface has been collected in the past, issues like discontinuous sample distributions and/or differences in analysis techniques, have limited their helpfulness for applied research (Harting et al. 2014). We present new results from an analysis of a vast nationwide data set of lithological properties from the shallow subsurface of the Northern Netherlands that were collected, measured and analysed in a systematic and uniform way (Bosch et al. 2015). The dataset, which is managed by the Geological Survey of the Netherlands, contains a large number of parameters that were collected from undisturbed sediments taken from 75 continuous boreholes with a depth up-to 40m. Since the cores have an excellent geographical distribution and the samples are measured on in situ samples with a standard set of techniques, the data overcomes the data density and measurement inconsistencies known from previous studies. In this presentation we will focus on the lithological properties of the Boxtel, Drente, Drachten and Peelo Formations in the northern Netherlands that were deposited under periglacial, glacial, eolian and glacial to glacio-lacustrine conditions respectively. We found that Weichselian (Boxtel Fm.) and early Saalian (Drachten Fm.) eolian sediments are nearly identical, suggesting comparable sedimentary processes during deposition. In several cores we newly identified a Elsterian till belonging to the Peelo Formation, which is very different from the Saalian till belonging to the Gieten Member of the Drente Formation. This difference may be related to the sediment sources available during formation of the till. Finally, we describe a new type of bedded till from the Saalian Gieten Member, strongly contrasting with the well-known massive till. With these examples, we show that using our lithofacies based approach (Harting et al. 2014) is a strong and powerful tool to study and better understand heterogeneity in these sediments (Bosch et al. 2015).
Download (.pdf)
Harting_INQUAPeribaltic_Abstract.pdf
Harting_INQUA-Peribaltic2015_LithologicalPropertiesGlacialSedimentsNNL.pdf
Society has an increasing demand from the subsurface, which in the Dutch shallow subsurface (upper 30 to 40 meters) mainly focuses on natural aggregate resources, groundwater, infrastructure and dike safety. This stimulates the demand for... more
Society has an increasing demand from the subsurface, which in the Dutch shallow subsurface (upper 30 to 40 meters) mainly focuses on natural aggregate resources, groundwater, infrastructure and dike safety. This stimulates the demand for knowledge about the composition and heterogeneity of the subsurface and its physical and chemical properties, including the uncertainties involved. Physical and chemical properties of sediments in the subsurface have been under investigation for decades; however, the usefulness of this data for applied research and the understanding of these properties is limited. This  is due to several factors: studies consist mainly of separately collected datasets, targeted at a limited amount of parameters, focused on a small number of geological units, distributed unevenly with depth and usually collected from clustered drillings with limited spatial extent or are analysed with different techniques and methods, often on disturbed samples. These factors result in a heterogeneous and biased dataset not suitable to function as a reference dataset or to statistically determine regional characteristics of geological units.

To overcome these shortcomings, the Geological Survey of the Netherlands is establishing a nation-wide reference dataset for physical and chemical properties. In 2006, a drilling campaign was started using cone penetration tests, cored drillings and geophysical well logs, choosing the sites for a good geographical distribution. The lithological properties of the undisturbed cores are visually described and interpreted for lithostratigraphy and inferred sedimentary environment based on lithofacies. The location of the samples in the cores are chosen based on this description and interpretation, resulting in an evenly distributed dataset of in situ samples with respect to geological units as well as an adequate number of samples suitable for statistical analysis. Analyses are uniformly performed for grain size distribution, permeability (both high and low permeable lithologies) and geochemical methods (X-Ray Fluorescence, Thermo-Gravimetric Analysis, Total Carbon, Total Sulphur and Total Organic Carbon). These analyses result in a large number of lithological, hydrological and geochemical parameters, i.e. clay content, sand median, vertical and horizontal permeability and CaCO3-content.

We present the results from the analysis of lithological properties for the Northern Netherlands. Besides geology, these properties can be applied directly in studies concerning (amongst others) groundwater, natural aggregates and dike safety. We demonstrate the use of sedimentary environments based on lithofacies as a useful tool for comparison between lithostratigraphic units and lithofacies. These lithofacies match distinct parts of the marine, fluvial, glacial, eolian or organogenic environment, i.e. tidal channel sand, floodbasin clay and subglacial till. This results in lithological properties illustrating the heterogeneity within a geological unit and between equal depositional environments in different lithostratigraphic units.
The acquired data have so far been used in several applied studies, i.e. improving parameterisation of 3D models leading to increased accuracy in groundwater models and dike safety studies concerning dike failure due to undermining. Recently, grain size distributions measured with different methods were recalibrated into a homogeneous dataset using this reference set, which greatly enlarged the dataset to be incorporated in the parameterisation of a 3D voxel model.
Download (.pdf)
The NAM is preparing a new "Winningsplan", to be submitted in 2016. For this new Winningsplan, a new generation of Ground Motion Prediction Equations (GMPEs) will be developed. The overall scope is to reduce uncertainties in the hazard... more
The NAM is preparing a new "Winningsplan", to be submitted in 2016. For this new Winningsplan, a new generation of Ground Motion Prediction Equations (GMPEs) will be developed. The overall scope is to reduce uncertainties in the hazard and risk analysis by improvement of input data, such as Groningen-specific data, and better GMPEs. In the current GMPE, only one value for shear wave velocity (Vs) is used for the entire Groningen field (Vs30 = 200 mIs). The shallow subsurface of Groningen, consisting of Holocene and Pleistocene sediments is heterogeneous, resulting in variations of shear wave velocity. It is expected that part of the uncertainties in the seismic hazard and risk analysis can be reduced by including Groningen-specific information and knowledge of the subsurface to improve quantification of the site response caused by earthquakes

Deltares has built a geological model for the Groningen field (+ 5 km buffer) for the purpose of the construction of Vs30 maps and as input for the calculations of site amplification. These results will feed into the new GMPEs. The Geological model for the ~ite response at the Groningen Field (GSG-model) is, among other data sources, based on the beta version of GeoTOP (a 3D geological model of the Netherlands), provided by TNO Geological Survey of the Netherlands. The GSG-model built by Deltares consists of a map defining geological areas and voxel stacks containing stratigraphy and lithological class with depth. Additionally, a state-of-the-art Vs30 map was derived for the Groningen field + 5 km buffer, taking into account Groningen-specific Vs relations and the geology from the GSG-model.

This report describes the method for the construction and the results of version 1 of the GSGmodel, the quality checks performed on the model and recommendations for future versions. When more data becomes available, updates of the GSG-model are anticipated.
Research Interests:
Download (.pdf)