Understanding The Mechanisms Controlling The Concentration-Discha
Understanding The Mechanisms Controlling The Concentration-Discha
2020
Recommended Citation
Porter, Veronica, "Understanding the Mechanisms Controlling the Concentration-Discharge Relationships
of Streams in Northeastern and North Midwestern United States", Open Access Master's Thesis, Michigan
Technological University, 2020.
https://doi.org/10.37099/mtu.dc.etdr/1098
By
Veronica Porter
A THESIS
Submitted in partial fulfillment of the requirements for the degree of
MASTER OF SCIENCE
In Forest Ecology and Management
Acknowledgements....................................................................................... vii
iii
2.6.3 Seasonal Pattern ................................................................................26
2.7 CONCLUSION ..............................................................................................27
2.8 REFERENCES ...............................................................................................28
2.9 TABLES AND FIGURES..............................................................................34
2.9.1 Tables ................................................................................................34
2.9.2 Figures...............................................................................................37
38
2.10 APPENDIX 2.A SUPPLEMENTAL TABLES .............................................47
2.11 APPENDIX 2.B SUPPLEMENTAL FIGURES ............................................62
iv
3.9 TABLES AND FIGURES..............................................................................90
3.9.1 Tables ................................................................................................90
3.9.2 Figures...............................................................................................91
3.10 APPENDIX 3.A: SUPPLEMENTAL TABLES ..........................................101
4 Conclusion ...........................................................................................105
4.1 REFERENCES .............................................................................................107
v
Preface
The chapters of this thesis are composed of two multi-author, journal articles that
are in preparation to be submitted for publication in the scientific journals: Hydrological
Processes (Chapter 2) and Journal of Great Lakes Research (Chapter 3). For these articles,
I primarily led the research (study design, data analysis, and interpretation) and wrote the
manuscripts. The data used for Chapter 2 was collected by numerous researchers and
technicians at Sleepers River Research Watershed and the data in Chapter 3 was collected
by USGS as part of their water watch program in addition to multiple public datasets. Both
articles were written in collaboration with Fengjing Liu who also acquired the funding,
guided the studies, assisted with the end-member mixing analysis, and edited the
manuscripts. Stephen Sebestyen and Jamie Shanley aided in data analysis questions and
draft revisions of Chapter 2.
vi
Acknowledgements
I want to thank Fengjing Liu for the intellectual, analytical and writing advice and
securing funding for the project. I want to thank the funding sources that made this research
possible: the McIntire-Stennis Cooperative Forestry Research Program and Michigan
Technological University’s College of Forest Resources and Environmental Science,
including Ecosystem Science Center and the Graduate Student Government. Thanks to
Jamie Shanley, and Stephen Sebestyen for all their guidance, analytical and technical
support, and willingness to collaborate by sharing their valuable datasets. I would like to
thank my committee members Stephen Sebestyen and Veronica Webster for their
contributions and time. Thanks to Joseph Shannon who patiently answered questions and
selflessly provided coding help in R. Lastly, thanks to my family and friends who have
been supportive and encouraging throughout these last two years.
vii
Abstract
Concentration-discharge (C-Q) relationships, concentration expressed as a power-
law function of stream discharge (C = aQb), are an easy and efficient way to gain insight
into solute export processes within a catchment. Despite decades of research, we have
limited understanding of the controls on C-Q relationships, consequently making it difficult
to confidently apply findings to unstudied catchments. Additionally, there is limited
knowledge on how changes in climate and land uses may affect solute export by modifying
solute availability and stream water chemistry at catchment scales.
The goal in Chapter 2 was to close some of those knowledge gaps by understanding
factors controlling variation in stream water chemistry between timescales through end-
member mixing analysis and concentration-discharge analysis. Over a decadal timescale,
major nutrients and some trace elements exhibited a positive relationship with discharge
(positive b for flushing), while major ions, particularly geogenic ions, displayed a
contrasting pattern of dilution (negative b). I found a key factor controlling export,
represented by the slope of the C-Q power-law relation, is the vertical distribution of
solutes between shallow to deeper flowpaths. A negative b-coefficient corresponds to an
increase in solute concentrations with depth, while a positive b is caused by a decrease.
When solute concentrations varied greatly between deep and shallow flowpaths, the
magnitude of the b-coefficient was higher than when there was little variation in
concentrations. Interannual variability was much lower for conservative solutes (resulting
from mixing rather than chemical equilibrium over a short time period, e.g., days) than for
non-conservative solutes, making the export pattern of conservative solutes, mostly
geogenic ones, more predictable and useful in forensic hydrology. During snowmelt, the
magnitudes of b-coefficients of most solutes were higher on the falling than rising limbs
due to greater contributions of saturation overland flow. These findings indicate that
seasonal stream chemistry will likely shift under new climate conditions as precipitation
patterns change with a decrease in snow abundance and earlier snowmelt.
viii
The objective in Chapter 3 was to examine patterns in C-Q relationships throughout
138 catchments within the Great Lakes Basin (GLB) and the catchment physical and
environmental factors (e.g., land-use and precipitation pattern) controlling the C-Q
behaviors. Major geogenic ions were consistently diluted throughout the GLB during
periods of high flow, while major nutrients and trace elements were enriched. Snow
abundance and land-use were frequently correlated to the R2 values of C-Q relationships
and the pattern and intensity of changes in solute concentrations. In more forested areas,
the strength of C-Q relations (R2) for specific conductance (SC) increased, while that of
sulfate and ammonium decreased. There was also increased dilution of SC, calcium, and
magnesium and increased enrichment of phosphorus in catchments with more forested
area. The results indicate that under future conditions of increasing urban areas and
decreasing snow abundance, the C-Q power-law relationships may become weakened for
major geogenic solutes, suggesting that solute export for multiple solutes will become more
chemostatic and the buffering capacity for nutrients such as NO3- and contaminants will be
reduced.
ix
1 Chapter 1: Concentration-Discharge Relationships
1.1 INTRODUCTION
Analyzing the relationship between solute concentrations and changes in stream
discharge is a valuable analytical method used to study hydrological processes and solute
export processes at a catchment scale. Many studies have shown that concentration-
discharge (C-Q) analysis is a useful tool in understanding the internal processes controlling
streamflow generation, solute distribution and availability, export flowpaths, and timing of
export (Duncan, Welty, Kemper, Groffman, & Band 2017; Hofmeister et al., 2019). Shifts
in C-Q relations across spatial-temporal scales and varying landscapes reflect dominant
flowpaths and solute availability controlling stream water chemistry creating an effective
way to understand integrated biological and hydrological processes within catchments
(Hofmeister et al., 2019; Zimmer, Pellerin, Burns, & Petrochenkov 2019). Concentration-
discharge relationships can be one tool to improve the understanding of hydrologic
mechanisms and predict potential impacts of climate change on the quality and quantity of
future water resources. (Hunsaker & Johnson, 2017; Litt, Gardner, Ogden, & Lyons 2015)
One benefit of the C-Q analytical method is that it does not require high-frequency
sampling or data collection from multiple water sources, making it a relatively inexpensive
but effective method to gain an understanding of hydrological and biogeochemical
processes occurring within a catchment. To calculate C-Q relationships, data collection
involves simultaneously measuring stream discharge and the concentration of a solute
either by in situ measurements or grab samples typically at the outlet of the catchment. The
concentration-discharge relationship is represented by a power-law function such that:
𝐶𝐶 = aQb (1)
1
where ℇ is the error (Hofmeister et al., 2019; Zimmer et al., 2019). The strength of the C-
Q relation is measured by the coefficient of determination (R²) which indicates how much
of variation in solute concentration is explained by variation in stream discharge. The
significance of the relationship is measured by the p-value of the regression with a value
of less than 0.05, 0.01 or 0.001.
When using the logarithmic scale, the slope (b) indicates the trend of solute
concentration with flat slopes (-0.1< b ≤ 0.1) indicating chemostatic behavior, positive
slopes (b > 0) indicating a flushing response, and negative slopes (-0.1 ≥ b) occurring with
dilution (Evans & Davies, 1998; Herndon, Dere, Sullivan, Norris, Reynolds, & Brantley
2015; Hofmeister et al., 2019; Zimmer et al., 2019). Solutes with positive slopes are
theorized as transport-limited, those with negative slopes are source limited, and those with
no slope are dominated by neither process (Hofmeister et al. 2019, Godsey et al. 2009,
Zimmer et al. 2019).
The slope may also change between the rising and falling limbs of snowmelt or
event hydrographs, creating a hysteresis pattern. Clockwise hysteresis occurs when the
rising limb has higher solute concentrations than the falling during the same values of
stream discharge, often related to mineral weathering in the vadose or saturated zone
(Evans & Davies, 1998). An anti-clockwise hysteresis is the reverse where the rising limb
has lower solute concentrations than the falling (Evans & Davies, 1998). Evans & Davies
(1998) argued that one can infer dominant mixing components from the direction and shape
of the hysteresis. This should be done cautiously since changes in stream chemistry does
not always reflect the timing of inflows (Chanet & Hornberger, 2003). Hysteretic flushing
patterns are often associated with transient (short-term) water sources and flowpaths that
intersect surface soils with higher organic carbon content and higher hydraulic conductivity
than deeper soils (Bishop, Seibert, Nyberg, & Rodhe, 2011; Duncan et al., 2017; Kendall,
Shanley, & McDonald, 1999; McIntosh et al., 2017). Additionally, hysteretic patterns can
reveal information on the moisture content of soils since well-drained soils are less likely
to exhibit a dilution effect compared to less well-drained soils (McIntosh et al., 2017)
2
Thus, C-Q analysis provides a diagnostic tool to quickly and efficiently assess water
quality changes in response to stream flow changes by reflecting the internal catchment
structures (e.g., hydrologic pathways) influencing hydrological, biogeochemical, and
ecological functions (Liu, Conklin, & Shaw, 2017). However, research investigating C-Q
relationships is still young. The general patterns of C-Q relationships are readily observed,
yet, the underlying causes of variability over timescales and drivers (e.g., water and energy
inputs) are still debated (Chanet, Rice, & Hornberger, 2002; McIntosh et al., 2017). This
knowledge gap makes extrapolating findings from C-Q analysis to unstudied catchments
difficult. Teasing out the factors contributing to this unknown variability is made more
difficult with climate change and changes in land use, which can exacerbate the natural
complexity of these systems. However, understanding these interactions between natural
processes and human activities is necessary in balancing environmental conservation and
human development (Hunsaker & Johnson, 2017). Thus, the goal of this thesis is to bridge
this knowledge gap by furthering our understanding of factors controlling C-Q
relationships by focusing on the influences of climate, the critical zone, and land-use with
the aim to increase the applicability and interpretation of C-Q relationships as a tool to
solve water management issues.
1.2 REFERENCES
Bishop K H, Seibert J, Nyberg L, & Rodhe A. (2011), Water storage in a till catchment.
II: Implications of transmissivity feedback for flow paths and turnover times,
Hydrological Processes, 25(25), pp. 3950-3959. doi:10.1002/hyp.8355
3
Evans C, & Davies T D. (1998). Causes of concentration/discharge hysteresis and its
potential as a tool for analysis of episode hydrochemistry. Water Resources
Research, 34(1), pp. 129–137. doi:10.1029/97WR01881
Liu F, Conklin M H, and Shaw G D. (2017). Insights into hydrologic and hydrochemical
processes based on concentration-discharge and end-member mixing analyses in
the mid-Merced River Basin, Sierra Nevada, California, Water Resources
Research. 53, pp. 832– 850. doi:10.1002/2016WR019437
4
Zimmer M A, Pellerin B, Burns D A, & Petrochenkov G. (2019). Temporal variability in
nitrate-discharge relationships in large rivers as revealed by high-frequency data.
Water Resources Research, 55, pp. 973 – 989. doi:10.1029/2018WR023478
5
2 Chapter 2: Controls on Long-Term, Annual, and Seasonal
Solute Export Patterns in a Forested Headwater Watershed
in Northeastern U.S.
2.1 Abstract
The objective of this study was to compare responses in stream water chemistry to
changes in discharge over decadal, annual, and seasonal timescales to better understand
solute export patterns using concentration-discharge (C-Q) relationships and end-member-
mixing analysis (EMMA). C-Q power-law relationships (C=aQb) were assessed using 15
solutes for water years 1992-2015 at W-9 of Sleepers River Research Watershed (0.405
km2), Vermont. Over a decadal timescale, NO3-, DOC, and Fe exhibited a positive
relationship with discharge (flushing), while SC, Ca2+, Mg2+, Na+, K+, SO42-, Si, Sr, Ba,
and Mn displayed a negative relationship (dilution). Cl- and Al were chemostatic. The
power-law relationship is primarily controlled by vertical distribution of solutes from
shallow to deeper flowpaths. A negative b-coefficient corresponds to an increase with
depth while a positive b-coefficient corresponds to a decrease. The b-coefficient decreased
significantly (p < 0.05) with an increase in concentration ratio between deeper and shallow
end-members for geogenic solutes. Geogenic solutes that are conservative, distinct, and
temporally relatively constant in the three end-members identified by EMMA, showed
consistent annual C-Q patterns. However, C-Q patterns for the remaining solutes did
change between years, more frequently switching from chemodynamic to chemostatic
behavior than remaining chemodynamic and switching the b-coefficient sign. During
snowmelt with greater contributions of saturation overland flow, the power-law relations
became stronger and the b-coefficients of geogenic solutes were frequently more negative
on the falling than rising limbs. As snowmelt timing and precipitation intensities are likely
to change in the future, concentration-discharge relationships in streamflow could change,
particularly at the seasonal scale.
6
2.2 INTRODUCTION
For decades, hydrologists have studied processes governing solute export through
concentration-discharge (C-Q) analysis to better understand export behavior in relation to
catchment characteristics, precipitation patterns, biotic, and abiotic processes (Bieroza,
Heathwaite, Bechmann, Kyllmar, & Jordan, 2018; Evans & Davies, 1998; Godsey et al.,
2009; Hall 1970; Hall 1971; Johnson, Likens, Bormann, Fisher, & Pierce, 1969; Moater,
Abbott, Minaudo, Curie, & Pinay, 2017; Pinder & Jones, 1969). Limited by chemical data,
few studies have focused on the long-term patterns of C-Q relationships. Thus, we have
minimal understanding of how the C-Q relationship changes in response changes in
precipitation, namely increases in rainfall intensity and declining snow abundance.
Understanding export process dynamics is particularly important since these processes will
likely change under predicted climate conditions. These climatic changes will modify
many of the dominant factors controlling solute export such as snow cover extent,
snowmelt timing, winter temperatures, duration of seasons, and stream discharge (Hidalgo
et al., 2009; Derksen & Brown, 2012; Pierce & Cayan, 2013; Demaria, Roundy, Wi, &
Palmer, 2016). These changes will continue to occur as concentrations of greenhouse gases
rise and further modify water budget partitioning, solute export, and nutrient cycling
processes, which will have important implications for water quality and watershed
management (Markstrom & Hay, 2009; Maher, Shanley, Bailey, & Mitchell, 2010; Krause
et al., 2017).
7
et al., 2014; Winnick, Carroll, & Williams, 2017). By understanding current linkages
between solute export, activation of hydrologic flowpaths, and end-member contributions,
we can better identify the implications of these changes on stream water chemistry and
better address future water quality vulnerabilities.
8
McDonnell, 1999; Ohte et al. 2004; Sebestyen et al., 2008, Shanley et al., 2015; Smith
1997). Kendall et al. (1999) showed that 94% of the variance in stream chemistry could be
explained by a two-component mixture of groundwater high in base cations and an O-
horizon/overland flow water high in DOC. Smith (1997) also used EMMA to find three
end-members, including precipitation (rain and meltwater), soil water, and groundwater.
Both studies selected the chemical species as conservative tracers based on assumptions
and their behavior in the catchment (Kendall et al, 1999; Smith 1997). Two of these studies
occurred before diagnostic tools of mixing models (DTMM) were developed, which are
tools to select conservative tracers and determine the number of end-members (Hooper,
2003). Sebestyen et al. (2008) used DTMM but performed EMMA only over the 2004
snowmelt period. A recent study has demonstrated that combining end-member mixing
analysis and concentration-discharge relationships strongly enhances watershed
hydrologic and biogeochemical studies and provides physically meaningful results of end-
member mixing analysis (Liu, Conklin, & Shaw, 2017).
9
2.3 SITE DESCRIPTION
Sleepers River Research Watershed (44° 28'N, 72° 9'W), a long-term research
watershed located in rural north-eastern Vermont, was established in 1957 by the
Agricultural Research Service (ARS) of the United States Department of Agriculture
(USDA). Watershed 9 (W-9) is a 0.405 km2 forested sub-basin included in Sleepers River
Research Watershed. Elevation ranges from 524 to 679 m and is dominated by south-facing
slopes ranging from 0 to 90% with a mean of 13 (Shanley, Sundquist, & Kendall, 1995).
Temperature ranges from an average of -10°C in winter to 12°C in summer with an annual
mean of 4.6°C (Shanley, Kram, Hruska, & Bullen, 2004). The area is impacted by acid
rain, but the high buffering capacity from carbonates in bedrock and till limits the negative
impacts (Hornbeck, Bailey, Buso, & Shanley, 1997) There are also chronically high
nitrogen inputs typical in northeastern United States with the average total N input from
1978 to 1998 of 13.2 kg ha-1a-1 (Campbell et al. 2004). Mean annual precipitation of 1100
mm falls evenly over seasons throughout the year with 25-35% of precipitation falling as
snow, and snow typically accumulating from December until late March or April (Shanley
et al., 2015; Armfield et al., 2019).
Soils belong to Spodosol, Inceptisol, and Histosol soil orders that range from 50 to
70 cm depth on top of up to three meters of dense till and silt (Shanley et al., 2008; Mayer
et al., 2010). Soil hydraulic conductivities exponentially decrease with depth while water
levels in upland areas vary over several meters depending on antecedent wetness conditions
(Kendall et al., 1999; Shanley et al., 2003). The forest overstory at W-9 is dominated by
sugar maple, yellow birch, and white ash with an understory of striped maple, sugar maple,
10
hobblebush, balsam fir, and red spruce (Sebestyen et al., 2008). The forest was last
harvested in 1960 when there was a selective harvest, and the forest was partially logged
in 1929 (Sebestyen et al., 2008).
2.4 METHODOLOGY
Soil water and groundwater samples were collected from numerous zero-tension
lysimeters, piezometers, and wells throughout the catchment. The soil water was collected
from three sets of 2-3 nested shallow soil zero-tension lysimeters (names starting with LYS
in Table 2.A.1) along a hillslope transect (Figure 2.1). The shallowest lysimeters (SH) are
just under the O-horizon (0.05-0.1 m), the middle depth (MID) is at the top of the B-horizon
(~0.3 m), and the deepest (DP) is at the bottom of the B-horizon (0.7-1 m). Samples from
the lysimeters were collected after a rainfall or melt event. The inlets of the lysimeters at
the middle and deep depths were frequently under the top of the water table. Six
piezometers (T1-T6) were deployed at two locations in the riparian zone along the stream.
Groundwater samples were also collected from 36 groundwater wells (BW-##) at varying
depths and sampling frequency at various landscape positions, including riparian zone,
hillslope, hillslope bench, and hillslope hollow (Figure 2.1; Table 2.A.1). There were
between 1 to 30 samples collected at 32 of the wells for each solute. There were only 4
wells where samples were collected more frequently with 30 to 166 samples collected for
each solute (Table 2.A.1). Precipitation was measured using a weighing bucket gage at the
11
meteorological station (R29) in a clearing adjacent to the stream gauge. Precipitation
chemistry samples were collected from a polyethylene bucket in a wet-only precipitation
collector near the weighing bucket after each event or weekly.
To look at seasonal variability in solute export mechanisms, each water year was
divided into four periods: the dormant season (Oct. 1 – beginning of snowmelt), the rising
limb of snowmelt (the beginning of snowmelt to date of peak streamflow during
snowmelt), the falling limb of snowmelt (peak snowmelt to the return to baseflow
conditions), and the growing season (return to base flow to Sept. 30). For each water year,
the snowmelt period was visually identified as the date ranges during the rising and falling
limbs using stream discharge, precipitation, and temperature data. The rising limb of the
snowmelt period started when there was a steady, persistent increase in stream discharge
that was coupled with consecutive mean daily temperatures greater than 0°C. The falling
limb began after peak streamflow and continued until streamflow returned to near baseflow
conditions.
12
Years that had less than 20 data points were discarded for the annual C-Q analysis
because there were too few samples, but these points were included in the multi-annual C-
Q relationship analysis. Similarly, seasons were excluded in the seasonal C-Q analysis if
they had less than 5 data points. The rising and falling limbs of snowmelt frequently had
the fewest data points, ranging from 5 to 52 with a mean of 16 and a standard deviation of
13.
13
indicates greater variability in the solute concentrations and that discharge is not the sole
factor influencing concentrations.
14
2.5 RESULTS
Solute concentrations in stream water varied slightly over seasons and events, as
demonstrated by Ca2+, Cl-, NO3-, and DOC in Figures 2.B.1c-1e. The SC and median
15
concentrations of all solutes except NO3- were highest during the growing and dormant
period. Most solutes had lowest median concentrations during the snowmelt period, except
DOC, which was lowest during the dormant period, and NO3-, which was lowest during
the growing and dormant season (Table 2.A.3). Unlike other solutes, median DOC, NO3-,
and Fe concentrations in stream water were higher during high flow conditions (e.g., NO3-
= 0.9 mg L-1, DOC = 3.6 mg L-1, and Fe = 17.1 μg L-1 at the 95th flow percentile) compared
to low flow conditions (e.g., NO3- = 0.5 mg L-1, DOC = 1.3 mg L-1, Fe = 8.7 μg L-1 at the
5th flow percentile) (Table 2.A.4).
The solute concentrations in groundwater often varied with the geographic location
of the wells (Table 2.A.6). Riparian wells had the highest median concentrations of SC,
16
Ca2+, Mg2+, K+, SO42-, Sr, and Ba, while the hillslope hollow wells had highest median
concentrations of Cl- and Al (Table 2.A.1). The hillslope wells had on average higher
median concentrations of Na+ (1.1 mg L-1) and Si (11.0 mg L-1) while the hillslope bench
wells had the highest median concentrations of NO3-.
17
The annual C-Q relationships were consistent with the multi-annual ones for most
solutes (Figures 2.3 and 2.5). However, SC, K+, Cl-, NO3-, SO42-, DOC, Si, Al had at least
one year where the C-Q pattern changed. These eight solutes shifted between a
chemodynamic behavior to a chemostatic behavior. Al was the only solute that had
significant annual C-Q relationships for all three C-Q patterns. CVC/CVQ mean values for
all the solutes ranged from 0.14 to 1.51 with standard deviations ranging from 0.02 to 0.55
(Table 2.2). Most CVC/CVQ values were less than 0.50, while DOC, Al, and Fe were
between 0.50 and 1.0 (Table 2.2).
Some solutes had individual seasonal export patterns that differed from the multi-
annual pattern. The median b-coefficient for Si remained chemostatic for all periods except
the growing season where it was diluted. SC, K+, NO3-, Ba, and Mn also had periods during
which the median b-coefficient showed evidence of chemostasis (Figure 2.6). Conversely,
during the seasons where Al had significant C-Q relationships, the medians of all four
periods showed evidence of flushing. Aside from those six solutes, all others had seasonal
median b-coefficients consistent with their respective multi-annual C-Q relationship.
When concentration was plotted against discharge during the snowmelt period, few
solutes showed consistent visible hysteresis patterns (Figure 2.B.2). NO3- and DOC were
the only solutes that seemed to show weak hysteretic behavior with NO3- exhibiting a
clockwise pattern and DOC exhibiting a much narrower and unclear clockwise pattern.
18
2.5.3 End-Member Mixing Analysis
19
Table 2.4). Although only 6.1% of stream samples were outsiders, this model was rejected
because NSE was so poor for Ca2+ (0.44) and Mg2+ (0.23). A triangle of UP SH, T2, and
BW18 (M2), which tightly bound 87% of stream samples, significantly improved NSE
values for Ca2+, Mg2+, and Na+ (0.94-1.00) and marginally for SO42- and Sr. The remaining
trials, we supplemented either one or two of UP SH, T2, and BW18 end-members to
identify the best combination. Substituting UP SH with precipitation changed the modeling
results only slightly (e.g., M2 vs M9 and M7 vs M3). Any trials substituting BW18 with
nearby neighbors worsened the modeling (e.g., M4, M6, M8, and M10), indicating that
BW18 was necessary. Substituting T4 for either T1 or T2 produced poor NSE values for
Ca2+ and Mg2+ even if both UP SH and BW18 were used (e.g., M5). In summary, four best,
competing models were identified, which were composed of UP SH/precipitation, T1/T2,
and BW18 (M2, M3, M7, and M9).
The end-members in these four competing models were further evaluated from a
hydrologic perspective. UP SH, a soil lysimeter placed 10 cm below the soil surface in an
upslope location, represented soil water that was derived from infiltration of snowmelt and
rainfall(Kendall et al., 1999; Shanley et al., 2002). Both UP SH and precipitation worked
well in EMMA because they had similar ionic concentrations, particularly of Ca2+, Mg2+,
and Na+ (Table S1), and plotted very closely to each other in the mixing diagram (Figure
2.7). With hindsight, the mean fractional contribution of corresponding end-members
between alternative models had a difference < 3% with a R2 of almost 1.0 when either UP
SH or precipitation was used (M3 vs M7 and M2 vs M9).
T1 and T2 were co-located in a riparian zone and differed only in screened depth
(T1’s screen centered at 2.2 m while T2’s at 1.5 m). Comparing the results of M2 and M3,
NSE was almost identical for Ca2+, Mg2+, and Na+ and slightly better for SC and Sr from
M3 (Table 2.4). The total number of samples outside the triangle (therefore negative
solutions for at least one end-member) was about the same, but the number of samples with
a negative contribution of SOF from M3 was much higher than that from M2 (219 vs 147),
while that number of hillslope hollow (characterized by BW18) was much lower from M3
than from M2 (40 vs 122). Examination of ionic concentrations in stream samples outside
20
the triangle of M3 revealedthat samples below the UP SH-T1 side and above the UP SH-
BW18 side had abnormally high or low concentrations for one or two tracers compared to
samples collected within the day or a few days apart if there were no samples collected the
same day. However, abnormal ionic concentrations were hardly identified for those
unbounded on the T1-BW18 side. Interestingly, these unbounded points, with negative
contributions of SOF, occurred during very low flows (< 0.006 m3s-1), regardless of
whether T1 or T2 was used. Having more samples with negative contributions of SOF (zero
after corrected following Liu et al., 2004) during low flows was more plausible than having
more samples with negative contributions of hillslope hollow groundwater (which turned
out to dominate baseflow; see below). Therefore, M3 was superior to M2.
21
2.6 DISCUSSION
By combining the results from EMMA and the C-Q analysis, this study determines
the solute export patterns at multi-timescales and the key factors controlling annual or
seasonal shifts in export at W-9. The multi-annual C-Q relationships provide the “big
picture” of long-term solute concentration variation while the C-Q relationships at the
annual and seasonal timescales reveal increasingly more specific patterns in stream water
chemistry. In the multi-annual C-Q relationships, major geogenic ions and major nutrients
have contrasting patterns (Figure 2.3). The concentrations of major ions become diluted
during periods of high streamflow while the concentrations of major nutrients and some
trace elements enrich with increasing streamflow. This contrast is caused by the
distribution of solute concentrations throughout the soil profile (Table 2.A.1; Table 2.A.4;
Figure 2.2) and, subsequently, the conservative behavior and contrasting solute
concentrations in contributing end-members, consistent with Evans & Davies (1998).
22
increase in solute concentrations with depth creates a negative b, while a decrease creates
a positive b. This field observation-based finding agrees with a hypothetical numerical
modeling result by Chen, Tague, Melack, & Keller (2020). However, if source pools are
not large enough in volume, as suggested by Ackerer et al. (2020), the C-Q power-law may
be very weak. Source limitation may explain why Ba and Mn do not exhibit flushing
patterns despite having higher median concentrations in the shallow soils.
Contrasting C-Q patterns and solute vertical distribution throughout the soil profile
revealed which major flowpaths became active during high flow events. The flushing
patterns of NO3- and DOC and their high concentrations in surface soils indicate that a
significant amount of water flowed through or over shallow soils during periods of high
flow, as supported by EMMA. Several prior studies at SRRW indicated that SOF greatly
shapes NO3- and DOC export (Dunne and Black, 1970; Kendall et al., 1999; Shanley et al.,
2002; Sebestyen et al., 2008; Sebestyen et al., 2014; Sebestyen et al., 2019). SOF begins
heavily contributing to streamflow during storm events once precipitation wets the
nutrient-rich shallow soils (Dunne and Black, 1970; Kendall et al., 1999).
It is worth noting that though Si-Q had a significant power-law relation at multi-
annual scale based on p value, its b was extremely low (-0.05) like that of Cl- (Figure 2.4).
Silica has a significant biogenic cycling component and the dissolution of opal phytoliths
contained in leaf litter confound the shallow to deep concentration profile. When mixed,
this shallow source counters dilution to make the b-coefficient less negative. Additionally,
over the 24-year-long sampling period, SO42- concentrations in streamflow were
systematically lower after 2008 than before 2008 and were consistently under-estimated
using the EMMA results before 2008 but over-estimated after 2008 (Figure 2.10). This
decreasing trend is a result of the declining intensity of acid rain deposition in the
northeastern United States and has been observed at Hubbard Brook Experimental Forest
(Fuss, Driscoll, & Campbell, 2015). Using data from 2008 to 2015 in simulating stream
chemistry by the EMMA results, NSE becomes 0.56 and 0.55 for SO42- and Si,
respectively, significantly higher than their values for the entire period (M3 results in Table
4). Si appears to have become at least quasi-conservative after 2008.
23
2.6.2 Yearly Pattern
The variation of C-Q power-law coefficients over years is stronger for some solutes
than others (Figures 2.4 and 2.5). A key factor controlling this variability is whether the
solute behaves conservatively and whether solute concentrations in end-members are
relatively stable over time (Evans and Davis, 1998). The DTMM determined Ca2+, Mg2+,
and Na+ were conservative (Table 2.3). Using their median concentrations in EMMA, their
concentrations in streamflow were very well simulated (e.g., Ca2+ in Figure 2.10),
suggesting their concentrations in end-members were relatively constant over time.
Chloride and NO3- were among the solutes that were determined non-conservative
and had they some of the most interannual variation, with multiple insignificant C-Q annual
relationships. Both R2 and |b| values for Si became relatively high and stable after 2008,
consistent with Si behaving quasi-conservatively after 2008. DOC was assumed
conservative by Kendall et al. (1999), but this study’s EMMA results significantly
underestimated DOC concentrations (Figure 2.10), showing that DOC concentrations in
our SOF end-member were too low and also varied dramatically over time. A doubling of
DOC concentrations in SOF improved the mean values but not the overall temporal pattern.
Trace elements (other than Sr) were deemed to be non-conservative because their
concentrations in streamflow were not well simulated using the EMMA results (Table 2.4).
Their annual C-Q relationships had greater variability as well (Figures 2.4 and 2.5). Since
Al and Fe can be co-transported with DOC and colloidal materials, it is not surprising that
all three solutes have high interannual variability since all are controlled by similar
processes (Chorover, Derry, & McDowell, 2017). Non-conservative behavior degrades the
C-Q relationship and creates more variable concentrations due to factors such as instream
processes, biotic processes, plant uptake, anthropogenic additions, and colloid formation
and transport. The latter process tends to mobilize nutrients and trace metals (Mulholland
and Hill, 1997).
Furthermore, both NO3- and DOC concentrations are mainly delivered to the stream
via shallow flowpaths, with the amount of seasonal quick flow as one of the main factors
24
explaining seasonal variation of stream NO3- loading and the single factor explaining DOC
loading (Kendall et al., 1999; Shanley et al., 2002; Sebestyen et al., 2008, Sebestyen,
Boyer, & Shanley, 2009). Hence, the C-Q relationships deteriorate from variability in
production and uptake/degradation of NO3- and DOC due to weather and wetting/drying
cycles.
Like many areas of the US, there is high nitrogen deposition in SRRW. Nitrate can
be quickly transported to the stream during storms and snowmelt without being
transformed or immobilized by biological uptake (Campbell et al., 2004; Sebestyen et al.,
2008; Sebestyen et al., 2019). New N inputs may be highly variable and thus further
weaken the C-Q relationship. Prior studies from SRRW show that NO3- is source limited,
meaning that NO3- becomes progressively depleted over long events such as snowmelt
(Sebestyen et al., 2008), further eroding the C-Q relation.
For those solutes with high interannual variation, b-coefficients more frequently
switched from flushing or diluting to chemostatic behaviors than from flushing to diluting
(or vice versa). Only Al had significant (p < 0.05) annual C-Q relationships where the b-
coefficient was positive one year and negative in another. This result reveals that, for most
solutes, flushing or diluting patterns are relatively stable, rarely switching between
chemodynamic patterns.
It is important to note that neither the b-coefficient nor CVC/CVQ were necessarily
good indicators of mixing and strength of the C-Q relationship, as is seen in the multi-
annual and annual C-Q relationships. Results for Ca2+, Mg2+, and Na+ showed that strong
mixing does not require a b near -1 (Godsey, 2009). Here, these major ions had high R2
values and small p values, but |b| was only 0.11-0.15 (Tables 2.1 and 2.2), which would be
considered chemostatic by the thresholds in Godsey et al. (2009). Meanwhile, CVC/CVQ is
also low for these major ions and higher for NO3- and DOC (Tables 2.1 and 2.2). DOC has
a b of 0.26 and CVC/CVQ of 0.54, but it cannot be predicted by conservative mixing and
its annual C-Q varies significantly (Figure 2.5). The EMMA results showed that C-Q
25
relationships with low b-coefficients can still be strong when there is strong mixing among
distinct end-members with unique concentrations.
During the snowmelt period, the C-Q relationships showed that the median
intensity of solute export for most solutes changed between the rising and falling limb of
snowmelt. Most solutes had smaller median |b| coefficients on the rising limb than the
falling limb of snowmelt (Figure 2.6). Only Si and Al had larger median |b| coefficients on
the rising limb compared to the falling limb. The EMMA results indicate that the
proportional contributions of hillslope hollow groundwater decreased and that of SOF and
riparian groundwater increased during the falling limb of snowmelt (Figure 2.8). The
difference in solute concentrations of major ions (except NO3-) was significant between
SOF (characterized by UP SH) and either riparian groundwater (T1) or hillslope hollow
groundwater (BW18) (Figure 2.2). The increase in fractional contribution of SOF during
the falling limb (Figure 2.8) increases the |b| value for all solutes except Si and Al (Figure
2.6). Because solute concentrations in groundwater from riparian zone and hillslope hollow
are very close (Figure 2.2), relative changes in fractional contributions between riparian
and hillslope hollow groundwater do not change much stream chemistry. Thus, a hysteresis
was not observed for most solutes. Similar to DOC, NO3- is an exception because its
concentration was much higher in SOF based on the above discussion and also somewhat
different between riparian groundwater and hillslope groundwater (Figure 2.2), resulting
26
in a subtle hysteresis (Figure 2.B.2). Certain studies have found hysteretic patterns in data
from individual years (Kendall et al., 1999; Sebestyen et al., 2014). Since our study
combined the data from all years to study snowmelt hysteresis, it is possible that the
hysteresis seen within individual years becomes obscured when looking at all data due to
interannual variation.
Changing C-Q patterns between seasons also means that timing of solute export
and nutrient availability might shift in a future climate as snow melt occurs earlier, dormant
season ends earlier, and the growing season becomes longer. Coupled with results from
existing studies on precipitation patterns and influence of storm events on stream nutrient
loadings, the observed seasonal patterns from my study of C-Q analysis and EMMA
reinforce the theory that climate change effects on the frequency and magnitude of storm
events will drive future patterns of stream nutrient loadings (Sebestyen et al., 2009). With
a warmer climate, solute export during snowmelt will continue to occur earlier in the year.
Additionally, during the growing season, there will be greater proportional contributions
of SOF as more intense and precipitation events are projected in the future (Hayhoe et al.
2007), causing more water to flow through quick flowpaths. Increased storm intensity will
also increase DOC export since DOC had a flushing pattern, consistent with Sebestyen et
al. (2008).
The observed seasonal variation in C-Q relationships also has implications for
sampling regime and data analysis. If sampling is not conducted representatively
throughout the year, researchers can potentially draw inaccurate conclusions from the
constructed relationships. This seasonal variation particularly poses a challenge for
collecting representative NO3- and DOC samples from surficial soils.
2.7 CONCLUSION
The power-law relationship between concentration and discharge in streamflow is
primarily controlled by vertical distribution of solutes in source waters, an increase with
depth generating a negative b-coefficient (dilution) while a decrease with depth generating
a positive b-coefficient (flushing). The strength (R2) and direction (b-coefficient) of the
27
power-law depends on (1) conservative behavior of solutes, (2) distinctiveness of solute
concentrations among source waters, (3) temporal variability of solute concentrations in
source waters. Solutes that are conservative with distinct and relatively constant
concentrations in source waters, show a strong and persistent power-law relation over
variable timescales such as Ca2+, Mg2+, Na+, and Sr of this study. The strength of the power-
law cannot be solely measured by the magnitude of b-coefficient and CVQ/CVC. A |b|-
coefficient as low as 0.15, like those of Ca2+, Mg2+, Na+, and Sr in this study, could still be
considered a strong power-law relationship and interpreted primarily as a result of mixing.
The power-law relationships of some solutes such as nutrients in this study could change
over timescales due to non-conservative behavior, variable source stores, and/or
intermittent flow pathways. With a change in snowmelt timing and precipitation intensities
in the future, concentration-discharge relationship in streamflow could change, particularly
at the seasonal scale.
2.8 REFERENCES
Ackerer J, Jeannot B, Delay F, Weill S, Lucas Y, Fritz B, … & Chabaux F. (2020).
Crossing hydrological and geochemical modeling to understand the spatiotemporal
variability of water chemistry in a headwater catchment. Hydrology and Earth
System Science. 24, pp. 3111–3133. doi:10.5194/hess-24-3111-2020
Armfield J R, Perdrial J N, Gagnon A, Ehrenkranz J, Perdrial N, Cincotta M, … & Ryan
P. (2019). Does stream water composition at Sleepers River in Vermont reflect
dynamic changes in soils during recovery from acidification? Frontiers in Earth
Science. 6, pp. 1–13. doi:10.3389/feart.2018.00246
Bieroza M Z, Heathwaite A L, Bechmann M, Kyllmar K, & Jordan P. (2018). The
concentration-discharge slope as a tool for water quality management. Science of the
Total Environment. 630, pp. 738–749. doi:10.1016/j.scitotenv.2018.02.256
Bishop K H, Grip H, & O'Neill A. (1990). The origins of acid runoff in a hillslope during
storm events. Journal of Hydrology. 116(1-4), pp. 35-61. doi:10.1016/0022-
1694(90)90114-D
28
Bishop K H (1991). Episodic increases in stream acidity, catchment flow pathways and
hydrograph separation, dissertation thesis. pp. 246, Cambridge University,
Cambridge, England.
29
Evans C, & Davies T D. (1998). Causes of concentration/discharge hysteresis and its
potential as a tool for analysis of episode hydrochemistry. Water Resources
Research, 34(1), pp. 129–137.
Field C B, Barros V R, Mach K J, Mastrandrea M D, Van Aalst M, Adger W N, …. &
Yohe G W. (2014). Climate Change 2014: Impacts, Adaptation, and Vulnerability.
Part A: Global and Sectoral Aspects. Contribution of Working Group II to the Fifth
Assessment Report of the Intergovernmental Panel on Climate Change. Technical
Summary, pp. 35–94.
Frisbee M D, Phillips F M, Campbell A R, Liu F, and Sanchez S A. (2011), Streamflow
generation in a large, alpine watershed in the southern Rocky Mountains of
Colorado: Is streamflow generation simply the aggregation of hillslope runoff
responses? Water Resources Research, 47, W06512, doi:10.1029/2010WR009391
Fuss C B, Driscoll C T, & Campbell J L. (2015). Recovery from chronic and snowmelt
acidification: Long-term trends in stream and soil water chemistry at the Hubbard
Brook Experimental Forest, New Hampshire, USA. Journal Geophysical Research,
120, pp. 2360–2374, doi:10.1002/2015JG003063
Godsey S E, Kirchner J W, & Clow D W. (2009). Concentration-discharge relationships
reflect chemostatic characteristics of US catchments. Hydrological Processes, 23,
pp. 1844–1864. doi:10.1002/hyp.7315
Hall LM. 1959. The geology of the St. Johnsbury Quadrangle, Vermont and New
Hampshire. Bulletin No. 13. Montpelier: Vermont Development Commission. pp.
105.
Hall F R. (1970). Dissolved solids-discharge relationships 1: mixing models, Water
Resources Research. 6(3), pp. 845-850.
Hall, F. R. (1971). Dissolved solids-discharge relationships 2: applications to field data,
Water Resources Research. 7(3), pp. 591-601.
Hayhoe K, Wake C P, Huntington T G, Luo L, Schwartz M D, Sheffield J, … & Wolfe
D. (2007). Past and future changes in climate and hydrological indicators in the US
Northeast. Climate Dynamics. 28, pp. 381–407.
Herndon E M, Dere A L, Sullivan P L, Norris D, Reynolds B, & Brantley S L. (2015).
Landscape heterogeneity drives contrasting concentration-discharge relationships in
shale headwater catchments. Hydrology and Earth System Sciences, 19, pp. 3333–
3347. doi:10.5194/hess-19-3333-2015
Hidalgo H G, Das T, Dettinger M D, Cayan D R, Pierce D W, Barnett T P, …& Nozawa
T. (2009). Detection and attribution of streamflow timing changes to climate change
in the western United States. 22, pp. 3838–3855. doi.org/10.1175/2009JCLI2470.1
Hofmeister K L, Nave L E, Drevnick P, Veverica T, Knudstrup R, Heckman K A, …&
Walter M T. (2019). Seasonal dynamics and exports of elements from a first-order
stream to a large inland lake in Michigan. Hydrological Processes, 33, pp. 1476–
1491. doi:10.1002/hyp.13416
30
Hooper, R. P. (2003). Diagnostic tools for mixing models of stream water chemistry.
Water Resources Research, 39(3), pp. 1–13. doi:10.1029/2002WR001528
Hornbeck J, Bailey S, Buso D, & Shanley J. (1997). Streamwater chemistry and nutrient
budgets for forested watersheds in New England: Variability and management
implications. Forest Ecology and Management. 93, pp. 73–89. doi:
10.1016/S0378- 1127(96)03937-0
Johnson N M, Likens G E, Bormann F H, Fisher D W, & Pierce R S. (1969). A working
model for the variation in stream water chemistry at the Hubbard Brook
Experimental Forest, New Hampshire, Water Resources Research. 5(6), pp. 1353-
1363.
Kendall K A, Shanley J B, & McDonnell J J. (1999). A hydrometric and geochemical
approach to test the transmissivity feedback hypothesis during snowmelt. Journal of
Hydrology, 219, pp. 188–205.
Krause S, Lewandowski J, Grimm N B, Hannah D M, Pinay G, McDonald K, … & Turk
V. (2017). Ecohydrological interfaces as hot spots of ecosystem processes. Water
Resources Research, 53, pp. 6359–6376. doi:10.1002/2016WR019516
Liu F, Bales R C, Conklin M H, & Conrad M E. (2008). Streamflow generation from
snowmelt in semi-arid, seasonally snow-covered, forested catchments, Valles
Caldera, New Mexico. Water Resources Research, 44, W12443. pp. 1–13.
doi:10.1029/2007WR006728
Liu F, Hunsaker C, & Bales R C. (2013), Controls of streamflow generation in small
catchments across snow-rain transition in the Southern Sierra, California.
Hydrological Processes, 27, pp. 1959–1972. doi:10.1002/hyp.9304
Liu F, Williams M W, & Caine N. (2004). Source waters and flow paths in an alpine
catchment, Colorado Front Range, United States. Water Resources Research, 40,
W09401. pp. 1–16. doi:10.1029/2004WR003076
Liu F, Conklin M H, & Shaw G D. (2017). Insights into hydrologic and hydrochemical
processes based on concentration-discharge and end-member mixing analyses in the
mid-Merced River Basin, Sierra Nevada, California, Water Resoures Research, 53,
pp. 832– 850, doi:10.1002/2016WR019437
Liu F, Lerch R N, Yang J, and Peters G. (2020). Determining streamflow pathways using
geochemical tracers in a claypan watershed, Hydrological Processes, 34, pp. 2494–
2509. doi:10.1002/hyp.13743
Markstrom S L, & Hay L E. (2009). Integrated watershed scale response to climate
change for selected basins across the United States. Water Resources. 11(2), pp. 8–
10.
Maher K. (2011). The role of fluid residence time and topographic scales in determining
chemical fluxes from landscapes. Earth and Planetary Science Letters, 312, pp. 48 –
58.
31
Mayer B, Shanley J B, Bailey S W, & Mitchell M J. (2010). Identifying sources of
stream water sulfate after a summer drought in the Sleepers River watershed
(Vermont, USA) using hydrological, chemical, and isotopic techniques. Applied
Geochemistry, 25, pp. 747–754. doi:10.1016/j.apgeochem.2010.02.007
Moater F, Abbott B, Minaudo C, Curie F, & Pinay G. (2017). Elemental properties,
hydrology, and biology interact to shape concentration-discharge curves for carbon,
nutrients, sediment, and major ions. Water Resources Research, 53, pp. 1270–1287.
doi:10.1002/2016WR019635
Mulholland P J, & Hill W R. (1997). Seasonal patterns in streamwater nutrient and
dissolved organic carbon concentrations: Separating catchment flow path and in-
stream effects. Water Resources Research. 33 (6), pp. 1297 – 1306.
Ohte N, Sebestyen S D, Kendall C, Shanley J B, Wankel S D, Doctor D H, & Boyer E
W. (2004). Tracing sources of nitrate in snowmelt runoff using a high-resolution
isotopic technique. Geophysical Research Letters. 31, L21506.
doi:10.1029/2004GL020908
Pierce D W, & Cayan D R. (2013). The uneven response of different snow measures to
human-induced climate warming. Journal of Climate, 26, 4148–4167.
doi:10.1175/JCLI-D-12-00534.1
Pinder G F & Jones J F. (1969). Determination of the ground-water component of peak
discharge from the chemistry of total runoff. Water Resources Research. 5(2), pp.
438-445. doi:10.1029/WR005i002p00438
Rantz S E. (1982), Measurement and computation of streamflow: Volume 1:
Measurement of stage and discharge, Water-Supply Paper. 2175, pp. 284, U.S.
Geological Survey, Washington, D. C.
Sebestyen S D, Shanley J B, Boyer E W, Kendall C, and Doctor DH. (2014). Coupled
hydrological and biogeochemical processes controlling variability of nitrogen
species in streamflow during autumn in an upland forest. Water Resources Research,
50, pp. 1569–1591. doi:10.1002/2013WR013670
Sebestyen S D, Boyer E W, & Shanley J B. (2009). Responses of stream nitrate and DOC
loadings to hydrological forcing and climate change in an upland forest of the
northeastern United States. Journal of Geophysical Research, 114, G02002.
doi:10.1029/2008JG000778
Sebestyen S D, Boyer E W, Shanley J B, Kendall C, Doctor D H, Aiken G R, & Ohte N.
(2008). Sources, transformations, and hydrological processes that control stream
nitrate and dissolved organic matter concentrations during snowmelt in an upland
forest. Water Resources Research. 44, W12410. doi:10.1029/2008WR006983
Sebestyen S D, Ross D S, Shanley J B, Elliott E M, Kendall C, Campbell J L, …&
Williard K W J. (2019), Unprocessed atmospheric nitrate in waters of the Northern
Forest Region in the U.S. and Canada. Environmental Science and Technology, 53,
pp. 3620-3633. doi:10.1021/acs.est.9b01276
32
Shanley J B, Hjerdt K N, McDonnell J J, & Kendall C. (2003). Shallow water table
fluctuations in relation to soil penetration resistance. Ground Water, 41(7), pp. 964–
972.
Shanley J B, Kendall C, Smith T E, Wolock D M, & McDonnell J J. (2002). Controls on
old and new water contributions to stream flow at some nested catchments in
Vermont, USA. 16, pp. 589–609. doi:10.1002/hyp.312
Shanley J B, Krám P, Hruška J, & Bullen T D. (2004). A biogeochemical comparison of
two well-buffered catchments with contrasting histories of acid deposition. Water,
Air, and Soil Pollution 4, pp. 325–342.
Shanley J B, Mayer B, Mitchell M J, & Bailey S W. (2008). Seasonal and event
variations in δ34S values of stream sulfate in a Vermont forested catchment:
Implications for sulfur sources and cycling. Science of the Total Environment, 404,
pp. 262–268. doi:10.1016/j.scitotenv.2008.03.020
Shanley J B, Sebestyen S D, MacDonnell J J, McGlynn B L, & Dunne T. (2015). Water’s
Way at Sleepers River watershed – revisiting flow generation in a post-glacial
landscape, Vermont USA. Hydrological Processes. 29, pp. 3447 – 3459.
doi:10.1002/hyp.10377
Shanley J B, Sundquist E T, Kendall C, 1995. Water, Energy, and Biogeochemical
Budget research at Sleepers River Research Watershed, Vermont. US Geological
Survey Open- File Report, pp. 94-475, 22
Smith T E. 1997. Comparison of runoff flow paths in four basins at Sleepers River
Research Watershed, Vermont. Hydrology and Water Resources. University of New
Hampshire: Durham, NH, USA; 197.
Thompson S E, Basu N B, Lascurain J, Aubeneau A, & Rao P S C. (2011). Relative
dominance of hydrologic versus biogeochemical factors on solute export across
impact gradients. Water Resources Research, 47, W00J05,
doi:10.1029/2010WR009605
Winnick M J, Carroll R W H, Williams K H, Maxwell R M, Dong W, & Maher K.
(2017). Snowmelt controls on concentration-discharge relationships and the balance
of oxidative and acid-base weathering fluxes in an alpine catchment, East River,
Colorado. Water Resources Research, 53(3). pp. 2507 – 2523.
doi:10.1002/2016WR019724
Zimmer M A, Pellerin B, Burns D A, & Petrochenkov G. (2019). Temporal variability in
nitrate-discharge relationships in large rivers as revealed by high-frequency data.
Water Resources Research, 55, pp. 973–989.doi:10.1029/2018WR023478
33
2.9 TABLES AND FIGURES
2.9.1 Tables
Table 2.1. Multi-annual C-Q coefficients, CVC/CVQ, and root mean squared error
(RMSE) for the 15 solutes.
Table 2.2. Summary statistics of annual C-Q relationships from years with significant
relationships (p < 0.05).
34
Table 2.3. The results of diagnostic tools of mixing models using concentrations of major
ions (n = 2420)
35
Table 2.4. Evaluation of alternative end-member mixing models using Nash-Sutcliffe Efficiency (NSE) and percent of samples outside
the triangle.
Model End- Out-
Code members siders SC Ca2+ Mg2+ Na+ K+ Cl- NO3- SO42- DOC Si Sr Al Fe
(μS (mg L- (mg L- (mg L- (mg L- (mg L- (mg L- (mg L- (mg L- (mg L- (μg L- (μg L- (μg L-
(%) cm-1) 1) 1) 1) 1) 1) 1) 1) 1) 1) 1) 1) 1)
UP SH;
BW39A;
M1 BW6 6.1 0.64 0.44 0.23 0.98 0.01 -0.09 -55.9 -0.32 -1.11 -2.34 0.79 -10.66 -3.13
UP SH;
M2 T2; BW18 13.1 0.26 0.96 0.94 1 0.06 -0.01 0.05 0.26 -0.11 -0.99 0.82 -4.57 -0.71
UP SH;
M3 T1; BW18 12.7 0.44 0.95 0.93 1 -0.67 -0.03 0.07 0.27 -0.08 -0.3 0.86 -4.48 -0.7
UP SH;
M4 T2; BW6 11.7 0.41 0.72 0.62 0.99 0 -0.05 -23.49 -0.12 -1.01 -1.51 0.73 -9.57 -2.59
UP SH;
36
M5 T4; BW18 5.1 0.62 0.43 0.22 0.98 -1.56 -0.01 -69.48 0.01 -0.6 0.03 0.68 -7.17 -1.2
UP SH;
M6 T1; BW26 9.2 0.52 0.41 0.19 0.98 -0.45 0 0.07 0.38 -0.19 -0.64 0.26 -5.57 -1.17
PPT; T1;
M7 BW18 12.1 0.4 0.96 0.94 1 0.44 -0.06 0.03 0.27 -0.81 -0.7 0.85 -0.09 -0.16
UP SH;
T1;
M8 8/9/2011 9.4 0.92 0.89 1 -0.3 -0.56 -0.09 -0.29 -3.54 -0.7 0.75 -6.15 -3.59
PPT; T2;
M9 BW18 12.2 0.18 0.96 0.94 1 0.39 -0.07 -0.01 0.26 -0.79 -0.9 0.83 -0.19 -0.22
PPT; T2;
M10 BW17 10.9 -1.15 0.89 0.85 1 0 -0.12 -0.17 0.17 -0.66 -0.92 -0.65
Note: the simulation cannot be done for Ba and Mn because their concentrations in hillslope (both bench and hollow) samples were
not available.
2.9.2 Figures
37
Figure 2.1. Map of sampling locations including W-9 stream gauge, precipitation gauge (R29), soil lysimeters, piezometers (T-#), and
groundwater wells (BW-##).
38
Figure 2.2. Depth profiles of Ca2+, NO3-, and DOC concentrations in soil lysimeters, piezometers, and selected groundwater wells
using all data collected during the 1992-2015 water years. Note that depth increases toward the right of the graphs and top of the page
within each sample group. Solute concentrations in precipitation, high (5%) and low (95%) flows of all years are shown for
comparison. Data for all groundwater wells are summarized in Table S1.
39
Figure 2.3. C-Q relationships at the multi-annual timescales for water years 1992-2015. All relationships were calculated on the log10-
log10 scale with red indicating insignificant C-Q relationships (p < 0.05).
40
Figure 2.4. Temporal variation of the R2 values from the yearly C-Q relationships. All relationships were calculated on the log10-log10
scale with red indicating insignificant C-Q relationships (p < 0.05).
41
Figure 2.5. Temporal variation of the b-coefficients from the yearly C-Q relationships. All relationships were calculated on the log10-
log10 scale. Red indicates an insignificant C-Q relationship (p < 0.05), and the light red shaded section indicates the range of b where
the relationship is considered chemostatic.
42
Figure 2.6. Boxplots of the b-coefficients from the seasonal C-Q relationships during the four periods: dormant, rising limb of
snowmelt, falling limb of snowmelt, and growing season. All relationships were calculated on the log10-log10 scale, and only those
with significant relationships (p < 0.05) were included in the statistics. The red line indicates the b-coefficient from the multiannual C-
Q relationships.
43
Figure 2.7. Mixing diagram using the first two principal components (PC1 and PC2) constructed from Ca2+, Mg2+, and Na+. Triangles
show examples of alternative models with varying end-members. Sites with labels are discussed in the text. The final selected end-
members were UP SH, T1, and BW18.
44
Figure 2.8. Boxplots of fractional contributions of saturation overland flow, riparian groundwater, and hillslope groundwater from
hollows across the four seasons
45
Figure 2.9. Variation of b-coefficients with geogenic solute concentration (SC, Ca2+, Mg2+, Na+, K+, Cl-, SO42-, Si, and Sr) ratios
between T1 and UP SH (a) and BW18 and UP SH (b).
46
Figure 2.10. Measured (blue) vs modelled (red) stream concentrations for Ca2+, SO42-, NO3-, and DOC using the EMMA results. Note
that an additional curve (orange) was shown for modelled NO3- results with hypothetically increased concentrations in saturation
overland flow (see text).
2.10 APPENDIX 2.A SUPPLEMENTAL TABLES
Supplemental Table 2.A.1: Summary of solute concentrations in all types of samples collected during water years 1992-2015,
including 25%, 50%, 75% quartiles, the sample period, and the number of samples.
SC Ca2+ Mg2+ Na+ K+ Cl- NO3- SO42- DOC Si Sr Al Ba Fe Mn
Station pH (μS (mg (mg (mg (mg (mg (mg (mg (mg (mg (μg (μg (μg (μg (μg
cm-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1)
Median 7.6 108 21.6 1 0.6 1.1 0.3 0.7 6.3 1.6 5.3 87.3 15.4 8.1 8.6 1.4
25% 7.3 92.9 17.6 0.8 0.5 0.9 0.3 0.5 5.1 1.2 4.6 70.9 7.1 6.7 4.2 0.6
75% 7.7 122.2 24.8 1.2 0.7 1.3 0.4 1.1 7.6 2.6 6 103.5 30 9.9 17 3.3
W-9
Date 1991- 1991- 1991- 1991- 1991- 1991- 1991- 1991- 1991- 1991- 1991- 2008- 1993- 2008- 2008- 2008-
Range 2015 2007 2015 2015 2015 2015 2015 2015 2015 2015 2015 2015 2015 2015 2015 2015
n 1760 693 2336 2336 2336 2336 2336 2336 2336 1837 2321 988 1431 988 987 986
Median 4.7 12.6 0.1 0 0.1 0.1 0.1 0.9 0.8 0.8 0.1 0.7 6 0.7 5.1 2
25% 4.4 6.9 0.1 0 0 0 0.1 0.5 0.4 0.5 0 0.4 2.7 0.4 2.6 1
75% 5.2 20.4 0.2 0 0.1 0.1 0.2 1.6 1.6 1.4 0.2 1.3 10 2.4 11.2 3
PPT
Date 1991- 1991- 1991- 1991- 1991- 1991- 1991- 1991- 1991- 1991- 1991- 1999- 1991- 2007- 1991- 2007-
Range 2020 2020 2020 2020 2020 2020 2020 2020 2020 2020 2020 2018 2020 2017 2020 2020
47
n 1090 967 1795 1777 1736 1730 1722 1756 1794 1428 1400 378 1286 265 1511 518
Median 6.9 45 6.8 0.5 0.7 0.2 0.3 0.2 3.9 1.3 7 35.1 31.6 2.7 14.7 0.2
25% 6.5 32 4.6 0.4 0.5 0.1 0.2 0.1 3 0.7 6.3 26.4 7.1 2.3 1.8 0.1
LYS
75% 7.2 57.8 7.3 0.6 0.8 0.3 0.4 2.3 4.8 3.5 7.7 41.5 81.9 5 38.5 0.5
LO
Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1992- 2008- 1999- 2008-
SH
Range 2013 2002 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013
n 71 62 95 104 104 94 111 90 111 103 103 74 67 47 60 46
Median 6.9 52.4 6.4 0.5 0.8 0.2 0.3 0.1 4.1 0.7 7.1 37 6.3 2.5 1.6 0.2
25% 6.6 43.3 5.8 0.5 0.7 0.1 0.2 0 3.8 0.5 6.5 31.6 4 2 0.6 0.1
LYS
75% 7.2 56.3 7.3 0.6 0.8 0.3 0.3 5.2 5 0.9 7.8 40 12 7.3 7.2 0.3
LO
Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1992- 2008- 1999- 2008-
DP
Range 2013 2002 2013 2013 2013 2013 2013 2012 2013 2013 2013 2013 2013 2013 2013 2013
n 78 56 93 103 103 91 114 72 114 107 103 77 69 52 47 49
Median 5.5 15.8 1.2 0.2 0.1 0.5 0.2 0.8 1 8.5 3.9 7.1 267.4 16.4 47 2.9
25% 5.2 12.2 1 0.1 0.1 0.2 0.2 0.4 0.7 6.8 2.7 6.2 207.4 13 29.4 1.8
LYS
75% 5.7 21.8 1.6 0.2 0.2 0.8 0.5 1.2 1.9 10.4 5.2 8.6 345 20 66.4 4.9
MI
Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 2003- 1992- 2003- 2003- 2008-
SH
Range 2016 2002 2015 2015 2015 2015 2016 2016 2016 2016 2015 2015 2015 2015 2015 2015
n 56 19 61 66 63 63 76 72 75 74 61 57 51 50 56 49
SC Ca2+ Mg2+ Na+ K+ Cl- NO3- SO42- DOC Si Sr Al Ba Fe Mn
Station pH (μS (mg (mg (mg (mg (mg (mg (mg (mg (mg (μg (μg (μg (μg (μg
cm-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1)
Median 6.4 15.5 2.1 0.2 0.2 0.1 0.3 0.4 2.4 3.9 5.4 18.4 81.3 8.3 15.4 1.4
25% 6.2 15.4 1.9 0.2 0.2 0.1 0.2 0.2 1.8 2.2 4.5 15.3 30 4.8 8.8 1
LYS
75% 6.5 19.5 2.4 0.2 0.3 0.2 0.3 2 2.6 4.4 6 20.1 90.5 10.1 37.1 1.6
MI
Date 1992- 1993- 199 - 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 2003- 1992- 2003- 2003- 2008-
MI
Range 2011 1998 2011 2011 2011 2011 2011 2011 2011 2011 2011 2011 2011 2011 2011 2011
n 9 5 12 13 13 13 15 15 15 13 12 10 9 8 9 7
Median 6.6 24.2 3 0.2 0.3 0.1 0.3 0.3 2.7 1.6 5.3 20.9 30.7 2.7 3.1 0.5
25% 6.3 19.3 2.6 0.2 0.3 0 0.2 0.2 2.1 1.3 4.8 19.5 26.4 1.9 1.8 0.4
LYS
75% 6.7 25.7 3.3 0.2 0.4 0.1 0.3 1.1 3.3 1.7 6.2 25.4 35.1 6.9 6.5 0.8
MI
Date 1992- 1993- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1992- 2008- 1999- 2008
DP
Range 2015 1999 2015 2015 2015 2015 2015 2015 2015 2015 2015 2015 2015 2015 2015 2015
n 23 17 29 32 32 26 34 30 34 30 30 25 22 18 25 18
Median 5.5 16.1 0.9 0.2 0.1 1 0.2 1 1.4 9.6 3.5 6.2 340.4 5.7 83.6 4.4
25% 5.3 12 0.7 0.2 0.1 0.7 0.2 0.5 0.8 5.5 1.9 4.8 219 5.1 65 1.9
LYS
75% 6 28.2 1.3 0.3 0.2 1.4 0.4 1.9 2 11.9 5.2 8.1 437.8 7.7 110 7.8
UP
Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1992- 2008- 1999- 2008-
SH
48
Range 2016 2002 2015 2015 2015 2015 2016 2016 2016 2016 2015 2015 2015 2015 2015 2015
n 82 40 89 99 98 99 119 114 120 106 92 71 72 58 72 58
Median 6.4 15.2 1.2 0.3 0.2 0.1 0.3 0.2 2.8 1 4.9 10.1 24.1 6.5 1.7 1.6
25% 6 12.8 1.1 0.2 0.2 0.1 0.2 0.1 2.6 0.8 4.7 8.4 17.7 4 0.9 1.6
LYS
75% 6.6 19.8 1.5 0.3 0.3 0.2 0.4 2.7 3.2 1.4 6.1 10.6 31.9 9.6 2.9 2
UP
Date 1992- 1993- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1992- 2008- 1999- 2008-
MI
Range 2011 1999 2011 2011 2011 2011 2011 2011 2011 2011 2011 2011 2011 2011 2011 2011
n 15 15 23 26 25 22 27 24 27 23 23 14 15 8 12 8
Median 6.6 24.3 2.6 0.3 0.6 0.2 0.3 0.2 3.3 0.5 5.7 13.9 5.7 1.1 3.7 0.4
25% 6.4 18.6 2.4 0.3 0.5 0.2 0.2 0.1 2.8 0.4 5.4 12.3 3.3 0.9 1.2 0.2
LYS
75% 7 27.1 2.9 0.3 0.6 0.3 0.3 4.8 3.7 0.6 6.1 15.8 24.8 1.4 7.8 0.6
UP
Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1992- 2008- 1999- 2008-
DP
Range 2016 1999 2015 2015 2015 2015 2016 2016 2016 2016 2015 2015 2015 2015 2015 2015
n 53 35 68 72 72 70 80 48 80 68 70 56 45 37 37 35
SC Ca2+ Mg2+ Na+ K+ Cl- NO3- SO42- DOC Si Sr Al Ba Fe Mn
Station pH (μS (mg (mg (mg (mg (mg (mg (mg (mg (mg (μg (μg (μg (μg (μg
cm-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1)
Median 7.6 216.5 28.8 6.3 4.7 2.1 0.4 1.3 9.7 1.6 11.1 69.9
25% 7.5 209.8 27.9 4.9 4.6 2.1 0.4 1.3 8.6 1.3 10.7 64.1
BW- 75% 7.7 223.3 29.7 7.7 4.8 2.1 0.4 1.3 10.8 1.9 11.6 75.7
1 Date
Range 1992 1992 1992 1992 1992 1992 1992 1992 1992 1992 1992 1992
n 2 2 2 2 2 2 2 1 2 2 2 2
Median 8 373 43.3 4 3.3 2.6 0.6 0.2 9.8 3.6 10.8 74.7
25% 7.9 373 36 3.9 2.8 2.2 0.5 0.2 9.2 2.8 10.6 64.5
BW- 75% 8.1 373 50.5 4.2 3.8 2.9 0.6 0.3 10.5 4.3 11 85
2 Date
Range 1992 1992 1992 1992 1992 1992 1992 1992 1992 1992 1992 1992
n 2 1 2 2 2 2 2 2 2 2 2 2
Median 7.58 232 36.7 3.64 1.52 2.15 0.39 1.49 9.09 2.91 7.62 198.9 48.03 1.5
25% 7.45 232 31.84 2.4 0.94 1.74 0.39 0.99 7.17 2.25 7.2 198.9 46.95 1.5
BW- 75% 7.71 232 43.37 3.67 1.68 2.34 0.51 6.62 9.46 4.84 7.72 198.9 49.11 1.5
3 Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992-
49
Range 1992 1992 2003 2003 2003 2003 2003 2003 2003 2003 2003 2003 1992 2003
n 2 1 3 3 3 3 3 3 3 3 3 1 2 1
Median 8.18 187.9 33.26 1.58 0.93 1.65 0.42 6.36 9.26 0.04 9.44 151.9 54.5 3.7
25% 8.18 181.8 33.23 1.55 0.88 1.63 0.38 3.64 9.19 0.04 9.34 151.9 54.5 3.7
BW- 75% 8.18 193.8 33.28 1.6 0.98 1.68 0.45 9.07 9.32 0.05 9.54 151.9 54.5 3.7
4 Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992-
Range 1992 1996 2003 2003 2003 2003 2003 2003 2003 2003 2003 2003 1992 2003
n 1 6 2 2 2 2 2 2 2 2 2 1 1 1
Median 7.54 150 24.14 0.9 0.88 1.45 0.46 6.04 9.18 0.58 8.22 110.4 36.97
25% 7.54 150 24.05 0.87 0.88 1.44 0.45 3.2 9.13 0.35 8.18 110.4 36.97
BW- 75% 7.54 150 24.22 0.92 0.88 1.46 0.48 8.89 9.22 0.81 8.26 110.4 36.97
5 Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992-
Range 1992 1992 2003 2003 2003 2003 2003 2003 2003 2003 2003 2003 1992
n 1 1 2 2 2 2 2 2 2 2 2 1 1
SC Ca2+ Mg2+ Na+ K+ Cl- NO3- SO42- DOC Si Sr Al Ba Fe Mn
Station pH (μS (mg (mg (mg (mg (mg (mg (mg (mg (mg (μg (μg (μg (μg (μg
cm-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1)
Median 8 213 38.7 2.5 1 1.6 0.3 11.2 8.7 0.4 9.1 185.1 70.4 25
25% 8 169 38 2.4 1 1.5 0.3 4.5 8.6 0.3 8.7 179.6 70.4 25
BW- 75% 8 219.5 39.4 2.5 1 1.7 0.4 12.2 9.4 0.4 9.7 186.9 70.4 25
6 Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1993- 1992- 2002- 1992-
Range 1993 1993 2004 2004 2004 2004 2004 2004 2004 2004 2004 2004 1993 2003
n 3 3 9 9 9 9 9 9 9 4 9 6 2 1
Median 176.5 30.8 1.6 1 1.4 0.4 6.7 9.4 0.6 7.3 130.7 77.4
25% 175.2 30.3 1.6 1 1.4 0.4 4.5 9 0.6 7.2 130.7 77.4
BW- 75% 176.8 31.3 1.6 1 1.4 0.4 9 9.8 0.6 7.4 130.7 77.4
7 Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992-
Range 1996 2003 2003 2003 2003 2003 2003 2003 1992 2003 2003 1992
n 6 2 2 2 2 2 2 2 1 2 1 1
Median 7.7 128 21.3 1.3 0.6 1.2 0.3 2 7.6 0.5 6.4 92.5 44
25% 7.7 120.5 20 1.2 0.6 1.1 0.3 1.4 7.1 0.2 6.2 92.5 38.3
BW- 75% 7.8 137.5 22.5 1.3 0.7 1.2 0.3 2.3 8.2 0.5 7.1 92.5 88.2
8 Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992-
50
Range 1993 1993 2003 2003 2003 2003 2003 2003 2003 1993 2003 2003 1993
n 5 4 6 6 6 6 6 6 6 5 6 1 5
Median 8.1 232 39.7 1.8 0.9 2 0.5 2.8 10.2 0.4 7.6 153.5 45.5
25% 8 228.5 39.4 1.7 0.9 1.9 0.4 2.4 10.1 0.3 7.5 153.5 36.8
BW- 75% 8.1 235.5 40.3 1.9 0.9 2.1 0.5 7.9 10.4 0.5 8.2 153.5 54.2
9 Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992-
Range 1993 1993 2003 2003 2003 2003 2003 2003 2003 1993 2003 2003 1993
n 2 2 3 3 3 3 3 3 3 2 3 1 2
Median 122.5 20.9 1.2 0.5 1.1 0.3 1.3 6.5 1.9 5.5 86.5 58.3 15.4
25% 113.1 17.2 1 0.5 1.1 0.2 0.5 5.5 1.7 5.3 74.6 58.3 14
BW- 75% 137.5 28.6 1.7 0.7 1.2 0.3 1.5 7.2 2.1 7 104.5 58.3 21.5
10 Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1999-
Range 1999 2003 2003 2003 2003 2003 2003 2003 2000 2003 2003 1992 2003
n 11 8 8 8 7 15 12 15 15 8 7 1 7
SC Ca2+ Mg2+ Na+ K+ Cl- NO3- SO42- DOC Si Sr Al Ba Fe Mn
Station pH (μS (mg (mg (mg (mg (mg (mg (mg (mg (mg (μg (μg (μg (μg (μg
cm-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1)
Median 97.5 15.5 0.7 0.6 0.8 0.3 1.5 5.5 0.4 5.8 72.6 51.3 15.9
25% 80.5 13.3 0.6 0.5 0.7 0.3 0.6 5.1 0.4 5.6 65.9 51.3 10.7
BW- 75% 111.5 20 1 0.7 1.2 0.3 3.4 6.3 0.5 6.9 78.8 51.3 27.6
11 Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1992- 1999- 1999-
Range 1999 2003 2003 2003 2003 2003 2003 2003 2000 2003 2003 1992 2003
n 6 8 8 8 7 15 15 15 12 8 7 1 7
Median 145 22.4 1.1 0.7 1.2 0.3 1.6 9 0.4 6.3 98.6 99.6 25.6
25% 141.9 22.3 1 0.6 1.2 0.3 0.4 8.6 0.4 6 96.6 99.6 16.8
BW- 75% 149.4 24.8 1.1 0.7 1.4 0.4 2.5 9.2 0.6 6.4 100.8 99.6 34
12 Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999-
Range 1999 2003 2003 2003 2003 2003 2003 2003 2000 2003 2003 1992 1999
n 11 7 7 7 7 15 15 15 19 7 6 1 5
Median 93.9 11.2 0.6 0.5 0.7 0.3 1.7 5.4 0.7 5.6 55.4 13.7
25% 80.7 9.4 0.5 0.4 0.7 0.2 0.5 5.2 0.5 5 41.5 10.8
BW- 75% 107.8 12.3 0.6 0.5 0.8 0.4 2.1 5.7 0.9 5.8 63.3 20.9
14 Date 1996- 1999- 1999- 1999- 1999- 1999- 1999- 1996- 1999- 1999- 1999-
51
Range 1999 2000 2000 2000 1999 2000 2000 2000 2000 2000 2000 2000
n 10 6 6 6 5 12 12 12 18 6 6 6
Median 7.8 108 16.4 0.9 0.6 1.1 0.3 1.7 6.6 0.5 5.5 68.8 87.8 12.6
25% 7.8 88 13.6 0.6 0.5 0.9 0.2 0.9 5.7 0.4 5.3 62.6 84 8.7
BW- 75% 7.8 132 21.1 1.2 0.7 1.4 0.3 3 7.1 0.5 6.6 89.4 91.7 15.6
15 Date 1993- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1992- 1999-
Range 1993 1999 2004 2004 2004 2004 2004 2004 2004 2003 2004 2004 1993 2004
n 1 7 13 13 13 11 19 19 19 17 13 11 2 8
Median 112 18.7 0.7 0.4 0.8 0.3 0.8 5.1 0.4 4.3 103.8 11.8
25% 104 16.2 0.6 0.3 0.8 0.2 0.6 4.9 0.4 4 90.8 10.7
BW- 75% 125 20.6 3.9 0.5 0.8 0.4 2.3 6 0.5 6 111 15.8
16 Date 1999- 1999- 1999- 1999- 1999- 1999- 1999- 1999- 1999- 1999-
Range 1999 2000 2000 2000 1999 2000 2000 2000 2000 2000 2000 2000
n 5 7 7 7 5 13 12 13 12 7 7 5
SC Ca2+ Mg2+ Na+ K+ Cl- NO3- SO42- DOC Si Sr Al Ba Fe Mn
Station pH (μS (mg (mg (mg (mg (mg (mg (mg (mg (mg (μg (μg (μg (μg (μg
cm-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1)
Median 198 32.5 1.1 0.7 1.9 0.3 0.5 9.9 0.9 7 132.5
25% 198 32.5 1.1 0.7 1.9 0.3 0.5 9.9 0.9 7 132.5
BW- 75% 198 32.5 1.1 0.7 1.9 0.3 0.5 9.9 0.9 7 132.5
17 Date
Range 1992 1992 1992 1992 1992 1992 1992 1992 1992 1992 1992
n 1 1 1 1 1 1 1 1 1 1 1
Median 133.5 27.6 1.1 0.6 1.3 0.4 0.9 8.3 0.5 5 120 62.6 9.2
25% 130.3 20.4 0.8 0.5 1 0.3 0.2 7.8 0.4 4.9 101.5 62.6 6.5
BW- 75% 157 29.8 1.2 0.7 1.5 0.4 10.5 9 0.5 6.4 139.6 62.6 16.1
18 Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1999-
Range 1999 2004 2004 2004 2004 2004 2004 2004 2003 2004 2004 1992 2000
n 6 14 14 14 12 19 17 19 15 14 13 1 6
Median 146 25.3 1 0.6 1.2 0.3 0.6 7.6 0.5 5.2 112.4 61.5 4
25% 141.5 23.3 0.8 0.5 1 0.3 0.2 7.1 0.5 5.1 105.6 61.5 2
BW- 75% 193 28.4 1.2 0.7 1.3 0.4 9.1 7.8 0.6 5.7 127.5 61.5 8.7
19 Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1999-
52
Range 1999 2004 2004 2004 2004 2004 2004 2004 2003 2004 2004 1992 2004
n 7 15 15 15 14 21 21 21 16 15 14 1 9
Median 143.7 24.4 1.2 0.6 1 0.4 1 7.2 0.5 5 102.8 67.5 6.2
25% 138.8 21.9 1 0.6 0.9 0.4 0.5 6.2 0.4 4.8 98.1 67.5 2.7
BW- 75% 147.1 28.8 1.4 0.7 1.4 0.4 1.5 8.3 0.6 5.9 113.1 67.5 12.5
21 Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1999-
Range 1999 2003 2003 2003 2003 2003 2003 2003 2003 2003 2003 1992 2003
n 11 11 11 11 9 17 17 17 13 11 10 1 8
Median 135 22.7 0.8 0.5 1.2 0.4 0.1 8 0.5 5.2 90.5 8
25% 128.5 21 0.7 0.5 1.1 0.3 0 7.7 0.5 4.9 88.1 8
BW- 75% 168 26.9 0.9 0.6 1.3 0.4 0.4 8.5 0.8 5.5 109.7 8
22 Date 1996- 1999- 1999- 1999- 1999- 1999- 1999- 1999- 1996- 1999- 1999-
Range 1999 2003 2003 2003 2003 2003 2003 2003 2000 2003 2003 1999
n 9 5 5 5 5 13 13 13 16 5 5 1
SC Ca2+ Mg2+ Na+ K+ Cl- NO3- SO42- DOC Si Sr Al Ba Fe Mn
Station pH (μS (mg (mg (mg (mg (mg (mg (mg (mg (mg (μg (μg (μg (μg (μg
cm-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1)
Median 180 25.4 1 0.6 0.9 0.4 7.2 8.7 1.1 5.3 84 61.3
25% 180 20.9 0.8 0.5 0.6 0.3 4.6 7.1 0.9 4.7 74.2 61.3
BW- 75% 180 26.4 1.1 0.6 1.2 0.4 9.2 9.2 1.4 5.5 93.9 61.3
23 Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992-
Range 1992 2003 2003 2003 2003 2003 2003 2003 2003 2003 2003 1992
n 1 3 3 3 3 3 3 3 2 3 2 1
Median 144 24.1 1.4 0.6 0.9 0.3 0.4 7 0.8 6.1 38
25% 144 24.1 1.4 0.6 0.9 0.3 0.4 7 0.8 6.1 38
BW- 75% 144 24.1 1.4 0.6 0.9 0.3 0.4 7 0.8 6.1 38
24 Date
Range 1992 1992 1992 1992 1992 1992 1992 1992 1992 1992 1992
n 1 1 1 1 1 1 1 1 1 1 1
Median 7.46 7.5 163.9 26.4 1.2 0.8 1.8 0.5 4.2 9.1 0.5 7.3
25% 7.46 7.5 159.3 25 1.2 0.8 1.8 0.4 3 8.9 0.5 7.2
BW- 75% 7.46 7.5 168.8 27.8 1.2 0.9 1.8 0.5 5.4 9.2 0.5 7.5
25 Date 1992-
53
Range 1992 1996 1992 1992 1992 1992 1992 1992 1992 1992 1992 1992
n 1 7 2 2 2 2 2 2 2 2 2 2
Median 146.5 23.7 1.7 0.6 1.3 0.3 0.8 8.4 0.6 6.4 97.3 56.7 13.8
25% 101.1 17.9 1.1 0.6 1.1 0.3 0.5 7.2 0.5 5.4 73.3 56.7 11
BW- 75% 194.4 32.7 2 0.8 2.1 0.4 1.9 10 0.9 6.8 141.2 56.7 18.2
26 Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1999-
Range 1999 2004 2004 2004 2004 2004 2004 2004 2004 2004 2004 1992 2000
n 10 12 12 12 9 18 18 18 16 12 11 1 7
Median 7.7 91.9 17.3 1.2 1.1 0.8 0.3 0.1 5.9 0.8 11 31 90.9 37.6
25% 7.6 54.6 9.4 0.7 0.9 0.7 0.3 0 4.2 0.5 9.3 23.8 66.9 32.8
BW- 75% 7.8 107.3 27.4 1.6 1.9 0.9 0.3 0.3 6.5 1.5 17.2 51.6 150 57.3
27 Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1992- 1999-
Range 1993 1999 2000 2000 2000 1999 2000 2000 2000 2000 2000 2000 1993 2000
n 4 18 22 22 22 10 18 13 18 26 22 7 5 7
SC Ca2+ Mg2+ Na+ K+ Cl- NO3- SO42- DOC Si Sr Al Ba Fe Mn
Station pH (μS (mg (mg (mg (mg (mg (mg (mg (mg (mg (μg (μg (μg (μg (μg
cm-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1)
Median 7.8 128.1 21.5 1.1 1.1 1 0.3 0.2 6.6 0.4 9.7 95.7 7.9 4.3 6.2 0.4
25% 7.5 121.7 20.5 1.1 1 0.9 0.3 0.1 6 0.3 9.1 86.8 3.4 3.8 1.6 0.2
BW- 75% 7.9 134.1 22.6 1.2 1.2 1.1 0.3 8.1 7.1 0.5 10.3 101.2 42.4 5.3 22.4 0.6
27A Date 1992- 1993- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1992- 2008- 1999- 2008-
Range 2013 2000 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013
n 40 28 49 57 57 55 61 52 61 56 55 39 29 20 27 20
Median 92.9 22.3 0.7 0.5 1 0.4 9.2 7.5 1.6 5.5 83 45.3 1.4
25% 82.4 17.5 0.6 0.5 0.7 0.3 8.7 7.1 1.5 4.6 67.1 45.3 1.4
BW- 75% 99.3 26.9 0.9 0.6 1.4 0.4 11.5 9 1.7 6.5 118 45.3 1.4
29 Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 2003-
Range 1996 2004 2004 2004 2004 2004 2004 2004 2003 2004 2004 1992 2003
n 6 6 6 6 6 6 6 6 2 6 5 1 1
Median 207 28.3 0.8 0.6 1.6 0.4 11.4 9.3 0.5 6 108 46.7 127
25% 207 26.5 0.8 0.6 1.5 0.4 11.2 8.8 0.4 5.7 104.8 46.7 64.9
BW- 75% 207 29.3 0.8 0.7 1.9 0.4 12 9.7 0.7 6.6 112.3 46.7 189
32 Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 2002- 2003-
54
Range 1992 2004 2004 2004 2004 2004 2004 2004 2003 2004 2004 1992 2004
n 1 7 7 7 7 7 7 7 3 7 6 1 2
Median 125.6 22.4 0.9 0.5 0.9 0.3 0.1 6.1 0.5 5.2 92.7 8.4
25% 107.8 18.4 0.6 0.4 0.7 0.2 0.1 5.6 0.4 4.7 78.8 4.4
BW- 75% 139.4 39.4 1.4 1 1 0.4 0.3 6.4 0.7 9.2 104.6 11
34 Date 1996- 1996- 1996- 1996- 1999- 1999- 1999- 1999- 1996- 1996- 1999- 1999-
Range 1999 2004 2004 2004 2004 2004 2004 2004 2003 2004 2004 2003
n 16 22 22 22 12 21 19 21 25 22 14 8
Median 6.7 29.8 3.7 0.6 0.8 0.3 0.3 0.1 4.8 0.7 8.5 23 10.8 1.3 22.8 2.2
25% 6.6 26.3 2.5 0.5 0.7 0.3 0.2 0.1 4.4 0.6 6.8 17.8 8.3 1 6.2 1.5
BW- 75% 6.8 39 4.9 1 1.1 0.3 0.3 6 5.9 1 13.5 32.1 11.7 1.6 57.1 2.6
36 Date 2008- 1996- 1996- 1996- 1996- 1999- 1999- 1999- 1999- 1996- 1996- 1999- 2008- 2008- 1999- 2008-
Range 2012 1999 2012 2012 2012 2012 2012 2011 2012 2012 2012 2012 2012 2012 2012 2012
n 8 13 29 29 29 18 25 16 25 30 29 21 9 9 19 9
SC Ca2+ Mg2+ Na+ K+ Cl- NO3- SO42- DOC Si Sr Al Ba Fe Mn
Station pH (μS (mg (mg (mg (mg (mg (mg (mg (mg (mg (μg (μg (μg (μg (μg
cm-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1)
Median 7.5 146 28.1 1.2 0.8 1.4 0.3 0.6 7.3 0.4 6.7 128.7 1.7 9.4 5 0.2
25% 7.4 119.7 26.9 1.2 0.7 1.2 0.3 0.4 6.6 0.4 6.4 117.3 1.1 9 1.2 0.1
BW- 75% 7.8 154.7 29.3 1.3 0.8 1.5 0.4 10.5 8.2 0.5 7.1 137.4 4.4 10.8 11.4 0.2
37 Date 1998- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1992- 1999- 1992- 2008- 1999- 2008-
Range 2013 2000 2013 2013 2013 2013 2013 2013 2013 2013 2013 2013 2012 2013 2013 2013
n 21 19 41 41 41 40 50 48 50 41 41 40 16 20 28 19
Median 132 21.1 1.4 0.8 0.9 0.3 0.1 6 0.4 8 88.8 83.6 29.3
25% 125.8 16.7 1.1 0.7 0.8 0.2 0.1 5.8 0.3 7.8 77.7 83.6 25
BW- 75% 139.8 24.3 1.6 0.9 1.1 0.3 1.2 6.5 0.5 8.7 101 83.6 35.5
38 Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1999-
Range 1999 2003 2003 2003 2003 2003 2003 2003 2003 2003 2003 1992 2003
n 6 12 12 12 10 18 18 18 15 12 11 1 8
Median 7.8 84.1 15.6 1 1.2 0.6 0.3 0.1 5.2 0.4 12 82.9 9.4 3.2 4.8
25% 7.4 80.1 14.5 0.9 1.1 0.6 0.3 0.1 4.6 0.3 11.5 77.5 5.5 3 1
BW- 75% 8 88.2 17 1.1 1.4 0.7 0.4 0.5 6.1 0.5 12.8 87.7 32 3.5 19.5
39 Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1992- 1999-
55
Range 1993 1999 2000 2000 2000 1999 2000 2000 2000 2000 2000 2000 1993 2000
n 2 7 10 10 10 8 15 7 15 15 10 7 3 7
Median 7.8 84.1 15.6 1 1.2 0.6 0.3 0.1 5.2 0.4 12 82.9 9.4 3.2 4.8 0.2
25% 7.4 80.1 14.5 0.9 1.1 0.6 0.3 0.1 4.6 0.3 11.5 77.5 5.5 3 1 0.1
BW- 75% 8 88.2 17 1.1 1.4 0.7 0.4 0.5 6.1 0.5 12.8 87.7 32 3.5 19.5 0.4
39A Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1992- 2008- 1999- 2008-
Range 2016 2002 2015 2015 2015 2015 2017 2017 2017 2017 2015 2015 2015 2015 2015 2015
n 126 69 149 158 159 147 176 145 176 166 157 122 105 87 98 79
Median 7.7 106.4 20.5 1.6 1.3 0.8 0.3 0.2 5.3 0.4 14.1 66 48.6 10.6
25% 7.6 98.5 15 1.1 0.9 0.7 0.3 0.2 4.8 0.3 10 64.8 25.2 6
BW- 75% 7.9 120.9 26.6 2.4 1.9 1 0.4 0.3 5.5 0.7 19.2 68.9 79.6 13.9
40 Date 1993- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1993- 1992- 1999- 1992- 1999-
Range 1993 1999 2000 2000 2000 1999 2000 2000 2000 2000 2000 2000 1993 2000
n 7 23 25 25 25 13 22 22 22 29 25 7 7 7
SC Ca2+ Mg2+ Na+ K+ Cl- NO3- SO42- DOC Si Sr Al Ba Fe Mn
Station pH (μS (mg (mg (mg (mg (mg (mg (mg (mg (mg (μg (μg (μg (μg (μg
cm-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1)
Median 7.5 97.3 16.5 1 1.1 0.6 0.3 0.2 5.6 0.5 10.1 57.6 2.3 3.3 2.7 0.1
25% 7.1 84.2 10.2 0.8 1 0.4 0.3 0.1 4.7 0.4 9.2 50.5 1.1 3 1.9 0.1
BW- 75% 7.9 119 25.8 1.2 1.2 0.9 0.3 7.2 5.9 0.7 11.4 79.2 31.8 4.4 6.4 0.1
40A Date 1992- 1993- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1992- 2008- 1999- 2008-
Range 2013 2000 2013 2013 2013 2013 2013 2011 2013 2013 2013 2013 2013 2013 2012 2013
n 36 34 59 67 67 55 60 50 61 67 65 39 29 20 25 19
Median 7.9 161.3 28.4 1.6 1.2 2 0.3 0.5 9.4 0.4 9.2 124.3 5.5 5.5 8.2 0.1
25% 7.7 156 27 1.5 1.1 1.8 0.3 0.3 8.1 0.3 8.5 114.8 2.9 4.9 2 0
75% 8 165.3 29.6 1.7 1.3 2 0.4 6.2 9.7 0.6 10 131.7 21.7 7.9 18.8 0.1
T-1
Date 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1992- 2008- 1999- 2008-
Range 2012 2000 2012 2012 2012 2012 2012 2012 2012 2012 2012 2012 2012 2012 2012 2012
n 35 33 56 65 65 54 62 62 62 64 63 40 27 19 28 16
Median 7.9 172 31.4 1.5 1.1 1.4 0.3 0.4 7.5 0.5 9.3 123.1 4.5 5.9 3 0.1
25% 7.7 167.7 29.3 1.5 1 1.3 0.3 0.3 6.2 0.4 8.3 117.7 2.2 5.3 1 0.1
75% 8 176.9 34.6 1.7 1.2 1.6 0.4 1 8 0.7 10.2 129.9 20.2 8.8 14.4 0.1
T-2
Date 1993- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1992- 2008- 1999- 2008-
56
Range 2012 2000 2012 2012 2012 2012 2012 2012 2012 2012 2012 2012 2012 2012 2012 2012
n 33 35 52 61 61 49 57 57 58 63 59 35 25 18 25 17
Median 7.6 164.3 27.7 1.3 0.8 1.1 0.3 0.3 5.7 0.8 7.9 113 7.1 7.4 7.2 0.7
25% 7.2 150 23.4 1.2 0.7 1 0.3 0.2 5.2 0.6 7.2 101 3.1 6.4 4.2 0.3
75% 7.9 172.6 31.6 1.6 1 1.3 0.4 0.6 7.2 1.1 9.2 127.1 23.1 8.9 14.8 1.2
T-3
Date 1993- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1999- 1992- 2008- 1999- 2008-
Range 2016 2000 2015 2015 2015 2015 2016 2016 2016 2016 2015 2015 2015 2015 2015 2015
n 52 37 65 73 74 64 78 77 79 84 72 47 34 28 36 28
Median 8.1 202 29.9 3.1 1.7 2.9 0.3 21.7 16.8 1.3 7.9 257 62.3 2
25% 8 196.5 29.3 3.1 1.6 2.9 0.3 21 16.2 0.8 7.7 253 51 1.4
75% 8.2 208 30.5 3.2 1.8 3.1 0.3 22 17 1.6 8.1 260.5 117.6 3
T-4
Date 1993- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 2003- 1992- 2003-
Range 1993 1996 2004 2004 2004 2004 2004 2004 2004 2003 2004 2004 1993 2004
n 4 9 30 30 30 30 30 29 30 6 30 25 5 15
SC Ca2+ Mg2+ Na+ K+ Cl- NO3- SO42- DOC Si Sr Al Ba Fe Mn
Station pH (μS (mg (mg (mg (mg (mg (mg (mg (mg (mg (μg (μg (μg (μg (μg
cm-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1) L-1)
Median 8.1 197.7 30.6 3.2 2 2.3 0.4 11.1 15.3 0.6 8.7 203.5 50.6 23.6
25% 8 191.3 27.9 3.1 1.9 2.1 0.3 0.4 13.8 0.3 8 195.3 47.3 23.6
T- 75% 8.2 203 31.6 3.2 2.1 2.5 0.4 21.3 16.4 0.9 9.6 205.8 51.9 23.6
5 Date 1993- 1992- 1992- 1992- 1992- 1992- 1992- 1993- 1992- 1992- 1992- 2003- 1992-
Range 1993 1996 2004 2004 2004 2004 2004 2004 2004 2003 2004 2004 1993 2003
n 4 9 8 7 8 8 8 6 8 6 7 3 4 1
Median 8 143.1 13.6 0.8 0.5 0.7 0.3 7.3 6.3 2.1 4.3 57.5 43 17.5
25% 7.8 135.6 12.2 0.6 0.4 0.6 0.3 4.5 5.7 1.8 4 55.2 31.7 14
T- 75% 8 214 21.2 1.8 1.4 1.1 0.4 8.1 8.9 2.8 5 65.1 70.4 28.1
6 Date 1993- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 1992- 2003- 1992- 2003-
Range 1993 1996 2004 2004 2004 2004 2004 2004 2004 2003 2004 2004 1993 2004
n 8 13 33 33 33 33 33 33 33 9 33 24 8 24
57
Supplemental Table 2.A.2. Annual total precipitation and average air temperature for the
water years 1992-2015.
Supplemental Table 2.A.3. Median solute concentrations in stream water during the four
periods during water years 1992-2015.
Rising Falling
Species
Dormant limb limb Growing
SC (μs) 110.2 89.3 83.3 130.3
Ca2+ (mg L-1) 22.6 17.1 16.8 23.0
Mg2+ (mg L-1) 1.1 0.8 0.8 1.0
Na+ (mg L-1) 0.7 0.5 0.5 0.7
K+ (mg L-1) 1.1 0.9 0.9 1.2
Cl- (mg L-1) 0.3 0.3 0.3 0.3
NO3- (mg L-1) 0.6 1.6 0.8 0.6
SO42- (mg L-1) 6.8 6.1 5.7 5.8
DOC (mg L-1) 1.3 1.9 1.7 2.0
Si (mg L-1) 5.3 4.8 5.0 5.4
Sr (μg L-1) 91.4 63.7 64.2 89.7
Al (μg L-1) 13.4 13.7 12.7 16.4
Ba (μg L-1) 7.8 5.5 7.2 8.5
Fe (μg L-1) 7.9 6.2 4.9 10.9
Mn (μg L-1) 2.1 0.8 0.8 1.3
58
Supplemental Table 2.A.4. Median solute concentrations at the 5th and 95th quantiles of
streamflow from 1992 to 2015
5th
Species
Quartile 95th Quartile
SC (μs) 149.2 82.1
Ca2+ (mg L-1) 31.2 15.0
Mg2+ (mg L-1) 1.5 0.7
Na+ (mg L-1) 1.0 0.5
K+ (mg L-1) 1.6 0.9
Cl- (mg L-1) 0.3 0.3
NO3- (mg L-1) 0.5 0.9
SO42- (mg L-1) 9.1 4.3
DOC (mg L-1) 1.3 3.6
Si (mg L-1) 5.8 4.5
Sr (μg L-1) 131.1 62.9
Al (μg L-1) 56.9 21.1
Ba (μg L-1) 11.1 6.0
Fe (μg L-1) 8.7 17.1
Mn (μg L-1) 22.0 1.3
59
Supplemental Table 2.A.5. Median solute concentrations in soil lysimeters and groundwater wells deeper that 1 meter
growing 49.1 3.3 0.3 0.6 0.2 0.2 0.1 3.6 0.8 6.4 30.2 5.7 2.6 1.6 0.3
dormant 114.1 19.3 1.1 1.1 0.8 0.3 0.3 5.9 0.5 10.2 87.2 5.4 4.0 2.5 0.1
rising limb 121.8 24.4 1.4 1.1 1.0 0.3 0.4 6.6 0.4 9.8 94.0 66.9 3.0 11.5 0.1
WELLS
falling limb 121.0 19.4 1.0 0.9 0.8 0.3 0.3 6.3 0.5 8.1 85.8 8.0 4.8 15.5 0.2
growing 144.0 23.6 1.2 1.1 1.2 0.3 1.4 7.0 0.5 9.6 100.0 17.9 3.6 3.6 0.1
Supplemental Table 2.A.6. Median solute concentrations across soil lysimeters and
groundwater well types
61
2.11 APPENDIX 2.B SUPPLEMENTAL FIGURES
62
Supplemental Figure 2.B.1. Timeseries of precipitation, air temperature, selected solute (Ca2+, Cl-, NO3-, DOC) concentrations in
streamflow, and stream discharge from water years: 1992-2015. The shaded grey areas indicate the snowmelt period.
63
Supplemental Figure 2.B.2 Solute concentrations in streamflow during the rising (blue) and falling (green) limbs of the snowmelt
period using data from all water years (1992-2015), showing a mild hysteresis for NO3- and a slight hysteresis for DOC.
3 Chapter 3: Synoptic Diagnosis of Catchment Hydrologic
and Biogeochemical Processes and Climate and Land Use
Change Effects in the Great Lakes Basin
3.1 Abstract
The goal of this study was to investigate the hydrologic and biogeochemical cycles
present in 138 catchments to better understand how land-use, snow abundance, and
surficial geology control stream chemistry using C-Q relationships (C=aQb) throughout
the US portion of the Great Lakes Basin (GLB). Multi-decadal C-Q relationships for
specific conductance and 15 solutes (Ca2+, Mg2+, Na+, K+, Cl-, SO42-, P, NH4+, NO3-, NO2-
Al, Fe, Mn, Si, and Sr) were created using data from downloaded from the USGS website
and the Water Quality Portal. Specific conductance (SC) was measured daily over varying
periods in the past 5 decades by the US Geological Survey (USGS) at 98 locations while
data for the other 15 solutes was collected less frequently, ranging from 1 to 11 samples
per year over multiple years to decades at 67 different catchments. SC, Ca2+, Mg2+, Na+,
K+, Cl-, SO42-, Si, and Sr concentrations often diluted (b < -0.1) during high flows,
becoming the dominant pattern for 46% - 97% of the sampled sites measuring each solute.
P, NH4+, NO3-, NO2- Al, and Fe concentrations often exhibited a flushing response (b >
0.1), increasing with flow for 70% - 100% of the sampled sites. The enrichment patterns
indicate that those solutes likely have higher concentrations in shallow flowpaths that are
activated during high flow conditions while dilution patterns indicate that those have higher
concentrations in deeper flowpaths. No measured solute was predominately chemostatic
throughout the studied catchments and solutes were instead frequently chemodynamic.
Land use and precipitation characteristics were frequently correlated to the strength (R2)
and slope (b-coefficient) of the C-Q relationships. These results indicate that land cover
and precipitation patterns are potential factors impacting the intensity and the pattern of
solute export by rerouting dominant flowpaths and changing biogeochemical processes.
Primarily, as urban areas expand and the precipitation regimes shift away from snow
dominated systems to rain dominated systems, the C-Q power-law relationship may
weaken, and the buffering capacity of nutrients and other contaminants may be reduced.
64
3.2 INTRODUCTION
Shifting precipitation patterns resulting from climate change coupled with rapid
land-use changes highlight the importance of water protection and sustainable use,
requiring that we understand processes controlling water quality and streamflow
generation. Natural hydrologic processes vary spatially and temporally, making each
watershed unique with differing hydrologic controls on inputs and outputs (Meybeck &
Moatar, 2012; Williams et al., 2018). This natural heterogeneity often makes it difficult to
apply findings between watersheds, and this heterogeneity can be further influenced by
human activities. Anthropogenic modifications to land use and atmospheric carbon dioxide
concentrations are quickly shifting hydrologic processes (flowpaths, stormflow,
infiltration, etc.) (Hung, James, Carbone, 2018; Rose, 2003) and climate variables
(precipitation, evapotranspiration, evaporation, storm event frequency and duration)
(Bieroza, Heathwaite, Bechmann, Kyllmar, Jordan, 2018; Williams et al., 2018), making
it even more critical to understand potential water quality vulnerabilities.
Under warmer conditions, northern latitudes see more rapid changes to snow
precipitation patterns than lower latitudes. Northern latitudes with will continue to see
earlier onset of snowmelt compared to historic trends (D’Orgeville, Peltier, Erler, & Gula,
2014; Demaria et al., 2016). Catchments that are currently snow-dominated areas will shift
to rain-dominated systems as temperatures rise and snow cover decreases (D’Orgeville,
Peltier, Erler, & Gula, 2014; Demaria et al., 2016). Snow cover extent has significantly
decreased over the last 100 years especially during the spring snowmelt season (Derksen
& Brown, 2012). Additionally, Demaria et al. (2016) project that future winter precipitation
will increase while snow water equivalent (SWE) will decrease in the northeastern and
upper Midwestern United States, including the Great Lakes Basin. Several models predict
that the number of days with snowpack and the snow depth will decrease while snowfall
intensity will increase (Demaria et al., 2016). Peak streamflows are predicted to increase
5% to 25% over historical patterns in the Great Lakes Region. The center of volume or the
hydrograph centroid, the day on which 50% of annual streamflow has occurred, is
estimated to happen 22 to 27 days earlier (Demaria et al., 2016).
65
In addition to changes in climate, anthropogenic changes in land-use through
expanding agriculture and urbanization can alter certain hydrologic, geomorphic, ecologic,
aquatic, and water quality features (Hung et al. 2018; Liu, Wang, & Li, 2013; Zambrano,
Pacheco-Muñoz, & Fernández, 2017). Urbanization increases impervious surfaces while
decreasing forest cover, reducing infiltration capacities and often increasing risks of
flooding (Zambrano et al., 2017). Changes in landscape characteristics coupled with
changing meteorological inputs increase runoff and severity and frequency of flood events,
subsequently impacting water quality (Liu et al., 2013; Zambrano et al., 2017).
Additionally, the expansive croplands in the Great Lakes area have extensive tile drainage
networks that change hydrological flowpaths while the accompanied fertilizer additions
change solute availability (Guo et al. 2017). These changes fuel rising water quality
concerns as the human population increases, and approximately 50% of the global
population is living in communities faced with complex water quality issues (Liu et al.,
2013).
66
temporal variance of C-Q relationships, with a particular focus on catchments within the
Great Lakes Basin.
The Great Lakes region includes land areas in the province of Ontario, Canada and
eight US states – Minnesota, Wisconsin, Michigan, Indiana, Illinois, Ohio, Pennsylvania,
and New York. However, since the same stream chemistry and landform data is not
collected in Canada, only the United States portion of the Great Lakes Basin is included in
this study. Water is an important resource within the Great Lakes region since the area
contains 84% of North America’s fresh water and 21% of the world’s surface fresh water,
and its population relies heavily on the Great Lakes for transportation, drinking water,
recreation, and power regeneration (Pathak, Kalra, & Ahmad, 2017). The 764,000 km2 area
has over 35 million inhabitants and contributes one fifth of the gross national product of
the US and one half of that of Canada (D’Orgeville et al., 2014; Quinn, 1981). The region
also contains a gradient of urban development from intensely urbanized areas around the
large urban centers of Detroit, Chicago, Cleveland, Green Bay, and Milwaukee to
catchments covered by cropland to those dominated by forests, creating an ideal landscape
to study C-Q relations. The area was glaciated resulting in varying types and depth of
surficial deposits with generally shallower deposits in the northern areas and deeper
deposits throughout the lower basin. There is also a gradient of snow inputs that is primarily
controlled by latitude and lake proximity (Eichenlaub, 1970). These variations in climate
and landscape create ideal conditions to study variation in C-Q relationships.
The objective of this study was to answer related questions by first investigating
how C-Q patterns respond with variation in snow abundance and urbanization amount as a
proxy of future urbanization. This objective was pursued by investigating multi-decadal C-
Q relationships throughout watersheds within the United States portion of Great Lakes
Basin using stream chemistry and discharge data from public stream chemistry datasets.
High frequency specific conductance and discharge measurements and low-frequency,
long-term water chemistry dataset for 15 additional solutes was used to create C-Q
relationships at decadal timescales. The questions posed are
Q1. What factors control C-Q relations in the Great Lakes Basin?
67
Q2. How might expected reductions in snowfall and increases in urbanization impact
future water quality and abundance?
Overall, about 35% of the land is managed for agriculture, 52% is forested and
often used for timber production, and 7% is urbanized with varying intensities, the highest
being around the larger urban centers of Detroit, MI, Milwaukee, WI, and Chicago, IL
(Wilcox et al., 2007). Due to the longer growing season and more favorable soil types,
more land is used for agriculture in the lower basin while most land in higher latitudes is
managed for timber production by private, commercial, or government stakeholders.
Annual average precipitation and air temperature varies throughout the four sub-
basins: Superior Basin (80.9 cm and 2.93°C), Michigan-Huron (84.5 cm and 6.49°C), Erie
(91.9 cm and 9.19°C), and Ontario (92.3 cm and 7.41°C) (Croley and Lewis, 2006). Winter
precipitation occurs as snow throughout the northern basin with snowfall amount and
timing varying with latitude, cyclone frequency, moisture, lake proximity, and altitude
(Eichenlaub, 1970). The mean seasonal snowfall within the Great Lakes Basin is on
average between 10% to 50% of the total regional winter precipitation (Eichenlaub, 1970;
Wright, Posselt, & Steiner, 2013).
68
3.4 METHODOLOGY
Most of the 138 sites were located within Michigan and Wisconsin (MI = 53; WI =
39) and 1 to 16 in the surrounding states (IN = 3, IN/IL=1, IN/MI=1, IN/OH=2, MI/OH=1,
MI/WI=1, MN = 5, NY = 15, OH = 16, PA = 1), which roughly reflected the contributing
areas composing the Great Lakes Basin (Table 3.A.1). The sampling frequency varied
between the solutes. SC had much higher sampling frequency with daily measurements
over 1 to 20 years with a median of 4 and a mean of 5 years. The other solutes had lower
sampling frequency that ranged from 1 to 356 samples per year with a mean of 8.9 samples
69
per year over 1 to 53 water years with a mean of 15 years. Mean stream discharge varied
from 0.029 m3s-1 to 5,832.1 m3 s-1 with a mean of 61.7 m3 s-1 and a standard deviation of
496.1 m3s-1. Watershed area ranged from 3.77 km2 to 163.85 km2 with a mean of 75.91
km2 and a standard deviation of 31.56 km2.
HUC 8 12-digit watershed boundaries from USGS were used for the watershed
boundaries for all 138 sites. Seven watersheds had 2 or more sampling locations. These
watersheds were delineated into smaller watersheds to contain only one gauge using the
ArcGIS “watershed tool” and 10m DEMs, using the USGS stream gauge locations as pour
points. The hydrologic interpretation of a pour point is that it is the outlet of the newly
created catchment. Those split watersheds were considered nested and all contributing area
to each USGS stream gauge was considered in following calculations (Table 3.A.2). If
discharge was not measured where the chemical concentrations were sampled (n=9), the
discharge measured at the nearest USGS stream gauge was used, which ranged between
approximately 3 m to 11.64 km from the water chemistry sampling location.
If there were multiple measurements in a day for a chemical species, only the first
observation was selected. Multiple daily entries were not used because there was only mean
daily stream discharge data available. This deleted 40,273 entries from 28,015 days out of
final 124,972. Most of the duplicates were analytical duplicates where 2 to 6 samples were
collected at the same time with a mean of 2. Only 6,498 days had entries that were deleted
because multiple samples were collected on the same day but not at the same time. The
number of samples deleted on days with multiple samples ranged from 1 to 14 samples per
day with a mean of 2.4 and a standard deviation of 1. A median of 4 entries were deleted
from C-Q relationships with multiple, non-analytical duplicates. The mean standard
deviations of the concentrations of these duplicates ranged from 7.0 x 10-4 to 6.25 x 102
with aluminum, iron, and manganese having the highest standard deviations. The mean
standard deviation of duplicates from the other 12 solutes (removing Al, Fe, and Mn) was
2.7 and ranged from 7.0 x10-4 to 13.7.
70
3.4.2 Concentration-Discharge Analysis
All concentration-discharge (C-Q) analysis was calculated on log10-log10 scales
across multi-annual timescales for all solutes using the equation: log(C) = log(a) +
b*log(Q) (Hofmeister et al., 2019; Zimmer et al., 2019), where C is the solute
concentration, Q is stream discharge, and a and b are constants with a as a coefficient with
units of concentration and b as an exponent (Hofmeister et al., 2019; Liu et al., 2017;
Zimmer et al., 2019). The b-coefficients were classified as chemostatic when concentration
varied little with changes in stream discharge (-0.1 ≤ b ≤ 0.1) or chemodynamic when the
concentration varied significantly with stream discharge (|b| > 0.1) (Godsey et al., 2009;
Herndon et al., 2015). Chemodynamic patterns were further classified into dilution when
b was negative and mobilization when b was positive (Figure 3.2). Whether the best-fit
slope was significantly different from a reference slope of zero was determined using
Student’s t-test. The strength of the C-Q relationship was measured by the coefficient of
determination (R²) and a relationship was considered significant if the p-value was less
than 0.05. The physical meaning of a-coefficient is unit-flow stream concentration.
Here 𝜇𝜇 is the mean and 𝜎𝜎 is the standard deviation of the solute concentration (C) and the
stream discharge (Q) (Zimmer et al., 2019). Chemostatic behavior is CVC/CVQ < 1 and
chemodynamic behavior is > 1 but the sign of the b-coefficient must be considered to
determine whether the chemodynamic pattern is flushing or dilution. CVC/CVQ >1
indicates there is greater variability in the solute concentrations and that discharge is not
the sole influencing factor, and solute concentrations could be influenced by other abiotic
and biotic factors such as in-stream processes, biotic processes, plant uptake,
71
anthropogenic additions, and colloidal transport or formation which greatly control
concentrations of nutrients and trace metals (Mulholland & Hill, 1997).
The b-coefficients and R2 values from the multi-annual C-Q relationships were
correlated against 9 independent variables: the watershed area, the weighted mean depth
of overlaying glacial deposits, the 30-year-average precipitation, the average yearly
accumulated SWE (in), the percent of total precipitation delivered as snow, the baseflow
index (BFI) at the point where discharge was measured, the percent of urbanized area
within the catchment, the percent of cropland area within the catchment, and the percent of
forested area within the catchment.
72
(2004-2019) and averaged. BFI describes the ratio of base flow to total stream flow as a
percent and was recorded from a geospatial layer that estimated the BFI throughout
ungauged streams in the US by extrapolating the BFI values calculated at USGS stream
gauges using an automated hydrograph separation computer program called the BFI
program (Wahl & Wahl, 1988; Wahl & Wahl, 1995) (water.usgs.gov/GIS/metadata/
usgswrd/XML/bfi48 grd.xml). Land-use data was derived from the 2016 version of the 30-
meter pixel National Land Cover Database 2016 (NLCD) (www.mrlc.gov/data), using all
intensities of developed area described by NLCD (Developed-Open Space, Developed-
Low intensity, Developed-Medium intensity, Developed-High intensity). The proportion
of cropland and forested areas was similarly calculated using pixels that were described as
cropland (“Cultivated Crops”) and forested (“Evergreen Forest”, “Deciduous Forest”, and
“Mixed Forest”).
3.5 RESULTS
73
Si displayed all three C-Q patterns, but no solute was persistently dominated by
chemostatic C-Q behavior (Figures 3.5 and 3.6). It is worth noting that 22% of catchments
(7 catchments) had a significant dilution pattern for NO3-. Most solutes that were
dominated by dilution had mean CVC/CVQ values that were < 1, while those exhibiting
flushing patterns were often ~1 or > 1 (Figure 3.7).
74
cropland was also positively correlated with the b-coefficients of NO3-, and Si (ρ = 0.4).
Forested area was negatively correlated with the b-coefficients from SC, Ca2+, Mg2+ (ρ = -
0.4 to -0.5) and positively correlated with P (ρ = 0.3).
3.6 DISCUSSION
By combining information from multiple public datasets, I analyzed stream
chemistry data from 138 diverse catchments across the GLB to identify which physical and
environmental catchment characteristics were correlated with specific C-Q attributes and
strengths. In the multi-annual C-Q relationships, solute behavior was predominantly
chemodynamic among the sites, as geogenic solutes were diluted and biogeochemically
reactive and anthropogenically-enriched solutes were often enriched with an increase in
streamflow. The opposite chemodynamic patterns between the geogenic and biogenic
solutes are likely caused by contrasting distributions of solutes from shallow soil water to
deeper groundwater (Chen et al., 2020; Herndon et al., 2015) and the mixing from two or
more endmembers with distinct solute concentrations (Evans & Davies, 1998). Results
from Spearman’s Rank-Order Correlation indicate that land-use and precipitation
characteristics, especially urbanization, agriculture, forestation, and snow abundance, are
frequently significantly correlated to the b-coefficients and R2 values from the C-Q power-
law relationships, emphasizing the influences of anthropogenic modifications on
streamflow generation processes and flowpath dynamics.
The combination of R2 values and b-coefficients from the C-Q power-law relations
reveals relevant information on many aspects of hydrologic and hydrochemical processes
at the catchment scale (Figure 3.2), as described in the method section. In the following
sections, a holistic view is applied to provide perspective on some hydrologic and
biogeochemical processes useful for managing water resources and evaluating climate and
land use change effects that can be inferred from C-Q relationships. Examples are given to
demonstrate how to use R2 values and b-coefficients to gain insights into hydrologic and
hydrochemical characteristics of a catchment. Limited by data availability, however, every
application cannot be mentioned in this study. The demonstration below includes (1) strong
75
and consistent dilution patterns of major ions in forested catchments, (2) disconnection of
hydrologic flowpaths accompanied by transport of pollutants in urban areas, and (3)
persistent flushing of nutrients, specially NO3-, in agriculture areas.
Additionally, major ions and nutrients had contrasting CVC/CVQ values with the
former having low values (<1) and the latter having higher mean values (~1 or >1) (Figure
3.7; Table 3.1). When interpreted, this pattern indicates that the variation in nutrient
concentrations is similar to the amount of variation in discharge while there is much less
76
variation in geogenic solute concentrations compared to variation in discharge, indicating
major ions are chemostatic or near-chemostatic while nutrients and trace elements are
chemodynamic or near chemodynamic. This interpretation is not always a good
determinate of C-Q behavior. For example, Ca2+ and Mg2+ have CVC/CVQ values that are
less than 1 even when the b-coefficients are relatively small (-0.2 < b < -0.6), but both have
high R2 values indicating a strong C-Q power-law relationship (Figure 3.7). Additionally,
Ca2+, Mg2+, and Na+ are often conservative, indicating that changes in concentrations (b-
coefficients) are a product of mixing between contributing endmembers even when b is
less than 1 (Liu et al., 2004; Liu et al., 2020; Sebestyen et al., 2008). Meanwhile, major
nutrients had lower median R2 values compared to major ions (major nutrients R2: 0.1 –
0.27; major ions R2: 0.36 – 0.52) and are not conservative but have other processes
controlling concentrations. Thus, although the CVC/CVQ values for nutrients where close
to 1, CVC/CVQ values are not an accurate indicator of how concentration and discharge are
coupled since the observed variation in concentrations could match that of discharge but
due to other active processes independent of discharge. However, comparing CVC/CVQ
values could be useful in identifying locations or moments where or when there is variable
solute export.
77
processes (ie, colloidal formation and transport, biotic uptake, and microbial
transformations), which deteriorate the C-Q relationship. Additionally, nutrients and trace
elements may have variable solute stores throughout the catchment or limited solute
sources that deplete or become hydrologically connected during the snowmelt event. If
changes in solute concentrations are uncoupled with changes in stream discharge, this
would deteriorate the C-Q relationships for those solutes. An example would be if a solute
begins a storm or melt event with enrichment but as the event progresses becomes depleted
and cannot match the initial response, changing export patterns mid-event and weakening
R2.
Snow abundance and baseflow index had similar effects on solute export. Both
were positively correlated with the b-coefficient of major ions and Si (Figure 3.8b).
Snowpack releases large amounts of meltwater during spring snowmelt and, as that new
water flows through the subsurface, it accumulates solutes from chemical weathering
reactions. This release of new water that is low in ionic concentrations strongly dilutes
concentrations of major ions and Si that are highest in groundwater, leading to more
negative b-coefficients for geogenic solutes in areas with high snow abundance. Under
future conditions when there is less snow precipitation but continued snow cover, one can
expect less dilution of major ions and subsequently higher b-coefficients since there is less
new water released during snowmelt. However, less snowmelt accompanied by freezing
conditions can create deeper ground frost which can impede infiltration and result in more
runoff (Shanley et al., 2002). Thus, there is uncertainty in how hydrologic processes and
solute export will respond to climate change with the predicted responses dependent on
winter temperatures.
78
3.6.3 Land-Use Impacts on Solute Transport and Hydrologic
Flowpaths
Unlike C-Q patterns often observed in urban and agriculture areas, C-Q patterns for
nutrients can be diluted in forested catchments as seen in a study on C-Q patterns of a
headwater stream in northern Michigan (Hofmeister et al., 2019). Hofmeister et al. (2019)
observed dilution C-Q patterns for NO3-, PO43-, and SO42- which they attributed to restricted
solute availability by biological activity and large inputs of snowmelt and precipitation
diluting concentrations. Thus, transport from uplands may be supply limited. During
events, the subsurface flow low in NO3- dilutes the higher NO3- stream concentrations,
creating the observed dilution pattern. The spatial distribution of NO3- b-coefficients
79
throughout the GLB also supports this idea of limited NO3- with frequent NO3--Q dilution
or chemostatic behavior in the northern basin where catchments are forested. However,
nearly 80% of the studied catchments measuring NO3- showed enrichment, indicating that
NO3- is no longer source limited in most areas in the GLB.
80
contrasting concentrations creating stronger C-Q relationships (Chen et al. 2020; Herndon
et al. 2015). Fertilizer application on cropland elevates nutrient solute concentrations in
surface soils, creating solute rich dominant shallow flowpaths and extenuating distinction
in solute concentration between endmembers. This distinction between endmember
concentrations simultaneously increase the C-Q R2 value and |b| coefficient when those
endmembers mix during events.
Catchments with more cropland also had stronger enrichment for NO3- with over
75% of the study sites exhibiting enrichment, but cropland is accompanied by more
chemostatic relations for SC, Ca2+, and Mg2= (Figure 3.5, Figure 3.8b; Figure 3.10). Rapid
flushing of NO3- rich shallow near-stream groundwater or soil water create stronger
enrichment behavior (Bowes et al., 2015; Chen, Wu, & Hong, 2012; Jaing et al, 2010;
Outram, Cooper, Sunnenberg, Hiscock, & Lovett, 2016). Flushing is enhanced when
groundwater tables rise into nutrient-rich soil layers and mobilize solutes into nearby
streams along natural flowpaths or subsurface tile drainage (Williams et al., 2018). Tile
drainage is a common agricultural practice in the Midwest with 20.6 million hectares of
land drained by surface and subsurface tile drainage networks in the states creating the
Great Lakes Basin (Fausey, Brown, Belcher, & Kanwar, 1995). Jiang, Frankenberger,
Bowling, & Sun (2014) studied nitrate export from catchments near Lake Erie dominated
by agriculture (including USGS site: 4193500). They found nitrate generally increased and
the enrichment was especially clear in areas with tile drainage, supporting the findings of
Williams et al. (2018) who found that larger drainage networks created higher solute
concentrations. Antecedent conditions also control the magnitude of NO3- export with more
export occuring when groundwater levels are already high and soil is saturated, allowing
for faster and more contributions from shallow soils (Outram et al., 2016). Tile drainage
can delay the response of the groundwater table and transport which might complicate the
hydrologic response of SC, Ca2+, Mg2+, and Si and create the observed chemostatic
behaviors. Nitrification of nitrogen fertilizers can create nitric acid (HNO3) which
accelerates chemical weathering and raise concentrations of SC, Ca2+, Mg2+, and Si in
81
shallow soils and decreasing the disparity between solute concentrations in shallow
flowpaths compared to deeper ones.
82
from mixing between groundwater and lake water and/or denitrification (Knights et al.,
2017). The three observed flushing patterns in the GLB could be a sign of the growing
issue of NO3- groundwater contamination that will become increasingly important if left
unaddressed and nitrogen fertilizers are continually applied.
Groundwater contamination is not the only future concern and high stream nutrient
loading is a continued threat despite recent decreases as loads are predicted to increase due
to increased precipitation and stream discharge under future warming scenarios
(Culbertson et al., 2016). If fertilizer application does not remain the same and instead
matches increased P plant uptake, then TP loading will increase under both climate change
scenarios and exacerbate eutrophication (Culbertson et al., 2016). These pressing issues
highlight the importance of studying C-Q patterns, especially since controlling factors are
changing so rapidly through land-use and climate changes.
3.7 CONCLUSION
The relationships between environmental and physical catchment characteristics and
solute export behavior were investigated using concentration-discharge analysis and
Spearman’s Rank Order Correlation Analysis. Results indicate that the long-term export
patterns of Mg, Na, Cl, SO42-, P, NH4+, Al, Fe, and Sr are relatively consistent across the
Great Lakes Basin despite changes in land-use and precipitation patterns. CVC/CVQ values
do not accurately summarize the concentration-discharge pattern and can under-represent
endmember mixing, resulting in misclassifying solute export as chemostatic rather than
chemodynamic. Among the studied environmental factors, snow abundance, potentially
through greater snowmelt, increased the dilution of solutes and strengthened the
concentration-discharge power-law relationship. Fertilizer additions and rapid subsurface
flow increased the b-coefficients from certain solutes and increased the power-law
relationship by exaggerating the gradient of solute concentrations throughout the soil
profile. Urbanization circumvents shallow flowpaths increasing some major ions b-
coefficients and deteriorating the concentration-discharge relationship by decoupling the
response in concentration to discharge through inconstant additions and flowpaths. The
83
observed impacts of cropland and urbanization on solute concentrations might change with
the proximity of the cropland to the sampling location, with more prominent effects when
the cropland and urban areas are close to the sampling location. This could be an interesting
research question for future studies and would provide useful information for land
managers prescribing best management practices. Lastly, forested areas have strong
dilution of major ions and less negative or positive b-coefficients for NO3- and P. These
results indicate that snow abundance and land-use are key factors influencing
concentration-discharge patterns and that certain patterns can be expected in those
environments, which permits characterization of solute export patterns across similar
catchments.
3.8 REFERENCES
Argerich A, Johnson S L, Sebestyen S D, Rhoades C C, Greathouse E, Knoepp J D, …&
Ice G G. (2013). Trends in stream nitrogen concentrations for forested reference
catchments across the USA. Environmental Research Letters, 8,
doi:10.1088/1748-9326/8/1/014039
Chen N, Wu J, & Hong H. (2012). Effect of storm events on riverine nitrogen dynamics
in a subtropical watershed, southeastern China. Science of the Total Environment,
431, pp. 357–365. doi:10.1016/j.scitotenv.2012.05.072
84
Chen X, Tague C L, Melack J M, & Keller A A. (2020). Sensitivity of nitrate
concentration-discharge patterns to soil nitrate distribution and drainage
properties in the vertical dimension. Hydrological Processes, 34, 2477–2493.
doi:10.1002/hyp.13742
Croley T E, & Lewis C F M. (2006). Warmer and drier climates that make terminal Great
Lakes. Journal of Great Lakes Research, 32, pp. 852–869.
Derksen C, & Brown R. (2012). Spring snow cover extent reductions in the 2008 – 2012
period exceeding climate model projections. Geophysical Research Letters, 39,
L19504. doi:10.1029/2012GL053387
Dyer J L, & Mote T L. (2006). Spatial variability and trends in observed snow depth over
North America. Geophysical Research Letters. 33, L16503,
doi:10.1029/2006GL027258
Eichenlaub V L. (1970). Lake effect snowfall to the lee of the Great Lakes: Its role in
Michigan. Bulletin of the American Meteorological Society, 51, pp. 403–412.
Elser J, & Bennett, E. (2011). Phosphorus cycle: A broken biogeochemical cycle. Nature,
478, 29–31, doi:10.1038/478029a
85
Godsey S E, Kirchner J W, & Clow D W. (2009). Concentration-discharge relationships
reflect chemostatic characteristics of US catchments. Hydrological Processes, 23,
pp. 1844–1864. doi:10.1002/hyp.7315
Hidalgo, H. G., and Coauthors. (2009). Detection and attribution of streamflow timing
changes to climate change in the western United States. Journal of Climate, 22,
pp. 3838–3855, doi:10.1175/ 2009JCLI2470.1.
86
Larson G, & Schaetzl R. (2001). Origin and evolution of the Great Lakes. Journal of
Great Lakes Research, 27, pp. 518–546. doi:10.1016/S0380-1330(01)70665-X
Liu F, Conklin M H, & Shaw G D. (2017). Insights into hydrologic and hydrochemical
processes based on concentration-discharge and end-member mixing analyses in
the mid-Merced River Basin, Sierra Nevada, California. Water Resources
Research, 53, pp. 832–850. doi:10.1002/2016WR019437
Liu F, Williams M W, & Caine N. (2004). Source waters and flow paths in an alpine
catchment, Colorado Front Range, United States. Water Resources Research, 40,
W09401. doi:10.1029/2004WR003076
Liu Z, Wang Y, Li Z, & Peng J. (2013). Impervious surface impact on water quality in
the process of rapid urbanization in Shenzhen, China. Environmental Earth
Sciences, 68, pp. 2365–2373. doi:10.1007/s12665-012-1918-2
Meybeck M, & Moatar F. (2012). Daily variability of river concentrations and fluxes:
Indicators based on the segmentation of the rating curve. Hydrological Processes,
26, pp. 1188–1207. doi:10.1002/hyp.8211
87
Pathak P, Kalra A, & Ahmad S. (2017). Temperature and precipitation changes in the
Midwestern United States: Implications for water management. International
Journal of Water Resources Development, 33(6), pp. 1003–1019.
doi:10.1080/07900627.2016.1238343
Quinn F H. (1981). Secular changes in annual and seasonal Great Lakes precipitation,
1854-1979, and their implications for the Great Lakes water resources studies.
Water Resources Research, 17(6), 1619–1624.
88
University of Milwaukee. (2005). Stormwater Monitoring Program 200-2004. (Report
No. M03002P17) Milwaukee Metropolitan Sewerage District (MMSD).
89
3.9 TABLES AND FIGURES
3.9.1 Tables
Table 3.1. Mean values of multi-annual C-Q relationships by solute
Mean
Mean R2
Species Sig/Total Mean b (sd) CVC/CVQ
(sd)
(sd)
SC (us) 95/98 0.44 (0.28) -0.15 (0.13) 0.21 (0.11)
Ca2+ (mg L-1) 55/65 0.36 (0.26) -0.15 (0.11) 0.18 (0.08)
Mg2+ (mg L-1) 59/65 0.53 (0.23) -0.19 (0.08) 0.19 (0.07)
Na+ (mg L-1) 48/59 0.52 (0.23) -0.28 (0.17) 0.37 (0.22)
K+ (mg L-1) 27/57 0.22 (0.14) -0.04 (0.18) 0.42 (0.66)
Cl- (mg L-1) 55/64 0.32 (0.20) -0.23 (0.26) 0.46 (0.57)
NO3- (mg L-1) 32/48 0.27 (0.19) 0.31 (0.54) 0.71 (0.63)
NH4+ (mg L-1) 15/40 0.10 (0.07) 0.21 (0.20) 1.09 (0.36)
NO2- (mg L-1) 24/57 0.18 (0.10) 0.15 (0.33) 0.91 (0.44)
P (mg L-1) 44/65 0.24 (0.20) 0.32 (0.24) 0.80 (0.34)
SO42- (mg L-1) 51/57 0.32 (0.22) -0.19 (0.11) 0.31 (0.17)
Si (mg L-1) 26/38 0.18 (0.12) 0.06 (0.34) 0.41 (0.26)
Sr (μg L-1) 36/42 0.59 (0.28) -0.27 (0.13) 0.41 (0.55)
Al (μg L-1) 27/45 0.40 (0.20) 0.54 (0.25) 1.23 (1.41)
Fe (μg L-1) 34/53 0.27 (0.22) 1.11 (2.92) 2.27 (3.63)
Mn (μg L-1) 24/53 0.19 (0.20) 0.18 (0.81) 1.47 (1.49)
90
3.9.2 Figures
Figure 3.1. Map showing study catchments throughout the US portion of the Great Lakes
Basin.
91
92
Figure 3.2. A flow chart used for diagnosing catchment hydrologic and biogeochemical characteristics/behaviors based on
Concentration-Discharge (C-Q) relationships. Note that some behaviors may exist at various power-law strengths (those crossing the
dashed boxes).
93
Figure 3.3. Distribution of catchment characteristics among the 138 study sites with the dashed vertical line indicating the median
value.
94
Figure 3.4. Frequency of R2 values from the significant (p < 0.05) multi-annual C-Q relationships created for 16 solutes. The vertical
dashed line indicates the median of the R2 values
Figure 3.5. Prevalence of dilution, chemostatic, and enrichment C-Q patterns among the
sites with significant (p < 0.05) multi-annual power-law relationships for 16 solutes.
95
96
Figure 3.6. Frequency of b-coefficients from the significant (p < 0.05) multi-annual C-Q relationships created for 16 solutes. The
vertical dashed line indicates the median while the shaded red area indicated the chemostatic range (b≤ ±0.1).
97
Figure 3.7. Relationship between C-Q |b|-coefficients and CVC/CVQ from the significant (p<0.05) multi-annual C-Q relationships for
16 analytes
98
Figure 3.8. 𝜌𝜌 from Spearman Rank Correlation Analysis for significant relationships (p<0.05) between the R2 and b-coefficients from
the C-Q relationships and catchment characteristics among SC and the 15 solutes. (a) illustrates the correlations between the
catchment characteristics and the R2 values from all significant C-Q relationships. (b) illustrates the correlations between the
catchment characteristics and the b-coefficients from all significant C-Q relationships. Blue indicates positive correlation while green
indicates negative correlation. Note that the results of solutes that do not have a significant relationship with any factors were not
shown.
99
Figure 3.9. Spatial distribution of SC-Q b-coefficients and R2 values with SWE (a) and land cover (b) on background.
100
Figure 3.10. Spatial distribution of NO3--Q b-coefficients and R2 values with SWE (a) and land cover (b) on background.
3.10 APPENDIX 3.A: SUPPLEMENTAL TABLES
Supplemental Table 3.A.1. Site locations for the 138 locations used for this study. *
indicates that the site measured water chemistry parameters and not specific conductance.
** indicates that the site measured specific conductance and additional water chemistry
parameters. Sites without a * are sites that only measures specific conductance and
stream discharge.
Elevation
Site No Site ID STATES Latitude Longitude
(ft)
4014500** USGS-04014500** MN 47.33742 -91.20072 755
4015475 USGS-04015475 MN 47.527146 -92.12267 1483
4024000** USGS-04024000** MN 46.70328 -92.4188 1158
4024090 USGS-04024090 MN 46.517444 -92.48214 984
4024093 USGS-04024093 MN 46.5155 -92.46269 928
4024430* USGS-04024430* WI 46.63333 -92.09389 659
4026561 USGS-04026561 WI 46.394722 -90.59 1266
4027000** USGS-04027000** WI 46.48661 -90.6963 735
4040000** USGS-04040000** MI 46.720774 -89.20709 856
4043097 USGS-04043097 MI 46.734655 -88.44319 702
4043126 USGS-04043126 MI 46.721043 -88.33013 1401
4043140 USGS-04043140 MI 46.751043 -88.3618 1004
4043150 USGS-04043150 MI 46.804097 -88.31707 797
4043238 USGS-04043238 MI 46.782222 -87.8775 1106
4043244 USGS-04043244 MI 46.785759 -87.85235 1033
4043275 USGS-04043275 MI 46.713539 -87.84069 1532
4044724* 21MICH-030077* MI 46.34078 -86.85015 853
4045500* USGS-04045500* MI 46.574583 -85.26964 748
4059000** USGS-04059000** MI 45.908917 -87.21353 801
4059500** USGS-04059500** MI 45.754967 -87.20208 702
4063700** USGS-04063700** WI 45.76361 -88.46361 1460
4067500** USGS-04067500** MI,WI 45.32583 -87.66333 666
4072150** USGS-04072150** WI 44.53339 -88.12969 659
SOKAOGON_WQX-
4074538* WI 45.4883 -88.96372 1572
SC2@55*
SOKAOGON_WQX-
4074548* WI 45.47941 -88.99789 1575
SC3@M*
4080798 USGS-04080798 WI 44.524444 -89.33778 1112
4081000 USGS-04081000 WI 44.329167 -88.99583 804
4085108 USGS-04085108 WI 44.371944 -88.09222 656
4085427** USGS-04085427** WI 44.10617 -87.71603 676
4085813 USGS-04085813 WI 43.63833 -87.89842 840
4085845 USGS-04085845 WI 43.696662 -87.82092 715
4086149 USGS-04086149 WI 43.527219 -88.22232 961
4086200 USGS-04086200 WI 43.550274 -88.18843 961
4086265 USGS-04086265 WI 43.425 -88.08093 886
4086340 USGS-04086340 WI 43.482776 -88.06093 853
4086360 USGS-04086360 WI 43.472777 -87.98981 860
4086500 USGS-04086500 WI 43.323056 -87.97861 853
4086600* WI_MMSD-RI-01S* WI 43.28028 -87.9425 725
101
WI_MMSD-RI-04S &
4087000** WI 43.10001 -87.90897 679
USGS-04087000**
4087030 USGS-04087030 WI 43.172778 -88.10389 781
4087050 USGS-04087050 WI 43.206674 -88.03842 758
4087070 USGS-04087070 WI 43.123611 -88.04361 738
4087088 USGS-04087088 WI 43.054722 -88.04611 741
4087119 USGS-04087119 WI 43.043833 -88.00511 696
4087120 USGS-04087120** WI 43.04556 -87.99972 696
WI_MMSD-RI-11M &
4087142** WI 43.0339 -87.93425 627
USGS-04087142**
4087170 USGS-04087170 WI 43.024444 -87.89833 627
WI_MMSD-OC-05S &
4087204** WI 42.925 -87.87 673
USGS-04087204**
4087214* WI_MMSD-RR-04S* WI 42.945 -88.01417 778
4087220* WI_MMSD-RR-05S* WI 42.87361 -87.99583 725
4087233 USGS-04087233 WI 42.815556 -87.99472 728
4087240 USGS-04087240 WI 42.751389 -87.82361 682
4087242 USGS-04087242 WI 42.733909 -87.78452 620
INSTOR_WQX-1899 &
4092750** IL,IN 41.649202 -87.46865 587
USGS-04092750**
INSTOR_WQX-1918 &
4095090** IN 41.622259 -87.17587 600
USGS-04095090**
4095300* INSTOR_WQX-1924* IN 41.716707 -86.85975 620
4099000* USGS-04099000* IN,MI 41.800883 -85.7561 784
4101500** USGS-04101500** MI 41.829214 -86.25973 702
4107850 USGS-04107850 MI 42.482256 -85.79836 728
4108500* USGS-04108500* MI 42.593641 -85.9842 650
21MICH_WQX-030077 &
4108660** MI 42.65086 -86.1067 705
USGS-04108660**
4116000* USGS-04116000* MI 42.971977 -85.06917 745
4119400 USGS-04119400 MI 43.024188 -86.02644 610
4121650 USGS-04121650 MI 43.693634 -85.46754 948
4121970 USGS-04121970 MI 43.434746 -85.66532 860
4122500* USGS-04122500* MI 43.945006 -86.27869 686
4126802 USGS-04126802 MI 44.903329 -85.96231 689
4126970 USGS-04126970 MI 44.65667 -85.43673 892
4135000* USGS-04135000* MI 45.093901 -83.49998 653
21MICH-350061 & USGS-
4137500** MI 44.4364 -83.43386 636
04137500**
4142000** USGS-04142000** MI 44.07252 -84.01999 741
4144032 USGS-04144032 MI 42.881698 -83.98468 833
4145000* USGS-04145000* MI 43.254748 -84.10553 600
4149000* USGS-04149000* MI 43.308359 -83.95358 591
4157005 USGS-04157005 MI 43.42197 -83.95192 597
4159492 USGS-04159492 MI 43.15086 -82.62465 728
4160398 USGS-04160398 MI 42.858642 -82.53797 623
4160625 USGS-04160625 MI 42.768366 -82.51214 584
4160800 USGS-04160800 MI 42.720031 -83.35355 994
4160900 USGS-04160900 MI 42.660309 -83.39022 961
4161000 USGS-04161000 MI 42.633366 -83.22438 856
4161540 USGS-04161540 MI 42.688366 -83.14299 823
4161820** USGS-04161820** MI 42.614478 -83.02659 643
4163030 USGS-04163030 MI 42.537813 -83.00576 610
102
4164000 USGS-04164000 MI 42.577812 -82.95159 600
4164980 USGS-04164980 MI 42.642811 -82.93326 594
4165500 USGS-04165500 MI 42.59587 -82.90881 600
4165557 USGS-04165557 MI 42.561424 -82.8452 574
21MICH-500233 & USGS-
4165559** MI 42.596423 -82.82603 581
04165559**
4166100 USGS-04166100 MI 42.447536 -83.29743 656
4166500 USGS-04166500 MI 42.373092 -83.25465 617
4167150 USGS-04167150 MI 42.330593 -83.24798 610
4168400 USGS-04168400 MI 42.308371 -83.25271 607
4170000 USGS-04170000 MI 42.578921 -83.62661 958
4174500 USGS-04174500 MI 42.286982 -83.73383 860
4174518 USGS-04174518 MI 42.264761 -83.68827 833
4176500** USGS-04176500** MI 41.9606 -83.53105 633
4178000* 21OHIO_WQX-510220* IN,OH 41.38561 -84.80163 830
4183038 USGS-04183038 IN 41.183083 -84.86969 738
4183500 USGS-04183500 IN,OH 41.198939 -84.7444 725
4184500* 21OHIO_WQX-500330* MI,OH 41.6595 -84.24911 728
4185000* 21OHIO_WQX-300020* OH 41.5045 -84.42967 712
4186500* USGS-04186500* OH 40.94866 -84.26606 732
4188496 USGS-04188496 OH 40.97922 -83.65327 801
4189000* 21OHIO_WQX-500040* OH 41.05589 -83.68799 778
4193500* USGS-04193500* OH 41.50005 -83.71271 640
4195820* 21OHIO_WQX-S02P08* OH 41.49116 -83.22465 594
4196500* 21OHIO_WQX-500860* OH 40.85061 -83.25631 843
4196800* 21OHIO_WQX-500850* OH 40.92283 -83.34881 804
4197100* 21OHIO_WQX-U03S02* OH 41.02228 -83.10964 833
4198000* USGS-04198000* OH 41.30783 -83.15881 656
21OHIO_WQX-501030
4199000** OH 41.30089 -82.60823 643
&USGS-04199000**
4199155* 21OHIO_WQX-K01G31* OH 41.34838 -82.51378 666
21OHIO_WQX-501260 &
4199500** OH 41.38199 -82.31683 712
USGS-04199500**
4200500 USGS-04200500 OH 41.380324 -82.10459 719
4208000** USGS-04208000** OH 41.39533 -81.62985 732
4212100 USGS-04212100 OH 41.718934 -81.22788 686
4213152 USGS-04213152 PA 42.073889 -80.23472 577
4213500** USGS-04213500** NY 42.46333 -78.93417 948
4214500* USGS-04214500* NY 42.85472 -78.755 620
4216000 21NYDECA-01031002 NY 42.87778 -78.91639 571
4217000* USGS-04217000* NY 42.9975 -78.18861 899
4226000 USGS-04226000 NY 42.682472 -77.82869 728
4231600 USGS-04231600 NY 43.141722 -77.61631 541
4232040* USGS-04232040* NY 43.05417 -77.49083 518
4232046* USGS-04232046* NY 43.10611 -77.46194 459
4232047* USGS-04232047* NY 43.12111 -77.47639 459
4232050* USGS-04232050* NY 43.12803 -77.527 446
4249000** USGS-04249000** NY 43.45167 -76.50528 302
4250750* USGS-04250750* NY 43.81333 -76.07472 545
4260500* USGS-04260500* NY 43.98556 -75.92472 423
4263000* USGS-04263000* NY 44.59944 -75.37889 315
4268000* USGS-04268000* NY 44.83889 -74.97889 269
40851325 USGS-040851325 WI 44.501111 -87.93611 663
103
40851385 USGS-040851385 WI 44.528611 -88.01 584
40869416 USGS-040869416 WI 43.097222 -87.96722 676
40871488 USGS-040871488 WI 42.989722 -87.95194 663
41482663 USGS-041482663 MI 42.941944 -83.84611 784
Supplemental Table 3.A.2. The order of the 7 sets of nested catchments that were split
from the respected HUC 8 watersheds
Nested Elevation
Site No Site ID State Latitude Longitude
order (ft)
1 4024090 USGS-04024090 MN 46.51744 -92.48214 984
2 4024093 USGS-04024093 MN 46.5155 -92.46269 928
1 4043140 USGS-04043140 MI 46.75104 -88.3618 1004
2 4043150 USGS-04043150 MI 46.8041 -88.31707 797
WI_MMSD-RI-04S &
1 4087000** WI 43.10001 -87.90897 679
USGS-04087000**
2 4087170 USGS-04087170 WI 43.02444 -87.89833 627
1 4087050 USGS-04087050 WI 43.20667 -88.03842 758
2 4087070 USGS-04087070 WI 43.12361 -88.04361 738
1 4087119 USGS-04087119 WI 43.04383 -88.00511 696
2 4087120 USGS-04087120** WI 43.04556 -87.99972 696
WI_MMSD-RI-11M &
3 4087142** WI 43.0339 -87.93425 627
USGS-04087142**
1 4087240 USGS-04087240 WI 42.75139 -87.82361 682
2 4087242 USGS-04087242 WI 42.73391 -87.78452 620
1 4161820** USGS-04161820** MI 42.61448 -83.02659 643
2 4164000 USGS-04164000 MI 42.57781 -82.95159 600
3 4165500 USGS-04165500 MI 42.59587 -82.90881 600
21MICH-500233 &
4 4165559** MI 42.59642 -82.82603 581
USGS-04165559**
104
4 Conclusion
Concentration-Discharge analysis continues to be a promising analytical method to
understand and monitor solute export processes. These studies illustrated some of the
hydrogeochemical and hydrological controls on multi-annual and annual concentration-
discharge patterns within Northeastern and North Midwestern United States and
throughout catchments with different land-covers.
The sites in the Great Lakes Basin offered valuable insight into catchment scale
changes in hydrologic inputs and processes over times, such as, solute availability and
dominant hydrological pathways (Chapter 3). Most solutes displayed significant (p < 0.05)
chemodynamic behavior, indicating that stream discharge is a significant driving factor in
controlling stream chemistry. Snow abundance also had a positive correlation with
concentration-discharge strength of SC but a negative correlation with nutrients and trace
elements, forecasting weaker concentration-discharge relationships for major ions under
future climate conditions. Land-use, especially cropland and urbanization, is correlated
with concentration-discharge strength and direction and intensity of solute export, related
to fertilizer additions and drainage networks facilitating rapid subsurface flow. A major
105
finding was that urbanization deteriorates the concentration-discharge relationship, likely
due to anthropogenic solute additions and disruption of natural flowpaths, decoupling the
concentration response from change in discharge.
The low frequency sampling resolution of the data used to study change in stream
water chemistry in the Great Lakes Basin presents some limitations to the interpretation of
the study results by limiting the interpretation to general long-term patterns (Chapter 3).
The results exclude many of the likely present temporal variations in stream chemistry, and
the observed trends should be viewed as overall patterns in stream water chemistry. The
impact of land-use was studied using the percent of watershed area without considering
proximity to the sampling location. Future related research could investigate how
proximity of cropland, urban areas, or forested areas to the sampling locations controls the
impact those land-cover types have on stream concentrations. Another potential future
study could investigate if there are thresholds that, once crossed, the studied land-uses have
a significantly stronger impact on stream chemistry. These future studies could provide
valuable information for sustainable development and land-use that minimizes the negative
effects of necessary structures and practices.
Results from these studies improve our understanding of current hydrologic processes
relating to streamflow generation and solute export while highlighting this importance of
solute vertical distribution and land-use. Climate scientists predict an average rise of 1.5 to
4.5ºC globally accompanied with changes in seasonal distribution, intensity, and duration
of precipitation (Mortsch & Quinn, 1996; D’Orgeville et al., 2014). High temperature
changes in winter will affect freeze thaw frequencies, ice formation, and the heat balance
of water bodies (D’Orgeville et al., 2014; Mortsch & Quinn, 1996). These changes will
modify biogeochemical and biotic processes, impacting solute availability, and negatively
impact the concentration-discharge behaviors positively associated with snow abundance.
Additionally, streamflow at higher latitudes is expected to increase by 10-40% over the
next 100 years, which could lead to a proportional increase in solute loadings while
maintaining similar concentrations in streams with chemostatic behavior (IPCC, Core
Writing Team, 2007). By better understanding current hydrologic and biogeochemical
106
processes and how these are modified by existing climate and land-use change, scientists
and land managers can predict how ecosystems might change under further modified
conditions.
4.1 REFERENCES
D’Orgeville, M., Peltier, W. R., Erler, A. R., & Gula, J. (2014). Climate change impacts
on Great Lakes Basin precipitation extremes. Journal of Geophysical Research,
119(18), 10.799-10.812. https://doi.org/10.1002/2014JD021855
Mortsch, L. D., & Quinn, F. H. (1996). Climate change scenarios for Great Lakes Basin
ecosystem studies. Limnology and Oceanography, 41(5), 903–911.
https://doi.org/10.4319/lo.1996.41.5.0903
IPCC, Core Writing Team. (2007). Climate Change 2007: Synthesis Report. Contribution
of Working Groups I, II and III to the Fourth Assessment Report of the
Intergovernmental Panel on Climate Change. Pachauri RK and Reisinger A (eds).
IPCC: Geneva, Switzerland; 104.
107