SEARCH · Search USGS

Results for “Solutions Journal”

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

SourceResearch connectionFrom yearThrough year

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

At least 577 records · Page 32Linked to original sources

USGS2022

Microbial endophytes and compost improve plant growth in two contrasting types of hard rock mining waste

The re-vegetation of mining wastes with native plants is a comparatively low-cost solution for mine reclamation. However, re-vegetation fails when extreme pH values, low organic matter, or high concentrations of phytotoxic elements inhibit plant establishment and growth. Our aim was to determine whether the combined addition of municipal waste compost and diazotrophic endophytes (i.e., microorganisms that fix atmospheric N 2 and live within plants) could improve plant growth, organic matter accumulation, and phytostabilization of trace element contaminants in two types of hard rock mine waste. We grew a widespread native perennial grass, Bouteloua curtipendula , for one month in alkaline waste rock (porphyry copper mine) and tailings (Ag-Pb-Au mine, amended with dolomite) sourced from southeastern Arizona, United States. B. curtipendula tolerated elevated concentrations of multiple phytotoxic trace elements in the tailings (Mn, Pb, Zn), stabilizing them in roots without foliar translocation. Adding compost and endophyte seed coats improved plant growth, microbial biomass, and organic matter accumulation despite stark differences in the geochemical and physical characteristics of the mining wastes. The widespread grass B. curtipendula is a potential candidate for re-vegetating mine wastes when seeded with soil additives to increase pH and with microbial and organic amendments to increase plant growth.

International Journal of Phytoremediation↗

USGS2007

Transient storage assessments of dye-tracer injections in rivers of the Willamette Basin, Oregon

Rhodamine WT dye-tracer injections in rivers of the Willamette Basin yield concentration-time curves with characteristically long recession times suggestive of active transient storage processes. The scale of drainage areas contributing to the stream reaches studied in the Willamette Basin ranges from 10 to 12,000 km2. A transient storage assessment of the tracer studies has been completed using the U.S. Geological Survey's One-dimensional Transport with Inflow and Storage (OTIS) model, which incorporates storage exchange and decay functions along with the traditional dispersion and advection transport equation. The analysis estimates solute transport of the dye. It identifies first-order decay coefficients to be on the order of 10-5/sec for the nonconservative Rhodamine WT. On an individual subreach basis, the first-order decay is slower (typically by an order of magnitude) than the transient storage process, indicating that nonconservative tracers may be used to evaluate transient storage in rivers. In the transient storage analysis, a dimensionless parameter (As/A) expresses the spatial extent of storage zone area relative to stream cross section. In certain reaches of Willamette Basin pool-and-riffle, gravel-bed rivers, this parameter was as large as 0.5. A measure of the storage exchange flux was calculated for each stream subreach in the simulation analysis. This storage exchange is shown subjectively to be higher at higher stream discharges. Hyporheic linkage between streams and subsurface flows is the probable physical mechanism contributing to a significant part of this inferred active transient storage. Hyporheic linkages are further suggested by detailed measurements of river discharge with an Acoustic Doppler Current Profiler system delineating zones in two large rivers where water alternately enters and leaves the surface channels through graveland-cobble riverbeds. Measurements show patterns of hyporheic exchange that are highly variable in time and space.

Oregon↗

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↗

USGS2017

Benchmarking computational fluid dynamics models of lava flow simulation for hazard assessment, forecasting, and risk management

Numerical simulations of lava flow emplacement are valuable for assessing lava flow hazards, forecasting active flows, designing flow mitigation measures, interpreting past eruptions, and understanding the controls on lava flow behavior. Existing lava flow models vary in simplifying assumptions, physics, dimensionality, and the degree to which they have been validated against analytical solutions, experiments, and natural observations. In order to assess existing models and guide the development of new codes, we conduct a benchmarking study of computational fluid dynamics (CFD) models for lava flow emplacement, including VolcFlow, OpenFOAM, FLOW-3D, COMSOL, and MOLASSES. We model viscous, cooling, and solidifying flows over horizontal planes, sloping surfaces, and into topographic obstacles. We compare model results to physical observations made during well-controlled analogue and molten basalt experiments, and to analytical theory when available. Overall, the models accurately simulate viscous flow with some variability in flow thickness where flows intersect obstacles. OpenFOAM, COMSOL, and FLOW-3D can each reproduce experimental measurements of cooling viscous flows, and OpenFOAM and FLOW-3D simulations with temperature-dependent rheology match results from molten basalt experiments. We assess the goodness-of-fit of the simulation results and the computational cost. Our results guide the selection of numerical simulation codes for different applications, including inferring emplacement conditions of past lava flows, modeling the temporal evolution of ongoing flows during eruption, and probabilistic assessment of lava flow hazard prior to eruption. Finally, we outline potential experiments and desired key observational data from future flows that would extend existing benchmarking data sets.

Journal of Applied Volcanology↗

USGS1991

An analytical model for in situ extraction of organic vapors

This paper introduces a simple convective-flow model that can be used as a screening tool and for conducting sensitivity analyses for in situ vapor extraction of organic compounds from porous media. An assumption basic to this model was that the total mass of volatile organic chemicals (VOC) exists in three forms: as vapors, in the soil solution, and adsorbed to soil particles. The equilibrium partitioning between the vapor-liquid phase was described by Henry's law constants (K(H)) and between the liquid-soil phase by soil adsorption constants (K(d)) derived from soil organic carbon-water partition coefficients (K(oc)). The model was used to assess the extractability of 36 VOCs from a hypothetical site. Most of the VOCs appeared to be removable from soil by this technology, although modeling results suggested that rates for the alcohols and ketones may be very slow. In general, rates for weakly adsorbed compounds (K(oc) < 100 mL/g) were significantly higher when K(H) was greater than 10-4 atm??m3??mol-1. When K(oc) was greater than about 100 mL/g, the rates of extraction were sensitive to the amount of organic carbon present in the soil. The air permeability of the soil material (k) was a critical factor. In situ extraction needs careful evaluation when k is less than 10 millidarcies to determine its applicability. An increase in the vacuum applied to an extraction well accelerated removal rates but the diameter of the well had little effect. However, an increase in the length of the well screen open to the contaminated zone significantly affected removal rates, especially in low-permeability materials.This paper introduces a simple convective-flow model that can be used as a screening tool and for conducting sensitivity analyses for in situ vapor extraction of organic compounds from porous media. An assumption basic to this model was that the total mass of volatile organic chemicals (VOC) exists in three forms: as vapors, in the soil solution, and adsorbed to soil particles. The equilibrium partitioning between the vapor-liquid phase was described by Henry's law constants (KH) and between the liquid-soil phase by soil adsorption constants (Kd) derived from soil organic carbon-water partition coefficients (Koc). The model was used to assess the extractability of 36 VOCs from a hypothetical site. Most of the VOCs appeared to be removable from soil by this technology, although modeling results suggested that rates for the alcohols and ketones may be very slow. In general, rates for weakly adsorbed compounds (Koc < 100 mL/g) were significantly higher when KH was greater than 10-4atm-m3-mol-1. When Koc was greater than about 100 mL/g, the rates of extraction were sensitive to the amount of organic carbon present in the soil. The air permeability of the soil material (k) was a critical factor. In situ extraction needs careful evaluation when k is less than 10 millidarcies to determine its applicability. An increase in the vacuum applied to an extraction well accelerated removal rates but the diameter of the well had little effect. However, an increase in the length of the well screen open to the contaminated zone significantly affected removal rates, especially in low-permeability materials.

Journal of Hazardous Materials↗

USGS2026

Computational electromagnetic geophysics for groundwater system studies: A review on established practices and recent advances

Identifying effective solutions for locating groundwater resources and ensuring the quality of drinking water is increasingly urgent, given the challenges posed by climate change and population growth. This review investigates electromagnetic geophysical imaging techniques, in both time- and frequency-domain, that can provide valuable insights for groundwater assessment. We explore computational electromagnetic methods used to evaluate electromagnetic data and several recent hydrogeophysical case studies. As open-source frameworks for modeling electromagnetic geophysical problems become available, a broader range of researchers can interpret their data with computationally advanced software. We provide an overview of documented open-source codes for evaluating electromagnetic data and analyze various hydrological targets in relation to their electromagnetic surveying technique and the computational method applied. Furthermore, we evaluate the potential of advanced computational techniques, including three-dimensional modeling, non-deterministic inversion and machine learning, to couple geophysical with numerical groundwater modeling and apply it in groundwater system studies. Despite obstacles such as complexity and resource demands, our findings indicate that the quantification and integration of predictive uncertainties from both electromagnetic and hydrological data and simulations would significantly improve the reliability of hydrogeophysical models. This can lead to a deeper understanding of groundwater systems and improved management practices.

Journal of Hydrology↗

USGS2025

Rupture continuity through intermittent pauses in Cascadia slow slip events

Cascadia slow slip events (SSEs) are often envisioned as smooth, continuous ruptures, progressively activating tremor asperities as they propagate. Macroscopically, geodetic inversions and spatiotemporal maps of tremor epicenters show steady, uniform migration. In detail tremor is more chaotic and discontinuous. Larger long-term SSEs observed in daily geodetic solutions are inferred to exhibit intermittent pauses that reflect temporary re-locking of the fault, but this temporal resolution limits tests for similar re-locking on shorter timescales. We use temporal measurements of the areal growth and radiated energy of tremor clusters to investigate SSE intermittence. We find that ruptures mirror tremor pauses. Areal growth rate, however, does not reset, and removing the pauses results in smoother and more similar growth measurements among all SSEs. The rupture similarity occurs regardless of size or location and hints at an underlying uniformity and lack of predeterminism in eventual SSE size. Epicentral uncertainty precludes quantifying early rupture stages, but for larger events areal growth follows a power-law and slows with increasing size. Temporal correlations in tremor energy with inferred SSE propagation velocities and tremor rates suggest its use as a proxy for slip velocity. We find that tremor energy is tidally modulated at daily and sub-daily frequencies, and this modulation is continuous through pauses, suggesting a memory of slip state is sustained through them. We argue these pauses reflect unsteady propagation of the slip front, marked by rapid re- and un-locking, and excluding them removes rupture complexity to reveal a diffusive-like slip process and underlying universality in growth.

California, Oregon, Washington↗

USGS1996

Uncertainty of groundwater vulnerability assessments for agricultural regions in Hawaii: Review

There are important challenges associated with assessing potential groundwater vulnerability hazards that may result from regional scale applications of agrochemicals. The increasing availability of Geographic Information System (GIS) software to those involved in assisting with landuse decisions has resulted in the widespread production of multicolored risk management maps for many environmentally sensitive issues. Soil-based GIS's have recently been coupled to various solute-leaching models to make near-surface groundwater vulnerability assessments for guidance in pesticide regulation in several states. In general, these assessments rest on soil, climatic, and chemical data that are extremely sparse and contain considerable uncertainty. It is also important to acknowledge the uncertainty associated with the transport/fate processes that are not accounted for by the modeling approach used to make the assessment. In this paper, we review the results from a series of papers that have focused on characterization of uncertainty in pesticide mobility estimates, using the attenuation and retardation indices (AF and RF), for the Pearl Harbor Basin on the Hawaiian island of Oahu. Relative to data error uncertainties, we discuss the impacts of: (i) soil, climatic, and chemical data base uncertainties, (ii) reductions in data base uncertainties, (iii) extrapolation of soil data base information based on soil taxonomy and soil survey, and (iv) importing information from outside the region of interest. Relative to model error uncertainties, we compare pesticide leaching estimates from the simple AF and RF mobility indices with simulations from the EPA's Pesticide Root Zone Model (PRZM) and field observations. Finally, we outline a Regional Integrated Risk Assessment approach for characterizing regional scale groundwater vulnerability for near-surface nonpoint sources.

Hawaii↗

USGS2013

Seismic evidence for a slab tear at the Puerto Rico Trench

The fore-arc region of the northeast Caribbean plate north of Puerto Rico and the Virgin Islands has been the site of numerous seismic swarms since at least 1976. A 6 month deployment of five ocean bottom seismographs recorded two such tightly clustered swarms, along with additional events. Joint analyses of the ocean bottom seismographs and land-based seismic data reveal that the swarms are located at depths of 50–150 km. Focal mechanism solutions, found by jointly fitting P wave first-motion polarities and S/P amplitude ratios, indicate that the broadly distributed events outside the swarm generally have strike- and dip-slip mechanisms at depths of 50–100 km, while events at depths of 100–150 km have oblique mechanisms. A stress inversion reveals two distinct stress regimes: The slab segment east of 65°W longitude is dominated by trench-normal tensile stresses at shallower depths (50–100 km) and by trench-parallel tensile stresses at deeper depths (100–150 km), whereas the slab segment west of 65°W longitude has tensile stresses that are consistently trench normal throughout the depth range at which events were observed (50–100 km). The simple stress pattern in the western segment implies relatively straightforward subduction of an unimpeded slab, while the stress pattern observed in the eastern segment, shallow trench-normal tension and deeper trench-normal compression, is consistent with flexure of the slab due to rollback. These results support the hypothesis that the subducting North American plate is tearing at or near these swarms. The 35 year record of seismic swarms at this location and the recent increase in seismicity suggest that the tear is still propagating.

Puerto Rico↗

USGS1997

Continuous monitoring of surface deformation at Long Valley Caldera, California, with GPS

Continuous Global Positioning System (GPS) measurements at Long Valley Caldera, an active volcanic region in east central California, have been made on the south side of the resurgent dome since early 1993. A site on the north side of the dome was added in late 1994. Special adaptations for autonomous operation in remote regions and enhanced vertical precision were made. The data record ongoing volcanic deformation consistent with uplift and expansion of the surface above a shallow magma chamber. Measurement precisions (1 standard error) for “absolute” position coordinates, i.e., relative to a global reference frame, are 3–4 mm (north), 5–6 mm (east), and 10–12 mm (vertical) using 24 hour solutions. Corresponding velocity uncertainties for a 12 month period are about 2 mm/yr in the horizontal components and 3–4 mm/yr in the vertical component. High precision can also be achieved for relative position coordinates on short (less than 10 km) baselines using broadcast ephemerides and observing times as short as 3 hours, even when data are processed rapidly on site. Comparison of baseline length changes across the resurgent dome between the two GPS sites and corresponding two-color electronic distance measurements indicates similar extension rates within error (∼2 mm/yr) once we account for a random walk noise component in both systems that may reflect spurious monument motion. Both data sets suggest a pause in deformation for a 3.5 month period in mid-1995, when the extension rate across the dome decreased essentially to zero. Three dimensional positioning data from the two GPS stations suggest a depth (5.8±1.6 km) and location (west side of the resurgent dome) of a major inflation center, in agreement with other geodetic techniques, near the top of a magma chamber inferred from seismic data. GPS systems similar to those installed at Long Valley can provide a practical method for near real-time monitoring and hazard assessment on many active volcanoes.

Journal of Geophysical Research B: Solid Earth↗

USGS2003

Sequestration of phosphorus by acid mine drainage floc

Solubilization and transport of phosphorus (P) to the water environment is a critical environmental issue. Flocs resulting from neutralizing acid mine drainage (AMD) were tested as a possible low-cost amendment to reduce the loss of soluble P from agricultural fields and animal wastewater. Flocs were prepared by neutralizing natural and synthetic solutions of AMD with limestone, lime, ammonium hydroxide, and sodium hydroxide. Phosphorus sequestration was tested in three distinct environments: water, soil, and manure storage basins. In water, flocs prepared from AMD adsorbed 10 to 20 g P kg −1 dry floc in equilibrium with 1 mg L −1 soluble P. Similar results were observed for both Fe-based and Al-based synthetic flocs. A local soil sample adsorbed about 0.1 g P kg −1 , about two orders of magnitude less. The AMD-derived flocs were mixed with a high-P soil at 5 to 80 g floc kg −1 soil, followed by water and acid (Mehlich-1) extractions. All flocs performed similarly. About 70% of the water-extractable P was sequestered by the floc when applied at a rate of 20 g floc kg −1 soil, whereas plant-available P only decreased by about 30%. Under anaerobic conditions simulating manure storage basins, all AMD flocs reduced soluble P by greater than 95% at a rate of 0.2 g floc g −1 rainbow trout ( Oncorhynchus mykiss ) manure. These findings indicate that AMD flocs could be an effective agent for preventing soluble P losses from soil and manure to the water environment, while at the same time decreasing the costs associated with AMD treatment.

Journal of Environmental Quality↗

USGS2024

Evaluation of coal mine drainage and associated precipitates for radium and rare earth element concentrations

Coal mine drainage (CMD) and associated metal-rich precipitates have recently been proposed as unconventional sources of rare earth elements (REEs). However, the potential occurrence of radium (Ra), a known carcinogen, with the REE-bearing phases has not been investigated. We hypothesized that Ra may occur in solids that are precipitated from CMD as a “radiobarite” solid solution ((Ba,Sr,Ra)SO 4 ) and/or adsorbed with hydrous metal oxides. REEs have been documented to sorb or co-precipitate with iron (Fe), manganese (Mn), and aluminum (Al) oxyhydroxide in CMD solids. Likewise, Ra has been documented to sorb to hydrous Fe and Mn oxides especially where sulfate (SO 4 ) and/or barium (Ba) concentrations are insufficient to precipitate radiobarite. Thus, we conducted the first-ever survey of Ra concentrations in corresponding CMD water and solid samples in the United States. Samples were analyzed from 4 untreated and 9 treated CMD sites in both the bituminous and anthracite coal regions of Pennsylvania across a range of pH and SO 4 concentrations. The dissolved Ra in CMD was relatively low (<0.5 Bq/L), consistent with radiobarite solubility; however, CMD solids were largely composed of amorphous Fe, Al, and Mn oxyhydroxide and silicate minerals. Ra was associated with Mn-enriched CMD solids, upwards of 875 Bq/kg. Total REE + yttrium (Y) content in the CMD solids was enriched upwards of 3600 mg/kg and was significantly correlated with Al content. These preliminary results suggest that REE extraction may target Al-rich solids to avoid Ra in Mn-rich solids.

Pennsylvania↗

USGS1996

Exchangeable hydrogen explains the pH of spodosol Oa horizons

The chemistry of extremely acid Oa horizons does not conform to traditional pH, Al, and base saturation relationships. Results from two separate studies of northeastern U.S. forested soils were used to investigate relationships between pH in water or dilute salt solutions and other soil characteristics. In Oa horizons with pH below 4, soil pH in dilute CaCl 2 solution was correlated with exchangeable H + measured either by titration ( r = −0.88, P = 0.0001, n = 142) or by electrode ( r = −0.89, P = 0.0001, n = 45). Exchangeable H + expressed as a percentage of the cation-exchange capacity (CEC) was linear with pH and showed similar slopes for data from both studies. For all samples, pHw = 4.21 − 1.80 × H + /CEC ( R 2 = 0.69, n = 194). The reciprocal of the H + /CEC ratio is base saturation with Al added to the bases. Because of the low pH, exchangeable Al does not appear to behave as an acid. Exchangeable H + remains an operationally defined quantity because of the difficulty in separating exchange and hydrolysis reactions. In a variety of neutral-salt extractants, concentrations of H + were correlated with 0.1 M BaCl 2 -exchangeable H + ( r > 0.91, P = 0.0001, n = 26) regardless of the strength of the extract. Nine successive extractions with 0.33 m M CaCl 2 removed more H + than was removed by single batch extractions with either 1 M KCl or 0.1 M BaCl 2 (average H + of 70, 43, and 49 mmol kg −1 , respectively for 26 samples). The data showed little difference in the chemical behavior of Oa horizons from a variety of geographical sites and vegetation types.

Maine, New Hampshire, New York, Vermont↗

USGS2014

Real-time inversions for finite fault slip models and rupture geometry based on high-rate GPS data

We present an inversion strategy capable of using real-time high-rate GPS data to simultaneously solve for a distributed slip model and fault geometry in real time as a rupture unfolds. We employ Bayesian inference to find the optimal fault geometry and the distribution of possible slip models for that geometry using a simple analytical solution. By adopting an analytical Bayesian approach, we can solve this complex inversion problem (including calculating the uncertainties on our results) in real time. Furthermore, since the joint inversion for distributed slip and fault geometry can be computed in real time, the time required to obtain a source model of the earthquake does not depend on the computational cost. Instead, the time required is controlled by the duration of the rupture and the time required for information to propagate from the source to the receivers. We apply our modeling approach, called Bayesian Evidence-based Fault Orientation and Real-time Earthquake Slip, to the 2011 Tohoku-oki earthquake, 2003 Tokachi-oki earthquake, and a simulated Hayward fault earthquake. In all three cases, the inversion recovers the magnitude, spatial distribution of slip, and fault geometry in real time. Since our inversion relies on static offsets estimated from real-time high-rate GPS data, we also present performance tests of various approaches to estimating quasi-static offsets in real time. We find that the raw high-rate time series are the best data to use for determining the moment magnitude of the event, but slightly smoothing the raw time series helps stabilize the inversion for fault geometry.

Journal of Geophysical Research↗

USGS2012

Effects of sulfate ligand on uranyl carbonato surface species on ferrihydrite surfaces

Understanding uranium (U) sorption processes in permeable reactive barriers (PRB) are critical in modeling reactive transport for evaluating PRB performance at the Fry Canyon demonstration site in Utah, USA. To gain insight into the U sequestration mechanism in the amorphous ferric oxyhydroxide (AFO)-coated gravel PRB, U(VI) sorption processes on ferrihydrite surfaces were studied in 0.01 M Na 2 SO 4 solutions to simulate the major chemical composition of U-contaminatedgroundwater (i.e., [SO 4 2- ]~13 mM L -1 ) at the site. Uranyl sorption was greater at pH 7.5 than that at pH 4 in both air- and 2% pCO 2 -equilibrated systems. While there were negligible effects of sulfate ligands on the pH-dependent U(VI) sorption (<24 h) in both systems, X-ray absorption spectroscopy (XAS) analysis showed sulfate ligand associated U(VI) surface species at the ferrihydrite–water interface. In air-equilibrated systems, binary and mono-sulfate U(VI) ternary surface species co-existed at pH 5.43. At pH 6.55–7.83, a mixture of mono-sulfate and bis-carbonato U(VI) ternary surface species became more important. At 2% pCO 2 , there was no contribution of sulfate ligands on the U(VI) ternary surface species. Instead, a mixture of bis-carbonato inner-sphere (38%) and tris-carbonato outer-sphere U(VI) ternary surface species (62%) was found at pH 7.62. The study suggests that the competitive ligand (bicarbonate and sulfate) coordination on U(VI) surface species might be important in evaluating the U solid-state speciation in the AFO PRB at the study site where pCO 2 fluctuates between 1 and 2 pCO 2 %.

Journal of Colloid and Interface Science↗

USGS2021

Event scale relationships of DOC and TDN fluxes in throughfall and stemflow diverge from stream exports in a forested catchment

Aquatic fluxes of carbon and nutrients link terrestrial and aquatic ecosystems. Within forests, storm events drive both the delivery of carbon and nitrogen to the forest floor and the export of these solutes from the land via streams. To increase understanding of the relationships between hydrologic event character and the relative fluxes of carbon and nitrogen in throughfall, stemflow and streams, we measured dissolved organic carbon (DOC) and total dissolved nitrogen (TDN) concentrations in each flow path for 23 events in a forested watershed in Vermont, USA. DOC and TDN concentrations increased with streamflow, indicating their export was limited by water transport of catchment stores. DOC and TDN concentrations in throughfall and stemflow decreased exponentially with increasing precipitation, suggesting that precipitation removed a portion of available sources from tree surfaces during the events. DOC and TDN fluxes were estimated for 76 events across a 2-year period. For most events, throughfall and stemflow fluxes greatly exceeded stream fluxes, but the imbalance narrowed for larger storms (>30 mm). The largest 10 stream events exported 40% of all stream event DOC whereas those same 10 events contributed 14% of all throughfall export. Approximately 2–5 times more DOC and TDN was exported from trees during rain events than left the catchment via streams annually. The diverging influence of event size on tree versus stream fluxes has important implications for forested ecosystems as hydrological events increase in intensity and frequency due to climate change.

Vermont↗

USGS1999

Estimates of bottom roughness length and bottom shear stress in South San Francisco Bay, California

A field investigation of the hydrodynamics and the resuspension and transport of particulate matter in a bottom boundary layer was carried out in South San Francisco Bay (South Bay), California, during March-April 1995. Using broadband acoustic Doppler current profilers, detailed measurements of turbulent mean velocity distribution within 1.5 m above bed have been obtained. A global method of data analysis was used for estimating bottom roughness length z o and bottom shear stress (or friction velocities u * ). Field data have been examined by dividing the time series of velocity profiles into 24-hour periods and independently analyzing the velocity profile time series by flooding and ebbing periods. The global method of solution gives consistent properties of bottom roughness length z o and bottom shear stress values (or friction velocities u * ) in South Bay. Estimated mean values of z o and u * for flooding and ebbing cycles are different. The differences in mean z o and u * are shown to be caused by tidal current flood-ebb inequality, rather than the flooding or ebbing of tidal currents. The bed shear stress correlates well with a reference velocity; the slope of the correlation defines a drag coefficient. Forty-three days of field data in South Bay show two regimes of z o (and drag coefficient) as a function of a reference velocity. When the mean velocity is >25–30 cm s −1 , the ln z o (and thus the drag coefficient) is inversely proportional to the reference velocity. The cause for the reduction of roughness length is hypothesized as sediment erosion due to intensifying tidal currents thereby reducing bed roughness. When the mean velocity is <25–30 cm s −1 , the correlation between z o and the reference velocity is less clear. A plausible explanation of scattered values of z o under this condition may be sediment deposition. Measured sediment data were inadequate to support this hypothesis, but the proposed hypothesis warrants further field investigation.

Journal of Geophysical Research C: Oceans↗

USGS2024

High resolution identification and quantification of diffuse deep groundwater discharge in mountain rivers using continuous boat-mounted helium measurements

Discharge of deeply sourced groundwater to streams is difficult to locate and quantify, particularly where both discrete and diffuse discharge points exist, but diffuse discharge is one of the primary controls on solute budgets in mountainous watersheds. The noble gas helium is a unique identifier of deep groundwater discharge because groundwater with long residence times is commonly enriched in helium. In this study, a portable mass spectrometer was used to measure longitudinal variation in dissolved helium concentrations in two mountainous rivers at high spatial resolution not feasible with traditional sampling techniques. Helium profiles were then simulated using a mass-balance model to quantify longitudinal variation in groundwater discharge to the receiving rivers. Results indicate helium concentrations were enriched by multiple orders of magnitude above atmospheric equilibrium in both rivers and that this persisted for up to 18 km below observed pulse inputs in the Colorado River. Helium mass-balance models match observed longitudinal patterns with the exception of sharp initial increases in helium observed in the rivers. Increased longitudinal groundwater discharge rates correspond to mapped geologic structures in both watersheds that likely transport deep geothermal water. Models show variable sensitivity to spatial assignment of input variables representing the groundwater source, illustrating the importance of collecting data from discrete groundwater discharges where possible. The methodology shows promise for field experiments designed to assess air–water exchange rates and to quantify total groundwater discharge from a combination of discrete and diffuse sources.

Colorado, Utah↗