The Geological Survey of the Netherlands aims at building a 3D geological voxel model of the upper 30 m of the subsurface of the Netherlands in order to provide a sound basis for subsurface related questions on, amongst others, groundwater extraction and management, land subsidence studies, aggregate resources and infrastructural issues. The Province of Zeeland (SW Netherlands, covering an area of approximately 70 by 75 km) was chosen as the starting point for this model due to an excellent dataset of 23,000 stratigraphically interpreted borehole descriptions.The modelling procedure involved a number of steps. The first step is a geological schematisation of the borehole descriptions into units that have uniform sediment characteristics, using lithostratigraphical, lithofacies and lithological criteria. During the second modelling step, 2D bounding surfaces are constructed. These surfaces represent the top and base of the lithostratigraphical units and are used to place each voxel (100 by 100 by 0.5 metres) in the model within the correct lithostratigraphical unit. The lithological units in the borehole descriptions are used to perform a final 3D stochastic interpolation of lithofacies, lithology (clay, sand, peat) and if applicable, sand grain-size class within each lithostratigraphical unit. After this step, a three-dimensional geological model is obtained. The use of stochastic techniques such as Sequential Gaussian Simulation and Sequential Indicator Simulation, allowed us to compute probabilities for lithostratigraphy, lithofacies and lithology for each voxel, providing a measure of model uncertainty.The procedures described above resulted in the first fully 3D regional-scale lithofacies model of the shallow subsurface in the Netherlands. The model provides important new insights on spatial connectivity of sediment units of, for example, sandy Holocene tidal channel systems. Our results represent a major step forward towards a fully 3D voxel model of the Netherlands.

A regional shear-wave velocity (V_S) model has been developed for the Groningen gas field in the Netherlands as the basis for seismic microzonation of an area of more than 1000 km². The V_S model, extending to a depth of almost 1 km, is an essential input to the modelling of hazard and risk due to induced earthquakes in the region. The detailed V_S profiles are constructed from a novel combination of three data sets covering different, partially overlapping depth ranges. The uppermost 50 m of the V_S profiles are obtained from a high-resolution geological model with representative V_S values assigned to the sediments. Field measurements of V_S were used to derive representative V_S values for the different types of sediments. The profiles from 50 to 120 m are obtained from inversion of surface waves recorded (as noise) during deep seismic reflection profiling of the gas reservoir. The deepest part of the profiles is obtained from sonic logging and V_P–V_S relationships based on measurements in deep boreholes. Criteria were established for the splicing of the three portions to generate continuous models over the entire depth range for use in site response calculations, for which an elastic half-space is assumed to exist below a clear stratigraphic boundary and impedance contrast encountered at about 800 m depth. In order to facilitate fully probabilistic site response analyses, a scheme for the randomisation of the V_S profiles is implemented.

Over the last ten to twenty years, geological surveys all over the world have been entangled in a process of digitisation. Their paper archives, built over many decades, have largely been replaced by electronic databases. The systematic production of geological map sheets is being replaced by 3D subsurface modelling, the results of which are distributed electronically. In the Netherlands, this transition is both being accelerated and concluded by a new law that will govern management and utilisation of subsurface information. Under this law, the Geological Survey of the Netherlands has been commissioned to build a key register for the subsurface: a single national database for subsurface data and information, which Dutch government bodies are obliged to use when making policies or decisions that pertain to, or can be affected by the subsurface. This requires the Survey to rethink and redesign a substantial part of its operation: from data acquisition and interpretation to delivery. It has also helped shape our view on geological surveying in the future.The key register, which is expected to start becoming operational in 2015, will contain vast quantities of subsurface data, as well as their interpretation into 3D models. The obligatory consultation of the register will raise user expectations of the reliability of all information it contains, and requires a strong focus on confidence issues. Building the necessary systems and meeting quality requirements is our biggest challenge in the upcoming years. The next step change will be towards building 4D models, which represent not only geological conditions in space, but also processes in time such as subsidence, anthropogenic effects, and those associated with global change.

The urbanised peat-rich coastal-deltaic plain of the Netherlands is severely subsiding due to human-induced phreatic groundwater level lowering, as this causes peat layers to compress and oxidise. To determine the potential susceptibility of this area to future subsidence by peat compression and oxidation, the effects of lowering present-day phreatic groundwater levels were quantitatively evaluated using a subsidence model. Input were a 3D geological subsurface voxel-model, modelled phreatic groundwater levels, and functions for peat compression and oxidation. Phreatic groundwater levels were lowered by 0.25 and 0.5m, and the resulting peat compression and oxidation over periods of 15 and 30 years were determined. The model area comprised the major cities Amsterdam and Rotterdam, and their surrounding agricultural lands.The results revealed that for these scenarios agricultural areas may subside between 0.3 and 0.8m; potential subsidence in Amsterdam and Rotterdam is considerably lower, less than 0.4m. This is due to the presence of several metres thick anthropogenic brought-up soils overlying the peat below the urban areas, which has already compressed the peat to a depth below groundwater level, and thus minimises further compression and oxidation. In agricultural areas peat is often situated near the surface, and is therefore highly compressible and prone to oxidation. The averaged subsidence rates for the scenarios range between 7 and 13mma⁻¹, which is corresponds to present-day rates of subsidence in the peat areas of the Netherlands. These results contrast with the trend of coastal-deltaic subsidence in other deltas, with cities subsiding faster than agricultural areas. This difference is explained by the driver of subsidence: in other deltas, subsidence of urban areas is mainly due to deep aquifer extraction, whereas in the Netherlands subsidence is due to phreatic groundwater level lowering.