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 1,423 records · Page 79Linked to original sources

Hydrogeology of aquifers in Cretaceous and younger rocks in the vicinity of Onslow and southern Jones counties, North Carolina

Unconsolidated sediments in Onslow and Jones Counties, North Carolina overlie crystalline basement rocks and range in thickness from about 700 ft to more than 1,800 ft, thickening toward the east. This material is composed of permeable sand and limestone interlayered with relatively impermeable clay and silt beds. Sediments are divided into two groups: aquifers in Quaternary-, and Tertiary-aged rocks and aquifers in Cretaceous-aged rocks. Aquifers in the Cretaceous rocks provide most of the groundwater for public supplies and are the focus of this report. The aquifers in Cretaceous rocks are the Peedee, Black Creek, upper Cape Fear, and lower Cape Fear aquifers, which are composed of beds or groups of beds of sand and gravel. Each aquifer is overlain by a clay and silt bed, that impedes the flow of water between aquifers. The thickness of Cretaceous hydrogeologic units ranges from about 700 ft to more than 1, 300 ft. Hydrogeologic units are correlated using 60 geophysical logs and accompanying drillers ' logs along with water level and water quality data. Three hydrogeologic sections demonstrate the continuity of the aquifers and confining units, show water levels and chloride concentration in water from test intervals, and delineate where chloride concentration in water exceeds 250 mg/L within each aquifer. Maps of each aquifer in Cretaceous rocks show altitude of its top, thickness, sand percentage, and the transition from freshwater to saltwater. Maps of the confining units show thickness and sand percentage of each.

North Carolina↗

Multiple technologies applied to characterization of the porosity and permeability of the Biscayne aquifer, Florida

Research is needed to determine how seepage-control actions planned by the Comprehensive Everglades Restoration Plan (CERP) will affect recharge, groundwater flow, and discharge within the dual-porosity karstic Biscayne aquifer where it extends eastward from the Everglades to Biscayne Bay. A key issue is whether the plan can be accomplished without causing urban flooding in adjacent populated areas and diminishing coastal freshwater flow needed in the restoration of the ecologic systems. Predictive simulation of groundwater flow is a prudent approach to understanding hydrologic change and potential ecologic impacts. A fundamental problem to simulation of karst groundwater flow is how best to represent aquifer heterogeneity. Currently, U.S. Geological Survey (USGS) researchers and academic partners are applying multiple innovative technologies to characterize the spatial distribution of porosity and permeability within the Biscayne aquifer.

Open-File Report↗

Hydrogeology, aquifer characteristics, and ground-water flow of the surficial aquifer system, Broward County, Florida

The surficial aquifer system, in which an unconfined groundwater flow system exists, comprises the sediments from land surface to the intermediate confining unit (formerly called the Floridan aquiclude) in Broward County, Florida. These sediments have hydraulic conductivities that range more than seven orders of magnitude from about 0.001 ft/d to more than 10,000 ft/d. The sediments are grouped into (1) the Biscayne aquifer, (2) a semiconfining unit, (3) a gray limestone aquifer in west Broward County, and (4) basal sand or clayey sand. The system is about 160 ft thick in west Broward County and more than 350 ft thick in east Broward. Transmissivities, locally variable, show a clear areal trend from greater than 300,000 sq ft/day in southeast Broward County to less than 75,000 sq ft/day in the northwest part of the county. Very high transmissivity is associated with the Biscayne aquifer. Transmissivity of the gray limestone aquifer ranged from about 20,000 to 88,000 sq ft/day. Topography, water levels, geologic framework, and water quality in west Broward County suggest past and present groundwater movement is south or southeast. Drainage and urban development have greatly altered circulation patterns in east Broward County , eliminating the original coastal groundwater ridge, causing rapid surface runoff and short groundwater flow paths, and saltwater intrusion. (USGS)

Florida↗

Hydrologic framework and characterization of the Little Colorado River alluvial aquifer near Leupp, Arizona

The Little Colorado River alluvial aquifer near Leupp, Arizona, was investigated as a possible source of irrigation water for the Leupp and Birdsprings Chapters of the Navajo Nation. The physical, chemical, and hydraulic characteristics of the alluvial aquifer were studied using geophysical surveys, installation of observation wells, water-level measurements, chemical analyses, groundwater pumping simulations, and review of previous investigations. Geophysical surveys and well borings revealed that the aquifer ranges in thickness from near 0 feet around its periphery to about 100 feet in its thickest parts. Water levels were monitored in nine alluvial wells within the study area and compared with earlier measurements collected in 1998 and 1999. Comparison of those earlier water levels with water levels collected as part of this study showed a decline of between about 5 and 12 feet has occurred in the last 23 years. The water chemistry of the aquifer was analyzed for salinity hazard, sodium-adsorption hazard, and specific-ion toxicity. The sodium-adsorption hazard and specific-ion toxicity of alluvial aquifer water were found to be low. However, the salinity hazard was high enough in most areas that it could negatively affect salt-sensitive crops. Well-field pumping scenarios conducted for this study demonstrated that using groundwater from the alluvial aquifer for irrigated agriculture is theoretically possible but may be economically challenging owing to the hydraulic properties of the aquifer.

Arizona↗

Hydrogeologic framework and hydrologic conditions of the Piney Point aquifer in Virginia

The Piney Point aquifer in Virginia is newly described and delineated as being composed of six geologic units, in a study conducted by the U.S. Geological Survey in cooperation with the Virginia Department of Environmental Quality (VA DEQ). The eastward-dipping geologic units include, in stratigraphically ascending order, the sand of the Nanjemoy Formation Woodstock Member, interbedded limestone and sand of the Piney Point Formation, silty and clayey sand of the Gosport Formation equivalent sediments, silty sand of the Oligocene-age sediments, silty fine-grained sand of the Old Church Formation, and silty sand of the Calvert Formation, Newport News unit and basal Plum Point Member. Identification of geologic units is based on typical sediment lithologies of geologic formations. Fine-grained sediments that compose confining units positioned immediately above and below the Piney Point aquifer are also described. The Piney Point aquifer is one of several confined aquifers within the Virginia Coastal Plain and includes a highly porous and solution-channeled indurated limestone within the Piney Point Formation from which withdrawals are made. The limestone is relatively continuous laterally across central parts of the Northern Neck, Middle Peninsula, and York-James Peninsula. Other geologic units are of variable extent. The configurations of most of the geologic units are further affected by newly identified faults that are aligned radially from the Chesapeake Bay impact crater and create constrictions or barriers to groundwater flow. Some geologic units are also truncated beneath the lower Rappahannock River by a resurge channel associated with the impact crater. Groundwater withdrawals from the Piney Point aquifer increased from approximately 1 million gallons per day (Mgal/d) during 1900 to 7.35 Mgal/d during 2004. As a result, a water-level cone of depression in James City and northern York Counties was estimated to be as low as 70 feet (ft) below the National Geodetic Vertical Datum of 1929 (NGVD 29) by 2005. Withdrawals decreased to 5.01 Mgal/d by 2009 as withdrawals were shifted toward other sources, and by 2015 water levels had recovered to approximately 50 ft below NGVD 29. The mean estimated transmissivity of the Piney Point aquifer in York and James City Counties is 16,300 feet squared per day (ft 2 /d), but farther north it is only 925 ft 2 /d. The mean well specific capacity in York and James City Counties is 11.4 gallons per minute per foot (gal/min/ft). Farther north in Virginia, mean specific capacity is only 2.26 gal/min/ft, and in Maryland it is 0.99 gal/min/ft. The northward decrease in specific capacity probably reflects the northward decrease in transmissivity, which results from poor development of the solution-channeled limestone. An aquifer test in northern York County induced vertical leakage to the solution-channeled limestone from overlying silty sand and a change in response of the aquifer to pumping from a single layer to two layers. Transmissivity of the limestone of approximately 19,800 ft 2 /d was distinguished from the silty sand of approximately 2,500 ft 2 /d. Most of the water in the Piney Point aquifer is slightly alkaline with moderate concentrations primarily of sodium and bicarbonate that are slightly undersaturated with respect to calcite. Iron concentrations are generally less than 0.3 milligrams per liter (mg/L). Mixing of freshwater with seawater has elevated chloride concentrations to the southeast to as much as 7,120 mg/L. Information on the Piney Point aquifer can benefit water-resource management in siting production wells, predicting likely well yield, and anticipating water-level response to withdrawals. Models that vertically discretize individual geologic units can potentially be used to evaluate groundwater flow in greater detail by representing lateral flow and vertical leakage among the geologic units. Because groundwater withdrawals are made primarily from the limestone and sand of the Piney Point Formation, the VA DEQ has considered regarding the limestone and sand singly as a regulated aquifer apart from the other geologic units. Under current policy in Virginia, if only the limestone and sand were regarded as a regulated aquifer, a greater amount of drawdown would be allowed than is allowed for the Piney Point aquifer consisting of six geologic units. Some production wells intercept multiple geologic units, and the units can undergo water-level decline and vertical leakage induced by pumping from the limestone and sand. Whether the other geologic units are to be regarded as regulated aquifers is an additional consideration for the VA DEQ.

Virginia↗

Effects of surface-water and groundwater inflows and outflows on the hydrology of the Tsala Apopka Lake Basin in Citrus County, Florida

The U.S. Geological Survey, in cooperation with the Southwest Florida Water Management District, initiated a study to quantify the inflows and outflows in the Floral City, Inverness, and Hernando pools of the Tsala Apopka Lake Basin in Citrus County, Florida. This study assesses hydrologic changes in pool stages, groundwater levels, spring flows, and streamflows caused by the diversion of streamflow from the Withlacoochee River to the Tsala Apopka Lake Basin through water-control structures. A surface-water/groundwater flow model was developed using hydraulic parameters for lakes, streams, the unsaturated zone, and the underlying surficial and Upper Floridan aquifers estimated using an inverse modeling calibration technique. After calibration, the model was used to assess the relation between inflows and outflows in the Tsala Apopka Lake Basin and changes in pool stages. Simulation results using the calibrated surface-water/groundwater flow model showed that leakage rates from the pools to the Upper Floridan aquifer were largest at the deep lake cells and that these leakage rates to the Upper Floridan aquifer were the highest in the model area. Downward leakage to the Upper Floridan aquifer occurred beneath most of the extent of the Floral City, Inverness, and Hernando pools. These leakage rates depended on the lakebed leakance and the difference between lake stages and heads in the Upper Floridan aquifer. Leakage rates were higher for the Floral City pool than for the Inverness pool, and higher for the Inverness pool than for the Hernando pool. Lakebed leakance was higher for the Floral City pool than for the Hernando pool, and higher for the Hernando pool than for the Inverness pool. Simulation results showed that the average recharge rate to the surficial aquifer was 10.3 inches per year for the 2004 to 2012 simulation period. Areas that recharge the surficial aquifer covered about 86 percent of the model area. Simulations identified areas along segments of the Withlacoochee River and within land-surface depressions that receive water from the surficial aquifer. Recharge rates were largest in physiographic regions having a deep water table. Simulated heads in the Upper Floridan aquifer indicated the general flow directions in the active flow model area were from the northeast toward the southwest and then westward toward the coast, and from the southeast toward the northwest and then westward toward the coast, consistent with flow directions inferred from the estimated potentiometric surface map for May 2010. The largest inflow in the water budget of the Upper Floridan aquifer was downward leakage from the overlying hydrogeologic unit. The largest outflow in the water budget of the Upper Floridan aquifer was spring flow. The calibrated surface-water and groundwater flow model was used to simulate hydrologic scenarios that included changes in rainfall rates, projected increases in groundwater pumping rates for 2025 and 2035, no flow for the 2004–12 period through the eight water-control structures in the Tsala Apopka Lake Basin, and the removal of the Inglis Dam and the Inglis Bypass Spillway on Lake Rousseau. Scenario simulation results were compared to annual average calibrated water levels and flows from 2004 to 2012. Simulated declines in the Tsala Apopka Lake pool stages under the 10-percent lower rainfall scenario were about 0.8, 0.3, and 1.3 feet (ft) for the Floral City, Inverness, and Hernando pools, respectively. Simulated groundwater levels under the same scenario declined up to 5.4 ft in the surficial aquifer and up to 2.9 ft in the Upper Floridan aquifer. Under the projected increases in groundwater pumping rates for 2035 that represented an increase of 36 percent from average 2004 to 2012 pumping rates, the simulated declines in the Floral City, Inverness, and Hernando pool stages were, in downstream order, 0.02, 0.06, and 0.04 ft. The largest drawdown under the projected increases in groundwater pumping rates for 2035 was 2.1 ft in the surficial aquifer and about 1.8 ft in the Upper Floridan aquifer. A scenario of decreased rainfall by 10 percent caused greater declines in water levels and pool stages than projected increases in groundwater pumping rates. The simulation with no flow through the eight Tsala Apopka Lake water-control structures resulted in simulated declines in average pool stage of 1.8, 1.9, and 0.5 ft in the Floral City, Inverness, and Hernando pools, respectively. The simulated removal of the two water-control structures in Lake Rousseau caused flow to increase at Rainbow Springs by 28 cubic feet per second, an increase of 4.7 percent from the average calibrated flow for 2004 to 2012.

Florida↗

Ground-water resources and simulation of flow in aquifers containing freshwater and seawater, Island County, Washington

Aquifers in Island County, Washington, that are intruded by seawater from Puget Sound contain chloride concentrations that exceed 100 mg/L. Chloride concentrations exceeded 100 mg/L in 24 % of the wells that were drilled below sea level and sampled in August 1981, but most of the chloride concentrations did not exceed 1,000 mg/L. Groundwater occurs in glacial deposits that have a maximum thickness of 3,000 ft; the deposits were divided into five aquifers and five confining units. Four overlapping digital models were calibrated, using time-averaged data, to simulate three-dimensional steady flow of fresh groundwater in multiple aquifers containing freshwater and seawater separated by a sharp interface. Model simulations indicate that most of the recharge is discharged from aquifers C and D as springs below sea levels, and only a small fraction of the recharge moves downward below aquifer C. Simulations also indicate that aquifers beneath the islands are not recharged by groundwater that moves from the mainland through aquifers beneath Puget Sound except in the area of northeast Camano Island. Between Whidbey and Camano Islands, the freshwater-seawater interface intersects the bottom of Puget Sound and prevents movement of fresh groundwater between the two islands. (USGS)

Water-Resources Investigations Report↗

Land subsidence along the California Aqueduct in west-central San Joaquin Valley, California, 2003–10

Extensive groundwater withdrawal from the unconsolidated deposits in the San Joaquin Valley caused widespread aquifer-system compaction and resultant land subsidence from 1926 to 1970—locally exceeding 8.5 meters. The importation of surface water beginning in the early 1950s through the Delta-Mendota Canal and in the early 1970s through the California Aqueduct resulted in decreased groundwater pumping, recovery of water levels, and a reduced rate of compaction in some areas of the San Joaquin Valley. However, drought conditions during 1976–77, 1987–92, and drought conditions and operational reductions in surface-water deliveries during 2007–10 decreased surface-water availability, causing pumping to increase, water levels to decline, and renewed compaction. Land subsidence from this compaction has reduced freeboard and flow capacity of the California Aqueduct, Delta-Mendota Canal, and other canals that deliver irrigation water and transport floodwater. The U.S. Geological Survey, in cooperation with the California Department of Water Resources, assessed more recent land subsidence near a 145-kilometer reach of the California Aqueduct in the west-central part of the San Joaquin Valley as part of an effort to minimize future subsidence-related damages to the California Aqueduct. The location, magnitude, and stress regime of land-surface deformation during 2003–10 were determined by using data and analyses associated with extensometers, Global Positioning System surveys, Interferometric Synthetic Aperture Radar, spirit-leveling surveys, and groundwater wells. Comparison of continuous Global Positioning System, shallow-extensometer, and groundwater-level data indicated that most of the compaction in this area took place beneath the Corcoran Clay, the primary regional confining unit. The integration of measurements strengthens confidence in individual measurement methods and provides the information at spatial and temporal scales that water managers need to design and implement groundwater sustainability plans in compliance with California’s Sustainable Groundwater Management Act. Measurements of land-surface deformation during 2003–10 indicated that the parts of the California Aqueduct closest to the Coast Ranges in the west-central part of the San Joaquin Valley were fairly stable or minimally subsiding on an annual basis; some areas show seasonal periods of subsidence and uplift that resulted in little or no longer-term elevation loss. Many groundwater levels in these areas did not reach historical lows during 2003–10, indicating that deformation nearest the Coast Ranges was likely primarily elastic. Land-surface deformation measurements indicated that some parts of the California Aqueduct that traverse farther from the Coast Ranges toward the valley center subsided. Some parts of the California Aqueduct subsided locally, but generally the California Aqueduct is within part of a 12,000-square-kilometer area affected by 25 millimeters or more of subsidence during 2008–10, with maxima in Madera County, south of the town of El Nido near the San Joaquin River and the Eastside Bypass (540 millimeters), and in Tulare County, west of the town of Pixley (345 millimeters). Interferometric Synthetic Aperture Radar-derived subsidence maps for various periods during 2003–10 show that the area of maximum active subsidence (that is, the largest rates of subsidence) shifted from its historical (1926–70) location southwest of the town of Mendota to these areas nearer the valley center. Calculations indicated that the subsidence rate doubled in 2008 in parts of the study area. Water levels declined during 2007–10 in many shallow and deep wells in the most rapidly subsiding areas, where water levels in many deep wells reached their historical lows, indicating that subsidence measured during this period was largely inelastic. Continued groundwater-level and land-subsidence monitoring in the San Joaquin Valley is important because (1) operational- and drought-related reductions in surface-water deliveries since 1976 have resulted in increased groundwater pumping and associated water-level declines and land subsidence, (2) land use and associated pumping continue to change throughout the valley, and (3) subsidence management is stipulated in the Sustainable Groundwater Management Act. The availability of surface water remains uncertain; even during record-setting precipitation years, such as 2010–11, water deliveries fell short of requests and groundwater pumping was required to meet the irrigation demand. In some areas, the infrastructure is not available to supply surface water, and groundwater is the only source of water. Because of the expected continued demand for water and the limitations and uncertainty of surface-water supplies, groundwater pumping and associated land subsidence remains a concern. Spatially detailed information on land subsidence is needed to minimize future subsidence-related damages to the California Aqueduct and other infrastructure in the San Joaquin Valley, as well as alterations to natural resources such as stream gradients, water depths, and water temperatures. The integration of data on land-surface elevation, subsurface deformation, and water levels—particularly continuous measurements—enables the analysis of aquifer-system response to groundwater pumping, which in turn, enables estimation of the preconsolidation head and calculation of aquifer-system storage properties. This information can be used to improve numerical model simulations of groundwater flow and aquifer-system compaction and allow for consideration of land subsidence in the evaluation of water resource management alternatives and compliance with the Sustainable Groundwater Management Act.

California↗

Detection and measurement of land subsidence and uplift using Global Positioning System surveys and interferometric synthetic aperture radar, Coachella Valley, California, 2010–17

Groundwater has been a major source of agricultural, recreational, municipal, and domestic supply in the Coachella Valley of California since the early 1920s. Pumping of groundwater resulted in groundwater-level declines as large as 50 feet (ft) or 15 meters (m) by the late 1940s. Because of concerns that the declines could cause land subsidence, the Coachella Valley Water District (CVWD) and the U.S. Geological Survey (USGS) have cooperatively investigated subsidence in the Coachella Valley since 1996. Importation of Colorado River water to the southern Coachella Valley began in 1949, resulting in a reduction in groundwater pumping and a recovery of groundwater levels during the 1950s through the 1970s. Since the late 1970s, the demand for water in the valley increased to the point that groundwater levels again declined in response to increased pumping and, consequently, increased the potential for land subsidence caused by aquifer-system compaction. Several management actions to increase recharge or to reduce reliance on groundwater have been implemented since as early as 1973 to address overdraft in the Coachella Valley. The implementation of three particular projects has markedly improved groundwater conditions in some of the historically most overdrafted areas of the valley: (1) groundwater substitution with surface-water imports since 2006 using Colorado River water through the Mid-Valley Pipeline project, which was expanded through 2017; (2) budget-based, tiered rates since 2009; and (3) managed aquifer recharge at the Thomas E. Levy Groundwater Replenishment Facility since 2009. Global Positioning System (GPS) surveying and interferometric synthetic aperture radar (InSAR) methods were used to determine the location, extent, and magnitude of the vertical land-surface changes in the Coachella Valley during 2010–17, updating 1993–2010 information presented in previous USGS reports. The GPS measurements taken at 24 geodetic monuments in August 2010 and September 2015 indicated that the land-surface elevation was stable at 17 monuments but changed at seven monuments during the 5-year period. Subsidence ranged from 0.17 to 0.43 ±0.09 ft (52 to 132 ±28 millimeters, or mm) at three monuments, and uplift ranged from 0.11 to 0.18 ±0.09 ft (33 to 54 ±28 mm) at four monuments between 2010 and 2015. At two of the monuments that subsided, the subsidence rates decreased between 2010 and 2015 from those computed between 2005 and 2010. Data prior to 2010 were not available for the third monument that subsided; thus, the 2010–15 subsidence rate could not be compared to an earlier period. At three of the monuments that uplifted between 2010 and 2015, data collected in 2005 and 2010 indicated stability. Data prior to 2010 were not available for the fourth monument that uplifted; thus, the 2010–15 uplift rate could not be compared to an earlier period. InSAR analyses for December 28, 2014–June 27, 2017, indicated that the land surface uplifted as much as about 0.20 ft (60 mm) near the Whitewater River Groundwater Replenishment Facility in the northern Coachella Valley and subsided as much as about 0.26 ft (80 mm) in the La Quinta area and less in Palm Desert, Indian Wells, and other localized areas in the southern Coachella Valley. These areas were identified as subsidence areas in previous reports covering periods during 1993–2010. The comparison of 2014–17 subsidence rates with those derived for 1995–2010 generally indicated a substantial slowing of subsidence, however. Analyses of deformation in the northern Coachella Valley were not included in the previous reports, so a comparison to deformation during the earlier period could not be made. Water levels in wells near the subsiding geodetic monuments, in and near the three subsiding areas shown by InSAR, and throughout the valley generally indicated seasonal fluctuations and longer-term stability or rising groundwater levels since about 2010. These results mark a reversal in trends of groundwater-level declines during the preceding decades. This trend reversal provides new insights into aquifer-system mechanics. Although many areas have stopped subsiding, and a few have even uplifted, the few areas that did subside during 2010–17—albeit at a slower rate—indicate a mixed aquifer-system response. Subsidence when groundwater levels are stable or recovering indicates that residual compaction may have occurred. At the same time, coarse-grained materials and thin aquitards may have expanded as groundwater levels recovered. The continued valley-wide stabilization and recovery of groundwater levels since 2010 likely is a result of various projects designed to increase recharge or to reduce reliance on groundwater.

California↗

Ground-water flow and solute transport at a municipal landfill site on Long Island, New York. Part 3, simulation of solute transport

A solute transport model representing a 2.3-sq mi area surrounding and downgradient from a municipal landfill site in the Town of Brookhaven, N.Y. was used to simulate migration of a conservative solute (chloride) in the upper glacial aquifer. Aquifer values used in the model were: hydraulic conductivity, 200 ft/day; effective porosity, 0.30; longitudinal dispersivity, 100 ft; transverse dispersivity, 20 ft. Average concentration of chloride was set at 875.0 mg/L in leachate and 10 mg/L in recharge and in ambient groundwater. Entry of leachate into the aquifer was assumed to have begun in 1977. Chloride concentrations in the simulated plume after 6 years of travel matched reasonably well the chloride data collected in October-December 1982. After 12 years of travel, the simulated plume extended 6,200 ft and was 2,600 ft wide. Maximum predicted concentration at the site boundary was 160 mg/L. Additional simulations were made to test the model 's ability to predict the effect of several remedial strategies on the movement of solutes. These included capping the landfill with an impermeable surface, removal of contaminated groundwater through four recovery wells, and a combination of the first two actions.

New York↗

Geospatial compilation of historical water-level altitudes in the Chicot and Evangeline aquifers 1977-2013 and Jasper aquifer 2000-13 in the Gulf Coast aquifer system, Houston-Galveston Region, Texas

The U.S. Geological Survey (USGS) in cooperation with the Harris-Galveston Subsidence District, City of Houston, Fort Bend Subsidence District, Lone Star Groundwater Conservation District, and Brazoria County Groundwater Conservation District has produced a series of annual reports depicting groundwater-level altitudes in the Chicot, Evangeline, and Jasper aquifers of the Gulf Coast aquifer system in the Houston-Galveston region, Texas. To produce these annual reports, contours of equal water-level altitudes are created from water levels measured between December and March of each year from groundwater wells screened completely within one of these three aquifers. Information obtained from maps published in the annual series of USGS reports and geospatial datasets of water-level altitude contours used to create the annual series of USGS reports were compiled into a comprehensive geodatabase. The geospatial compilation contains 88 datasets from previously published contour maps showing water-level altitudes for each primary aquifer of the Gulf Coast aquifer system, 37 for the Chicot (1977–2013), 37 for the Evangeline aquifer (1977–2013), and 14 for the Jasper aquifer (2000–13). Maps were georeferenced and digitized where existing geographic information system (GIS) data were unavailable (1977–89, 1991, 1995–99). Existing GIS data available for 1990, 1992–94, and 2000–13 were included in the geodatabase. The feature classes were organized into three feature datasets by principal aquifer: Chicot, Evangeline, and Jasper aquifers.

Texas↗

Simulation of the flow system of Barton Springs and associated Edwards Aquifer in the Austin area, Texas

A digital model of two-dimensional groundwater flow was used to estimate the hydraulic properties of the Edwards Aquifer in a 151 sq mi area near Austin, Texas. The transmissivity, hydraulic conductivity, and specific yield were estimated for the part of the aquifer that discharges at Barton Springs in Austin. The aquifer is composed of the Edwards and overlying Georgetown Limestones of Cretaceous age and ranges in thickness from about 100 ft to about 450 ft. More than 60 years of discharge measurements and 5 years of gaged discharge for Barton Springs were used to adjust springflow for the simulations. Barton Springs accounts for about 96% of springflow from the study area and 90% of the total discharge. The remaining discharge was pumpage from wells which was entered in the model. Four years of gaged recharge were used in the simulations. The potentiometric surfaces used by the models were constructed from water level measurements in as many as 75 wells. The transmissivity was calibrated through steady-state simulations that used the mean value of recharge and mean potentiometric surface to represent average conditions for the aquifer. The transmissivities vary from about 100 sq ft/day in the western part of the aquifer to > 1 million sq ft/day near Barton Springs. Specific yield was calibrated through transient-state simulations for 5 consecutive months using time-dependent data for recharge, discharge, and water levels. The mean specific yield for the aquifer is 0.014 and ranges from 0.008 to 0.064. Additional aquifer properties used in the simulations include storage coefficient, altitudes of the base and top of the aquifer, and hydraulic conductivity. A simulation for the year 2000 using projected pumping rates for municipal, industrial, agricultural, and domestic supplies indicates that the aquifer would be dewatered in the southwestern part of the study area and have large declines in the southeastern part of the study area. Another simulation of projected conditions using potential recharge enhancement predicts a rise in the potentiometric surface of about 50 feet in the southwestern part of the aquifer and moderate water-level declines in the southeastern part of the aquifer.

Texas↗

Impact of land use and land cover change on groundwater recharge and quality in the southwestern US

Humans have exerted large‐scale changes on the terrestrial biosphere, primarily through agriculture; however, the impacts of such changes on the hydrologic cycle are poorly understood. The purpose of this study was to test the hypothesis that the conversion of natural rangeland ecosystems to agricultural ecosystems impacts the subsurface portion of the hydrologic cycle by changing groundwater recharge and flushing salts to underlying aquifers. The hypothesis was examined through point and areal studies investigating the effects of land use/land cover (LU/LC) changes on groundwater recharge and solute transport in the Amargosa Desert (AD) in Nevada and in the High Plains (HP) in Texas, US. Studies use the fact that matric (pore‐water‐pressure) potential and environmental‐tracer profiles in thick unsaturated zones archive past changes in recharging fluxes. Results show that recharge is related to LU/LC as follows: discharge through evapotranspiration (i.e., no recharge; upward fluxes <0.1 mm yr −1 ) in natural rangeland ecosystems (low matric potentials; high chloride and nitrate concentrations); moderate‐to‐high recharge in irrigated agricultural ecosystems (high matric potentials; low‐to‐moderate chloride and nitrate concentrations) (AD recharge: ∼130–640 mm yr −1 ); and moderate recharge in nonirrigated (dryland) agricultural ecosystems (high matric potentials; low chloride and nitrate concentrations, and increasing groundwater levels) (HP recharge: ∼9–32 mm yr −1 ). Replacement of rangeland with agriculture changed flow directions from upward (discharge) to downward (recharge). Recent replacement of rangeland with irrigated ecosystems was documented through downward displacement of chloride and nitrate fronts. Thick unsaturated zones contain a reservoir of salts that are readily mobilized under increased recharge related to LU/LC changes, potentially degrading groundwater quality. Sustainable land use requires quantitative knowledge of the linkages between ecosystem change, recharge, and groundwater quality.

Global Change Biology↗

Hydrogeologic framework and extent of saltwater intrusion in Kings, Queens, and Nassau Counties, Long Island, New York

In 2016, the U.S. Geological Survey began a multiyear cooperative study with the New York State Department of Environmental Conservation to evaluate the sustainability of Long Island’s sole-source aquifer system through hydrogeologic mapping, compilation of groundwater chloride concentrations, and groundwater flow modeling. In the initial phase of the islandwide study, the hydrogeologic framework and extent of saltwater intrusion in aquifers in Kings, Queens, and Nassau Counties on western Long Island, N.Y., were investigated. The aquifer system underlying western Long Island has been under stress from pumping of public, irrigation (golf course), and industrial supply wells. Saltwater intrusion has occurred from surrounding embayments (East River, Long Island Sound, Jamaica Bay, and the Atlantic Ocean) due to pumping. Eighteen boreholes were drilled and cores taken during 2019–21 to collect hydrogeologic, geochemical, and geophysical data to delineate the complex subsurface hydrogeology and extent of saltwater intrusion within the study area. Evaluation of the new cores, reexamination of legacy core descriptions, and analysis of borehole geophysical logs was used to refine the previously published hydrogeologic framework of Pleistocene and Cretaceous unconsolidated sediments in the area, including delineation of a previously undefined hydrogeologic unit between the Magothy aquifer and the Raritan confining unit, herein named the “upper Raritan aquifer.” The upper Raritan aquifer was first recognized in southeastern Nassau County from an analysis of about 50 closely spaced boreholes with high-resolution core descriptions and gamma-ray (gamma) logs. Further analysis of borehole logs across the study area indicated that the upper Raritan aquifer was also present in Kings and Queens Counties. Nuclear magnetic resonance (NMR) logging was used for the first time on Long Island to provide estimates of the hydraulic properties of the major aquifer and confining units. Unlike other geophysical logs that record responses to the rock matrix and fluid properties and are strongly dependent on mineralogy, NMR logs record responses to the presence of hydrogen protons in the formation fluid to determine water fraction and pore-size distribution. NMR log analysis provided estimates of the clay-bound, capillary-bound, and mobile water fractions and hydraulic conductivity of aquifers and confining units penetrated by five wells in Nassau County. Pumpage for public-supply and industrial wells on Long Island began in the 1870s with small, localized suppliers of populated areas in Kings and Queens Counties. By 1904–16, pumpage for public water supply in Kings County averaged 21 million gallons per day, and averaged 37 million gallons per day in Queens County, mostly from the upper glacial aquifer. Saltwater intrusion was reported as early as the beginning of the 20th century and included the upper glacial-Jameco-Magothy and Lloyd-North Shore aquifer systems. By 1936, pumping in central Kings County created a major cone of depression in the water table extending to the south shore of much of Kings County and into southwestern Queens County. Saltwater intrusion has caused the shutdown of public-supply wells in Kings, Queens, and Nassau Counties. A large saltwater intrusion wedge in the Lloyd aquifer was indicated in southern Queens County in the early part of the 20th century, and the saltwater interface may have been onshore predevelopment. Most of Kings and Queens Counties are intruded with saltwater in both the upper glacial-Jameco-Magothy and Lloyd-North Shore aquifers systems. Saltwater increased during the 20th century and continues to increase to the present (2023) in the Lloyd-North Shore aquifer system in Great Neck and Manhasset Neck in northern Nassau County. A major wedge of saltwater intrusion in the upper glacial-Jameco-Magothy aquifer in southwestern Nassau County appears to be increasing.

New York↗

Water levels and water quality in the Sparta-Memphis aquifer (middle Claiborne aquifer) in Arkansas, spring-summer 2009

The U.S. Geological Survey in cooperation with the Arkansas Natural Resources Commission and the Arkansas Geological Survey has monitored water levels in the Sparta Sand of Claiborne Group and Memphis Sand of Claiborne Group (herein referred to as the Sparta Sand and the Memphis Sand, respectively) since the 1920s. Groundwater withdrawals have increased while water levels have declined since monitoring was initiated. Herein, aquifers in the Sparta Sand and Memphis Sand will be referred to as the Sparta-Memphis aquifer throughout Arkansas. During the spring of 2009, 324 water levels were measured in wells completed in the Sparta-Memphis aquifer and used to produce a regional potentiometric-surface map. During the summer of 2009, 64 water-quality samples were collected and measured for specific conductance, temperature, and pH from wells completed in the Sparta-Memphis aquifer. The regional direction of groundwater flow in the Sparta-Memphis aquifer is generally to the south-southeast in the northern half of Arkansas and to the east and south in the southern half of Arkansas, away from the outcrop area except where affected by large groundwater withdrawals. The highest and lowest water-level altitudes measured in the Sparta-Memphis aquifer were 325 feet above and 157 feet below National Geodetic Vertical Datum of 1929, respectively. Eight depressions (generally represented by closed contours) are located in the following counties: Bradley; Ashley; Calhoun; Cleveland; Columbia; Arkansas, Jefferson, Lincoln, and Prairie; Cross and Poinsett; and Union. Two large depressions shown on the 2009 potentiometric-surface map, centered in Jefferson and Union Counties, are the result of large withdrawals for industrial, irrigation, or public supply. The depression centered in Jefferson County deepened and expanded in recent years into Arkansas and Prairie Counties. The area enclosed within the 40-foot contour on the 2009 potentiometric-surface map has expanded south to the Drew County line and moved west from the intersection of Arkansas, Jefferson, and Lincoln Counties when compared with the 2007 potentiometric-surface map. To the north, east, and west, the 40-foot contour is comparable to the 2007 potentiometric-surface map. The lowest water-level altitude measurement during 2009 in the center of the depression in Union County represents a rise of 42 feet since 2003. The area enclosed by the lowest altitude contour, 140 feet below National Geodetic Vertical Datum of 1929, on the 2009 potentiometric-surface map is about half the area on the 2007 potentiometric-surface map. In the depression in western Poinsett and Cross Counties, the 140-foot contour extended north to the Poinsett-Craighead County line and south across Cross County about two-thirds of the distance to the St. Francis County line. A water-level difference map was constructed using water-level measurements made during 2005 and 2009 from 309 wells. The difference in water level between 2005 and 2009 ranged from -74.6 to 60.2 feet. Areas with a general rise in water levels occur in central Columbia County, southern Jefferson County, and most of Union County. In the area around west-central Union County, water levels rose as much as 60.2 feet with water levels in 18 wells rising 20 feet or more, representing an average annual rise of 5 feet or more. Water levels generally declined throughout most of the rest of Arkansas. Hydrographs were constructed using a minimum of 25 years of water-level measurements at each of 206 wells. During the period 1985–2009, mean annual water levels rose in Calhoun, Columbia, Lafayette, and Union Counties, about 1.3 feet per year (ft/yr), 0.2 ft/yr, 0.1 ft/yr, and 0.6 ft/yr, respectively. Mean annual water-level declines between 0.0 and 2.3 ft/yr occurred in all other counties. In western Arkansas County, water-level altitudes in a continuously monitored well declined 60 feet during the irrigation season (April to September). Specific conductance ranged from 43 microsiemens per centimeter at 25 degrees Celsius (μS/cm) in Ouachita County to 1,230 μS/cm in Phillips County. The mean specific conductance was 392 μS/cm. Although there is a regional increase in specific conductance to the east and south, specific conductance values greater than 700 μS/cm occurred in samples from wells in Arkansas, Ashley, Monroe, Phillips, and Union Counties.

Arkansas↗

Assessing the vulnerability of public-supply wells to contamination&mdash;Edwards aquifer near San Antonio, Texas

This fact sheet highlights findings from the vulnerability study of a public-supply well field in San Antonio, Texas. The well field consists of six production wells that tap the Edwards aquifer. Typically, one or two wells are pumped at a time, yielding an average total of 20-21 million gallons per day. Water samples were collected from public-supply wells in the well field and from monitoring wells installed along general directions of flow to the well field. Samples from the well field contained some constituents of concern for drinking-water quality, including nitrate; the pesticide compounds atrazine, deethylatrazine, and simazine; and the volatile organic compounds tetrachloroethene (also called perchloroethene, or PCE), chloroform, bromoform, and dibromochloromethane. These constituents were detected in untreated water at concentrations much less than established drinking-water standards, where such standards exist. Overall, the study findings point to four primary factors that affect the movement and fate of contaminants and the vulnerability of the public-supply well field in San Antonio, Texas: (1) groundwater age (how long ago water entered, or recharged, the aquifer), (2) fast pathways for flow of groundwater through features formed or enlarged by dissolution of bedrock, (3) recharge characteristics of the aquifer, and (4) natural geochemical processes within the aquifer. A computer-model simulation of groundwater flow and transport was used to estimate the traveltime (or age) of water particles entering public-supply well W4 in the well field. Modeled findings show that almost half of the water reaching the public-supply well is less than 2 years old. Such a large percentage of very young water indicates that (1) contaminants entering the aquifer may be transported rapidly to the well, (2) there is limited time for chemical reactions to occur in the aquifer that may attenuate contaminants, and (3) should recharge water become contaminated with pathogenic microorganisms (which have limited survival times in aquifers), the microorganisms may be able to persist to the well. Features formed or enlarged by dissolution of bedrock allow most of the water reaching the well field to travel rapidly from the recharge zone to the supply wells along fast pathways rather than through the aquifer matrix. Supporting evidence includes (1) geophysical logging and flowmeter measurements in public-supply well W4 and in nearby monitoring wells showing that most of the flow volume into and out of the wells occurs in three horizontal zones, thought to be dissolution-enlarged bedding planes; and (2) fluctuations in groundwater chemistry that can be correlated to individual precipitation events. Analysis of water samples collected from shallow, intermediate, and deep zones of the Edwards aquifer at public-supply well W4 and from nearby monitoring wells reveal that water in the vicinity of the selected well field is notably well mixed throughout the sampled thickness of the Edwards aquifer, showing little of the chemical variation with depth that is commonly seen in other aquifers. Contaminants were found at all depths, and they did not enter the well through a specific horizon. The well-mixed nature of the Edwards aquifer is caused by the recharge characteristics of the area combined with fast flow paths through karst features. Constituents of concern in the Edwards aquifer for the long-term sustainability of the groundwater resource include the nutrient nitrate and anthropogenic contaminants such as atrazine, PCE, and chloroform. A scenario of hypothetical contaminant loading in the aquifer recharge zone was evaluated by using results from groundwater-flow-model particle tracking to assess the response of the aquifer to potential contamination. Results indicate that the concentrations at public-supply well W4 would begin to respond to contaminant loading in the recharge zone within 1 year because of short traveltimes through fast flow paths. Within 10 years, contaminant concentrations in the public-supply well would be equal to 90 percent of the input concentration for a contaminant (such as nitrate) that does not degrade in the oxic conditions of the Edwards aquifer.

Texas↗

Multilevel groundwater monitoring of hydraulic head and temperature in the eastern Snake River Plain aquifer, Idaho National Laboratory, Idaho, 2007-08

During 2007 and 2008, the U.S. Geological Survey, in cooperation with the U.S. Department of Energy, collected quarterly depth-discrete measurements of fluid pressure and temperature in six boreholes located in the eastern Snake River Plain aquifer of Idaho. Each borehole was instrumented with a multilevel monitoring system consisting of a series of valved measurement ports, packer bladders, casing segments, and couplers. Hydraulic heads (head) and water temperatures in boreholes were monitored at 86 hydraulically-isolated depth intervals located 448.0 to 1,377.6 feet below land surface. The calculation of head is most sensitive to fluid pressure and the altitude of the pressure transducer at each port coupling; it is least sensitive to barometric pressure and water temperature. An analysis of errors associated with the head calculation determined the accuracy of an individual head measurement at +/- 2.3 feet. Many of the sources of measurement error are diminished when considering the differences between two closely-spaced readings of head; therefore, a +/- 0.1 foot measurement accuracy was assumed for vertical head differences (and gradients) calculated between adjacent monitoring zones. Vertical head and temperature profiles were unique to each borehole, and were characteristic of the heterogeneity and anisotropy of the eastern Snake River Plain aquifer. The vertical hydraulic gradients in each borehole remained relatively constant over time with minimum Pearson correlation coefficients between head profiles ranging from 0.72 at borehole USGS 103 to 1.00 at boreholes USGS 133 and MIDDLE 2051. Major inflections in the head profiles almost always coincided with low permeability sediment layers. The presence of a sediment layer, however, was insufficient for identifying the location of a major head change in a borehole. The vertical hydraulic gradients were defined for the major inflections in the head profiles and were as much as 2.2 feet per foot. Head gradients generally were downward in boreholes USGS 133, 134, and MIDDLE 2050A, zero in boreholes USGS 103 and 132, and exhibited a reversal in direction in borehole MIDDLE 2051. Water temperatures in all boreholes ranged from 10.2 to 16.3 degrees Celsius. Boreholes USGS 103 and 132 are in an area of concentrated volcanic vents and fissures, and measurements show water temperature decreasing with depth. All other measurements in boreholes show water temperature increasing with depth. A comparison among boreholes of the normalized mean head over time indicates a moderately positive correlation.

Idaho↗

Towards reproducible environmental modeling for decision support: A worked example

A fully worked example of decision-support-scale uncertainty quantification (UQ) and parameter estimation (PE) is presented. The analyses are implemented for an existing groundwater flow model of the Edwards aquifer, Texas, USA, and are completed in a script-based workflow that strives to be transparent and reproducible. High-dimensional PE is used to history-match simulated outputs to corresponding state observations of spring flow and groundwater level. Then a hindcast of a historical drought is made. Using available state observations recorded during drought conditions, the combined UQ and PE analyses are shown to yield an ensemble of model results that bracket the observed hydrologic responses. All files and scripts used for the analyses are placed in the public domain to serve as a template for other practitioners who are interested in undertaking these types of analyses.

Texas↗