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.

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.

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.

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.

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.

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.

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.

Stream temperature data are useful for deciphering watershed processes important for aquatic ecosystems. Accurately extracting signal trends from stream temperature is essential for predicting responses of environmental and ecological indicators to change. Missing data periods are common for various reasons, and pose a challenge for scientists using temperature signal analysis to support stream research and ecological management objectives. However, the sensitivity of estimated temperature signal patterns to missing data has not been thoroughly evaluated, despite the potentially large impact on interpretation. In this study, we explored the effects of simulated missing daily data on the characterization of annual water temperature signals measured at headwater sites in the Pacific Northwest and Mid-Atlantic regions of the USA. For each site, we used linear regressions of sine-waves fitted to complete (365-d) and partial (7–357 consecutive missing data points) annual datasets of daily mean water temperature and computed three thermal parameters (mean, phase, and amplitude), which together can indicate thermally and ecologically influential watershed processes (e.g., depth and magnitude of groundwater discharge). Expected values (derived from complete datasets) ranged from 7.0 to 12.6 °C, 205 to 254 d, and 1.9 to 9.5 °C for annual mean, phase, and amplitude, respectively. While annual phase and amplitude could be accurately estimated (i.e., within 95–99% confidence intervals of expected values) with up to approximately two months of consecutively missing data, annual mean temperature required more complete datasets. We found that datasets with less than seven weeks of consecutively missing data enabled estimation of all annual signal parameters with reasonable accuracy (>75% probability of being within the 95–99% confidence intervals of expected values). Imputation of missing data expanded this range to approximately 20 weeks, with the greatest improvements in parameter estimation between 9 and 27 weeks of imputed missing data. However, caution should be exercised when applying this technique. For example, imputation improved the accuracy of parameter estimation for most sites, but accuracy decreased for some sites exhibiting strong groundwater influence. The timing of consecutive missing data points within a year had inconsistent effects on annual thermal parameter estimates among regions, years, and individual parameters. Utilizing sites with more than approximately seven consecutive weeks of missing data or 20 weeks of imputed data increases the probability of mischaracterization of annual stream thermal regimes. Understanding this limitation is vital for identifying the potential of streams to serve as climate refugia for ecological indicator species and effective future management of stream systems.

Freshwater lenses underlying small ocean islands exhibit spatial variability and temporal fluctuations in volume, influencing ecologic management. For example, The Palmyra Atoll National Wildlife Refuge harbors one of the few surviving native stands of Pisonia grandis in the central Pacific Ocean, yet these trees face pressure from groundwater salinization, with little basic groundwater data to guide decision making. Adding to natural complexity, the geology of Palmyra was heavily altered by dredge and fill activities. Our study based at this atoll combines geophysical and hydrological field measurements from 2008 to 2019 with groundwater modeling to study the drivers of observed freshwater lens dynamics. Electromagnetic induction (EMI) field data were collected on the main atoll islands over repeat transects in 2008 following 'strong' La Niña conditions (wet) and in 2016 during 'very strong' El Niño conditions (dry). Shallow monitoring wells were installed adjacent to the geophysical transects in 2013 and screened within the fresh/saline groundwater transition zone. Temporal EMI and monitoring well data showed a strong contraction of the freshwater lens in response to El Niño conditions, and indicated a thicker lens toward the ocean side, an opposite spatial pattern to that observed for many other Pacific islands. On an outer islet where a stand of mature Pisonia trees exist, EMI surveys revealed only a thin (3 m from land surface) layer of brackish groundwater during El Niño. Numerical groundwater simulations were performed for a range of permeability distributions and climate conditions at Palmyra. Results revealed that the observed atypical lens asymmetry is likely due to more efficient submarine groundwater discharge on the lagoon side as a result of lagoon dredging and filling with high-permeability material. Simulations also predict large decreases (40%) in freshwater lens volume during dry cycles and highlight threats to the Pisonia trees, yielding insight for atoll ecosystem management worldwide.

Groundwater discharge zones in streams are important habitats for aquatic organisms. The use of discharge zones for thermal refuge and spawning by fish and other biota renders them susceptible to potential focused discharge of groundwater contamination. Currently, there is a paucity of information about discharge zones as a potential exposure pathway of chemicals to stream ecosystems. Using thermal mapping technologies to locate groundwater discharges, shallow groundwater and surface water from three rivers in the Chesapeake Bay Watershed, USA were analyzed for phytoestrogens, pesticides and their degradates, steroid hormones, sterols and bisphenol A. A Bayesian censored regression model was used to compare groundwater and surface water chemical concentrations. The most frequently detected chemicals in both ground and surface water were the phytoestrogens genistein (79%) and formononetin (55%), the herbicides metolachlor (50%) and atrazine (74%), and the sterol cholesterol (88%). There was evidence suggesting groundwater discharge zones could be a unique exposure pathway of chemicals to surface water systems, in our case, metolachlor sulfonic acid (posterior mean concentration = 150 ng/L in groundwater and 4.6 ng/L in surface water). Our study also demonstrated heterogeneity of chemical concentration in groundwater discharge zones within a stream for the phytoestrogen formononetin, the herbicides metolachlor and atrazine, and cholesterol. Results support the hypothesis that discharge zones are an important source of exposure of phytoestrogens and herbicides to aquatic organisms. To manage critical resources within the Chesapeake Bay Watershed, more work is needed to characterize exposure in discharge zones more broadly across time and space.

Zones of strong groundwater upwelling to streams enhance thermal stability and moderate thermal extremes, which is particularly important to aquatic ecosystems in a warming climate. Passive thermal tracer methods used to quantify vertical upwelling rates rely on downward conduction of surface temperature signals. However, moderate to high groundwater flux rates (>−1.5 m d−1) restrict downward propagation of diurnal temperature signals, and therefore the applicability of several passive thermal methods. Active streambed heating from within high-resolution fiber-optic temperature sensors (A-HRTS) has the potential to define multidimensional fluid-flux patterns below the extinction depth of surface thermal signals, allowing better quantification and separation of local and regional groundwater discharge. To demonstrate this concept, nine A-HRTS were emplaced vertically into the streambed in a grid with ∼0.40 m lateral spacing at a stream with strong upward vertical flux in Mashpee, Massachusetts, USA. Long-term (8–9 h) heating events were performed to confirm the dominance of vertical flow to the 0.6 m depth, well below the extinction of ambient diurnal signals. To quantify vertical flux, short-term heating events (28 min) were performed at each A-HRTS, and heat-pulse decay over vertical profiles was numerically modeled in radial two dimension (2-D) using SUTRA. Modeled flux values are similar to those obtained with seepage meters, Darcy methods, and analytical modeling of shallow diurnal signals. We also observed repeatable differential heating patterns along the length of vertically oriented sensors that may indicate sediment layering and hyporheic exchange superimposed on regional groundwater discharge.

Analytical solutions that use diurnal temperature signals to estimate vertical fluxes between groundwater and surface water based on either amplitude ratios (Ar ) or phase shifts (Δϕ) produce results that rarely agree. Analytical solutions that simultaneously utilize Ar and Δϕ within a single solution have more recently been derived, decreasing uncertainty in flux estimates in some applications. Benefits of combined (Ar Δϕ) methods also include that thermal diffusivity and sensor spacing can be calculated. However, poor identification of either Ar or Δϕ from raw temperature signals can lead to erratic parameter estimates from Ar Δϕ methods. An add-on program for VFLUX 2 is presented to address this issue. Using thermal diffusivity selected from an Ar Δϕ method during a reliable time period, fluxes are recalculated using an Ar method. This approach maximizes the benefits of the Ar and Ar Δϕ methods. Additionally, sensor spacing calculations can be used to identify periods with unreliable flux estimates, or to assess streambed scour. Using synthetic and field examples, the use of these solutions in series was particularly useful for gaining conditions where fluxes exceeded 1 m/d.

Quantifying flow of groundwater through streambeds often is difficult due to the complexity of aquifer-scale heterogeneity combined with local-scale hyporheic exchange. We used fiber-optic distributed temperature sensing (FO-DTS), seepage meters, and vertical temperature profiling to locate, quantify, and monitor areas of focused groundwater discharge in a geomorphically simple sand-bed stream. This combined approach allowed us to rapidly focus efforts at locations where prodigious amounts of groundwater discharged to the Quashnet River on Cape Cod, Massachusetts, northeastern USA. FO-DTS detected numerous anomalously cold reaches one to several m long that persisted over two summers. Seepage meters positioned upstream, within, and downstream of 7 anomalously cold reaches indicated that rapid groundwater discharge occurred precisely where the bed was cold; median upward seepage was nearly 5 times faster than seepage measured in streambed areas not identified as cold. Vertical temperature profilers deployed next to 8 seepage meters provided diurnal-signal-based seepage estimates that compared remarkably well with seepage-meter values. Regression slope and R2 values both were near 1 for seepage ranging from 0.05 to 3.0 m d−1. Temperature-based seepage model accuracy was improved with thermal diffusivity determined locally from diurnal signals. Similar calculations provided values for streambed sediment scour and deposition at subdaily resolution. Seepage was strongly heterogeneous even along a sand-bed river that flows over a relatively uniform sand and fine-gravel aquifer. FO-DTS was an efficient method for detecting areas of rapid groundwater discharge, even in a strongly gaining river, that can then be quantified over time with inexpensive streambed thermal methods.

The distribution of shallow frozen ground is paramount to research in cold regions, and is subject to temporal and spatial changes influenced by climate, landscape disturbance and ecosystem succession. Remote sensing from airborne and satellite platforms is increasing our understanding of landscape-scale permafrost distribution, but typically lacks the resolution to characterise finer-scale processes and phenomena, which are better captured by integrated surface geophysical methods. Here, we demonstrate the use of electrical resistivity imaging (ERI), electromagnetic induction (EMI), ground penetrating radar (GPR) and infrared imaging over multiple summer field seasons around the highly dynamic Twelvemile Lake, Yukon Flats, central Alaska, USA. Twelvemile Lake has generally receded in the past 30 yr, allowing permafrost aggradation in the receded margins, resulting in a mosaic of transient frozen ground adjacent to thick, older permafrost outside the original lakebed. ERI and EMI best evaluated the thickness of shallow, thin permafrost aggradation, which was not clear from frost probing or GPR surveys. GPR most precisely estimated the depth of the active layer, which forward electrical resistivity modelling indicated to be a difficult target for electrical methods, but could be more tractable in time-lapse mode. Infrared imaging of freshly dug soil pit walls captured active-layer thermal gradients at unprecedented resolution, which may be useful in calibrating emerging numerical models. GPR and EMI were able to cover landscape scales (several kilometres) efficiently, and new analysis software showcased here yields calibrated EMI data that reveal the complicated distribution of shallow permafrost in a transitional landscape. Copyright © 2016 John Wiley & Sons, Ltd.

Despite decades of research into air and stream temperature dynamics, paired air-water annual temperature signals have been underutilized to characterize watershed processes. Annual stream temperature dynamics are useful in classifying fundamental thermal regimes and can enhance process-based interpretation of stream temperature controls, including deep and shallow groundwater discharge, when paired with air signals. In this study, we investigated multi-scale variability in annual paired air-water temperature patterns using sine-wave linear regressions of multi-year daily temperature data from streams of various sizes. A total of 311 sites from two spatially intensive regional datasets (Shenandoah National Park and Olympic Experimental State Forest) and a spatially extensive national dataset spanning the contiguous United States (U.S. Geological Survey gages) were evaluated. We calculated three annual air-water thermal metrics (mean ratio, phase lag, and amplitude ratio) to deduce the influence of groundwater and other watershed processes on stream thermal regimes at multiple spatial scales. Site-specific values of the three annual air-water thermal metrics ranged from 0.69 to 5.29 (mean ratio), −9 to 40 days (phase lag), and 0.29 to 1.12 (amplitude ratio). Regional patterns in the annual thermal metrics revealed persistent yet spatially variable influences of shallow groundwater discharge and high levels of thermal variability within watersheds, indicating the importance of local hydrogeological controls on stream temperature. Furthermore, annual thermal metric patterns from the regional datasets were generally concordant with the national dataset suggesting the utility of these annual thermal metrics for analysis at multiple scales. Analysis of the national dataset showed that previously defined thermal regimes based on water temperature alone could be further refined using air-water metrics and these metrics were related to physiographic watershed characteristics such as contributing area, elevation, and slope. This research demonstrates the importance of spatial scale and heterogeneity for inferring hydrological process in streams and provides guidance for the interpretation of annual air-water temperature metrics that can be efficiently applied to the growing database of multi-year temperature records. Results from this research can aid in the prediction of future thermal habitat suitability for coldwater-adapted species at ecologically and management-relevant spatial scales with readily available data.

The characterization of pore-space connectivity in porous media at the sediment/water interface is critical in understanding contaminant transport and reactive biogeochemical processes in zones of groundwater and surface-water exchange. Previous in situ studies of dual-domain (i.e., mobile/less-mobile porosity) systems have been limited to solute tracer injections at scales of meters to hundreds of meters and subsequent numerical model parameterization using fluid concentration histories. Pairing fine-scale (e.g., sub-meter) geoelectrical measurements with fluid tracer data over time alleviates dependence on flowpath-scale experiments, enabling spatially targeted characterization of shallow sediment/water interface media where biogeochemical reactivity is often high. The Dual-Domain Porosity Apparatus is a field-tested device capable of variable rate-controlled downward flow experiments. The Dual-Domain Porosity Apparatus facilitates inference of dual-domain parameters, i.e., mobile/less-mobile exchange rate coefficient and the ratio of less mobile to mobile porosity. The Dual-Domain Porosity Apparatus experimental procedure uses water electrical conductivity as a conservative tracer of differential loading and flushing of pore spaces within the region of measurement. Variable injection rates permit the direct quantification of the flow-dependence of dual-domain parameters, which has been theorized for decades but remains challenging to assess using existing experimental methodologies.

Brook trout (Salvelinus fontinalis) spawn in fall, and overwintering egg development can benefit from stable, relatively warm temperatures in groundwater seepage zones. However, eggs also are sensitive to dissolved oxygen concentration, which may be reduced in discharging groundwater. We investigated a 2-km reach of the coastal Quashnet River, Cape Cod, Massachusetts, USA, to relate preferred fish spawning habitat to geology, geomorphology, and groundwater discharge. Thermal reconnaissance methods were used to locate zones of rapid groundwater discharge, which were predominantly found along the center channel of a wider stream valley section. Pore-water chemistry and temporal vertical groundwater flux were measured at a subset of these zones during field campaigns over several seasons. Seepage zones in open valley sub-reaches generally showed suboxic conditions and higher dissolved solutes compared to the underlying glacial outwash aquifer. These discharge zones were cross-referenced with preferred brook trout redds, evaluated during 10 yr of observation, all of which were associated with discrete alcove features in steep cut banks where stream meander bends intersect the glacial valley walls. Seepage in these repeat spawning zones was generally stronger and more variable than open valley sites, with higher dissolved oxygen and reduced solute concentrations. The combined evidence indicates that regional groundwater discharge along the broader valley bottom is predominantly suboxic due to the influence of near-stream organic deposits; trout show no obvious preference for these zones when spawning. However, the meander bends that cut into sandy deposits near the valley walls generate strong, oxic seepage zones that are utilized routinely for redd construction and the overwintering of trout eggs. In similar coastal systems with extensive valley peat deposits, specific use of groundwater discharge points by brook trout may be limited to morphologies such as cut banks where groundwater flowpaths can short circuit buried organic material and remain oxygen rich.