Interdisciplinary, landscape-oriented studies from an archaeological viewpoint in the Low Countries mainly concentrate on cultural and economic research questions. Focal points are the physical setting of settlements and cemeteries, land use patterns and subsistence economy. As a result, the collected data are rather site-based and one-dimensional. As a counterweight, this study aims to look beyond the boundaries of settlements and cemeteries by offering a regional and diachronic perspective on the development of the landscape, vegetation and habitation of Twente (the Netherlands) since the Late Glacial. A detailed search for existing pollen data yielded 125 sites containing information from a wide variety of sampling contexts. A series of six evidence-based regional vegetation maps have been constructed by analysing relations between pollen data, soil data and topography. The maps serve as first-stage generalised models that predict regional trends, allow subsequent testing and place site-specific archaeological data in a wider context. The method developed is applicable to other regions. A comparison with contemporary habitation patterns, based on archaeological and historical data, reveals spatio-temporal trends in human influence on vegetation and in physical factors influencing site location. Five maps have been ‘translated’ into artist impressions.

The history of the Amstel and Vecht rivers in the central part of the Netherlands during the Holocene has been debated for a considerable time (Bos et al., 2009; De Bont, 1994, 2000, 2008; De Bont & Kleij, 2012; De Gans, 2011; De Gans & Bunnik, 2011, 2012). Many publications in different disciplines have made their contribution to solve issues such as initiation and cessation of river flow, different and opposing river flow directions, and different chronologies.

This paper describes the background of five palaeogeographical maps between the Marsdiep and the Weser River, and discusses the natural and anthropogenic processes driving the coastal changes during the last part of the Holocene. Before 2500 BC, during the first half of the Holocene, tidal basins were formed in the lower lying Pleistocene valley system as a result of the Holocene sea-level rise. The tidal basins were filled during the second half of the Holocene and on the deposits from the Pleistocene in the hinterland large coastal peat bogs developed. These peat bogs were vulnerable and sensitive to marine ingressions when the peat surface subsided due to drainage, compaction and erosion. During the Subatlantic (450 BC to present), the different ingression systems in the coastal area between Marsdiep and Weser had their own histories in timing and evolution. The ingressions were naturally caused by lateral migration of coastal barrier and tidal-inlet systems or by changes in the natural drainage system in the hinterland. From the Late Iron Age onwards, humans started to be the major cause of ingressions. By reclaiming and cultivating the seaward margins of coastal peat bogs, these areas subsided significantly and were flooded by high storm surges. When coastal areas were embanked during the historical period, the situation for the lower lying peat lands became more dramatic. When the sea dikes breached, the peat land was flooded, leading to casualties and huge material damages and loss of land. Drowning of the peat lands of the Jade and Dollard in the 14th and 15th centuries are examples of such catastrophes.

The prehistoric Oer-IJ tidal system in the coastal area of Noord-Holland, between Castricum, Uitgeest, Velsen and Amsterdam, was the successor of the Haarlem tidal system when this was silted up around 3000 BC and large peat bogs formed in the Zaanstreek and Haarlemmermeer areas. Since then the Oer-IJ has been the natural outlet to the sea, draining the peat hinterlands. About 800 BC the Oer-IJ system was connected to the fluvial system of the Utrechtse Vecht and became the northern branch of the river Rhine. During the Late Iron Age, when the Flevo lakes in the IJsselmeer region and the Utrechtse Vecht were connected with the Wadden Sea, the Oer-IJ lost its discharge function. The tidal area silted up and was closed between 200 and 100 BC by a barrier ridge. The settlement history of the Oer-IJ system and the archaeological heritage in the subsurface is closely related to the geological and hydrological development throughout the ages. The shape and location of the continuously migrating outlet determined the opportunities for human settlements and activities. The best locations were beach ridges, higher, silted-up salt marshes and marginal zones of the peatland. In the Late Iron Age the sand flats also became habitable since tidal activity had stopped. In the Early Roman period there was no direct connection from the harbour of Castellum Flevum at Velsen to the North Sea, but ships could navigate from the Oer-IJ channel between Velsen and Amsterdam, through the Flevo lakes and the Utrechtse Vecht to the Wadden Sea, and to the Roman border (Limes) along the Oude Rijn. Here the data used for the palaeogeographic landscape reconstruction of the Oer-IJ are presented and explained, and the most important landscape-forming processes, which led to the emergence and closure of the Oer-IJ, are described. The landscape reconstructions give a new perspective on the migration of the main tidal channel and the formation of the tidal-inlet system near Castricum, which was the result of the progradation of the beach ridges south and north of the Castricummerpolder (Binnendelta according to De Roo, 1953). The geological and archaeological observations in the Binnendelta prove that the Oer-IJ was closed from the open sea in the early Late Iron Age.

Aggregate resource assessments, derived from three subsequent generations of voxel models, were compared in a qualitative way to illustrate and discuss modelling progress. We compared the models in terms of both methodology and usability. All three models were produced by the Geological Survey of the Netherlands. Aggregate is granular mineral material used in building and construction, and in this case consists of sand and gravel. On each occasion ever-increasing computer power allowed us to model at a higher resolution and use more geological information to constrain interpolations. The two oldest models, built in 2005 and 2007, were created specifically for aggregate resource assessments, the first as proof of concept, the second for an online resource information system. The third model was derived from the ongoing multipurpose systematic 3D modelling programme GeoTOP. We used a study area of 40 × 40 km located in the central Netherlands, which encompasses a section of the Rhine-Meuse delta and adjacent glacial terrains to the north. Aggregate resource assessments rely on the extent to which the occurrence and grain size of sand and gravel are resolved, and on proper representation of clay and peat layers (overburden and intercalations) that affect exploitability. Average model properties (e.g. total aggregate content) are about the same in all three models, except for a difference resulting from converting older lithological classifications to the current one. This difference illustrates that data selection and preparation are paramount, especially when dealing with quality issues. Generally speaking the results of the aggregate resource assessments are consistent and satisfactory for all three models, provided that they are judged at the appropriate scale. However, the assessments based on GeoTOP best approach the desired scale of use for the aggregates industry; in that sense progress was significant and each model was a better fit for the purpose.

Since 2003 extensive archaeological research has been conducted in Amsterdam, the Netherlands, connected with the initial phase of the new underground system (Noord/Zuidlijn). Research has mainly focused on two locations, Damrak and Rokin, in the centre of Medieval Amsterdam. Both sites are situated around the (former) River Amstel, which is of vital importance for the origin and development of the city of Amsterdam. Information on the Holocene evolution of the river, however, is relatively sparse. This project has provided new evidence combining archaeological and geological data, and allowed the reconstruction of six consecutive landscape phases associated with the development of the River Amstel. The course of the present-day Amstel is the result of a complex interaction of processes that started with an early prehistoric tidal gully within the Wormer Member of the Naaldwijk Formation, including Late Neolithic (2400–2000 BC) occupation debris in its fill that was subsequently eroded. Next, this system developed into a later prehistoric Amstel river course that was part of the Angstel–Vecht–Oer-IJ system (1020–350 BC), meandering through a peat-dominated landscape. Later on the processes included intensive reclamation of land, drainage and canalisation, although the Amstel was also strongly influenced by natural storm tides. After intense land reclamation, starting around 1200 AD, the meandering Amstel from Nes to Kalverstraat, which was originally 150 m wide, became the rather straight 20–50 m wide tamed, canalised river of today.

The ice caps that covered large parts of the continents of the northern hemisphere during the Pleistocene glaciations drained huge quantities of meltwater. In several places the erosive power of the meltwater rivers has led to the formation of ice-marginal valleys (IMVs). A much-debated question is whether sediments deposited in IMVs by proglacial and extraglacial streams can be distinguished on the basis of their heavy-mineral content. This question was assessed by an inventory of the heavy-mineral assemblages from the middle part of the Toruń-Eberswalde IMV in northwest Poland, two sandurs that supplied sediment from the north and the pre-Wisła river system that supplied sediment from the south; all these streams fed the IMV. The largely similar heavy-mineral compositions and sediments concentrations of the middle part of the IMV and sandurs suggest that the sediment in the IMV was supplied almost entirely by the streams on the sandurs but also that some sediments were eroded from the Miocene subsoil of the IMV itself and for a small part from the south by the pre-Wisła river system. The only heavy mineral in the pre-Wisła sediments for which the percentage is significantly different from those in the sediments of the sandurs and the IMV terrace is epidote. The difference, however, is not seen in the sediments of the IMV so it can be concluded that the sediment supply to the middle part of this IMV by streams from the south was insignificant. This is in contrast with what was hitherto commonly assumed.

Various trace fossils are known from the Upper Triassic (Rhaetian) shale deposits of western Europe, especially from Great Britain. Here we present a unique specimen, unknown from Rhaetian shales in western Europe thus far to our knowledge. The specimen consists of a string of small knobs collected from the dark-coloured Rhaetian shales from the eastern Netherlands, deposited in a marine, near-coastal environment. The specimen represents the first described trace fossil from these shales. The identity of this specimen appears enigmatic. However, SEM-EDS analysis showed that the string of knobs is pyritised and does not contain phosphorus nor did the sediment directly around the specimen, suggesting a non-coprolitic origin of the specimen. Eggs and larvae are also excluded as possibilities. The specimen closely resembles several trace fossils identified as burrows, which is why we favour this interpretation. The rare presence of trace fossils reinforces the hypothesis that the Dutch Rhaetian shales were deposited under a stresses regime with low oxygen conditions.

In Germany, mosasaur remains are very rare and only incompletely known. However, the earliest records date back to the 1830s, when tooth crowns were found in the chalk of the Isle of Rügen. A number of prominent figures in German palaeontology and geosciences of the 19th and 20th centuries focused on these remains, including, among others, Friedrich von Hagenow, Hermann von Meyer, Andreas Wagner, Hanns Bruno Geinitz and Josef Pompeckj. Most of these works were only short notes, given the scant material. However, the discovery of fragmentary cranial remains in Westphalia in 1908 led to a more comprehensive discussion, which is also of historical importance, as it illustrates the discussions on the highly controversial and radical universal phylogenetic theory proposed by Gustav Steinmann in 1908. This theory saw the existence of continuous lines of descent, evolving in parallel, and did not regard higher taxonomic units as monophyletic groups but as intermediate paraphyletic stages of evolution. In this idea, nearly all fossil taxa form part of these lineages, which extend into the present time, and natural extinction occurs very rarely, if ever. In Steinmann's concept, mosasaurs were not closely related to squamates but formed an intermediate member in a anagenetic chain from Triassic thalattosaurs to extant baleen whales. The newly found specimen led Josef Pompeckj to write a vehement rebuttal to Steinmann's theory, published in 1910, showing that his conclusions were conjectural and speculative, being based on convergence and not supported by scientific evidence. This particular specimen, housed in Göttingen, later also inspired a piece of palaeoart by Franz Roubal under the instructions of Othenio Abel.

With the exception of a vertebra from the Campanian of former East Prussia (now Russian Federation), and a possible vertebra from the Cenomanian of Dresden, Saxony, all datable material – today partly lost – originated from the northern part of present-day Germany and stratigraphically from the Campanian–Maastrichtian. The purported record from the Cenomanian of Bavaria (southeastern Germany) was most probably an error, based on Upper Jurassic crocodilian material.

Airborne electromagnetics (AEM) is an effective technique for geophysical investigations of the shallow subsurface and has successfully been applied in various geological settings to analyse the depositional architecture of sedimentary systems for groundwater and environmental purposes. However, interpretation of AEM data is often restricted to 1D inversion results imaged on resistivity maps and vertical resistivity sections. The integration of geophysical data based on AEM surveys with geological data is often missing and this deficiency can lead to uncertainties in the interpretation process. The aim of this study is to provide an improved methodology for the interpretation of AEM data and the construction of more realistic 3D geological subsurface models. This is achieved by the development of an integrated workflow and 3D modelling approaches based on combining different geophysical and geological data sets (frequency-domain helicopter-borne electromagnetic data (HFEM), time-domain helicopter-borne electromagnetic data (HTEM), three 2D reflection seismic sections and 488 borehole logs). We used 1D inversion results gained from both HFEM and HTEM surveys and applied a 3D resistivity gridding procedure based on geostatistical analyses and interpolation techniques to create continuous 3D resistivity grids. Subsequently, geological interpretations have been performed by combining with, and validation against, borehole and reflection seismic data. To verify the modelling results and to identify uncertainties of AEM inversions and interpretation, we compared the apparent resistivity values of the constructed 3D geological subsurface models with those of AEM field measurements. Our methodology is applied to a test site near Cuxhaven, northwest Germany, where Neogene sediments are incised by a Pleistocene tunnel valley. The Neogene succession is subdivided by four unconformities and consists of fine-grained shelf to marginal marine deposits. At the end of the Miocene an incised valley was formed and filled with Pliocene delta deposits, probably indicating a palaeo-course of the River Weser or Elbe. The Middle Pleistocene (Elsterian) tunnel valley is up to 350 m deep, 0.8–2 km wide, and incised into the Neogene succession. The unconsolidated fill of the Late Miocene to Pliocene incised valley probably formed a preferred pathway for the Pleistocene meltwater flows, favouring the incision. Based on the 3D AEM resistivity model the tunnel-valley fills could be imaged in high detail. They consist of a complex sedimentary succession with alternating fine- and coarse-grained Elsterian meltwater deposits, overlain by glaciolacustrine (Lauenburg Clay Complex) and marine Holsteinian interglacial deposits. The applied approaches and results show a reliable methodology, especially for future investigations of similar geological settings. The 3D resistivity models clearly allow a distinction to be made between different lithologies and enables the detection of major bounding surfaces and architectural elements.

A partial quadrate, here assigned to the globidensine mosasaur Prognathodon cf. saturator, is recorded from the basal Valkenburg Member (Maastricht Formation) in the type area of the Maastrichtian Stage (southeast Netherlands, northeast Belgium). Occurrences of Prognathodon saturator in this area are extremely rare, which might well be linked to the offshore habitat preferred by this species.

Around 800 AD the peatlands around Amstelland were drained by two rivers. The ‘northern’ Amstel discharged surplus water directly into Lake Almere, which from the 12th century onwards was enlarged to form the Sudersee (in Dutch: Zuiderzee). The ‘southern’ Amstel was a near-abandoned westernmost arm of the River Vecht discharging the poorly drained borderland around the river Vecht into Lake Almere/the Sudersee. As part of the reclamation of peatlands associated with agricultural activities between the end of the 10th and the mid-13th centuries, the rivers were connected via a canal, thus creating the river Amstel, as it is known today.

A new genus and species of euphaneropid, Ciderius cooperi, is recorded from the Lower Silurian Fish Bed Formation of the Midland Valley of Scotland on the basis of articulated material. Euphaneropids constitute an enigmatic group which is known from Devonian deposits in Scotland and Canada. The new find adds to our understanding of this group, in particular with regard to the morphology of some common anatomical elements. The paired head stains are here shown to contain clefts which can be interpreted as optic fissures, indicating that these are the remains of eyes. The anterior head stain is reconstructed and demonstrated to be a barrel-shaped object of an uncertain nature. Paired mineralisations situated in the posterior part of the cranial region appear to represent remains of the parachordals, while fossilised blood vessels might be preserved in the form of black lines which mark a greatly elongated branchial region, similar to euphaneropids. Serially repetitive rows of short horizontal stripes on the posterior half of the body are unique for the new taxo, but their interpretation is problematic.

The Amstel river is located in the peat-covered coastal plain of the Netherlands and gives Amsterdam its name (Amstel dam). It is a small secondary branch of the repeatedly bifurcating Rhine delta system. Historically, the Amstel debouched into the peat-fringed former Oer-IJ estuary, which was connected to the North Sea, but after the closure of this inlet the estuary was transformed into an inland sea (IJ) due to erosion of the adjacent peat. The Amstel river was active between 3000 BP and 1122 AD after which time the supply water from the Rhine was stopped due to the construction of a dam far upstream near Wijk bij Duurstede. On the basis of borehole data from various sources, four cross-sections were constructed in the Amstel branch to study the unknown lithology and lithostratigraphy of the Amstel sediments in the Amsterdam area. The deposits show the Amstel was a low-energy river which carried mainly clay. The cross-sections reveal that the Amstel in its downstream part was flanked by two lithologically identical layers of overbank clay, intercalated by a peat layer. The lowermost overbank clay was deposited from 3000 BP to about 1000 AD. The intercalated peat layer is estimated to have developed between the 11th and 12th centuries AD, indicating a decreased sediment supply in the Amstel, and rise of water level in the downstream river caused by Zuiderzee influences such as storms and tide. The uppermost overbank clay was deposited during major storm surges such as those documented in 1164 and 1170 AD, and was derived from the brackish Zuiderzee; it has been traced upstream along the Amstel for over 10 km. Near the mouth of the Amstel channel in the Oer-IJ estuary its bottom has been scoured by estuarine processes to a lower level. On the basis of archaeological and geological data it is argued that the Amstel channel of medieval Amsterdam had a water depth of about 6 m before the construction of a dam in the 13th century. Soil scientists, historical geographers and historians have argued that the Amstel once consisted of two separate rivers: a northern Oer-IJ connected channel draining from the Amsterdam Stopera to the north, and a southern peat draining channel draining from the Amsterdam-Watergraafsmeer to the south. The relatively straight stretch of the present-day Amstel now positioned within the urban area has been hypothesised to be man-made between the 11th and 13th centuries AD. In this paper, on the basis of geological arguments such as channel depths, overbank clays, peat composition and other characteristics, it is concluded that the Amstel had a natural channel in the Amsterdam area.

Estimates of particle-size made by operators in the field and laboratory represent a vast and relatively untapped data archive. The wide spatial distribution of particle-size estimates makes them ideal for constructing geological models and soil maps. This study uses a large data set from the Netherlands (n = 4837) containing both operator estimates of particle-size and complete particle-size distributions measured by laser granulometry. Operator estimates are inaccurate and imprecise relative to measured laser data; only 16.68% of samples were successfully classified using the Dutch classification scheme for fine sediment. Operator estimates of sediment particle-size encompass the same range of percentage values as particle-size distributions measured by laser analysis. However, the distributions measured by laser analysis show that most of the sand percentage values lie between 0 and 1, so the majority of the variability in the data is lost because operator estimates are made to the nearest 1% at best, and more frequently to the nearest 5%. Operator estimates made by three technicians trained by the Geological Survey of the Netherlands are found not to be influenced by bias, rather they exhibit very similar levels of accuracy and precision. This study compares five different methods of modelling complete particle-size distributions from sparse data: (i) a four-part Pearson's probability distribution function, (i) a log-linear interpolation, (iii) a logit-linear interpolation, (iv) a logistic probability distribution function and (v) a logit constrained cubic-spline (logit-CCS) interpolation. The logit-CCS interpolation performed best across all the samples used, although the performance of all models was very similar for normal Gaussian, skewed and peaked distributions. Predictions for bimodal distributions using the Pearson's, logit-linear and logistic models are markedly less accurate than both log-linear and logit-CCS interpolation models. Although the logit-CCS interpolation model produces the best predictions of continuous particle-size distributions, the low accuracy and precision of operator estimates does not warrant the use of such a complex algorithm. Given this, it is suggested that a standard log-linear interpolation is the most effective means of modelling complete particle-size distributions from sparse data. Interpolation-based models are recommended over probability distribution functions because they allow for a greater degree of flexibility and will always honour the available input data.

We report here a large mosasaur skull, preserved three-dimensionally in a concretion recovered from Ziebach County, South Dakota, USA. This fossil represents the first articulated mosasaur skull from the Trail City Member of the Fox Hills Formation and the first definitive occurrence of Mosasaurus hoffmanni Mantell, 1829 from that area and the northernmost occurrence in the Western Interior Seaway, greatly extending the paleobiogeographic range of this taxon. The age of this specimen is determined to be between 68.3 and 67.6 Ma based on the associated invertebrate fauna. Although previous authors have suggested synonomy of the North American Mosasaurus maximus Cope, 1869 and the European M. hoffmanni, this is the most comprehensive analysis to date and is based on comparisons with Mosasaurus specimens recovered across the Northern Hemisphere, allowing an emended diagnosis of the species M. hoffmanni. Minor differences are considered individual variation or to reflect ontogenetic stage, including slender dentaries in some individuals, range of development of the C-shaped notch of the coronoid and differences in the shape of the supratemporal fenestra.