Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Interfaces”

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 757 records · Page 42Linked to original sources
USGS2010

Assessment of basin-scale hydrologic impacts of CO2 sequestration, Illinois basin

Idealized, basin-scale sharp-interface models of CO2 injection were constructed for the Illinois basin. Porosity and permeability were decreased with depth within the Mount Simon Formation. Eau Claire confining unit porosity and permeability were kept fixed. We used 726 injection wells located near 42 power plants to deliver 80 million metric tons of CO2/year. After 100 years of continuous injection, deviatoric fluid pressures varied between 5.6 and 18 MPa across central and southern part of the Illinois basin. Maximum deviatoric pressure reached about 50% of lithostatic levels to the south. The pressure disturbance (>0.03 MPa) propagated 10-25 km away from the injection wells resulting in significant well-well pressure interference. These findings are consistent with single-phase analytical solutions of injection. The radial footprint of the CO2 plume at each well was only 0.5-2 km after 100 years of injection. Net lateral brine displacement was insignificant due to increasing radial distance from injection well and leakage across the Eau Claire confining unit. On geologic time scales CO2 would migrate northward at a rate of about 6 m/1000 years. Because of paleo-seismic events in this region (M5.5-M7.5), care should be taken to avoid high pore pressures in the southern Illinois basin. ?? 2010 Elsevier Ltd.

International Journal of Greenhouse Gas Control↗
USGS2007

Wave field features of shallow vertical discontinuity and their application in non-destructive detection

The geotechnical integrity of critical infrastructure can be seriously compromised by the presence of fractures or crevices. Non-destructive techniques to accurately detect fractures in critical infrastructure such as dams and highways could be of significant benefit to the geotechnical industry. This paper investigates the application of shallow seismic and georadar methods to the detection of a vertical discontinuity using numerical simulations. The objective is to address the kinematical analysis of a vertical discontinuity, determine the resulting wave field characteristics, and provide the basis for determining the existence of vertical discontinuities based on the recorded signals. Simulation results demonstrate that: (1) A reflection from a vertical discontinuity produces a hyperbolic feature on a seismic or georadar profile; (2) In order for a reflection from a vertical discontinuity to be produced, a reflecting horizon below the discontinuity must exist, the offset between source and receiver (x0) must be non-zero, on the same side of the vertical discontinuity; (3) The range of distances from the vertical discontinuity where a reflection event is observed is proportional to its length and to x0; (4) Should the vertical crevice (or fracture) pass through a reflecting horizon, dual hyperbolic features can be observed on the records, and this can be used as a determining factor that the vertical crevice passes through the interface; and (5) diffractions from the edges of the discontinuity can be recorded with relatively smaller amplitude than reflections and their ranges are not constrained by the length of discontinuity. If the length of discontinuity is short enough, diffractions are the dominant feature. Real-world examples show that the shallow seismic reflection method and the georadar method are capable of recording the hyperbolic feature, which can be interpreted as vertical discontinuity. Thus, these methods show some promise as effective non-destructive detection methods for locating vertical discontinuities (e.g., fractures or crevices) in infrastructure such as dams and highway pavement. ?? 2007 Elsevier B.V. All rights reserved.

Journal of Applied Geophysics↗
USGS2012

Geophysical investigations of geology and structure at the Martis Creek Dam, Truckee, California

A recent evaluation of Martis Creek Dam highlighted the potential for dam failure due to either seepage or an earthquake on nearby faults. In 1972, the U.S. Army Corps of Engineers constructed this earthen dam, located within the Truckee Basin to the north of Lake Tahoe, CA for water storage and flood control. Past attempts to raise the level of the Martis Creek Reservoir to its design level have been aborted due to seepage at locations downstream, along the west dam abutment, and at the base of the spillway. In response to these concerns, the U.S. Geological Survey has undertaken a comprehensive suite of geophysical investigations aimed at understanding the interplay between geologic structure, seepage patterns, and reservoir and groundwater levels. This paper concerns the geologic structure surrounding Martis Creek Dam and emphasizes the importance of a regional-scale understanding to the interpretation of engineering-scale geophysical data. Our studies reveal a thick package of sedimentary deposits interbedded with Plio-Pleistocene volcanic flows; both the deposits and the flows are covered by glacial outwash. Magnetic field data, seismic tomography models, and seismic reflections are used to determine the distribution and chronology of the volcanic flows. Previous estimates of depth to basement (or the thickness of the interbedded deposits) was 100 m. Magnetotelluric soundings suggest that electrically resistive bedrock may be up to 2500 m deep. Both the Polaris Fault, identified outside of the study area using airborne LiDAR, and the previously unnamed Martis Creek Fault, have been mapped through the dam area using ground and airborne geophysics. Finally, as determined by direct-current resistivity imaging, time-domain electromagnetic sounding, and seismic refraction, the paleotopography of the interface between the sedimentary deposits and the overlying glacial outwash plays a principal role both in controlling groundwater flow and in the distribution of the observed seepage.

California↗
USGS2013

Electrical signatures of ethanol-liquid mixtures: implications for monitoring biofuels migration in the subsurface

Ethanol (EtOH), an emerging contaminant with potential direct and indirect environmental effects, poses threats to water supplies when spilled in large volumes. A series of experiments was directed at understanding the electrical geophysical signatures arising from groundwater contamination by ethanol. Conductivity measurements were performed at the laboratory scale on EtOH–water mixtures (0 to 0.97 v/v EtOH) and EtOH–salt solution mixtures (0 to 0.99 v/v EtOH) with and without a sand matrix using a conductivity probe and a four-electrode electrical measurement over the low frequency range (1–1000 Hz). A Lichtenecker–Rother (L–R) type mixing model was used to simulate electrical conductivity as a function of EtOH concentration in the mixture. For all three experimental treatments increasing EtOH concentration resulted in a decrease in measured conductivity magnitude (|σ|). The applied L–R model fitted the experimental data at concentration ≤ 0.4 v/v EtOH, presumably due to predominant and symmetric intermolecular (EtOH–water) interaction in the mixture. The deviation of the experimental |σ| data from the model prediction at higher EtOH concentrations may be associated with hydrophobic effects of EtOH–EtOH interactions in the mixture. The |σ| data presumably reflected changes in relative strength of the three types of interactions (water–water, EtOH–water, and EtOH–EtOH) occurring simultaneously in EtOH–water mixtures as the ratio of EtOH to water changed. No evidence of measurable polarization effects at the EtOH–water and EtOH–water–mineral interfaces over the investigated frequency range was found. Our results indicate the potential for using electrical measurements to characterize and monitor EtOH spills in the subsurface.

Journal of Contaminant Hydrology↗
USGS2013

Complex resistivity signatures of ethanol in sand-clay mixtures

We performed complex resistivity (CR) measurements on laboratory columns to investigate changes in electrical properties as a result of varying ethanol (EtOH) concentration (0% to 30% v/v) in a sand–clay (bentonite) matrix. We applied Debye decomposition, a phenomenological model commonly used to fit CR data, to determine model parameters (time constant: τ, chargeability: m, and normalized chargeability: m n ). The CR data showed a significant (P ≤ 0.001) time-dependent variation in the clay driven polarization response (~ 12 mrad) for 0% EtOH concentration. This temporal variation probably results from the clay–water reaction kinetics trending towards equilibrium in the sand–clay–water system. The clay polarization is significantly suppressed (P ≤ 0.001) for both measured phase (ϕ) and imaginary conductivity (σ″) with increasing EtOH concentration. Normalized chargeability consistently decreases (by up to a factor of ~ 2) as EtOH concentration increases from 0% to 10% and 10 to 20%, respectively. We propose that such suppression effects are associated with alterations in the electrical double layer (EDL) at the clay–fluid interface due to (a) strong EtOH adsorption on clay, and (b) complex intermolecular EtOH–water interactions and subsequent changes in ionic mobility on the surface in the EDL. Changes in the CR data following a change of the saturating fluid from EtOH 20% to plain water indicate strong hysteresis effects in the electrical response, which we attribute to persistent EtOH adsorption on clay. Our results demonstrate high sensitivity of CR measurements to clay–EtOH interactions in porous media, indicating the potential application of this technique for characterization and monitoring of ethanol contamination in sediments containing clays.

Journal of Contaminant Hydrology↗
USGS2012

The impact of antecedent fire area on burned area in southern California coastal ecosystems

Frequent wildfire disasters in southern California highlight the need for risk reduction strategies for the region, of which fuel reduction via prescribed burning is one option. However, there is no consensus about the effectiveness of prescribed fire in reducing the area of wildfire. Here, we use 29 years of historical fire mapping to quantify the relationship between annual wildfire area and antecedent fire area in predominantly shrub and grassland fuels in seven southern California counties, controlling for annual variation in weather patterns. This method has been used elsewhere to measure leverage: the reduction in wildfire area resulting from one unit of prescribed fire treatment. We found little evidence for a leverage effect (leverage = zero). Specifically our results showed no evidence that wildfire area was negatively influenced by previous fires, and only weak relationships with weather variables rainfall and Santa Ana wind occurrences, which were variables included to control for inter-annual variation. We conclude that this is because only 2% of the vegetation burns each year and so wildfires rarely encounter burned patches and chaparral shrublands can carry a fire within 1 or 2 years after previous fire. Prescribed burning is unlikely to have much influence on fire regimes in this area, though targeted treatment at the urban interface may be effective at providing defensible space for protecting assets. These results fit an emerging global model of fire leverage which position California at the bottom end of a continuum, with tropical savannas at the top (leverage = 1: direct replacement of wildfire by prescribed fire) and Australian eucalypt forests in the middle (leverage ∼ 0.25).

California↗
USGS2018

Evaluating long-term patterns of decreasing groundwater discharge through a lake-bottom permeable reactive barrier

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.

Journal of Environmental Management↗
USGS2018

Regeneration of Salicaceae riparian forests in the Northern Hemisphere: A new framework and management tool

Human activities on floodplains have severely disrupted the regeneration of foundation riparian shrub and tree species of the Salicaceae family ( Populus and Salix spp.) throughout the Northern Hemisphere. Restoration ecologists initially tackled this problem from a terrestrial perspective that emphasized planting. More recently, floodplain restoration activities have embraced an aquatic perspective, inspired by the expanding practice of managing river flows to improve river health (environmental flows) . However, riparian Salicaceae species occupy floodplain and riparian areas, which lie at the interface of both terrestrial and aquatic ecosystems along watercourses . Thus, their regeneration depends on a complex interaction of hydrologic and geomorphic processes that have shaped key life-cycle requirements for seedling establishment . Ultimately, restoration needs to integrate these concepts to succeed. However, while regeneration of Salicaceae is now reasonably well-understood, the literature reporting restoration actions on Salicaceae regeneration is sparse, and a specific theoretical framework is still missing. Here, we have reviewed 105 peer-reviewed published experiences in restoration of Salicaceae forests, including 91 projects in 10 world regions, to construct a decision tree to inform restoration planning through explicit links between the well-studied biophysical requirements of Salicaceae regeneration and 17 specific restoration actions, the most popular being planting (in 55% of the projects), land contouring (30%), removal of competing vegetation (30%), site selection (26%), and irrigation (24%). We also identified research gaps related to Salicaceae forest restoration and discuss alternative, innovative and feasible approaches that incorporate the human component.

Journal of Environmental Management↗
USGS2011

Polymethylene-interrupted fatty acids: Biomarkers for native and exotic mussels in the Laurentian Great Lakes

Freshwater organisms synthesize a wide variety of fatty acids (FAs); however, the ability to synthesize and/or subsequently modify a particular FA is not universal, making it possible to use certain FAs as biomarkers. Herein we document the occurrence of unusual FAs (polymethylene-interrupted fatty acids; PMI-FAs) in select freshwater organisms in the Laurentian Great Lakes. We did not detect PMI-FAs in: (a) natural seston from Lake Erie and Hamilton Harbor (Lake Ontario), (b) various species of laboratory-cultured algae including a green alga ( Scenedesmus obliquus ), two cyanobacteria ( Aphanizomenon flos-aquae and Synechococystis sp.), two diatoms ( Asterionella formosa , Diatoma elongatum ) and a chrysophyte ( Dinobryon cylindricum ) or, (c) zooplankton ( Daphnia spp., calanoid or cyclopoid copepods) from Lake Ontario, suggesting that PMI-FAs are not substantively incorporated into consumers at the phytoplankton–zooplankton interface. However, these unusual FAs comprised 4-6% of total fatty acids (on a dry tissue weight basis) of native fat mucket ( Lampsilis siliquoidea ) and plain pocketbook ( L. cardium ) mussels and in invasive zebra ( Dreissena polymorpha ) and quagga ( D. bugensis ) mussels. We were able to clearly partition Great Lakes' mussels into three separate groups (zebra, quagga, and native mussels) based solely on their PMI-FA profiles. We also provide evidence for the trophic transfer of PMI-FAs from mussels to various fishes in Lakes Ontario and Michigan, further underlining the potential usefulness of PMI-FAs for tracking the dietary contribution of mollusks in food web and contaminant-fate studies.

Michigan, Ontario, Prince Edward Island, Wisconsin↗
USGS2013

Great Lakes rivermouth ecosystems: scientific synthesis and management implications

At the interface of the Great Lakes and their tributary rivers lies the rivermouths, a class of aquatic ecosystem where lake and lotic processes mix and distinct features emerge. Many rivermouths are the focal point of both human interaction with the Great Lakes and human impacts to the lakes; many cities, ports, and beaches are located in rivermouth ecosystems, and these human pressures often degrade key ecological functions that rivermouths provide. Despite their ecological uniqueness and apparent economic importance, there has been relatively little research on these ecosystems as a class relative to studies on upstream rivers or the open-lake waters. Here we present a synthesis of current knowledge about ecosystem structure and function in Great Lakes rivermouths based on studies in both Laurentian rivermouths, coastal wetlands, and marine estuarine systems. A conceptual model is presented that establishes a common semantic framework for discussing the characteristic spatial features of rivermouths. This model then is used to conceptually link ecosystem structure and function to ecological services provided by rivermouths. This synthesis helps identify the critical gaps in understanding rivermouth ecology. Specifically, additional information is needed on how rivermouths collectively influence the Great Lakes ecosystem, how human alterations influence rivermouth functions, and how ecosystem services provided by rivermouths can be managed to benefit the surrounding socioeconomic networks.

Great Lakes↗
USGS2015

Application of the FluEgg model to predict transport of Asian carp eggs in the Saint Joseph River (Great Lakes tributary)

The Fluvial Egg Drift Simulator (FluEgg) is a three-dimensional Lagrangian model that simulates the movement and development of Asian carp eggs until hatching based on the physical characteristics of the flow field and the physical and biological characteristics of the eggs. This tool provides information concerning egg development and spawning habitat suitability including: egg plume location, egg vertical and travel time distribution, and egg-hatching risk. A case study of the simulation of Asian carp eggs in the Lower Saint Joseph River, a tributary of Lake Michigan, is presented. The river hydrodynamic input for FluEgg was generated in two ways — using hydroacoustic data and using HEC-RAS model data. The HEC-RAS model hydrodynamic input data were used to simulate 52 scenarios covering a broad range of flows and water temperatures with the eggs at risk of hatching ranging from 0 to 93% depending on river conditions. FluEgg simulations depict the highest percentage of eggs at risk of hatching occurs at the lowest discharge and at peak water temperatures. Analysis of these scenarios illustrates how the interactive relation among river length, hydrodynamics, and water temperature influence egg transport and hatching risk. An improved version of FluEgg, which more realistically simulates dispersion and egg development, is presented. Also presented is a graphical user interface that facilitates the use of FluEgg and provides a set of post-processing analysis tools to support management decision-making regarding the prevention and control of Asian carp reproduction in rivers with or without Asian carp populations.

Michigan↗
USGS2026

The Great Lakes aquatic tissue analysis repository (GLATAR)

Empiricists and modellers use information on energy density, proximate composition, stable isotopes, fatty acids, thiamine, and bioaccumulative tracers (e.g., PCBs and mercury) to understand the state and inter-relationships of aquatic food webs. Data exist in many published and unpublished sources, but are not consolidated in an easily accessible database that would serve as a vital resource to i) provide basic estimates of these diet-derived measures of body composition, ii) understand sources of variation in the underlying data, iii) facilitate exploration and development of data proxies, and iv) assist in study design. We designed GLATAR (Great Lakes Aquatic Tissue Analysis Repository, glatar.org ) to address this need. GLATAR is an open-access, searchable database linked to a web-based toolbox to visualise and generate user-defined summaries on diet-derived ecological metrics, with a focus on taxa of importance to the Great Lakes. GLATAR currently contains over 50,000 records on energy density, chemical tracers, and proximate body composition from 67 species of fish and 72 invertebrate taxa. We hope this user-friendly interface will entice others to upload their published and unpublished data to the repository, enriching the breadth of data accessible to researchers. In this way, GLATAR will become a ‘living’ and interactive resource for empiricists and modellers working in freshwater ecosystems as they make critical decisions related to growth, production, and consumption across a diverse group of economically and ecologically important aquatic species.

Journal of Great Lakes Research↗
USGS2014

Alternative waste residue materials for passive in situ prevention of sulfide-mine tailings oxidation: A field evaluation

Novel solutions for sulfide-mine tailings remediation were evaluated in field-scale experiments on a former tailings repository in northern Sweden. Uncovered sulfide-tailings were compared to sewage-sludge biosolid amended tailings over 2 years. An application of a 0.2 m single-layer sewage-sludge amendment was unsuccessful at preventing oxygen ingress to underlying tailings. It merely slowed the sulfide-oxidation rate by 20%. In addition, sludge-derived metals (Cu, Ni, Fe, and Zn) migrated and precipitated at the tailings-to-sludge interface. By using an additional 0.6 m thick fly-ash sealing layer underlying the sewage sludge layer, a solution to mitigate oxygen transport to the underlying tailings and minimize sulfide-oxidation was found. The fly-ash acted as a hardened physical barrier that prevented oxygen diffusion and provided a trap for sludge-borne metals. Nevertheless, the biosolid application hampered the application, despite the advances in the effectiveness of the fly-ash layer, as sludge-borne nitrate leached through the cover system into the underlying tailings, oxidizing pyrite. This created a 0.3 m deep oxidized zone in 6-years. This study highlights that using sewage sludge in unconventional cover systems is not always a practical solution for the remediation of sulfide-bearing mine tailings to mitigate against sulfide weathering and acid rock drainage formation.

Skellefte�;Stockholm↗
USGS2008

A seepage meter designed for use in flowing water

Seepage meters provide one of the most direct means to measure exchange of water across the sediment-water interface, but they generally have been unsuitable for use in fluvial settings. Although the seepage bag can be placed inside a rigid container to minimize velocity head concerns, the seepage cylinder installed in the sediment bed projects into and disrupts the flow field, altering both the local-scale fluid exchange as well as measurement of that exchange. A low-profile seepage meter designed for use in moving water was tested in a seepage meter flux tank where both current velocity and seepage velocity could be controlled. The conical seepage cylinder protrudes only slightly above the sediment bed and is connected via tubing to a seepage bag or flowmeter positioned inside a rigid shelter that is located nearby where current velocity is much slower. Laboratory and field tests indicate that the net effect of the small protrusion of the seepage cylinder into the surface water flow field is inconsequentially small for surface water currents up to 65 cm s-1. Current velocity affects the variability of seepage measurements; seepage standard deviation increased from ???2 to ???6 cm d-1 as current velocity increased from 9 to 65 cm s-1. Substantial bias can result if the shelter is not placed to minimize hydraulic gradient between the bag and the seepage cylinder.

Journal of Hydrology↗
USGS2011

Subsurface transport of orthophosphate in five agricultural watersheds, USA

Concentrations of dissolved orthophosphate (ortho P) in the unsaturated zone, groundwater, tile drains, and groundwater/stream water interfaces were assessed in five agricultural watersheds to determine the potential for subsurface transport. Concentrations of iron oxides were measured in the aquifer material and adsorption of ortho P on oxide surfaces was assessed by geochemical modeling. Attenuation of ortho P in these aquifers was attributed primarily to sorption onto iron oxides, and in one location onto clay minerals. Only one location showed a clear indication of phosphorus transport to a stream from groundwater discharge, although groundwater did contribute to the stream load elsewhere. Subsurface ortho P movement at a site in California resulted in a plume down gradient from orchards, which was attenuated by a 200 m thick riparian zone with natural vegetation. Iron oxides had an effect on phosphorus movement and concentrations at all locations, and groundwater chemistry, especially pH, exerted a major control on the amount of phosphorus adsorbed. Groundwater pH at a site in Maryland was below 5 and that resulted in complete sequestration of phosphorus and no movement toward the stream. Geochemical modeling indicated that as the surfaces approached saturation, groundwater concentrations of ortho P rise rapidly.

Journal of Hydrology↗
USGS2014

Quantity and quality of groundwater discharge in a hypersaline lake environment

Geophysical and geochemical surveys were conducted to understand groundwater discharge to Great Salt Lake (GSL) and assess the potential significance of groundwater discharge as a source of selenium (Se). Continuous resistivity profiling (CRP) focusing below the sediment/water interface and fiber-optic distributed temperature sensing (FO-DTS) surveys were conducted along the south shore of GSL. FO-DTS surveys identified persistent cold-water temperature anomalies at 10 separate locations. Seepage measurements were conducted at 17 sites (mean seepage rate = 0.8 cm/day). High resistivity anomalies identified by the CRP survey were likely a mirabilite (Na 2 SO 4 ·10H 2 O) salt layer acting as a semi-confining layer for the shallow groundwater below the south shore of the lake. Positive seepage rates measured along the near-shore areas of GSL indicate that a ∼1-m thick oolitic sand overlying the mirabilite layer is likely acting as a shallow, unconfined aquifer. Using the average seepage rate of 0.8 cm/day over an area of 1.6 km 2 , an annual Se mass loading to GSL of 23.5 kg was estimated. Determination of R/Ra values (calculated 3 He/ 4 He ratio over the present-day atmospheric 3 He/ 4 He ratio) <1 and tritium activities of 1.2–2.0 tritium units in groundwater within and below the mirabilite layer indicates a convergence of regional and local groundwater flow paths discharging into GSL. Groundwater within and below the mirabilite layer obtains its high sulfate salinity from the dissolution of mirabilite. The δ 34 S and δ 18 O isotopic values in samples of dissolved sulfate from the shallow groundwater below the mirabilite are almost identical to the isotopic signature of the mirabilite core material. The saturation index calculated for groundwater samples using PHREEQC indicates the water is at equilibrium with mirabilite. Water samples collected from GSL immediately off shore contained Se concentrations that were 3–4 times higher than other sampling sites >25 km offshore from the study site and may be originating from less saline groundwater seeps mixing with the more saline water from GSL. Additional evidence for mixing with near shore seeps is found in the δD and δ 18 O isotopic values and Br:Cl ratios. Geochemical modeling for a water sample collected in the vicinity of the study area indicates that under chemically reducing conditions, arsenic- (As) bearing minerals could dissolve while Se-bearing minerals will likely precipitate out of solution, possibly explaining why the shallow groundwater below and within the mirabilite salt layer contains low concentrations of Se (0.9–2.3 μg/L).

Utah↗
USGS2014

Practical limitations on the use of diurnal temperature signals to quantify groundwater upwelling

Groundwater upwelling to streams creates unique habitat by influencing stream water quality and temperature; upwelling zones also serve as vectors for contamination when groundwater is degraded. Temperature time series data acquired along vertical profiles in the streambed have been applied to simple analytical models to determine rates of vertical fluid flux. These models are based on the downward propagation characteristics (amplitude attenuation and phase-lag) of the surface diurnal signal. Despite the popularity of these models, there are few published characterizations of moderate-to-strong upwelling. We attribute this limitation to the thermodynamics of upwelling, under which the downward conductive signal transport from the streambed interface occurs opposite the upward advective fluid flux. Governing equations describing the advection&ndash;diffusion of heat within the streambed predict that under upwelling conditions, signal amplitude attenuation will increase, but, counterintuitively, phase-lag will decrease. Therefore the extinction (measurable) depth of the diurnal signal is very shallow, but phase lag is also short, yielding low signal to noise ratio and poor model sensitivity. Conversely, amplitude attenuation over similar sensor spacing is strong, yielding greater potential model sensitivity. Here we present streambed thermal time series over a range of moderate to strong upwelling sites in the Quashnet River, Cape Cod, Massachusetts. The predicted inverse relationship between phase-lag and rate of upwelling was observed in the field data over a range of conditions, but the observed phase-lags were consistently shorter than predicted. Analytical solutions for fluid flux based on signal amplitude attenuation return results consistent with numerical models and physical seepage meters, but the phase-lag analytical model results are generally unreasonable. Through numerical modeling we explore reasons why phase-lag may have been over-predicted by the analytical models, and develop guiding relations of diurnal temperature signal extinction depth based on stream diurnal signal amplitude, upwelling magnitude, and streambed thermal properties that will be useful in designing future experiments.

Journal of Hydrology↗
USGS2019

Multi-scale preferential flow processes in an urban streambed under variable hydraulic conditions

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. <3 h) of the hyporheic zone with stream water, followed by longer-term (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.

Journal of Hydrology↗