To provide quantitative information on the ground acceleration necessary to break speleothems, laboratory measurements on samples of stalagmite have been performed to study their failure in bending. Due to their high natural frequencies, speleothems can be considered as rigid bodies to seismic strong ground motion. Using this simple hypothesis and the determined mechanical properties (a minimum value of 0.4 MPa for the tensile failure stress has been considered), modelling indicates that horizontal acceleration ranging from 0.3 m/s2 to 100 m/s2 (0.03 to 10g) are necessary to break 35 broken speleothems of the Hotton cave for which the geometrical parameters have been determined. Thus, at the present time, a strong discrepancy exists between the peak accelerations observed during earthquakes and most of the calculated values necessary to break speleothems. One of the future research efforts will be to understand the reasons of the defined behaviour. It appears fundamental to perform measurements on in situ speleothems.

A palynological study of the Chalk Group (Campanian-Danian) in the Meer borehole (northern Belgium), which penetrated the thickest succession known in the Campine Basin, has revealed diverse, well-preserved organic-walled dinoflagellate cyst assemblages. The succession contains numerous chronostratigraphically significant dinocyst events, which are based mainly on the highest consistent occurrences of index species. At least 35 bio-events have enabled a subdivision into nine intervals, at stage or substage level, within the Campanian to Danian interval, as based on comparison with coeval assemblages elsewhere in northwest Europe, inclusive of stratotypes of stages and stage boundaries. Bio-events allow correlation of the section studied with the Campanian Exochosphaeridium? masureae, Areoligera coronata and Samlandia mayi zones, the Maastrichtian Pervosphaeridium tubuloaculeatum, Deflandrea galeata and Hystrichostrogylon coninckii zones, and the Danian Damassadinium californicum Zone. In addition, a correlation with other zonal schemes for the southern North Sea Basin and with conventional northwest European belemnite zones is presented. Comparisons with Boreal and Tethyan realms confirm that most bio-events may also be useful for interregional and global correlation. The Campanian-Danian dinocyst biostratigraphy of the Meer borehole is put alongside geophysical well logs and an ecozonation, in order to check the validity of lithostratigraphical correlations across the Campine Basin. This first, detailed correlation attempt shows that sensitivities to facies change associated with differences in accommodation space and sediment supply appear to be at the base of slight, yet consistent, shifts between the local lithological succession and the standard lithostratigraphical scheme of the Maastricht type area.

Accurate dating is necessary to get insight in the temporal variations in sediment deposition in floodplains. The interpretation of such dates is however dependent on the fluvial architecture of the floodplain. In this study we discuss the fluvial architecture of three contrasting Belgian catchments (Dijle, Geul and Amblève catchment) and how this influences the dating possibilities of net floodplain sediment storage. Although vertical aggradation occurred in all three floodplains during the last part of the Holocene, they differ in the importance of lateral accretion and vertical aggradation during the entire Holocene. Holocene floodplain aggradation is the dominant process in the Dijle catchment. Lateral reworking of the floodplain sediments by river meandering was limited to a part of the floodplain, resulting in stacked point bar deposits. The fluvial architecture allows identifying vertical aggradation without erosional hiatuses. Results show that trends in vertical floodplain aggradation in the Dijle catchment are mainly related to land use changes. In the other two catchments, lateral reworking was the dominant process, and channel lag and point bar deposits occur over the entire floodplain width. Here, tracers were used to date the sediment dynamics: lead from metal mining in the Geul and iron slag from ironworks in the Amblève catchment. These methods allow the identification of two or three discrete periods, but their spatial extent and variations is identified in a continuous way. The fluvial architecture and the limitation in dating with tracers hampered the identification of dominant environmental changes for sediment dynamics in both catchments. Dating methods which provide only discrete point information, like radiocarbon or OSL dating, are best suited for fluvial systems which contain continuous aggradation profiles. Spatially more continuous dating methods, e.g. through the use of tracers, allow to reconstruct past surfaces and allow to reconstruct reworked parts of the floodplain. As such they allow a better reconstruction of past sedimentation rates in systems with important lateral reworking.

The Late weichselian Pleniglacial wind regime in the eastern Netherlands is reconstructed by means of landform and sedimentological analysis. This analysis involves aeolian and fluvial landforms in the Dinkel river valley in the Twente region. The aeolian deposits considered here date from the Last Glacial Maximum (approximately 22 ka) to the start of the Belling Interstadial at 14.7 ka.

A major event in this period is the formation of a cryoturbation level caused by permafrost degradation, overlain by an erosional hiatus dated between 21 and 17 ka. Both features are attributed to a period of warmer and moister climate, causing permafrost degradation and erosion by surficial runoff. Thereafter aeolian activity prevailed under relatively arid conditions. A deflation surface was formed, the Beuningen Gravel Bed. This deflation surface is present in many Weichselian sections in the Netherlands and the adjacent parts of Belgium and Germany. The deflation occurred concurrently with deposition of coversand at other places.

The morphology of the coversand-landscape in the Dinkel valley was controlled by the relief of the pre-existing floodplain and the wind pattern. Coversand ridges consisting of low dunes accumulated near the margins of the active channel belt. Relatively thick sand sheets occur in the leesides of the ridges, thin sand sheets are found at greater distance.

Mainly westerly sand-transporting winds operated during winter and summer. In winter aeolian deposition occurred by frequent and strong easterly winds also. On the smallest, local scale, the pattern of deposition was determined by the topography and moisture of the receiving surface.

Coversand deposition came to an end with the formation of a sand sheet under relatively warm and less arid conditions. Coversand deposition continued into the Belling Interstadial; colonization of the coversand surface by vegetation probably has been delayed by nutrient-poor conditions.

Mosasaur taxa currently known from Campanian-Maastrichtian strata in the Liège-Limburg (SE Netherlands, NE Belgium) and Mons (southern Belgium) basins are listed and briefly discussed and their stratigraphic ranges indicated. Recently published and/or ongoing work on coleoid and ammonoid cephalopods in these areas allows tie points between NW Europe and the United States (Western Interior, Gulf Coast, Atlantic Seaboard) to be established. Future studies need to refine the resultant, rather crude, scheme. The ultimate aim is a more robust picture of mosasaur taxonomy, biostratigraphy and palaeobiogeography, and a detailed evaluation of migratory patterns across the Atlantic.

Field mapping of fluvial terraces, aerial photographs, ground penetrating radar and seismic data from gas and oil exploration were used at four different locations to detect young tectonic and halokinetic movements in the North-German-Basin.

i) The course of the Rivers Weser and Aller follow precisely a shallow Tertiary graben on the northwestern flank of the Verden salt diapir. Recent local depressions and vegetation anomalies on the alluvial plain have the same orientation as the strike direction of the faults at subsurface depth. Apparently, the river follows tectonic lines, and thus the river sediments can be used for the interpretation of recent crustal movements.

ii) The Wedehof diapir, in contrast, is topped by a local topographic high which follows exactly the shape of the underlying salt. Either the diapir formed an obstacle for the advance of the continental glaciers or one has to assume halokinetic uplift of more than 50 m during the post-Saalian Pleistocene. Either way, the Wedehof diapir shows control of the modern surface morphology by halokinesis.

iii) The course of the river Hunte, in contrast, outside the area of salt diapirism, shows anomalies of incision and terrace width over a local updoming caused by tectonic inversion of distinct blocks in the basin. The confluence of several tributaries of the Hunte lies exactly over the updoming of Barnstorf. Thus, the rivers do not avoid the local high, but focus in this area, which is characterised by a graben on top of the domestructure, as visible in seismic profiles. Again, tectonism controls river development.

iv) The last case study is from Lake Plön, where seismic profiles reveal that linear shorelines of the lake parallel the flanks of two local graben structures of Tertiary age. It is apparent that the Weichselian glaciers that formed the lake and the surrounding moraines interacted with the existing grabens.

The Tertiary morphology in the North German basin was apparently draped by Quaternary glacial deposits, but rivers and lakes that dominate the topography of the modern landscape still reflect the geodynamic centers of Tertiary tectonism and halokinesis. Faults from the depth of the Tertiary penetrate the Quaternary strata and allow upward fluid migration, which becomes visible on aerial photographs as linear vegetation anomalies.

Carinodens is among the smallest known mosasaurs, and one of the most elusive at that. This enigmatic taxon is known from only two dentary fragments and just over a dozen of isolated teeth. Based on this material, an attempt is made to improve our understanding of the dietary habits of this mosasaur, by comparing tooth and jaw morphology to extant analogues, and by testing a biomechanical model with ‘feeding’ experiments using an artificial mosasaur jaw equipped with a force gauge. Carinodens appears to have been a durophagous mosasaur, capable of crushing small molluscs and arthropods, although its dietary habits may not necessarily have been limited to hard-shelled food.

Dendrochronology has been used to date sub-fossil wood excavated from former wetland woods located in Zwolle-Stadshagen, the Netherlands. Tree-ring analysis was performed on oak and ash to gain insight into changes in the growth dynamics of the trees resulting from changes in external site conditions, particularly hydrology. In addition population dynamics of oak and ash were studied by observing the temporal distribution and mean age of the oak and ash trees in the period when the woodland existed. The results showed that woodlands with oak and ash at Zwolle-Stadshagen existed at least for a period of about 700 years, from ca. 150 BC to AD 580 (ca. 2200 - 1400 cal. BP). The oak and ash were slowly growing and showed alternating phases of normal and depressed growth. Growth depressions occurred with a frequency of about 20 to 40 years. Almost all of these depressions were synchronous with those in oaks from other - mainly wetland areas - in NW Europe. Consequently, a regional (climate) factor must be considered responsible for this phenomenon. By taking into account the specific growth conditions in wetland woods it is very likely that this factor is linked to site hydrology. Around AD 300 (ca. 1650 cal. BP) a remarkable shift in population dynamics occurred in both oak and ash with many new trees establishing. This indicates a period of striking changes in site hydrology preceding the collapse of the woodland starting around AD 530 (ca. 1420 cal. BP).

The Late Glacial landscape of the Netherlands was a landscape with changing geomorphology and vegetation. Glacial, eolian and fluvial processes in the time before the Late Glacial initially had formed the main landscape types that still existed during the Late Glacial. In these landscape types, geomorphological processes persisted, particularly during intervals when the vegetation cover was discontinuous. Vegetation development initiated soil formation and stabilised the substratum. On the other hand, the abiotic landscape influenced vegetation development, and particularly vegetation patterns.

The Late Glacial vegetation patterns, changing in time, have been reconstructed on the basis of a palynological database containing the data from over 250 pollen diagrams from the Netherlands. Patterns of vegetation and abiotic landscape appear to compare to each other in many instances, indicating the close interrelationship between vegetation and the abiotic landscape.

The investigation of foraminiferal assemblages from a series of Oligocene borehole sections allowed paleoenvironment and paleoclimate reconstructions for the Rupelian and Chattian (Lower and Upper Oligocene) Stages in their type region, the southern North Sea Basin. A striking feature coinciding with the Rupelian-Chattian (R-C) unconformity is the major change in paleotemperature and paleobathymetry. The shallow marine to restricted marine subtropical fauna at the base of the Chattian is in strong contrast with the deeper marine and cooler upper Rupelian assemblages. This study suggests that the early Chattian transgression is genetically related to a widespread major warming pulse, known as the Late Oligocene Warming Event.

Earthquakes associated with gas production have been recorded in the northern part of the Netherlands since 1986. The Huizinge earthquake of 16 August 2012, the strongest so far with a magnitude of ML = 3.6, prompted reassessment of the seismicity induced by production from the Groningen gas field. An international research programme was initiated, with the participation of many Dutch and international universities, knowledge institutes and recognised experts.

The prime aim of the programme was to assess the hazard and risk resulting from the induced seismicity. Classic probabilistic seismic hazard and risk assessment (PSHA) was implemented using a Monte Carlo method. The scope of the research programme extended from the cause (production of gas from the underground reservoir) to the effects (risk to people and damage to buildings). Data acquisition through field measurements and laboratory experiments was a substantial element of the research programme. The existing geophone and accelerometer monitoring network was extended, a new network of accelerometers in building foundations was installed, geophones were placed at reservoir level in deep wells, GPS stations were installed and a gravity survey was conducted.

Results of the probabilistic seismic hazard and risk assessment have been published in production plans submitted to the Minister of Economic Affairs, Winningsplan Groningen 2013 and 2016 and several intermediate updates. The studies and data acquisition further constrained the uncertainties and resulted in a reduction of the initially assessed hazard and risk.

The results of a detailed seismic mapping campaign of 13 horizons in the northwestern German North Sea, covering Late Permian to Palaeogene sedimentary successions, are presented. Based on the interpretation of four 3D and two 2D seismic surveys, thickness and depth maps of prominent stratigraphic units were constructed. These maps provide an overview of key structural elements, the sedimentation and erosion, and give insights into the evolution of the German Central Graben. The base of the Zechstein Group reaches a maximum depth of 7800 m within the German Central Graben. Lateral thickness variations in the Zechstein reflect the extensive mobilisation of Zechstein salt. Complex rift-related structures, with the Central Graben as the main structural element, were found not later than the Early Triassic. Up to 3000-m thick Triassic sediments are preserved in the eastern German Central Graben of which 1800 m consist of Keuper sediments. The Lower Buntsandstein unit shows increasing thicknesses towards the southeastern study area, likely related to distinct lateral subsidence. As a consequence of uplift of the North Sea Dome, Middle Jurassic sediments were eroded in large parts of the northwestern German North Sea and are only preserved in the German Central Graben. The NNW–SSE oriented John Basin is another important structural element, which shows maximum subsidence during the Late Jurassic. In most parts of the study area Lower Cretaceous sediments are absent due to either erosion or non-deposition. Lower Cretaceous deposits are preserved in the Outer Rough Basin in the northwest and within the German Central Graben. Upper Cretaceous sediments are found at depths between 1500 and 3600 m, reaching a maximum thickness of approximately 1600 m on the Schillgrund High. Contraction and inversion of pre-existing Mesozoic faults during the Late Cretaceous is distinct at the Schillgrund Fault, i.e. the eastern border fault of the Central Graben. The Palaeogene is predominantly a period of strong basin subsidence. Within 37 Myrs, up to 1400 m of Palaeogene sediments were deposited in the northwesternmost part of the study area. Detailed mapping of salt structures enables a reconstruction of halokinetic movements over time and a deciphering of the influence of the Zechstein salt on the sedimentary evolution during the Mesozoic and Cenozoic. Increasing sediment thicknesses in rim-synclines indicate that most of the salt structures in the German Central Graben had their main growth phase during the Late Jurassic.

We report here a new elasmosaurid from the early Maastrichtian at Bentiaba, southern Angola. Phylogenetic analysis places the new taxon as the sister taxon to Styxosaurus snowii, and that clade as the sister of a clade composed of (Hydrotherosaurus alexandrae (Libonectes morgani + Elasmosaurus platyurus)). The new taxon has a reduced dorsal blade of the scapula, a feature unique amongst elasmosaurids, but convergent with cryptoclidid plesiosaurs, and indicates a longitudinal protraction-retraction limb cycle rowing style with simple pitch rotation at the glenohumeral articulation. Morphometric phylogenetic analysis of the coracoids of 40 eosauropterygian taxa suggests that there was a broad range of swimming styles within the clade.

The Glueckstadt Graben is a prominent structure of the Central European Basin System, where the sedimentary patterns are extensively affected by Permian salt movements. The relations of the sedimentary patterns to salt structures have been analyzed through present-day distributions of sediments. In addition, a three-dimensional backward modelling approach has been applied to determine the original salt distribution in response to the unloading due to sequential backstripping of the stratigraphic layers. The results of the modelling reveal the thickness distribution of the Permian salt for 5 time intervals from the end of the Triassic to present day. Spatial agreement has been found between the development of the depleted zone of the Permian salt through time and the observed distribution of the maximum subsidence for the different stratigraphic units above the salt. The sedimentation centres for each time interval are always located above the zone of reduced or depleted Permian salt. In the central part of the Glueckstadt Graben, the depletion occurred already in the Triassic and perfectly correlates with the thickest Triassic. During the Jurassic, Cretaceous and Tertiary, the areas of depleted Permian salt shifted towards the basin flanks, and the same occurred with the centres of maximum sediment deposition. Thus, the results of the modelling strongly support the conclusion that salt withdrawal has played a major role during the Meso-Cenozoic evolution of the Glueckstadt Graben and that the progressive depletion of the Permian salt layer, from the central part towards the margins, created the large part of the accommodation space for sedimentation in addition to tectonic subsidence.

Furthermore, our study has several important implications for salt behaviour in different tectonic settings. In general, the results of modelling indicate a good correlation between the main phases of salt movements and tectonic events in the area under consideration. During the Triassic, the first stage of diapirism in the Glueckstadt Graben occurred within the central part of the basin. Regional extension may have triggered reactive diapirism and caused the formation of the deep primary rim synclines. Once the salt structures had reached the critical size, buoyancy forces supported their continued growth until the Jurassic when extension-induced regional stresses once more affected the Glueckstadt Graben. The results of the modelling indicate very little salt activity during the late Early Cretaceous-early Late Cretaceous when the area of the Glueckstadt Graben was tectonically silent. Therefore, our study supports the concept of tectonically induced salt movements which can be interrupted during the absence of tectonic forces. Salt movements were reactivated in the marginal troughs by compressional forces during the latest Late Cretaceous-Early Cenozoic. Paleogene-Neogene salt withdrawal led to the growth of N-S oriented salt structures mainly at the margins of the basin. This phase of salt tectonics correlates temporally with almost W-E extension. This indicates a renewed change in tectonic regime after Late Cretaceous-Early Cenozoic compression.

Micromorphology is used to examine and compare a Late Jurassic diamictite from northeast Scotland with a Pleistocene diamict from southern Ontario, Canada in order to test if a statistical difference between diamicts can be recognized and used to separate differing types of diamicts/diamictites. The diamictites from Scotland have been ascribed to various depositional agencies occurring in several distinctly differing terrestrial and marine palaeoenvironments. In contrast, the Pleistocene diamicton is regarded as a subglacial till. Both diamicts appear remarkably similar visually and contain many corresponding features such as macrostructures, and exotic and fractured subangular to subrounded clasts. Micromorphology is used to re-examine these diamicts/diamictites at the microscopic level to detect if the palaeoenvironments within which they were deposited can be ascertained. In this paper a quantitative assessment of microstructures using micromorphology is developed. Comparative statistical analyses of these diamicts, using micromorphological features, reveals that the Jurassic diamictites are non-glacigenic, non-terrestrial and most likely deposited within a marine environment as a result of subaquatic debris mass movement, while, in contrast, the Pleistocene diamicts were most likely subglacial tectomicts deposited beneath the active base of the Laurentide Ice Sheet.

Exploration in a mature basin requires a detailed classification and standardisation of rock stratigraphy to adequately comprehend the depositional history and prospect architecture. The pre-Quaternary Stratigraphic Nomenclature of the Netherlands compiled by Van Adrichem Boogaert & Kouwe in 1993 provided a consistent framework for use by the Dutch geological community. Over the past twenty years, new biostratigraphic techniques and continued exploration in the Netherlands have provided additional stratigraphic information. Based on this information the Late Jurassic lithostratigraphy in particular, shows significant inaccuracies. The Callovian-Ryazanian strata from the northern offshore of the Netherlands' territorial waters, termed the Central Graben Subgroup and Scruff Group, reveal a complex sedimentary history. The combination of non-marine to shallow marine lateral facies changes, repetitive log and facies characteristics in time, sea-level and climate change, salt tectonics and structural compartmentalisation hamper straightforward seismic interpretation and log correlation. Recognition of three genetic sequences by Abbink et al. in 2006 enabled an improved reconstruction of the geological history. Further improvements in refinement and reliability of the stratigraphy together with new information on the facies and ages of the successions and about the subsequent tectonostratigraphic development of the northern Dutch offshore area form the basis of the present revision. As a result, earlier lithostratigraphic models have been changed and new lithostratigraphic relationships and names are introduced in this paper.

Sediments in north-central Texas, ranging in age from >117 to 85 Ma, represent a variety of terrestrial and marine depositional settings. Isotopic analyses of wood fragments found throughout the section allow correlation to the standard secular marine δ13C curve because of characteristic peaks at the Aptian-Albian and Cenomanian-Turonian boundaries. Consistency of the north-central Texas δ13C curve with the marine standard facilitates correlation among non-marine and marine environments on a global scale. Radiometrically dated ammonite zones recognised in Texas provide calibration for the Cenomanian and Turonian portions of the section. Cenomanian and Turonian sediments in north-central Texas preserve the oldest (96 Ma) and the youngest (<85 Ma) well-documented Coniasaurus, a dolichosaur also known from the southern North Sea Basin during that interval. Haasiasaurus, the oldest known well-documented early mosasaur, is found at ‘Ein Yabrud, Israel (98 Ma), followed by other poorly dated Cenomanian taxa from the eastern Mediterranean region, and then by Dallasaurus turneri and Russellosaurus coheni in Texas (92 Ma) and Tethysaurus (90.5 Ma) in Morocco. Neither shifts in δ13C nor large-scale sea level change seem to have influenced dolichosaur or mosasaur evolution in substantial ways during the Cenomanian and Turonian stages.

A first trench has been excavated for paleoseismological analysis in the German part of the Roer Valley graben, which has experienced several historical earthquakes with a maximum intensity up to VIII on the MSK-scale.The trench has exposed the Rurrand fault as a complex fault zone with at least five separate, SW-dipping, normal fault strands displacing an early Pleistocene terrace of the Rhine river by more than 7 m. The major part of the observed deformation was produced during or after deposition of an overlying unit of stratified loess of middle Weichselian to probably Saalian age. The faulting history is shown to be episodic, with different fault strands active at different times. Growth faulting that would be indicative of continuous, aseismic fault motion has not been observed. Our stratigraphic control is not sufficient to constrain the timing and to provide evidence of the coseismic nature for each observed fault displacement. However, two units of structureless, gravelly loess are interpreted as the result of extensive solifluction triggered by two large surface-rupturing events. This is suggested by the position of these units, which is controlled by the main faults, and by their remarkably young age (< 400 cal. BC), indicated by radiocarbon and OSL datings and by the presence of historic brick fragments. At least two faults show moderate activity that is even younger. Our interpretation is not in agreement with earlier hypotheses that ongoing vertical movements of circa 1 mm/a in the German part of the Lower Rhine graben are the result of aseismic fault creep, but is in line with the results of similar investigations on the southwestern border fault of the Roer Valley Graben in Belgium, which demonstrates the need for further paleoseismological research in this region. The Rurrand fault is presently experiencing aseismic slip on its superficial portion, induced by extensive groundwater lowering for mining purposes. This ongoing deformation seems to be expressed in the trench as diffuse bundles of anastomosing cracks extending up to, and in some cases even into the plough zone, rather than as sharp fault planes which are typical of older, tectonic fault movements.

Tylosaurus kansasensis sp. nov. is described herein on the basis of thirteen specimens collected from the Smoky Hill Chalk (upper Coniacian) of western Kansas, USA. The new species, originally designated Tylosaurus n. sp., co-occurred with T. nepaeolicus and exhibits a number of primitive characters that place it in a basal position in the mosasaur phylogeny. Among the key differences separating this species from other tylosaurines are a shortened, more rounded pre-dental process of the premaxilla, a distinctive quadrate lacking an infrastapedial process, and a parietal foramen located adjacent to the frontal-parietal suture.

A resource assessment methodology has been developed to designate prospective high permeable clastic aquifers and to assess the amount of potential geothermal energy in the Netherlands. It builds from the wealth of deep subsurface data from oil and gas exploration and production which is publicly and digitally available. In the resource assessment various performance indicator maps have been produced for direct heat applications (greenhouse and spatial heating). These maps are based on detailed mapping of depth, thickness, porosity, permeability, temperature and transmissivity (methodology presented in other papers in this NJG issue). In the resource assessment analysis 14 lithostratigraphic units (clastic aquifers) have been considered, ranging in age from the Permian to the Cenozoic. Performance maps have been made which include a) the expected doublet power (MWth) to be retrieved; b) the number of houses or hectares that can be heated from one doublet; and c) a potential indicator map, which provides insight in subsurface suitability for specific applications from a techno-economic perspective. To obtain a nationwide overview of the resource potential in terms of recoverable geothermal energy, a progressive filtering approach was used from total heat content of the reservoirs (Heat In Place – HIP) via the heat that can potentially be recovered (Potential Recovery Heat – PRH) to energy maps taking into account a techno-economic performance evaluation (Recoverable Heat – RH). Results show that the HIP is approximately 820,000 PJ which is significantly more than previous estimates of around 90,000 PJ. This considerable increase in geothermal energy potential is the result of accurate geological mapping of key reservoir properties and the development of state-of-the-art techno-economic performance assessment tools that performs Monte Carlo simulation. Moreover, for the previous estimates boundary conditions were set with the aim to compare the geothermal potential between different EU countries (Rijkers & Van Doorn, 1997). Taking into account techno-economic aspects, the RH is in the order of 85,000 PJ. This is equivalent to ~70% of the ultimate recoverable gas of the Slochteren Gas field. In total over 400 maps have been created or used as input for the resource assessment. Together, they provide comprehensive information for geothermal energy development from various stakeholder perspectives. The maps can be interactively assessed in the web-based portal ThermoGIS (www.thermogis.nl). This application complements existing subsurface information systems available in the Netherlands and supports the geothermal community in assessing the feasibility of a geothermal system on a regional scale.