USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Hydrology”

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

Ecohydrological response of a forested headwater catchment to a flash drought in the Southeastern U.S.

Flash droughts differ from traditionally defined droughts in their rapidity of intensification and often associated high vapor-pressure deficit. These droughts can lead to declines in streamflow and water table depth and induce water stress to vegetation at a greater rate than droughts that manifest over longer periods. However, little is known regarding the response of forested environments to flash drought because most studies of impacts have been conducted in agricultural settings. In this study we investigated water-use patterns of riparian trees using sap flow methods and examined the role of groundwater as a source of moisture over three periods that were delimited by antecedent soil moisture conditions. For a longer-term perspective we also examine monthly streamflow over the 35-year record. We observed that trees at only one monitoring plot showed a decrease in water use relative to evaporative demand during a flash drought. Total reverse sap flow (flow toward the roots rather than the canopy) greatly increased during the flash drought period, suggesting the likely occurrence of hydraulic redistribution to the excessively dry soils. Over the drought period groundwater became a more dominant source of moisture for sustaining forest water use. Monthly mean streamflow during the flash drought approached levels observed in past multiyear droughts. This is the first study, to our knowledge, to specifically investigate the response of multiple water budget components to flash drought in a humid forest. As more studies are conducted, a better understanding of the range of expected responses are likely to emerge.

Georgia↗
USGS2025

Ranking river basins for stream temperature research and monitoring in the contiguous United States

There is a need to prioritize research and data collection in river basins by integrating information from environmental, ecological, and socioeconomic datasets to maintain acceptable water quality for human uses and ecosystem health. Multiple anthropogenic and natural stressors are responsible for driving changes in stream temperatures that can alter ecosystems and degrade water quality. These stressors are variable spatially and temporally, which can be challenging for prioritizing monitoring and research. In this study, an evaluation of variables related to stream temperature was performed for 163 candidate basins of the contiguous United States to highlight potential focal areas. Thirty variables were combined in six components: anthropogenic forcing, natural forcing, climate change, ecological sensitivity, socioeconomic sensitivity, and data availability. The components were incorporated into three themes: vulnerability of streams to thermal change, vulnerability of ecosystems, and vulnerability of communities. By evaluating each theme separately, patterns of vulnerability and potential resiliency were identified as well as consistency in ranks within the themes. For example, results of the national scale ranking indicated the Southern California Coastal basin (within California-Nevada region) was the highest-ranking priority in vulnerability of streams to thermal change and vulnerability of ecosystems. The analysis also identified vulnerable basins with gaps in monitoring. For example, the Missouri-Oahe basin (within the Northern High Plains region) was the highest-ranking priority for vulnerable communities with limitations in data availability. The ranking approach provides insight into basins that are resilient and are ideal candidates for monitoring and research.

Journal of Hydrology↗
USGS2025

A joint Gaussian process model of geochemistry, geophysics, and temperature for groundwater TDS in the San Ardo Oil Field, California, USA

Decline in availability of fresh groundwater has expanded interest in brackish groundwater resources; however, the distribution of brackish groundwater is poorly understood. Water resources in sedimentary basins across the United States often overlie oil and gas development. Mapping of groundwater total dissolved solids (TDS) using data from oil well geophysical logs has become an important technique for identifying fresh and brackish groundwater. Existing geophysical log analysis methods use porosity and temperature to relate formation resistivity to TDS. Typically, natural geothermal gradients are used to estimate temperature at the location of collected resistivity. However, in thermally enhanced oil fields, steam is injected into the subsurface to mobilize high viscosity oil, creating variable temperature distributions. Furthermore, TDS derived from resistivity also depends on the fractions of dominant ions. Typically, chloride and bicarbonate fractions must be determined. It is also necessary to model TDS across many geologic units with heterogenous porosity distributions. Collectively, each quantity used to estimate TDS (resistivity, porosity, temperature, bicarbonate fraction) varies in space and time, and available data points are rarely collocated. Here, we present a new method of mapping groundwater TDS that continuously models each quantity together with a joint Gaussian process. This method enables mapping fresh and brackish water with practically available data. We apply this method to the San Ardo Oil Field in Monterey County, California, where steam injection occurs. In some areas of the aquifer system overlying the oil zone, the temperature is ∼75 °C, roughly twice the natural background value. Groundwater TDS is typically <1,500 mg/L in the aquifer and increases with depth to ∼9,000 mg/L in the oil-producing zone. A low-permeability clay layer delineates the fresh and brackish water, likely by inhibiting surface recharge from penetrating the deeper zones, allowing higher-TDS connate water to remain in place. Weaker lateral TDS trends may be controlled by recharge patterns associated with the Salinas River. Our model reveals with high certainty that groundwater has freshened in one localized part of the oil-producing zone and suggests with less certainty that more widespread freshening has also occurred. The lowering of TDS was possibly from decades of low-TDS steam injection and the associated fluid production and disposal operations.

California↗
USGS2025

Tracking persistent declines in suspended sediment in the Lower Mississippi and Atchafalaya Rivers, 1992–2021: Harnessing WRTDSplus to characterize longitudinally varying trends and explore connections to streamflow

Suspended sediment (SS) continues a century-long decline in the Lower Mississippi and Atchafalaya Rivers, United States. In this study, we use the WRTDSplus model to estimate concentrations and loads for total, fine (<0.0625 millimeter (mm)), and coarse (≥0.0625 mm) SS for 11 sites. This extension of the Weighted Regressions on Time, Discharge, and Season (WRTDS) model allows a fourth explanatory variable in the model formulation. We incorporated hysteresis terms for most models based on a residual analysis, which allowed for the identification of decreased flushing over time at some sites. Total, fine, and coarse SS concentrations and loads decreased at all sites over two trend periods (water years (WY) 1992–2021 and WY 2012–2021). Declines were largely due to changes in fine SS (mud and silt) but decreases in coarse SS (sands) were also widespread. On average, recent declines are more severe in the Lower Mississippi River below the Old River Control Complex (ORCC, −3.7 mg per liter per year (mg/L/yr)) compared to the Atchafalaya River (−2.0 mg/L/yr), although there is longitudinal variability within each river. The reach below the ORCC is a net SS sink, leading to complex temporal changes for the sites in this area. Streamflows (low, moderate, and high) have increased over these periods, with the last decade being particularly wet. Increasing streamflow and decreasing SS, with little evidence of amelioration, may influence spillway operations during floods, sediment diversion construction and operation, coastal restoration efforts, and aquatic health.

Louisiana, Mississippi↗
USGS2025

A spatiotemporal deep learning approach for predicting daily air-water temperature signal coupling and identification of key watershed physical parameters in a montane watershed

Seasonal shifts from runoff to groundwater dominance influence daily headwater stream temperatures, especially where local groundwater input is strong. This input buffers temperature during hot periods, supporting cold-water habitats. Recent studies use air–water temperature signal metrics to identify zones of strong stream–groundwater connectivity. While Previous studies used air–water signal ratios as proxies for groundwater influence but were limited to specific sites and periods, without dynamic forecasting. This study is the first to forecast daily A r as a spatiotemporal signal using a Graph Convolutional Network–Long Short-Term Memory (GCN-LSTM) model. The model was trained using hydroclimate data (air temperature, precipitation, shortwave radiation, streamflow) and watershed physical features (e.g., sand content, slope). Results showed high predictive skill, achieving R 2 (NSE, RMSE) of 0.86 (0.73, 0.0004) for one-day-ahead to 0.52 (0.50, 0.0009) for seven-days ahead forecasts. Prior studies often have not explicitly incorporated spatial hydrogeologic drivers, but this model explicitly incorporates them to assess their impact on A r forecasting and stream-groundwater connectivity. Feature analysis identified mean sand, elevation, slope, clay, and TWI as key predictors of A r . Stronger groundwater signals appeared in hillslopes, elevations, and tributaries, highlighting watershed influence on streamflow. However, limitations include reliance on historical air–water temperature patterns for training and limited representation of extreme climate conditions. Despite these limitations, unlike previous studies relying on measured in-situ stream and air temperature, this study forecasts A r directly from climate and physiographic features after training, avoiding in-situ data requirements. Findings aiding predictions of stream ecosystem resilience.

New York↗
USGS2025

Multi-lithofacies alluvial characterization via airborne electromagnetic-borehole fusion using ordinary interval kriging and geologic constraints

Airborne electromagnetic (AEM) survey provides extensive spatial coverage and detailed resolution at the near surface and can be used to develop hydrogeological models. However, utilization of AEM data is not straightforward because AEM resistivity is an indirect measurement for inferring sediment types. This study develops an ordinary interval kriging (OIK) algorithm and a resistivity-to-multi-lithofacies (R2ML) data fusion workflow for multi-lithofacies alluvial characterization. OIK utilizes irregular interval data to construct three-dimensional (3D) resistivity fields from one-dimensional inverted AEM resistivity models. The R2ML workflow maps the resistivity field generated from OIK into a multi-facies lithological model, incorporating geologic constraints derived from well logs and geological observations. The numerical and real-world cases demonstrate that OIK is computationally efficient, accounts for 3D anisotropy, and minimizes the smoothing effect, thereby preserving resistivity contrasts and reducing interpolation uncertainty. The methodology is applied to lithologic characterization of the Mississippi River Valley alluvial aquifer (MRVA) in the Shellmound area, Mississippi, U.S. A frequency-domain AEM survey was conducted to support groundwater studies for the managed aquifer recharge (MAR) to the MRVA. The resulting lithological model, including four types of lithofacies—clay, very fine sand, fine-medium sands, and graveliferous sands, illustrates the geomorphological processes of the MRVA and implies potential MAR. The alignment between the lithological model and existing geological and hydrogeological investigations demonstrates that OIK and R2ML workflow effectively capture the subsurface architecture of the MRVA. The methods have broad applicability for characterizing alluvial aquifers through AEM-borehole data fusion, supporting sustainable groundwater management.

Mississippi↗
USGS2026

Groundwater spatial variability within an atoll island: Assessing shallow aquifer heterogeneity with geophysical and physicochemical measurements

This study examines the spatial variability of shallow groundwater on Dhigelaabadhoo Island using electromagnetic induction surveys, groundwater monitoring, and sediment analyses. The research reveals how variations in island morphology—such as differences in elevation, reef flat width, and sediment composition—affect the spatial distribution of groundwater lenses and the overall aquifer dynamics. Saltwater intrusion is especially pronounced in low elevated areas, with narrow reef flat plate and areas where higher hydraulic conductivity—driven by the presence of coarser sediments—is observed, whereas regions characterized by finer sediments, higher elevation, and wider reef flat plates tend to support more symmetric and less saline groundwater lenses. The geophysical investigations reveal that tidal oscillations alter groundwater movement by markedly changing water levels and conductivity, thereby underscoring the critical need to account for temporal dynamics in atoll coastal aquifer systems and the importance of integrating tidal dynamics into the aquifer zone. The findings highlight the significant role of intrinsic morphological and external hydrodynamic factors in shaping groundwater distribution on atoll islands, offering critical insights for sustainable freshwater resource management.

Dhigelaabadhoo Island↗
USGS2026

Groundwater drought in the United States: Spatial and temporal variability

Many communities and ecosystems in the United States that are dependent on groundwater are potentially adversely affected by groundwater drought. We computed yearly groundwater-drought metrics and mean groundwater levels at well locations across the conterminous United States (CONUS), using data from wells and remotely sensed and modeled Gravity Recovery and Climate Experiment Drought Monitor Data Assimilation (GRACE-DADM). We also modeled the probability of low or high human impact at each well location. The spatial distribution of groundwater-drought duration and severity from 2001 to 2020 for 1,510 wells shows longer maximum duration and higher maximum severity events in drier regions like the Southwest than in wetter regions like the Northeast. Based on 613 wells in CONUS from 1981 to 2020, there are many significant decreases in drought duration and severity in the Northeast and many significant increases in annual-mean groundwater levels. In contrast, there are many significant increases in drought metrics and decreases in mean water levels in parts of the Southeast. There are major differences in trends from 2001 to 2020 between well-based and GRACE-DADM-based groundwater metrics in some CONUS regions and a very low correlation between trends at individual locations across CONUS. A potential reason for this disparity is the low GRACE-DADM resolution (∼12 km) and the potential for a large amount of groundwater variation at the local scale. Also, GRACE-DADM represents shallow, unconfined aquifers which may not match the screened interval of the monitoring wells we evaluated. Large spatial gaps in long-term, high frequency, and quality-assured groundwater-well monitoring data present a challenge for understanding groundwater-drought variability across CONUS. Remote sensing tools such as GRACE can help but cannot fully replace well monitoring, as highlighted by our study results. Substantially more long-term monitoring wells would more accurately represent groundwater-drought trends and spatial variability across CONUS, particularly in western regions.

conterminous United States↗
USGS2026

Reservoir releases and land cover interact to drive event-scale nitrate export in a large agricultural basin

Understanding the drivers of nitrate export in rivers is critical for developing effective nutrient management strategies. However, few studies have explored event-scale drivers of export in large river basins with human modifications like reservoirs. Here, we analyzed nitrate concentration-discharge (C-Q) relationships from 215 events at the outlet of the Kansas River Basin, USA (155,690 km 2 ) from 2014 to 2022 to (1) characterize event-scale nitrate export behaviors in a large agricultural basin, and (2) determine how different event characteristics are linked to these behaviors. We found that C-Q behaviors varied greatly, with 60% of events producing nitrate enrichment (n = 130) and 40% of events producing nitrate dilution (n = 85). These behaviors were correlated with complex spatial interlinkages between climate and land cover: across most of the basin, nitrate enrichment was correlated with drier antecedent conditions, but in wetter areas with higher proportions of urban and forested land cover, enrichment was more strongly correlated with precipitation magnitude/intensity. This difference in hydroclimatic controls on nitrate export might be related to differential distributions of nitrate sources within these land covers. Upstream of major reservoirs, however, neither variable was strongly correlated with C-Q behavior, suggesting that nitrate attenuation within reservoirs decouples event-scale concentration signals in upstream waters from those downstream. Reservoir outflows had variable impacts on C-Q behavior, reflecting reservoir-specific variations in nitrate attenuation efficiency. Together, these results identify specific complex interactions between hydroclimate, land cover, reservoir positioning, and individual reservoir properties that control event-scale nitrate export from large basins.

Colorado, Kansas, Nebraska↗
USGS2026

Drone–based camera discharge derived from the entropy–based probability concept and image velocimetry algorithms

In June and August 2023, the U.S. Geological Survey (USGS) conducted small Uncrewed Aircraft System (sUAS or drone) flights at five streamgages in the Upper Colorado River Basin, Colorado, USA. The objectives of this research were (1) compare maximum surface velocities and discharges derived from Space-time Image Velocimetry (STIV) and Particle Image Velocimetry (PIV) and (2) evaluate drone-based camera algorithms that rely on the entropy-based probability concept and mid-section method relative to standard discharge measurements. To evaluate the scalability and transferability of the methods, river reaches were selected with varying top widths, depths, and discharges and colocated with existing USGS streamgages. Percent errors between drone-based camera algorithms and hydroacoustic measurements were substantially smaller for large rivers such as the Colorado River, Colorado when compared to small rivers such as the Blue River, Colorado. A Wilcoxon signed-rank test suggests there was no difference between the discharges generated by STIV and PIV algorithms and hydroacoustic measurements. This implies that, if the appropriate algorithms are properly parameterized, drone-based camera discharges are consistent with standard discharge measurements when siting is conducted appropriately.

Colorado↗
USGS2000

Precipitation areal-reduction factor estimation using an annual-maxima centered approach

The adjustment of precipitation depth of a point storm to an effective (mean) depth over a watershed is important for characterizing rainfall-runoff relations and for cost-effective designs of hydraulic structures when design storms are considered. A design storm is the precipitation point depth having a specified duration and frequency (recurrence interval). Effective depths are often computed by multiplying point depths by areal-reduction factors (ARF). ARF range from 0 to 1, vary according to storm characteristics, such as recurrence interval; and are a function of watershed characteristics, such as watershed size, shape, and geographic location. This paper presents a new approach for estimating ARF and includes applications for the 1-day design storm in Austin, Dallas, and Houston, Texas. The approach, termed 'annual-maxima centered,' specifically considers the distribution of concurrent precipitation surrounding an annual-precipitation maxima, which is a feature not seen in other approaches. The approach does not require the prior spatial averaging of precipitation, explicit determination of spatial correlation coefficients, nor explicit definition of a representative area of a particular storm in the analysis. The annual-maxima centered approach was designed to exploit the wide availability of dense precipitation gauge data in many regions of the world. The approach produces ARF that decrease more rapidly than those from TP-29. Furthermore, the ARF from the approach decay rapidly with increasing recurrence interval of the annual-precipitation maxima. (C) 2000 Elsevier Science B.V.The adjustment of precipitation depth of a point storm to an effective (mean) depth over a watershed is important for characterizing rainfall-runoff relations and for cost-effective designs of hydraulic structures when design storms are considered. A design storm is the precipitation point depth having a specified duration and frequency (recurrence interval). Effective depths are often computed by multiplying point depths by areal-reduction factors (ARF). ARF range from 0 to 1, vary according to storm characteristics, such as recurrence interval; and are a function of watershed characteristics, such as watershed size, shape, and geographic location. This paper presents a new approach for estimating ARF and includes applications for the 1-day design storm in Austin, Dallas, and Houston, Texas. The approach, termed 'annual-maxima centered,' specifically considers the distribution of concurrent precipitation surrounding an annual-precipitation maxima, which is a feature not seen in other approaches. The approach does not require the prior spatial averaging of precipitation, explicit determination of spatial correlation coefficients, nor explicit definition of a representative area of a particular storm in the analysis. The annual-maxima centered approach was designed to exploit the wide availability of dense precipitation gauge data in many regions of the world. The approach produces ARF that decrease more rapidly than those from TP-29. Furthermore, the ARF from the approach decay rapidly with increasing recurrence interval of the annual-precipitation maxima.

Journal of Hydrology↗
USGS2000

Aquifer response to stream-stage and recharge variations. I. Analytical step-response functions

Laplace transform step-response functions are presented for various homogeneous confined and leaky aquifer types and for anisotropic, homogeneous unconfined aquifers interacting with perennial streams. Flow is one-dimensional, perpendicular to the stream in the confined and leaky aquifers, and two-dimensional in a plane perpendicular to the stream in the water-table aquifers. The stream is assumed to penetrate the full thickness of the aquifer. The aquifers may be semi-infinite or finite in width and may or may not be bounded at the stream by a semipervious streambank. The solutions are presented in a unified manner so that mathematical relations among the various aquifer configurations are clearly demonstrated. The Laplace transform solutions are inverted numerically to obtain the real-time step-response functions for use in the convolution (or superposition) integral. To maintain linearity in the case of unconfined aquifers, fluctuations in the elevation of the water table are assumed to be small relative to the saturated thickness, and vertical flow into or out of the zone above the water table is assumed to occur instantaneously. Effects of hysteresis in the moisture distribution above the water table are therefore neglected. Graphical comparisons of the new solutions are made with known closed-form solutions.Laplace transform step-response functions are presented for various homogeneous confined and leaky aquifer types and for anisotropic, homogeneous unconfined aquifers interacting with perennial streams. Flow is one-dimensional, perpendicular to the stream in the confined and leaky aquifers, and two-dimensional in a plane perpendicular to the stream in the water-table aquifers. The stream is assumed to penetrate the full thickness of the aquifer. The aquifers may be semi-infinite or finite in width and may or may not be bounded at the stream by a semipervious streambank. The solutions are presented in a unified manner so that mathematical relations among the various aquifer configurations are clearly demonstrated. The Laplace transform solutions are inverted numerically to obtain the real-time step-response functions for use in the convolution (or superposition) integral. To maintain linearity in the case of unconfined aquifers, fluctuations in the elevation of the water table are assumed to be small relative to the saturated thickness, and vertical flow into or out of the zone above the water table is assumed to occur instantaneously. Effects of hysteresis in the moisture distribution above the water table are therefore neglected. Graphical comparisons of the new solutions are made with known closed-form solutions.

Journal of Hydrology↗
USGS2000

Aquifer response to stream-stage and recharge variations. II. Convolution method and applications

In this second of two papers, analytical step-response functions, developed in the companion paper for several cases of transient hydraulic interaction between a fully penetrating stream and a confined, leaky, or water-table aquifer, are used in the convolution integral to calculate aquifer heads, streambank seepage rates, and bank storage that occur in response to streamstage fluctuations and basinwide recharge or evapotranspiration. Two computer programs developed on the basis of these step-response functions and the convolution integral are applied to the analysis of hydraulic interaction of two alluvial stream-aquifer systems in the northeastern and central United States. These applications demonstrate the utility of the analytical functions and computer programs for estimating aquifer and streambank hydraulic properties, recharge rates, streambank seepage rates, and bank storage. Analysis of the water-table aquifer adjacent to the Blackstone River in Massachusetts suggests that the very shallow depth of water table and associated thin unsaturated zone at the site cause the aquifer to behave like a confined aquifer (negligible specific yield). This finding is consistent with previous studies that have shown that the effective specific yield of an unconfined aquifer approaches zero when the capillary fringe, where sediment pores are saturated by tension, extends to land surface. Under this condition, the aquifer's response is determined by elastic storage only. Estimates of horizontal and vertical hydraulic conductivity, specific yield, specific storage, and recharge for a water-table aquifer adjacent to the Cedar River in eastern Iowa, determined by the use of analytical methods, are in close agreement with those estimated by use of a more complex, multilayer numerical model of the aquifer. Streambank leakance of the semipervious streambank materials also was estimated for the site. The streambank-leakance parameter may be considered to be a general (or lumped) parameter that accounts not only for the resistance of flow at the river-aquifer boundary, but also for the effects of partial penetration of the river and other near-stream flow phenomena not included in the theoretical development of the step-response functions.Analytical step-response functions, developed for several cases of transient hydraulic interaction between a fully penetrating stream and a confined, leaky, or water-table aquifer, are used in the convolution integral to calculate aquifer heads, streambank seepage rates, and bank storage that occur in response to stream-stage fluctuations and basinwide recharge or evapotranspiration. Two computer programs developed on the basis of these step-response functions and the convolution integral are applied to the analysis of hydraulic interaction of two alluvial stream-aquifer systems. These applications demonstrate the utility of the analytical functions and computer programs for estimating aquifer and streambank seepage rates and bank storage.

Journal of Hydrology↗
USGS2000

Flowmetering of drainage wells in Kuwait City, Kuwait

A heat-pulse flowmeter was used in six drainage wells in Kuwait City for flow profiling under both ambient and pumping conditions. The data collected were used in: (a) estimating the cross-flow among the screened intervals under ambient conditions; (b) estimating the relative transmissivity adjacent to the individual screen zones; and (c) determination of the hydraulic heads at the far boundaries of the large-scale aquifer zones. These inferences were cross-checked against known hydrogeology of the aquifer-aquitard system in the study area, and the calibration results of numerical flow modeling. The major conclusions derived from the flow measurements were: (a) the presence of natural downward cross-flow under ambient condition supported the hypothesis that the upper part of the Kuwait Group aquifer in the study area was divided into a series of permeable units (aquifers), separated by confining or semi-confining beds (aquitards); (b) the head differences between the different screened zones, derived through modeling of the flowmeter data of the wells, provided additional confirmation for the division of the upper part of the Kuwait Group aquifer into compartments in the study area; (c) flowmeter data indicated that the second and third aquifers were contributing most of the water to the well bores, compared with the uppermost (first) and the lowermost (fourth) aquifers; and (d) inflow to the wells during pumping was associated with discrete sub-intervals in the screened zones, controlled by local aquifer heterogeneity, and possibly clogging of screens and gravel pack.

Kuwait City↗
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↗
USGS2000

Concentrations and characteristics of organic carbon in surface water in Arizona: Influence of urbanization

Dissolved (DOC) and total (TOC) organic carbon concentrations and compositions were studied for several river systems in Arizona, USA. DOC composition was characterized by ultraviolet and visible absorption and fluorescence emission (excitation wavelength of 370 nm) spectra characteristics. Ephemeral sites had the highest DOC concentrations, and unregulated perennial sites had lower concentrations than unregulated intermittent sites, regulated sites, and sites downstream from wastewater-treatment plants (p < 0.05). Reservoir outflows and wastewater-treatment plant effluent were higher in DOC concentration (p < 0.05) and exhibited less variability in concentration than inflows to the reservoirs. Specific ultraviolet absorbance values at 254 nm were typically less than 2 m-1(milligram DOC per liter)-1 and lower than values found in most temperate-region rivers, but specific ultraviolet absorbance values increased during runoff events. Fluorescence measurements indicated that DOC in desert streams typically exhibit characteristics of autochthonous sources; however, DOC in unregulated upland rivers and desert streams experienced sudden shifts from autochthonous to allochthonous sources during runoff events. The urban water system (reservoir systems and wastewater-treatment plants) was found to affect temporal variability in DOC concentration and composition. (C) 2000 Elsevier Science B.V.Dissolved (DOC) and total (TOC) organic carbon concentrations and compositions were studied for several river systems in Arizona, USA. DOC composition was characterized by ultraviolet and visible absorption and fluorescence emission (excitation wavelength of 370 nm) spectra characteristics. Ephemeral sites had the highest DOC concentrations, and unregulated perennial sites had lower concentrations than unregulated intermittent sites, regulated sites, and sites downstream from wastewater-treatment plants (p<0.05). Reservoir outflows and wastewater-treatment plant effluent were higher in DOC concentration (p<0.05) and exhibited less variability in concentration than inflows to the reservoirs. Specific ultraviolet absorbance values at 254 nm were typically less than 2 m-1(milligram DOC per liter)-1 and lower than values found in most temperate-region rivers, but specific ultraviolet absorbance values increased during runoff events. Fluorescence measurements indicated that DOC in desert streams typically exhibit characteristics of autochthonous sources; however, DOC in unregulated upland rivers and desert streams experienced sudden shifts from autochthonous to allochthonous sources during runoff events. The urban water system (reservoir systems and wastewater-treatment plants) was found to affect temporal variability in DOC concentration and composition.The influence of urbanization, becoming increasingly common in arid regions, on dissolved organic carbon (DOC) concentrations in surface water resources was studied. DOC concentration and composition, seasonal watershed runoff events, streamflow variations, water management practices, and urban infrastructure in several Arizona watersheds were monitored. Ephemeral sites had the highest DOC levels, and unregulated perennial sites and lower concentrations than unregulated intermittent sites, regulated sites, and sites downstream from wastewater treatment plants. Reservoir outflows and wastewater treatment plant effluent had higher and less variable DOC concentrations than inflows to reservoirs. UV absorbance values, fluorescence measurements, and other indicators suggest that urban water systems (reservoirs and wastewater treatment plants) affect temporal variability in DOC concentration and composition.

Journal of Hydrology↗
USGS2000

Geochemical modeling of iron, sulfur, oxygen and carbon in a coastal plain aquifer

Fe(III) reduction in the Magothy aquifer of Long Island, NY, results in high dissolved-iron concentrations that degrade water quality. Geochemical modeling was used to constrain iron-related geochemical processes and redox zonation along a flow path. The observed increase in dissolved inorganic carbon is consistent with the oxidation of sedimentary organic matter coupled to the reduction of O2 and SO4/2- in the aerobic zone, and to the reduction of SO4/2- in the anaerobic zone; estimated rates of CO2 production through reduction of Fe(III) were relatively minor by comparison. The rates of CO2 production calculated from dissolved inorganic carbon mass transfer (2.55 x 10-4 to 48.6 x 10-4 mmol 1-1 yr-1) generally were comparable to the calculated rates of CO2 production by the combined reduction of O2, Fe(III) and SO4/2- (1.31 x 10-4 to 15 x 10-4 mmol 1-1 yr-1). The overall increase in SO4/2- concentrations along the flow path, together with the results of mass-balance calculations, and variations in ??34S values along the flow path indicate that SO4/2- loss through microbial reduction is exceeded by SO4/2- gain through diffusion from sediments and through the oxidation of FeS2. Geochemichal and microbial data on cores indicate that Fe(III) oxyhydroxide coatings on sediment grains in local, organic carbon- and SO4/2- -rich zones have localized SO4/2- -reducing zones in which the formation of iron disulfides been depleted by microbial reduction and resulted in decreases dissolved iron concentrations. These localized zones of SO4/2- reduction, which are important for assessing zones of low dissolved iron for water-supply development, could be overlooked by aquifer studies that rely only on groundwater data from well-water samples for geochemical modeling. (C) 2000 Elsevier Science B.V.Fe(III) reduction in the Magothy aquifer of Long Island, NY, results in high dissolved-iron concentrations that degrade water quality. Geochemical modeling was used to constrain iron-related geochemical processes and redox zonation along a flow path. The observed increase in dissolved inorganic carbon is consistent with the oxidation of sedimentary organic matter coupled to the reduction of O2 and SO42- in the aerobic zone, and to the reduction of SO42- in the anaerobic zone; estimated rates of CO2 production through reduction of Fe(III) were relatively minor by comparison. The rates of CO2 production calculated from dissolved inorganic carbon mass transfer (2.55??10-4 to 48.6??10-4mmol l-1yr-1) generally were comparable to the calculated rates of CO2 production by the combined reduction of O2, Fe(III) and SO42- (1.31??10-4 to 15??10-4mmol l-1yr-1). The overall increase in SO42- concentrations along the flow path, together with the results of mass-balance calculations, and variations in ??34S values along the flow path indicate that SO42- loss through microbial reduction is exceeded by SO42- gain through diffusion from sediments and through the oxidation of FeS2. Geochemical and microbial data on cores indicate that Fe(III) oxyhydroxide coatings on sediment grains in local, organic carbon- and SO42--rich zones have been depleted by microbial reduction and resulted in localized SO42--reducing zones in which the formation of iron disulfides decreases dissolved iron concentrations. These localized zones of SO42- reduction, which are important for assessing zones of low dissolved iron for water-supply development, could be overlooked by aquifer studies that rely only on groundwater data from well-water samples for geochemical modeling.

Journal of Hydrology↗
USGS2000

Modeling regional salinization of the Ogallala aquifer, Southern High Plains, TX, USA

Two extensive plumes (combined area > 1000 km2) have been delineated within the Ogallala aquifer in the Southern High Plains, TX, USA. Salinity varies within the plumes spatially and increases with depth; Cl ranges from 50 to >500 mg 1-1. Variable-density flow modeling using SUTRA has identified three broad regions of upward cross-formational flow from the underlying evaporite units. The upward discharge within the modeled plume area is in the range of 10-4-10-5 m3 day-1, and the TDS concentrations are typically >3000 mg 1-1. Regions of increased salinity, identified within the Whitehorse Group (evaporite unit) underlying the Ogallala aquifer, are controlled by the structure and thickness variations relative to the recharge areas. Distinct flow paths, on the order of tens of km to >100 km in length, and varying flow velocities indicate that the salinization of the Ogallala aquifer has been a slow, ongoing process and may represent circulation of waters recharged during Pleistocene or earlier times. On-going pumping has had negligible impact on the salinity distribution in the Ogallala aquifer, although simulations indicate that the velocity distribution in the underlying units may have been affected to depths of 150 m after 30 years of pumping. Because the distribution of saline ground water in this region of the Ogallala aquifer is heterogeneous, careful areal and vertical characterization is warranted prior to any well-field development. (C) 2000 Elsevier Science B.V.Two extensive plumes (combined area >1000 km2) have been delineated within the Ogallala aquifer in the Southern High Plains, TX, USA. Salinity varies within the plumes spatially and increases with depth; Cl ranges from 50 to >500 mg l-1. Variable-density flow modeling using SUTRA has identified three broad regions of upward cross-formational flow from the underlying evaporite units. The upward discharge within the modeled plume area is in the range of 10-4-10-5 m3 day-1, and the TDS concentrations are typically >3000 mg l-1. Regions of increased salinity, identified within the Whitehorse Group (evaporite unit) underlying the Ogallala aquifer, are controlled by the structure and thickness variations relative to the recharge areas. Distinct flow paths, on the order of tens of km to >100 km in length, and varying flow velocities indicate that the salinization of the Ogallala aquifer has been a slow, ongoing process and may represent circulation of waters recharged during Pleistocene or earlier times. On-going pumping has had negligible impact on the salinity distribution in the Ogallala aquifer, although simulations indicate that the velocity distribution in the underlying units may have been affected to depths of 150 m after 30 years of pumping. Because the distribution of saline ground water in this region of the Ogallala aquifer is heterogeneous, careful areal and vertical characterization is warranted prior to any well-field development.

Journal of Hydrology↗