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 523 records · Page 29Linked to original sources

USGS2011

The fate and transport of nitrate in shallow groundwater in northwestern Mississippi, USA

Agricultural contamination of groundwater in northwestern Mississippi, USA, has not been studied extensively, and subsurface fluxes of agricultural chemicals have been presumed minimal. To determine the factors controlling transport of nitrate-N into the Mississippi River Valley alluvial aquifer, a study was conducted from 2006 to 2008 to estimate fluxes of water and solutes for a site in the Bogue Phalia basin (1,250 km 2 ). Water-quality data were collected from a shallow water-table well, a vertical profile of temporary sampling points, and a nearby irrigation well. Nitrate was detected within 4.4 m of the water table but was absent in deeper waters with evidence of reducing conditions and denitrification. Recharge estimates from 6.2 to 10.9 cm/year were quantified using water-table fluctuations, a Cl – tracer method, and atmospheric age-tracers. A mathematical advection-reaction model predicted similar recharge to the aquifer, and also predicted that 15% of applied nitrogen is leached into the saturated zone. With current denitrification and application rates, the nitrate-N front is expected to remain in shallow groundwater, less than 6–9 m deep. Increasing application rates resulting from intensifying agricultural demands may advance the nitrate-N front to 16–23 m, within the zone of groundwater pumping.

Mississippi↗

USGS2012

Superposed local and regional paleostresses: fault-slip analysis of Neogene extensional faulting near coeval caldera complexes, Yucca Flat, Nevada

Numerous reduced stress tensors are computed by multiple inversions of 906 temporally and spatially partitioned fault-slip data from the Yucca Flat region in the southwest Nevada volcanic field to constrain the Neogene paleostress and faulting history and to investigate how the regional tectonic stress field was affected by local caldera magmatism. Perturbed, shallow (<400 m), pre-11 Ma paleostress configurations, determined west and northwest of present (post-11 Ma) Yucca Flat basin, existed during mild extensional faulting and are attributed to superposition of transient caldera-magmatic stresses on the regional stress field. Northwest of Yucca Flat a progressive shift in least principal stress (σ 3 ) directions near known calderas located 5–15 km to the west occurred under a normal-slip stress state during caldera development between about 15 and 13 Ma. A brief (∼0.5 m.y.) change to a strike-slip stress state occurred at about 13 Ma and was accompanied by small-offset, quasi-conjugate strike-slip faulting. This stress state was most distinct, relative to a normal-slip state, near calderas where stress solutions and fault relations indicate closer affinities to a reverse-slip state. Inferred 11.6–11.45 Ma paleostress tensors indicate radial tension associated with either initial caldera collapse or local post-collapse topographic modification of the stress field. Post-11 Ma normal-slip stress tensors are associated with normal- and oblique-slip faults that accommodated subsidence and eastward extension of Yucca Flat basin away from the caldera complexes. These tensors do not indicate stress modifications due to residual caldera-related effects and thus were used to infer post-11 Ma regional stress changes. The stress field has rotated as much as 65° clockwise since 11 Ma during extensional development of Yucca Flat basin, with most of the rotation and extension occurring before about 8.5 Ma. Results suggest that shallow magmatism and caldera development can strongly alter extensional tectonic stress fields, fault patterns, and slip directions in the uppermost crust out to distances of roughly two magma chamber radii away from a magma body.

Journal of Geophysical Research↗

USGS2003

Constraints on the composition of ore fluids and implications for mineralising events at the Cleo gold deposit, Eastern Goldfields Province, Western Australia

The Cleo gold deposit, 55 km south of Laverton in the Eastern Goldfields Province of Western Australia, is characterised by banded iron-formation (BIF)-hosted ore zones in the gently dipping Sunrise Shear Zone and high-grade vein-hosted ore in the Western Lodes. There is evidence that gold mineralisation in the Western Lodes (which occurred at ca 2655 Ma) post-dates the majority of displacement along the Sunrise Shear Zone, but it remains uncertain if the ore in both structures formed simultaneously or separately. Overall, the Pb, Nd, Sr, C. O and S isotopic compositions of ore-related minerals from both the Western Lodes and ore zones in the Sunrise Shear Zone are similar. Early low-salinity aqueous-carbonic fluids and late high-salinity fluids with similar characteristics are trapped in inclusions in quartz veins from both the Sunrise Shear Zone and the Western Lodes. The early CO2, CO2-H2O, and H2O- dominant inclusions are interpreted as being related to ore formation, and to have formed from a single low-salinity aqueous-carbonic fluid as a result of intermittent fluid immiscibility. Homogenisation temperatures indicate that these inclusions were trapped at approximately 280??C and at approximately 4 km depth, in the deeper epizonal range. Differences between the ore zones are detected in the trace-element composition of gold samples, with gold from the Sunrise Shear Zone enriched in Ni, Pb, Sn, Te and Zn, and depleted In As, Bi, Cd, Cu and Sb, relative to gold from the Western Lodes. Although there are differences in gold composition between the Sunrise Shear Zone and Western Lodes, and hence the metal content of ore fluids may have varied slightly between the different ore zones, no other systematic fluid or solute differences are detected between the ore zones. Given the fact that the ore fluids in each zone have very similar bulk properties, the considerable differences in gold grade, sulfide mineral abundance, and ore textures between the two ore zones most likely result from different gold-deposition mechanisms. The association of ore zones in the Sunrise Shear Zone with pyrite-replaced BIF suggests that wall-rock sulfidation was the most significant mechanism of gold precipitation, through the destabilisation of gold-bisulfide complexes. The Western Lodes, however, do not exhibit any host-rock preference and multistage veins commonly contain coarse-grained gold. Fluid-inclusion characteristics and breccia textures in veins in the Western Lodes suggest that rapid pressure changes, brought about by intermittent release of overpressured fluids and concomitant phase separation, are likely to have caused the destabilisation of gold-thiocomplexes, leading to formation of higher-grade gold ore zones.

Australian Journal of Earth Sciences↗

USGS2006

Preferential flow estimates to an agricultural tile drain with implications for glyphosate transport

Agricultural subsurface drains, commonly referred to as tile drains, are potentially significant pathways for the movement of fertilizers and pesticides to streams and ditches in much of the Midwest. Preferential flow in the unsaturated zone provides a route for water and solutes to bypass the soil matrix and reach tile drains faster than predicted by traditional displacement theory. This paper uses chloride concentrations to estimate preferential flow contributions to a tile drain during two storms in May 2004. Chloride, a conservative anion, was selected as the tracer because of differences in chloride concentrations between the two sources of water to the tile drain, preferential and matrix flow. A strong correlation between specific conductance and chloride concentration provided a mechanism to estimate chloride concentrations in the tile drain throughout the storm hydrographs. A simple mixing analysis was used to identify the preferential flow component of the storm hydrograph. During two storms, preferential flow contributed 11 and 51% of total storm tile drain flow; the peak contributions, 40 and 81%, coincided with the peak tile drain flow. Positive relations between glyphosate [N-(phosphonomethyl)glycine] concentrations and preferential flow for the two storms suggest that preferential flow is an important transport pathway to the tile drain. ?? ASA, CSSA, SSSA.

Journal of Environmental Quality↗

USGS2014

The behavioural response of adult Petromyzon marinus to damage-released alarm and predator cues

Using semi‐natural enclosures, this study investigated (1) whether adult sea lamprey Petromyzon marinus show avoidance of damage‐released conspecific cues, damage‐released heterospecific cues and predator cues and (2) whether this is a general response to injured heterospecific fishes or a specific response to injured P. marinus . Ten replicate groups of 10 adult P. marinus , separated by sex, were exposed to one of the following nine stimuli: deionized water (control), extracts prepared from adult P. marinus , decayed adult P. marinus (conspecific stimuli), sympatric white sucker Catostomus commersonii , Amazon sailfin catfish Pterygoplichthys pardalis (heterospecific stimuli), 2‐phenylethylamine (PEA HCl) solution, northern water snake Nerodia sipedon washing, human saliva (predator cues) and an adult P. marinus extract and human saliva combination (a damage‐released conspecific cue and a predator cue). Adult P. marinus showed a significant avoidance response to the adult P. marinus extract as well as to C. commersonii , human saliva, PEA and the adult P. marinus extract and human saliva combination. For mobile P. marinus , the N. sipedon washing induced behaviour consistent with predator inspection. Exposure to the P. pardalis extract did not induce a significant avoidance response during the stimulus release period. Mobile adult female P. marinus showed a stronger avoidance behaviour than mobile adult male P. marinus in response to the adult P. marinus extract and the adult P. marinus extract and human saliva combination. The findings support the continued investigation of natural damage‐released alarm cue and predator‐based repellents for the behavioural manipulation of P. marinus populations in the Laurentian Great Lakes.

Journal of Fish Biology↗

USGS2023

Free long wave transformation in the nearshore zone through partial reflections

Long waves play an important role in coastal inundation and shoreline and dune erosion, requiring a detailed understanding of their evolution in nearshore regions and interaction with shorelines. While their generation and dissipation mechanisms are relatively well understood, there are fewer studies describing how reflection processes govern their propagation in the nearshore. We propose a new approach, accounting for partial reflections, which leads to an analytical solution to the free wave linear shallow-water equations at the wave-group scale over general varying bathymetry. The approach, supported by numerical modeling, agrees with the classic Bessel standing solution for a plane sloping beach but extends the solution to arbitrary alongshore uniform bathymetry profiles and decomposes it into incoming and outgoing wave components, which are a combination of successively partially reflected waves lagging each other. The phase lags introduced by partial reflections modify the wave amplitude and explain why Green’s law, which describes the wave growth of free waves with decreasing depth, breaks down in very shallow water. This reveals that the wave amplitude at the shoreline is highly dependent on partial reflections. Consistent with laboratory and field observations, our analytical model predicts a reflection coefficient that increases and is highly correlated with the normalized bed slope (bed slope relative to wave frequency). Our approach shows that partial reflections occurring due to depth variations in the nearshore are responsible for the relationship between the normalized bed slope and the amplitude of long waves in the nearshore, with direct implications for determining long-wave amplitudes at the shoreline and wave runup.

Journal of Physical Oceanography↗

USGS1987

Analysis of an anisotropic coastal aquifer system using variable-density flow and solute transport simulation

The groundwater system in southern Oahu, Hawaii consists of a thick, areally extensive freshwater lens overlying a zone of transition to a thick saltwater body. This system is analyzed in cross section with a variable-density groundwater flow and solute transport model on a regional scale. The simulation is difficult, because the coastal aquifer system has a saltwater transition zone that is broadly dispersed near the discharge area, but is very sharply defined inland. Steady-state simulation analysis of the transition zone in the layered basalt aquifer of southern Oahu indicates that a small transverse dispersivity is characteristic of horizontal regional flow. Further, in this system flow is generally parallel to isochlors and steady-state behavior is insensitive to the longitudinal dispersivity. Parameter analysis identifies that only six parameters control the complex hydraulics of the system: horizontal and vertical hydraulic conductivity of the basalt aquifer; hydraulic conductivity of the confining "caprock" layer; leakance below the caprock; specific yield; and aquifer matrix compressibility. The best-fitting models indicate the horizontal hydraulic conductivity is significantly greater than the vertical hydraulic conductivity. These models give values for specific yield and aquifer compressibility which imply a considerable degree of compressive storage in the water table aquifer.

Hawaii↗

USGS2013

Uncertainty in simulated groundwater-quality trends in transient flow

In numerical modeling of groundwater flow, the result of a given solution method is affected by the way in which transient flow conditions and geologic heterogeneity are simulated. An algorithm is demonstrated that simulates breakthrough curves at a pumping well by convolution-based particle tracking in a transient flow field for several synthetic basin-scale aquifers. In comparison to grid-based (Eulerian) methods, the particle (Lagrangian) method is better able to capture multimodal breakthrough caused by changes in pumping at the well, although the particle method may be apparently nonlinear because of the discrete nature of particle arrival times. Trial-and-error choice of number of particles and release times can perhaps overcome the apparent nonlinearity. Heterogeneous aquifer properties tend to smooth the effects of transient pumping, making it difficult to separate their effects in parameter estimation. Porosity, a new parameter added for advective transport, can be accurately estimated using both grid-based and particle-based methods, but predictions can be highly uncertain, even in the simple, nonreactive case.

Hydrogeology Journal↗

USGS2017

Leaching and sorption of neonicotinoid insecticides and fungicides from seed coatings

Seed coatings are a treatment used on a variety of crops to improve production and offer protection against pests and fungal outbreaks. The leaching of the active ingredients associated with the seed coatings and the sorption to soil was evaluated under laboratory conditions using commercially available corn and soybean seeds to study the fate and transport of these pesticides under controlled conditions. The active ingredients (AI) included one neonicotinoid insecticide (thiamethoxam) and five fungicides (azoxystrobin, fludioxonil, metalaxyl, sedaxane thiabendazole). An aqueous leaching experiment was conducted with treated corn and soybean seeds. Leaching potential was a function of solubility and seed type. The leaching of fludioxonil, was dependent on seed type with a shorter time to equilibrium on the corn compared to the soybean seeds. Sorption experiments with the treated seeds and a solution of the AIs were conducted using three different soil types. Sorption behavior was a function of soil organic matter as well as seed type. For most AIs, a negative relationship was observed between the aqueous concentration and the log K oc . Sorption to all soils tested was limited for the hydrophilic pesticides thiamethoxam and metalaxyl. However, partitioning for the more hydrophobic fungicides was dependent on both seed type and soil properties. The mobility of fludioxonil in the sorption experiment varied by seed type indicating that the adjuvants associated with the seed coating could potentially play a role in the environmental fate of fludioxonil. This is the first study to assess, under laboratory conditions, the fate of pesticides associated with seed coatings using commercially available treated seeds. This information can be used to understand how alterations in agricultural practices (e.g., increasing use of seed treatments) can impact the exposure (concentration and duration) and potential effects of these chemicals to aquatic and terrestrial organisms.

Journal of Environmental Science and Health - Part↗

USGS2016

The removal kinetics of dissolved organic matter and the optical clarity of groundwater

Concentrations of dissolved organic matter (DOM) and ultraviolet/visible light absorbance decrease systematically as groundwater moves through the unsaturated zones overlying aquifers and along flowpaths within aquifers. These changes occur over distances of tens of meters (m) implying rapid removal kinetics of the chromophoric DOM that imparts color to groundwater. A one-compartment input-output model was used to derive a differential equation describing the removal of DOM from the dissolved phase due to the combined effects of biodegradation and sorption. The general solution to the equation was parameterized using a 2-year record of dissolved organic carbon (DOC) concentration changes in groundwater at a long-term observation well. Estimated rates of DOC loss were rapid and ranged from 0.093 to 0.21 micromoles per liter per day (μM d −1 ), and rate constants for DOC removal ranged from 0.0021 to 0.011 per day (d −1 ). Applying these removal rate constants to an advective-dispersion model illustrates substantial depletion of DOC over flow-path distances of 200 m or less and in timeframes of 2 years or less. These results explain the low to moderate DOC concentrations (20–75 μM; 0.26–1 mg L −1 ) and ultraviolet absorption coefficient values ( a 254 < 5 m −1 ) observed in groundwater produced from 59 wells tapping eight different aquifer systems of the United States. The nearly uniform optical clarity of groundwater, therefore, results from similarly rapid DOM-removal kinetics exhibited by geologically and hydrologically dissimilar aquifers.

California, Colorado, Connecticut, Georgia, Illino↗

USGS2024

Improving crop-specific groundwater use estimation in the Mississippi Alluvial Plain: Implications for integrated remote sensing and machine learning approaches in data-scarce regions

Study region The Mississippi Alluvial Plain (MAP) in the United States (US). Study focus Understanding local-scale groundwater use, a critical component of the water budget, is necessary for implementing sustainable water management practices. The MAP is one of the most productive agricultural regions in the US and extracts more than 11 km 3 /year for irrigation activities. Consequently, groundwater-level declines in the MAP region pose a substantial challenge to water sustainability, and hence, we need reliable groundwater pumping monitoring solutions to manage this resource appropriately. New hydrological insights for the region We incorporate remote sensing datasets and machine learning to improve an existing lookup table-based model of groundwater use previously developed by the U.S. Geological Survey (USGS). Here, we employ Distributed Random Forests, an ensemble machine learning algorithm to predict annual and monthly groundwater use (2014–2020) throughout this region at 1-km resolution, using pumping data from existing flowmeters in the Mississippi Delta. Our model compares favorably with the existing USGS model, with higher R 2 (0.51 compared to 0.42 in the previous model), and lower root mean square error (RMSE) and mean absolute error (MAE)— 0.14 m and 0.09 m, respectively in our model, compared to 0.15 m and 0.1 m in the previous model. Therefore, this work advances our ability to predict groundwater use in regions with scarce or limited in-situ groundwater withdrawal data availability.

Journal of Hydrology Regional Studies↗

USGS2017

A probabilistic approach to remote compositional analysis of planetary surfaces

Reflected light from planetary surfaces provides information, including mineral/ice compositions and grain sizes, by study of albedo and absorption features as a function of wavelength. However, deconvolving the compositional signal in spectra is complicated by the nonuniqueness of the inverse problem. Trade-offs between mineral abundances and grain sizes in setting reflectance, instrument noise, and systematic errors in the forward model are potential sources of uncertainty, which are often unquantified. Here we adopt a Bayesian implementation of the Hapke model to determine sets of acceptable-fit mineral assemblages, as opposed to single best fit solutions. We quantify errors and uncertainties in mineral abundances and grain sizes that arise from instrument noise, compositional end members, optical constants, and systematic forward model errors for two suites of ternary mixtures (olivine-enstatite-anorthite and olivine-nontronite-basaltic glass) in a series of six experiments in the visible-shortwave infrared (VSWIR) wavelength range. We show that grain sizes are generally poorly constrained from VSWIR spectroscopy. Abundance and grain size trade-offs lead to typical abundance errors of ≤1 wt % (occasionally up to ~5 wt %), while ~3% noise in the data increases errors by up to ~2 wt %. Systematic errors further increase inaccuracies by a factor of 4. Finally, phases with low spectral contrast or inaccurate optical constants can further increase errors. Overall, typical errors in abundance are <10%, but sometimes significantly increase for specific mixtures, prone to abundance/grain-size trade-offs that lead to high unmixing uncertainties. These results highlight the need for probabilistic approaches to remote determination of planetary surface composition.

Journal of Geophysical Research E: Planets↗

USGS1973

A calorimetric determination of the standard enthalpies of formation of huntite, CaMg3 (CO3)4 , and artinite, Mg2(OH)2 CO3 * 3H2O, and their standard Gibbs free energies of formation

The enthalpies of formation, ΔH° f, of huntite, CaMg 3 (CO 3 ) 4 , and artinite, Mg 2 (OH) 2 CO 3 * 3H 2 O, have been determined by HCl solution calorimetry using a constant-volume isoperibol reaction calorimeter. For the reaction CaO(c) + 3MgO(c) + 4CO 2 (g) = CaMg 3 (CO 3 ) 4 (c), the enthalpy change at 298.15 K, ΔH° 298 , is -123,203±145 cal mol -1 . For the reaction 2MgO(c) + 4H 2 O(l) + CO 2 (g)=Mg 2 (OH) 2 CO 3 * 3H 2 O(c), we obtained -45,132±100 cal mol -1 . These results combined with the standard enthalpies of formation of CaO, MgO, H 2 O, and CO 2 lead to ΔH° 298 (huntite) = -1,082,600±375 cal mol -1 and ΔH° 298 (artinite) = -698,043±170 cal mol -1 . Using recently determined values for the standard entropies of huntite, CaMg 3 (CO 3 ) 4 , and artinite, Mg 2 (OH) 2 CO 3 * 3H 2 O, and of Mg, Ca, C, O 2 , and H 2 , we calculate ΔG° f,298 (huntite) = -1,004,707±390 cal mol -1 and ΔG° f,298 (artinite) = -613,924±180 cal mol -1 .

Journal of Research of the U.S. Geological Survey↗

USGS2015

Peclet number as affected by molecular diffusion controls transient anomalous transport in alluvial aquifer-aquitard complexes

This study evaluates the role of the Peclet number as affected by molecular diffusion in transient anomalous transport, which is one of the major knowledge gaps in anomalous transport, by combining Monte Carlo simulations and stochastic model analysis. Two alluvial settings containing either short- or long-connected hydrofacies are generated and used as media for flow and transport modeling. Numerical experiments show that 1) the Peclet number affects both the duration of the power-law segment of tracer breakthrough curves (BTCs) and the transition rate from anomalous to Fickian transport by determining the solute residence time for a given low-permeability layer, 2) mechanical dispersion has a limited contribution to the anomalous characteristics of late-time transport as compared to molecular diffusion due to an almost negligible velocity in floodplain deposits, and 3) the initial source dimensions only enhance the power-law tail of the BTCs at short travel distances. A tempered stable stochastic (TSS) model is then applied to analyze the modeled transport. Applications show that the time-nonlocal parameters in the TSS model relate to the Peclet number, P e . In particular, the truncation parameter in the TSS model increases nonlinearly with a decrease in P e due to the decrease of the mean residence time, and the capacity coefficient increases with an increase in molecular diffusion which is probably due to the increase in the number of immobile particles. The above numerical experiments and stochastic analysis therefore reveal that the Peclet number as affected by molecular diffusion controls transient anomalous transport in alluvial aquifer–aquitard complexes.

Journal of Contaminant Hydrology↗

USGS1996

Chemical and isotopic characteristics of fluids along the cameroon volcanic line, cameroon

Results of the chemical and isotopic analysis of the water and gases discharged from volcanic crater lakes and soda springs located along the Cameroon Volcanic Line were used to characterize and infer their genetic relationships. Variations in the solute compositions of the waters indicate the dominant influence of silicate hydrolysis. Na+ (40-95%) constitutes the major cation in the springs while Fe2+ + Mg2+ (70%) dominate in the CO2-rich lakes. The principal anion is HCO3 (>90%), except in the coastal springs where Cl predominates. Lakes Nyos and Monoun have Fe-Mg-Ca-HCO3 type signatures; the soda springs are essentially Na-HCO3 type, while all other lakes show similar ionic compositions to dilute surface waters. Dissolved gases show essentially CO2 (>90%), with small amounts of Ar and N2, while CH4 constitutes the principal component in the non-gassy lakes. Active volcanic gases are generally absent, except in the Lobe spring with detectable H2S. Stable isotope ratio evidence indicates that the bicarbonate waters are essentially of meteoric origin. CO2 (??13C = -2 to -8???) and He (3He/4He = 1 to 5.6Ra) infer a mantle contribution to the total CO2. CH4 has a biogenic source, while Ar and N2 are essentially atmospheric in origin, but mixing is quite common.

Journal of African Earth Sciences↗

USGS2013

Source mechanisms of persistent shallow earthquakes during eruptive and non-eruptive periods between 1981 and 2011 at Mount St. Helens, Washington

Shallow seismicity between 0 and 3-km depth has persisted at Mount St. Helens, Washington (MSH) during both eruptive and non-eruptive periods for at least the past thirty years. In this study we investigate the source mechanisms of shallow volcano-tectonic (VT) earthquakes at MSH by calculating high-quality hypocenter locations and fault plane solutions (FPS) for all VT events recorded during two eruptive periods (1981–1986 and 2004–2008) and two non-eruptive periods (1987–2004 and 2008–2011). FPS show a mixture of normal, reverse, and strike-slip faulting during all periods, with a sharp increase in strike-slip faulting observed in 1987–1997 and an increase in normal faulting in 1998–2004. FPS P-axis orientations show a ~ 90° rotation with respect to regional σ1 (N23°E) during 1981–1986 and 2004–2008, bimodal orientations (~ N-S and ~ E-W) during 1987–2004, and bimodal orientations at ~ N-E and ~ S-W from 2008–2011. We interpret these orientations to likely be due to pressurization accompanying the shallow intrusion and subsequent eruption of magma as domes during 1981–1986 and 2004–2008 and the buildup of pore pressure beneath a seismogenic volume (located at 0–1 km) with a smaller component due to the buildup of tectonic forces during 1987–2004 and 2008–2011.

Washington↗

USGS2000

Estimating formation properties from early-time oscillatory water levels in a pumped well

Hydrologists often attempt to estimate formation properties from aquifer tests for which only the hydraulic responses in a pumped well are available. Borehole storage, turbulent head losses, and borehole skin, however, can mask the hydraulic behavior of the formation inferred from the water level in the pumped well. Also, in highly permeable formations or in formations at significant depth below land surface, where there is a long column of water in the well casing, oscillatory water levels may arise during the onset of pumping to further mask formation responses in the pumped well. Usually borehole phenomena are confined to the early stages of pumping or recovery, and late-time hydraulic data can be used to estimate formation properties. In many instances, however, early-time hydraulic data provide valuable information about the formation, especially if there are interferences in the late-time data. A mathematical model and its Laplace transform solution that account for inertial influences and turbulent head losses during pumping is developed for the coupled response between the pumped borehole and the formation. The formation is assumed to be homogeneous, isotropic, of infinite areal extent, and uniform thickness, with leakage from an overlying aquifer, and the screened or open interval of the pumped well is assumed to fully penetrate the pumped aquifer. Other mathematical models of aquifer flow can also be coupled with the equations describing turbulent head losses and the inertial effects on the water column in the pumped well. The mathematical model developed in this paper is sufficiently general to consider both underdamped conditions for which oscillations arise, and overdamped conditions for which there are no oscillations. Through numerical inversion of the Laplace transform solution, type curves from the mathematical model are developed to estimate formation properties through comparison with the measured hydraulic response in the pumped well. The mathematical model is applied to estimate formation properties from a singlewell test conducted near Waialua, Oahu, Hawaii. At this site, both the drawdown and recovery showed oscillatory water levels in the pumped well, and a step-drawdown test showed that approximately 86% of the drawdown is attributed to turbulent head losses. Analyses at this site using late-time drawdown data were confounded by the noise present in the measured water levels due primarily to nearby irrigation wells and ocean tides. By analyzing the early-time oscillatory recovery data at the Waialua site, upper and lower bounds were placed on the transmissivity, T, storage coefficient, S, and the leakance of the confining unit, K′/B′. The upper and lower bounds on T differ by a factor of 2. Upper and lower bounds on S and K′/B′ are much larger, because drawdown stabilized relatively quickly after the onset of pumping.

Journal of Hydrology↗

USGS2025

Mitigating flood risks in urban estuaries: Tidal dynamics, shoreline hardening, nature-based solutions, and floodgates in San Francisco Bay

Hydrodynamic models are valuable tools for understanding the primary factors influencing daily and peak water levels and for guiding discussions on potential adaptation strategies for managing flood risk in coastal areas. This analysis uses the Delft3D San Francisco Bay-Delta Community Model to simulate water levels and incorporates the effects of a number of adaptation measures in the urban San Francisco Bay estuary, California. In particular, we examine the influence of shoreline hardening, nature-based solutions, and subregional floodgates on regional water levels. The result shows that under present conditions, tidal amplification is responsible for generating a wide distribution of extreme water levels across San Francisco Bay. Tidal amplification is found to decrease under sea level rise, thereby producing a relative damping effect on extremes. A comparison of different shoreline scenarios demonstrates that hard frontal shorelines result in higher tidal amplification, whereas restored (soft) shorelines lower amplification. The current shoreline configuration has both hard and soft characteristics and results in an intermediate tidal response. In some areas, wetland restoration reduces extreme water levels by as much as 20 cm, whereas hard-shoreline addition elevates them by as much as 10 cm for 1.5 m of sea level rise. Furthermore, local floodgates can significantly reduce high water levels without major adverse effects elsewhere in San Francisco Bay. These findings point toward the justification for a range of adaptive measures across political boundaries, weighing hard and soft options in addressing the mounting danger of sea level rise.

California↗