Camden M.
Hatley
a,
Brooklyn
Armijo
a,
Katherine
Andrews
a,
Christa
Anhold
a,
Jesse B.
Nippert
b and
Matthew F.
Kirk
*a
aDepartment of Geology, Kansas State University, 108 Thompson Hall, Manhattan 66506, Kansas, USA. Tel: +1-785-532-6724E-mail: mfkirk@ksu.edu
bDivision of Biology, Kansas State University, Manhattan 66506, Kansas, USA
First published on 23rd November 2022
Intermittent headwater streams are highly vulnerable to environmental disturbances, but effective management of these water resources requires first understanding the mechanisms that generate streamflow. This study examined mechanisms governing streamflow generation in merokarst terrains, a type of carbonate terrain that covers much of the central United States yet has received relatively little attention in hydrological studies. We used high-frequency sampling of precipitation, stream water, and groundwater during summer 2021 to quantify the contributions to streamflow from different water sources and characterize their short-term dynamics in a 1.2 km2 merokarst catchment at the Konza Prairie Biological Station (Kansas, USA). Mixing calculations using stable water isotopes and dissolved ions indicate that streamflow is overwhelmingly contributed by groundwater discharge from thin (1–2 m) limestone aquifers, even during wet periods, when soil water and surface runoff are generally expected to be more important. Relationships between hydraulic heads in the aquifers and their contributions to streamflow differed early in the study period compared to later, after a major storm occurred, suggesting there is a critical threshold of groundwater storage that the bedrock needs to attain before fully connecting to the stream. Furthermore, contributions from each limestone unit varied during the study period in response to differences in their hydrogeological properties and/or their stratigraphic position, which in turn impacted both the length of streamflow and its composition. Taken together, we interpret that the subsurface storage threshold and variation in aquifer properties are major controllers of flow intermittency in merokarst headwater catchments.
Environmental significanceIntermittent streamflow accounts for up to 50% of global discharge, but we currently lack a full understanding of the mechanisms which generate this streamflow in many physiographic settings. Characterizing streamflow generation for these settings will allow us to better manage and protect their water resources. In this study, we identify mechanisms governing streamflow generation in an intermittent headwater catchment in a merokarst terrain, an understudied subset of carbonate terrains that underlies much of the west- and south-central USA. Our results demonstrate how thin carbonate units within merokarst bedrock can serve as important hydrologic buffers and that the properties of individual units affect the timing and significance of their contribution to streamflow, and therefore the nature of intermittency in the stream. |
Merokarst terrains are a subset of carbonate terrains in which thin, relatively high permeability carbonate units are interbedded with low permeability mudstones or other siliciclastic units. A high proportion of carbonate terrains are merokarst, including significant parts of central North America.11,12 However, they remain largely understudied compared to other carbonate terrains and may be particularly sensitive to environmental change.13,14 The rapid transport of water generally seen in carbonate terrains amplifies the sensitivity of the environment to disturbances,15 exacerbating the responsiveness of headwater catchments therein.
Our research focuses on a headwater catchment within the Konza Prairie Biological Station (hereafter referred to as Konza Prairie) in northeastern Kansas, USA (Fig. 1). Konza Prairie is a part of the broader Flint Hills physiographic province of eastern Kansas and north-central Oklahoma and is classified as a merokarst terrain.14 Previous studies conducted at Konza Prairie reveal key pieces of information about merokarst streamflow generation and flow intermittency. Costigan et al.9 used long-term precipitation and discharge data (1987–2011) to characterize the overall surface water hydrology of Konza Prairie and observed that headwater discharge patterns were disconnected from precipitation patterns. They concluded that threshold-driven processes were causing this disconnect and ultimately leading to intermittent flow behavior. They offer the gradual closing of soil macropores and the filling of groundwater storage zones as two possible mechanisms that might be giving rise to the threshold dynamic, but their results did not allow them to test these hypotheses.
Fig. 1 Map of Konza Prairie Biological Station and watershed N04D. Locations of wells included in this study are indicated with open red circles. Each of these wells also contains a pressure transducer used to monitor groundwater levels. In the well identifiers, EIS 1 is used for wells screened in the Lower Eiss Limestone and EIS 2 is used for wells screened in the Upper Eiss Limestone. Wells 4-6 MOR, 4-6 EIS 1, and 4-6 EIS 2 are marked by a single circle because they are located within 1 m of each other. Approximate locations where select bedrock units outcrop in the stream channel are identified based on previous research.20 Limestone and shale are abbreviated Ls and Sh, respectively, in the legend. |
More recently, Sullivan et al.16 employed a principal-component analysis to identify potential end members that contribute solutes to streams on Konza Prairie. They concluded that the stream water is a mixture of four different sources: surface runoff, soil water, groundwater flowing through the carbonate aquifers, and groundwater flowing along the carbonate-mudstone contacts. Sullivan et al.16 did not quantify the relative contributions to streamflow from these endmembers in their study, but data collected by Keen et al.17 suggest that the groundwater inputs are dominant. They analyzed stable isotope compositions of water from soil cores as well as archived spring and surface water samples collected during 2010–2017. Water isotope compositions of stream samples closely resembled those of groundwater discharging at the spring.
These analyses provide valuable insight into the potential water sources that contribute to streamflow generation, but they do not reveal the short-term dynamics of these water sources that might help give rise to flow intermittency. Additionally, they do not capture separate contributions to streamflow from the different stacked carbonate units within Konza Prairie bedrock stratigraphy. The different carbonate units at Konza Prairie vary chemically, meaning that shifting contributions from these aquifers to the stream can impact stream chemistry.11,18 Furthermore, each aquifer can have significantly different hydrogeologic properties19 and outcrop patterns, suggesting that they might impart different flow pathways on water moving through them.
In this study, we seek to build on these previous analyses by investigating the relative importance of different water sources, their variation on short time scales, and how this variation might influence stream flow intermittency. To achieve this goal, we employed high frequency sampling of surface and groundwater at Konza Prairie during the 2021 growing season (April–July) to quantify streamflow source contributions and characterize their short-term dynamics. We collected stream samples on a weekly to hourly basis depending on flow conditions and applied mixing analyses to capture small timescale changes to source contributions. We established the endmembers for such mixing analyses using groundwater sampled from numerous monitoring wells along with long-term precipitation, stream, and groundwater chemistry records that are available in the extensive Konza Prairie data catalog. Lastly, we utilized sensors measuring high frequency discharge, precipitation, and groundwater elevation data to provide context for any changes to source contributions and aid in interpretations of threshold-driven behavior that might be controlling flow intermittency.
Konza Prairie receives an average of 835 mm of precipitation annually.23 Most of the precipitation falls during the growing season which runs from April–September24 (ESI Fig. SI1A†). Discharge is measured at a gauging station at the outlet of N04D and generally peaks during May–June, coinciding with the peak of precipitation. However, streamflow typically dries out by July, despite the persistence of frequent precipitation over the remainder of the growing season25 (ESI Fig. SI1B†). This drying is likely due to the watershed achieving maximum rates of evapotranspiration, which is estimated to account for 75% of the total annual water budget in Konza Prairie.26 Despite this tendency for the stream to dry out by early summer, flow can return at any later month given enough rainfall.
The geology of N04D is characterized by an alternating sequence of nearly horizontal 1–2 meter-thick limestone units (abbreviated ‘Ls’) and 2–4 meter-thick mudstones/shales11 (Fig. 2). The bedrock is Permian in age and the lowest (oldest) unit is the Cottonwood Ls, which crops out immediately below the gauging station (Fig. 1). The next two limestone units above the Cottonwood Ls are the Morrill Ls and Eiss Ls, with the Eiss Ls being further divided into the hydrogeologically distinct Upper Eiss Ls and Lower Eiss Ls members. The Morrill Ls and Eiss Ls crop out in the lower fourth of the catchment and have 31 groundwater monitoring wells completed within them (Fig. 1). The wells are made of 5 cm diameter PVC and have 61 cm long screens with gravel packs and bentonite seals.11 Pressure transducers are installed in three of the wells that are completed in the Morrill Ls, one of the wells completed in the Upper Eiss Ls, and one of the wells completed in the Lower Eiss Ls. These transducers record the pressure and temperature of water and atmosphere within the wells every five minutes, which we use to evaluate hydraulic heads.
Fig. 2 Konza Prairie bedrock units that crop out approximately within the lower two-thirds of watershed N04D. Wells included in this study are screened within the Upper Eiss Limestone Member and Lower Eiss Limestone Member of the Bader Formation, and the Morrill Limestone Member of the Beattie Formation, Council Grove Group, Permian System. Limestone and shale are abbreviated Ls and Sh, respectively, in the figure. Unit thickness ranges are shown in gray font beneath the name of each unit. This stratigraphic column is based on Fig. 2a of Macpherson.11 |
Groundwater flow within these limestone aquifers is mostly accommodated by fractures and solution-enlarged pores, meaning that hydraulic conductivity within each unit likely has extreme spatial heterogeneity. Indeed, well productivity varies considerably from site to site within the same unit11 and recent geophysical work has shown a large degree of heterogeneity in the distribution of water within each unit.14 Previous studies have measured the hydraulic conductivities of the Morrill and Eiss Limestones via rising-head slug tests in the wells and reported values of 10−5–10−3 m d−1 for the Morrill Ls, 10−8–10−5 m d−1 for the Lower Eiss Ls, and 10−5–10−4 m d−1 for the Upper Eiss Ls.19,20 These values fall in line with those typical of karstic and non-karstic limestones and dolomites,27 but given the aforementioned heterogeneity they should not be taken as measurements that can be accurately generalized across the entire spatial extent of each limestone.
Streamflow in the merokarst of Konza Prairie is predominantly generated by groundwater discharge from the limestone units via springs along the hillslopes and within the streambed. Because the limestones are interbedded with low permeability mudstones/shales, vertical movement of water between the limestones is likely restricted and thus each limestone unit acts as a semi-isolated aquifer predominately hosting horizontal groundwater flow as shown in the Graphical Abstract. Lithological differences between these aquifers can cause groundwater contained within each unit to have unique chemical compositions and flow properties.11,18 Surface runoff can also contribute to streamflow during storms, but it is likely limited in importance due to the presence of soil macropores and bedrock fractures which route precipitation into the subsurface instead.28,29
While the above hydrogeologic properties form the basic understanding of streamflow generation at N04D, water flow paths through the catchment can in fact be highly variable both in space and time. For instance, whereas each limestone unit is generally thought of as an isolated aquifer with limited flow across its neighboring mudstones, recent studies have shown that significant vertical transport can occur in localized areas. Barry20 first found evidence of this cross-unit transport when a tracer dye injected in the Morrill Ls appeared in significant concentrations within groundwater discharging from the Cottonwood Ls below it. They posited that this flow was accommodated by a buried collapse feature. Sullivan et al.14 followed up on this finding with a geophysical survey that revealed many such areas of enhanced vertical flow, suggesting that this phenomenon might be occurring extensively throughout the catchment. As further evidence of the complicated water flow paths at N04D, the stream gauge at the outlet of the watershed generally records zero-flow between July and February of the following year, but a portion of the stream above the gauge experiences flow nearly perennially.11 This perennial portion begins where the outcropping unit in the stream transitions from Upper Eiss Ls to Lower Eiss Ls and continues until the outcropping unit transitions from the Stearns Shale to the Morrill Ls. These observations suggests that a temporary flow path occurs during this portion of the year in which water discharges from the Upper Eiss Ls, flows some distance in the channel, then reinfiltrates into the subsurface and recharges the Morrill Ls.
Surface water was collected from the stream immediately downstream of the flume at the catchment outlet. The frequency at which we collected stream samples varied based on the conditions in the watershed: during low flow conditions towards the end of the sampling period this frequency was roughly weekly, during the peak of the wet season this frequency was roughly daily, and during selected storms this frequency was roughly hourly. We sampled hourly during five storms, which fell on May 8–9, May 26–27, June 11, June 24–25, and July 15. In addition to the stream, for each of these events we also collected precipitation with a ball-in-funnel type collector.30 These events were chosen to represent a range of storm intensities and antecedent watershed conditions. Sampling for these events began at the onset of precipitation and continued until approximately two hours after precipitation ceased. In total, we collected 44 stream samples over the 93 days sampling period.
For each stream sample, we measured water temperature, pH, specific conductivity, and dissolved oxygen using field probes placed directly into the stream. Probes were recalibrated for all parameters before each sampling trip. Immediately upstream of the probes, we collected and filtered water through 0.45 μm syringe filters into two 30 mL high density polyethylene bottles for major cation and anion analysis, a 10 mL glass vial for stable water isotope analysis, and a 60 mL amber glass bottle for analysis of non-purgeable organic carbon (NPOC) and total dissolved nitrogen (TN). The filtration equipment and polyethylene bottles that we used were either brand new or were acid washed prior to sampling. Amber glass bottles were baked overnight at 120 °C prior to sampling. Every piece of sampling equipment was rinsed with sample water three times before being filled up with sample. Within a day of sampling, we preserved cation and NPOC/TN samples with three drops of concentrated trace element grade HNO3 and analytical grade HCl, respectively. All samples were stored on ice during extended field work outings and at 2 °C in a refrigerator upon return to the laboratory.
For each of the hourly-sampled storms, we left out a rain collector over the entire duration of the storm and, when volumes were sufficient, divided the resulting precipitation into two 10 mL aliquots for replicate analysis of stable water isotopes.
We collected groundwater samples on six different days throughout the sampling period. For each groundwater sampling event, we collected water from the three wells in the 4-6 well nest (see Fig. 1 for well locations). Each of these wells is screened in one of the three limestone aquifers which crop out in the lower reach of the stream: the Morrill Ls (well name 4-6 MOR), Lower Eiss Ls (4-6 EIS 1), and Upper Eiss Ls (4-6 EIS 2). On three of the sampling trips, we also collected groundwater from the 3-5 and 3-5-1 wells which draw from the Morrill Ls (3-5 and 3-5-1 MOR). Before sampling each well, we purged a volume of water equal to two times the standing water volume via bailing. We then bailed water into two glass beakers, rinsing each beaker three times with sample water before filling them up for collection and analysis at each well. We used one beaker for field probe measurements and the other for collection into sample bottles, as described for stream samples. At the 3-5 and 3-5-1 MOR wells, we only collected water for analysis of major ions and not stable water isotopes, NPOC, or TN. The 4-6 well completed in the Lower Eiss Ls often did not yield enough water to fill all four collection bottles, and at these times, we also only collected water for major ion analysis.
The δ18O and δD of water isotope samples were measured at 0.1 and 0.5‰ precision using a Picarro L-i2130 water analyzer in the Stable Isotope Mass Spectrometry Laboratory at the Kansas State University Division of Biology. Results are presented as per mil (‰) relative to Vienna Standard Mean Ocean Water (VSMOW).
To quantify contributions to streamflow from soil water, we relied on soil water chemistry data from Tsypin and Macpherson,29 who collected eight samples from the A horizon (∼17 cm deep) and nine samples from the B horizon (∼150 cm deep) in N04D during the summer of 2010. We recognize that there is some uncertainty introduced by establishing an endmember from one season of field data, but we are not aware of other available soil water chemistry data sets for our study site. Nevertheless, the fact that both Tsypin and Macpherson29 and the present study were conducted during the summer suggests that their soil data may be generally applicable for this analysis.
XM = XAfA + XBfB +…+ Xnfn | (1) |
fA + fB +…+ fn = 1 | (2) |
Selecting which tracers to employ and which water sources to use as endmembers is the key consideration for mixing analysis. While endmembers can be statistically determined from the stream chemistry data alone,35 they can also be chosen ahead of time based on previous understanding of the system in question.36 Based on the work of Sullivan et al.,16 we assumed that N04D stream water is composed of precipitation (incorporating precipitation falling directly into the stream channel and surface runoff), soil water, and multiple groundwater endmembers, each sourced from different limestones.
Having more than three endmembers gives rise to significant mathematical complexity.35 To avoid this complexity, we separated our mixing analysis into two stages. In the first stage, we used stable water isotope ratios (δ18O and δD) as tracers in a two-endmember mixing calculation that determined the fractional contributions of precipitation and subsurface water in streamflow, an analysis known as isotope hydrograph separation.37 The δ18O and δD values used as the endmember compositions varied throughout the sampling period. For calculations at each stream sampling time point, we used the isotope ratios from the most recent precipitation sample as the precipitation endmember and the mean isotope ratios from the most recent groundwater samples as the subsurface endmember. We conducted the mixing calculations separately for the two isotopic tracers and then averaged the results. As a last step in this stage of the analysis, we calculated discharge-weighted averages of the fractional contribution of precipitation to streamflow across each individual storm and the entire study period to understand the overall importance of runoff to streamflow.
Some error is likely introduced to these calculations by the fact that the precipitation endmember is established from a single composite sample of precipitation collected over the course of a storm. Isotopic ratios in precipitation are known to vary throughout storms,38 so our use of a single, volume-averaged endmember means our calculation of the fractional contribution of precipitation in streamflow is not entirely accurate on the hourly time scale of some of our stream samples. While this creates uncertainty in the distribution of precipitation inputs over time, we maintain that the proportions of precipitation in streamflow averaged across the duration of each storm are accurate because our precipitation end-member values were measured on composite samples. Given that the error introduced by within-storm variability in precipitation stable water isotopic ratios is typically only 0.5–2‰ for δ18O values,39–41 and that previous studies have reported a general insignificance of surface runoff at N04D,42,43 we reason that the uncertainty is low in the distribution of precipitation inputs over time and thus does not impact the conclusions we make about overall behaviors in the catchment.
After accounting for the contribution of precipitation to each stream sample using stable water isotopes, the second stage of our analysis was to then use major ion concentrations as tracers to calculate the contributions from the remaining three endmembers: the soil water and two groundwater sources. We first adjusted the ion concentrations in each stream sample to remove the influence of precipitation via a rearrangement of eqn (1):
(3) |
The flume at the catchment outlet recorded flow for most of the sampling period, with the exception of a dry period between July 5 and July 14 (Fig. 3). Flow was modest at the beginning of the sampling period, averaging 0.004 m3 s−1 between April 15 and May 15. A large storm on May 16 (96.5 mm as measured at the catchment outlet) caused discharge to spike to its maximum for the year of 4.69 m3 s−1, which then dropped to an average of 0.077 m3 s−1 on the following day. An equivalent water depth of 39.6 mm exited the catchment via the stream in the 24 hours following this large storm, and 74.2 mm exited by the beginning of the next major storm on May 26. After the May 16 storm, discharge generally decreased over the remainder of the sampling period and dried out completely on July 5. Occasional storms, such as the one on May 26, brought small spikes in discharge during this drying period. After the stream completely dried out, the large storm on July 15 caused a week-long return of modest flow. Across the whole sampling period, 1.49 × 105 m3 of water exited N04D via the stream, which equates to 124.2 mm of water depth over the whole catchment.
The storms that we selected for hourly sampling are marked in Fig. 3 and covered a range of hydrological conditions. Total precipitation that fell for these events ranged from 2.80 mm to 80.0 mm and their spacing throughout the sampling period meant that stream flow levels were different for each. One event coincided with the low-flow period before the large May 16 storm, three events occurred over the progressive drying period after the May 16 storm, and one event covered the large July 15 storm that produced flow after the prolonged dry period.
Fig. 4 Groundwater head in meters above sea level (m.a.s.l.) within each of the sampled wells from February 23, 2021 to August 7, 2021. Daily total precipitation is also shown. Well IDs refer to locations in Fig. 1, with EIS 1 corresponding to the Lower Eiss Ls and EIS 2 corresponding to the Upper Eis LS. The vertical dashed lines outline the sampling period of this study. Sharp downward spikes in the groundwater head of 4-6 EIS 1 occurred in response to groundwater removal from the well during sampling. These spikes are largely not visible in the head profiles of the other wells, reflecting differences in hydraulic conductivities between limestone units. |
The wells in the Morrill Ls (3-5 MOR, 3-5-1 MOR, and 4-6 MOR) all responded similarly to one another. As with 4-6 EIS 2, heads in each of these wells increased as a result of the three moderate storms in March and April. Unlike 4-6 EIS 2, however, these increases were relatively drawn out rather than instantaneous spikes, and the increases associated with each storm did not fully recess to pre-storm levels prior to the onset of the next storm. Heads in these wells were still in the process of recessing from these March and April storms when the large May 16 storm occurred. As a result of the May 16 storm, water levels instantaneously spiked upwards similarly to 4-6 EIS 2, but the magnitudes of their increases were not as large (0.938 m in 3-5 MOR, 0.563 m in 3-5-1 MOR, and 0.339 m in 4-6 MOR). Following this storm, heads resumed their drawn out recessions which continued for well over a month and did not land on low, stable values until late July or early August. As with head in 4-6 EIS 2, moderate storms which occurred during this recession phase produced small spikes in head, but in the Morrill wells these interruptions were lower in magnitude and were again followed by slower recessions.
These differences in groundwater head behavior indicate a difference in how water moves through the Upper Eiss Ls and Morrill Ls, despite previous slug tests resulting in similar hydraulic conductivities for them.19,20 One possibility is that the slug tests do not capture the spatial heterogeneity in fractures and solution conduits that causes the Upper Eiss Ls to have an overall higher hydraulic conductivity than the Morrill Ls. This would explain why groundwater head fluctuations in March and April were flashier for the Upper Eiss Ls than the Morrill Ls, and why the recessions following these storms took longer in the Morrill Ls. This hypothesis is complicated, however, by the fact that wells in both the Upper Eiss Ls and the Morrill Ls seem to have similar increasing pressure responses to storms occurring after April. Although the spikes produced in the Morrill Ls wells were lower than those in 4-6 EIS 2, they were instantaneous spikes nonetheless and not the relatively drawn-out peaks like those in March and April. This similarity in response times to the onset of storms in May and beyond suggests that overall hydraulic conductivities in the two units might not be significantly different.
An alternative explanation is that the differences in behavior are caused by the simple stratigraphic order of the two units. The Upper Eiss Ls is positioned above the Morrill Ls, so it would naturally be more impacted by the addition of recharge from the surface. This relationship could create a stronger pressure response in the groundwater head in 4-6 EIS 2 than in the underlying Morrill Ls wells, consistent with observed head responses to storms throughout the study period. Similarly, the slower recessions for the Morrill wells can also be explained by stratigraphic order. While horizontal groundwater flow is favored in these limestone aquifers, the presence of spatially discrete vertical connections between different units causes higher units to ‘leak’ into lower ones.16,20 Thus, the Morrill Ls might be receiving continuous input of groundwater from the overlying Upper Eiss Ls as it recesses following a storm, effectively extending the amount of time it takes to return to a low, stable head.
Compared to the wells discussed thus far, changes to groundwater head in the 4-6 EIS 1 well (Lower Eiss Ls) occurred very slowly. Whereas the wells in the Upper Eiss Ls and Morrill Ls had near instantaneous pressure responses following the May 16 storm, 4-6 EIS 1 did not achieve its maximum head until May 22, six days after the others. Further, the well took several days to recover the water that was removed during each sampling (Fig. 4). This is in accordance with previous studies concluding that the Lower Eiss Ls has a significantly lower conductivity than the adjacent limestone units.11,19
Fig. 5 Box plots of select major ion concentrations in all surface and groundwater samples. Hinges mark the interquartile range of each group of samples, and stems extend to 1.5 times the interquartile range. Individual points mark outlier samples that had concentrations greater or less than 1.5 times the interquartile range. Extreme outlier points are not shown for K+ and Cl− to show finer variation in concentrations. Complete data are available in the ESI (Tables SI3 and SI4†). |
Some samples of the stream water contained Cl−, Na+, and K+ concentrations that are 3–50 times higher than the median (Fig. 5). These outliers all came from samples collected during the May 8–May 9 storm, which was the first major rain event to occur after N04D was burned on April 12. Previous studies in various grassland and non-grassland environments have shown that fires can mobilize large amounts of nutrients from soil and vegetation that can later be leached into surface waters.44–46 Runoff during the May 8–May 9 storm could have carried these recently mobilized nutrients into the stream, offering one potential explanation for the spikes in Cl−, Na+, and K+ levels. Why other major ions such as Ca2+ and Mg2+ did not spike in concentration as well is unclear.
Values of pH were generally higher in stream water samples (avg. 7.69) than in the groundwater samples (avg. 6.99), and vice versa for alkalinity (ESI Tables SI3 and SI4†). This difference can be attributed to CO2 outgassing from the water after it is discharged from the subsurface.47 Stream water NPOC values averaged 3.00 mg L−1 and did not vary significantly with discharge (ESI Table SI3†). NPOC was not measured in groundwater due to use of dye tracers for previous research.20 For both surface water and groundwater samples, the concentrations of NO2−, NO3−, Br−, NH4+, and TN were frequently below the detection limits of our methods and thus excluded from the remainder of the analysis.
Fig. 6 Stable water isotopes in stream water, groundwater, and precipitation samples collected at N04D. Solid line denotes the global meteoric water line (GMWL) (δD = 8 × δ18O + 10), and the dashed line denotes the linear regression of all plotted samples. Isotope data are available in the ESI (Tables SI3–SI5†). |
The δ18O and δD values of the stream water samples generally fell between the relatively depleted values of the groundwater and the relatively enriched values of the precipitation, but positioned closer to the former than the latter (mean δ18O of the stream samples = −5.4‰). The isotopic ratios of the stream water only deflected closer to the precipitation values at the peak of major storms. The maximum δ18O value of our stream water samples was −4.4‰, which occurs at the peak of the May 8–May 9 storm. These results suggest that, during our sampling events, the stream water was comprised mostly of groundwater with somewhat minor amounts of more recent precipitation (i.e., surface runoff) mixed in, consistent with findings of Keen et al.17
Date range | Event endmember | Pre-event endmember | Fraction event | |||
---|---|---|---|---|---|---|
Start | End | δ 18O | δD | δ 18O | δD | — |
5/6 | 5/14 | −1.8 | −6.5 | −5.7 | −36 | 0.181 |
5/22 | 6/1 | −1.25 | −10 | −5.7 | −35 | 0.018 |
6/8 | 6/14 | 0 | 12 | −5.6 | −35 | 0.002 |
6/22 | 6/27 | −1.3 | −7 | −5.6 | −35 | 0.021 |
7/15 | 7/16 | −4.4 | −22 | −5.6 | −36 | 0.500 |
The second stage of the mixing analysis was a three-endmember calculation between contributions from soil water, groundwater from the Morrill Ls as represented in the 3-5-1 MOR well, and groundwater from the Upper Eiss Ls as represented in the 4-6 EIS 2 well (Fig. 7). We selected 3-5-1 MOR to represent the Morrill Ls because it had a stable chemistry over the entire sampling period, displaying little variation in all chemical parameters (Fig. 5 and ESI Table SI4†). The 3-5 MOR and 4-6 MOR wells receive significant amounts of recharge directly from the stream,11,20 making them unsuitable for stream water mixing calculations and leaving 3-5-1 MOR as the better representation of ‘pure’ groundwater held in the Morrill Ls. We did not include the Lower Eiss Ls as an endmember because it has a much lower hydraulic conductivity than the other aquifers and has been shown to have little interaction with the stream.11,19
Fig. 7 Bivariate plot of Ca2+/Mg2+vs. SO42− for the three endmembers selected for mixing analysis, with stream water sample chemistries adjusted for the influence of precipitation as described in the text. The average values of Ca2+/Mg2+ and SO42− that were used as the endmember compositions are marked with red diamonds, and error bars show one standard deviation from the mean. Dashed lines denote the mixing envelope. Endmember compositions (Table SI2) and an identical plot with raw stream chemistry data (Fig. SI2) are available in the ESI.† |
To calculate the fractional contributions from the endmembers, we used SO42− concentration and the mass ratio of Ca2+/Mg2+ as tracers. We selected these tracers based on the criteria laid out in the methods: SO42− and Mg2+ concentrations both displayed high variation between the groundwater wells selected as endmembers and low variation within each well (Fig. 5). Using the concentration ratio Ca2+/Mg2+ instead of the concentration of either ion made the tracer more robust to the effects of evaporation, sorption, and mineral precipitation in the stream channel. Indeed, long-term data from N04D has shown that there is little variation in Ca2+/Mg2+ values between stream samples collected on the same day from different points along the stream reach.52 Because we already accounted for the contributions to streamflow from precipitation in the first stage of the analysis, we adjusted the SO42− and Ca2+/Mg2+ values of all stream samples using eqn (3) to remove the influence of precipitation on these concentrations. All adjusted stream values were successfully bounded by a mixing envelope created by our selected endmembers in a SO42−vs. Ca2+/Mg2+ plot (Fig. 7). Endmember compositions (Table SI2†) and an identical plot without adjusted stream compositions (Fig. SI2) are available in the ESI.†
A major assumption that we made in the mixing analysis of this study is that groundwater contributions to the stream were only coming from the Upper Eiss Ls and the Morrill Ls. While we provided justification above for excluding the Lower Eiss Ls, N04D contains a total of six different limestone formations that outcrop above the flume,14 each of which could potentially contribute water to streamflow. However, the Eiss and Morrill members are the stratigraphically lowest of these limestones and have much more contact area with the stream channel than the higher units in the steeper upper reaches of the catchment. Several dozen springs drain the limestone units of N04D20 primarily in the upper reaches of the watershed, providing a potential path for water to reach the stream from these upper limestones without needing direct contact with the channel. However, these springs only produce water ephemerally and, based on long-term monitoring of stream and groundwater chemistry at the site,52 their contribution to streamflow is likely very small. Specifically, the stream water chemistry adjacent to the Upper Eiss Ls is generally very similar to the groundwater from the Upper Eiss Ls in terms of major ions and temperature. Thus, we reason that much of the water from the upper limestone units likely makes its way into the Eiss or Morrill members eventually.
Fig. 8 Fractional contributions to streamflow from precipitation, the two groundwater endmembers, and soil water through the sampling period. Stream discharge is marked with the light blue line for reference to flow conditions. Dotted lines indicate extended periods of time between samplings when short-term dynamics are not captured. Mixing calculation data are available in the ESI (Table SI6†). |
While groundwater dominated streamflow for most of the sampling period, our mixing calculations indicate that contributions from each aquifer varied in relative importance through the study period (Fig. 8). The Upper Eiss Ls contributed the most water to streamflow among all sources for most of the sampling period, especially during high flow conditions (daily average stream discharge > 0.02 m3 s−1) when its fractional contribution could exceed 0.75. Given that the Upper Eiss Ls is stratigraphically higher and outcrops over a larger surface area than the Morrill Ls, it follows that it generally contributed more to streamflow across the sampling period. However, as the stream dried following the May 16 storm, contributions from the Morrill Ls slowly overtook those from the Upper Eiss Ls and exceeded them by late June when daily average stream discharge had fallen below 0.01 m3 s−1.
This pattern of dominance by the Upper Eiss Ls during high flow conditions and by the Morrill Ls during prolonged dry periods might be related to the differences between the units in pressure response and recession time that we discussed in Section 4.2 (Fig. 4). The Upper Eiss Ls generally has higher magnitude spikes in groundwater head and more rapid rates of recovery towards pre-storm levels than the Morrill Ls, which we hypothesized might be due to enhanced fracturing and dissolution in the Upper Eiss Ls, giving it a greater overall hydraulic conductivity. According to this hypothesis, the Upper Eiss Ls dominates streamflow under wet conditions by quickly routing recharge into the channel, but in the absence of further additions of precipitation this output quickly dwindles. In comparison, the slightly slower responding Morrill Ls releases its water in more of a steady trickle, which has a lower volumetric flow rate than the discharge from the Upper Eiss Ls but can be sustained for longer and thus becomes proportionally more important as the Upper Eiss Ls depletes.
Even if the overall hydraulic conductivities of the two limestone units are not significantly different, this pattern in their fractional contributions can also be explained by their stratigraphic order. As the stratigraphically lower unit, much of the storm response in the Morrill Ls relies on vertical transport of groundwater from the overlying Upper Eiss Ls and through the Stearns Shale that separates them. Though some spatially discrete structures exist throughout the catchment and may be able to accomodate rapid vertical transport,14,20 vertical flow between limestone units will likely still be overall slower than horizontal flow within each individual limestone unit. Thus, when a storm creates a pressure response in the groundwater within the Upper Eiss Ls, much of the resulting flow will be directed horizontally towards the stream and relatively less will move vertically into the Morrill Ls. As a result, discharge from the Upper Eiss Ls will be more important to streamflow than discharge from the Morrill Ls in the immediate aftermath of storms. During dry periods, however, the slower vertical movement of water into the Morrill Ls sustains an elevated head for a longer period of time, allowing it to continue discharging into the stream after the Upper Eiss Ls has drained. The net results of either hydraulic conductivity or stratigraphic order controlling the observed behavior in streamflow generation are the same, but further studies would be helpful in evaluating the extent to which either hypothesis is supported.
Surface runoff only comprised a significant proportion of the total stream discharge during the peak of storms occurring when the stream was dry (July 15 storm) or under low-flow conditions (May 8–May 9 storm) (Fig. 8). The May 8–May 9 and May 26–May 27 storms had very similar amounts of precipitation (22.1 mm and 22.5 mm, respectively), but were comprised of an order of magnitude difference in proportions of surface runoff (0.181 and 0.018). The difference between these storms was the flow conditions in the stream prior to precipitation. In the May 8–May 9 storm, the stream was nearly dry (discharge = 0.0016 m3 s−1) and in the May 26–May 27 storm flow was over an order of magnitude larger (discharge = 0.0272 m3 s−1). Thus, the importance of surface runoff in the stream compared to inputs from subsurface sources seems to diminish drastically during high flow conditions.
We hypothesize that this behavior reflects a ‘fill and spill’ dynamic to the hydrology of N04D: when the watershed is dry, precipitation infiltrating into the subsurface merely ‘fills’ the subsurface storage units but does not yet push groundwater into the stream, causing the precipitation that does run over the surface to comprise a relatively large proportion of streamflow. Once these subsurface units surpass some critical thresholds of storage, however, they ‘spill’ into the stream and dwarf the surface runoff in terms of fractional contribution.
Fill and spill models were originally conceived to explain subsurface stormflow on the hillslope scale,53 but subsequent studies have shown that fill and spill dynamics can operate across entire watersheds and are heavily controlled by sediment and rock properties.54,55 Gutierrez-Jurado et al.8 modeled streamflow generation in a non-karstic catchment and showed that the gradual filling of groundwater storage zones can be a major controller on flow intermittency. The hydrogeologic regime of merokarst terrains is very different than that of the sandy aquifer that Gutiérrez-Jurado et al. modeled. However, Costigan et al.9 hypothesized that such a fill and spill dynamic was operating in N04D and concluded that streamflow and precipitation only became synchronized once the watershed attained a threshold amount of subsurface water storage. Our results on the dynamics of surface runoff agree that fill and spill behavior characterizes streamflow generation at N04D, and further analyses on groundwater inputs to the stream as discussed below seem to add support to this notion.
The SO42− and Mg2+ concentrations in groundwater samples from the 4-6 EIS 2 well in early May were over one standard deviation higher than their long-term averages, potentially reflecting an increased residence time of groundwater in the subsurface. SO42− concentration may have increased in response to dissolution of gypsum or anhydrite whereas Mg2+ levels may have been high as a result of exchange with Ca2+ in clay minerals or potentially chlorite dissolution.11,57 The lack of increased Ca2+ concentrations in these same samples is likely due to the fact that the groundwater generally equilibrated with calcite, so an increase in contact time will not result in further dissolution without some accompanying change in temperature or pH.
Regardless of the mineral sources of these ions, the increased contact time with bedrock for these early May samples implies that the Upper Eiss groundwater was not being extensively flushed into the stream at this time. This implied lack of a hydrologic connection suggests that the groundwater units had not yet reached the critical storage threshold required to ‘spill’ water into the stream. As such, this might explain why our calculated fractional contributions for this early period are different from what is predicted in the F–Q relationships and it also means that the results of our mixing calculations (Fig. 8) might not be accurate for this same period. The May 16 storm input a massive amount of water into the system, however, making it likely that the subsurface units then reached their storage thresholds, fully connected to the stream, and were able to spill water into it according to the established F–Q relationships. Based on these results, we interpret that fractional contributions to streamflow from the groundwater units at N04D can be reasonably estimated from the F–Q relationships so long as the aquifers have surpassed some threshold water storage amount.
Establishing relationships between groundwater contributions and hydraulic head in each limestone unit allows us to better understand the threshold that governs subsurface connection to the stream. In the Upper Eiss Ls, volumetric discharge of water to the stream (fractional contribution multiplied by stream discharge) varies linearly with groundwater hydraulic head in the 4-6 EIS 2 well (Fig. 10A). In the Morrill Ls, a linear relationship between discharge into the stream and hydraulic head in 3-5-1 MOR only holds for samples collected after the large May 16 storm (Fig. 10B). The samples collected prior to May 16 all plot well below the linear trendline of the post-May 16 samples, indicating that discharge into the stream is much lower than predicted by the relatively high groundwater heads in 3-5-1 MOR at this time.
One interpretation of this discrepancy is that it is further evidence of a disconnect between the groundwater and the stream during these earlier sampling points. After the May 16 storm, the behavior switches to the predictable linear relationship, possibly indicating that a storage threshold was achieved. The clear difference in fractional contribution values for the samples collected from the two periods might also be related to the fact that our mixing model may be inaccurate for the pre-May 16 samples, but in this case, we would also expect there to be a similar difference in the samples for the Upper Eiss Ls. The lack of a similar behavior switch in the Upper Eiss Ls might instead suggest that the threshold dynamics are different for these two aquifers, and that the Upper Eiss Ls requires much less addition of water to connect to the stream than the Morrill Ls does. A higher overall hydraulic conductivity in the Upper Eiss Ls or its higher stratigraphic position with respect to the Morrill Ls might cause this difference in threshold behavior, much as it causes the differences in hydraulic head responses to storms as discussed previously.
Another interpretation of the discrepancy between discharge into the stream and well hydraulic head for the Morril Ls is that head, or at least head in a single well such as 3-5-1 MOR, is not a suitable measure for determining if the storage threshold has been attained in the aquifers. The wells only measure hydraulic head at a single point for each limestone unit which can crop out along several dozens of meters of stream reach. The extreme spatial heterogeneity in fractures and solution conduits in these limestones means that head measured at one point may be significantly different from head measured at another point which exists along a more conductive flow path. Studies employing greater well control or geophysical techniques to constrain groundwater head across a greater spatial extent could be used to quantify the storage thresholds that control groundwater connection to the stream.
We interpret that groundwater discharge dominates headwater streamflow generation in the merokarst of N04D, during both high and low flow conditions (Fig. 8). Contributions from surface runoff are only important during some storms, and contributions from soil water are volumetrically insignificant across the entire wet season. These behaviors are in contrast to the widespread shallow-and-deep hypothesis of streamflow generation, which states that contributions to streamflow from shallow units (i.e. surface runoff and soil water) become increasingly more important relative to deep units (i.e. groundwater) under wet conditions.58,59 While studies conducted at non-karstic catchments of similar size to N04D generally support the shallow-and-deep hypothesis,60,61 we hold that the karstic properties of N04D prevent the development of such a regime. Thin soils place an upper limit on the significance of shallow flow paths, while soil macropores and bedrock fractures rapidly route water to limestone aquifers ensuring that deep flow paths remain dominant even during very wet conditions. Given that groundwater flow paths are critical in sustaining streamflow, our finding that there are groundwater storage thresholds which govern subsurface connections to the stream then implies that these thresholds are major controllers on stream intermittency. Significant, sustained streamflow can only occur when enough precipitation has fallen in a short enough amount of time to surpass the groundwater storage thresholds.
The thinly bedded nature of the limestones at N04D means that its hydrological behavior also deviates from what is typically seen in massive karst systems. Confining mudstone layers largely restrict vertical movement of water between limestones and the unique depositional history of each individual limestone bed means that each one might possess unique hydrogeological properties. Whether caused by differences in hydraulic conductivity or simple stratigraphic position, our results show that each of our studied limestone aquifers, which are vertically separated by only a couple of meters from one another, respond to recharge and interact with the stream in different ways (Fig. 4). This heterogeneity expands the capacity of the subsurface to act as a hydrological buffer, in that each limestone unit releases recharge into the stream at a different rate to effectively extend the amount of time that the stream can flow. Consistent with this observation, streams across the Flint Hills region of Kansas have greater proportions of baseflow in total streamflow than streams in adjacent, non-merokarstic regions,62 suggesting that the merokarst bedrock does improve hydrological buffering. Because the sensitivity of a catchment to different lengths of drought is strongly influenced by the hydraulic conductivity of the underlying rocks/sediment,63 the existence of a spectrum of conductivities might make merokarst terrains such as N04D more resilient to drought overall. This effect should not be overstated, as flow at N04D is still highly intermittent, but it might act to mitigate some of the risk that climate change poses to intermittent streams worldwide.6,64
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d2va00191h |
This journal is © The Royal Society of Chemistry 2023 |