Transport of nitrate-containing groundwater to coastal areas through buried tunnel valleys, Denmark

Nitrogen impact on the aquatic environment, including coastal areas, is too high in many countries worldwide, particularly in areas with intensive agriculture. Efficient mitigation initiatives demand that important pathways and the fate of nitrate in the hydrological cycle are known. In this study, we focus on groundwater nitrate contamination in two near-shore catchment areas in north-west Denmark. Groundwater in the area is mainly located in buried tunnel valleys, which are subsurface structures eroded by meltwater during Pleistocene glaciations in former glaciated areas. Ground-water samples from the aquifers inside the buried valleys reveal the presence of up to 120 mg/l nitrate down to 10 m below sea level and about 1 km down from the stream outlet towards the coast. We interpret the complex tunnel-valley infill to be responsible for the spatial heterogeneity of the groundwater geochemistry, where sandy geological windows create localised hydraulic pathways and complex redox structures. Groundwater and stream water chemistry in the study area clearly demonstrate the role of groundwater in nitrate transport within the catchment as well as the direct pathway to the coast bypassing the stream and riverine systems. Our results show that the buried tunnel valleys potentially contribute to submarine groundwater discharge and therefore could be responsible for a hitherto unaccounted input of nitrogen to the marine environment.


Introduction
Eutrophication of sea waters is a global phenomenon, due to nitrogen (N) enrichment that mainly stems from inland riverine systems (e.g.DEFRA 2016).Anoxic dead zones, which are considered a key stressor on marine ecosystems, can be the consequence (Diaz & Rosenberg 2008).Nitrogen pollution to coastal zones around the world originates from different terrestrial sources, such as waste water, atmospheric fixation and deposition, livestock manure and fertiliser application (Erisman et al. 2011).The implementation of Danish national action plans since the mid-1980s for the reduction of N losses has successfully reduced the impact of N on the aquatic environment (Dalgaard et al. 2014;Hansen et al. 2017).However, Denmark still has a relatively high level of agricultural N loss, surpassed only by a few areas in Europe, such as in the Netherlands and the Po Valley in Italy (De Vries et al. 2011).
To achieve the goals of the EU water framework directives, losses of N must be reduced, especially to coastal areas.This is true for the case study presented in this paper where mitigation actions are needed to considerably reduce the N impact on Skive Fjord, Denmark, according to the recently published Water Plans for 2021-2027(Danish Environmental Protection Agency 2023).
Nitrate (NO 3 -) in groundwater, surface water, the riparian zone and the sea floor can be reduced or retained (Lutz et al. 2020).Especially the longterm capacity of riparian zones in regulating groundwater nitrate fluxes is www.geusbulletin.orgnot well understood (Hill et al. 2014).Nitrate removal in some riparian zones may be absent or minimal because deeper groundwater flow paths may not interact with nitrate-reducing minerals and organic matter (Hill 2019).However, N flux to the marine environment might not be limited to output via streams and rivers.Several studies suggest that transport of N to the marine environment via submarine groundwater discharge (SGD) can be substantial (e.g.Bishop et al. 2017;Duque et al. 2019;Taniguchi et al. 2019;Szymczycha et al. 2020), and its contribution to the local N budget may be much greater in geologically complex settings where fractures or preferential flow pathways can create hotspots or point sources of N fluxes (Santos et al. 2021).A quantification of the amount of SGD under such geological settings, therefore, is challenging because in addition to direct measurements of groundwater flow, it requires detailed geological and hydrological models for both onshore and offshore (e.g.Haider et al. 2015).Likewise, assessments of the subsurface redox conditions are needed to determine where nitrate reduction takes place (Hansen et al. 2021;Kim et al. 2021).Therefore, modelling of geology, hydrology and geochemistry, along with calculations of the magnitude of the SGD and the amount of N transported to the coast is needed to evaluate the impact on the coastal environment.
Because the landscape and the near-surface geology of Denmark are highly affected by the ice sheets and their meltwater during the Pleistocene (e.g.Houmark-Nielsen 2011), the complex geological setting generally calls for advanced hydrogeological and redox modelling in order to map the groundwater flow paths and the fate of N (Kim et al. 2019(Kim et al. , 2021)).Buried tunnel valleys are structures that have a high impact on the geological architecture of the subsurface (Jørgensen & Sandersen 2006;Sandersen & Jørgensen 2017).Most of these structures were formed as tunnel valleys eroded by high-pressure meltwater underneath the Pleistocene glaciers.After formation, the valleys were filled with predominantly glacial tills and meltwater sediments (Sandersen & Jørgensen 2003).Being up to 400 m deep and on average 1-1.5 km wide, buried valleys have been found to significantly impact the groundwater flow (e.g.Sandersen & Jørgensen 2003;Andersen et al. 2013).Several valleys are found close to the coastline (see Sandersen & Jørgensen 2022) and a few have been interpreted to cross the coastline (e.g.Haider et al. 2015;Andersen et al. 2016).In one example, it has been shown that groundwater in inferred buried onshore channels connects with shoreface sediments carrying nitrate to the coastal zone (Andersen et al. 2007).In combination, these studies show a potential for buried valleys in enabling or enhancing preferential flow of groundwater across the coast, bypassing nitrate reduction zones in, for example, the riparian zones and thus transporting nitrate directly to the marine environment.
In this paper, we present results from a multidisciplinary study on nitrate transport in the uppermost 40 m of the subsurface in a near-shore area on the Salling peninsula, north-west Denmark.Here, large and deep as well as narrow and shallow buried tunnel valleys close to the shore are mapped in high detail.Our study hypothesises that these mainly sand-filled and gravel-filled tunnel valleys play an important role as preferential pathways of nitrate to the coastal ecosystem, thereby increasing the risk of eutrophication.

Topography
The study area comprises the two adjacent catchments of the streams Hulebro Baek (study site 1; 11.3 km 2 ) and Hagens Møllebaek (study site 2; 27.6 km 2 ) located on the Salling peninsula in the north-western part of Jylland (Fig. 1).The topography is highest to the north and in the south-west with elevations up to 40-50 m above sea level (a.s.l.), sloping down to below 20 m a.s.l. in the central and southern parts of the study area (Fig. 1a).The two streams cut through the terrain in an overall south-easterly direction towards Skive Fjord (Fig. 1) -a fjord that is generally 3 to 4 km deep, but up to 5 m deep in places and that constitutes a southern branch of the larger Limfjorden farther north.Limfjorden has connections to the Kattegat to the east and the North Sea to the west (Fig. 1).

Geological setting
The pre-Quaternary sedimentary succession relevant for this study consists of Palaeogene clays with thin layers of Miocene mica sand, overlain by silt (Rasmussen et al. 2010).The overlying Quaternary succession varies in thickness from less than 10 m to more than 120 m and is dominated by tills and meltwater deposits.The study area is characterised by a clay-dominated moraine landscape from the Weichselian glaciation with occurrences of postglacial freshwater deposits in low-lying areas.Several ice advances have previously covered the area, resulting in a complex Quaternary succession (Houmark-Nielsen 2011).
The valley system previously mapped in the study area comprises two valleys with a N-S-and NNE-SSW-orientation and an older generation of valleys with a WNW-ESE-orientation (Fig. 2a).The valleys are mainly eroded into Palaeogene clay-dominated sediments (down to 100 to 150 m below sea level [b.s.l.]).According to borehole data, the buried valleys are www.geusbulletin.orgfilled with successions of Quaternary clay tills, meltwater sand and meltwater clays.All valleys are now completely buried and not directly recognisable in the present-day landscape (Sandersen & Jørgensen 2016b).

Climate and agriculture
The agricultural area in Denmark spans about 61% (2022) of the country's total area of 42 951 km 2 .The climate is costal temperate with a mean temperature of 9.1°C, and an average precipitation of 781.7 mm per year as reported by the Danish Meteorological Institute for the years 2011-2020(DMI 2023)).The study area is dominated by agricultural fields (c.88%), whilst buildings (c.10%) and nature (c.2%) occupy smaller parts of the area (Levin et al. 2014).Both sandy and clay-rich soil types are found within the study area (Adhikari et al. 2013) with a relatively low content of organic matter (<6%; Adhikari et al. 2014).The crop types are mainly cereals and to a lesser extent permanent grass, and dairy, pig and arable farming types are abundant (Danish Agricultural Agency 2018).

Drillings and lithological descriptions
In this study, a drilling campaign operated by the Ejlskov company was carried out to collect lithology and water chemistry samples at seven locations in November 2021.The lithology samples were collected using a Geoprobe direct push method (DT-22), and the groundwater samples were collected using a modified Geoprobe  (Sandersen & Jørgensen 2016a, b, 2017).(b) Buried valleys mapped using the tTEM method.Interpreted from tTEM 3D grid slices, with mean resistivity slices at 2 m a.s.l.(see resistivity scale in figure legend).(c) A south-north-orientated cross-section through the 3D tTEM resistivity grid superimposed with geological interpretations from boreholes.BV: buried valley.

Geophysical surveys
Two types of geophysical TEM (Transient Electromagnetic) data were used in this study: (1) tTEM (towed TEM) soundings data (Auken et al. 2019; Aarhus University 2021) with a line spacing of 25 to 60 m shown in black on Fig. 1d, and (2) TEM40 (ground-based TEM) soundings data (Christiansen et al. 2006).The TEM40 data are not included in Fig. 1d.The TEM method uses a transmitter coil to emit a current, which gives rise to a magnetic field.The current is then rapidly turned off, resulting in a decrease of the magnetic field, which, in turn, gives rise to a secondary magnetic field decaying over time.The rate of decay is measured in a receiver coil, which provides the resistivity of the subsurface.Clay deposits are characterised by slow decays due to the conductive nature of the material, whereas sand or gravel layers, which are characterised by a high resistivity, will have fast decays.
The two systems have different spatial resolutions and depths of investigation: tTEM provides high-resolution images of the uppermost 70 m (Auken et al. 2019), whereas the depths of investigation for the TEM40 system range from 200 to 300 m.As a result of this and in combination with the generally denser data coverage, highly detailed mapping of the near-surface geology is possible (Sandersen et al. 2021).All TEM systems are based on the time-domain electromagnetic method and measurement of the conductivity of the subsurface layers (Danielsen et al. 2003).The method therefore provides information of the resistivity distribution of the sedimentary succession (Sørensen & Auken 2004).

Stream discharge and water chemistry
The stream discharge and nitrate concentrations were extracted from the public database ODA (https://odaforalle.au.dk).For streams, water chemistry monitoring focuses only on nutrient compounds, such as nitrate, nitrite and phosphate.In both catchments, discharge and stream chemistry were monitored at the outlet of the catchment (Fig. 1a).In Hulebro Baek, stream discharge has been measured daily since 1989, and nitrate concentration has been monitored at monthly (dry season) to bi-weekly (wet season) intervals since 1984.In Hagens Møllebaek, discharge and stream water chemistry were monitored since 2020 (daily) and 2017 (bi-weekly to monthly), respectively.

Groundwater chemistry
The groundwater chemistry data were extracted from the Jupiter database.In this study, the focus is on major cations (Ca 2+ , Mg 2+ and Na + ), anions (Cl -, NO 3 -and SO 4 2-), pH and the groundwater alkalinity.Groundwater chemistry information was available from 48 boreholes (Fig. 1c) within a depth range from 3 to 157 m b.l.s.(below land surface).The groundwater chemistry in these boreholes has been monitored periodically since 1951, and the data availability varied from 1 to 5 measurements per borehole.This study includes all the available measurements.Except for nitrate, concentrations of major cations and anions are expressed in mM.Stoichiometric ratios of these elements are used to infer underlying processes.Nitrate, however, is shown in mg/l because the standard for groundwater and drinking water is 50 mg/l of nitrate according to Danish and European legislation (vs.0.8 mM), and nitrate concentrations below 1 mg/l (vs.0.016 mM) are considered reduced conditions.Thus, it is more intuitive to express nitrate in mg/l.

Visualisation and interpretation of data
The tTEM soundings are visualised with low resistivities being blue to green and high resistivities being red to purple (see Fig. 2).The DOI (depth of investigation) shown in Fig. 2c illustrates the maximum depth at which the tTEM models can be considered reliable (Christiansen & Auken 2012).Based on the tTEM data, a 3D grid and a series of horizontal 2 m mean resistivity slices were made.The map in Fig. 2b shows a selected resistivity slice at 2 m a.s.l., and Fig. 2c shows a selected vertical cross-section through the 3D resistivity grid and selected boreholes.The TEM40 data were used as single soundings to support the geological interpretations of the tTEM data and are not visualised on the profile.
The borehole data and the geophysical data were imported into the 3D geological modelling tool Geo-Scene3D, enabling joint interpretations of the geological and geophysical data within the 3D model space (see Figs 2b and c).

The pre-Quaternary succession
The tTEM mapping has confirmed the widespread presence of low resistive sediments with resistivities below 10-12 ohm-metre, in boreholes described as Palaeogene clay.On top of this clay, sporadic sandy and clayey Neogene sediments are found.The elevation of the top of the pre-Quaternary sediments varies from 20 m a.s.l.

The Quaternary succession and buried valleys
The tTEM survey has added more detail to the uppermost 50 to 70 m of the subsurface and has especially enhanced the picture of the buried valleys in the area (Sandersen & Kallesøe 2021).The deep and broad buried valleys that were earlier mapped primarily with TEM40 and boreholes (see Fig. 2a) were confirmed with the tTEM survey.In addition to these valleys, several narrow and shallower valleys are mapped with tTEM (Figs 2b and c).
A large N-S-orientated buried valley is seen on the horizontal mean resistivity slice 2 m a.s.l. in Fig. 2b (labelled BV-E) as generally high resistivities (mainly red to yellow colours) surrounded by low resistivities (blue colours).The valley is clearly seen to the right on the cross-section in Fig. 2c, where sandy and clayey Quaternary sediments represent the infill of the c. 2 km wide structure eroded to depths below 100 m b.s.l.The valley is surrounded by the Palaeogene clays showing very low resistivities.The presence of the valley below the tTEM DOI is confirmed by TEM40 and boreholes.The data show that the deep parts of the valley (below c. 10 m b.s.l.) are divided into two parts, one dominated by clay and the other by sand, represented by low and high resistivities, respectively (Fig. 2c).According to borehole data, the sediments consist of meltwater sands and clays.Between c. 10 m b.s.l. and c. 30 m a.s.l., the tTEM shows a 10-20 m thick low-resistive layer.The clay, generally interpreted as meltwater clay in borehole data, has a wavy appearance and is apparently confined to the valley.Above the clay, there is a succession of meltwater sand or gravel irregularly covered with clayey till.According to previous investigations, the buried valley structure continues southwards to the coast of Skive Fjord (Fig. 2a).
A WNW-ESE-orientated buried valley in the central part of the study area is also around 2 km wide eroded no deeper than c. 60 m b.s.l.(labelled BV-C in Fig. 2).The infill is dominated by high resistivities representing meltwater sand and gravel.Above sea level, clay layers are found occasionally.Along sections of the northern valley flank, a 500-700 m wide high-resistivity structure can be observed from the terrain surface and down to c. 10 m b.s.l.The structure is split into a NW-SE segment and a WSW-ENE segment (BV-D in Fig. 2b).At elevations higher than the resistivity slice shown in Fig. 2b, another narrow high-resistivity structure is present above and with the same orientation as the deep, broad valley.Both valleys are labelled BV-D on the cross-section in Fig. 2c.According to previous investigations, the buried valley structure continues eastwards to the coast of Skive Fjord (Fig. 2a).
To the south, two narrow high-resistive structures with an overall WSW-ENE orientation can be seen (Fig. 2b).The northern of the two, BV-B, is 200-400 m wide, eroded down to 20 m b.s.l. and generally filled with high resistivity layers corresponding to sand.The deepest parts of the valley are eroded down to 50-60 m b.s.l.into the low-resistive Palaeogene clays.The southern valley, BV-A (Fig. 2b and c), appears to be two to three separate erosions, just 200 to 300 m wide, which merge eastwards to create a single valley trace.The tTEM soundings map the southernmost valleys (BV-A) all the way to the coastline of Skive Fjord, whereas the valley BV-B, just to the north, is mapped as close as 800 m from the coast.
Most of the buried valleys in Denmark are formed as buried valleys eroded by meltwater flow underneath the Pleistocene ice sheets (Jørgensen & Sandersen 2006).The valleys mapped in this study area are no exception.The tTEM mapping has revealed an additional system of younger, narrower and shallower buried valleys in the upper part of the subsurface (Fig. 2b; valleys BV-A, BV-B and BV-D).As seen on the cross-section in Fig. 2c, these valleys are filled with high-resistive sediments dominated by meltwater sand and gravel, which are only sporadically covered by low-resistive clayey sediments.
All of the buried valleys in the area are several kilometres long and have been mapped very close to the shore of Skive Fjord to the east.

Nitrate concentration-discharge relationships
The nitrate concentration-discharge (C-Q) relationships of both catchments show positive slopes up to a specific discharge of 1 mm/day and then negative slopes at higher discharge, likely due to dilution at extreme rainfall events (Fig. 3).However, the nitrate concentrations at baseflow were much higher in Hulebro Baek (c.2.0 mg/l) than in Hagens Møllebaek (c. 5 mg/l).In addition, the positive slope of the C-Q relationship was much steeper in Hulebro Baek compared to that of Hagens Møllebaek (Fig. 3).At Hulebro Baek, the long-time series of nitrate data (1984-2020) reveals a clear reduction in stream nitrate concentrations, likely due to N regulation of agriculture since the 1980s.

Groundwater chemistry
The boreholes holding groundwater chemistry data were divided into groups representing the individual buried tunnel valleys (Fig. 4).In all the valleys, groundwater is reduced (i.e.nitrate concentration <1 mg/l) at www.geusbulletin.orgaround 30 m b.l.s., roughly corresponding to sea level, except for groundwater in BV-C.At shallow depths (<30 m b.l.s.), nitrate concentrations in the BV-E valley showed the largest variability (0-100 mg/l) compared to concentrations elsewhere (Fig. 4).Groundwater nitrate in BV-C and BV-A increased with increasing depth (Fig. 4).For instance, nitrate concentrations in groundwater from valley BV-C displayed 120 mg/l at around 30 m b.l.s. or at sea level, which is the highest concentration amongst all the groundwater chemistry measurements.
Groundwater from valley BV-D showed moderately high nitrate concentrations, which decreased with increasing depth.

Nitrate pathways from land to coast via streams
The C-Q relationships provide valuable information about stream water and solute transport at the catchment scale (Godsey et al. 2009;Kim et al. 2017Kim et al. , 2021;;Bieroza et al. 2018).The differing C-Q relationships in the two catchments indicate that nitrate is transported via different pathways from the agriculture fields to the stream outlets.The lower nitrate concentration at baseflow and the steeper slope of C-Q in Hagens Møllebaek compared to those of Hulebro Baek can be interpreted as a limited contribution of groundwater to outlet (Fig. 3).
During baseflow, the stream Hagens Møllebaek may be fed by the groundwater in the riparian zone; thus, it likely contains very low nitrate.The steep slope of the C-Q relationships of Hagens Møllebaek implies that nitrate is primarily transported via shallow subsurface pathways such as drains, interflow or overland flow in this system.
The surface geology of the Hagens Møllebaek catchment, which is predominantly clayey till, supports this interpretation (Fig. 1b).
In Hulebro Baek, groundwater appears to be significant in transporting nitrate to the stream.Valley BV-C intersects with the Hulebro Baek outlet (Fig. 2), and the groundwater in this buried valley is oxic with nitrate concentrations exceeding 100 mg/l (Fig. 4).This oxic, nitrate-containing groundwater may discharge into the stream year-round, resulting in higher nitrate concentration during baseflow and also the gentle slope of the C-Q relationship shown in Fig. 3.In addition, this catchment is covered with sandy till, and the subsurface contains meltwater sand and gravel, which explains the deeper percolation of nitrate containing groundwater directly discharging to the stream (Figs 1b and 2).Thus, in this catchment, shallow subsurface flow appears to play a minor role in delivering nitrate to the stream.

Nitrate pathways from land to coast bypassing the stream outlets: the role of buried tunnel valleys
The groundwater and stream chemistry observations of Hulebro Baek clearly demonstrate an important role of groundwater for nitrate transport within the catchment.This implies that groundwater may be a direct pathway to the near-coast waters, bypassing the stream as well.
Along valley BV-C, for instance, groundwater containing up to 120 mg/l of nitrate was found far from any stream outlet and very near the coast, and this groundwater is present 2-10 m b.s.l.The biogeochemical evolution of nitrate discharging via groundwater to the sea floor is spatially complex and temporally dynamic (e.g.Kroeger & Charette 2008); however, various studies have documented that SGD should be part of the nitrate budget of the coastal ecosystems (Andersen et al. 2007;Luijendijk et al. 2020;Brookfield et al. 2021).
The spatial heterogeneity of SGD is one of the biggest challenges in quantifying its contribution to biogeochemistry of coastal ecosystems (Kroeger & Charette 2008;Luijendijk et al. 2020).In this study, we found a high degree of spatial heterogeneity.In valleys BV-C and BV-D, groundwater chemistry showed a signature of pyrite oxidation by oxygen as in Eq. ( 1), coupled with carbonate dissolution by sulfuric acid as in Eq. ( 2).This results in SO 4 2-/(Ca 2+ + Mg 2+ ) of c. 0.5 and HCO 3 -/( Ca 2+ + Mg 2+ ) of c. 1 (Fig. 6; Appelo & Postma 2005): The pyrite content, however, may be insufficient to reduce nitrate.These two buried valleys are primarily filled with sand without protective clay layers (Fig. 2), and therefore, it is possible that the infill of these buried valleys is more weathered and thus depleted of reduced compounds, that is, pyrite or reactive organic matter, than those of the other buried valleys.
In valley BV-E, which is the largest and deepest buried valley in the study catchment, the groundwater chemistry varies widely.The nitrate containing groundwater of BV-E shows similar signatures as those of valleys BV-C and BV-D, except for two cases of high nitrate concentrations (SO 4 2-/(Ca 2+ + Mg 2+ ) of c. 0.1), which probably represent young oxic groundwater with a relatively low sulphate content due to a decrease in atmospheric deposition and absence of pyrite oxidation (Fig. 5).We conclude that the upper oxic part of valley BV-E does not contain reduced compounds reactive enough to remove all the oxygen and initiate nitrate reduction.
However, many, mostly deeper, groundwater samples were nitrate-free in valley BV-E (Fig. 4).This groundwater shows a high variability: SO 4 2-/(Ca 2+ + Mg 2+ ) ranges between 0.1 and 0.5, and HCO 3 -/(Ca 2+ + Mg 2+ ) ranges between 1 and 3 (Fig. 6).Complete pyrite oxidation by nitrate reduction is defined in Eq. ( 3).Followed by carbonate dissolution, this will result in SO These stoichiometric ratios were not observed in our catchments; thus, we conclude that pyrite oxidation may not be responsible for the nitrate-free conditions.Instead, we attribute these conditions to the complex geological structure of valley BV-E.In BV-E, the clay layers protect the groundwater from nitrate contamination, but because they are deformed, sandy geological windows most likely result in localised hydraulic pathways resulting in complex redox structures.The nitrate-free groundwater in valley BV-E may be well protected by the clay layers, whilst the oxic groundwater in BV-E may be related to sandy geological windows.The increase of HCO 3-/(Ca 2+ + Mg 2+ ) above 1, which is positively correlated to Na + /Cl -, may support this hypothesis (Fig. 6).We attribute this positive correlation to cation exchange between Ca 2+ in groundwater and Na + in the salt-saturated marine sediments.The sediments of the protected areas may not be in equilibrium with the fresh groundwater because they have not been exposed to influx of enough fresh groundwater.
Especially in areas where buried valleys are eroded into high-lying, clay-dominated substrata and filled with coarse-grained sediments, groundwater flow will predominantly take place within the buried valleys.Several deep and shallow buried valleys are mapped within the study area.Infill ranges from disturbed layers of heterogeneous sandy, gravelly and clayey layers to an apparently homogeneous infill of sand and gravel (Fig. 2).Because the Palaeogene clays outside the valleys can be found very close to the terrain surface, the groundwater flow will be focused in the hydraulically transmissive infill of the buried valleys.To the east in the study area where the valleys are very close to the coast, the groundwater flow will be towards the coast within the buried valleys.The focused flow results in high local fluxes implying that the groundwater emanates through the seabed with relatively high velocities so that nitrate reduction in the seabed may not be fast enough to remove all nitrate before it enters the sea.If the valleys terminate at the coast, for instance, due to later erosion, groundwater discharge will take place along the shore.However, if the valleys continue underneath the seabed, SGD may take place at some point offshore, depending on the transmissivity of the seabed.Due to a lack of offshore data, the extension of the buried valleys underneath Skive Fjord is unknown.However, the onshore geophysical mapping and the general knowledge of buried valleys in Denmark make it highly likely that the valleys have continuations under the fjord.
For SGD to occur, the valley aquifers need to be 'leaky'.Where, on the other hand, deep groundwater www.geusbulletin.orgwithin the buried valley has a clayey seal that prevents salinisation, fresh groundwater within the valley infill could potentially form an overlooked resource of offshore groundwater (referred to as offshore fresh groundwater, OFG; Micallef et al. 2020;Zamrsky et al. 2022).The map in Fig. 7 shows the current knowledge of the location of buried valleys in Denmark based on interpretations of TEM data primarily in combination with borehole data (Sandersen & Jørgensen 2022).Although only c. 50% of the onshore area is mapped with TEM, Fig. 7 shows many valleys close to the coast.Based on comparisons between onshore and offshore tunnel valleys in the Danish North Sea, Andersen et al. (2012) concluded that valley morphologies on sides of the shoreline were comparable, and as today's coastline does not constitute a natural limit for the extension of the buried valleys, it is highly likely that many of the Danish onshore valleys extend offshore.
Buried valleys in Denmark are infilled with mainly tills and meltwater sediments.Rough calculations based on borehole lithologies show that around 40-75% of the infill is composed of meltwater sand and gravel, with the remaining part of the infill consisting of silty and clayey meltwater sediments and tills (Sandersen & Jørgensen 2003).Because the tunnel valleys are eroded deeply into the subsurface and generally have heterogeneous infill, the impact on the architecture of the subsurface is generally significant (Sandersen & Jørgensen 2017).The coarse-grained meltwater sediments often constitute important, deep-seated aquifers, but due to the erosive nature of the valleys and the typically heterogenic infill, the groundwater resources are generally vulnerable to contaminants from the surface (Sandersen & Jørgensen 2003).

Structural analogies and perspectives
The buried valleys represent geological structures that are expected to expand beyond the coastline, and the aquifers contained within them therefore constitute a part of the present-day hydrological cycle.This implies a continuous flow of fresh groundwater and a flux of dissolved compounds to coastal areas through the valley sediments.Potential flow of nitrate-containing groundwater directly to the coastline and the sea floor, by-passing the riparian zone, is currently not considered in Danish nitrogen budgets in the water plans for the protection of the aquatic environment.Nitrate is the major form of N in groundwater.However, it might be converted to reduced forms when discharging through the seabed depending on the balance between flow rates and nitrate reduction rates (Andersen et al. 2007).Measurement of the total N content in the bottom water  column below the halocline of the inner fjord has shown that this could be a considerable source of N for biological production creating oxygen deficiency in the summer period (Gertz et al. 2018).Direct contributions of nitrate to the near-coastal environment via groundwater could be significant and might therefore be important for the choice of the right mitigation measures in agricultural areas to reduce the N impact on the marine environment.
Buried tunnel valleys are common both onshore and offshore in former glaciated areas (van der Vegt et al. 2012).Since many are found close to the coastline, their inferred continuations are obvious targets for both SGD and OFG investigations.Other geological structures, which in equal ways may enable or enhance groundwater flow across the coastline, could be fault zones, graben structures and offshore continuations of large aquifers in porous rocks and sediments (e.g.Edmunds et al. 2001;Hinsby et al. 2001;Cantarero et al. 2019;Santos et al. 2021).Outside Denmark, the number of investiga tions of both SGD and OFG is increasing (see Taniguchi et al. 2019;Micallef et al. 2020).Pertaining to the conditions in Denmark, widespread near-surface chalk and limestone aquifers could also be potential targets for further investigations.In relation to OFG, the buried tunnel valleys potentially represent a renewable groundwater resource that has not previously been considered.This contrasts with the potential OFG resources described by Zamrsky et al. (2022) who regard most fresh offshore groundwater as likely non-renewable resources.This is because, at a global scale, they are mostly located far from the shore and thus interpreted as being detached from inland aquifers.

Conclusions
The near-surface geology in the study area is dominated by buried tunnel valleys.New detailed geological mapping has revealed near-surface valleys that are narrower and less deeply eroded in comparison to the previously known valleys.All buried valleys are found very close to the shore of Skive Fjord, indicating that the valleys most likely extend offshore.
Groundwater and stream water chemistry observations in the study area clearly demonstrate the important role of groundwater in nitrate transport within the catchment, and that groundwater in buried valleys may constitute a direct pathway to the coast, bypassing the stream and riverine systems.Groundwater samples from the valley infills revealed that nitrate is present below the sea level within the buried valley infill.Samples with as high as 50-120 mg/l of nitrate were found far from the stream outlet down to 10 m b.s.l.close to the coastline.The complex geological structure of the buried tunnel valleys appears to be responsible for the spatial heterogeneity of the groundwater geochemistry in this area, where sandy geological windows result in localised hydraulic pathways, resulting in complex redox structures.In addition to this, the buried valleys generally control the groundwater flow towards the coast.
Consequently, there is a high risk that N will discharge to the coastal waters of Limfjorden via buried valleys.As such, this flux will represent a N contribution to the environment that is usually not included in estimates of N transported to the inner Danish marine environment.Based on our findings, we therefore suggest that direct groundwater discharge to the coast is included in the N management plans in areas such as these and in comparable geological settings.We also stress the importance of producing subsurface models, which include detailed delineation of geology and geochemistry.Otherwise, calculations of the flux of N to the marine environment will not encompass all sources.
The number of investigations of both SGD and OFG is rising internationally, and with the large number of buried tunnel valleys found in Denmark, there appears to be a potential for further studies, here.To calculate and evaluate actual N fluxes to the shallow marine environment, joint investigations should include detailed onshore and offshore mapping of buried tunnel valleys, groundwater modelling and geochemical investigations of both groundwater and sea water.
Apart from the buried tunnel valleys in former glaciated areas, fault zones, graben structures and offshore continuations of porous sediments with high flow rates and presumably low nitrate reduction rates would also be potential targets for further investigations.
www.geusbulletin.orgdata, and PS performed joint interpretations of geophysical and geological data.PS wrote the first draft, and all co-authors contributed with text and/or comments to subsequent versions of the manuscript.
(shorter screen, i.e. 0.25 m)   and a peristaltic pump.Older borehole data have been extracted from the freely accessible National Borehole Archive Jupiter hosted by GEUS(GEUS 2022; Hansen & Pjetursson 2011).

Fig. 4
Fig. 4 Nitrate concentration in groundwater.On the map (a), the average of all existing data of each borehole is shown in colour.Nitrate concentrations of groundwater are shown as a function of sampling depth (b) and elevation (c).

Fig. 3
Fig. 3 Nitrate concentration-discharge (C-Q) relationships of the study catchments.The sampling year of the stream chemistry is shown in colour.

Fig. 5
Fig. 5 Major cation and anion concentration of groundwater in buried tunnel valleys.Valleys BV-A, BV-C, BV-D and BV-E are shown as triangles, circles, diamond, and squares, respectively.The colours of the markers represent nitrate concentrations.