USGS Science⌕ Search

SEARCH · USGS Science

Results for “Hydrological Processes”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 523 records · Page 29Linked to original sources

Sediment source apportionment in Laurel Hill Creek, PA, using Bayesian chemical mass balance and isotope fingerprinting

A Bayesian chemical mass balance (CMB) approach was used to assess the contribution of potential sources for fluvial samples from Laurel Hill Creek in southwest Pennsylvania. The Bayesian approach provides joint probability density functions of the sources' contributions considering the uncertainties due to source and fluvial sample heterogeneity and measurement error. Both elemental profiles of sources and fluvial samples and 13 C and 15 N isotopes were used for source apportionment. The sources considered include stream bank erosion, forest, roads and agriculture (pasture and cropland). Agriculture was found to have the largest contribution, followed by stream bank erosion. Also, road erosion was found to have a significant contribution in three of the samples collected during lower-intensity rain events. The source apportionment was performed with and without isotopes. The results were largely consistent; however, the use of isotopes was found to slightly increase the uncertainty in most of the cases. The correlation analysis between the contributions of sources shows strong correlations between stream bank and agriculture, whereas roads and forest seem to be less correlated to other sources. Thus, the method was better able to estimate road and forest contributions independently. The hypothesis that the contributions of sources are not seasonally changing was tested by assuming that all ten fluvial samples had the same source contributions. This hypothesis was rejected, demonstrating a significant seasonal variation in the sources of sediments in the stream.

Pennsylvania↗

Water's Way at Sleepers River watershed – revisiting flow generation in a post-glacial landscape, Vermont USA

The Sleepers River Research Watershed (SRRW) in Vermont, USA, has been the site of active hydrologic research since 1959 and was the setting where Dunne and Black demonstrated the importance and controls of saturation-excess overland flow (SOF) on streamflow generation. Here, we review the early studies from the SRRW and show how they guided our conceptual approach to hydrologic research at the SRRW during the most recent 25 years. In so doing, we chronicle a shift in the field from early studies that relied exclusively on hydrometric measurements to today's studies that include chemical and isotopic approaches to further elucidate streamflow generation mechanisms. Highlights of this evolution in hydrologic understanding include the following: (i) confirmation of the importance of SOF to streamflow generation, and at larger scales than first imagined; (ii) stored catchment water dominates stream response, except under unusual conditions such as deep frozen ground; (iii) hydrometric, chemical and isotopic approaches to hydrograph separation yield consistent and complementary results; (iv) nitrate and sulfate isotopic compositions specific to atmospheric inputs constrain new water contributions to streamflow; and (v) convergent areas, or ‘hillslope hollows’, contribute disproportionately to event hydrographs. We conclude by summarizing some remaining challenges that lead us to a vision for the future of research at the SRRW to address fundamental questions in the catchment sciences.

Vermont↗

A comparison of high-resolution specific conductance-based end-member mixing analysis and a graphical method for baseflow separation of four streams in hydrologically challenging agricultural watersheds

Quantifying the relative contributions of different sources of water to a stream hydrograph is important for understanding the hydrology and water quality dynamics of a given watershed. To compare the performance of two methods of hydrograph separation, a graphical program [baseflow index (BFI)] and an end-member mixing analysis that used high-resolution specific conductance measurements (SC-EMMA) were used to estimate daily and average long-term slowflow additions of water to four small, primarily agricultural streams with different dominant sources of water (natural groundwater, overland flow, subsurface drain outflow, and groundwater from irrigation). Because the result of hydrograph separation by SC-EMMA is strongly related to the choice of slowflow and fastflow end-member values, a sensitivity analysis was conducted based on the various approaches reported in the literature to inform the selection of end-members. There were substantial discrepancies among the BFI and SC-EMMA, and neither method produced reasonable results for all four streams. Streams that had a small difference in the SC of slowflow compared with fastflow or did not have a monotonic relationship between streamflow and stream SC posed a challenge to the SC-EMMA method. The utility of the graphical BFI program was limited in the stream that had only gradual changes in streamflow. The results of this comparison suggest that the two methods may be quantifying different sources of water. Even though both methods are easy to apply, they should be applied with consideration of the streamflow and/or SC characteristics of a stream, especially where anthropogenic water sources (irrigation and subsurface drainage) are present.

Hydrological Processes↗

Groundwater – The disregarded component in lake water and nutrient budgets. Part 2: effects of groundwater on nutrients

Lacustrine groundwater discharge (LGD) transports nutrients from a catchment to a lake, which may fuel eutrophication, one of the major threats to our fresh waters. Unfortunately, LGD has often been disregarded in lake nutrient studies. Most measurement techniques are based on separate determinations of volume and nutrient concentration of LGD: Loads are calculated by multiplying seepage volumes by concentrations of exfiltrating water. Typically low phosphorus (P) concentrations of pristine groundwater often are increased due to anthropogenic sources such as fertilizer, manure or sewage. Mineralization of naturally present organic matter might also increase groundwater P. Reducing redox conditions favour P transport through the aquifer to the reactive aquifer-lake interface. In some cases, large decreases of P concentrations may occur at the interface, for example, due to increased oxygen availability, while in other cases, there is nearly no decrease in P. The high reactivity of the interface complicates quantification of groundwater-borne P loads to the lake, making difficult clear differentiation of internal and external P loads to surface water. Anthropogenic sources of nitrogen (N) in groundwater are similar to those of phosphate. However, the environmental fate of N differs fundamentally from P because N occurs in several different redox states, each with different mobility. While nitrate behaves essentially conservatively in most oxic aquifers, ammonium's mobility is similar to that of phosphate. Nitrate may be transformed to gaseous N 2 in reducing conditions and permanently removed from the system. Biogeochemical turnover of N is common at the reactive aquifer-lake interface. Nutrient loads from LGD were compiled from the literature. Groundwater-borne P loads vary from 0.74 to 2900 mg PO 4 -P m −2 year −1 ; for N, these loads vary from 0.001 to 640 g m −2 year −1 . Even small amounts of seepage can carry large nutrient loads due to often high nutrient concentrations in groundwater. Large spatial heterogeneity, uncertain areal extent of the interface and difficult accessibility make every determination of LGD a challenge. However, determinations of LGD are essential to effective lake management.

Hydrological Processes↗

Groundwater: The disregarded component in lake water and nutrient budgets, Part 1: Effects of groundwater on hydrology

Lake eutrophication is a large and growing problem in many parts of the world, commonly due to anthropogenic sources of nutrients. Improved quantification of nutrient inputs is required to address this problem, including better determination of exchanges between groundwater and lakes. This first of a two-part review provides a brief history of the evolution of the study of groundwater exchange with lakes, followed by a listing of the most commonly used methods for quantifying this exchange. Rates of exchange between lakes and groundwater compiled from the literature are statistically summarized for both exfiltration (flow from groundwater to a lake) and infiltration (flow from a lake to groundwater), including per cent contribution of groundwater to lake-water budgets. Reported rates of exchange between groundwater and lakes span more than five orders of magnitude. Median exfiltration is 0.74 cm/day, and median infiltration is 0.60 cm/day. Exfiltration ranges from near 0% to 94% of input terms in lake-water budgets, and infiltration ranges from near 0% to 91% of loss terms. Median values for exfiltration and infiltration as percentages of input and loss terms of lake-water budgets are 25% and 35%, respectively. Quantification of the groundwater term is somewhat method dependent, indicating that calculating the groundwater component with multiple methods can provide a better understanding of the accuracy of estimates. The importance of exfiltration to a lake budget ranges widely for lakes less than about 100 ha in area but generally decreases with increasing lake area, particularly for lakes that exceed 100 ha in area. No such relation is evident for lakes where infiltration occurs, perhaps because of the smaller sample size. Copyright © 2014 John Wiley & Sons, Ltd.

Hydrological Processes↗

Variably-saturated groundwater modeling for optimizing managed aquifer recharge using trench infiltration

Spreading-basin methods have resulted in more than 130 million cubic meters of recharge to the unconfined Navajo Sandstone of southern Utah in the past decade, but infiltration rates have slowed in recent years because of reduced hydraulic gradients and clogging. Trench infiltration is a promising alternative technique for increasing recharge and minimizing evaporation. This paper uses a variably saturated flow model to further investigate the relative importance of the following variables on rates of trench infiltration to unconfined aquifers: saturated hydraulic conductivity, trench spacing and dimensions, initial water-table depth, alternate wet/dry periods, and number of parallel trenches. Modeling results showed (1) increased infiltration with higher hydraulic conductivity, deeper initial water tables, and larger spacing between parallel trenches, (2) deeper or wider trenches do not substantially increase infiltration, (3) alternating wet/dry periods result in less overall infiltration than keeping the trenches continuously full, and (4) larger numbers of parallel trenches within a fixed area increases infiltration but with a diminishing effect as trench spacing becomes tighter. An empirical equation for estimating expected trench infiltration rates as a function of hydraulic conductivity and initial water-table depth was derived and can be used for evaluating feasibility of trench infiltration in other hydrogeologic settings

Utah↗

Effects of urbanization and stormwater control measures on streamflows in the vicinity of Clarksburg, Maryland, USA

Understanding the efficacy of revised watershed management methods is important to mitigating the impacts of urbanization on streamflow. We evaluated the influence of land use change, primarily as urbanization, and stormwater control measures on the relationship between precipitation and stream discharge over an 8-year period for five catchments near Clarksburg, Montgomery County, Maryland, USA. A unit-hydrograph model based on a temporal transfer function was employed to account for and standardize temporal variation in rainfall pattern, and properly apportion rainfall to streamflow at different time lags. From these lagged relationships, we quantified a correction to the precipitation time series to achieve a hydrograph that showed good agreement between precipitation and discharge records. Positive corrections appeared to include precipitation events that were of limited areal extent and therefore not captured by our rain gages. Negative corrections were analysed for potential causal relationships. We used mixed-model statistical techniques to isolate different sources of variance as drivers that mediate the rainfall–runoff dynamic before and after management. Seasonal periodicity mediated rainfall–runoff relationships, and land uses (i.e. agriculture, natural lands, wetlands and stormwater control measures) were statistically significant predictors of precipitation apportionment to stream discharge. Our approach is one way to evaluate actual effectiveness of management efforts in the face of complicating circumstances and could be paired with cost data to understand economic efficiency or life cycle aspects of watershed management. Published 2015. This article is a U.S. Government work and is in the public domain in the USA.

Maryland↗

Rethinking the longitudinal stream temperature paradigm: region-wide comparison of thermal infrared imagery reveals unexpected complexity of river temperatures

Prevailing theory suggests that stream temperature warms asymptotically in a downstream direction, beginning at the temperature of the source in the headwaters and leveling off downstream as it converges to match meteorological conditions. However, there have been few empirical examples of longitudinal patterns of temperature in large rivers due to a paucity of data. We constructed longitudinal thermal profiles (temperature versus distance) for 53 rivers in the Pacific Northwest (USA) using an extensive dataset of remotely sensed summertime river temperatures and classified each profile into one of five patterns of downstream warming: asymptotic (increasing then flattening), linear (increasing steadily), uniform (not changing), parabolic (increasing then decreasing), or complex (not fitting other classes). We evaluated (1) how frequently profiles warmed asymptotically downstream as expected, and (2) whether relationships between river temperature and common hydroclimatic variables differed by profile class. We found considerable diversity in profile shape, with 47% of rivers warming asymptotically, and 53% having alternative profile shapes. Water temperature did not warm substantially over the course of the river for coastal parabolic and uniform profiles, and for some linear and complex profiles. Profile classes showed no clear geographical trends. The degree of correlation between river temperature and hydroclimatic variables differed among profile classes, but there was overlap among classes. Water temperature in rivers with asymptotic or parabolic profiles was positively correlated with August air temperature, tributary temperature and velocity, and negatively correlated with elevation, August precipitation, gradient, and distance upstream. Conversely, associations were less apparent in rivers with linear, uniform, or complex profiles. Factors contributing to the unique shape of parabolic profiles differed for coastal and inland rivers, where downstream cooling was influenced locally by climate or cool water inputs, respectively. Potential drivers of shape for complex profiles were specific to each river. These thermal patterns indicate diverse thermal habitats that may promote resilience of aquatic biota to climate change. Without this spatial context, climate change models may incorrectly estimate loss of thermally suitable habitat.

California, Idaho, Oregon, Washington↗

Climate-water quality relationships in Texas reservoirs

Water temperature, dissolved oxygen, and concentrations of salts in surface water bodies can be affected by the natural environment, local human activities such as surface and ground water withdrawals, land use, and energy extraction, and variability and long-term trends in atmospheric conditions including temperature and precipitation. Here, we quantify the relationship between 121 indicators of mean and extreme temperature and precipitation and 24 water quality parameters in 57 Texas reservoirs using observational data records covering the period 1960 to 2010. We find that water temperature, dissolved oxygen, pH, specific conductance, chloride, sulfate, and phosphorus all show consistent correlations with atmospheric predictors, including high and low temperature extremes, dry days, heavy precipitation events, and mean temperature and precipitation over time scales ranging from one week to two years. Based on this analysis and published future projections for this region, we expect climate change to increase water temperatures, decrease dissolved oxygen levels, decrease pH, increase specific conductance, and increase levels of sulfate, chloride in Texas reservoirs. Over decadal time scales, this may affect aquatic ecosystems in the reservoirs, including altering the risk of conditions conducive to algae occurrence, as well as affecting the quality of water available for human consumption and recreation.

Texas↗

Wide-area estimates of evapotranspiration by red gum ( Eucalyptus camaldulensis ) and associated vegetation in the Murray-Darling River Basin, Australia

Floodplain red gum forests ( Eucalyptus camaldulensis plus associated grasses, reeds and sedges) are sites of high biodiversity in otherwise arid regions of southeastern Australia. They depend on periodic floods from rivers, but dams and diversions have reduced flood frequencies and volumes, leading to deterioration of trees and associated biota. There is a need to determine their water requirements so environmental flows can be administered to maintain or restore the forests. Their water requirements include the frequency and extent of overbank flooding, which recharges the floodplain soils with water, as well as the actual amount of water consumed in evapotranspiration (ET). We estimated the flooding requirements and ET for a 38 134 ha area of red gum forest fed by the Murrumbidgee River in Yanga National Park, New South Wales. ET was estimated by three methods: sap flux sensors placed in individual trees; a remote sensing method based on the Enhanced Vegetation Index from MODIS satellite imagery and a water balance method based on differences between river flows into and out of the forest. The methods gave comparable estimates yet covered different spatial and temporal scales. We estimated flood frequency and volume requirements by comparing Normalized Difference Vegetation Index values from Landsat images with flood history from 1995 to 2014, which included both wet periods and dry periods. ET during wet years is about 50% of potential ET but is much less in dry years because of the trees' ability to control stomatal conductance. Based on our analyses plus other studies, red gum trees at this location require environmental flows of 2000 GL yr −1 every other year, with peak flows of 20 000 ML d −1 , to produce flooding sufficient to keep them in good condition. However, only about 120–200 GL yr −1 of river water is consumed in ET, with the remainder flowing out of the forest where it enters the Murray River system.

Murray-Darling River Basin↗

Water-quality response to a high-elevation wildfire in the Colorado Front Range

Water quality of the Big Thompson River in the Front Range of Colorado was studied for 2 years following a high‐elevation wildfire that started in October 2012 and burned 15% of the watershed. A combination of fixed‐interval sampling and continuous water‐quality monitors was used to examine the timing and magnitude of water‐quality changes caused by the wildfire. Prefire water quality was well characterized because the site has been monitored at least monthly since the early 2000s. Major ions and nitrate showed the largest changes in concentrations; major ion increases were greatest in the first postfire snowmelt period, but nitrate increases were greatest in the second snowmelt period. The delay in nitrate release until the second snowmelt season likely reflected a combination of factors including fire timing, hydrologic regime, and rates of nitrogen transformations. Despite the small size of the fire, annual yields of dissolved constituents from the watershed increased 20–52% in the first 2 years following the fire. Turbidity data from the continuous sensor indicated high‐intensity summer rain storms had a much greater effect on sediment transport compared to snowmelt. High‐frequency sensor data also revealed that weekly sampling missed the concentration peak during snowmelt and short‐duration spikes during rain events, underscoring the challenge of characterizing postfire water‐quality response with fixed‐interval sampling.

Colorado↗

Effect of permafrost thaw on the dynamics of lakes recharged by ice-jam floods: case study in Yukon Flats, Alaska

Large river floods are a key water source for many lakes in fluvial periglacial settings. Where permeable sediments occur, the distribution of permafrost may play an important role in the routing of floodwaters across a floodplain. This relationship is explored for lakes in the discontinuous permafrost of Yukon Flats, interior Alaska, using an analysis that integrates satellite-derived gradients in water surface elevation, knowledge of hydrogeology, and hydrologic modeling. We observed gradients in water surface elevation between neighboring lakes ranging from 0.001 to 0.004. These high gradients, despite a ubiquitous layer of continuous shallow gravel across the flats, are consistent with limited groundwater flow across lake basins resulting from the presence of permafrost. Permafrost impedes the propagation of floodwaters in the shallow subsurface and constrains transmission to “fill-and-spill” over topographic depressions (surface sills), as we observed for the Twelvemile-Buddy Lake pair following a May 2013 ice-jam flood on the Yukon River. Model results indicate that permafrost table deepening of 1–11 m in gravel, depending on watershed geometry and subsurface properties, could shift important routing of floodwater to lakes from overland flow (fill-and-spill) to shallow groundwater flow (“fill-and-seep”). Such a shift is possible in the next several hundred years of ground surface warming, and may bring about more synchronous water level changes between neighboring lakes following large flood events. This relationship offers a potentially useful tool, well-suited to remote sensing, for identifying long-term changes in shallow groundwater flow resulting from thawing of permafrost.

Alaska↗

Urban base flow with low impact development

A novel form of urbanization, low impact development (LID), aims to engineer systems that replicate natural hydrologic functioning, in part by infiltrating stormwater close to the impervious surfaces that generate it. We sought to statistically evaluate changes in a base flow regime because of urbanization with LID, specifically changes in base flow magnitude, seasonality, and rate of change. We used a case study watershed in Clarksburg, Maryland, in which streamflow was monitored during whole-watershed urbanization from forest and agricultural to suburban residential development using LID. The 1.11-km 2 watershed contains 73 infiltration-focused stormwater facilities, including bioretention facilities, dry wells, and dry swales. We examined annual and monthly flow during and after urbanization (2004–2014) and compared alterations to nearby forested and urban control watersheds. We show that total streamflow and base flow increased in the LID watershed during urbanization as compared with control watersheds. The LID watershed had more gradual storm recessions after urbanization and attenuated seasonality in base flow. These flow regime changes may be because of a reduction in evapotranspiration because of the overall decrease in vegetative cover with urbanization and the increase in point sources of recharge. Precipitation that may once have infiltrated soil, been stored in soil moisture to be eventually transpired in a forested landscape, may now be recharged and become base flow. The transfer of evapotranspiration to base flow is an unintended consequence to the water balance of LID.

Hydrological Processes↗

Synthesis of soil-hydraulic properties and infiltration timescales in wildfire-affected soils

We collected soil-hydraulic property data from the literature for wildfire-affected soils, ash, and unburned soils. These data were used to calculate metrics and timescales of hydrologic response related to infiltration and surface runoff generation. Sorptivity ( S ) and wetting front potential (Ψ f ) were significantly different (lower) in burned soils compared with unburned soils, whereas field-saturated hydraulic conductivity ( K fs ) was not significantly different. The magnitude and duration of the influence of capillarity during infiltration was greatly reduced in burned soils, causing faster ponding times in response to rainfall. Ash had large values of S and K fs but moderate values of Ψ f , compared with unburned and burned soils, indicating ash has long ponding times in response to rainfall. The ratio of S 2 /K fs was nearly constant (~100 mm) for unburned soils but more variable in burned soils, suggesting that unburned soils have a balance between gravity and capillarity contributions to infiltration that may depend on soil organic matter, whereas in burned soils the gravity contribution to infiltration is greater. Changes in S and K fs in burned soils act synergistically to reduce infiltration and accelerate and amplify surface runoff generation. Synthesis of these findings identifies three key areas for future research. First, short timescales of capillary influences on infiltration indicate the need for better measurements of infiltration at times less than 1 min to accurately characterize S in burned soils. Second, using parameter values, such as Ψ f , from unburned areas could produce substantial errors in hydrologic modeling when used without adjustment for wildfire effects, causing parameter compensation and resulting underestimation of K fs . Third, more thorough measurement campaigns that capture soil-structural changes, organic matter impacts, quantitative water repellency trends, and soil-water content along with soil-hydraulic properties could drive the development of better techniques for numerically simulating infiltration in burned areas.

Hydrological Processes↗

Spatial and temporal patterns of dissolved organic matter quantity and quality in the Mississippi River Basin, 1997–2013

Recent studies have found insignificant or decreasing trends in time-series dissolved organic carbon (DOC) datasets, questioning the assumption that long-term DOC concentrations in surface waters are increasing in response to anthropogenic forcing, including climate change, land use, and atmospheric acid deposition. We used the weighted regressions on time, discharge, and season (WRTDS) model to estimate annual flow-normalized concentrations and fluxes to determine if changes in DOC quantity and quality signal anthropogenic forcing at 10 locations in the Mississippi River Basin. Despite increases in agriculture and urban development throughout the basin, net increases in DOC concentration and flux were significant at only 3 of 10 sites from 1997 to 2013 and ranged between −3.5% to +18% and −0.1 to 19%, respectively. Positive shifts in DOC quality, characterized by increasing specific ultraviolet absorbance at 254 nm, ranged between +8% and +45%, but only occurred at one of the sites with significant DOC quantity increases. Basinwide reductions in atmospheric sulfate deposition did not result in large increases in DOC either, likely because of the high buffering capacity of the soil. Hydroclimatic factors including annual discharge, precipitation, and temperature did not significantly change during the 17-year timespan of this study, which contrasts with results from previous studies showing significant increases in precipitation and discharge over a century time scale. Our study also contrasts with those from smaller catchments, which have shown stronger DOC responses to climate, land use, and acidic deposition. This temporal and spatial analysis indicated that there was a potential change in DOC sources in the Mississippi River Basin between 1997 and 2013. However, the overall magnitude of DOC trends was not large, and the pattern in quantity and quality increases for the 10 study sites was not consistent throughout the basin.

Hydrological Processes↗

Comparing catchment hydrologic response to a regional storm using specific conductivity sensors

A better understanding of stormwater generation and solute sources is needed to improve the protection of aquatic ecosystems, infrastructure, and human health from large runoff events. Much of our understanding of water and solutes produced during stormflow comes from studies of individual, small headwater catchments. This study compared many different types of catchments during a single large event to help isolate landscape controls on streamwater and solute generation, including human-impacted land cover. We used a distributed network of specific electrical conductivity sensors to trace storm response during the post-tropical cyclone Sandy event of October 2012 at 29 catchments across the state of New Hampshire. A citizen science sensor network, Lotic Volunteer for Temperature, Electrical Conductivity, and Stage, provided a unique opportunity to investigate high-temporal resolution stream behavior at a broad spatial scale. Three storm response metrics were analyzed in this study: (a) fraction of new water contributing to the hydrograph; (b) presence of first flush (mobilization of solutes during the beginning of the rain event); and (c) magnitude of first flush. We compared new water and first flush to 64 predictor attributes related to land cover, soil, topography, and precipitation. The new water fraction was positively correlated with low and medium intensity development in the catchment and riparian buffers and with the precipitation from a rain event 9 days prior to Sandy. The presence of first flush was most closely related (positively) to soil organic matter. Magnitude of first flush was not strongly related to any of the catchment variables. Our results highlight the potentially important role of human landscape modification in runoff generation at multiple spatial scales and the lack of a clear role in solute flushing. Further development of regional-scale in situ sensor networks will provide better understanding of stormflow and solute generation across a wide range of landscape conditions.

New Hampshire↗

Tracer-based evidence of heterogeneity in subsurface flow and storage within a boreal hillslope

Runoff from boreal hillslopes is often affected by distinct soil boundaries, including the frozen boundary and the organic – mineral boundary (OMB), where highly porous and hydraulically-conductive organic material overlies fine-grained mineral soils. Viewed from the surface, ground cover appears as a patchwork on sub-meter scales, with thick, moss mats interspersed with lichen-covered, silty soils with gravel inclusions. We conducted a decameter-scale subsurface tracer test on a boreal forest hillslope in interior Alaska to quantify locations and mechanisms of transport and storage in these soils, focusing on the OMB. A sodium bromide tracer was added as a slug addition to a pit and sampled at 40 down-gradient wells, screened primarily at the OMB and within a 7 by 12 m well field. We maintained an elevated head in the injection pit for 8.5 h to simulate a storm. Tracer breakthrough velocities ranged from < 0.12 to 0.93 m hr-1, with the highest velocities in lichen-covered soils. After 12 hours and cessation of the elevated head, the tracer coalesced and was only detected in thick mosses at a trough in the OMB. By 24 hours, approximately 17% of the tracer mass could be accounted for. The majority of the mass loss occurred between 4 and 12 hours, while the tracer was in contact with lichen-covered soils, which is consistent with tracer transport into deeper flow paths via preferential flow through discrete gravelly areas. Slow breakthroughs suggest that storage and exchange also occurred in shallow soils, likely related to saturation and drainage in fine-grained mineral soils caused by the elevated hydraulic head. These findings highlight the complex nature of storage and transmission of water and solutes from boreal hillslopes to streams, and are particularly relevant given rapid changes to boreal environments related to climate change, thawing permafrost and increasing fire severity.

Hydrological Processes↗

Heat as a groundwater tracer in shallow and deep heterogeneous media: Analytical solution, spreadsheet tool, and field applications

Groundwater flow advects heat, and thus, the deviation of subsurface temperatures from an expected conduction‐dominated regime can be analysed to estimate vertical water fluxes. A number of analytical approaches have been proposed for using heat as a groundwater tracer, and these have typically assumed a homogeneous medium. However, heterogeneous thermal properties are ubiquitous in subsurface environments, both at the scale of geologic strata and at finer scales in streambeds. Herein, we apply the analytical solution of Shan and Bodvarsson (2004), developed for estimating vertical water fluxes in layered systems, in 2 new environments distinct from previous vadose zone applications. The utility of the solution for studying groundwater‐surface water exchange is demonstrated using temperature data collected from an upwelling streambed with sediment layers, and a simple sensitivity analysis using these data indicates the solution is relatively robust. Also, a deeper temperature profile recorded in a borehole in South Australia is analysed to estimate deeper water fluxes. The analytical solution is able to match observed thermal gradients, including the change in slope at sediment interfaces. Results indicate that not accounting for layering can yield errors in the magnitude and even direction of the inferred Darcy fluxes. A simple automated spreadsheet tool (Flux‐LM) is presented to allow users to input temperature and layer data and solve the inverse problem to estimate groundwater flux rates from shallow (e.g., <1 m) or deep (e.g., up to 100 m) profiles. The solution is not transient, and thus, it should be cautiously applied where diel signals propagate or in deeper zones where multi‐decadal surface signals have disturbed subsurface thermal regimes.

Hydrological Processes↗