Till Plain Research Articles

Geospatial analysis of groundwater arsenic and fluoride in Quaternary aquifers of southern Ontario, Canada

The geological setting of an area plays a critical role in the transfer and ultimate distribution of hydrochemical constituents present in groundwater. In southern Ontario, Canada, the present physiography was significantly influenced by glacial processes during the Quaternary period. The heterogeneous nature and complex pattern of shallow subsurface glacial overburden sediments, likely affect the fate of different groundwater constituents. In this study, arsenic (As) and fluoride (F−) concentrations from 515 water wells, that are constructed within overburden sediment, were analyzed with the physiographic map of southern Ontario along with other related variables. Geospatial mapping and several spatial statistical analyses were performed to examine the possible geological influence on As and F− distribution and water-well susceptibility at a regional scale. Key findings suggest four physiographic settings were significant variables influencing the distribution of As and F− in differently constructed bored/dug and drilled wells. Bored/dug wells in Bevelled Till Plains and bored/dug wells in Undrumlinized and Bevelled Till Plains were found to be relatively susceptible to As and F− contamination respectively. In contrast, bored/dug and drilled wells in Drumlinized Till Plains and Drumlins and drilled wells in Sand Plains seemed to be relatively safe from F− and As respectively. The statistical regression analyses suggested that other variables, such as the application of phosphate fertilizer and the textures of till, influenced the spatial distribution of As and F− as well as which types of wells (bored/dug or drilled) were impacted. The geospatial mapping and statistical cluster analysis indicated that the possible sources of elevated As and F− in drilled wells are the clasts of underlying bedrock. The relationship between physiographic settings and impacted overburden wells in southern Ontario provides planners with an approach to water-well susceptibility assessments at the regional scale, which in turn can guide further local analysis for water resource management.

Estimating natural soil drainage classes in the Wisconsin till plain of the Midwestern U.S.A. based on lidar derived terrain indices: Evaluating prediction accuracy of multinomial logistic regression and machine learning algorithms

Accurate spatial prediction of natural soil drainage condition is not only important for agriculture and hydrological modeling but also for installing subsurface drainage and onsite waste disposal systems. For this research, 154 sites were selected based on a stratified random sampling method. For each site, drainage class was identified based on visual examination of soil cores. A digital elevation model developed from lidar data was used to derive seven terrain indices. Terrain indices were used to predict drainage class using four prediction models: multinomial logistic regression and three machine learning algorithms (random forest, C5.0, and artificial neural network). Based on 30% random hold-back validation data, all digital soil mapping (DSM) models provided similar results. The artificial neural network provided relatively higher overall accuracy of 70%, a kappa coefficient of 0.59, and a Brier score of 0.34. The rest of the DSM models as well as SSURGO provided an accuracy that ranged between 64 and 66%, a kappa coefficient that ranged between 0.52 and 0.54, and a Brier score that ranged between 0.41 and 0.46.

Origin of glacial curvilineations by subglacial meltwater erosion: Evidence from the Stargard drumlin field, Poland

AbstractDeciphering the origin of subglacial landforms is crucial for understanding processes operating under ice sheets. Due to the inaccessibility of the ice/bed interface limiting real‐time observations, the origin of these landforms is typically interpreted from the morphological and sedimentological record exposed after ice retreat, which makes the conclusions conjectural. Here, we focus on a recently discovered new type of subglacial landform called glacial curvilineations (GCLs), whose formation generates some controversy. We use LiDAR imagery and geological data to constrain the origin of GCLs found in a tunnel valley in the Stargard drumlin field, Poland. The tunnel valley is about 3800 m long, 520 m wide, and it occupies about 2 km2. A total of 61 individual GCLs have been mapped there; they are up to 550 m long and typically 0.5–1.5 m high, but some reach a height of up to 5 m. The GCLs are composed of mostly massive till represented by either a single or several diamicton units occasionally separated by thin layers of meltwater sand. Till macrofabrics in the GCLs and in the nearby till plain are distinctly consistent, and have high and uniform clustering strengths. At all locations, the dominant till fabric direction broadly corresponds with the ice movement direction given by the drumlin orientations, but it is oblique to the orientation of the tunnel valley and independent of the orientations of the individual GCLs. The lack of correspondence between the orientation of the GCLs and the till fabrics in them suggests that the GCLs consist of an antecedent till deposited during the drumlinizing ice advance and subsequently carved by subglacial meltwater erosion in the tunnel valley. These results substantiate the hypothesis that GCLs are erosional features postdating the material constituting them, and emphasize the role of subglacial meltwater flows as a landforming agent.

Harmonized landform regions in the glaciated Central Lowlands, USA

ABSTRACT Many distinct glacial episodes in the past ~1 million years in the Central Lowlands of North America left behind a patchwork of glaciated landscapes of different ages and formed through different glacial, paraglacial, and proglacial processes. Herein, we synthesize and reconcile diverse data sources across nine states in the Central Lowlands to create a generalized landform regions map, incorporating information from surficial geology, age-constrained glacial boundaries, soil properties and parent material, and topography. The resulting map presents regions of till plains, moraines, outwash plains, and glacial lakes with ages ranging from 10.2 cal ka BP to greater than 500 cal ka BP. This new map improves on regional perspectives of the glaciated Central Lowlands by modifying and reconciling boundaries to agree with a variety of spatial data sources with finer detail. This refined map enables a more accurate spatial analysis of the landform regions by reducing noise from imprecise boundaries.

Permafrost Landscape History Shapes Fluvial Chemistry, Ecosystem Carbon Balance, and Potential Trajectories of Future Change

AbstractIntensifying permafrost thaw alters carbon cycling by mobilizing large amounts of terrestrial substrate into aquatic ecosystems. Yet, few studies have measured aquatic carbon fluxes and constrained drivers of ecosystem carbon balance across heterogeneous Arctic landscapes. Here, we characterized hydrochemical and landscape controls on fluvial carbon cycling, quantified fluvial carbon fluxes, and estimated fluvial contributions to ecosystem carbon balance across 33 watersheds in four ecoregions in the continuous permafrost zone of the western Canadian Arctic: unglaciated uplands, ice‐rich moraine, and organic‐rich lowlands and till plains. Major ions, stable isotopes, and carbon speciation and fluxes revealed patterns in carbon cycling across ecoregions defined by terrain relief and accumulation of organics. In previously unglaciated mountainous watersheds, bicarbonate dominated carbon export (70% of total) due to chemical weathering of bedrock. In lowland watersheds, where soil organic carbon stores were largest, lateral transport of dissolved organic carbon (50%) and efflux of biotic CO2 (25%) dominated. In watersheds affected by thaw‐induced mass wasting, erosion of ice‐rich tills enhanced chemical weathering and increased particulate carbon fluxes by two orders of magnitude. From an ecosystem carbon balance perspective, fluvial carbon export in watersheds not affected by thaw‐induced wasting was, on average, equivalent to 6%–16% of estimated net ecosystem exchange (NEE). In watersheds affected by thaw‐induced wasting, fluvial carbon export approached 60% of NEE. Because future intensification of thermokarst activity will amplify fluvial carbon export, determining the fate of carbon across diverse northern landscapes is a priority for constraining trajectories of permafrost region ecosystem carbon balance.

Estimating soil organic carbon levels in cultivated soils from satellite image using parametric and data-driven methods

ABSTRACT Soil organic carbon (SOC) is one of the key soil components for cultivated soils. SOC is regularly monitored and mapped to improve the quality, health, and productivity of the soil. However, traditional SOC-level monitoring is expensive for land managers and farmers. Estimating SOC using satellite imagery provides an easy, efficient, and cost-effective way to monitor surface SOC levels. The objective of this study was to estimate the surface SOC distribution in selected soils of Major Land Resource Areas (MLRA), 102A (Rolling Till Plain, Brookings County, SD), and 103 (Central Iowa and Minnesota Till Prairies, Lac qui Parle County, MN), using satellite imagery with different resolutions (Landsat 8 and PlanetScope). The dominant soils in the study area are Haplustolls, Calciustolls, and Endoaquolls, which are formed in silty sediments, local silty alluvium, and till. Landsat 8 and PlanetScope spectral bands were used to develop SOC prediction models. Parametric and data-driven methods were employed to predict the SOC. Multiple linear regression and Linear Spatial Mixed Model (LSMM) were used on the Landsat 8 and PlanetScope data. In addition to the parametric models, Regression Trees and Random Forest were also employed on both satellite data. The results showed that reduced LSMM provided the lowest RMSE, which are 0.401 and 0.367 for Landsat 8 and PlanetScope, respectively. Furthermore, the random forest has the highest RPD and RPIQ for Landsat 8 (RPD of 2.67 and RPIQ of 2.49) and PlanetScope (RPD of 2.85 and RPIQ of 3.7). In all cases, models obtained from PlanetScope are better than those obtained from Landsat 8.

De-icing landsystem model for the Universidad Glacier (34° S) in the Central Andes of Chile during the past ~660 years

Reconstructing latest Holocene (<1000 years) glacial landscape development in the Central Andes of Chile (30–35° S) is key for understanding the response of the cryosphere during periods of negative glacier mass balance, such as the current one. The excellently preserved glacial landscape produced during the latest ice advance and retreat cycle is of particular interest for examining the detailed response of glaciers to deglaciation. To establish a conceptual model of glacier behavior under warm and dry climatic conditions, we reconstructed and dated the recent glacial history of the Universidad Glacier (34° S) through detailed geomorphological mapping and 10Be cosmogenic surface exposure dating. Our mapping describes a landsystem that spans from the current ice front to ~3 km down-valley, where a mosaic of glacial landforms includes mounded relief; sinkholes; debris-filled stripes; moraine belts; flutings; and a prominent basal till plain. Our 10Be ages suggest that the Universidad Glacier has fluctuated in its forefield since the 13th – 15th centuries CE. We propose that the glacier evolved from a clean glacier to a debris-covered glacier, to an ice-cored moraine, and finally, to a massive dead-ice topography. This deglacial evolution intermittent and potentially reset by multiple standstills and/or re-advances during the overall retreat. The implication is that phases of active ice were followed by stagnation associated with progressive melting of dead ice under the supraglacial debris layer. Similar geomorphic features and processes are recorded in the present-day Universidad Glacier ablation zone, denoting a recurrent reconditioning over time analogous to the glacier's evolution during the latest Holocene.

Exploration of Glacial Landforms by Object-Based Image Analysis and Spectral Parameters of Digital Elevation Model

Glacial landforms are a significant element of landscape in many regions of Earth. The increasing availability of high-resolution digital elevation models (DEMs) provides an opportunity to develop automated methods of glacial landscape exploration and classification. In this study, we aimed to: 1) identify glacial landforms based on high-resolution DEM datasets; 2) determine relevant geomorphometric and spectral parameters and object-based features for the mapping of glacial landforms; and 3) develop an accurate workflow for glacial landform classification based on DEM. The developed methodology included the extraction of secondary features from DEM, feature selection with the Boruta algorithm, object-based image analysis, and random forest supervised classification. We applied the workflow for three study sites: one in Svalbard and two in Poland. It allowed the identification of six categories of glacial landforms: till plains, end moraines, hummocky moraines, outwash/glaciolacustrine plains, valleys, and kettle holes. The majority of relevant secondary features represented DEM spectral parameters calculated from 2-D Fourier analysis. The supervised classification models with the highest performance exhibited up to 96% overall accuracy with regard to a groundtruth dataset. This study showed that glacial landforms can be identified using novel image-processing methodology and spectral parameters of high-resolution DEM. The complete classification workflow developed herein provides a solution for the transparent generation of thematic maps of glacial landforms that may be reproducible and transferrable to various glacial regions worldwide.

Status and short-term environmental changes of lakes in the area of Devil's Bay, Vega Island, Antarctic Peninsula

AbstractClimate change affects various components of the polar environment, including lacustrine systems in many regions of the Antarctic Peninsula. However, the recent status of these life-preserving habitats remains poorly known. We performed field geomorphological mapping and limnological characterization of lakes distributed across Devil's Bay in the northern area of Vega Island, James Ross Archipelago, and described some significant changes in their physical and chemical properties during the 2013 summer season. Fifty lakes were described in the area, even though the area is relatively small (~12 km2). Six major lake types were distinguished based on their geomorphological settings and their origins, and their hydrochemical properties were compared with the neighbouring lake districts. The geomorphic characteristics of the lakes range from stable lakes on till plains (21 out of 50 lakes), to kettle lakes on fresh moraines with an unstable shoreline, to very dynamic ice-contact lakes with fluctuating water levels and inundated areas. The ice-contact lake dynamics were observed using ‘Lake Payer’ as an example. At the interface between hydrology and geosciences, this study provides new and original data from a yet largely unexplored region and emphasizes the importance of cross-disciplinary research in polar sciences.

Longitudinal development of clast shape characteristics from different material sources in Hrbye River, Central Svalbard

The sediment transport in polar regions is highly changeable and it is getting faster in connection with a climate change. This study describes the Hrbye River catchment located in the northern Billefjorden, Central Svalbard. The Czech Arctic Station and AMUPS - Adam Mickiewicz University Polish Polar Station are located in near this locality Petunia Bay. The material for this study was sampled in August 2016, during the summer research campaign of Czech Arctic Station together with a cooperation between Masaryk University in Brno and the University of Oslo via Norway Grants. The catchment area is 60 km2. The area of interest lies around the 10 km long Hrbye River in its braidplain, which is 2.3 km wide and 4.5 km long. In the Hrbye Glacier forefield, 27 sediment sampling localities were selected and defined into seven groups: (i) esker complex; (ii) debris stripes; (iii) till plain; (iv) hummocky moraine; (v) post-LIA braidplain; (vi) LIA moraine; (vii) LIA braidplain. Three main petrological types of rocks were studied (SVP – sandstone, VAP – limestone, ORT – orthogneiss). Lithology and roundness of the clasts were evaluated in order to study clast shape properties from various glacial sediments. The results show the dominant role of lithology on the clast shape modification in the Hrbye Glacier forefield.

Twenty Years of Tallgrass Prairie Reconstruction and Restoration at Pawnee Prairie Natural Area, Missouri

In 1996 the Missouri Department of Conservation purchased Pawnee Prairie, a 190-ha mix of remnant tallgrass prairie and formerly row-cropped prairie with varying degrees of Festuca arundinacea invasion and past cattle grazing intensities on rolling terrain in the central dissected till plains ecological section. Management actions implemented over the following 20 y included prescribed fire, herbicide treatments of invasive nonnative species, and seeding of local ecotype prairie seed. Concurrently, four vegetation monitoring transects were sampled for plant species composition and cover five times between 1996 and 2017. Each of the transects increased significantly over time in the following per-quadrat means: % native plant species cover, plant species conservatism, and cover-weighted plant species conservatism. At the site level, native grasses increased by 22%, nonnative grasses declined by 76%, native forbs increased by 91%, nonnative forbs declined by 94%, and native sedges declined by 37%. In 1996 the top species in importance value across all transects included weedy native species (e.g., Dicanthelium lanuginosum) and nonnative species (e.g., Daucus carota). By 2017 the top species had transitioned to characteristic prairie species (e.g., Schizachyrium scoparium). Ordination results documented compositional trends across all transects toward greater native species richness, cover, and species conservatism values. At Pawnee Prairie, 20 y of sustained prairie reconstruction and restoration practices applied across an area of differing land use histories resulted in significant gains in the natural quality of the site's vegetation, including a greater abundance of prairie flora matrix species and some conservative species.

New State Records and Habitat Characteristics for Kloosia Kruseman (Diptera: Chironomidae) from Wadeable Streams in Missouri

Abstract First records of occurrence of the midge genus Kloosia Kruseman in Missouri are reported based upon aquatic macroinvertebrate community samples collected during April 2012 and October 2015 from reaches of Grindstone Creek in Dekalb County and South Blackbird Creek in Schuyler County in the Central Dissected Till Plains. Select physical and water quality characteristics from the reaches are also provided.

Geomorphological evidences of paleosurge activity in Lake Viedma Lobe, Patagonia, Argentina

The valley of Lake Viedma represents one of the main ice flow axes during the Quaternary glaciations of Patagonia. Geomorphological mapping of the eastern sector of the lake, between the shoreline and the end moraine systems attributed to the Last Glaciation, allowed us to recognize a large number of ice-marginal and subglacial landforms. The landform assemblage of this area matches with previous models of surging glacier landsystems around the world. We identified the three geomorphological zones that characterise surging glacier landsystems, with some mixing of the zones, possibly related to multiple surges. The outer zone includes push moraines, thrust-block moraines, overridden moraines and hummocky moraines. In the middle zone, overridden moraines and hummocky moraines were developed. Finally, the inner zone is composed of megaflutes, crevasse-squeeze ridges, eskers and till eskers in a till plain. We interpreted these geomorphological features as evidence of surge activity of the Lake Viedma Lobe. A combination of topographic setting, soft-bedrock and the presence of a proglacial lake in contact with ice could have contributed to the development of surge activity in the Viedma paleoglacier. It could have been triggered by a change in the glacial thermal regime of the Viedma paleoglacier from polar during the Last Glacial Maximum to polythermal or temperate during the beginning of the Late Glacial.

Subsurface heterogeneity in the geological and hydraulic properties of the hummocky Paris Moraine, Guelph, Ontario

The advance and retreat of ice lobe margins of the Laurentide Ice Sheet formed moraines that are a prevalent feature throughout southwestern Ontario. In contrast to the well-studied stratified moraine complexes, recessional and end moraines have largely been ignored in the context of hydrogeological studies. Recent urban growth has led to development pressures on these moraines and a need to better understand their hydrogeology. This study presents data sets from the Paris Moraine near Guelph, Ontario, to examine its geomorphology, internal composition, and the corresponding hydraulic properties of these ice-marginal features. The moraine’s geomorphic elements were mapped using high-resolution Global Positioning System transects, aerial photograph analysis, and ground truthing. Nine continuous sediment cores were recovered to determine the nature and distribution of subsurface sedimentary units and their relation to the regional stratigraphic framework. Cores were described in detail using standard sedimentological techniques, and significant sediment heterogeneity was observed in cross sections. Grain-size analyses of over 150 samples provide site-specific estimates of saturated hydraulic conductivity. In addition, saturated hydraulic conductivity was measured on 104 samples using the falling head permeameter method. This study found that different scales of sediment heterogeneity occur across the moraine and the associated till plain and outwash. In contrast, the hydraulic conductivity varies much less. It is expected that certain sedimentary units at specific depths will impact groundwater flow at the centimetre to hundreds of metres scale, which is significant in environmental site assessments or for understanding contaminant hydrogeological problems.

Karst sinkhole mapping using GIS and digital terrain models

Danish glacial landscape elements such as basal till plains, hummocky moraine areas and outwash plains contain a variety of small and large depressions. They were probably formed in glacial, late-glacial or Holocene time and may represent dead-ice holes or degraded pingos, or sinkholes formed by interaction between pre-Quaternary chalk or limestone bedrock and the thin glacial co ver. The aim of this study is to map terrain depressions that might potentially be karst sinkholes by analysing digital terrain models in the geographic information system (GIS). The incentive to apply the technique for mapping of sinkholes came from an accidental acquaintance with a farmer, Jens Kirk, whose farmland is located near Thisted. Jens Kirk told us that the front end of his tractor had suddenly sunk into the ground during routine farming work, and this incident was our inspiration to start the project described here.

Orthophosphorus Contributions to Total Phosphorus Concentrations and Loads in Iowa Agricultural Watersheds.

Phosphorus (P) is delivered to streams as episodic particulate P and more continuous soluble P (orthophosphorus [OP]), and it is important to determine the proportion of each P form in river water to more effectively design remedial measures. In this study, we evaluated the annual mean ratios of OP to total P (TP) concentrations and loads in 12 Iowa rivers and found systematic variation in the ratios. The OP/TP ratios were >60% in two tile-drained watersheds of the Des Moines Lobe and in a shallow fractured bedrock watershed in northeast Iowa, whereas in southern and western Iowa, OP contributions to TP were <30%. Higher OP/TP ratios were associated with greater row crop intensity in the watershed and a greater proportion of baseflow in the river. Orthophosphorus contributions from croplands would be greater in watersheds characterized by widespread tile drainage and well-drained soils, whereas cropland TP export would be dominated by particulate P in dissected till plains with poorly drained soils. Understanding the dominant form and transport pathway of P from agricultural areas in a watershed is seen as an important first step in determining appropriate conservation practices to reduce P loads.

Hydro-geomorphic controls of greenhouse gas fluxes in riparian buffers of the White River watershed, IN (USA)

Riparian ecosystems are defined by the nature and regularity of the interactions between a given river system and its floodplains, and past studies have often presented vegetation cover as the exclusive expression of these interactions. There has been to our knowledge, no systematic attempt at linking greenhouse gases (GHG) fluxes and types of riparian buffers. The present study was conducted to investigate the intensity and seasonality of carbon dioxide (CO2), methane (CH4) and nitrous oxide (N2O) fluxes in riparian buffers in three common hydro-geomorphic settings (HGM) across the White River watershed (Indiana, USA). These classes included riparian sites located: (i) in till plain depressions near 1st order streams (HGM-1), (ii) in incised narrow valleys with thin alluvium layers above glacial till (HGM-2), and (iii) along 3rd–4th order streams in broad floodplains with thick alluvial and glacial outwash deposits (HGM-3). For each class, 3 sites were selected and GHG fluxes were measured during the wet (May) and dry seasons (August). Strong relationships were found between GHG fluxes, soil properties and environmental factors, but these relationships varied with season and gas species, making it challenging to rely on these relationships for GHG fluxes upscaling. Analysis of variance and discriminant analysis showed that the HGM-defined riparian buffers were distinct in terms of GHG flux intensity. Regardless of season, the HGM-1 sites emitted CO2 at rates 1.6 times higher than at the other sites, likely due to difference in soil C quality. During the wet season, N2O emission was significantly higher at the HGM-3 than at the other sites (0.88 vs 0.27mgNm−2d−1), and was negatively related with the gradient of the adjacent channel (r2: 0.69). The riparian buffers acted as CH4 sinks, with the HGM-2 sites exhibiting CH4 uptake rates significantly greater than the other riparian types (−0.80 vs −0.34mgCH4-Cm−2d−1). The consistency of these results underscores the potential of an HGM-based monitoring approach to derive watershed-scale GHG budgets for riparian buffers.

Sukcesja osadowa i etapy formowania ozu gostynińskiego, Równina Kutnowska, centralna Polska

Artykuł nie zawiera streszczenia w języku polski.

Geochronological (OSL) and geomorphological investigations at the presumed Frankfurt ice marginal position in northeast Germany

Abstract The Weichselian Frankfurt ice marginal position in northeast Germany has been critically discussed in the past owing to weak morphological evidence and a lack of clear sedimentological records. This study aims to contribute to this discussion with new geochronological and geomorphological results. Apart from very few cosmogenic exposure ages, the time frame is to date still based on long distance correlation with radiocarbon chronologies. We selected a study site in a key position regarding the classic location of the Frankfurt ice marginal position and the recently described arcuate ridge structures on the Barnim plateau. For the first time we present Optically Stimulated Luminescence (OSL) ages of quartz from glaciofluvial deposits for this Weichselian phase. Our results indicate an advance of the Scandinavian Ice Sheet (SIS) at around 34.1 ± 4.6 ka. This is in agreement with OSL ages from sandur deposits at the Brandenburg ice marginal position located farther south and could also be correlated with the Klintholm advance in Denmark. The subsequent meltdown phase lasted until around 26.3 ± 3.7 ka. During the meltdown phase a minor oscillation of the SIS caused the formation of the recently described arcuate ridges on the Barnim till plain. Recalculated surface exposure ages of glacigenic boulders with an updated global production rate indicate a landscape stabilization phase at around 22.7 ± 1.6 ka, which is in agreement with our ages. A phase of strong aeolian activity has been dated with OSL to 1 ± 0.1 ka; this may have been triggered by human activities that are documented in this region for the medieval period.

Miocene deposits at Silkeborg, Jylland, and their influence on hydrology

A motorway was constructed in 2010–2016 through the suburbs of the city of Silkeborg (Fig. 1). The Danish Road Directorate wished to climate-proof the motorway against adverse future climate changes. The directorate collaborated with the Geological Survey of Denmark and Greenland (GEUS) to study the hydrological conditions. Studies of historical and projected climate-change-driven variations in groundwater levels in relation to urbanised hydrological fluxes were conducted by Kidmose et al. (2013, 2015). During the construction of the motorway, Miocene and Quaternary deposits were exposed in the slopes of the Gudenå valley and late-glacial glaciofluvial deposits were found in the valley floor. This paper focuses on the Miocene sediments and their influence on the local hydrological conditions.&#x0D; At Silkeborg the Gudenå valley is c. 35 m deep (Fig. 1). The surrounding terrain is a till plain. In the slope of the valley, glaciofluvial sand is found below the till. Miocene deposits are found below the glaciofluvial sand. The floor of the Gudenå valley is covered by c. 15 m thick glaciofluvial deposits, which rest on Miocene deposits. In borehole no. DGU 87.907 49 m of Miocene deposits belonging to the Vejle Fjord Formation are recorded, consisting primarily of marine clay with minor occurrences of sandy deposits. About 12 km south of Silkeborg lower Miocene deposits are seen in outcrops and boreholes (Fig. 2). Here the fluvial Addit Member of the Billund Formation (Rasmussen et al. 2010) is separated from the underlying marine Vejle Fjord Formation by a sharp erosional contact (Rasmussen 2014).