Articles | Volume 26, issue 23
https://doi.org/10.5194/hess-26-6227-2022
https://doi.org/10.5194/hess-26-6227-2022
Research article
 | 
12 Dec 2022
Research article |  | 12 Dec 2022

Disentangling scatter in long-term concentration–discharge relationships: the role of event types

Felipe A. Saavedra, Andreas Musolff, Jana von Freyberg, Ralf Merz, Stefano Basso, and Larisa Tarasova
Abstract

Relationships between nitrate concentrations and discharge rates (CQ) at the catchment outlet can provide insights into sources, mobilization and biogeochemical transformations of nitrate within the catchment. Nitrate CQ relationships often exhibit considerable scatter that might be related to variable hydrologic conditions during runoff events at sampling time, corresponding to variable sources and flow paths despite similar discharge (Q) rates. Although previous studies investigated the origins of this scatter in individual or in a few catchments, the role of different runoff event types across a large set of catchments is not yet fully understood.

This study combines a hydrological runoff event classification framework with low-frequency nitrate samples in 184 catchments to explore the role of different runoff events in shaping long-term CQ relationships and their variability across contrasting catchments. In most of the catchments, snow-impacted events produce positive deviations of concentrations, indicating an increased nitrate mobilization compared to the long-term pattern. In contrast, negative deviations occur mostly for rainfall-induced events with dry antecedent conditions, indicating the occurrence of lower nitrate concentrations (C) in river flows than their long-term pattern values during this type of event. Pronounced differences in event runoff coefficients among different event types indicate their contrasting levels of hydrologic connectivity that in turn might play a key role in controlling nitrate transport due to the activation of faster flow paths between sources and streams. Using long-term, low-frequency nitrate data, we demonstrate that runoff event types shape observed scatter in long-term CQ relationships according to their level of hydrologic connectivity. In addition, we hypothesize that the level of biogeochemical attenuation of catchments can partially explain the spatial variability of the scatter during different event types.

1 Introduction

Diffuse nutrient inputs in catchments are a challenge for water quality management (Paerl 1997; Stumpf et al., 2016). An excess of nutrients, such as nitrate, harms ecosystems by creating favorable conditions for eutrophication in water bodies and leading to biodiversity loss (GEA, 2017; EEA, 2019; Weitere et al., 2021). Fertilizer application on agricultural land remains the main source of nitrate contamination in human-impacted catchments, despite regulations of the past decades that stimulated a reduction of fertilizer application in Europe (Grinsven et al., 2012). Moreover, due to long-lasting legacy effects, a delay in the reduction of riverine nitrate concentration (C) was reported in many catchments (Tesoriero et al., 2013; Meter and Basu, 2017; Bieroza et al., 2018; Chang et al., 2021).

Long-term concentration–discharge (CQ) relationships are a valuable tool for analyzing water quality gradients and trends, and for developing water management strategies (Bowes et al., 2014). The shape of CQ relationships encodes export patterns and reflects the temporally varying quantities of critical substances such as nutrients delivered to streams (Godsey et al., 2009; Meybeck and Moatar, 2012; Rose et al., 2018). Depending on the slope of the log–log linear dependency of concentrations from discharge (Q), three different export patterns (Godsey et al., 2009) can be defined: dilution (negative slope), enrichment (positive slope) and neutral (no relationship between C and Q or slope close to 0). Differences in long-term CQ-relationships among catchments can be associated with differences in availability and spatial distribution of solute sources (Musolff et al., 2017; Dupas et al., 2019; Zhi et al., 2019; Casquin et al., 2021), their hydrologic connectivity (Seibert et al., 2009; Dupas et al., 2016; Covino, 2017) and biogeochemical processes within the soil and stream that can retain or permanently remove nitrate from stream water (Mulholland et al., 2008; Dupas et al., 2016; Moatar et al., 2017; Benettin et al., 2020).

Biogeochemical processes that affect nutrient cycles in soil and water might add variability to long-term CQ relationships. The effectiveness of the denitrification process, which removes nitrate from the soil, depends on periodic environmental factors such as temperature and soil moisture and the availability of electron donors (Korom et al., 2012; Ortmeyer et al., 2021). Instream removal processes are also more efficient during low flows and higher temperatures, adding more variability to the low-flow portion of the long-term CQ relationships (Dehaspe et al., 2021; Moatar et al., 2017). Moreover, the availability of nitrate sources is balanced by fertilizer application and mineralization of organic nitrogen compounds and hence varies in time, adding temporal variability to CQ relationships. Timing of fertilizer application is often unknown, and the mineralization processes depend on chemical soil conditions and environmental factors (e.g., soil moisture and temperature) that mediate communities of microorganisms (Curtin et al., 2012; Guntiñas et al., 2012). Average residence times of nitrate in agricultural catchments can last for decades, producing a legacy in soil (Meter et al., 2016; Puckett et al., 2011; Tesoriero et al., 2013; Vervloet et al., 2018) that can buffer the periodic effect of biogeochemical processes which reduces the variability in the concentration of nitrate (Basu et al., 2011; Bieroza et al., 2018; Thompson et al., 2011).

The scatter of CQ relationships might also be related to hydrologic conditions at the time of sampling (Knapp et al., 2020, Musolff et al., 2021), which are investigated for a large number of catchments only by a few recent studies (Minaudo et al., 2019; Pohle et al., 2021). Minaudo et al. (2019) showed that in most of the 219 French catchments, nitrate samples taken during baseflow conditions exhibit an enrichment export pattern, while during runoff events, a neutral or opposite pattern (dilution) prevails, generating scatter in the combined long-term CQ relationships. The cause of this scatter can be also traced to a variety of responses observed at the event scale in several studies with high-frequency data in single or a few catchments (e.g., Bowes et al., 2015; Lloyd et al., 2016; Koenig et al., 2017; Gorski and Zimmer, 2021).

Our study relies on low-frequency nitrate data, which are often used to build long-term CQ relationships (e.g., Cartwright, 2020; Diamond and Cohen, 2018). However, studies with high-frequency data found large variability in the CQ patterns during events (event CQ relationship; e.g., Knapp et al., 2020; Dupas et al., 2016; Vaughan et al., 2017) that might add scatter to the long-term CQ relationship. Disparate event CQ relationships in a catchment over time are mainly attributed to varying dominant flow sources (e.g., groundwater, shallow subsurface flow), antecedent wetness conditions (Inamdar et al., 2006; Knapp et al., 2020; Vaughan et al., 2017), time of fertilizer application (Bowes et al., 2015; Dupas et al., 2016; Outram et al., 2016), biogeochemical cycling (Heathwaite and Bieroza, 2021) and runoff event characteristics or types (Butturini et al., 2006; Bauwe et al., 2015; Chen et al., 2020; Knapp et al., 2020). For example, Winter et al. (2022) showed that in a few catchments located in Central Germany, runoff events generated by rainfall with dry antecedent conditions export lower nitrate concentrations due to lower hydrologic connectivity but exhibit a high variability of event CQ slopes. In contrast, Knapp et al. (2020) showed that using high-frequency concentration and discharge observations from one small forested catchment located in Switzerland during larger runoff events with dry antecedent conditions the slopes of the event CQ relationships are more positive due to the accumulation of nitrate in the soil during dry periods by atmospheric deposition and the subsequent mobilization by event water. Moreover, in several catchments in USA and Europe, snow-induced events were found to export high nitrate concentration (Koenig et al., 2017; Inamdar et al., 2006; Casson et al., 2014). Similarly, in the previously mentioned Central German catchments, Winter et al. (2022) found high nitrate concentrations and flat event CQ slopes during snow-impacted events, indicating that sufficient nitrate sources are available and most of the relevant flow paths are activated and connected to the stream during such events.

It was shown that hydrologic connectivity as a portion of the catchment connected to the stream via surface or subsurface pathways increases according to the wetness state of the catchment (Blume and van Meerveld, 2015; Jencso et al., 2009) and modulates export of nutrients at different scales. At seasonal scale, nutrient transport to streams can be increased with higher hydrologic connectivity in catchments with abundant sources (Martin et al., 2004; Veith et al., 2020; Guillemot et al., 2021). At event scale, the activation of different flow paths during different levels of hydrologic connectivity evaluated using shallow wells or models can partially explain changes in nitrate concentration during events (von Freyberg et al., 2014; Ocampo et al., 2006; Stieglitz et al., 2003). However, at the larger scale, such observations are not available.

At catchment scale, soil moisture or discharge rates are often used as proxy of hydrologic connectivity (e.g., Bracken et al., 2013; Jencso et al., 2009). The event runoff coefficient (i.e., a volumetric ratio of quick flow and input precipitation or snowmelt), which represents how efficiently streamflow responds to catchment water inputs, can also be considered as its proxy (e.g., Blume et al., 2007; Outram et al., 2016; von Freyberg et al., 2014). Higher runoff coefficients are associated with wetter antecedent catchment states, indicating that such conditions favor a more efficient rainfall–runoff response (Tarasova et al., 2018; Outram et al., 2016) and possibly activation of more surface and subsurface hydrologic flow pathways that facilitate fast transport of water and nutrients from the landscape to the stream (Blume and van Meerveld, 2015; Hardie et al., 2011; Stieglitz et al., 2003).

New approaches to characterize and classify runoff events according to hydrologic conditions offer a possibility to efficiently aggregate information about the antecedent wetness state of catchments and characteristics of inducing events (e.g., rainfall, snowmelt) and to distinguish events with contrasting hydrological responses for a large number of catchments (Tarasova et al., 2020). Such classification of event types combined with concentration of nitrate in stream water might unravel scatter in long-term CQ relationships as exemplified in Fig. 1. In Fig. 1a, biweekly nitrate data are associated with the event type at the time of stream water sample collection. When these data are plotted in the log–log CQ space (Fig. 1b), some event types exhibit positive (higher concentration) or negative (lower concentration) deviations from the long-term CQ relationship. Our study aims for the first time to investigate the presence of systematic deviations in long-term CQ relationships produced by different runoff event types for a large dataset of catchments.

https://hess.copernicus.org/articles/26/6227/2022/hess-26-6227-2022-f01

Figure 1(a) Time series of daily discharge and biweekly grab sample nitrate concentrations during event and no-event conditions in the Naab River at the gauge of Unterköblitz, Bavaria over a period of 5 years. Event types are differentiated by colors (see Fig. 3 for details). (b) Double logarithmic plot of CQ pairs for samples (from 2000 to 2012) taken during different event types and no-event conditions. Dashed black lines show the long-term CQ relationships (same line in each subplot) obtained from linear regression in a double logarithmic plot of CQ values for all available samples.

Download

We hypothesize that these deviations are related to the differences in nitrate transport during these event types and we investigate such deviations from the long-term CQ relationships in 184 German catchments. Specifically, our goal is to examine the effect of runoff event types on the observed scatter in CQ relationships by addressing the following research questions:

  1. Do samples collected during different event types deviate differently from the long-term CQ relationships observed at the catchment outlets?

  2. Which climatic and landscape characteristics explain differences in the observed CQ deviations among German catchments?

  3. Which are the potential mechanisms that explain the direction and magnitude of CQ deviations for different event types?

Understanding the nature of nitrate deviations from the long-term CQ relationships might provide useful information for water quality managers to reduce the risk of extreme nitrate loads to water bodies, as well as improve sampling campaigns to better capture nitrate CQ scatter.

2 Methods

2.1 Study catchments and data

In this study, we analyzed low-frequency (biweekly to monthly) nitrate concentration data from 184 mesoscale catchments in Germany for the period from 2000 to 2015. The data were obtained from the water quality and quantity database of Germany (Musolff, 2020; Ebeling et al., 2021) in combination with a recently developed classification framework of runoff events (Tarasova et al., 2020). Similar to Ebeling et al. (2021), we exclude the data prior to the 2000s to avoid impacts of improved wastewater treatment technologies in Germany. In total, we considered 33 713 nitrate samples.

Sizes of study catchments range from 95 to 23 615 km2 (with a median size of 704 km2) and cover all four main German natural regions: the North German Plain, Central Uplands, South German Scarplands and Alpine Foreland (Fig. 2a). The climate varies from temperate oceanic to temperate continental from west to east. Mean annual precipitation ranges from 567 mm in the lowland northeastern catchments up to 1379 mm in the alpine catchments in the south. The predominant land use in the study catchments is agriculture, with a median coverage among catchments of 50 % and a range from 13 % to 84 %. The median portion of catchment area covered by forest is 41 % of the catchment area (Fig. 2b).

https://hess.copernicus.org/articles/26/6227/2022/hess-26-6227-2022-f02

Figure 2(a) Study area and stations of nitrate concentration measurements in stream water. Gray lines show catchment boundaries. Catchment outlets (points) are color-coded according to the long-term export pattern (dilution, neutral and enrichment). Blue lines show the main rivers. The background color map corresponds to the elevation. Purple labels indicate German natural regions. (b) Area, fraction of agriculture, fraction of forest and mean annual precipitation of study catchments grouped according to export patterns (dilution, neutral and enrichment). Red lines show medians of boxplots and significance of median differences between adjacent boxplots was estimated using the Kruskal–Wallis test (displayed as  for p< 0.05 and  for p< 0.01).

The runoff event classification framework of Tarasova et al. (2020) considers runoff events identified from daily discharge data in catchments with no major flow regulations. The location of the discharge stations does not always coincide with water quality stations in the dataset of Ebeling et al. (2021). Both datasets are linked by pairing stations that are located on the same stream and differ less than 20 % in their drainage areas. These were considered as identical outlets, similar to Guillemot et al. (2021). The mean overlap between drainage areas of the corresponding outlets from the two datasets is 95 % with a standard deviation of 5 %.

2.2 Identification and classification of hydrological events

Runoff events and corresponding precipitation events were separated using an automated time-series approach developed by Tarasova et al. (2018). The method was applied to daily discharge and precipitation data obtained from the REGNIE dataset (Rauthe et al., 2013). The method includes baseflow separation, precipitation attribution (i.e., corresponding inducing events (rainfall and/or snowmelt) are linked to runoff events) and an iterative procedure to adjust site-specific thresholds for the refinement of multi-peak events. The median event duration is 12 d with a standard deviation of 7.7 d. The shortest event duration is 1 d, however, 95 % of the identified events exhibit a duration of 3 or more days. Each identified runoff event was then classified in the first place by considering the nature of inducing events (rainfall, mixture of rainfall and snowmelt or rain-on-snow) (Fig. 3a) using the proportions of rainfall and snowmelt in the total volume of precipitation events (Table S1 in the Supplement). In the second step, we considered the antecedent wetness state (wet or dry) by accounting for the catchment-averaged soil moisture state prior to the event. Catchment-average snow water equivalent and soil moisture were simulated by the mesoscale hydrological model (Samaniego et al., 2010; Kumar et al., 2013) and provided in Zink et al. (2017). Additionally, the classification considers spatial organization of soil moisture within the catchment using the spatial coefficient of variation of soil moisture, classifying events as uniform or patchy, with the latter corresponding to highly variable soil moisture within the catchment. A more detailed description of the classification framework is provided in Tarasova et al. (2020).

https://hess.copernicus.org/articles/26/6227/2022/hess-26-6227-2022-f03

Figure 3(a) Hierarchical scheme for event classification (modified from Tarasova et al., 2020; classification criteria are provided in Table S1). Colored dots located next to the five different event types indicate their markers. (b) Cobs and Qobs are observed concentration and discharge, Cfit is the nitrate concentration estimated from fitting the long-term CQ relationship with a linear relation in log–log space, and res is the residual value. (c) CQ plots for three different catchments attributed to different long-term nitrate export patterns based on the logClogQ slope b, i.e., dilution (b<0.1, the Würm River in Pforzheim), neutral (b 0, the Wupper River in Opladen) and enrichment (b> 0.1, the Naab River in Unterköblitz).

Each nitrate sample was linked to either no event (No.event), or to one of the five event types (Fig. 3a): rain-on-snow (Rain.on.snow), mixture of rainfall and snowmelt (Mix), rainfall during wet antecedent conditions (Rain.wet), rainfall during dry antecedent conditions with spatial uniform distribution of soil moisture (Rain.dry.uniform) and rainfall with dry antecedent conditions with heterogeneous spatial distribution of soil moisture (Rain.dry.patchy). Note that we simplified the event types to increase the number of nitrate samples of each event type.

2.3 Long-term CQ export patterns

For each catchment, the long-term CQ relationship was derived as a linear regression between nitrate concentration (C) and discharge (Q) in the log–log space (Fig. 3b). Based on the slope of the long-term CQ relationships (b), we grouped all study catchments according to three different long-term CQ export patterns: dilution (b< 0.1) refers to a limitation of sources during high flows, enrichment (b> 0.1) is related to a transport limitation with abundant sources or solute uptake during low flows (Moatar et al., 2017) and neutral (b 0) indicates no monotonic relationship between C and Q. As stated by Ebeling et al. (2021), this latter group exhibits largely invariable concentration with low ratios of coefficients of variation (CVc/CVq). Three different catchments are shown as an example of each export pattern in Fig. 3c.

2.4 Quantifying the deviations from long-term CQ relationship

For each catchment, we want to quantify whether samples taken at a specific event type show systematic deviations from the long-term CQ regression compared to all samples. We quantified the deviation of each grab sample from the long-term CQ relationship for each catchment by computing the corresponding residual concentration from the long-term CQ linear regression line (Fig. 3b). Resulting residuals were subsequently grouped according to the hydrological event type at the time of sampling.

Due to the variable number of grab samples attributed to different event types (Fig. S1 in the Supplement), for each catchment, we performed a bootstrapping procedure that can explicitly handle unbalanced data by iteratively comparing two random subgroups of samples with the same size (undersampling method, e.g., Branco et al., 2015). The procedure is implemented in the following way for each catchment: n nitrate samples of a certain event type and the same number of nitrate samples from all samples (general pattern) are chosen randomly with replacement (i.e., each data point can be chosen more than once, following bootstrapping procedure). The difference of median residuals of an event type and residuals of the general pattern is then the measure of deviation of a corresponding event type from the long-term CQ relationship (Δres). We obtained this measure 10 000 times to robustly compute its distribution (Fig. S2) and median value (Δres50). The number of samples n was chosen for each catchment and event type according to the number of nitrate samples available for the corresponding event type. For each catchment, event types with less than 10 nitrate samples are excluded from the analysis. The median number of nitrate samples among all study catchments and event types is 27.

For each catchment and from all the iterations, we obtain the median deviations between event types and the general pattern (Δres50). In order to evaluate the persistence of CQ deviations across catchments, we tested the significance of Δres50 across catchments for each event type using the non-parametric Kruskal–Wallis test (Kruskal and Wallis, 1952) at the significance level α= 0.05.

Low-frequency datasets such as the one used in our study might contain samples collected during different phases of the event hydrograph (e.g., falling or rising limb). This might hamper the interpretability of the results due to possible bias in observed nitrate concentration linked to the time of sampling and the hysteresis effect revealed in high-frequency observations (e.g., Lloyd et al., 2016; Vaughan et al., 2017). In fact, Pohle et al. (2021) showed systematic differences in nitrate concentration between samples collected during rising and falling limbs for numerous catchments in Scotland. To understand the potential effect of the hysteresis on the deviations from long-term CQ (Δres50), we repeat the bootstrapping procedure described above considering samples collected during the rising limb, falling limb and near the event peak (near-to-peak). The rising limb of a runoff event starts at the beginning of the event and finishes 1 d before the day of the peak discharge. The falling limb starts 1 d after the day of the peak discharge and finishes at the end of the runoff event. In addition, we defined near-to-peak as samples collected from 1 d before to 1 d after the day of the peak discharge. Of the total samples taken during runoff event types, 34 % correspond to the rising limb, 55 % to the falling limb and 30 % to near-to-peak. Notice that the definition of near-to-peak samples allows some overlap with the other two groups of samples to use a more balanced number of samples than considering samples collected on the day of the peak of discharge only (11 % of the samples were collected during the day of the peak discharge).

2.5 Catchment descriptors and relationships to CQ deviations

In order to explore the differences of deviations from the long-term CQ relationships across the catchments, we examined the Spearman rank correlation of median residuals for each catchment with various catchment descriptors. Here, we only examine catchment descriptors that were previously identified as primary controls of the nitrate CQ export patterns in Germany (Ebeling et al., 2021). This includes topographic descriptors (median topographic wetness index, median slope and area); land cover descriptors (fraction of agriculture, forest and artificial surface); soil and aquifer descriptors (median soil depth and fraction of sedimentary aquifer); nitrate sources descriptors (nitrate surplus, agricultural horizontal heterogeneity, nitrate vertical ratio); and hydrometeorological descriptors (aridity index, mean annual potential evapotranspiration, precipitation and temperature) (Table S2). Detailed derivations of the above-mentioned catchment descriptors are provided in Ebeling et al. (2021).

3 Results

3.1 Frequency of runoff event types

Stream water samples taken during runoff event conditions account for 58 % of all samples. These samples are classified to one of the five event types as follows: 18 % – Rain.dry.patchy, 11 % – Rain.dry.uniform, 15 % – Rain.wet, 7 % – Rain.on.snow and 7 % – Mix.

On average across catchments, the fraction of samples taken during each event type vary at different discharge rates. Above median discharge rate, 74 % of all samples correspond to an event and the event types Rain.wet, Rain.on.snow and Mix occur more frequently (Fig. 4a). In contrast, only 49 % of samples below median discharge rate were taken during an event and most of these grab samples correspond to Rain.dry.patchy and Rain.dry.uniform types.

https://hess.copernicus.org/articles/26/6227/2022/hess-26-6227-2022-f04

Figure 4(a) Mean fraction of samples linked to each event type according to each catchment decile of discharge, and (b) seasonal distribution of mean fraction of samples linked to each event type in the study catchments.

Download

The frequency of event types also varies seasonally (Fig. 4b). In winter, most of the grab samples were taken during Rain.on.snow, Mix and Rain.wet event types. In the spring months, Rain.dry events become more frequent than Rain.on.snow, Mix and Rain.wet event types. During summer, most of the samples were taken either under No.event conditions or during Rain.dry.uniform and Rain.dry.patchy events. In autumn, the frequency of grab samples taken during Rain.wet, Rain.on.snow and Mix event types increases.

3.2 Long-term CQ relationships and deviations during event types

We computed long-term nitrate CQ relationships for the 184 catchments, obtaining slopes (b) from 0.6 to 1.48, with a mean of 0.13. In total, 88 study catchments exhibit neutral patterns, 80 catchments are characterized by enrichment patterns and only 16 catchments show dilution patterns. Across all catchments, the median R2 value of the long-term CQ relationship was low (0.14), indicating the presence of considerable scatter in the regressions.

We explored the residuals (res) of all nitrate data from all catchments and found that 65 % and 68 % of the samples taken during Rain.on.snow and Mix event types, respectively have positive residual values, indicating that concentrations were higher than the long-term log–log linear CQ regressions. In contrast, 69 % and 60 % of the samples during Rain.dry.patchy and Rain.dry.uniform events, respectively have negative residuals values. We found a less clear picture for samples taken during Rain.wet events and No.event conditions with 53 % and 56 % of positive residuals, respectively (Fig. S3).

https://hess.copernicus.org/articles/26/6227/2022/hess-26-6227-2022-f05

Figure 5Median deviations of nitrate concentrations from the long-term CQ relationships (Δres50). (a) Δres50 values of different event types for each catchment. On the right-hand side of each map, boxplots show the distribution of Δres50 values across catchments for each event type (box limits represent the interquartile range and whiskers correspond to the 5th and 95th percentiles). (b) Heatmap of Δres50 values averaged across different groups of catchments, considering all nitrate data for each event type and No.event. The first three columns of the heatmap correspond to one of the long-term export patterns (i.e., dilution, slope b< 0, neutral, slope b 0, and enrichment, slope b> 0) and the fourth column corresponds to all study catchments. Bold font and  indicates significant differences (Kruskal–Wallis test, p< 0.05) between median deviations across catchments for each event type and median deviation across catchments of all nitrate samples.

We found strong differences in median deviations from the long-term CQ relationships (Δres50) among different event types (Fig. 5a). Rain.on.snow and Mix event types often have more positive Δres50 values (79 % and 93 % of the study catchments correspondingly) (Fig. 5b) when compared across catchments. Contrastingly, Rain.dry.patchy and Rain.dry.uniform event types show negative values of Δres50 more often (96 % and 61 % of the study catchments), with Rain.dry.patchy events showing stronger deviations. Contrasting behavior between snow-impacted events (i.e., Mix and Rain.on.snow) and rainfall events with dry antecedent wetness conditions (Rain.dry.patchy and Rain.dry.uniform) occurs across most of the study catchments independently of their long-term export pattern (Fig. 5b). For Rain.wet events, deviations can be negative as well as positive (52 % and 48 % of study catchments, respectively) with a median of Δres50 across catchments close to zero (Fig. 5a). For samples that were taken during No.event conditions, the Δres50 value is slightly positive in 85 % of all catchments.

The sign of CQ deviations are in line with observed nitrate concentration during different event types (Fig. S4). Negative residuals during Rain.dry.patchy and Rain.dry.uniform events coincide with lower nitrate concentrations for most of the catchments, independent of the long-term CQ pattern. Similarly, during Rain.on.snow and Mix events, positive CQ deviations correspond to nitrate concentrations higher than median for most of the catchments with a neutral or enrichment CQ pattern. For catchment with the dilution export pattern, nitrate concentration for Rain.on.snow and Mix events is similar to the average, however higher discharge generates positive residuals in this case.

We analyzed the influence of the sampling time within runoff events separating samples taken during the rising limb, near-to-peak and falling limb. Although there are certain data limitations for a few groups of samples (gray tiles in Fig. S5b), we are able to reproduce the analyses for most of the cases. Similar to the case when using all samples (Fig. 5b), the values of Δres50 for samples taken during the rising limb, near-to-peak and falling limb are mostly positive for Rain.on.snow and Mix events and negative for Rain.dry.patchy and Rain.dry.uniform. Our results confirm that the time of sampling during runoff events does not affect our findings regarding median CQ deviations for different types of runoff events.

Although the sign of CQ deviation is consistent across catchments for most of the event types, the magnitude of deviation varies across catchments (Fig. 5a). The variability of Δres50 expressed as interquartile ranges across catchments (boxplots in Fig. 5a) is the lowest for the samples taken during No.event conditions (0.03) and Rain.wet events (0.06). The largest variability was detected for Rain.dry.patchy events (0.1), followed by Mix (0.09) and Rain.dry.uniform events (0.09).

https://hess.copernicus.org/articles/26/6227/2022/hess-26-6227-2022-f06

Figure 6Spearman rank correlation coefficient between deviations of nitrate concentrations from the long-term CQ relationships (Δres50) of a particular event type across study catchments and catchment descriptors. Significant correlations are indicated by bold font and  for p< 0.05 and  for p< 0.01.

Download

3.3 Variability of CQ deviations across German catchments

We analyzed the spatial variability of CQ deviations for different event types (Fig. 5a) computing Spearman rank correlations between deviations and catchment descriptors. We found significant correlations between Δres50 for each event type and catchment descriptors. Topographic properties (i.e., median slope and topographic wetness index) have the strongest correlation to the Δres50 values of almost all event types (Fig. 6). Specifically, flatter catchments (low median topographic slope) with greater soil depths that are mostly located in the North German Plain and Alpine Foreland tend to exhibit more positive residuals for Rain.wet, Rain.on.snow and Mix events, and more negative residuals for Rain.dry.patchy events and samples are taken during No.event conditions (Fig. 5a). Catchments with these characteristics often show high agricultural land cover (Fig. S6), however the fraction of agriculture shows less significant correlations with Δres50 than topographic descriptors. Moreover, in catchments with larger fractions of water-impacted soils (e.g., stagnosols, semi-terrestrial, semi-subhydric, subhydric and moor soils), we found more positive residuals for snow-impacted events (Rain.on.snow, Mix) and more negative residuals for Rain.dry.patchy events. These catchments are often located in Central East or North-West Germany.

Correlations between Δres50 and fraction of agriculture in the catchments are less significant than those with topographic descriptors (Fig. 6). Instead, we observed strong correlations between Δres50 and the fraction of forest (p< 0.01). Forested catchments show less positive Δres50 values for Rain.on.snow events and less negative values for Rain.dry.patchy events. However, we also noticed that the fraction of forest is positively correlated with topographic slope and negatively correlated with soil depth and the fraction of agriculture (Fig. S6).

Nutrient source descriptors were also significantly correlated with Δres50. Horizontal heterogeneity of agricultural sources correlates negatively with Rain.dry.pachy residuals and the vertical concentration ratio of nitrate correlates negatively with Δres50 values of Rain.dry.patchy and No.event conditions. Nitrate surplus is significantly related only to Rain.wet residuals.

3.4 Relationship between hydrologic connectivity and event type variations in residuals

We examined event runoff coefficients corresponding to different catchments and event types to link the relation between hydrologic connectivity for these event types and corresponding deviations of their samples from the long-term CQ relationships (Fig. 7a). Catchment median event runoff coefficients exhibit a coefficient of variation of 41 % across catchments. Nevertheless, variability of median runoff coefficients across event types for single catchments is larger in most of the cases, with coefficients of variation from 12 % to 118 % and a median value of 67 % across catchments. We found that event types with significantly higher median runoff coefficients also exhibit significant differences in Δres50 values (Fig. 7b and c). Only Mix and Rain.on.snow events have similar runoff coefficients and similar Δres50 values.

https://hess.copernicus.org/articles/26/6227/2022/hess-26-6227-2022-f07

Figure 7(a) Relationship between Δres50 for each catchment and event type and median runoff coefficient (rc); runoff coefficient is not defined for No.event. (b) Variability of runoff coefficients (rc) for each event type and (c) median residuals for each event type. Significance of median differences between adjacent boxplots was estimated using the Kruskal–Wallis test (displayed as  for p< 0.001).

Download

4 Discussion

4.1 Direction and magnitude of CQ deviations for different event types

We found systematic differences in the direction and magnitude of deviations of nitrate concentrations (Δres50) from the long-term CQ relationships during different types of runoff events despite the large variety of study catchments (Fig. 5). In the following paragraphs, we will discuss potential mechanisms that can explain the variability of CQ deviations across event types.

Positive deviations for nitrate concentrations during snow-impacted events (i.e., higher nitrate concentration compared to the general CQ pattern) are in line with previous studies that have shown an increase of nitrate concentration in stream water during snow-impacted events in forested and agricultural catchments (Inamdar et al., 2006; Casson et al., 2014; Koenig et al., 2017). This is in line with Winter et al. (2022), who showed using high-frequency data that snow-induced events export the highest nitrate concentration compared to other event types in six German catchments with mixed land use. Our results clearly show that snowmelt does not generate lower concentration of nitrate compared to the long-term CQ relationship, although this might be expected due to lower nitrate concentration in snowfall than in stream water from agricultural catchments (Johannsen et al., 2008). Instead, higher concentration indicates enhanced nitrate transport from soil sources with no source limitation during these types of events. We argue that during snow-impacted events, hydrologic connectivity is high between nitrate sources and streams due to elevated wetness conditions (Stieglitz et al., 2003), which is consistent with previously reported high nitrate concentration during the winter period (Martin et al., 2004; Ocampo et al., 2006; Yang et al., 2018). Due to excessive catchment wetness during snow-impacted events, a high amount of new water transported by faster and shallower pathways can reach the stream (a so-called inverse storage effect; Fang et al., 2019), mobilizing large amounts of nitrate available in the soil (Yang et al., 2018). In addition, during these events, the mobilized water is less affected by biogeochemical processes due to lower microbial activity induced by low temperature during snow-impacted events (Johannsen et al., 2008).

Furthermore, our analysis shows that Rain.dry.uniform and Rain.dry.patchy events generate lower nitrate concentrations compared to the other types of events or No.event conditions (Fig. S4), producing strong negative CQ deviations (Fig. 5). Along the same lines, Winter et al. (2022) showed that runoff events with dry antecedent conditions exhibit lower concentration compared to other event types in six German catchments with mixed land use. There are two possible explanations for the occurrence of this phenomenon. On the one hand, Rain.dry.uniform and Rain.dry.patchy events occur more often during the dry season, when nitrate concentrations are reported to be lower (House et al., 2001; Guillemot et al., 2021) due to a hydrological disconnection between agricultural sources and streams under dry conditions, as well as higher biogeochemical nitrate removal processes, including biotic uptake and denitrification (Mulholland et al., 2008; Rode et al., 2016; Lutz et al., 2020; Johannsen et al., 2008). On the other hand, during runoff events with dry antecedent conditions, nitrate concentrations can be diluted below pre-event concentration levels. This is shown by high-frequency observations in agricultural catchments that report more frequent negative event CQ slopes during the dry season (Winter et al., 2021, 2022; Zhang et al., 2020; Koenig et al., 2017). In such cases, nitrate concentration decreases compared to pre-event concentrations due to hydrologic disconnection between streams and agricultural land, and the growing importance of runoff generated from riparian zones (Fang et al., 2019; Grayson et al., 1997; McGlynn and Seibert, 2003), which are known to buffer nitrate inputs due to high denitrification potential (Ocampo et al., 2006; Cole et al., 2020; Sabater et al., 2003). Our results show that the combined effect of lower pre-event concentration and further decrease in concentrations due to runoff events magnifies the observed negative deviations of nitrate samples from the long-term CQ relationships. Nevertheless, the data available for this study do not allow us to quantify the contribution of individual effects of these two factors on the scatter of long-term CQ relationships of nitrate. On the contrary, studies in pristine headwaters and forested catchments found that rainfall events with dry antecedent conditions can mobilize large amounts of nitrate, increasing the concentration in streams (Knapp et al., 2020; Koenig et al., 2017). Since these findings are based on the observations in a single or only a few catchments with limited agricultural activity, different nitrate sources, such as atmospheric deposition or nitrate fixation and nitrate accumulation in soil between events, might be more relevant. Agriculture is a dominant land use type in the catchments used in this study (median fraction of agricultural land is 50 %), therefore a considerable nitrate accumulation in soil as the result of fertilization dominates over any other nitrate source (Häussermann et al., 2020; Lassaletta et al., 2014), explaining the discrepancy between our findings and the results from pristine headwaters and forested catchments on the role of rainfall events with dry antecedent conditions for nitrate mobilization.

Different from runoff events with dry antecedent conditions, we found that nitrate grab samples taken during No.event conditions exhibit slightly positive deviations, indicating higher concentrations compared to the long-term CQ relationships. No.event samples also exhibit higher nitrate concentrations (Fig. S4) compared to Rainfall events with dry antecedent conditions (i.e., Rain.dry.patchy and Rain.dry.uniform), with both groups of samples being collected during relatively low discharge conditions (Fig. 4a). This suggest that the lack of dilution during No.event conditions might produce more positive residuals.

4.2 The role of hydrologic connectivity between different event types

The hypothesized role of hydrologic connectivity on shaping nitrate deviations during runoff events is supported by the relation between event runoff coefficients and the deviation of nitrate concentrations from the long-term CQ relationships for different event types (Fig. 7). Higher runoff coefficients indicate a more efficient rainfall–runoff response, either due to the activation of stored water or the fast runoff of rainwater or snowmelt into the stream. Across all studied catchments, the highest runoff coefficients are consistently found for snow-impacted events (Rain.on.snow and Mix) (Fig. 7). High values of runoff coefficients were connected to highly positive residuals, indicating that compared to the CQ relationship, more nitrate was mobilized during high levels of hydrologic connectivity. Studies using high-frequency data show that during runoff events in wet seasons, when catchments are hydrologically more connected, shallow flow paths are activated transporting greater amounts of nitrate (Inamdar et al., 2006; Outram et al., 2016; Schwientek et al., 2013). Similarly, von Freyberg et al. (2014), Ocampo et al. (2006) and Stieglitz et al. (2003) showed that upland zones are more efficiently connected to riparian zones by shallow pathways during wet months, permitting effective transport of nutrients to the stream.

There are no significant differences in event runoff coefficients between two types of snow-impacted events (i.e., Mix and Rain.on.snow) (Fig. 7), despite possible differences in their characteristic snowmelt intensities (Tarasova et al., 2020). While the melting of the snowpack is only induced by temperature increase during Mix events, additional portions of snowpack might be melted by additional energy brought by rainfall during rain-on-snow events (Cohen et al., 2015). However, in both cases, event runoff coefficients are similarly high (Fig. 7b), suggesting similar hydrologic connectivity (Blume and van Meerveld, 2015). This, in turn, results in similar deviations of nitrate concentrations from the long-term CQ relationships for these two snow-impacted event types, reemphasizing the primary role of hydrologic connectivity on the observed deviations of CQ relationships.

In contrast to snow-impacted events, lower runoff coefficients typical for rainfall events with dry antecedent conditions (i.e., Rain.dry.patchy and Rain.dry.uniform) indicate that a small portion of event water can reach the stream, which means that distant zones from the stream network associated with longer pathways exhibit lower or lack of connection during these types of runoff events. Moreover, the overall dry antecedent conditions with a heterogeneous spatial distribution of soil moisture indicate a potential disconnection of runoff generation zones, and therefore the hydrologic connectivity might be lower during such events (i.e., Rain.dry.patchy) than during events with uniform spatial distribution of soil moisture (i.e., Rain.dry.uniform) as shown by lower event runoff coefficients of the former (Fig. 7). Differences in the connectivity of these two types of events are also in line with differences in residuals with the strongest negative residuals observed for Rain.dry.patchy events. Similarly, Outram et al. (2016) showed that lower event runoff coefficients during runoff events with dry antecedent conditions mobilize only a small quantity of nitrate due to the inactive subsurface pathways. Studies using high-frequency nitrate and discharge data showed that during dry periods, upland and riparian zones are usually disconnected (Ocampo et al., 2006; Stieglitz et al., 2003; von Freyberg et al., 2014). This evidence from single catchments are in line with our results across a large set of German catchments, suggesting a crucial role of hydrologic connectivity for nutrient transport.

During Rain.wet events, runoff coefficients are between those of snow-impacted events and rainfall events with dry antecedent conditions (Fig. 7) which indicate an intermediate level of hydrologic connectivity between sources and streams, and thus both positive and negative CQ deviations (Fig. 5).

4.3 Climatic and landscape controls of the variability of CQ deviations across Germany

Despite systematic differences of CQ deviations for different event types, we found considerable spatial variability in the magnitude of these deviations across German catchments (Fig. 5a). In the next paragraphs, we discuss how catchment characteristics may control deviations of concentrations taken during events from the long-term CQ relationship (Δres50).

The correlation of CQ deviations during snow-impacted events with topographic and soil properties (Fig. 6) indicates that in flatter catchments with thick soils and a high fraction of sedimentary aquifers, these types of events generate higher nitrate concentrations compared to the average long-term behavior. Contrarily, CQ deviations during Rain.dry.patchy events are more negative, indicating that nitrate concentrations during these events tend to be smaller than the long-term average in catchments with these characteristics. Previous studies have shown how these characteristics are able to promote nitrate removal in catchments. Deep sedimentary aquifers have a high potential of denitrification due to a great availability of electron donors, longer transit times and more anoxic conditions due to sufficient reduction capacity (Kunkel et al., 2004; Wendland et al., 2008; Knoll et al., 2020), generating a lower nitrate supply in deeper soils compared to shallow soil (Dupas et al., 2016). In addition, flat catchments (low topographic slope, higher topographic wetness index) tend to have a higher portion of riparian wetlands (Musolff et al., 2018) that can reduce nitrate concentrations from stream water. During snow-impacted events, fast flow pathways between nitrate sources and the stream are activated and nitrate stored in shallow soils can be mobilized, bypassing the denitrification attenuation of the soil and the riparian zone, which is also largely suppressed during low temperatures in winter (Johannsen et al., 2008; Lutz et al., 2020), resulting in higher nitrate CQ deviations. In contrast, Rain.dry.patchy events might mobilize water from connected near-stream source zones, where riparian wetlands from flatter areas contribute water with low nitrate concentration (Fovet et al., 2018; Lutz et al., 2020), generating more negative CQ deviations. In addition, deviations during these events might be associated with longer transit times due to thicker soil and less hydrologic connectivity (Yang et al., 2018) which can reduce nitrate concentrations in streams. Instead, steeper catchments with shallow soils during Rain.dry.patchy events show less nitrate attenuation due to shorter flow paths and less favorable conditions for denitrification, generating relatively higher stream water nitrate concentrations during these events and therefore decreasing the magnitude of CQ deviations.

We acknowledge that some catchment characteristics are highly correlated (Fig. S6). Flatter catchments often exhibit higher fractions of agriculture, therefore more diffuse source availability. Although the correlation of the fraction of agriculture and CQ deviations during Rain.on.snow events was less significant than topographic descriptors, a potential increment of diffuse sources in flatter catchments might also enhance the mechanism of nitrate bypassing the buffer capacity of catchments during Rain.on.snow events, generating higher CQ deviations. For most of the event types, we found that the fraction of agriculture itself is not sufficient to explain the differences in nitrate deviations from the long-term CQ relationships between catchments (Fig. 6). However, the vertical ratio of nitrate between topsoil and groundwater and the horizontal spatial distribution of agricultural land within the catchment (i.e., horizontal heterogeneity) were strongly correlated with CQ deviations for Rain.dry.patchy events (Fig. 6). During Rain.dry.patchy events, the mobilization of distant nitrate sources (horizontally and vertically) is reduced due to the low hydrologic connectivity, resulting in lower nitrate concentration of stream water and more negative deviations in catchments with top-loaded nitrate profiles as well as more distant agricultural lands from streams. The high spatial variability of agricultural nitrate sources, expressed as horizontal heterogeneity and vertical ratio of nitrate, and the temporal variability of sources possibly induced by elevated subsurface and riparian zone removal during different levels of hydrologic connectivity, promote deviations of nitrate concentrations from the long-term CQ relationships.

4.4 Implications of this study

In this study, we performed the first large-scale analysis of long-term nitrate CQ relationships, differentiating runoff event types. We show that flatter catchments with soil conditions favorable for denitrification or distant nitrate sources are prone to generate disproportional loads during runoff events with high levels of hydrologic connectivity, presenting an ecological risk for aquatic ecosystems. These findings can be instructive for implementing more effective water quality management strategies to prevent extreme nitrate loads from reaching water bodies in such catchments during events associated with high levels of hydrologic connectivity (i.e., snow-impacted events).

The connection between nitrate concentrations and different types of runoff events shown in our study indicates that possible changes in the occurrences of different event types due to the ongoing climate change might in turn affect the dynamics of nutrient exports in the catchments. With advancing climate change, air temperature is projected to increase further, leading to a substantial decline in seasonal snowpack accumulation and earlier snowmelt onset in Central Europe (Arias et al., 2021). Several studies reported a reduction in snow accumulation in Germany over the last decades (Fontrodona Bach et al., 2018; Chan et al., 2020; Taszarek et al., 2020), with a consistent reduction in the frequency of Rain.on.snow events (Cohen et al., 2015), suggesting that the corresponding positive deviations from the long-term nitrate CQ relationships are likely to occur less often in the future. Less frequent snow-impacted events would reduce nitrate mobilization from the soil under these critical event conditions. Consequently, more nitrate may remain in the soil sources. A fraction of this soil nitrate is expected to be removed by denitrification, whereas another fraction may last longer as soil nitrate legacy (Dupas et al., 2020; Meter et al., 2016), thus generating unknown long-term effects in the nitrate dynamics during future runoff events. On the other hand, higher temperatures lead to a decrease of soil moisture (Dai et al., 2004), propitiating dry conditions and reducing hydrologic connectivity. An increase in frequency of rainfall events with dry antecedent conditions observed in several German catchments (Winter et al., 2022) indicates that negative deviations might become even more frequent during warm seasons in the future.

By using low-frequency, long-term nitrate data, we were able to provide information about characteristic nitrate transport during different types of events and identify hydrologic connectivity associated with these types as a critical control of nitrate dynamics in German catchments. Our findings using low-frequency data are largely supported by the detailed analysis of high-frequency data in individual catchments from previous studies; but thanks to the large number of analyzed catchments, we were able to provide a more comprehensive analysis of systematic deviations of nitrate concentrations during events of different types and provide valuable insights on the origins of the scatter in CQ relationships. The abundance of low-frequency data worldwide and transferable nature of the applied event classification framework provide the means of further applications in contrasting environments to better understand long-term nitrate CQ relationships across contrasting environments. Moreover, our results suggest that sampling campaigns should be designed specifically to capture runoff events with different levels of hydrologic connectivity in order to better explain the scatter in long-term CQ relationships and better isolate the role of singular processes (i.e., nitrate uptake, denitrification).

Although the presence of the event-scale hysteresis effect might considerably affect nitrate concentration during rising and falling limbs of the event hydrograph in some catchments (Pohle et al., 2021), we found a similar direction of deviations from the long-term CQ relationships when we considered samples taken during rising limb, falling limb and near-to-peak (Fig. S6b). Hence, our results suggest that the variability, potentially added by the presence of hysteresis patterns, is lower than the deviations observed for different event types from the long-term CQ relationships. Increasing availability of high-frequency datasets coupled with new statistical modeling approaches might be used in the future to evaluate hysteresis-related effects in the existing long-term CQ datasets to further disentangle inter- and intra-event variability of nitrate dynamics at larger scales.

5 Conclusions

We analyzed for the first time the effect of different runoff event types on the scatter observed in concentration–discharge (CQ) relationships across 184 German catchments. Specifically, we examined the deviations of the concentration of nitrate samples collected during different runoff event types from the long-term CQ relationships. Our results highlight pronounced deviations in most of the catchments, regardless of their overall long-term CQ export patterns (dilution, neutral, or enrichment). Thus, scatter apparent in long-term CQ relationships can indeed be partially explained by different types of runoff event conditions.

We found that nitrate transport is enhanced during snow-impacted events compared to long-term CQ relationships. On the other hand, nitrate concentrations tend to be lower than the long-term CQ relationships when rainfall coincides with dry antecedent conditions. The CQ relationships during rainfall on wet antecedent conditions were not significantly different from the long-term relationships. We argue that hydrologic connectivity to the nitrate sources, here represented by the values of event runoff coefficients, is crucial to explain deviations from the long-term CQ relationship during different event types.

Finally, we found that flatter catchments with high denitrification potential (i.e., deep soils, presence of sedimentary aquifers), as well as catchments with agricultural areas located farther from the stream or with top-loaded nitrate profiles, exhibit an enhanced nitrate transport during snow-impacted events and lower nitrate concentrations during events induced by rainfall with dry antecedent conditions compared to the long-term CQ relationships. Catchments with these characteristics are prone to generate disproportional loads during snow-impacted events, exacerbating ecological risk for receiving water bodies. Findings from this study improve our understanding of the effects of runoff event types on nutrient dynamics and provide valuable insights for optimizing water quality management and monitoring.

Data availability

Water quality dataset is available online at: https://doi.org/10.4211/hs.a42addcbd59a466a9aa56472dfef8721 (Musolff, 2020). Runoff event classification is available upon request at: https://doi.org/10.5281/zenodo.3575024 (Tarasova, 2019). Catchment characteristics dataset is available at: https://doi.org/10.4211/hs.82f8094dd61e449a826afdef820a2c19 (Ebeling et al., 2021).

Supplement

The supplement related to this article is available online at: https://doi.org/10.5194/hess-26-6227-2022-supplement.

Author contributions

FAS contributed with formal analyses. LT and AM contributed to the conceptualization. FAS prepared the manuscript with contributions from all co-authors.

Competing interests

The contact author has declared that none of the authors has any competing interests.

Disclaimer

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Acknowledgements

We thank Camile Minaudo and the two anonymous reviewers for their valuable suggestions that helped to improve the original manuscript. Larisa Tarasova was supported by the German Research Foundation (“Deutsche Forschungsgemeinschaft”, DFG) in terms of the research group FOR 2416 “Space-Time Dynamics of Extreme Floods (SPATE)”. Jana von Freyberg was supported by the Swiss National Science Foundation SNSF (grant PR00P2_185931).

Financial support

This research has been supported by the Helmholtz-Zentrum für Umweltforschung (DYNAMO Cohort).

The article processing charges for this open-access publication were covered by the Helmholtz Centre for Environmental Research – UFZ.

Review statement

This paper was edited by Genevieve Ali and reviewed by Camille Minaudo and two anonymous referees.

References

Arias, P., Bellouin, N., Coppola, E., Jones, R., Krinner, G., Marotzke, J., Naik, V., Palmer, M., Plattner, G. K., Rogelj, J., and Rojas, M.: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, Technical Summary, https://www.ipcc.ch/report/ar6/wg1/ (last access: 8 December 2022), 2021. 

Basu, N. B., Thompson, S. E., and Rao, P. S. C.: Hydrologic and biogeochemical functioning of intensively managed catchments: A synthesis of top-down analyses, Water Resour. Res., 47, W00J15, https://doi.org/10.1029/2011WR010800, 2011. 

Bauwe, A., Tiemeyer, B., Kahle, P., and Lennartz, B.: Classifying hydrological events to quantify their impact on nitrate leaching across three spatial scales, J. Hydrol., 531, 589–601, https://doi.org/10.1016/j.jhydrol.2015.10.069, 2015. 

Benettin, P., Fovet, O., and Li, L.: Nitrate removal and young stream water fractions at the catchment scale, Hydrol. Process., 34, 2725–2738, https://doi.org/10.1002/hyp.13781, 2020. 

Bieroza, M. Z., Heathwaite, A. L., Bechmann, M., Kyllmar, K., and Jordan, P.: The concentration–discharge slope as a tool for water quality management, Sci. Total Environ., 630, 738–749, https://doi.org/10.1016/j.scitotenv.2018.02.256, 2018. 

Blume, T. and van Meerveld, H. J. (Ilja): From hillslope to stream: methods to investigate subsurface connectivity, WIREs Water, 2, 177–198, https://doi.org/10.1002/wat2.1071, 2015. 

Blume, T., Zehe, E., and Bronstert, A.: Rainfall—runoff response, event-based runoff coefficients and hydrograph separation, Hydrolog. Sci. J., 52, 843–862, https://doi.org/10.1623/hysj.52.5.843, 2007. 

Bowes, M. J., Jarvie, H. P., Naden, P. S., Old, G. H., Scarlett, P. M., Roberts, C., Armstrong, L. K., Harman, S. A., Wickham, H. D., and Collins, A. L.: Identifying priorities for nutrient mitigation using river concentration–flow relationships: The Thames basin, UK, J. Hydrol., 517, 1–12, https://doi.org/10.1016/j.jhydrol.2014.03.063, 2014. 

Bowes, M. J., Jarvie, H. P., Halliday, S. J., Skeffington, R. A., Wade, A. J., Loewenthal, M., Gozzard, E., Newman, J. R., and Palmer-Felgate, E. J.: Characterising phosphorus and nitrate inputs to a rural river using high-frequency concentration–flow relationships, Sci. Total Environ., 511, 608–620, https://doi.org/10.1016/j.scitotenv.2014.12.086, 2015. 

Bracken, L. J., Wainwright, J., Ali, G. A., Tetzlaff, D., Smith, M. W., Reaney, S. M., and Roy, A. G.: Concepts of hydrological connectivity: Research approaches, pathways and future agendas, Earth-Sci. Rev., 119, 17–34, https://doi.org/10.1016/j.earscirev.2013.02.001, 2013. 

Branco, P., Torgo, L., and Ribeiro, R.: A Survey of Predictive Modelling under Imbalanced Distributions, arXiv [cs], arXiv:1505.01658, 2015. 

Butturini, A., Gallart, F., Latron, J., Vazquez, E., and Sabater, F.: Cross-site Comparison of Variability of DOC and Nitrate cq Hysteresis during the Autumn–winter Period in Three Mediterranean Headwater Streams: A Synthetic Approach, Biogeochemistry, 77, 327–349, https://doi.org/10.1007/s10533-005-0711-7, 2006. 

Cartwright, I.: Concentration vs. streamflow (C-Q) relationships of major ions in south-eastern Australian rivers: Sources and fluxes of inorganic ions and nutrients, Appl. Geochem., 120, 104680, https://doi.org/10.1016/j.apgeochem.2020.104680, 2020. 

Casquin, A., Dupas, R., Gu, S., Couic, E., Gruau, G., and Durand, P.: The influence of landscape spatial configuration on nitrogen and phosphorus exports in agricultural catchments, Landscape Ecol., 36, 3383–3399, https://doi.org/10.1007/s10980-021-01308-5, 2021. 

Casson, N. J., Eimers, M. C., and Watmough, S. A.: Sources of nitrate export during rain-on-snow events at forested catchments, Biogeochemistry, 120, 23–36, https://doi.org/10.1007/s10533-013-9850-4, 2014. 

Chan, S. C., Kendon, E. J., Berthou, S., Fosser, G., Lewis, E., and Fowler, H. J.: Europe-wide precipitation projections at convection permitting scale with the Unified Model, Clim. Dynam., 55, 409–428, https://doi.org/10.1007/s00382-020-05192-8, 2020. 

Chang, S., Zhang, Q., Byrnes, D., Basu, N., and Van Meter, K.: Chesapeake legacies: The importance of legacy nitrogen to improving Chesapeake Bay water quality, Environ. Res. Lett., 16, 085002, https://doi.org/10.1088/1748-9326/ac0d7b, 2021. 

Chen, X., Parajka, J., Széles, B., Valent, P., Viglione, A., and Blöschl, G.: Impact of Climate and Geology on Event Runoff Characteristics at the Regional Scale, Water, 12, 3457, https://doi.org/10.3390/w12123457, 2020. 

Cohen, J., Ye, H., and Jones, J.: Trends and variability in rain-on-snow events, Geophys. Res. Lett., 42, 7115–7122, https://doi.org/10.1002/2015GL065320, 2015. 

Cole, L. J., Stockan, J., and Helliwell, R.: Managing riparian buffer strips to optimise ecosystem services: A review, Agriculture, Ecosystems and Environment, 296, 106891, https://doi.org/10.1016/j.agee.2020.106891, 2020. 

Covino, T.: Hydrologic connectivity as a framework for understanding biogeochemical flux through watersheds and along fluvial networks, Geomorphology, 277, 133–144, https://doi.org/10.1016/j.geomorph.2016.09.030, 2017. 

Curtin, D., Beare, M. H., and Hernandez-Ramirez, G.: Temperature and Moisture Effects on Microbial Biomass and Soil Organic Matter Mineralization, Soil Sci. Soc. Am. J., 76, 2055–2067, https://doi.org/10.2136/sssaj2012.0011, 2012. 

Dai, A., Trenberth, K. E., and Qian, T.: A Global Dataset of Palmer Drought Severity Index for 1870–2002: Relationship with Soil Moisture and Effects of Surface Warming, J. Hydrometeorol., 5, 1117–1130, https://doi.org/10.1175/JHM-386.1, 2004. 

Dehaspe, J., Sarrazin, F., Kumar, R., Fleckenstein, J. H., and Musolff, A.: Bending of the concentration discharge relationship can inform about in-stream nitrate removal, Hydrol. Earth Syst. Sci., 25, 6437–6463, https://doi.org/10.5194/hess-25-6437-2021, 2021. 

Diamond, J. S. and Cohen, M. J.: Complex patterns of catchment solute–discharge relationships for coastal plain rivers, Hydrol. Process., 32, 388–401, https://doi.org/10.1002/hyp.11424, 2018. 

Dupas, R., Jomaa, S., Musolff, A., Borchardt, D., and Rode, M.: Disentangling the influence of hydroclimatic patterns and agricultural management on river nitrate dynamics from sub-hourly to decadal time scales, Sci. Total Environ., 571, 791–800, https://doi.org/10.1016/j.scitotenv.2016.07.053, 2016. 

Dupas, R., Abbott, B. W., Minaudo, C., and Fovet, O.: Distribution of Landscape Units Within Catchments Influences Nutrient Export Dynamics, Frontiers in Environmental Science, 7, 43, https://doi.org/10.3389/fenvs.2019.00043, 2019. 

Dupas, R., Ehrhardt, S., Musolff, A., Fovet, O., and Durand, P.: Long-term nitrogen retention and transit time distribution in agricultural catchments in western France, Environ. Res. Lett., 15, 115011, https://doi.org/10.1088/1748-9326/abbe47, 2020. 

Ebeling, P., Kumar, R., and Musolff, A.: CCDB – catchment characteristics data base Germany, HydroShare [data set], https://doi.org/10.4211/hs.82f8094dd61e449a826afdef820a2c19, 2021. 

Ebeling, P., Kumar, R., Weber, M., Knoll, L., Fleckenstein, J. H., and Musolff, A.: Archetypes and Controls of Riverine Nutrient Export Across German Catchments, Water Resour. Res., 57, e2020WR028134, https://doi.org/10.1029/2020WR028134, 2021. 

EEA: The European environment-State and outlook 2020, in: Knowledge for Transition to a Sustainable Europe; Publications Office of the European Union, Luxembourg, Brussels, Belgium, https://doi.org/10.2800/96749, 2019. 

Fang, Z., Carroll, R. W. H., Schumer, R., Harman, C., Wilusz, D., and Williams, K. H.: Streamflow partitioning and transit time distribution in snow-dominated basins as a function of climate, J. Hydrol., 570, 726–738, https://doi.org/10.1016/j.jhydrol.2019.01.029, 2019. 

Fontrodona Bach, A., van der Schrier, G., Melsen, L. A., Klein Tank, A. M. G., and Teuling, A. J.: Widespread and Accelerated Decrease of Observed Mean and Extreme Snow Depth Over Europe, Geophys. Res. Lett., 45, 12,312-12,319, https://doi.org/10.1029/2018GL079799, 2018. 

Fovet, O., Humbert, G., Dupas, R., Gascuel-Odoux, C., Gruau, G., Jaffrezic, A., Thelusma, G., Faucheux, M., Gilliet, N., Hamon, Y., and Grimaldi, C.: Seasonal variability of stream water quality response to storm events captured using high-frequency and multi-parameter data, J. Hydrol., 559, 282–293, https://doi.org/10.1016/j.jhydrol.2018.02.040, 2018. 

GEA: Waters in Germany: Status and assessment, German Environment Agency, Dessau–Roßlau, 2017. 

Godsey, S. E., Kirchner, J. W., and Clow, D. W.: Concentration–discharge relationships reflect chemostatic characteristics of US catchments, Hydrol. Process., 23, 1844–1864, https://doi.org/10.1002/hyp.7315, 2009. 

Gorski, G. and Zimmer, M. A.: Hydrologic regimes drive nitrate export behavior in human-impacted watersheds, Hydrol. Earth Syst. Sci., 25, 1333–1345, https://doi.org/10.5194/hess-25-1333-2021, 2021. 

Grayson, R. B., Western, A. W., Chiew, F. H. S., and Blöschl, G.: Preferred states in spatial soil moisture patterns: Local and nonlocal controls, Water Resour. Res., 33, 2897–2908, https://doi.org/10.1029/97WR02174, 1997. 

Guillemot, S., Fovet, O., Gascuel-Odoux, C., Gruau, G., Casquin, A., Curie, F., Minaudo, C., Strohmenger, L., and Moatar, F.: Spatio-temporal controls of C–N–P dynamics across headwater catchments of a temperate agricultural region from public data analysis, Hydrol. Earth Syst. Sci., 25, 2491–2511, https://doi.org/10.5194/hess-25-2491-2021, 2021. 

Guntiñas, M. E., Leirós, M. C., Trasar-Cepeda, C., and Gil-Sotres, F.: Effects of moisture and temperature on net soil nitrogen mineralization: A laboratory study, Eur. J. Soil Biol., 48, 73–80, https://doi.org/10.1016/j.ejsobi.2011.07.015, 2012. 

Hardie, M. A., Cotching, W. E., Doyle, R. B., Holz, G., Lisson, S., and Mattern, K.: Effect of antecedent soil moisture on preferential flow in a texture-contrast soil, J. Hydrol., 398, 191–201, https://doi.org/10.1016/j.jhydrol.2010.12.008, 2011. 

Häussermann, U., Klement, L., Breuer, L., Ullrich, A., Wechsung, G., and Bach, M.: Nitrogen soil surface budgets for districts in Germany 1995 to 2017, Environmental Sciences Europe, 32, 109, https://doi.org/10.1186/s12302-020-00382-x, 2020. 

Heathwaite, A. L. and Bieroza, M.: Fingerprinting hydrological and biogeochemical drivers of freshwater quality, Hydrol. Process., 35, e13973, https://doi.org/10.1002/hyp.13973, 2021. 

House, W. A., Leach, D. V., and Armitage, P. D.: Study of dissolved silicon, and nitrate dynamics in a fresh water stream, Water Res., 35, 2749–2757, https://doi.org/10.1016/S0043-1354(00)00548-0, 2001. 

Inamdar, S. P., O'Leary, N., Mitchell, M. J., and Riley, J. T.: The impact of storm events on solute exports from a glaciated forested watershed in western New York, USA, Hydrol. Process., 20, 3423–3439, https://doi.org/10.1002/hyp.6141, 2006. 

Jencso, K. G., McGlynn, B. L., Gooseff, M. N., Wondzell, S. M., Bencala, K. E., and Marshall, L. A.: Hydrologic connectivity between landscapes and streams: Transferring reach- and plot-scale understanding to the catchment scale, Water Resour. Res., 45, W04428, https://doi.org/10.1029/2008WR007225, 2009. 

Johannsen, A., Dähnke, K., and Emeis, K.: Isotopic composition of nitrate in five German rivers discharging into the North Sea, Org. Geochem., 39, 1678–1689, https://doi.org/10.1016/j.orggeochem.2008.03.004, 2008. 

Knapp, J. L. A., von Freyberg, J., Studer, B., Kiewiet, L., and Kirchner, J. W.: Concentration–discharge relationships vary among hydrological events, reflecting differences in event characteristics, Hydrol. Earth Syst. Sci., 24, 2561–2576, https://doi.org/10.5194/hess-24-2561-2020, 2020. 

Knoll, L., Breuer, L., and Bach, M.: Nation-wide estimation of groundwater redox conditions and nitrate concentrations through machine learning, Environ. Res. Lett., 15, 064004, https://doi.org/10.1088/1748-9326/ab7d5c, 2020. 

Koenig, L. E., Shattuck, M. D., Snyder, L. E., Potter, J. D., and McDowell, W. H.: Deconstructing the Effects of Flow on DOC, Nitrate, and Major Ion Interactions Using a High-Frequency Aquatic Sensor Network, Water Resour. Res., 53, 10655–10673, https://doi.org/10.1002/2017WR020739, 2017. 

Korom, S. F., Schuh, W. M., Tesfay, T., and Spencer, E. J.: Aquifer denitrification and in situ mesocosms: Modeling electron donor contributions and measuring rates, J. Hydrol., 432–433, 112–126, https://doi.org/10.1016/j.jhydrol.2012.02.023, 2012. 

Kruskal, W. H. and Wallis, W. A.: Use of Ranks in One-Criterion Variance Analysis, J. Am. Stat. Assoc., 47, 583–621, https://doi.org/10.1080/01621459.1952.10483441, 1952. 

Kumar, R., Livneh, B., and Samaniego, L.: Toward computationally efficient large-scale hydrologic predictions with a multiscale regionalization scheme, Water Resour. Res., 49, 5700–5714, https://doi.org/10.1002/wrcr.20431, 2013. 

Kunkel, R., Bach, M., Behrendt, H., and Wendland, F.: Groundwater-borne nitrate intakes into surface waters in Germany, Water Sci. Technol., 49, 11–19, https://doi.org/10.2166/wst.2004.0152, 2004. 

Lassaletta, L., Billen, G., Grizzetti, B., Anglade, J., and Garnier, J.: 50 year trends in nitrogen use efficiency of world cropping systems: the relationship between yield and nitrogen input to cropland, Environ. Res. Lett., 9, 105011, https://doi.org/10.1088/1748-9326/9/10/105011, 2014. 

Lloyd, C. E. M., Freer, J. E., Johnes, P. J., and Collins, A. L.: Using hysteresis analysis of high-resolution water quality monitoring data, including uncertainty, to infer controls on nutrient and sediment transfer in catchments, Sci. Total Environ., 543, 388–404, https://doi.org/10.1016/j.scitotenv.2015.11.028, 2016. 

Lutz, S. R., Trauth, N., Musolff, A., Van Breukelen, B. M., Knöller, K., and Fleckenstein, J. H.: How Important is Denitrification in Riparian Zones? Combining End-Member Mixing and Isotope Modeling to Quantify Nitrate Removal from Riparian Groundwater, Water Resour. Res., 56, e2019WR025528, https://doi.org/10.1029/2019WR025528, 2020. 

Martin, C., Aquilina, L., Gascuel-Odoux, C., Molénat, J., Faucheux, M., and Ruiz, L.: Seasonal and interannual variations of nitrate and chloride in stream waters related to spatial and temporal patterns of groundwater concentrations in agricultural catchments, Hydrol. Process., 18, 1237–1254, https://doi.org/10.1002/hyp.1395, 2004. 

McGlynn, B. L. and Seibert, J.: Distributed assessment of contributing area and riparian buffering along stream networks, Water Resour. Res., 39, 1082, https://doi.org/10.1029/2002WR001521, 2003. 

Meter, K. J. V. and Basu, N. B.: Time lags in watershed-scale nutrient transport: an exploration of dominant controls, Environ. Res. Lett., 12, 084017, https://doi.org/10.1088/1748-9326/aa7bf4, 2017. 

Meter, K. J. V., Basu, N. B., Veenstra, J. J., and Burras, C. L.: The nitrogen legacy: emerging evidence of nitrogen accumulation in anthropogenic landscapes, Environ. Res. Lett., 11, 035014, https://doi.org/10.1088/1748-9326/11/3/035014, 2016. 

Meybeck, M. and Moatar, F.: Daily variability of river concentrations and fluxes: indicators based on the segmentation of the rating curve, Hydrol. Process., 26, 1188–1207, https://doi.org/10.1002/hyp.8211, 2012. 

Minaudo, C., Dupas, R., Gascuel-Odoux, C., Roubeix, V., Danis, P.-A., and Moatar, F.: Seasonal and event-based concentration–discharge relationships to identify catchment controls on nutrient export regimes, Adv. Water Resour., 131, 103379, https://doi.org/10.1016/j.advwatres.2019.103379, 2019. 

Moatar, F., Abbott, B. W., Minaudo, C., Curie, F., and Pinay, G.: Elemental properties, hydrology, and biology interact to shape concentration-discharge curves for carbon, nutrients, sediment, and major ions, Water Resour. Res., 53, 1270–1287, https://doi.org/10.1002/2016WR019635, 2017. 

Mulholland, P. J., Helton, A. M., Poole, G. C., Hall, R. O., Hamilton, S. K., Peterson, B. J., Tank, J. L., Ashkenas, L. R., Cooper, L. W., Dahm, C. N., Dodds, W. K., Findlay, S. E. G., Gregory, S. V., Grimm, N. B., Johnson, S. L., McDowell, W. H., Meyer, J. L., Valett, H. M., Webster, J. R., Arango, C. P., Beaulieu, J. J., Bernot, M. J., Burgin, A. J., Crenshaw, C. L., Johnson, L. T., Niederlehner, B. R., O'Brien, J. M., Potter, J. D., Sheibley, R. W., Sobota, D. J., and Thomas, S. M.: Stream denitrification across biomes and its response to anthropogenic nitrate loading, Nature, 452, 202–205, https://doi.org/10.1038/nature06686, 2008. 

Musolff, A. WQQDB – water quality and quantity data base Germany, HydroShare [data set], https://doi.org/10.4211/hs.a42addcbd59a466a9aa56472dfef8721, 2020. 

Musolff, A., Fleckenstein, J. H., Rao, P. S. C., and Jawitz, J. W.: Emergent archetype patterns of coupled hydrologic and biogeochemical responses in catchments: Emergence of Archetype CQ Patterns, Geophys. Res. Lett., 44, 4143–4151, https://doi.org/10.1002/2017GL072630, 2017. 

Musolff, A., Fleckenstein, J., Opitz, M., Büttner, O., Kumar, R., and Tittel, J.: Spatio-temporal controls of dissolved organic carbon stream water concentrations, J. Hydrol., 566, 205–215, https://doi.org/10.1016/j.jhydrol.2018.09.011, 2018. 

Musolff, A., Zhan, Q., Dupas, R., Minaudo, C., Fleckenstein, J., Rode, M., Dehaspe, J., and Rinke, K.: Spatial and Temporal Variability in Concentration–Discharge Relationships at the Event Scale, Water Resour. Res., 57, https://doi.org/10.1029/2020WR029442, 2021. 

Ocampo, C. J., Sivapalan, M., and Oldham, C.: Hydrological connectivity of upland-riparian zones in agricultural catchments: Implications for runoff generation and nitrate transport, J. Hydrol., 331, 643–658, https://doi.org/10.1016/j.jhydrol.2006.06.010, 2006. 

Ortmeyer, F., Begerow, D., Guerreiro, M. A., Wohnlich, S., and Banning, A.: Comparison of Denitrification Induced by Various Organic Substances—Reaction Rates, Microbiology, and Temperature Effect, Water Resour. Res., 57, https://doi.org/10.1029/2021WR029793, 2021. 

Outram, F. N., Cooper, R. J., Sünnenberg, G., Hiscock, K. M., and Lovett, A. A.: Antecedent conditions, hydrological connectivity and anthropogenic inputs: Factors affecting nitrate and phosphorus transfers to agricultural headwater streams, Sci. Total Environ., 545–546, 184–199, https://doi.org/10.1016/j.scitotenv.2015.12.025, 2016. 

Paerl, H. W.: Coastal eutrophication and harmful algal blooms: Importance of atmospheric deposition and groundwater as “new” nitrogen and other nutrient sources, Limnol. Oceanogr., 42, 1154–1165, https://doi.org/10.4319/lo.1997.42.5_part_2.1154, 1997. 

Pohle, I., Baggaley, N., Palarea-Albaladejo, J., Stutter, M., and Glendell, M.: A Framework for Assessing Concentration‐Discharge Catchment Behavior From Low‐Frequency Water Quality Data, Water Resour. Res., 57, e2021WR029692, https://doi.org/10.1029/2021WR029692, 2021. 

Puckett, L. J., Tesoriero, A. J., and Dubrovsky, N. M.: Nitrogen Contamination of Surficial Aquifers—A Growing Legacy, Environ. Sci. Technol., 45, 839–844, https://doi.org/10.1021/es1038358, 2011. 

Rauthe, M., Steiner, H., Riediger, U., Mazurkiewicz, A., and Gratzki, A.: A Central European precipitation climatology – Part I: Generation and validation of a high-resolution gridded daily data set (HYRAS), Meteorol. Z., 235–256, https://doi.org/10.1127/0941-2948/2013/0436, 2013. 

Rode, M., Halbedel née Angelstein, S., Anis, M. R., Borchardt, D., and Weitere, M.: Continuous In-Stream Assimilatory Nitrate Uptake from High-Frequency Sensor Measurements, Environ. Sci. Technol., 50, 5685–5694, https://doi.org/10.1021/acs.est.6b00943, 2016. 

Rose, L. A., Karwan, D. L., and Godsey, S. E.: Concentration–discharge relationships describe solute and sediment mobilization, reaction, and transport at event and longer timescales, Hydrol. Process., 32, 2829–2844, https://doi.org/10.1002/hyp.13235, 2018. 

Sabater, S., Butturini, A., Clement, J.-C., Burt, T., Dowrick, D., Hefting, M., Matre, V., Pinay, G., Postolache, C., Rzepecki, M., and Sabater, F.: Nitrogen Removal by Riparian Buffers along a European Climatic Gradient: Patterns and Factors of Variation, Ecosystems, 6, 0020–0030, https://doi.org/10.1007/s10021-002-0183-8, 2003. 

Samaniego, L., Kumar, R., and Attinger, S.: Multiscale parameter regionalization of a grid-based hydrologic model at the mesoscale, Water Resour. Res., 46, W05523, https://doi.org/10.1029/2008WR007327, 2010. 

Schwientek, M., Osenbrück, K., and Fleischer, M.: Investigating hydrological drivers of nitrate export dynamics in two agricultural catchments in Germany using high-frequency data series, Environ. Earth Sci., 69, 381–393, https://doi.org/10.1007/s12665-013-2322-2, 2013. 

Seibert, J., Grabs, T., Köhler, S., Laudon, H., Winterdahl, M., and Bishop, K.: Linking soil- and stream-water chemistry based on a Riparian Flow-Concentration Integration Model, Hydrol. Earth Syst. Sci., 13, 2287–2297, https://doi.org/10.5194/hess-13-2287-2009, 2009. 

Stieglitz, M., Shaman, J., McNamara, J., Engel, V., Shanley, J., and Kling, G. W.: An approach to understanding hydrologic connectivity on the hillslope and the implications for nutrient transport, Global Biogeochem. Cy., 17, 1105, https://doi.org/10.1029/2003GB002041, 2003. 

Stumpf, R. P., Johnson, L. T., Wynne, T. T., and Baker, D. B.: Forecasting annual cyanobacterial bloom biomass to inform management decisions in Lake Erie, J. Great Lakes Res., 42, 1174–1183, https://doi.org/10.1016/j.jglr.2016.08.006, 2016. 

Tarasova, L.: Classified runoff events, Zenodo [data set], https://doi.org/10.5281/zenodo.3575024, 2019. 

Tarasova, L., Basso, S., Zink, M., and Merz, R.: Exploring Controls on Rainfall-Runoff Events: 1. Time Series-Based Event Separation and Temporal Dynamics of Event Runoff Response in Germany, Water Resour. Res., 54, 7711–7732, https://doi.org/10.1029/2018WR022587, 2018. 

Tarasova, L., Basso, S., Wendi, D., Viglione, A., Kumar, R., and Merz, R.: A Process-Based Framework to Characterize and Classify Runoff Events: The Event Typology of Germany, Water Resour. Res., 56, e2019WR026951, https://doi.org/10.1029/2019WR026951, 2020. 

Taszarek, M., Kendzierski, S., and Pilguj, N.: Hazardous weather affecting European airports: Climatological estimates of situations with limited visibility, thunderstorm, low-level wind shear and snowfall from ERA5, Weather and Climate Extremes, 28, 100243, https://doi.org/10.1016/j.wace.2020.100243, 2020. 

Tesoriero, A. J., Duff, J. H., Saad, D. A., Spahr, N. E., and Wolock, D. M.: Vulnerability of Streams to Legacy Nitrate Sources, Environ. Sci. Technol., 47, 3623–3629, https://doi.org/10.1021/es305026x, 2013. 

Thompson, S. E., Basu, N. B., Lascurain Jr., J., Aubeneau, A., and Rao, P. S. C.: Relative dominance of hydrologic versus biogeochemical factors on solute export across impact gradients, Water Resour. Res., 47, W00J05, https://doi.org/10.1029/2010WR009605, 2011. 

van Grinsven, H. J. M., ten Berge, H. F. M., Dalgaard, T., Fraters, B., Durand, P., Hart, A., Hofman, G., Jacobsen, B. H., Lalor, S. T. J., Lesschen, J. P., Osterburg, B., Richards, K. G., Techen, A.-K., Vertès, F., Webb, J., and Willems, W. J.: Management, regulation and environmental impacts of nitrogen fertilization in northwestern Europe under the Nitrates Directive; a benchmark study, Biogeosciences, 9, 5143–5160, https://doi.org/10.5194/bg-9-5143-2012, 2012. 

Vaughan, M. C. H., Bowden, W. B., Shanley, J. B., Vermilyea, A., Sleeper, R., Gold, A. J., Pradhanang, S. M., Inamdar, S. P., Levia, D. F., Andres, A. S., Birgand, F., and Schroth, A. W.: High-frequency dissolved organic carbon and nitrate measurements reveal differences in storm hysteresis and loading in relation to land cover and seasonality, Water Resour. Res., 53, 5345–5363, https://doi.org/10.1002/2017WR020491, 2017. 

Veith, T. L., Preisendanz, H. E., and Elkin, K. R.: Characterizing transport of natural and anthropogenic constituents in a long-term agricultural watershed in the northeastern United States, J. Soil Water Conserv., 75, 319–329, https://doi.org/10.2489/jswc.75.3.319, 2020. 

Vervloet, L. S. C., Binning, P. J., Børgesen, C. D., and Højberg, A. L.: Delay in catchment nitrogen load to streams following restrictions on fertilizer application, Sci. Total Environ., 627, 1154–1166, https://doi.org/10.1016/j.scitotenv.2018.01.255, 2018. 

von Freyberg, J., Radny, D., Gall, H. E., and Schirmer, M.: Implications of hydrologic connectivity between hillslopes and riparian zones on streamflow composition, J. Contam. Hydrol., 169, 62–74, https://doi.org/10.1016/j.jconhyd.2014.07.005, 2014. 

Weitere, M., Altenburger, R., Anlanger, C., Baborowski, M., Bärlund, I., Beckers, L.-M., Borchardt, D., Brack, W., Brase, L., Busch, W., Chatzinotas, A., Deutschmann, B., Eligehausen, J., Frank, K., Graeber, D., Griebler, C., Hagemann, J., Herzsprung, P., Hollert, H., Inostroza, P. A., Jäger, C. G., Kallies, R., Kamjunke, N., Karrasch, B., Kaschuba, S., Kaus, A., Klauer, B., Knöller, K., Koschorreck, M., Krauss, M., Kunz, J. V., Kurz, M. J., Liess, M., Mages, M., Müller, C., Muschket, M., Musolff, A., Norf, H., Pöhlein, F., Reiber, L., Risse-Buhl, U., Schramm, K.-W., Schmitt-Jansen, M., Schmitz, M., Strachauer, U., von Tümpling, W., Weber, N., Wild, R., Wolf, C., and Brauns, M.: Disentangling multiple chemical and non-chemical stressors in a lotic ecosystem using a longitudinal approach, Sci. Total Environ., 769, 144324, https://doi.org/10.1016/j.scitotenv.2020.144324, 2021. 

Wendland, F., Blum, A., Coetsiers, M., Gorova, R., Griffioen, J., Grima, J., Hinsby, K., Kunkel, R., Marandi, A., Melo, T., Panagopoulos, A., Pauwels, H., Ruisi, M., Traversa, P., Vermooten, J. S. A., and Walraevens, K.: European aquifer typology: a practical framework for an overview of major groundwater composition at European scale, Environ. Geol., 55, 77–85, https://doi.org/10.1007/s00254-007-0966-5, 2008. 

Winter, C., Lutz, S. R., Musolff, A., Kumar, R., Weber, M., and Fleckenstein, J. H.: Disentangling the Impact of Catchment Heterogeneity on Nitrate Export Dynamics From Event to Long-Term Time Scales, Water Resour. Res., 57, e2020WR027992, https://doi.org/10.1029/2020WR027992, 2021. 

Winter, C., Tarasova, L., Lutz, S., Musolff, A., Kumar, R., and Fleckenstein, J.: Explaining the Variability in High-Frequency Nitrate Export Patterns Using Long-Term Hydrological Event Classification, Water Resour. Res., 58, e2021WR030938, https://doi.org/10.1002/essoar.10507676.1, 2022. 

Yang, J., Heidbüchel, I., Musolff, A., Reinstorf, F., and Fleckenstein, J. H.: Exploring the Dynamics of Transit Times and Subsurface Mixing in a Small Agricultural Catchment, Water Resour. Res., 54, 2317–2335, https://doi.org/10.1002/2017WR021896, 2018. 

Zhang, X., Yang, X., Jomaa, S., and Rode, M.: Analyzing impacts of seasonality and landscape gradient on event-scale nitrate-discharge dynamics based on nested high-frequency monitoring, J. Hydrol., 591, 125585, https://doi.org/10.1016/j.jhydrol.2020.125585, 2020. 

Zhi, W., Li, L., Dong, W., Brown, W., Kaye, J., Steefel, C., and Williams, K. H.: Distinct Source Water Chemistry Shapes Contrasting Concentration-Discharge Patterns, Water Resour. Res., 55, 4233–4251, https://doi.org/10.1029/2018WR024257, 2019.  

Zink, M., Kumar, R., Cuntz, M., and Samaniego, L.: A high-resolution dataset of water fluxes and states for Germany accounting for parametric uncertainty, Hydrol. Earth Syst. Sci., 21, 1769–1790, https://doi.org/10.5194/hess-21-1769-2017, 2017. 

Download
Short summary
Nitrate contamination of rivers from agricultural sources is a challenge for water quality management. During runoff events, different transport paths within the catchment might be activated, generating a variety of responses in nitrate concentration in stream water. Using nitrate samples from 184 German catchments and a runoff event classification, we show that hydrologic connectivity during runoff events is a key control of nitrate transport from catchments to streams in our study domain.