2.1 Site description

The Erlenbach catchment is a small (0.7 km²), steep catchment spanning an elevation range from 1100 to 1655 m above sea level in the northern Swiss pre-Alps. The underlying geologic formation is flysch, and the highly layered bedrock consists of limestone, claystone, marl, and shale, as well as conglomerate and calcareous sandstone (Zobrist, 2010). Groundwater chemistry is thus dominated by calcium, magnesium, and their counter-anion, bicarbonate. The bedrock is overlain by umbric Gleysols with high silt and clay content in the steeper areas (Schleppi et al., 1998; Xu et al., 2009), and by mollic Gleysols with a permanently reduced B_g horizon in the flatter areas (Hagedorn et al., 2000). The catchment landscape is characterized by interchanging slopes and plateaus, and the groundwater table is generally shallower under the plateaus than under the slopes (Rinderer et al., 2014). The soil and bedrock permeabilities are relatively low, resulting in highly saturated soils, particularly on the plateaus. In total, 53 % of the catchment is forested, and dry and wet meadows cover roughly 14 % and 33 % of the catchment area, respectively (van Meerveld et al., 2018). Coniferous forests cover the majority of the slopes, whereas meadows and partially forested areas can be found on the plateaus. Agricultural influence is limited to summer season cattle grazing in the upper part of the catchment. Average annual precipitation in the Erlenbach catchment is 2300 mm, of which up to 40 % falls as snow in the winter months (Stähli and Gustafsson, 2006), and about 20 % of incoming precipitation leaves the catchment as evapotranspiration (van Meerveld et al., 2018). In the years 2017 and 2018, stream discharge at the catchment outlet ranged from 0.2 to 2240 L s⁻¹, with an average value of 37.7 L s⁻¹.

2.2 Dataset

Streamwater and precipitation chemistry were analyzed semicontinuously at an automated field laboratory located at the Erlenbach catchment outlet (von Freyberg et al., 2017, 2018). Streamwater was pumped continuously from the stream to the field laboratory, and precipitation was collected with an open 45 cm diameter funnel. During periods without rain, the field laboratory received only streamwater for analysis. The field laboratory analyzed precipitation whenever more than 50 mL of precipitation accumulated in the rain sampler (corresponding to roughly 0.3 mm of precipitation), and the previously analyzed sample was streamwater. A new sampling and analysis cycle was initiated every 30 min; thus, during rainfall, streamwater and precipitation samples were each analyzed hourly. A drift correction standard was analyzed instead of streamwater every 4 h (every 6 h after 12 March 2018).

Before analysis, each precipitation or streamwater sample passed through a 0.2 µm tangential filter. One aliquot of each sample was automatically directed to an ion chromatograph (940 Professional IC Vario, Metrohm AG, Herisau, Switzerland, hereafter referred to as “IC”) for the analysis of major anions and cations (calcium, magnesium, sodium, potassium, chloride, nitrate, and sulfate). Another aliquot was injected into a continuous water sampler module (CWS, Picarro, Inc., Santa Clara, CA, USA) coupled to a cavity ring-down spectroscope (CRDS, Picarro, Inc., Santa Clara, CA, USA) for water isotope analysis (deuterium and oxygen-18). Further details on the sampling and analysis of isotopes and major ions in the field laboratory are described by von Freyberg et al. (2017, 2018).

In addition to the on-site analysis of major ions and stable water isotopes, aliquots of filtered rainwater (every sample) and streamwater (one sample per hour) were automatically collected into vials. Each vial contained 1 mL of ultrapure HNO₃ to stabilize the water sample. These acidified samples were collected approximately once per week and brought to the laboratory at ETH Zurich for subsequent analysis of a wide range of cations and trace elements using inductively coupled plasma mass spectrometry (Agilent 7900 ICP-MS, Agilent Technologies, Santa Clara, CA, USA, hereafter referred to as “ICP-MS”; more detail on the laboratory protocol can be found in the Supplement). Of the measured elements, we selected boron, strontium, barium, iron, manganese, copper, and chromium for further analysis in this study. Other elements were not included here because their concentrations were mostly below the analytical detection limits. Outlier removal of isotope, IC, and ICP-MS measurements was based on visual inspection.

Streamwater electrical conductivity (EC) was measured at 5 min resolution (s::can condu::lyser, S::CAN Messtechnik GmbH, Vienna, Austria). For the purpose of this study, only every second EC measurement was used to match the 10 min measurement frequency of river discharge (see below).

River discharge was measured at a concrete flume installed at the Erlenbach catchment outlet (Hegg et al., 2006). Precipitation rates were recorded with a heated tipping-bucket rain gauge (Joss-Tognini 15183, LAMBRECHT meteo GmbH, Göttingen, Germany) at a meteorological station located at 1216 m above sea level in the Erlenbach catchment. At the meteorological station, air temperature was measured with a ventilated thermometer (VT3, Meteolabor AG, Wetzikon, Switzerland), and groundwater level fluctuations were recorded in a fully screened well (these readings were relative rather than absolute because they have not been calibrated against manual water level measurements). Discharge, precipitation, air temperature, and groundwater levels were recorded at 10 min resolution. These data were aggregated to 30 min or 1 h intervals to match the frequency of the solute data, except for the analysis of the EC data, for which the 10 min resolution was used. For this purpose, we extracted those data that were closest to the sampling times of the streamwater samples from the 10 min discharge, air temperature, and groundwater level time series. For the precipitation samples, associated precipitation amounts were calculated as cumulative sums from the 10 min tipping-bucket recordings.

2.3 Dataset and event identification

The Erlenbach stream generally shows a fast and flashy response to rainfall events, resulting in very short durations of the rising limb of the storm hydrograph, in particular during short and intense events. Consequently, few samples were collected during the rising limb of each storm in spite of the high sampling frequency of the field laboratory. The falling limb of the streamflow hydrograph, on the other hand, was generally well captured. For this reason, we excluded all samples collected during periods of increasing discharge from our analysis (in Sect. 3.3 below, we show that excluding these periods has negligible effects on the 2-year cQ behavior). The relationship between solute concentrations and streamflow during hydrograph recession is informative in particular on outflow processes from shallow and deeper groundwater, as well as riparian water (Inamdar et al., 2006).

We excluded the winter months from the dataset to avoid ambiguities arising from rain-on-snow and snowfall events in our analysis, because solute responses to these events may differ from those during rain events in snow-free periods. Consequently, only data points between 1 May and 15 November in both 2017 and 2018 were analyzed (hereafter referred to as the “snow-free periods”). During these snow-free periods, we identified individual events based on the following criteria: (1) events had to have a substantial increase in discharge, i.e., peak discharge at least 20 L s⁻¹ above pre-event baseflow; (2) this discharge increase had to be triggered by rainfall; (3) the hydrograph had to recede by at least 75 % of the absolute increase in discharge before the event was cut off (e.g., due to a subsequent rain event or sampler failure); and finally, (4) only events were considered for which IC (i.e., major ions) measurements were available.

For each event, we analyzed the solute concentration and discharge data during the hydrograph recession, starting at the main discharge peak of an event until 95 % hydrograph recession to baseflow, or until the event was cut off by the start of a subsequent event (if discharge had receded by at least 75 %). We estimated cQ relationships for solute/event combinations for which at least five data points were available on the recession limb.

2.4 Assessment of source-area concentrations and input–output budgets

To identify likely streamwater sources, we quantified average solute concentrations of different compartments in the Erlenbach catchment: median solute concentrations in streamwater, precipitation, and their lower and upper quartiles were calculated from the snow-free periods in the time series. Solute concentrations in groundwater were obtained from sampling two pumping wells in the upper part of the catchment at three different instances between 2016 and 2017. Although these groundwater chemistry data may not provide a comprehensive picture of Erlenbach groundwater, they can nonetheless indicate which solutes are likely to be dominant therein. Concentrations in soil water of the neighboring Studibach catchment were obtained from sampling 18 suction lysimeters during eight baseflow snapshot campaigns in the snow-free season in 2016 and 2017 (detailed description of the campaigns are given in Kiewiet et al., 2019). We aggregated the data from six sites (three forested, three non-forested), spread over three elevations in the catchment (1361, 1502 and 1611 m above sea level), at which suction lysimeters were installed at 15, 30 and 50 cm depth. The lysimeters were emptied and set to a tension of 50 mbar the day prior to sampling.

We characterized the source-sink behavior of the Erlenbach catchment for individual solutes using a dimensionless solute flux index calculated from the high-frequency data of the snow-free periods. The solute flux index F relates solute fluxes in precipitation inputs I to those in streamwater outputs O:

where Q_tot and P_tot are the total streamflow and precipitation water fluxes, $c_{Q_i}$ and $c_{P_i}$ are the solute concentrations in streamflow and precipitation for all sampling times i, and the angled brackets indicate volume-weighted averages. The index is positive if the streamwater solute flux is larger than the precipitation solute flux and quantifies the fraction of the output flux O generated within the catchment. Conversely, if the index is negative, it quantifies the fraction of the input flux I retained within the catchment. If the input fluxes and output fluxes are exactly balanced, the flux index will be zero. Importantly, the solute flux indices as calculated here provide no information on long-term fluxes but relate the output fluxes and input fluxes during the snow-free periods between May and November from which the hydrologic events were extracted.

2.5 Quantification of concentration–discharge relationships

We estimated both the 2-year cQ behavior (i.e., using all recession data from the snow-free periods in 2017 and 2018) and the individual event-scale cQ behavior by fitting power-law relationships between concentration c and discharge Q to the data (Clow and Drever, 1996; Musolff et al., 2015):

This power-law relationship is identical to a linear relationship in double-logarithmic space:

where log ₁₀(a) and b are the intercept and slope of the cQ relationship, respectively. The cQ slopes and intercepts of the entire dataset will be referred to as “2-year” cQ slopes and intercepts, whereas “event scale” will refer to cQ relationships of individual hydrologic events. For the purpose of this study, we normalized discharge by the average discharge of the time series (log ₁₀(Q∕Q_mean)) to obtain 2-year intercepts that reflect the expected concentration at the mean discharge, rather than the arbitrary value of $Q=\mathrm{1}({\log }_{\mathrm{10}}Q=\mathrm{0})$. Centering the x axis in this way also has the benefit of making the slope and intercept estimates more independent from one another.

We calculated the relative standard errors of the event-scale slopes and intercepts, in order to exclude events for which the cQ relationships could not be well constrained. The cQ relationships were excluded from our analysis if the relative standard error of either the cQ slope or intercept for any event and solute exceeded 50 %, or if both the relative standard error of the cQ slope and intercept exceeded 25 %.

Given that our samples represent changes in chemistry during hydrograph recession, the meaning of the obtained cQ slopes can be interpreted as follows: a cQ slope of 1 (or −1) is obtained if changes in solute concentrations are proportional (or inversely proportional) to changes in discharge during recession. Consequently, a cQ slope between −1 and 1 indicates less than proportional changes in solute concentrations, and a cQ slope close to zero indicates solute concentrations that change relatively little, or that vary independently of discharge during recession. Conversely, cQ slopes greater than 1 indicate that solute concentrations decrease more than proportionally to discharge during recession. Examples of the relationships between cQ slopes, solute recessions, and hydrograph recessions are illustrated in Fig. 1.

2.6 Possible environmental controls of inter-event variability in concentration–discharge relationships

The 30 hydrologic events span wide ranges of storm durations, intensities, antecedent wetness conditions, and other potential controls, thus facilitating an investigation into how these environmental controls may influence the slopes and intercepts of the event-based cQ relationships. To this end, we quantified 15 different parameters for each event from the following five categories: (1) temperature and proximity to midyear as seasonality indicators; (2) relative input concentrations, which quantify the ratio between the volume-weighted average precipitation concentration during the event to the streamwater solute concentrations during pre-event baseflow; (3) groundwater levels, baseflow discharges, and antecedent precipitation as indicators of antecedent wetness conditions; (4) several measures of event magnitude and intensity; and (5) event and pre-event-water contributions determined from isotope hydrograph separation (following the approach presented by von Freyberg et al., 2018). Table 1 presents an overview of these 15 environmental controls, as well as their ranges in our dataset.

Table 1We assessed 15 environmental controls grouped into five different categories: seasonality indicators, relative input concentrations, antecedent wetness conditions, event characteristics, and event-water contributions. Minimum and maximum values indicate the ranges observed in the dataset. Relative input concentrations are specific for every solute (and not available for all solutes and events), whereas all other controls do not differ between individual solutes. Event-water contributions were only assessed for 22 out of 30 events. Groundwater levels are expressed as negative values so that the maximum corresponds to the wettest conditions, consistent with the other wetness indicators.

Download Print Version | Download XLSX

As a seasonality indicator, the proximity to midyear is calculated as a summer–winter index:

where doy is the day of year. This summer–winter index approaches 0 at the beginning and end of each calendar year, and approaches 1 at the beginning of July.

We used weighted rank correlation coefficients to quantify the dependence of event-scale cQ slopes and intercepts on the 15 environmental controls. The weights were the inverses of the standard errors of the individual cQ slopes and intercepts, to prevent highly uncertain points from substantially influencing the results. The statistical significance of these correlation coefficients (their p values) quantifies the probability of obtaining an equal or greater correlation if the null hypothesis were valid (i.e., if there were actually no relationship between the slope or intercept and the respective control).

3.1 Dataset

We extracted 30 events from the time series of 2017 and 2018 that fulfilled the criteria outlined in Sect. 2.3. While IC measurements were available for all 30 hydrologic events, not all events had the required five sample points for all of the cations analyzed by ICP-MS (i.e., boron, barium, iron, manganese, chromium, strontium, and copper). Furthermore, some events had to be excluded from further analysis for individual solutes due to high relative standard errors of cQ slopes and/or intercepts (see Sect. 2.5). This resulted in 24 to 30 usable events for major ions measured with the IC. More than 20 events were evaluated for all other solutes, except manganese and copper. The concentrations of these two solutes were low and variable, resulting in no clear power-law relationship with discharge for many events. Consequently, only data from 17 and 11 events were usable for manganese and copper, respectively.

We quantified most environmental controls for all 30 events. However, meaningful event-water contributions could only be calculated for 22 events from stable water isotope measurements. Ratios of precipitation concentration to the pre-event baseflow solute concentration were not available for all events and all solutes due to sporadic problems with the rain collector. Also, these ratios were not assessed for EC, because EC was not measured in precipitation. Table 1 provides an overview of the range of environmental controls covered by the selected events. Some small events lasted only a few hours, whereas other events were extended, multi-day storms, and antecedent wetness conditions ranged from relatively dry to very wet. Although all events took place between May and November, they spanned a wide range of air temperatures and weather conditions.

3.2 Characterization of solute contributions from different source areas

Streamwater concentrations at Erlenbach were dominated by calcium, sulfate, magnesium, and sodium, with median concentrations between 2 and 48 mg L⁻¹ (Table 2). Conversely, median streamwater concentrations of iron and copper were low at around 1 to 5 µg L⁻¹, and concentrations of manganese and chromium were even lower. All other solutes (strontium, barium, boron, chloride, potassium, and nitrate) were observed at intermediate concentrations in the Erlenbach streamwater.

Table 2Median solute concentrations and their upper and lower quartiles in streamwater, groundwater, and precipitation in the Erlenbach catchment and in soil water in the adjacent Studibach catchment. Groundwater solute data were collected at one to three different sampling times at two pumping wells located in the upper part of the catchment. These are probably not representative of groundwater concentrations throughout the catchment, but still provide a rough indication of which solutes dominate groundwater. Concentrations in precipitation and streamwater were obtained from the time series recorded at the Erlenbach outlet, excluding months with snow. N indicates the number of samples the calculations are based on. EC was not analyzed in precipitation or soil water. Concentrations greater than 10 µg L⁻¹ are rounded to the nearest integer. Flux index calculations (see Eq. 1) of the snow-free season are based on the same period as the analyzed time series and consequently do not represent long-term fluxes. Positive flux indices quantify the fraction of the output flux that was generated in the catchment, while negative flux indices quantify the fraction of the input flux retained in the catchment. A value of 0 is obtained if input fluxes and output fluxes balance. The ratio of precipitation to streamflow water fluxes during the snow-free period was 1.75 (compared to a ratio of 1.41 for the time from 1 January 2017 to 31 December 2018).

Download Print Version | Download XLSX

We can group different solutes based on their most important sources in the Erlenbach catchment as determined from their concentration ranges in groundwater, streamwater, and precipitation, as well as output–input flux indices (Table 2). The concentrations of calcium, strontium, barium, and boron in the groundwater samples were similar to, or higher than, those in streamwater. Concentrations of magnesium, sodium, and sulfate were higher in streamwater than in the groundwater samples, suggesting that other, unsampled, groundwaters with higher concentrations also contribute to streamflow. In general, groundwater concentrations of weathering products are likely to be highly heterogeneous due to spatially variable contributions from the geochemically complex flysch bedrock (Fischer et al., 2015; Kiewiet et al., 2019). Precipitation concentrations of weathering products are low, and flux indices close to 1, implying that they are primarily derived from within-catchment processes (Table 2). In this paper, we will refer to these solutes as groundwater-sourced.

The output–input flux index of chloride at Erlenbach was 0.33 for the snow-free periods, but was −0.13 if calculated from daily data over the 2 full years (not shown here), suggesting that precipitation is the most important source of chloride in the Erlenbach catchment on an annual basis (Zobrist, 2010). Chloride is likely to be relatively unreactive in the Erlenbach catchment.

Nitrate concentrations in Erlenbach streamwater can also be mainly attributed to atmospheric inputs as well as manure inputs from grazing. The output–input flux index of nitrate was −0.52 (even without taking manure inputs into account), indicating that nitrate was either taken up into vegetation or possibly volatilized as ammonia (i.e., the catchment acted as solute sink). Potassium is supplied to the stream through weathering of the bedrock, with some atmospheric inputs, which contribute more to streamwater fluxes of potassium than those of calcium or magnesium. Potassium is also known to be retained in the soils and cycled internally in forest stands (Hornung et al., 1986; Likens et al., 1994).

Iron was the most abundant trace metal in the Erlenbach streamwater, with median concentrations of around 5 µg L⁻¹. Streamwater iron concentrations were substantially higher than the concentrations measured in groundwater or precipitation, and very high in the soil water of the neighboring and geologically similar Studibach catchment, pointing toward the soil layer as a predominant source of streamwater iron (Table 2). Groundwater concentrations of manganese at Erlenbach were broadly similar to those measured in streamwater, but soil water concentrations (and also groundwater concentrations; not shown here) in the Studibach catchment were substantially higher. This indicates that manganese may be sourced both from groundwater and soil water. The storage of iron and manganese in and release from soil layers are likely controlled by the complexation of these trace metals with organic material (Bloomfield, 1953; Harter and Naidu, 1995), and the strong redox sensitivity of these elements (Drever, 1988; Basu et al., 2010; Koger et al., 2018). The concentrations of the trace metals chromium and copper were low in all analyzed compartments, and a distinct source cannot decisively be identified. In the neighboring Studibach catchment, Kiewiet et al. (2019) observed the highest concentrations of chromium at predominantly dry sites and highest concentrations of iron and manganese at predominantly wet sites. These observations suggest that concentrations of trace metals are probably not homogeneously distributed in the Erlenbach catchment.

3.3 High inter-event variability that differs from the 2-year behavior

We calculated 2-year cQ relationships based on the snow-free recession periods of the 2-year time series and event-scale cQ behavior from the recessions of the 30 extracted events. The 2-year cQ relationships were relatively narrow and well-defined for most solutes (Fig. 2), resulting in low uncertainties of 2-year slopes and intercepts (Table S1 in the Supplement). Furthermore, 2-year slopes and intercepts calculated for the whole time series were broadly similar to those determined from only hydrograph recession periods (Table S1), largely because recessions comprised nearly all of the 2-year data. Thus, our analysis can straightforwardly compare cQ behavior during individual recessions with the 2-year cQ behavior across 2 years of recessions.

Figure 2The 2-year cQ relationships from the entire recession time series (data points are shown in light blue and the dark blue line indicates the 2-year power-law fit), compared to the cQ behavior of individual events (colored lines; up to 20 events are shown for better visibility). Discharge is normalized by the average discharge of the time series. Most solutes exhibit a 2-year dilution pattern, whereas the trace metals iron, manganese, and chromium exhibit a 2-year mobilization pattern. The cQ relationships vary much more from event to event for some solutes (e.g., potassium, chloride, and nitrate) than for others (e.g., calcium, magnesium, sodium, and EC).

Download

Figure 2 shows that major ions and EC exhibited dilution behavior (slopes <0; lower concentrations at higher discharges) across the 2-year record, but with different degrees of variability. The data clouds in the cQ space were more scattered (e.g., chloride, potassium, and nitrate) compared to those of calcium, magnesium, and sodium. Iron, manganese, and chromium showed a mobilization behavior (2-year slope >0; higher concentrations at higher discharges), while the 2-year cQ relationship of copper indicated chemostasis (slopes ≈0).

Event-scale cQ relationships at Erlenbach were much more variable than the 2-year behavior. A comparison between different events revealed substantial inter-event variability in both slopes and intercepts for a number of solutes (colored lines in Fig. 2 and blue circles in Fig. 3). The extent of this variability differed dramatically among the solutes. Individual events followed the 2-year trend for calcium, magnesium, sodium, strontium, and EC, for example, whereas they deviated substantially from the 2-year trend for potassium, chloride, and nitrate.

Figure 3A comparison of event-scale and 2-year cQ intercepts (a) and slopes (b). Blue circles represent intercepts and slopes of individual events, light-blue diamonds represent the averages of all events, and red bars indicate the slopes and intercepts of the 2-year cQ relationships. If the red bar is close to the light-blue diamond, the 2-year slope or intercept value is a good approximation of the average event slope or intercept. In panel (a), event cQ intercepts are expressed relative to the 2-year average for better visual comparison of the solutes. Event-scale intercepts and slopes vary substantially for solutes in the lower half of the figure (from chloride to chromium), but vary little for most groundwater-sourced solutes. Because we normalized discharge by the mean discharge value, 2-year intercepts approximate the average event intercepts. The 2-year slopes also approximate the average event slopes for most groundwater-sourced solutes, but 2-year slopes are more negative than event slopes for other solutes, such as chloride, nitrate, and manganese.

Download

For most solutes, the averages of all event-scale intercepts (light-blue diamonds in Fig. 3a) were similar to the 2-year intercepts (red bars in Fig. 3a), because discharge was normalized by its mean value (see Sect. 2.5). The inter-event variability in intercept values was particularly high for chloride, nitrate, and some of the trace metals, possibly suggesting changes in sources, flow paths or reaction rates of these constituents between events.

Conversely, the slopes of the event-scale cQ relationships differed substantially from the 2-year behavior for some of the elements (Figs. 2 and 3). This was true in particular for three solutes with significant atmospheric inputs (chloride, potassium, and nitrate), which exhibited dilution behavior in the 2-year dataset (2-year slope <0), while their cQ slopes of individual events varied between dilution (event slope <0) and mobilization (event slope >0).

None of the solutes exhibited stronger dilution behavior, on average, during events (blue diamonds in Fig. 3b) than over the 2-year period (red bars in Fig. 3b). In the case of potassium, chloride, nitrate, and manganese, most individual events (blue circles in Fig. 3b) exhibited substantially weaker dilution, or greater mobilization, than the 2-year cQ relationships (red bars in Fig. 3b). Considered together, these results indicate that dilution processes were equally important, or more important, over the long term than during individual events.

Solutes with similar sources and similar chemical properties tend to exhibit similar 2-year cQ behavior and similar degrees of inter-event variability in slopes and intercepts. For example, calcium, magnesium, sodium, strontium, barium, boron, sulfate, and EC all exhibited well-defined 2-year dilution behavior and little variability in both cQ slopes and intercepts at the event scale. All solutes exhibiting this behavior were previously identified as groundwater-sourced in the Erlenbach catchment (Sect. 3.2). Average event slopes of these solutes were approximated well by the 2-year behavior (Fig. 3b), resulting in event patterns that fan out around the high-discharge section of the 2-year cQ trend (Fig. 2). This similarity between 2-year and event-scale behavior indicates that similar mechanisms controlled solute and water mobilization both during individual events and over the long term for these solutes.

Conversely, potassium, chloride, and nitrate have substantial atmospheric inputs at Erlenbach. These three solutes were characterized by 2-year dilution behavior, and high variability in event slopes ranging from negative to positive values, along with some inter-event variability in intercept values. These solutes exhibit event cQ patterns that fan out and stack up around the 2-year relationship. This large inter-event variability in cQ slopes may plausibly arise from varying degrees of atmospheric deposition, evapotranspiration, and dry deposition, resulting in temporally and spatially variable concentrations of chloride, potassium, and nitrate in shallow soil layers. A large variability in soil water concentrations was indeed observed in the neighboring Studibach catchment (Table 2). In the case of nitrate and potassium, the inter-event variability may also be affected by biological uptake and reaction processes. In summary, chloride, nitrate, and potassium showed a much more diverse behavior on the event scale and also – at least to some extent – on longer timescales than the purely groundwater-sourced solutes.

The trace metals iron, manganese, and chromium were mobilized during events, as reflected in positive event and 2-year cQ slopes. Interestingly, this mobilization behavior was more pronounced (i.e., slopes were steeper) on the event scale than on longer timescales. In contrast to the other metals, copper concentrations in streamwater suggested chemostatic behavior over long timescales, whereas some mobilization behavior was evident on the event scale (Fig. 2). These findings suggest that the trace metals were mainly sourced from the soil layers, and the observed cQ behavior of manganese and iron confirms common assumptions about them being co-located in the soil layers where they are bound to organic material or present as oxides.

The grouping of solutes based on cQ behavior aligns well with the previously identified source areas in the Erlenbach catchment (Sect. 3.2). These alignments were most obvious for the purely groundwater-sourced solutes because they could be clearly attributed to one main source. For most other solutes, the combination of different sources and the co-occurrence of reaction and mixing processes resulted in less clear patterns. A scatterplot of event slopes against event intercept values (Fig. 4) supports the grouping of solutes and illustrates the different degrees of inter-event variability. The low uncertainty of the calculated slopes and intercepts furthermore confirms that the observed variability among events is real-world behavior rather than noise.

Figure 4Scatterplots of event cQ slopes and intercepts of the 14 different solutes and EC (error bars indicate one standard error). Solutes from different dominant sources cluster and exhibit similar ranges of variability. (a) Groundwater-sourced solutes cluster closely around similar slopes and show little inter-event variability in both slopes and intercepts. (b) Intercepts and slopes of solutes with significant atmospheric input (i.e., chloride, nitrate, and potassium) vary substantially among events. (c) The slopes of trace metals are generally higher than those of the other solutes (indicating predominantly mobilization behavior) and are also highly variable among events. The uncertainty in the estimated slopes and intercepts is mostly smaller than the variability among events, indicating that the observed inter-event variability in slopes and intercepts reflects real-world behavior rather than sampling and measurement noise. To corroborate this, a version of the same figure including all cQ slopes and intercepts, also those from events excluded from further analysis due to high relative standard errors in slope and/or intercept values, is provided in the Supplement.

Download

Other likely modulators of the variability observed for different solutes are their ionic form and their possible occurrence as nanoparticulates. Cations are known to undergo exchange buffering through electrostatic binding to negatively charged sites in soils (Helling et al., 1964; Rhoades, 1982). This likely resulted in less variable cQ behavior on the event scale for the cationic solutes calcium, magnesium, sodium, strontium, and barium compared to the anionic solutes sulfate and boron. These anions are less buffered by ion exchange and showed a far more variable cQ behavior. Two other anions, chloride and nitrate, also exhibited highly variable cQ behavior at the event scale. It consequently seems likely that exchange buffering processes modulate the degree of observed inter-event variability. Nevertheless, the importance of cation exchange buffering on the variability of solute behavior cannot be well constrained with our dataset because most analyzed cations are groundwater-sourced species, while two of the three solutes with substantial atmospheric inputs, namely chloride and nitrate, are anions. The only cation with important atmospheric input (potassium) is biologically very active, potentially obscuring the effects of cation exchange buffering. With respect to the trace metals, a large part of their total concentration can occur as natural nanoparticulates (Gottselig et al., 2017), which may undergo different transport and adsorption processes compared to their dissolved forms. This likely had a strong effect on their cQ behavior at Erlenbach.

3.4 Environmental controls of the observed inter-event variability in concentration–discharge relationships

We used weighted rank correlation coefficients to assess how variations in cQ slopes and intercepts from event to event were related to seasonality indicators, relative input concentrations, antecedent wetness conditions, event characteristics, and event-water contributions. A heatmap (Fig. 5) illustrates how event-scale cQ slopes and intercepts depended on these different environmental controls. Individual examples of these relationships are shown in Fig. S2 in the Supplement for different solutes and drivers.

Figure 5Weighted rank correlation coefficients expressing the dependence of event-scale cQ intercepts (a) and slopes (b) on different environmental controls (seasonality indicators, relative input concentrations, antecedent wetness conditions, event characteristics, and importance of event-water contributions). Green fields indicate positive rank correlations, blue fields indicate negative correlations, and darker colors indicate stronger relationships. Only statistically significant (p<0.05) correlation coefficients are displayed. Gray fields for EC indicate relationships that could not be assessed because EC was not measured in precipitation. For solutes featuring dilution patterns, e.g., calcium in panel (b), positive correlations indicate relationships that become more chemostatic when the controlling variable increases. For solutes with mobilization patterns, e.g., iron in panel (b), positive correlations indicate enhanced mobilization when the controlling variable increases. For meanings of abbreviations, please refer to Table 1.

Download

The event-scale cQ intercepts for all solutes, regardless of their dominant source, generally responded similarly to each control (with the exception of chloride; Fig. 5a). Seasonality indicators and antecedent wetness conditions were the most important environmental controls on cQ intercepts for all solutes, with drier and warmer conditions resulting in higher concentrations in streamwater. By contrast, although the effects of each control on event-scale cQ slopes were broadly similar within groups of solutes that shared the same dominant source (precipitation, soil water, and groundwater), they often differed between these groups (Fig. 5b). In the following, we therefore discuss the behavior of the event-scale cQ slopes of these groups of solutes together, rather than discussing individual solutes.

We found that the event-scale cQ slopes of groundwater-sourced solutes were positively correlated with seasonality indicators, with weaker solute dilution (i.e., less negative cQ slopes) during warmer conditions. Antecedent wetness conditions and event characteristics were other important factors controlling the event-scale cQ behavior of most groundwater-sourced solutes. We found stronger dilution (i.e., more negative cQ slopes) during larger events associated with wetter antecedent conditions. cQ slopes were also more negative when event-water contributions were larger, consistent with stronger dilution of groundwater-sourced solutes by larger volumes of recent rainfall (i.e., event water). Most of the observed relationships were less clear for boron (which is primarily sourced from groundwater, but does not occur as a cation, but instead as either undissociated boric acid or as the borate anion, depending on pH), supporting the hypothesis that differences in ionic charge modulate event-scale cQ behavior. By contrast, however, sulfate (an anion) showed very similar dependencies as the cations in groundwater, contradicting this hypothesis. Streamwater electrical conductivity (EC) at Erlenbach is dominated by the two most abundant solutes, calcium and magnesium (along with their counter-ion, bicarbonate), which are groundwater-sourced. Consequently, the event-scale cQ behavior of EC depended on similar factors as that of the groundwater-sourced solutes.

For the solutes with significant atmospheric inputs (chloride, potassium, and nitrate), the effects of possible environmental controls on cQ slopes were less clear, but still indicated similar behavior among the solutes of this group. We found that the event-scale cQ slopes of chloride and nitrate were influenced by antecedent wetness conditions, with dilution behavior during wetter conditions and mobilization behavior during drier conditions (Fig. 5b). This can potentially be explained by stronger evapoconcentration of these solutes in the soil under drier conditions. We also observed a tendency toward more positive slopes with higher relative input concentrations for all three solutes, and during events in which the event-water contribution was larger. The observed dependencies of nitrate and potassium differed somewhat from those of chloride, likely due to the overprinting with reaction processes acting on these solutes, and due to overlapping patterns from atmospheric and groundwater sources in the case of potassium.

Among the trace metals, iron often exhibited contrasting behavior compared to the others (Fig. 5b). While manganese, chromium, and copper showed a tendency toward stronger mobilization under colder and wetter conditions, iron tended to be mobilized more during warmer and drier conditions. Given that both iron and manganese are found predominantly in wetter soils (in the neighboring Studibach catchment; Kiewiet et al., 2019), we would have expected these two metals to depend on similar drivers. However, Kiewiet et al. (2019) also observed higher concentrations of manganese in riparian areas compared to iron. We would expect these riparian zones to be activated first during a rain event, whereas soils farther from the stream would likely only start to contribute to streamflow later. This sequence of contributions from source areas dominated by different trace metals could potentially explain why the cQ slopes of the trace metals varied differently with seasonality, event characteristics, and event-water contributions. The high concentrations of manganese in riparian areas may also explain the very high cQ slopes observed for this element. As the relative contribution from riparian areas decreases later in the recession, manganese concentrations in streamwater drop rapidly, decreasing more than proportionally compared to discharge. As metal complexation and mobilization is known to depend on various factors such as pH and redox conditions in the soil layer (Gotoh and Patrick, 1972, 1974), further field measurements are necessary to better understand the mobilization of trace metals from soil layers during hydrologic events.