Detecting, quantifying and understanding groundwater discharge to streams
are crucial for the assessment of water, nutrient and contaminant
exchange at the groundwater–surface water interface. In lowland agricultural
catchments with significant groundwater discharge this is of particular
importance because of the risk of excess leaching of nutrients to streams.
Here we aim to combine hydraulic and tracer methods from point-to-catchment
scale to assess the temporal and spatial variability of groundwater
discharge in a lowland, groundwater gaining stream in Denmark. At the point-scale, groundwater fluxes to the stream were quantified based on vertical
streambed temperature profiles (VTPs). At the reach scale (0.15–2 km), the
spatial distribution of zones of focused groundwater discharge was
investigated by the use of distributed temperature sensing (DTS).
Groundwater discharge to the stream was quantified using differential
gauging with an acoustic Doppler current profiler (ADCP). At the catchment
scale (26–114 km
Groundwater and surface-water exchange dynamics are of great importance for a broad range of disciplines within the field of hydrology. For instance, groundwater discharge to streams governs the transfer of solutes and nutrients between sub-surface and surface-water environments (Boulton et al., 2010; Dahl et al., 2007; Gooseff, 2010; Kasahara and Hill, 2008; Krause et al., 2008) and is also a key parameter in controlling stream biodiversity (Malcolm et al., 2003; Hayashi and Rosenberry, 2002). Zones of groundwater recharge and discharge are particularly important in lowland groundwater-dominated streams as many lowland areas are intensively used for agriculture, which significantly increases the risk of transport of nutrients and pollutants to streams with the potential of severe consequences for stream ecology (Hoffmann and Baattrup-Pedersen, 2007; Kronvang et al., 2005); this is a pressing issue for instance in relation to lowering nutrient loads to rivers, lakes and seas (Danish Ministry of Environment, 2011; Griffith et al., 2006).
Controlled by a range of complex temporal and spatial processes governed by topography, catchment geology, hydrology and hydrometeorology (Brunke and Gonser, 1997; Winter, 1999), the exchange between groundwater and surface water is often spatially and temporally highly variable. Thus, the detection and quantification of groundwater–surface water dynamics present a challenge, particularly in lowland streams. In these streams the diffuse groundwater discharge along the stream channel reduces the sensitivity of thermal methods (Lowry et al., 2007; Krause et al., 2012), as well as tracer methods (Gonzales et al., 2009), and can cause low net increase in stream flow which also limits the available methods for detecting groundwater discharge (Briggs et al., 2011). At the same time due to the presence of focused, significant discharge zones (Lowry et al., 2007; Matheswaran et al., 2012) the spatial variability of groundwater discharge can be large (Krause et al., 2012). Therefore, there is a need to improve our understanding of these processes in lowland catchments across spatial scales in order to develop new approaches and simple tools to map and quantify them.
Different studies applying a range of hydraulic and tracer approaches have been summarised by, for instance, Kalbus et al. (2006). Groundwater fluxes at specific point locations have been measured in rivers by use of seepage metres, revealing large heterogeneity within metres (Landon et al., 2001; Langhoff et al., 2006; Rosenberry, 2008). This spatial heterogeneity has been confirmed by use of temperature as a natural tracer (Conant, 2004), where punctual vertical groundwater fluxes have been estimated from vertical sediment temperature profiles (VTPs) using the steady-state analytical solution to the 1-D conduction–convection equation (Schmidt et al., 2007; Jensen and Engesgaard, 2011) and streambed temperature time series (Hatch et al., 2006; Keery et al., 2007).
At the reach scale more integrated measures, such as differential flow gauging (McCallum et al., 2012; Briggs et al., 2011), have been applied to quantify net differences in stream discharge caused by groundwater recharge and discharge. The use of this method, however, is limited by the measurement uncertainty which prevents it from being applied for detecting small changes in groundwater discharge (Briggs et al., 2011). However, recent advances of acoustic Doppler current profiler (ADCP) instruments for stream discharge measurements open up new possibilities for a more detailed detection of net groundwater discharge with short measurement periods and with a high precision (Mueller and Wagner, 2009). Furthermore, distributed temperature sensing (DTS) has become a widely used method for reach-scale detection of groundwater discharge to streams by monitoring temperatures at the sediment–water interface along a fiber optic cable of several km length (Selker et al., 2006a; Tyler et al., 2009). Thereby, groundwater–surface water interactions can be detected over longer stream sections bridging the monitoring gap between point flux estimates and more integrated measures of net differences in stream discharge at the reach scale (Lowry et al., 2007).
However, most field studies presenting measurements of groundwater–surface water dynamics are carried out in stream sections of a few hundred metres (Conant, 2004; Lowry et al., 2007; Anibas et al., 2011), primarily due to the labour intensive work needed to extend measurements to quantify discharge fluxes beyond the kilometre scale. To obtain information about runoff sources at the catchment scale a common approach is stream hydrograph separations, often conducted by use of stable isotopes and chemical tracers (Sklash and Farvolden, 1979; Uhlenbrook and Hoeg, 2003). However, such hydrograph separations only reveal the integrated catchment response of the point-to-reach-scale groundwater–surface water dynamics, and have rarely been conducted in lowland agricultural catchments. The reason for this is that it can be difficult to clearly identify end members due to the damping of signals by the often constant, strong groundwater influence (Gonzales et al., 2009).
Since the large heterogeneity in groundwater–surface water interactions can be observed across scales, the necessity of combining the different hydraulic and tracer methods is widely recognised (Bencala et al., 2011; Kalbus et al., 2006; Lischeid, 2008; Scanlon et al., 2002) in order to avoid incorrect inferences regarding exchange processes based on observations at one spatial scale only (Schmadel et al., 2014). Hence, more recently point-to-reach-scale groundwater–surface water interactions have been studied by applying multiple methods covering different spatial scales such as groundwater head gradients and DTS (Krause et al., 2012); differential flow gauging, chemical tracers and DTS (Briggs et al., 2011); or chemical tracers and differential flow gauging (McCallum et al., 2012). However, either the studies did not detect small-scale spatial variability in groundwater discharge (Briggs et al., 2011; McCallum et al., 2012) or did not quantify discharge fluxes at the identified discharge zones (Krause et al., 2012). Furthermore, to our knowledge no study has so far combined point-to-reach-scale DTS, VTPs and differential gauging with catchment-scale tracer-based hydrograph separations.
Map of the study area and sampling sites.
The aim of this study was to combine hydraulic methods (ADCP, groundwater
head gradients) and tracer methods (hydrograph separations from electrical conductivity (EC) and
The study was carried out in the groundwater gaining lowland Holtum stream,
located in the Skjern river catchment in Jutland, western Denmark (Fig. 1a).
This glacial floodplain valley is characterised by thick sediment deposits
of sand and silt deposited during the latest Weichsel glacial period
(Houmark-Nielsen, 1989), and with podzols being the dominating soil layers. The
mean annual precipitation in the catchment is 950–1000 mm with an actual
evapotranspiration of 460–480 mm yr
The study catchment at Holtum stream was divided into four sub-catchments, and at each sub-catchment outlet a monitoring station was established (Fig. 1b). Three stations were located in the main stream network (stations 1, 2 and 4) and one station (station 3) was located in a tributary (Fig. 1b and c) which confluences with the main stream between stations 2 and 4. The point and reach-scale measurements were conducted between stations 2 and 4 and the catchment-scale investigations were carried out at each of the four stations.
Precipitation and discharge during 2011–2012. Upper part: hourly precipitation measured 6 km northwest of station 4. Lower part: the measured discharge at the catchment outlet (station 4).
Catchment characteristics and land use for each sub-catchment, with
mean annual discharge, catchment size, specific discharge, distance from the
source
Between stations 1 and 4, the stream flows from east to west with a mean
gradient of 1 ‰ receiving four main tributaries (Fig. 1b).
Between stations 2 and 4 there is a small inlet from a fishery,
constantly carrying a discharge of 0.07 m
For the point-to-reach-scale investigations, a 1 week campaign was carried out between stations 2 and 4 during the low-flow period of 9–15 June 2012 where point-scale VTPs and reach-scale DTS and ADCP measurements were conducted (Fig. 1c). The low-flow period was chosen to minimise the risk of surface discharge to the stream. The catchment-scale studies were conducted during three different rain events in 2012, one in spring, one in summer and one in autumn (Fig. 2). The decision of monitoring the three selected rain events were based on weather forecasts of upcoming large rain events, combined with antecedent medium to low stream-flow conditions. During the rain events, samples of stream water were collected at stations 1–4. Stream discharge at the catchment outlet and precipitation values during the investigation period are shown in Fig. 2. In addition, hydraulic heads were measured several times in piezometers installed in riparian zones/wetlands at stations 1, 2 and 4 (relative position of the screens shown on Fig. 3) to define hydraulic conditions at the stations. The different types of measurements are summarised in Table 2.
Manually interpolated cross-sectional contour maps of hydraulic
heads.
Deeper groundwater temperature in Denmark equals the annual average air
temperature of
Based on the VTP measurements, vertical groundwater fluxes were estimated by
fitting the steady-state analytical solution of the one-dimensional
conduction–convection equation (Bredehoeft and Papadopulos, 1965) to the
measured temperature data as described by Schmidt et al. (2007) and Jensen
and Engesgaard (2011):
Summary of sampling periods and data collection methods, with the scale covered by the method, the method/instrument and time of measurement.
As previous studies (Jensen and Engesgaard, 2011; Karan et al., 2013) in the
same area only detected moderate seasonal changes in streambed temperatures,
the steady-state conditions were assumed to be valid for the study period in
June. For each VTP,
During the June 2012 campaign DTS was used for reach-scale investigations of the groundwater discharge dynamics. A BruSteel fiber optic cable connected to a Sensornet Oryx-SR system was deployed along the middle of the stream on the sediment–water interface in three layouts (A, B and C) to cover the whole length of the stream section (Fig. 1c). To avoid damage of the fiber optic cable, no measurements were made between 1366 and 1530 m in the downstream direction from station 2 (Fig. 1c) due to remnants of a weir.
For each layout streambed temperature data were collected with double-ended
measurements of 10 min integration times and a 1.01 m spatial averaging
interval. Each installation was calibrated by running approximately 30 m
fiber optic cable through a calibration bath. The precision of the
installations is shown in Table 3. In each layout streambed temperature time
series of 22–23 h were collected with different starting times (Table 3),
but results are presented by aligning the measurements relative to the time
of day. Under the temperature conditions of the June campaign, low streambed
temperatures could indicate concentrated discharge zones. However, due to
different daily air temperatures, the decrease in streambed temperatures at
the potential concentrated discharge sites was not directly comparable
between the layouts. Hence, in order to compare streambed temperatures
measured at different days at different locations, the strength of the
groundwater temperature signal for each measurement location was calculated as
In the 2450 m long stream section between stations 2 and 4 (Fig. 1c), differential gauging of stream discharge was carried out during the June 2012 campaign for detection of the reach-scale variability of groundwater discharge. Stream discharge was measured with an ADCP StreamPro manufactured by Teledyne RD Instruments. The ADCP StreamPro has a four-beam 2 MHz transducer, a sampling frequency of 1 Hz and estimates discharge based on measured water velocities and a cross-sectional area. The ADCP is mounted on a platform and tethered across the stream, perpendicular to the main flow direction.
Discharge measurements were conducted for each 200 m in layout A and C and for each 150 m in layout B (Fig. 1c). The distances of 150 and 200 m between ADCP discharge measurements were chosen, based on a pilot study conducted in early spring 2012. This pilot study showed that at least 150 m between measurements were needed to obtain discharge increases larger than 5 %. The ADCP measurement procedure was optimised according to recent recommendations (Mueller and Wagner, 2009; Muste et al., 2004a, b) and a minimum of 10 discharge measurements with an average deviation less than or equal to 5 % were made at each location in order to minimise the uncertainty of the discharge estimates. A permanent gauging station was installed at the catchment outlet (station 4) for continuous discharge estimation based on the stage–discharge relation, continuous water stage measurements (OTT Thalimedes pressure transducer) and monthly current metre control measurements of discharge (Rantz, 1982; Herschy, 1999).
Time of DTS stream bed temperature measurements with the length and precision of each layout.
Different sources of runoff in the sub-catchments were studied based on the
analysis of the stream water EC and stable isotope fractions
The Teledyne samplers were programmed to collect 700 mL for every sampling,
and immediately after each sampling round of 24 samples, subsamples of 20 mL
were taken and sealed in plastic bottles and stored at a temperature of
4
To estimate the sources of runoff in the four different sub-catchments
during different events, a one-tracer, two-component hydrograph separation was
conducted (Sklash and Farvolden, 1979). The stream water was separated into
pre-event and event water fractions on the basis of the measured
An additional chemical one-tracer, two-component hydrograph separation was
conducted based on stream water EC. Thereby the fractions of sub-surface and
surface water can be estimated, where sub-surface water refers to the water
which has passed through the mineral soil, and surface water refers to water
which has not infiltrated the mineral soil (Genereux and Hooper, 1998). In
the case of an entirely groundwater-dominated stream network, the
sub-surface component will be equal to the groundwater component and surface
fractions will correspond to the rain component (Rodhe, 1998). Hence, any
discrepancies between the pre-event and sub-surface fractions can indicate
the likely presence of additional components (Wels et al., 1991). The same
mixing equation (Eq. 4) as used for the
Uncertainties in the pre-event water fractions inherent from uncertainties
in determination of the signatures used in Eq. (4) were calculated based on
the procedure by Genereux (1998). This method is based on an uncertainty
propagation technique using Gaussian error estimators, and was calculated at
the 0.05 confidence level. Uncertainties in EC and
In the majority of the piezometers installed at stations 1, 2 and 4, the
groundwater table was less than 2 m below ground during all measurements
conducted in the period December 2011 to June 2013. Due to the limited
fluctuations observed in groundwater levels at stations 1, 2 and 4
(
DTS, VTP and ADCP measurements from the campaign sampling in June 2012.
The results from the point-to-reach-scale investigations conducted during
the June 2012 campaign are summarised in Fig. 4. DTS measurements revealed a
number of focused groundwater discharge sites with their location generally
confirmed by the ADCP differential flow gauging. In layout A the ADCP
measurements showed only a slight net increase in stream discharge along the
first 400 m coinciding with no distinct temperature anomalies detected by
the DTS. However, at 600, 640, 705, 735, 800 and 825 m in the downstream
direction, colder streambed temperatures were detected by the DTS (Fig. 4a)
potentially indicating high groundwater discharge supported by an increase
in stream discharge of approximately 14 % along the layout (Fig. 4). In
the last 200 m of layout A the inflow of the tributary, which had an average
discharge of 0.23 m
Layout B revealed the largest spatial variability in groundwater discharge
of the three layouts, with both losing and gaining sections (Fig. 4).
The losing section was detected by the ADCP at the beginning of layout B
causing the stream discharge to decrease with approximately 13 % (Fig. 4).
However, ADCP measurements in the main stream suggested that stream
water is already recharging at the very last section of layout A, since only
an increase of 0.15 m
The most gradual net increase in discharge was observed along layout C by ADCP measurements and confirmed by several cold streambed temperature zones indicated by the DTS, suggesting more diffuse groundwater inflow compared to layout A and B. In layout C the most pronounced cold temperature anomalies were detected at the downstream end at 1900, 1980, 2285, 2380 and 2415 m (Fig. 4a). Due to a rain event on 9 June, the air temperature decreased and therefore the lowest streambed temperatures of all layouts were measured in layout C. The rain event also caused the stream discharge to be slightly higher on average during the first round of measurement compared to the second round (Fig. 4b). However, the event mainly occurred during the evening and night, and only the stream discharge pattern observed between the two most downstream ADCP measurements in layout C are suspected to be directly influenced by the rain event.
Summary of rainfall and runoff characteristics with rainfall intensity and duration, peak discharge, maximum discharge increase and number of rain samples for each precipitation event. September is divided into three sub-events.
Mean
To couple reach-scale groundwater discharge indications and point-scale
flux estimates VTP measurements were carried out at 18, 9, and 15 locations
in layout A, B and C, respectively, at the locations with the cold
temperature anomalies as shown by the DTS. The average root mean square error (RMSE) for the fitted
temperature profiles was 0.126
Correlation of the strength of the groundwater signal as recorded by DTS with upward groundwater fluxes estimated from VTPs. Measurements from each layout are separated by colour.
For layout A and B the comparison between estimated upward groundwater
fluxes and the strength of the groundwater signal (Eq. 3) at the
corresponding DTS locations showed a moderate correlation (significant on
the
Variability in tracer values measured during event sampling of
stream water.
Runoff and precipitation characteristics and pre-event fractions for
the September event.
The rainfall–runoff conditions during the three monitored events are
summarised in Table 4, and the signatures used for calculating the
hydrograph separations are shown in Table 5. No stream water samples were
collected prior to the September event. However, the May and September
events had similar antecedent conditions, and therefore the May pre-event
signatures were used as September pre-event and subsurface signatures of EC
and
The variability in stream water
The most significant event responses from all four sub-catchments were detected during the first part of the September event (Fig. 7). Station 1 showed the quickest and most pronounced response with the pre-event fraction reaching a minimum of 35 % (Fig. 7c) and a recovery time of approximately 9 h (recovery time is defined as the time it takes to reach pre-event concentrations). Stations 2 and 4 showed delayed and less pronounced event responses compared to station 1 with a minimum of 40 and 55 % pre-event water, respectively (Fig. 7d and f). Station 3 only showed a clear event response on 21 September, with the pre-event fraction being 70 % at the peak of the response (Fig. 7e). This response at station 3 was significantly delayed, approximately 15 h, compared to station 1, and showed a more gradually increasing response curve. Stations 2, 3 and 4 exhibited similar recovery times, approximately 24 h (Fig. 7c–f).
Generally, pre-event fractions were similar at all stations during peaks of the different events (Fig. 8a). Station 1 consistently showed the largest event responses and stations 2 and 4 reacted similarly but less pronounced than station 1. There was a tendency for station 4 to be damped in the pre-event responses as compared to station 2 (Fig. 8a). This is expected to be partly due to the inflow from the groundwater-dominated tributary between stations 2 and 4. Station 3 only showed modest peak response with min 70 % pre-event fractions during all events (Fig. 8a). The subsurface fractions showed similar responses at all stations as the pre-event fractions (Fig. 8b). However, with the exception of the September 1 event, the sub-surface fractions for stations 1, 2 and 4 varied significantly less than the pre-event fractions between events (Fig. 8). For instance, at station 4 the sub-surface fractions varied only between 80 and 90 %, whereas the pre-event fractions varied between 65 and 95 %. Calculated uncertainties at the peaks of the event and subsurface fractions depicted in Fig. 8 were all below 10 % (not shown).
The point-to-reach-scale investigation confirmed that the studied part of
the stream is groundwater dominated. On the reach scale, between station 2
and station 4, groundwater discharge to the stream resulted in approximately
30 % increase in total stream discharge. However, DTS and VTP measurements
showed that the spatial distribution of groundwater discharge in this
section is not homogeneous (Fig. 4), similarly to the DTS observations of
Lowry et al. (2007), Briggs et al. (2011) and the VTP-based flux estimations
of Schmidt et al. (2007) and Anibas et al. (2011). The large spatial
variability in groundwater discharge is most likely due to heterogeneity in
streambed hydraulic conductivity (Kalbus et al., 2006; Sebok et al., 2014),
which was also suggested by the streambed composition with interchanging
sand, gravel and clusters of macrophyte growth. The spatial heterogeneity
was also reflected at the point scale. Especially in layout C data showed
that even if the DTS streambed temperatures were higher than the mean, that
is
no high discharge was expected, upward fluxes up to 0.15 m day
DTS measurements have previously been used to locate and calculate
groundwater discharge to streams (Selker et al., 2006b; Briggs et al., 2011)
based on a temperature mixing approach combined with differential gauging
upstream and downstream of discharge sites. The DTS results from June 2012
also showed drops in streambed temperatures of 0.5–1
Consequently, our results suggest that a significant part of the groundwater
discharge along the studied 2.5 km long reach is concentrated in relatively
few focused zones. Hence, most likely the groundwater reaches the stream via
preferential flow paths governed by differences in streambed hydraulic
conductivity and hydraulic head conditions (Kalbus et al., 2006). Since
these focused high discharge zones will also carry the largest amounts of,
for instance, nutrients or potential contaminants (with a flux of up to
1.3 m day
So far, only few studies have endeavoured to confirm groundwater discharge sites indicated by the DTS with estimates of discharge based on either seepage metre data or vertical temperature profiles (Lowry et al., 2007; Sebok et al., 2013). This study shows that VTPs generally reflect the same spatial variability in groundwater discharge as the DTS (Fig. 4). There is a discrepancy of estimated fluxes and groundwater signal strength in the case of layout C, which is most likely caused by the reduced difference between the streambed temperatures and groundwater temperature (Fig. 4a). The DTS and VTP measurements of this study mostly complemented each other, confirming that cold streambed temperature anomalies correspond to locations of high upward groundwater fluxes. Thus, the combination of VTP and DTS measurements provides a useful tool for obtaining more robust groundwater discharge estimates in lowland groundwater-dominated streams, where the low temperature contrast between groundwater and surface water prevents discharge calculations by the method of mixing analysis.
The focused discharge locations detected by DTS and confirmed by VTPs agreed well with the net increases in stream discharge as measured by ADCP with the exception that DTS cannot identify losing stream sections. Contrary to the differential flow gauging of Briggs et al. (2011) where an acoustic Doppler velocimeter was used, the ADCP measurements here gave a good estimation of net groundwater discharge between measurement sections of 150–200 m spacing. In this study the combined ADCP and DTS methods made the detailed mapping of gaining and losing stream stretches possible, showing not only the net changes in discharge, but based on DTS also the approximate location of the focused discharge sites. However, a great logistical effort is required in order to map stream stretches longer than a few kilometres.
The discrepancy between the spatial resolutions of the methods is
illustrated when comparing the ADCP measurements to the DTS and VTP data.
Since the ADCP is expected to measure discharge within an uncertainty of
5 %, there exists a lower limit for measurement spacing during
differential gauging, since the change in discharge has to differ by more
than the 5 %. For this study, intervals of approximately 150–200 m were
close to the lower limit, especially for layout A and C, where the most
gradual increase in discharge was observed. Consequently, the ADCP method
was not capable of showing the same spatial variability in groundwater
discharge as the metre-scale DTS and the point-scale VTP measurements. For
this reason, it was also possible to still detect cold temperature anomalies
indicating groundwater discharge and relatively high upward fluxes of
0.43 m day
From the results of the hydrograph separations at the four different stations, it is clear that the most pronounced differences in runoff sources occur between station 1 and station 3 (Fig. 8a), with station 3 indicating a significantly larger and constant groundwater influence during events (maximum event water fraction was 30 %). The differences in forest cover (Table 1) could explain some of the differences in runoff sources during events since forest cover has previously been shown to significantly decrease surface runoff and enhance evaporation (Bosch and Hewlett, 1982; Brown et al., 2005). Also, the fact that station 3 was immediately surrounded by wetlands, while station 1 had a larger proportion of agricultural fields with tile drains in the near-stream area, may explain part of the larger proportion of event water observed at station 1, due to the importance of the riparian zones in terms of runoff processes (Tetzlaff et al., 2014; Vidon and Cuadra, 2010). The consistently high fractions of pre-event water observed at station 3 (Fig. 8a) suggest that the surrounding area has a shallow groundwater table as well as high hydraulic conductivity, allowing precipitation to seep to groundwater and preventing the presence of a zone of stored soil water which could otherwise have created a pre-event soil water component. This is supported by the fact, that no significant differences were seen between pre-event fractions and sub-surface fractions at station 3 (Fig. 8). Thus, most likely the assumption of two end members in the hydrograph separation was met, with pre-event water and sub-surface water representing the same groundwater component.
The large contribution from event water (maximum was 65 % event water) at station 1 could also be explained by the observed less strong groundwater gradients towards the stream, compared to the other three stations. Weaker groundwater gradients could potentially allow for a temporarily weakening of the groundwater discharge to the stream during large rain events, entailing a temporary dominance of surface and event water. Similar mechanisms were observed by Karan et al. (2014) where large rain events temporarily decreased groundwater discharge to Holtum stream. Also Gerecht et al. (2011) observed highly dynamic responses to rapid stage changes in terms of shifting between gaining and losing conditions in a groundwater influenced river. These observed differences in responses to large rain events between the studied catchments are of particular interest, in relation to being able to predict sensitive areas with the possibility of fast routing of nutrients and pollutants to streams. Catchments reacting similarly to station 1 would be more prone to fast routing of excess nutrients or pollutants than for instance catchments similar to station 3.
The discrepancies of around 10 % difference observed between subsurface
and pre-event factions at stations 1, 2 and 4 (Fig. 8) are similar to the
findings of Gonzales et al. (2009) for a lowland stream, and could indicate
the occurrence of a component which is not accounted for by either of the
two hydrograph separation methods (Wels et al., 1991; Hooper and Shoemaker,
1986). According to Karan et al. (2013), a shallow relatively young
groundwater component was discharging to the stream at station 4, supporting
the idea
that the stream-flow components could be divided into a deep groundwater
component discharging right beneath the stream channel, a shallow component
and a surface/event water component. However, there was no distinct
difference between the average EC and
Both the pre-event fractions as well as the subsurface fractions suggested that an event such as the one in April, with 15 mm rain and a resulting discharge increase of 30–50 %, constitutes a threshold below which runoff sources are not altered. These changes in contributing runoff sources between the sub-catchments are contrary to the findings of Gonzales et al. (2009). They found that their studied lowland stream system was at all times groundwater dominated, with minimum 90 % groundwater during events concluding that such consistently high influence of groundwater will most likely be found in the majority of similar lowland stream networks. However, our study illustrates that significant differences in event responses can exist among similar adjacent lowland catchments both in terms of the magnitude of event response and the response time (Fig. 7). These differences in catchment runoff sources during large rain events are important to take into account in water management practices, since a significant transport of phosphorus and nitrogen has been shown to be associated with storm events (Jordan et al., 2005; Kronvang and Bruhn, 1996; Stutter et al., 2008). Furthermore, the travel time (Flewelling et al., 2012) and origin (Clément et al., 2003) of discharging groundwater are decisive for the possibility of nitrate reduction.
Groundwater–surface water dynamics were studied in a groundwater gaining lowland stream in Denmark. The aim of this study was to combine hydraulic and tracer methods from point-to-catchment scale to assess the temporal and spatial variability of groundwater discharge and to assess the capability, limitations and synthesis of novel monitoring methods applied across the different spatial scales in terms of water management practices.
Significant groundwater discharge was observed, resulting in a total stream
discharge increase of approximately 30 % over a stream reach of 2400 m.
The groundwater discharge was found to be primarily confined in few distinct
zones, suggesting the presence of preferential flow paths. The major zones
of groundwater discharge were mapped by DTS and ADCP measurements and were
supported by point-scale VTP measurements indicating groundwater fluxes of
up to 1 m day
Based on this study it has been concluded, that despite a significantly groundwater influenced lowland catchment, there is still a high variability in the groundwater–surface water interaction. Hence, in relation to the growing demand of accurately estimating the transport of nutrients and other pollutants to streams, lakes and sea (e.g. European Commission, 2000), our study points to the challenges with variability in runoff sources in lowland streams. Our study emphasises the importance of considering the variability in groundwater discharge to streams across a range of scales. A strong focus should be not only on combining methods on the smaller scales, which has been subject to substantial investigations, but also on seeking to link to the catchment scale, where relatively simple hydrograph separations seem to be a useful tool even in lowland groundwater-dominated streams.
We would like to thank members of the Northern Rivers Institute, Aberdeen
University, for helpful discussions of data. We also thank Lars Rasmussen,
Jolanta Kazmierczak and Charlotte Ditlevsen for help in the field. This
study is part of the Hydrology Observatory, HOBE (