Search

The recent boom of commercially available drone technology capable of supporting geophysical imaging systems has led to several practitioners adopting the platforms to collect geophysical data. This has led to significant opportunities for studying hydrologic systems, as these platforms enable more efficient data collection and data collection in difficult-to-access locations. Recent examples from snow and stream studies show the viability of these platforms for studying complex systems in difficult locations and the added value of data provided by the technology. Small unoccupied aerial system (sUAS) ground-penetrating radar (GPR) surveys conducted over water bodies and snow have proven beneficial for determining bathymetry in dangerous swift water and depth of snowpack in mountainous regions. Challenges arise due to decoupling of the GPR antenna and the ground surface, which limits signal penetration. Additionally, multifrequency or multichannel GPR imaging remains a challenge due to the payload limitations of an sUAS. Drone-mounted thermal infrared data collection enables the generation of thermal orthomosaic images with unprecedented detail and spatial range for quantifying groundwater fluxes to surface water bodies and reduces issues with infrared reflection common in ground-based studies. Several drone-mounted electromagnetic studies also show utility in efficient, high-resolution mapping of near-surface soil conductivity. Overall, the use of drone platforms for studying hydrologic systems remains a frontier in science applications that needs further refinement before wide-scale adoption. The preliminary studies presented here hint that the development of robust geophysical data acquisition tools based on UAS will lead to an explosion of applications to near-surface hydrologic problems.

Groundwater/surface-water (GW/SW) exchange and hyporheic processes are topics receiving increasing attention from the hydrologic community. Hydraulic, chemical, temperature, geophysical, and remote sensing methods are used to achieve various goals (e.g., inference of GW/SW exchange, mapping of bed materials, etc.), but the application of these methods is constrained by site conditions such as water depth, specific conductance, bed material, and other factors. Researchers and environmental professionals working on GW/SW problems come from diverse fields and rarely have expertise in all available field methods; hence there is a need for guidance to design field campaigns and select methods that both contribute to study goals and are likely to work under site-specific conditions. Here, we present the spreadsheet-based GW/SW-Method Selection Tool (GW/SW-MST) to help practitioners identify methods for use in GW/SW and hyporheic studies. The GW/SW-MST is a Microsoft Excel-based decision support tool in which the user selects answers to questions about GW/SW-related study goals and site parameters and characteristics. Based on user input, the tool indicates which methods from a toolbox of 32 methods could potentially contribute to achieving the specified goals at the site described.

Groundwater discharge to rivers takes many forms, including preferential groundwater discharge points (PDPs) along riverbanks that are exposed at low flows, with multi-scale impacts on aquatic habitat and water quality. The physical controls on the spatial distribution of PDPs along riverbanks are not well-defined, rendering their prediction and representation in models challenging. To investigate the local riverbank sediment controls on PDP occurrence, we tested drone-based and handheld thermal infrared to efficiently map PDP locations along two mainstem rivers. Early in the study, we found drone imaging was better suited to locating tributary and stormwater inflows, which created relatively large water surface thermal anomalies in winter, compared to PDPs that often occurred at the sub-meter scale and beneath riparian tree canopy. Therefore, we primarily used handheld thermal infrared imaging from watercraft to map PDPs and larger seepage faces along 12-km of the fifth-order Housatonic River in Massachusetts, USA and 26-km of the Farmington River in Connecticut, USA. Overall, we mapped 31 riverbank PDPs along the Housatonic reach that meanders through lower permeability soils, and 104 PDPs along the Farmington reach that cuts through sandier sediments. Riverbank soil parameters extracted at PDP locations from the Soil Survey Geographic (SSURGO) database did not differ substantially from average bank soils along either reach, although the Farmington riverbank soils were on average 5× more permeable than Housatonic riverbank soils, likely contributing to the higher observed prevalence of PDPs. Dissolved oxygen measured in discharge water at these same PDPs varied widely, but showed no relation to measured sand, clay, or organic matter content in surficial soils indicating a lack of substantial near-surface aerobic reaction. The PDP locations were investigated for the presence of secondary bank structures, and commonly co-occurred with riparian tree root masses indicating the importance of localized physical controls on the spatial distribution of riverbank PDPs.

Exchanges of groundwater and surface-water are fundamental to a wide range of water-supply and water-quality management issues but challenging to map beyond the reach scale. Waterborne gradient self-potential (SP) measurements are directly sensitive to water flow through riverbed sediments and can be used to infer exchange locations, direction (gain versus loss), scale, and relative changes, but to date applications to river corridor hydrology are limited. Numerical modeling and field experiments were therefore performed herein, each emphasizing waterborne gradient SP logging for identifying and locating focused vertical groundwater discharge (surface-water gain) and recharge (surface-water loss) in a river. Two and three-dimensional numerical models were constructed to simulate the polarities, appearances, and peak amplitudes of streaming-potential and electric-field anomalies on a riverbed and in the surface-water that were attributable to steady-state vertical fluxes of groundwater through high-permeability conduits in the riverbed. Effects of varied hydraulic length-scale of exchange and surface-water depth were tested through numerical modeling. Modeling results aided in data acquisition and interpretation for three repeated field experiments performed along a 1.5–2.0 km reach of the Quashnet River in Cape Cod, Massachusetts, where focused, meter-scale groundwater discharges occur at discrete locations within otherwise ubiquitous and more diffuse groundwater upwelling conditions. Strong gradient SP anomalies were repeatedly measured in the Quashnet River at previously confirmed locations of focused groundwater discharge, showing the efficacy of waterborne gradient SP logging in identifying and characterizing groundwater/surface water exchange dynamics at multiple river network scales.

Groundwater discharge generates streamflow and influences stream thermal regimes. However, the water quality and thermal buffering capacity of groundwater depends on the aquifer source-depth. Here, we pair multi-year air and stream temperature signals to categorize 1729 sites across the continental United States as having major dam influence, shallow or deep groundwater signatures, or lack of pronounced groundwater (atmospheric) signatures. Approximately 40% of non-dam stream sites have substantial groundwater contributions as indicated by characteristic paired air and stream temperature signal metrics. Streams with shallow groundwater signatures account for half of all groundwater signature sites and show reduced baseflow and a higher proportion of warming trends compared to sites with deep groundwater signatures. These findings align with theory that shallow groundwater is more vulnerable to temperature increase and depletion. Streams with atmospheric signatures tend to drain watersheds with low slope and greater human disturbance, indicating reduced stream-groundwater connectivity in populated valley settings., Groundwater discharge generates streamflow and influences stream thermal regimes. Classifying more than 1700 streams across the US by using an empirically-based approach the study shows that the vulnerability of streams to stressors depends on the aquifer source-depth of groundwater discharge

The sediment–water interfaces (SWI) of streams serve as important biogeochemical hotspots in watersheds and contribute to whole-catchment reactive nitrogen budgets and water-quality conditions. Recently, the SWI has been identified as an important source of nitrous oxide (N2O) produced in streams, with SWI residence time among the principal controls on its production. Here, we conducted a series of controlled manipulations of SWI exchange in an urban stream that has high dissolved N2O concentrations and where we concurrently evaluated less-mobile porosity dynamics. Our experiments took place within isolated portions of two sediment types: a coarse sandy stream bed resulting from excess road-sand application in the watershed, and a coarse sand mixed with clay and organic particles. In these manipulation experiments we systematically varied SWI vertical-flux rates and residence times to evaluate their effect on the fate of nitrate and production rates of N2O. Our experiments demonstrate that the fate and transport of nitrate and N2O production are influenced by hydrologic flux rates through SWI sediments and associated residence times. Specifically, we show that manipulations of hydrologic flux systematically shifted the depth of the bulk oxic–anoxic interface in the sediments, and that nitrate removal increased with residence time. Our results also support the emerging hypothesis of a ‘Goldilocks’ timescale for the production of nitrous oxide, when transport and reaction timescales favor incomplete denitrification. Areal N2O production rates were up to threefold higher during an intermediate residence-time experiment, compared to shorter or longer residence times. In our companion study we documented that the studied sediments were dominated by a long-residence-time less-mobile porosity domain, which could explain why we observed N2O production even in bulk-oxic sediments. Overall, we have experimentally demonstrated that changes to SWI hydrologic residence times and SWI substrate associated with urbanization can change the biogeochemical function of the river corridor.

Manganese (Mn) plays a critical role in river-water quality because Mn-oxides serve as sorption sites for contaminant metals. The aim of this study is to understand the seasonal cycling of Mn in an alpine streambed that experiences large spring snowmelt events and the potential responses to changes in snowmelt timing and magnitude. To address this goal, annual variations in river-water/groundwater interaction and Mn(aq) transport were measured and modeled in the bed of East River, Colorado, USA. In observations and numerical models, oxygenated river water containing dissolved organic carbon (DOC) mixes with groundwater rich in Mn(aq) in the streambed. The mixing depth increases during spring snowmelt when river discharge increases, leading to a greater DOC supply to the hyporheic zone and net respiration of Mn-oxides, despite an enhanced supply of oxygen. As groundwater upwelling resumes during the subsequent baseflow period, Mn(aq)-rich groundwater mixes with oxygenated river water, resulting in net accumulation of Mn-oxides until the bed freezes in winter. To explore potential responses of Mn transport to different climate-induced hydrological regimes, three hydrograph scenarios were numerically modeled (historic, low-snow, and storm) for the Rocky Mountain region. In a warming climate, Mn(aq) export to the river decreases, and Mn(aq) oxidation is favored in the upper streambed sediments over more of the year. One important implication is that the streambed may have an increased sorption capacity for metals over more of the year, leading to potential changes in river-water quality.

In mountain headwater streams, the quality and resilience of summer cold-water habitat is generally regulated by stream discharge, longitudinal stream channel connectivity and groundwater exchange. These critical hydrologic processes are thought to be influenced by the stream corridor bedrock contact depth (sediment thickness), a parameter often inferred from sparse hillslope borehole information, piezometer refusal and remotely sensed data. To investigate how local bedrock depth might control summer stream temperature and channel disconnection (dewatering) patterns, we measured stream corridor bedrock depth by collecting and interpreting 191 passive seismic datasets along eight headwater streams in Shenandoah National Park (Virginia, USA). In addition, we used multi-year stream temperature and streamflow records to calculate several baseflow-related metrics along and among the study streams. Finally, comprehensive visual surveys of stream channel dewatering were conducted in 2016, 2019 and 2021 during summer low flow conditions (124 total km of stream length). We found that measured bedrock depths along the study streams were not well-characterized by soils maps or an existing global-scale geologic dataset where the latter overpredicted measured depths by 12.2 m (mean) or approximately four times the average bedrock depth of 2.9 m. Half of the eight study stream corridors had an average bedrock depth of less than 2 m. Of the eight study streams, Staunton River had the deepest average bedrock depth (3.4 m), the coldest summer temperature profiles and substantially higher summer baseflow indices compared to the other study steams. Staunton River also exhibited paired air and water annual temperature signals suggesting deeper groundwater influence, and the stream channel did not dewater in lower sections during any baseflow survey. In contrast, Paine Run and Piney River did show pronounced, patchy channel dewatering, with Paine Run having dozens of discrete dry channel sections ranging from 1 to greater than 300 m in length. Stream dewatering patterns were apparently influenced by a combination of discrete deep bedrock (20+ m) features and more subtle sediment thickness variation (1–4 m) depending on local stream valley hydrogeology. In combination, these unique datasets show the first large-scale empirical support for existing conceptual models of headwater stream disconnection based on spatially variable underflow capacity and shallow groundwater supply.

Spatially preferential flow processes occur at nested scales at the sediment-water interface (SWI), due in part to sediment heterogeneities, which may be enhanced in flashy urban streams with heavy road sand influence. However, several factors, including the flow-rate dependence of preferential hyporheic flow and discrete groundwater discharge zones are commonly overlooked in reach-scale models of groundwater/surface water exchange. Using a series of controlled-head tracer-injection experiments coupled with cm-scale geophysics within the highly reactive upper 30 cm of the hyporheic zone of an urban stream, we quantified the flow dependence of local less-mobile porosity volume, mass-transfer rate coefficient, and the resulting local residence time in the less-mobile pore space at three controlled downward fluid fluxes (0.8, 2, and 3 m/d). Experiments were performed in two adjacent streambed locations, representing different sediment bulk vertical permeability. Less-mobile porosity parameters were generally substantial and similar between the two streambed locations; though a more competent, thin, organic layer at ∼15 cm depth in one location strongly impacted tracer loading, flushing dynamics, and local residence times. Increased downward flux led to (1) a decrease in less-mobile porosity residence time in all experiments, and (2) an increase in less-mobile porosity fraction for most experiments. Additionally, at the larger stream reach-scale, surface electrodes for electrical resistivity measurement were installed along 22 m of the wetted stream channel. These surface electrode measurements were collected during a natural storm flow event, which revealed widespread, short-term, flushing (e.g. 60 h) flushing of the SWI with riparian zone groundwater. Flow dependence of preferential hyporheic zone flowpaths, like in the controlled tracer experiments, was also observed in these reach-scale electrical resistivity tomography measurements. Our findings reveal that the spatial and temporal dependence of preferential flow processes create highly dynamic SWI conditions that will affect the physical and coupled biogeochemical functions of the SWI in urbanized, sand-impacted streams.

Quantitative evaluation of groundwater/surface water exchange dynamics is universally challenging in large river systems, because existing methodology often does not yield spatially-distributed data and is difficult to apply in deeper water. Here we apply a combined near-surface geophysical and direct groundwater chemical toolkit to refine fresh groundwater discharge estimates to the Colorado River through a 4-km2 wetland that borders the town of Moab, Utah, USA. Preliminary characterization of raw electromagnetic imaging (EMI) data, collected by kayak and by walking, was used to guide additional direct-contact electrical measurements and installation of new monitoring wells. Chemical data from the wells strongly supported the EMI spatial characterization of preferential fresh groundwater discharge embedded in natural brine groundwaters and weighted to the southern wetland section. Inversion of the EMI data revealed sub-meter scale detail regarding bulk electrical conductivity zonation across approximately 15.5 km of transects, collected in only 3 days. This electrical detail indicates processes such as salinization of the unsaturated zone and direct discharge through the Colorado River sediments and a tributary creek bed. Overall, the study contributed to a substantial reduction in fresh groundwater discharge estimates previously made using sparse existing well data and a simplified assumption of diffuse fresh groundwater discharge below the entire wetland. EMI will likely become a widely used tool in systems with natural electrical contrast as groundwater/surface water hydrogeologists continue to recognize the prevalence of preferential groundwater discharge processes.

Heat has been widely applied to trace groundwater-surface water exchanges in inland environments, but it is infrequently applied in coastal sediment where head oscillations induce periodicity in water flux magnitude/direction and heat advection. This complicates interpretation of temperatures to estimate water fluxes. We investigate the convolution of thermal and hydraulic signals to assess the viability of using heat as a tracer in environments with tidal head oscillations superimposed on submarine groundwater discharge. We first generate sediment temperature and head time series for conditions ranging from no tide to mega-tidal using a numerical model (SUTRA) forced with periodic temperature and tidal head signals. We then analyze these synthetic temperature time series using heat tracing software (VFLUX2 and 1DTempPro) to evaluate if conventional terrestrial approaches to infer fluxes from temperatures are applicable for coastal settings. We consider high-frequency water flux variability within a tidal signal and averaged over tidal signals. Results show that VFLUX2 analytical methods reasonably estimated the mean discharge fluxes in most cases but could not reproduce the flux variability within tidal cycles. The model results further reveal that high-frequency time series of water fluxes varying in magnitude and direction can be accurately estimated if paired temperatures and hydraulic heads are analyzed using numerical models (e.g., 1DTempPro) that consider both dynamic hydraulic gradients and thermal signals. These results point to the opportunity to incorporate pressure sensors within heat tracing instrumentation to better assess sub-daily flux oscillations and associated reactive processes.

Releases of oil and gas (OG) wastewaters can have complex effects on stream-water quality and downstream organisms, due to sediment-water interactions and groundwater/surface water exchange. Previously, elevated concentrations of sodium (Na), chloride (Cl), barium (Ba), strontium (Sr), and lithium (Li), and trace hydrocarbons were determined to be key markers of OG wastewater releases when combined with Sr and radium (Ra) isotopic compositions. Here, we assessed the persistence of an OG wastewater spill in a creek in North Dakota using a combination of geochemical measurements and modeling, hydrologic analysis, and geophysical investigations. OG wastewater comprised 0.1 to 0.3% of the stream-water compositions at downstream sites in February and June 2015 but could not be quantified in 2016 and 2017. However, OG-wastewater markers persisted in sediments and pore water for 2.5 years after the spill and up to 7.2-km downstream from the spill site. Concentrations of OG wastewater constituents were highly variable depending on the hydrologic conditions. Electromagnetic measurements indicated substantially higher electrical conductivity under the bank adjacent to a seep 7.2 km downstream from the spill site. Geomorphic investigations revealed mobilization of sediment is an important contaminant transport process. Labile Ba, Ra, Sr, and ammonium (NH4) concentrations extracted from sediments indicated sediments are a long-term reservoir of these constituents, both in the creek and on the floodplain. Using the drivers of ecological effects identified at this intensively studied site we identified 41 watersheds across the North Dakota landscape that may be subject to similar episodic inputs from OG wastewater spills. Effects of contaminants released to the environment during OG waste management activities remain poorly understood; however, analyses of Ra and Sr isotopic compositions, as well as trace inorganic and organic compound concentrations at these sites in pore-water provide insights into potentials for animal and human exposures well outside source-remediation zones.

Aufeis are sheets of ice unique to cold regions that originate from repeated flooding and freezing events during the winter. They have hydrological importance associated with summer flows and winter insulation, but little is known about the seasonal dynamics of the unfrozen sediment layer beneath them. This layer may support perennial groundwater flow in regions with otherwise continuous permafrost. For this study, ground penetrating radar (GPR) were collected in September 2016 (maximum thaw) and April 2017 (maximum frozen) at the Kuparuk aufeis field on the North Slope of Alaska. Supporting surface nuclear magnetic resonance data were collected during the maximum frozen campaign. These point-in-time geophysical data sets were augmented by continuous subsurface temperature data and periodic Structure-from-Motion digital elevation models collected seasonally. GPR and difference digital elevation model data showed up to 6 m of ice over the sediment surface. Below the ice, GPR and nuclear magnetic resonance identified regions of permafrost and regions of seasonally frozen sediment (i.e., the active layer) underlain by a substantial lateral talik that reached >13-m thickness. The seasonally frozen cobble layer above the talik was typically 3- to 5-m thick, with freezing apparently enabled by relatively high thermal diffusivity of the overlying ice and rock cobbles. The large talik suggests that year-round groundwater flow and coupled heat transport occurs beneath much of the feature. Highly permeable alluvial material and discrete zones of apparent groundwater upwelling indicated by geophysical and ground temperature data allows direct connection between the aufeis and the talik below.

Fiber-optic distributed temperature sensing (FO-DTS) has proven to be a transformative technology for the hydrologic sciences, with application to diverse problems including hyporheic exchange, groundwater/surface-water interaction, fractured-rock characterization, and cold regions hydrology. FO-DTS produces large, complex, and information-rich datasets. Despite the potential of FO-DTS, adoption of the technology has been impeded by lack of tools for data processing, analysis, and visualization. New tools are needed to efficiently and fully capitalize on the information content of FO-DTS datasets. To this end, we present DTSGUI, a public-domain Python-based software package for editing, parsing, processing, statistical analysis, georeferencing, and visualization of FO-DTS data.

The hydrogeology below large surface water features such as rivers and estuaries is universally under-informed at the long reach to basin scales (tens of km+). This challenge inhibits the accurate modeling of fresh/saline groundwater interfaces and groundwater/surface water exchange patterns at management-relevant spatial extents. Here we introduce a towed, floating transient electromagnetic (TEM) system (i.e. FloaTEM) for rapid (up to 15 km/h) high resolution electrical mapping of the subsurface below large water bodies to depths often a factor of 10 greater than other towed instruments. The novel FloaTEM system is demonstrated at a range of diverse 4th through 6th-order riverine settings across the United States including 1) the Farmington River, near Hartford, Connecticut; 2) the Upper Delaware River near Barryville, New York; 3) the Tallahatchie River near Shellmound, Mississippi; and, 4) the Eel River estuary, on Cape Cod, near Falmouth, Massachusetts. Airborne frequency-domain electromagnetic and land-based towed TEM data are also compared at the Tallahatchie River site, and streambed geologic scenarios are explored with forward modeling. A range of geologic structures and pore water salinity interfaces were identified. Process-based interpretation of the case study data indicated FloaTEM can resolve varied sediment-water interface materials, such as the accumulation of fines at the bottom of a reservoir and permeable sand/gravel riverbed sediments that focus groundwater discharge. Bedrock layers were mapped at several sites, and aquifer confining units were defined at comparable resolution to airborne methods. Terrestrial fresh groundwater discharge with flowpaths extending hundreds of meters from shore was also imaged below the Eel River estuary, improving on previous hydrogeological characterizations of that nutrient-rich coastal exchange zone. In summary, the novel FloaTEM system fills a critical gap in our ability to characterize the hydrogeology below surface water features and will support more accurate prediction of groundwater/surface water exchange dynamics and fresh-saline groundwater interfaces.

Streams strongly influenced by groundwater discharge may serve as “climate refugia” for sensitive species in regions of increasingly marginal thermal conditions. The main goal of this study is to develop paired air and stream water annual temperature signal analysis techniques to elucidate the relative groundwater contribution to stream water and the effective groundwater flowpath depth. Groundwater discharge to streams attenuates surface water temperature signals, and this attenuation can be diagnostic of groundwater gaining systems. Additionally, discharge from shallow groundwater flowpaths can theoretically transfer lagged annual temperature signals from aquifer to stream water. Here we explore this concept using multi-year temperature records from 120 stream sites located across 18 mountain watersheds of Shenandoah National Park, VA, USA and a coastal watershed in Massachusetts, USA. Both areas constitute important cold-water habitat for native brook trout (Salvelinus fontinalis). Observed annual temperature signals indicate a dominance of shallow groundwater discharge to streams in the National Park, in contrast to the coastal watershed that has strong, apparently deeper, groundwater influence. The average phase lag from air to stream signals in Shenandoah National Park is 11 d; however, extended lags of approximately 1 month were observed in a subset of streams. In contrast, the coastal stream has pronounced attenuation of annual temperature signals without notable phase lag. To better understand these observed differences in signal characteristics, analytical and numerical models are used to quantify mixing of the annual temperature signals of surface and groundwater. Simulations using a total heat budget numerical model indicate groundwater-induced annual temperature signal phase lags are likely to show greater downstream propagation than the related signal amplitude attenuation. The measurement of multi-seasonal paired air and water temperatures offers great promise toward understanding catchment processes and informing current cold-water habitat management at ecologically-relevant scales.

Identifying and quantifying groundwater exchange is critical when considering contaminant fate and transport at the groundwater/surface-water interface. In this paper, areally distributed temperature and point seepage measurements are used to efficiently assess spatial and temporal groundwater discharge patterns through a glacial-kettle lakebed area containing a zero-valent iron permeable reactive barrier (PRB). Concern was that the PRB was becoming less permeable with time owing to biogeochemical processes within the PRB. Patterns of groundwater discharge over an 8-year period were examined using fiber-optic distributed temperature sensing (FO-DTS) and snapshot-in-time point measurements of temperature. The resulting thermal maps show complex and uneven distributions of temperatures across the lakebed and highlight zones of rapid seepage near the shoreline and along the outer boundaries of the PRB. Repeated thermal mapping indicates an increase in lakebed temperatures over time at periods of similar stage and surface-water temperature. Flux rates in six seepage meters permanently installed on the lakebed in the PRB area decreased on average by 0.021 md-1 (or about 4.5 percent) annually between 2004 and 2015. Modeling of diurnal temperature signals from shallow vertical profiles yielded mean flux values ranging from 0.39 to 1.15 md-1, with stronger fluxes generally related to colder lakebed temperatures. The combination of an increase in lakebed temperatures, declines in direct seepage, and observations of increased cementation of the lakebed surface provide in situ evidence that the permeability of the PRB is declining. The presence of temporally persistent rapid seepage zones is also discussed.

Groundwater/surface-water exchanges in streams are inexorably linked to adjacent aquifer dynamics. As surface-water temperatures continue to increase with climate warming, refugia created by groundwater connectivity is expected to enable cold water fish species to survive. The shallow alluvial aquifers that source groundwater seepage to headwater streams, however, may also be sensitive to seasonal and long-term air temperature dynamics. Depth to bedrock can directly influence shallow aquifer flow and thermal sensitivity, but is typically ill-defined along the stream corridor in steep mountain catchments. We employ rapid, cost-effective passive seismic measurements to evaluate the variable thickness of the shallow colluvial and alluvial aquifer sediments along a headwater stream supporting cold water-dependent brook trout ( Salvelinus fontinalis ) in Shenandoah National Park, VA, USA. Using a mean depth to bedrock of 2.6 m, numerical models predicted strong sensitivity of shallow aquifer temperature to the downward propagation of surface heat. The annual temperature dynamics (annual signal amplitude attenuation and phase shift) of potential seepage sourced from the shallow modeled aquifer were compared to several years of paired observed stream and air temperature records. Annual stream water temperature patterns were found to lag local air temperature by ∼8–19 d along the stream corridor, indicating that thermal exchange between the stream and shallow groundwater is spatially variable. Locations with greater annual signal phase lag were also associated with locally increased amplitude attenuation, further suggestion of year-round buffering of channel water temperature by groundwater seepage. Numerical models of shallow groundwater temperature that incorporate regional expected climate warming trends indicate that the summer cooling capacity of this groundwater seepage will be reduced over time, and lower-elevation stream sections may no longer serve as larger-scale climate refugia for cold water fish species, even with strong groundwater discharge.

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.

Dual-domain models are used to explain anomalous solute-transport behavior observed in diverse hydrologic settings and applications, from groundwater remediation to hyporheic exchange. To constrain such models, new methods are needed with sensitivity to both immobile and mobile domains. Recent experiments indicate that dual-domain transport of ionic tracers has an observable geoelectrical signature, appearing as a non-linear, hysteretic relation between paired bulk and fluid electrical conductivity. Here, we present a mechanistic explanation for this geoelectrical signature and evaluate assumptions underlying a previously published petrophysical model for bulk conductivity in dual-domain media. Pore network modeling of fluid flow, solute transport, and electrical conduction (1) verifies the geoelectrical signature of dual-domain transport, (2) reveals limitations of the previously used petrophysical model, and (3) demonstrates that a new petrophysical model, based on differential-effective media theory, closely approximates the simulated bulk/fluid conductivity relation. These findings underscore the potential of geophysically based calibration of dual-domain models.

Large-scale wetland restoration often focuses on repairing the hydrologic connections degraded by anthropogenic modifications. Of these hydrologic connections, groundwater discharge is an important target, as these surface water ecosystem control points are important for thermal stability, among other ecosystem services. However, evaluating the effectiveness of the restoration activities on establishing groundwater discharge connection is often difficult over large areas and inaccessible terrain. Unoccupied aircraft systems (UAS) are now routinely used for collecting aerial imagery and creating digital surface models (DSM). Lightweight thermal infrared (TIR) sensors provide another payload option for generation of sub-meter-resolution aerial TIR orthophotos. This technology allows for the rapid and safe survey of groundwater discharge areas. Aerial TIR water-surface data were collected in March 2019 at Tidmarsh Farms, a former commercial cranberry peatland located in coastal Massachusetts, USA (41&deg, 54&prime, 17&Prime, N 70&deg, 34&prime, W), where stream and wetland restoration actions were completed in 2016. Here, we present a 0.4 km2 georeferenced, temperature-calibrated TIR orthophoto of the area. The image represents a mosaic of nearly 900 TIR images captured by UAS in a single morning with a total flight time of 36 min and is supported by a DSM derived from UAS-visible imagery. The survey was conducted in winter to maximize temperature contrast between relatively warm groundwater and colder ambient surface environment, lower-density groundwater rises above cool surface waters and thus can be imaged by a UAS. The resulting TIR orthomosaic shows fine detail of seepage distribution and downstream influence along the several restored channel forms, which was an objective of the ecological restoration design. The restored stream channel has increased connectivity to peatland groundwater discharge, reducing the ecosystem thermal stressors. Such aerial techniques can be used to guide ecological restoration design and assess post-restoration outcomes, especially in settings where ecosystem structure and function is governed by groundwater and surface water interaction.

River to floodplain hydrologic connectivity is strongly enhanced by beaver- (Castor canadensis) engineered channel water diversions. The hydroecological impacts are wide ranging and generally positive, however, the hydrogeochemical characteristics of beaver-induced flowpaths have not been thoroughly examined. Using a suite of complementary ground- and drone-based heat tracing and remote sensing methodology we characterized the physical template of beaver-induced floodplain exchange for two alluvial mountain streams near Crested Butte, Colorado, USA. A flowpath-oriented perspective to water quality sampling allowed characterization of the chemical evolution of channel water diverted through floodplain beaver ponds and ultimately back to the channel in 'beaver pond return flows'. Subsurface return flow seepages were universally suboxic, while ponds and surface return flows showed a range of oxygen concentration due to in-situ photosynthesis and atmospheric mixing. Median concentrations of reduced metals: manganese (Mn), iron (Fe), aluminum (Al), and arsenic (As) were substantially higher along beaver-induced flowpaths than in geologically controlled seepages and upstream main channel locations. The areal footprint of reduced return seepage flowpaths were imaged with surface electromagnetic methods, indicating extensive zones of high-conductivity shallow groundwater flowing back toward the main channels and emerging at relatively warm bank seepage zones observed with infrared. Multiple-depth redox dynamics within one focused seepage zone showed coupled variation over time, likely driven by observed changes in seepage rate that may be controlled by pond stage. High-resolution times series of dissolved Mn and Fe collected downstream of the beaver-impacted reaches demonstrated seasonal dynamics in mixed river metal concentrations. Al time series concentrations showed proportional change to Fe at the smaller stream location, indicating chemically reduced flowpaths were sourcing Al to the channel. Overall our results indicated beaver-induced floodplain exchanges create important, and perhaps dominant, transport pathways for floodplain metals by expanding chemically-reduced zones paired with strong advective exchange.