USGS Science⌕ Search

SEARCH · USGS Science

Results for “groundwater aquifer”

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 685 records · Page 38Linked to original sources

A finite-element model for simulating hydraulic interchange of surface and ground water

A model was developed to be useful for predicting changes in streamflow as a result of groundwater pumping. The stream aquifer model is especially useful for simulating streams that flow intermittently owing to leakage to the aquifer or diversion for irrigation or streams that become perched owing to declining hydraulic head in the aquifer. The model couples the equation of two-dimensional groundwater flow with the kinematic equations of one-dimensional open-channel flow. Darcy 's law for vertical flow through a semipermeable streambed is used to couple the groundwater flow and streamflow equations. The equations of flow are approximated numerically by the finite-element method. A listing of the Fortran program that solves the equations of flow , and a description of data-input formats are given in the report. The program can simulate a variety of hydrologic characteristics including perched streams, streamflow diversions , springs, recharge from irrigated acreage, and evapotranspiration from the water table and phreatophytes. Time-dependent boundary conditions can be simulated. The program can be modified easily to simulate unconfined aquifers and aquifers with variable directions of anisotropy. (USGS)

Water-Resources Investigations Report↗

Potential effects of roadside dry wells on groundwater quality on the Island of Hawai'i — Assessment using numerical groundwater models

Widespread use of dry wells to dispose of roadside runoff has raised concern about the potential effects on the quality of groundwater on the Island of Hawai‘i. This study used semi-generic numerical models of groundwater flow and contaminant transport to assess the potential effect of dry wells on groundwater quality on the Island of Hawai‘i. The semi-generic models are generalized numerical groundwater-flow and solute-transport models that have a range of aquifer properties and regional groundwater gradients that are characteristic for the island. Several semi-generic models were created to study the effect of dry wells in different hydrogeologic conditions, such as different unsaturated-zone thicknesses or different aquifer characteristics. Results indicate that mixing of contaminated water from the surface with contaminant-free water in the saturated aquifer immediately reduces the contaminant concentration. The amount the concentration is reduced depends on the hydraulic properties of the aquifer in a given area, the thickness of the unsaturated zone, and whether the infiltration is focused in a small area of a dry well or spread naturally over a larger area. Model simulations indicate that focusing infiltration of contaminated runoff through a dry well can substantially increase contaminant concentrations in the underlying saturated aquifer relative to infiltration under natural conditions. Simulated concentrations directly beneath a dry well were nearly 8 times higher than the simulated concentrations directly beneath a broad infiltration area representing the natural condition. Where dry wells are present, contaminant concentrations in the underlying saturated aquifer are lower when the unsaturated zone is thicker and higher when the unsaturated zone is thinner. Contaminant concentrations decline quickly as the contaminant plume migrates, with the regional groundwater flow, away from the dry well. The differences among concentrations resulting from the various unsaturated-zone thicknesses also diminish with distance from the dry well. At a horizontal distance of about 700 ft downgradient from the dry well, all simulated maximum concentrations were less than 1 percent of the concentration in the infiltration water; at about 0.5 mi downgradient from the dry well, all simulated concentrations were equal to or less than 0.1 percent. Actual concentrations may be even lower than indicated by the models because of processes such as decay and reaction that were not simulated. Hydrologic and geologic differences from one location to the next also affect contaminant concentrations—simulations using models with properties representative of aquifers in the Hilo area resulted in lower overall concentrations than models with properties representative of aquifers in the Kona area. Results from this study can be used to assess how contaminants entering a dry well may affect receiving waters in a variety of situations on the Island of Hawai‘i. Better assessment would be obtained by using results from models having the most similar conditions (such as climate, hydraulic properties, regional groundwater gradient) to the dry well in question. The results of this study can help determine which dry wells are likely to have the greatest effect on nearby receiving waters and where more specific data and analyses may be needed.

Hawaii↗

Examples of deformation-dependent flow simulations of conjunctive use with MF-OWHM

The dependency of surface- and groundwater flows and aquifer hydraulic properties on deformation induced by changes in aquifer head is not accounted for in the standard version of MODFLOW. A new USGS integrated hydrologic model, MODFLOW-OWHM, incorporates this dependency by linking subsidence and mesh deformation with changes in aquifer transmissivity and storage coefficient, and with flows that also depend on aquifer characteristics and land-surface geometry. This new deformation-dependent approach is being used for the further development of the integrated Central Valley hydrologic model (CVHM) in California. Preliminary results from this application and from hypothetical test cases of similar systems show that changes in canal flows, stream seepage, and evapotranspiration from groundwater (ET gw ) are sensitive to deformation. Deformation feedback has been shown to also have an indirect effect on conjunctive surface- and groundwater use components with increased stream seepage and streamflows influencing surface-water deliveries and return flows. In the Central Valley model, land subsidence may significantly degrade the ability of the major canals to deliver surface water from the Delta to the San Joaquin and Tulare basins. Subsidence can also affect irrigation demand and ET gw , which, along with altered surface-water supplies, causes a feedback response resulting in changed estimates of groundwater pumping for irrigation. This modeling feature also may improve the impact assessment of dewatering-induced land subsidence/uplift (following irrigation pumping or coal-seam gas extraction) on surface receptors, inter-basin transfers, and surface infrastructure integrity.

Proceedings of the International Association of Hy↗

Groundwater-flow model of the northern High Plains aquifer in Colorado, Kansas, Nebraska, South Dakota, and Wyoming

The High Plains aquifer is a nationally important water resource underlying about 175,000 square miles in parts of eight states: Colorado, Kansas, Oklahoma, Nebraska, New Mexico, South Dakota, Texas, and Wyoming. Droughts across much of the Northern High Plains from 2001 to 2007 have combined with recent (2004) legislative mandates to elevate concerns regarding future availability of groundwater and the need for additional information to support science-based water-resource management. To address these needs, the U.S. Geological Survey began the High Plains Groundwater Availability Study to provide a tool for water-resource managers and other stakeholders to assess the status and availability of groundwater resources. A transient groundwater-flow model was constructed using the U.S. Geological Survey modular three-dimensional finite-difference groundwater-flow model with Newton-Rhapson solver (MODFLOW–NWT). The model uses an orthogonal grid of 565 rows and 795 columns, and each grid cell measures 3,281 feet per side, with one variably thick vertical layer, simulated as unconfined. Groundwater flow was simulated for two distinct periods: (1) the period before substantial groundwater withdrawals, or before about 1940, and (2) the period of increasing groundwater withdrawals from May 1940 through April 2009. A soil-water-balance model was used to estimate recharge from precipitation and groundwater withdrawals for irrigation. The soil-water-balance model uses spatially distributed soil and landscape properties with daily weather data and estimated historical land-cover maps to calculate spatial and temporal variations in potential recharge. Mean annual recharge estimated for 1940–49, early in the history of groundwater development, and 2000–2009, late in the history of groundwater development, was 3.3 and 3.5 inches per year, respectively. Primary model calibration was completed using statistical techniques through parameter estimation using the parameter estimation suite of software with Tikhonov regularization. Calibration targets for the groundwater model included 343,067 groundwater levels measured in wells and 10,820 estimated monthly stream base flows at streamgages. A total of 1,312 parameters were adjusted during calibration to improve the match between calibration targets and simulated equivalents. Comparison of calibration targets to simulated equivalents indicated that, at the regional scale, the model correctly reproduced groundwater levels and stream base flows for 1940–2009. This comparison indicates that the model can be used to examine the likely response of the aquifer system to potential future stresses. Mean calibrated recharge for 1940–49 and 2000–2009 was smaller than that estimated with the soil-water-balance model. This indicated that although the general spatial patterns of recharge estimated with the soil-water-balance model were approximately correct at the regional scale of the Northern High Plains aquifer, the soil-water-balance model had overestimated recharge, and adjustments were needed to decrease recharge to improve the match of the groundwater model to calibration targets. The largest components of the simulated groundwater budgets were recharge from precipitation, recharge from canal seepage, outflows to evapotranspiration, and outflows to stream base flow. Simulated outflows to irrigation wells increased from 7 percent of total outflows in 1940–49 to 38 percent of 1970–79 total outflows and 49 percent of 2000–2009 total outflows.

Colorado, Kansas, Nebraska, South Dakota, Wyoming↗

Effects of groundwater withdrawal on borehole flow and salinity measured in deep monitor wells in Hawai'i: implications for groundwater management

Water-resource managers in Hawai`i rely heavily on salinity profiles from deep monitor wells to estimate the thickness of freshwater and the depth to the midpoint of the transition zone between freshwater and saltwater in freshwater-lens systems. The deep monitor wells are typically open boreholes below the water table and extend hundreds of feet below sea level. Because of possible borehole-flow effects, there is concern that salinity profiles measured in these wells may not accurately reflect the salinity distribution in the aquifer and consequently lead to misinterpretations that adversely affect water-resource management. Steplike changes in salinity or temperature with depth in measured profiles from nonpumped deep monitor wells may be indicative of water moving within the well, and such changes are evident to some extent in all available profiles. The maximum vertical step length, or displacement, in measured profiles ranges from 7 to 644 feet. Vertical steps longer than 70 feet exceed the typical thickness of massive lava flows; they therefore cannot be attributed entirely to geologic structure and may be indicative of borehole flow. The longest vertical steps occur in monitor wells located in southern O'ahu, coinciding with the most heavily developed part of the aquifer. Although regional groundwater withdrawals have caused a thinning of the freshwater lens over the past several decades, the measured midpoint of the transition zone in most deep monitor wells has shown only inconsequential depth displacement in direct response to short-term variations in withdrawals from nearby production wells. For profiles from some deep monitor wells, however, the depth of the measured top of the transition zone, indicated by a specific-conductance value of 1,000 microsiemens per centimeter, has risen several hundred feet in response to withdrawals from nearby production wells. For these deep monitor wells, monitoring the apparent top of the transition zone may not provide an accurate indication of water quality in the adjacent aquifer. Hence, the measured midpoint in boreholes is a better proxy for freshwater-lens thickness. Brackish water transported upward in a deep monitor well can exit the borehole in the upper, freshwater part of the aquifer and affect the water quality in nearby production wells. Piezometers installed at different depths will provide the best information on aquifer salinity because they are unaffected by borehole flow. Despite the effects of borehole flow, monitoring the midpoint in deep monitor wells is still useful to identify long-term trends in the movement of the transition zone.

Hawaii↗

Enhancing drought resilience with conjunctive use and managed aquifer recharge in California and Arizona

Projected longer‐term droughts and intense floods underscore the need to store more water to manage climate extremes. Here we show how depleted aquifers have been used to store water by substituting surface water use for groundwater pumpage (conjunctive use, CU) or recharging groundwater with surface water (Managed Aquifer Recharge, MAR). Unique multi‐decadal monitoring from thousands of wells and regional modeling datasets for the California Central Valley and central Arizona were used to assess CU and MAR. In addition to natural reservoir capacity related to deep water tables, historical groundwater depletion further expanded aquifer storage by ~44 km3 in the Central Valley and by ~100 km3 in Arizona, similar to or exceeding current surface reservoir capacity by up to three times. Local river water and imported surface water, transported through 100s of km of canals, is substituted for groundwater (≤15 km3/yr, CU) or is used to recharge groundwater (MAR, ≤1.5 km3/yr) during wet years shifting to mostly groundwater pumpage during droughts. In the Central Valley, CU and MAR locally reversed historically declining water‐level trends, which contrasts with simulated net regional groundwater depletion. In Arizona, CU and MAR also reversed historically declining groundwater level trends in Active Management Areas. These rising trends contrast with current declining trends in irrigated areas that lack access to surface water to support CU or MAR. Use of depleted aquifers as reservoirs could expand with winter flood irrigation or capturing flood discharges to the Pacific (0 – 1.6 km3/yr, 2000–2014) with additional infrastructure in California. Because flexibility and expanded portfolio options translate to resilience, CU and MAR enhance drought resilience through multi‐year storage, complementing shorter term surface reservoir storage, and facilitating water markets.

Environmental Research Letters↗

Effects of groundwater levels and headwater wetlands on streamflow in the Charlie Creek basin, Peace River watershed, west-central Florida

The Charlie Creek basin was studied from April 2004 to December 2005 to better understand how groundwater levels in the underlying aquifers and storage and overflow of water from headwater wetlands preserve the streamflows exiting this least-developed tributary basin of the Peace River watershed. The hydrogeologic framework, physical characteristics, and streamflow were described and quantified for five subbasins of the 330-square mile Charlie Creek basin, allowing the contribution of its headwaters area and tributary subbasins to be separately quantified. A MIKE SHE model simulation of the integrated surface-water and groundwater flow processes in the basin was used to simulate daily streamflow observed over 21 months in 2004 and 2005 at five streamflow stations, and to quantify the monthly and annual water budgets for the five subbasins including the changing amount of water stored in wetlands. Groundwater heads were mapped in Zone 2 of the intermediate aquifer system and in the Upper Floridan aquifer, and were used to interpret the location of artesian head conditions in the Charlie Creek basin and its relation to streamflow. Artesian conditions in the intermediate aquifer system induce upward groundwater flow into the surficial aquifer and help sustain base flow which supplies about two-thirds of the streamflow from the Charlie Creek basin. Seepage measurements confirmed seepage inflow to Charlie Creek during the study period. The upper half of the basin, comprised largely of the Upper Charlie Creek subbasin, has lower runoff potential than the lower basin, more storage of runoff in wetlands, and periodically generates no streamflow. Artesian head conditions in the intermediate aquifer system were widespread in the upper half of the Charlie Creek basin, preventing downward leakage from expansive areas of wetlands and enabling them to act as headwaters to Charlie Creek once their storage requirements were met. Currently, the dynamic balance between wetland storage, rainfall-runoff processes, and groundwater-level differences in the upper basin allow it to generate approximately half of the streamflow from the Charlie Creek basin. Therefore, future development in the upper basin that would alter the hydraulic connectivity of wetlands during high flow conditions or expand recharging groundwater conditions could substantially affect streamflow in Charlie Creek. LIDAR (Light detection and ranging) based topographic maps and integrated modeling results were used to quantify the water stored in wetlands and other topographic depressions, and to describe the network of shallow stream channels connecting wetlands to Charlie Creek and its tributaries over distances of several thousand feet. Peak flows at all but one streamflow station were underpredicted in MIKE SHE simulations, possibly because the hydraulics of surface channels connecting wetlands to stream channels were not explicitly simulated in the model. Explicitly simulating the smaller channels connecting wetlands and stream channels should improve the ability of future watershed models to simulate peak flows in streams with headwater wetlands. The runoff potential was greater in the lower half of the Charlie Creek basin than in the upper half, and the streambed of Charlie Creek had greater potential to both directly gain streamflow from groundwater and lose streamflow to groundwater. Charlie Creek is more incised into the surficial aquifer in the lower basin than in the upper basin, and the streambed intersects the top of the intermediate aquifer system at two known locations. Groundwater levels in the intermediate aquifer system varied widely in the lower half of the basin from artesian conditions inducing upward flow toward the surficial aquifer and streams, to recharging conditions allowing downward flow and stream leakage. Recharge areas were greatest in May 2004 when rainfall was at a seasonal low and irrigation pumping was at a seasonal high. Recharge conditions

Scientific Investigations Report↗

Geogenic organic contaminants in the low-rank coal-bearing Carrizo-Wilcox aquifer of East Texas, USA

The organic composition of groundwater along the Carrizo-Wilcox aquifer in East Texas (USA), sampled from rural wells in May and September 2015, was examined as part of a larger study of the potential health and environmental effects of organic compounds derived from low-rank coals. The quality of water from the low-rank coal-bearing Carrizo-Wilcox aquifer is a potential environmental concern and no detailed studies of the organic compounds in this aquifer have been published. Organic compounds identified in the water samples included: aliphatics and their fatty acid derivatives, phenols, biphenyls, N-, O-, and S-containing heterocyclic compounds, polycyclic aromatic hydrocarbons (PAHs), aromatic amines, and phthalates. Many of the identified organic compounds (aliphatics, phenols, heterocyclic compounds, PAHs) are geogenic and originated from groundwater leaching of young and unmetamorphosed low-rank coals. Estimated concentrations of individual compounds ranged from about 3.9 to 0.01 μg/L. In many rural areas in East Texas, coal strata provide aquifers for drinking water wells. Organic compounds observed in groundwater are likely to be present in drinking water supplied from wells that penetrate the coal. Some of the organic compounds identified in the water samples are potentially toxic to humans, but at the estimated levels in these samples, the compounds are unlikely to cause acute health problems. The human health effects of low-level chronic exposure to coal-derived organic compounds in drinking water in East Texas are currently unknown, and continuing studies will evaluate possible toxicity.

Texas↗

Contribution of global groundwater depletion since 1900 to sea-level rise

Removal of water from terrestrial subsurface storage is a natural consequence of groundwater withdrawals, but global depletion is not well characterized. Cumulative groundwater depletion represents a transfer of mass from land to the oceans that contributes to sea-level rise. Depletion is directly calculated using calibrated groundwater models, analytical approaches, or volumetric budget analyses for multiple aquifer systems. Estimated global groundwater depletion during 1900–2008 totals ~4,500 km 3 , equivalent to a sea-level rise of 12.6 mm (>6% of the total). Furthermore, the rate of groundwater depletion has increased markedly since about 1950, with maximum rates occurring during the most recent period (2000–2008), when it averaged ~145 km 3 /yr (equivalent to 0.40 mm/yr of sea-level rise, or 13% of the reported rate of 3.1 mm/yr during this recent period).

Geophysical Research Letters↗

Fluid inclusions and biomarkers in the Upper Mississippi Valley zinc-lead district; implications for the fluid-flow and thermal history of the Illinois Basin

The Upper Mississippi Valley zinc-lead district is hosted by Ordovician carbonate rocks at the northern margin of the Illinois Basin. Fluid inclusion temperature measurements on Early Permian sphalerite ore from the district are predominantly between 90?C and I50?C. These temperatures are greater than can be explained by their reconstructed burial depth, which was a maximum of approximately 1 km at the time of mineralization. In contrast to the temperatures of mineral formation derived from fluid inclusions, biomarker maturities in the Upper Mississippi Valley district give an estimate of total thermal exposure integrated over time. Temperatures from fluid inclusions trapped during ore genesis with biomarker maturities were combined to construct an estimate of the district's overall thermal history and, by inference, the late Paleozoic thermal and hydrologic history of the Illinois Basin. Circulation of groundwater through regional aquifers, given sufficient flow rates, can redistribute heat from deep in a sedimentary basin to its shallower margins. Evidence for regional-scale circulation of fluids is provided by paleomagnetic studies, regionally correlated zoned dolomite, fluid inclusions, and thermal maturity of organic matter. Evidence for igneous acti vity contemporaneous with mineralization in the vicinity of the Upper Mississippi Valley district is absent. Regional fluid and heat circulation is the most likely explanation for the elevated fluid inclusion temperatures (relative to maximum estimated burial depth) in the Upper Mississippi Valley district. One plausible driving mechanism and flow path for the ore-forming fluids is groundwater recharge in the late Paleozoic Appalachian-Ouachita mountain belt and northward flow through the Reelfoot rift and the proto- Illinois Basin to the Upper Mississippi Valley district. Warm fluid flowing laterally through Cambrian and Ordovician aquifers would then move vertically upward through the fractures that control sphalerite mineralization in the Upper Mississippi Valley district. Biomarker reactant-product measurements on rock extracts from the Upper Mississippi Valley district define a relatively low level ofthermal maturity for the district, 0.353 for sterane and 0.577 for hopane. Recently published kinetic constants permit a time-temperature relationship to be determined from these biomarker maturities. Numerical calculations were made to simulate fluid heat flow through the fracture-controlled ore zones of the Thompson-Temperly mine and heat transfer to the adjacent rocks where biomarker samples were collected. Calculations that combine the fluid inclusion temperatures and the biomarker constraints on thermal maturity indicate that the time interval during which mineralizing fluids circulated through the Upper Mississippi Valley district is on the order of 200,000 years. Fluid inclusion measurements and thermal maturities from biomarkers in the district reflect the duration of peak temperatures resulting from regional fluid circulation. On the basis of thermal considerations, the timing of fluorite mineralization in southern Illinois, and the northward-decreasing pattern of fluorine enrichment in sediments, we hypothesize that the principal flow direction was northward through the Cambrian and Ordovician aquifers of the Illinois Basin. A basin-scale flow system would result in mass transport (hydrocarbon migration, transport of metals in solution) and energy (heat) transport, which would in turn drive chemical reactions (for example, maturation of organic matter, mineralization, diagenetic reactions) within the Illinois Basin and at its margins.

Bulletin↗

Stream seepage and groundwater levels, Wood River Valley, south-central Idaho, 2012-13

Stream discharge and water levels in wells were measured at multiple sites in the Wood River Valley, south-central Idaho, in August 2012, October 2012, and March 2013, as a component of data collection for a groundwater-flow model of the Wood River Valley aquifer system. This model is a cooperative and collaborative effort between the U.S. Geological Survey and the Idaho Department of Water Resources. Stream-discharge measurements for determination of seepage were made during several days on three occasions: August 27–28, 2012, October 22–24, 2012, and March 27–28, 2013. Discharge measurements were made at 49 sites in August and October, and 51 sites in March, on the Big Wood River, Silver Creek, their tributaries, and nearby canals. The Big Wood River generally gains flow between the Big Wood River near Ketchum streamgage (13135500) and the Big Wood River at Hailey streamgage (13139510), and loses flow between the Hailey streamgage and the Big Wood River at Stanton Crossing near Bellevue streamgage (13140800). Shorter reaches within these segments may differ in the direction or magnitude of seepage or may be indeterminate because of measurement uncertainty. Additional reaches were measured on Silver Creek, the North Fork Big Wood River, Warm Springs Creek, Trail Creek, and the East Fork Big Wood River. Discharge measurements also were made on the Hiawatha, Cove, District 45, Glendale, and Bypass Canals, and smaller tributaries to the Big Wood River and Silver Creek. Water levels in 93 wells completed in the Wood River Valley aquifer system were measured during October 22–24, 2012; these wells are part of a network established by the U.S. Geological Survey in 2006. Maps of the October 2012 water-table altitude in the unconfined aquifer and the potentiometric-surface altitude of the confined aquifer have similar topology to those on maps of October 2006 conditions. Between October 2006 and October 2012, water-table altitude in the unconfined aquifer rose by as much as 1.86 feet in 6 wells and declined by as much as 14.28 feet in 77 wells; average decline was 2.9 feet. A map of changes in the water‑table altitude of the unconfined aquifer shows that the largest declines were in tributary canyons and in an area roughly between Baseline and Glendale Roads. From October 2006 to October 2012, the potentiometric-surface altitude in 10 wells completed in the confined aquifer declined between 0.12 and 20.50 feet; average decline was 6.8 feet. A map of changes in the potentiometric-surface altitude of the confined aquifer shows that the largest declines were in the southwestern part of the Bellevue fan. Reduced precipitation prior to the October 2012 water-level measurements likely is partially responsible for 2006–12 water-table declines in the unconfined aquifer; the relative contribution of precipitation deficit and groundwater withdrawals to the declines is not known. Although the confined aquifer may not receive direct recharge from precipitation or streams, groundwater withdrawal from the confined aquifer induces flow from the unconfined aquifer. Declines in the confined aquifer are likely due to groundwater withdrawals and declines in the water table of the unconfined aquifer. A statistical analysis of five long-term monitoring wells (three completed in the unconfined aquifer, one in the confined aquifer, and one outside the aquifer system boundary) showed statistically significant declining trends in four wells.

Idaho↗

Evaluation of chemical and hydrologic processes in the eastern Snake River Plain Aquifer based on results from geochemical modeling, Idaho National Laboratory, eastern Idaho

Nuclear research activities at the U.S. Department of Energy (DOE) Idaho National Laboratory (INL) produced liquid and solid chemical and radiochemical wastes that were disposed to the subsurface resulting in detectable concentrations of some waste constituents in the eastern Snake River Plain (ESRP) aquifer. These waste constituents may affect the water quality of the aquifer and may pose risks to the eventual users of the aquifer water. To understand these risks to water quality the U.S. Geological Survey, in cooperation with the DOE, conducted geochemical mass-balance modeling of the ESRP aquifer to improve the understanding of chemical reactions, sources of recharge, mixing of water, and groundwater flow directions in the shallow (upper 250 feet) aquifer at the INL. Modeling was conducted using the water chemistry of 127 water samples collected from sites at and near the INL. Water samples were collected between 1952 and 2017 with most of the samples collected during the mid-1990s. Geochemistry and isotopic data used in geochemical modeling consisted of dissolved oxygen, carbon dioxide, major ions, silica, aluminum, iron, and the stable isotope ratios of hydrogen, oxygen, and carbon. Geochemical modeling results indicated that the primary chemical reactions in the aquifer were precipitation of calcite and dissolution of plagioclase (An 60 ) and basalt volcanic glass. Secondary minerals other than calcite included calcium montmorillonite and goethite. Reverse cation exchange, consisting of sodium exchanging for calcium on clay minerals, occurred near site facilities where large amounts of sodium were released to the ESRP aquifer in wastewater discharge. Reverse cation exchange acted to retard the movement of wastewater-derived sodium in the aquifer. Regional groundwater inflow was the primary source of recharge to the aquifer underlying the Northeast and Southeast INL Areas. Birch Creek (BC), the Big Lost River (BLR), and groundwater from BC valley provided recharge to the North INL Area, and the BLR and groundwater from BC and Little Lost River (LLR) valleys provided recharge to the Central INL Area. The BLR, groundwater from the BLR and LLR valleys and the Lost River Range, and precipitation provided recharge to the Northwest and Southwest INL Areas. The primary source of recharge west and southwest of the INL was groundwater inflow from BLR valley. Upwelling geothermal water was a small source of recharge at two wells. Aquifer recharge from surface water in the northern, central, and western parts of the INL indicated that the aquifer in these areas was a dynamic, open system, whereas the aquifer in the eastern part of the INL, which receives little recharge from surface water, was a relatively static and closed system. Sources of recharge identified from isotope ratios and geochemical modeling (major ion concentrations) were nearly identical for the North, Northeast, Southeast, and Central INL Areas, which indicated that both methods probably accurately identified the sources of recharge in these areas. Conversely, isotope ratios indicated that the BLR and groundwater from the LLR valley provided most recharge to the western parts of the Northwest and Southwest INL Areas, whereas geochemical modeling results indicated a smaller area of recharge from the BLR and groundwater from the LLR valley, a larger area of recharge from the Lost River Range, and recharge of groundwater from the BLR valley that extended to the west INL boundary. The results from geochemical modeling probably were more accurate because major ion concentrations, but not isotope ratios, were available to characterize groundwater from the BLR valley and the Lost River Range. Sources of recharge identified with a groundwater flow model (using particle tracking) and geochemical modeling were similar for the Northeast and Southeast INL Areas. However, differences between the models were that the geochemical model represented (1) recharge of groundwater from the Lost River Range in the western part of the INL, whereas the flow model did not, (2) recharge of groundwater from the BC and BLR valleys extending farther south and east, respectively, than the flow model, and (3) more recharge from the BLR in the Southwest INL Area than the flow model. Mixing of aquifer water beneath the INL included (1) mixing of regional groundwater and water from the BC valley in the Northeast and Southeast INL Areas and (2) mixing of surface water (primarily from the BLR) and groundwater across much of the North, Central, Northwest, and Southwest INL Areas. Localized recharge from precipitation mixed with groundwater in the Northwest and Southwest INL Areas, and localized upwelling geothermal water mixed with groundwater in the Central and Northeast INL Areas. Flow directions of regional groundwater were south in the eastern part of the INL and south-southwest at downgradient locations. Groundwater from the BC and LLR valleys initially flowed southeast before changing to south-southwest flow directions that paralleled regional groundwater, and groundwater from the BLR valley initially flowed south before changing to a southsouthwest direction. Wastewater-contaminated groundwater flowed south from the Idaho Nuclear Technology and Engineering Center (INTEC) infiltration ponds in a narrow plume, with the percentage of wastewater in groundwater decreasing due to dilution, dispersion, and (or) degradation from about 60‒80 percent wastewater 0.7‒0.8 mile (mi) south of the INTEC infiltration ponds to about 1.4 percent wastewater about 15.5 mi south of the INTEC infiltration ponds. Wastewater contaminated groundwater flowed southeast and then southwest from the Naval Reactors Facility industrial waste ditch, with the percentage of wastewater in groundwater decreasing from about 100 percent wastewater adjacent to the waste ditch to about 2 percent wastewater about 0.6 mi south of the waste ditch.

Idaho↗

Ground-water use in the coastal plain of Maryland, 1900-1980

This report presents groundwater withdrawal data from 1900 through 1980 for Maryland counties lying with the Coastal Plain physiographic province, as well as a summary section for the total Maryland Coastal Plain. The types of water use included are domestic, military, water supplier, industrial/commercial, and irrigation. The data were obtained from state and county reports, biannual pumpage reports submitted to the Maryland Water Resources Administration, communication with individual owners, and estimates based on existing published data. The amount of groundwater withdrawn from aquifers in the Maryland Coastal Plain in 1900 was approximately 26 million gallons per day (Mgal/d) compared to nearly 134 Mgal/d in 1980. Jurisdictions withdrawing more than 10 Mgal/d for most of the 80-year period were Anne Arundel and Baltimore Counties and Baltimore City. The greatest withdrawals for most of the early part of the period were for domestic and industrial/commercial uses; however, water-supplier use dominated after 1965. Groundwater use for irrigation became important in the Coastal Plain around 1960 and increased steadily from approximately 2 Mgal/d in 1960 to nearly 12 Mgal/d in 1980. (USGS)

Open-File Report↗

Nitrogen biogeochemistry of submarine groundwater discharge

To investigate the role of the seepage zone in transport, chemical speciation, and attenuation of nitrogen loads carried by submarine groundwater discharge, we collected nearshore groundwater samples (n = 328) and examined the distribution and isotopic signature (δ 15 N) of nitrate and ammonium. In addition, we estimated nutrient fluxes from terrestrial and marine groundwater sources. We discuss our results in the context of three aquifer zones: a fresh groundwater zone, a shallow salinity transition zone (STZ), and a deep STZ. Groundwater plumes containing nitrate and ammonium occurred in the freshwater zone, whereas the deep STZ carried almost exclusively ammonium. The distributions of redox-cycled elements were consistent with theoretical thermodynamic stability of chemical species, with sharp interfaces between water masses of distinct oxidation : reduction potential, suggesting that microbial transformations of nitrogen were rapid relative to dispersive mixing. In limited locations in which overlap occurs between distribution of nitrate with that of ammonium and dissolved Fe 2+ , changes in concentration and in δ 15 N suggest loss of all species. Concurrent removal of NO 3 − and NH 4 + , both in freshwater and the deep STZ, might occur through a range of mechanisms, including heterotrophic or autotrophic denitrification, coupled nitrfication : denitrification, anammox, or Mn oxidation of NH 4 + . Loss of nitrogen was not apparent in the shallow STZ, perhaps because of short water residence time. Despite organic Cpoor conditions, the nearshore aquifer and subterranean estuary are biogeochemically active zones, where attenuation of N loads can occur. Extent of attenuation is controlled by the degree of mixing of biogeochemically dissimilar water masses, highlighting the critical role of hydrogeology in N biogeochemistry. Mixing is related in part to thinning of the freshwater lens before discharge and to dispersion at the fresh : saline groundwater interface, features common to all submarine groundwater discharge zones.

Limnology and Oceanography↗

Development and application of a comprehensive simulation model to evaluate impacts of watershed structures and irrigation water use on streamflow and groundwater: The case of Wet Walnut Creek Watershed, Kansas, USA

This paper presents the results of a comprehensive modeling study of surface and groundwater systems, including stream-aquifer interactions, for the Wet Walnut Creek Watershed in west-central Kansas. The main objective of this study was to assess the impacts of watershed structures and irrigation water use on streamflow and groundwater levels, which in turn affect availability of water for the Cheyenne Bottoms Wildlife Refuge Management area. The surface-water flow model, POTYLDR, and the groundwater flow model, MODFLOW, were combined into an integrated, watershed-scale, continuous simulation model. Major revisions and enhancements were made to the POTYLDR and MODFLOW models for simulating the detailed hydrologic budget for the Wet Walnut Creek Watershed. The computer simulation model was calibrated and verified using historical streamflow records (at Albert and Nekoma gaging stations), reported irrigation water use, observed water-level elevations in watershed structure pools, and groundwater levels in the alluvial aquifer system. To assess the impact of watershed structures and irrigation water use on streamflow and groundwater levels, a number of hypothetical management scenarios were simulated under various operational criteria for watershed structures and different annual limits on water use for irrigation. A standard 'base case' was defined to allow comparative analysis of the results of different scenarios. The simulated streamflows showed that watershed structures decrease both streamflows and groundwater levels in the watershed. The amount of water used for irrigation has a substantial effect on the total simulated streamflow and groundwater levels, indicating that irrigation is a major budget item for managing water resources in the watershed. (C) 2000 Elsevier Science B.V.This paper presents the results of a comprehensive modeling study of surface and groundwater systems, including stream-aquifer interactions, for the Wet Walnut Creek Watershed in west-central Kansas. The main objective of this study was to assess the impacts of watershed structures and irrigation water use on streamflow and groundwater levels, which in turn affect availability of water for the Cheyenne Bottoms Wildlife Refuge Management area. The surface-water flow model, POTYLDR, and the groundwater flow model, MODFLOW, were combined into an integrated, watershed-scale, continuous simulation model. Major revisions and enhancements were made to the POTYLDR and MODFLOW models for simulating the detailed hydrologic budget for the Wet Walnut Creek Watershed. The computer simulation model was calibrated and verified using historical streamflow records (at Albert and Nekoma gaging stations), reported irrigation water use, observed water-level elevations in watershed structure pools, and groundwater levels in the alluvial aquifer system. To assess the impact of watershed structures and irrigation water use on streamflow and groundwater levels, a number of hypothetical management scenarios were simulated under various operational criteria for watershed structures and different annual limits on water use for irrigation. A standard `base case' was defined to allow comparative analysis of the results of different scenarios. The simulated streamflows showed that watershed structures decrease both streamflows and groundwater levels in the watershed. The amount of water used for irrigation has a substantial effect on the total simulated streamflow and groundwater levels, indicating that irrigation is a major budget item for managing water resources in the watershed.A comprehensive simulation model that combines the surface water flow model POTYLDR and the groundwater flow model MODFLOW was used to study the impacts of watershed structures (e.g., dams) and irrigation water use (including stream-aquifer interactions) on streamflow and groundwater. The model was revised, enhanced, calibrated, and verified, then applied to evaluate the hydrologic budget for Wet Wal

Journal of Hydrology↗

Simulation of the effects of groundwater withdrawals on water-level altitudes in the Sparta aquifer in the Bayou Meto-Grand Prairie area of eastern Arkansas, 2007-37

A groundwater-flow model of the Mississippi embayment was used to evaluate changes in water-level altitudes before (scenario 1) and after (scenario 2) the addition of wells that simulate potential future pumping from the Sparta aquifer in the Bayou Meto-Grand Prairie area of eastern Arkansas for the 30-year period from 2007 through 2037. Water-level altitudes at six model cell locations from the two different scenarios were compared for the period 2007 through 2037. Potential future pumping wells were added to the Mississippi Embayment Regional Aquifer Study model at a rate of 13 wells per year within areas of potential future pumping. Change maps for the Bayou Meto-Grand Prairie area were constructed for each scenario and water-level hydrographs were constructed for each scenario for each of the six model cell locations. The additional pumping from wells in the Sparta aquifer created greater water-level declines in the Bayou Meto-Grand Prairie area. In scenario 1, simulated water-level altitude declines range from 20 to 40 feet from 2007 through 2037. In scenario 2, the cone of depression in Lonoke County is the deepest, with a maximum water-level decline of approximately 102 feet. Water-level altitude declines range from 40 to 50 feet over most of the remainder of the Bayou Meto-Grand Prairie area in scenario 2. Simulated water-level altitudes across the Bayou Meto-Grand Prairie area and at all six model cell locations indicate substantial declines when additional wells pumping from the Sparta aquifer are introduced into the model from 2007 through 2037.

Arkansas↗

A preliminary evaluation of hydrologic conditions of the Lakeland Ridge area of Polk County, Florida

The Lakeland ridge are covers about 300 square miles in north-west Polk County. The rapid growth of this area has resulted in an increase in ground-water withdrawals such that in 1968 about 67.5 million gallons per day was pumped to satisfy the demands of municipal, irrigation and industrial users. Declines of water-levels of as much as 30 feet in the area cause pumps to lose suction and increase the hazard of upward migration of saline water. However, no widespread deterioration of water quality has been noted to date. Considerable hydrologic data are available for the Lakeland ridge area. However, in order to give water managers a proper base for decisions that would allow for optimum utilization of the groundwater resources, more detailed data are needed about the hydraulic characteristics of the aquifer, the quantity of groundwater withdrawals, and the quality of water at various depths in the Floridan aquifer.

Florida↗

An inset groundwater-flow model to evaluate the effects of layering configuration on model calibration and assess managed aquifer recharge near Shellmound, Mississippi

The U.S. Geological Survey has developed a high-resolution inset groundwater-flow model in the Mississippi Delta as part of an interdisciplinary collaboration coordinated by the Mississippi Alluvial Plain project to provide a tool that stakeholders can use to support water-resource management decisions. Groundwater withdrawals from the Mississippi River Valley alluvial (MRVA) aquifer have been vital to support agricultural production in the region, but substantial groundwater-level declines near Shellmound, Mississippi, have caused concerns for long-term sustainability of the aquifer. To better understand the subsurface and try to mitigate the long-term groundwater-level declines, stakeholders have undertaken actions including a Groundwater Transfer and Injection Pilot (GTIP) project using a riverbank filtration-based managed aquifer recharge approach. The pilot project consisted of extracting groundwater near the Tallahatchie River and reinjecting it into the aquifer 3 kilometers west where water levels have substantially declined. A high-resolution airborne electromagnetic (AEM) survey was also completed to collect electrical resistivity data to support the GTIP project and the development of the groundwater model. The inset groundwater-flow model was developed to (1) integrate the AEM data into the optimal layering configuration of the MRVA aquifer that the available observation data can support through calibration, and (2) assess the potential effect of the GTIP project on the groundwater levels. The AEM data were processed into three different layering configurations leading to the development of model A (18 layers), model B (16 layers), and model C (8 layers), all at a 100- x 100-meter cell spatial resolution using the U.S. Geological Survey modular finite-difference flow model 6 code with Newton-Raphson formulation. The model development process integrated recent advances in modeling, such as the incorporation of AEM data, the use of outputs from the soil-water-balance (SWB) model, and the Aquaculture and Irrigation Water-Use Model, and was facilitated by robust automation using the open-source python packages Modflow-setup and SFRmaker. Using Parameter Estimation ++ Iterative Ensemble Smoother, the three numerical groundwater-flow models (models A, B, and C) were calibrated against a set of observations, which included aquifer groundwater levels, streamflows, stream stage, and aquifer transmissivity. Results indicate that the detailed representation of MRVA aquifer layers in model A produced the best calibrated model by history matching, and the integration of data representing surficial connectivity played a key role in improving groundwater recharge and enhancing the ability of the model to match groundwater levels in the cone of depression. A forecast model simulated the managed aquifer recharge approach, and the results indicated that, given average irrigation and recharge conditions (2010–15), the GTIP project has the potential to induce groundwater-level increases of as much as 3 meters around the injection site, but a sustained increase would require repetition in subsequent years of water transfer at 2022 rates or above.

Mississippi↗