USGS Science⌕ Search

Geology topics

Groundwater & aquifers

Find USGS groundwater and aquifer publications. Search an aquifer name or place, compare source-supplied dates and open reports for hydrogeology, recharge and water-level studies.

At least 433 records · Page 24Linked to original sources

Water-quality comparison of the Gulf Coast aquifer system and Carrizo-Wilcox aquifer in Texas from National Water-Quality Assessment Project Principal Aquifer Surveys, 2013 and 2015

The U.S. Geological Survey’s National Water-Quality Assessment (NAWQA) Project assessed the quality of groundwater in aquifers that are important sources of drinking water in the United States. One major aquifer in Texas that was assessed by NAWQA in 2013 is the coastal lowlands aquifer system, which is often referred to in Texas as the “Gulf Coast aquifer system.” The coastal lowlands aquifer system supplies water for millions of people; self-supplied (private) well withdrawals in 2005 from this aquifer system were the sixth largest among all major aquifer systems in the Nation. A major aquifer in Texas that was assessed by NAWQA in 2015 is the Texas coastal uplands aquifer system; the Carrizo-Wilcox aquifer is one of several aquifers that compose this aquifer system in Texas. The rocks composing the Texas coastal uplands aquifer system extend east from Texas as part of the Mississippi embayment aquifer system and underlie areas of several States. The Texas coastal uplands aquifer system and Mississippi embayment aquifer system are often collectively referred to as the “Mississippi embayment-Texas coastal uplands aquifer system.” Self-supplied withdrawals from the Mississippi embayment-Texas coastal uplands aquifer system in 2005 were the eighth largest among all major aquifer systems in the Nation. The coastal lowlands aquifer system and Mississippi embayment-Texas coastal uplands aquifer system were assessed as part of the NAWQA Principal Aquifer Surveys (PAS), which were designed to evaluate constituent concentrations in water samples obtained from domestic and public-supply wells prior to any treatment. PAS assessments like these allow for the comparison of water-quality concentrations in untreated groundwater using preestablished benchmarks for the protection of human health and for aesthetic qualities such as taste, color, and odor. The use of preestablished benchmarks can provide a basis for comparison of groundwater quality among principal aquifers.

Texas↗

Groundwater hydrology and estimation of horizontal groundwater flux from the Rio Grande at selected locations in Albuquerque, New Mexico, 2009–10

The Albuquerque area of New Mexico has two principal sources of water: (1) groundwater from the Santa Fe Group aquifer system, and (2) surface water from the Rio Grande. From 1960 to 2002, pumping from the Santa Fe Group aquifer system caused groundwater levels to decline more than 120 feet while water-level declines along the Rio Grande in Albuquerque were generally less than 40 feet. These differences in water-level declines in the Albuquerque area have resulted in a great deal of interest in quantifying the river-aquifer interaction associated with the Rio Grande. In 2003, the U.S. Geological Survey, in cooperation with the Bureau of Reclamation, acting as fiscal agent for the Middle Rio Grande Endangered Species Collaborative Program, and the U.S. Army Corps of Engineers, began a study to characterize the hydrogeology of the Rio Grande inner valley alluvial aquifer in the Albuquerque area of New Mexico. The study provides hydrologic data in order to enhance the understanding of rates of water leakage from the Rio Grande to the alluvial aquifer, groundwater flow through the aquifer, and discharge of water from the aquifer to riverside drains. The study area extends about 20 miles along the Rio Grande in the Albuquerque area. Piezometers and surface-water gages were installed in paired transects at eight locations. Nested piezometers, completed at various depths in the alluvial aquifer, and surface-water gages, installed in the Rio Grande and riverside drains, were instrumented with pressure transducers. Water-level and water-temperature data were collected from 2009 to 2010. Water levels from the piezometers indicated that groundwater movement was usually away from the river towards the riverside drains. Annual mean horizontal groundwater gradients in the inner valley alluvial aquifer ranged from 0.0024 (I-25 East) to 0.0144 (Pajarito East). The median hydraulic conductivity values of the inner valley alluvial aquifer, determined from slug tests, ranged from 30 feet per day (ft/d) (Montaño) to 120 ft/d (Central) for paired transects, with a median hydraulic conductivity for all transects of 50 ft/d. Daily mean groundwater fluxes from the river through the inner valley alluvial aquifer computed using Darcy’s Law and the slug test results ranged from about 0.01 ft/d (Montaño West) to between 1.0 and 2.0 ft/d (Central East). Median annual groundwater fluxes from the river through the inner valley alluvial aquifer determined using the Suzuki-Stallman method was greatest at Alameda East (0.50 ft/d) and lowest at Alameda West (0.25 ft/d). The results from both methods agreed reasonably well. Seepage investigations conducted by measuring discharge in the east and west riverside drains provided information for computing changes in flow within the drains and for evaluating results from Darcy’s Law and Suzuki-Stallman method flux calculations. Discharge measured in the east riverside drain between the Barelas Bridge and the I-25 bridge indicated that the flow in the east riverside drain increased by an average of 56.5 cubic feet per day per linear foot (ft 3 /d/ft) of drain. Discharge measured in the west riverside drain between the Central bridge and the I-25 bridge indicated that flow increased between west drain miles 0 and 4, an average of 53.8 ft 3 /d/ft of drain, and that flow increased between west drain miles 7 and 10, an average of 44.9 ft 3 /d/ft of drain. In comparison to the seepage measurement results, the groundwater fluxes from the river through the inner valley alluvial aquifer calculated from Darcy’s Law (q slug ) and by the Suzuki-Stallman method (q heat ) would account for 20–36 percent or 53–95 percent, respectively, of the total flow in the east riverside drain and 22–31 percent or 19–26 percent, respectively, of the total flow in the west drain. These results indicate that the drains likely also receive water from outside the inner valley. The spatial variability of horizontal hydraulic gradients and groundwater fluxes can be primarily attributed to variability in the distances between the river and riverside drains throughout the study area and geologic heterogeneities in the alluvial aquifer. Temporal variability in the water levels, which control the horizontal hydraulic gradients and fluxes between the Rio Grande and the riverside drains, can be primarily attributed to seasonal fluctuations in river stage and irrigation practices.

New Mexico↗

Hydrogeology and simulation of groundwater flow in the Lucerne Valley groundwater basin, California

The Lucerne Valley is in the southwestern part of the Mojave Desert and is about 75 miles northeast of Los Angeles, California. The Lucerne Valley groundwater basin encompasses about 230 square miles and is separated from the Upper Mojave Valley groundwater basin by splays of the Helendale Fault. Since its settlement, groundwater has been the primary source of water for agricultural, industrial, municipal, and domestic uses. Groundwater withdrawal from pumping has exceeded the amount of water recharged to the basin, causing groundwater declines of more than 100 feet between 1917 and 2016 in the center of the basin. The continued withdrawal has resulted in an increase in pumping costs, reduced well efficiency, and land subsidence near Lucerne Lake. Although the volume of pumping has declined in recent years, there is concern that new agricultural growth and limits on imported water will continue to strain the sustainability of the groundwater system. To address these concerns, the U.S. Geological Survey entered into a cooperative agreement with the Mojave Water Agency to develop a better understanding of the Lucerne Valley hydrogeologic system and provide tools to help evaluate and manage the effects of future development in the Lucerne Valley. The objectives of this study were to (1) improve the understanding of the aquifer system, (2) improve the understanding of subsidence in the basin, and (3) incorporate the understanding into a groundwater-flow model that can be used to help manage the groundwater resources in the Lucerne Valley. The model developed for this study covers the period of 1942–2016 and can help evaluate various proposed water-management scenarios during different climatic and hydrologic conditions. The aquifer system consists of a shallow aquifer, a confining unit, and middle and lower aquifers. These layered water-bearing units were identified based on geologic units of the mostly unconsolidated sediments and hydrologic properties. These alluvial deposits consist of clay, silt, sand, and gravel; some places also contain clay and silty clay lacustrine deposits. Several faults act, at least in part, as barriers to groundwater flow on the eastern, southern, and western edges of the basin. Present-day natural recharge is primarily from the infiltration of runoff from the San Bernardino Mountains to the south; however, stable and radioactive isotopes show that groundwater from the middle of the Lucerne Valley was older than about 10,000 years and probably was recharged as infiltration from streams draining the mountains in the Mojave Desert to the north, which probably does not occur under present-day climatic conditions. The annual average natural recharge for 1942–2016, estimated by a Basin Characterization Model, was about 635 acre-feet per year; the average amount of treated wastewater effluent transferred to the Lucerne Valley for artificial recharge annually ranged from about 1,500 to 4,000 acre-feet per year during 1980–2016. Pumpage estimates for 1942–2016 ranged from about 3,000 acre-feet in 1942 to about 18,300 acre-feet in 1984. The total cumulative amount of groundwater removed from the basin by pumping between 1942 and 2016 was estimated to be about 700,000 acre-feet, which was about 10 times greater than the cumulative amount of recharge to the entire Lucerne Valley groundwater basin. Before groundwater development, the direction of groundwater flow was from the southern part of the basin northward to discharge areas near Lucerne Lake, where it discharged through springs along the Helendale Fault and by evapotranspiration. Since the early 1900s, groundwater-level declines have mostly eliminated the areas where natural discharge occurred and exceeded 100 feet in the middle of the basin between the early 1950s and mid-1990s, and as much as 25 feet near the margins from about the mid-1950s to 2000s. A decrease in the rate of pumping after the mid-1990s lessened the hydraulic stress on the middle and lower aquifers and enabled hydraulic heads in the middle of the basin to recover slightly as groundwater near the margins of the basin moved toward the pumping depression. Although trends in groundwater levels in the center of the basin have reversed since the mid-1990s, levels at the basin margins continue to decline as the movement of groundwater from the margins fills the pumping depression and gradually flattens the groundwater table throughout the basin. The long-term extraction of groundwater and associated dewatering of the fine-grained sediments present within the aquifer system has resulted in aquifer compaction and consequently land subsidence, primarily near Lucerne Lake. Analysis of interferometric synthetic aperture radar data shows that almost 11 inches of land subsidence has occurred south of Lucerne Lake between April 1992 and November 2009; less subsidence occurred elsewhere in the basin during this period. This differential land subsidence has caused fissures and cracks in the ground surface, which have buckled the pavement and undercut roads in several locations. The Lucerne Valley Hydrologic Model was developed using the finite-difference groundwater modeling software One Water Hydrologic Model to represent the hydrologic conditions and stresses during 1942–2016. The model has a uniform grid of approximately 92 acres per cell (2,000 feet by 2,000 feet) and has four layers representing the water-bearing units. The results from the calibrated model simulations indicated that groundwater pumpage exceeded recharge, resulting in an estimated net cumulative depletion of groundwater storage (discharge minus recharge) of about 465,000 acre-feet from 1942 to 2016. The model simulated as much as 7.5 feet (90 inches; 2,286 millimeters) of aquifer compaction, which indicates the extensive fine-grained deposits and measured subsidence near Lucerne Lake.

California↗

Effects of aquifer storage and recovery activities on water quality in the Little Arkansas River and Equus Beds Aquifer, south-central Kansas, 2011–14

The Equus Beds aquifer in south-central Kansas is aprimary water source for the city of Wichita. The Equus Beds aquifer storage and recovery (ASR) project was developed to help the city of Wichita meet increasing current (2016) and future water demands. The Equus Beds ASR project pumps water out of the Little Arkansas River during above-base flow conditions, treats it using drinking-water quality standards as a guideline, and recharges it into the Equus Beds aquifer for later use. Phase II of the Equus Beds ASR project currently (2016) includes a river intake facility and a surface-water treatment facility with a 30 million gallon per day capacity. Water diverted from the Little Arkansas River is delivered to an adjacent presedimentation basin for solids removal. Subsequently, waste from the surface-water treatment facility and the presedimentation basin is returned to the Little Arkansas River through a residuals return line. The U.S. Geological Survey, in cooperation with the city of Wichita, developed and implemented a hydrobiological monitoring program as part of the ASR project to characterize and quantify the effects of aquifer storage and recovery activities on the Little Arkansas River and Equus Beds aquifer water quality. Data were collected from 2 surface-water sites (one upstream and one downstream from the residuals return line), 1 residuals return line site, and 2 groundwater well sites (each having a shallow and deep part): the Little Arkansas River upstream from the ASR facility near Sedgwick, Kansas (upstream surface-water site 375350097262800), about 0.03 mile (mi) upstream from the residuals return line site; the Little Arkansas River near Sedgwick, Kans. (downstream surface-water site 07144100), about 1.68 mi downstream from the residuals return line site; discharge from the Little Arkansas River ASR facility near Sedgwick, Kansas (residuals return line site 375348097262800); 25S 01 W 07BCCC01 SMW–S11 near CW36 (MW–7 shallow groundwater well site 375327097285401); 25S01 W 07BCCC02 DMW–S10 near CW36 (MW–7 deep groundwater well site 375327097285402); 25S 01W 07BCCA01 SMW–S13 near CW36 (MW–8 shallow groundwater well site 375332097284801); and 25S 01W 07BCCA02 DMW–S14 near CW36 (MW–8 deep groundwater well site 375332097284802). The U.S. Geological Survey, in cooperation with the city of Wichita, assessed the effects of the ASR Phase II facility residuals return line discharges on stream quality of the Little Arkansas River by measuring continuous physicochemical properties and collecting discrete water-quality and sediment samples for about 2 years pre- (January 2011 through April 2013) and post-ASR (May 2013 through December 2014) Phase II facility operation upstream and downstream from the ASR Phase II facility. Additionally, habitat variables were quantified and macroinvertebrate and fish communities were sampled upstream and downstream from the ASR Phase II facility during the study period. To assess the effects of aquifer recharge on Equus Beds groundwater quality, continuous physicochemical properties were measured and discrete water-quality samples were collected before and during the onset of Phase II aquifer recharge in two (shallow and deep) groundwater wells. Little Arkansas River streamflow was about 10 times larger after the facility began operating because of greater rainfall. Residuals return line release volumes were a very minimal proportion (0.06 percent) of downstream streamflow volume during the months the ASR facility was operating. Upstream and downstream continuously measured water temperature and dissolved oxygen median differences were smaller post-ASR than pre-ASR. Turbidity generally was smaller at the downstream site throughout the study period and decreased at both sites after the ASR Phase II facility began discharging despite a median residuals return line turbidity that was about an order of magnitude larger than the median turbidity at the downstream site. Upstream and downstream continuously measured turbidity median differences were larger post-ASR than pre-ASR. Median post-ASR continuously measured nitrite plus nitrate and continuously computed total suspended solids and suspended-sediment concentrations were smaller than pre-ASR likely because of higher streamflows and dilution; whereas, median continuously computed dissolved and total organic carbon concentrations were larger likely because of higher streamflows and runoff conditions. None of the discretely measured water-quality constituents (dissolved and suspended solids, primary ions, suspended sediment, nutrients, carbon, trace elements, viral and bacterial indicators, and pesticides) in surface water were significantly different between the upstream and downstream sites after the ASR Phase II facility began discharging; however, pre-ASR calcium, sodium, hardness, manganese, and arsenate concentrations were significantly larger at the upstream site, which indicates that some water-quality conditions at the upstream and downstream sites were more similar post-ASR. Most of the primary constituents that make up dissolved solids decreased at both sites after the ASR Phase II facility began operation. Discretely collected total suspended solids concentrations were similar between the upstream and downstream sites before the facility began operating but were about 27 percent smaller at the downstream site after the facility began operating, despite the total suspended solids concentrations in the residuals return line being 15 times larger than the downstream site. Overall habitat scores were indicative of suboptimal conditions upstream and downstream from the ASR Phase II facility throughout the study period. Substrate fouling and sediment deposition mean scores indicated marginal conditions at the upstream and downstream sites during the study period, demonstrating that sediment deposition was evident pre- and post-ASR and no substantial changes in these habitat characteristics were noted after the ASR Phase II facility began discharging. Macroinvertebrate community composition (evaluated using functional feeding, behavioral, and tolerance metrics) generally was similar between sites during the study period. Fewer macroinvertebrate metrics were significant between the upstream and downstream sites post-ASR (6) than pre-ASR (14), which suggests that macroinvertebate communities were more similar after the ASR facility began discharging. Upstream-downstream comparisons in macroinvertebrate aquatic-life-support metrics had no significant differences for the post-ASR time period and neither site was fully supporting for any of the Kansas Department of Health and Environment aquatic-life-support metrics (Macroinvertebrate Biotic Index; Kansas Biotic Index with tolerances for nutrients and oxygen-demanding substances; Ephemeroptera, Plecoptera, and Trichoptera [EPT] richness; and percentage of EPT species). Overall, using macroinvertebrate aquatic life-support criteria from the Kansas Department of Health and Environment, upstream and downstream sites were classified as partially supporting before and after the onset of ASR facility operations. Fish community trophic status and tolerance groups generally were similar among sites during the study period. Fish community Little Arkansas River Basin Index of Biotic Integrity scores at the upstream and downstream sites were indicative of fair-to-good conditions before the facility began operating and decreased to fair conditions after the facility began operating. Groundwater physicochemical changes concurrent with the beginning of recharge operations at the Sedgwick basin were more pronounced in shallow groundwater. No constituent concentrations in the pre-recharge period in comparison to the post-recharge period increased to concentrations exceeding drinking water regulations; however, nitrate decreased significantly from a pre-recharge exceedance of the U.S. Environmental Protection Agency maximum contaminant level to a post recharge nonexceedance. Shallow groundwater chemical concentrations or rates of detection increased after artificial recharge began for the ions potassium, chloride, and fluoride; phosphorus and organic carbon species; trace elements barium, manganese, nickel, arsenate, arsenic, and boron; agricultural pesticides atrazine, metolachlor, metribuzin, and simazine; organic disinfection byproducts bromodichloromethane and trichloromethane; and gross beta levels. Additionally, water temperature, and pH were larger after recharge began; and total solids and slime-forming bacteria concentrations and densities were smaller. Total solids, nitrate, and selenium significantly decreased; and potassium, chloride, nickel, arsenic, fluoride, phosphorus and carbon species, and gross beta levels significantly increased in shallow groundwater after artificial recharge. Results of biological activity reaction tests indicated that water quality microbiology was different before and after artificial recharge began; at times, these differences may lead to changes in dominant bacterial populations that, in turn, may lead to formation and expansion in populations that may cause bioplugging and other unwanted effects. Calcite, iron (II) hydroxide, hydroxyapatite, and similar minerals, had shifts in saturation indices that generally were from undersaturation toward equilibrium and, in some cases, toward oversaturation. These shifts toward neutral saturation indices might suggest reduced weathering of the minerals present in the Equus Beds aquifer. Chemical weathering in the shallow parts of the aquifer may be accelerated because of the increased water temperatures and the system is more vulnerable to clogged pores and mineral dissolution as the equilibrium state is affected by recharge and withdrawal. When oversaturation is indicated for iron minerals, plugging of aquifer materials may happen.

Kansas↗

Inset groundwater-flow models for the Cache and Grand Prairie Critical Groundwater Areas, northeastern Arkansas

The water resources in the Mississippi alluvial plain, located in parts of Missouri, Kentucky, Tennessee, Mississippi, Louisiana, and Arkansas, supports a multibillion-dollar agricultural industry that relies heavily on pumping of groundwater for irrigation of crops and aquaculture. The primary source of groundwater for agricultural-related pumping is the Mississippi River Valley alluvial aquifer, which has declined in storage for decades; secondary groundwater sources include the middle Claiborne aquifer and Wilcox aquifer system. Two areas in northeastern Arkansas that lie within the Mississippi alluvial plain, part of the Cache and Grand Prairie regions, have been designated as Critical Groundwater Areas owing to decades of groundwater declines that resulted from past and current water use. The multidisciplinary Mississippi Alluvial Plain project, led by the U.S. Geological Survey, and funded by their Water Availability and Use Science Program, included objectives to develop numerical groundwater models in focus regions, including the part of the Cache and Grand Prairie regions of northeastern Arkansas. Two inset models were developed using the child model capabilities of MODFLOW 6, the U.S. Geological Survey’s Modular Hydrologic Model simulation software. Both models, called the Cache model and Grand Prairie model, simulated the groundwater system and surface-water/groundwater interactions for the Mississippi River Valley alluvial aquifer and underlying Tertiary-age aquifers and confining units to the Midway confining unit. Each model was spatially discretized into 500-meter x 500-meter orthogonal cells on a grid with 5-meter constant-thickness vertical layers that represented the Mississippi River Valley alluvial aquifer and increasing thickness layers for the aquifers and confining units below the alluvial aquifer. The Cache and Grand Prairie models were calibrated with the PEST++ iterative ensemble smoother Version 5 and employed high dimensional parameterization schemes of 13,740 and 30,436 parameters, respectively. The Cache mean absolute residual for groundwater-level observations within each model domain for the priority well was 1.58 meters. Grand Prairie mean absolute residuals for the alluvial aquifer and middle Claiborne aquifer groundwater-level observations were 2.71 and 10.78 meters, respectively. The groundwater budgets for the Cache and Grand Prairie models were characterized by substantial outflows to irrigation wells, which constituted about 52 and 54 percent of all outflows, with the primary source of water to those wells being releases from unconfined aquifer storage.

Arkansas↗

Potentiometric surfaces and water-level trends in the Cockfield (upper Claiborne) aquifer in southern Arkansas and the Wilcox (lower Wilcox) aquifer of northeastern and southern Arkansas, 2012

The Cockfield aquifer, located in southern Arkansas, is composed of Eocene-age sand beds found near the base of the Cockfield Formation of Claiborne Group. The Wilcox aquifer, located in northeastern and southern Arkansas, is composed of Paleocene-age sand beds found in the middle to lower part of the Wilcox Group. The Cockfield and Wilcox aquifers are primary sources of groundwater. In 2010, withdrawals from the Cockfield aquifer in Arkansas totaled 19.2 million gallons per day (Mgal/d), and withdrawals from the Wilcox aquifer totaled 36.5 Mgal/d. A study was conducted by the U.S. Geological Survey in cooperation with the Arkansas Natural Resources Commission and the Arkansas Geological Survey to measure water levels associated with the Cockfield aquifer and the Wilcox aquifer in northeastern and southern Arkansas. Water levels were measured at 43 wells completed in the Cockfield aquifer and 47 wells completed in the Wilcox aquifer in February and March 2012. Measurements from 2012 are presented as potentiometric-surface maps and in combination with measurements from 2006 as water-level difference maps. Trends in water-level change over time within the Cockfield and Wilcox aquifers were determined using the water-level difference maps and selected well hydrographs. The Cockfield aquifer study area in southern Arkansas is bounded on the east by the Mississippi River and on the west by the area that contains outcrops and subcrops of the Cockfield Formation. The northern boundary of the Cockfield aquifer study area is defined by the area that contains observation wells completed in the Cockfield aquifer and the southern boundary is the Louisiana State line. The Wilcox aquifer study area in northeastern Arkansas is bounded on the east by the Mississippi River and on the north by the Missouri State line. The southern and western boundaries are defined by areas containing observation wells completed in the Wilcox aquifer or by outcrop areas on or near Crowleys Ridge. The Wilcox aquifer study area in southern Arkansas is defined by observation wells completed in the Wilcox aquifer or by areas that contain outcrops of the Wilcox Group, or both. The potentiometric-surface map of the Cockfield aquifer shows the regional direction of groundwater flow was generally toward the east-southeast, except in areas of intense groundwater withdrawals such as southwestern Ashley County, where groundwater flows toward the town of Crossett. The highest water-level altitude measured was 350 feet (ft) above National Geodetic Vertical Datum of 1929 (NGVD 29) in central Columbia County. The lowest water-level altitude measured was 40 ft above NGVD 29 in southeastern Lincoln County. The water-level difference map for the Cockfield aquifer in Arkansas was constructed using 42 water-level measurements made during 2006 and 2012. The difference in water levels for the Cockfield aquifer ranged from 27.4 ft to -10.4 ft. The largest water-level rise was in Calhoun County, and the largest water-level decline was 10.4 ft in Union County. Of the 42 wells, 13 wells had a rise in water level, and the remaining 29 wells had a decline in water level. Hydrographs for 32 wells in the Cockfield aquifer with historical water-level data were evaluated using linear regression to calculate the annual rise or decline for each well. These data were aggregated by county and statistically evaluated for the range, mean, and median of water-level change in each county. Hydrographs for Bradley, Calhoun, Chicot, Columbia, and Union Counties indicated both rising and declining water levels. The mean annual water-level rise or decline for Calhoun County was 0.00 foot per year (ft/yr) or unchanged. The mean annual water-level for Ashley, Bradley, Chicot, Cleveland, Columbia, Lincoln, and Union Counties show declines ranging from -0.02 to -1.10 ft/yr. Two potentiometric-surface maps, one for the southern area and one for the northeastern area, were constructed to show the altitude of the water surface in the Wilcox aquifer. The direction of groundwater flow in the northeastern area was generally towards the south-southwest except for some areas immediately adjacent to the Mississippi River where the flow was more eastward towards the river. The highest water-level altitude was 219 ft in northern Mississippi County, and the lowest water-level altitude was 123 ft near West Memphis in Crittenden County. The direction of groundwater flow in the northern part of the southern area was generally towards the southwest. The direction of groundwater flow in the southern part was in all directions because of two cones of depression and two water-level mounds. The highest water-level altitude measured was 394 ft at the center of a water-level mound in eastern Hot Spring County and a water-level mound in southwestern Hempstead County. The lowest water-level altitude measured was 145 ft at the center of the cone of depression in Clark County. Water-level difference maps for the Wilcox aquifer in Arkansas were constructed using 47 water-level measurements made during 2006 and 2012. The difference in water levels for the Wilcox aquifer in the northeastern area ranged from 22.0 ft to -17.9 ft. The largest rise in water level occurred in Crittenden County, and the largest decline occurred in Lee County. Twenty-one wells had rising water levels, and 10 wells had declining water levels. The difference in water levels for the Wilcox aquifer in the southern area ranged from 18.1 ft to -4.2 ft. The largest rise and the largest decline in water level occurred in Nevada County. Twelve wells had rising water levels, and 4 wells had declining water levels. Linear regression analysis of long-term hydrographs was used to determine the mean annual water-level rise and decline in the Wilcox aquifer in the northeastern and southern areas of Arkansas. In the northeastern area, the mean annual water level declined in all seven counties. The mean annual declines ranged from -0.55 ft/yr in Craighead County to -1.46 ft/yr in St. Francis County. In the southern area, the annual rise and decline calculations for wells with over 20 years of records indicate rising and declining water levels in Clark, Hot Spring, and Nevada Counties. The mean annual water level declined in all counties except Hot Spring County.

Arkansas↗

Where groundwater seeps: Evaluating modeled groundwater discharge patterns with thermal infrared surveys at the river-network scale

Predicting baseflow dynamics, protecting aquatic habitat, and managing legacy contaminants requires explicit characterization and prediction of groundwater discharge patterns throughout river networks. Using handheld thermal infrared (TIR) cameras, we surveyed 47 km of stream length across the Farmington River watershed (1,570 km 2 ; CT and MA, USA), mapping locations of bank and waterline groundwater discharges based on their thermal signature. Using the observed groundwater discharge locations and predicted groundwater discharge rates from 6 variations of a numerical groundwater-flow model (MODFLOW-NWT), we compared 1) predicted groundwater-discharge rates in areas with and without observed groundwater discharge, 2) spatial patterns of observed and predicted groundwater discharge locations, and 3) density of observed groundwater discharge locations with predicted discharge rates. Five of six models reasonably predicted the spatial patterns of discharge locations along the 5th order mainstem, but fewer models predicted groundwater discharge patterns in smaller streams. Our results highlight 1) the feasibility of using TIR observations to evaluate groundwater models, 2) model parameters that influence discharge prediction accuracy (riverbed sediment and bedrock hydraulic conductivity and river-aquifer connections), and 3) current strengths and future opportunities for improved modeling of groundwater-discharge patterns.

Connecticut, Massachusetts↗

Groundwater hydrology and estimation of horizontal groundwater flux from the Rio Grande at selected locations in Albuquerque, New Mexico, 2003-9

The Albuquerque, New Mexico, area has two principal sources of water: groundwater from the Santa Fe Group aquifer system and surface water from the San Juan-Chama Diversion Project. From 1960 to 2002, groundwater withdrawals from the Santa Fe Group aquifer system have caused water levels to decline more than 120 feet in some places within the Albuquerque area, resulting in a great deal of interest in quantifying the river-aquifer interaction associated with the Rio Grande. In 2003, the U.S. Geological Survey in cooperation with the Bureau of Reclamation, the Middle Rio Grande Endangered Species Collaborative Program, and the U.S. Army Corps of Engineers began a detailed characterization of the hydrogeology of the Rio Grande riparian corridor in the Albuquerque, New Mexico, area to provide hydrologic data and enhance the understanding of rates of water leakage from the Rio Grande to the alluvial aquifer, groundwater flow through the aquifer, and discharge of water from the aquifer to the riverside drains. A simple conceptual model of flow indicates that the groundwater table gently slopes from the Rio Grande towards riverside drains and the outer boundaries of the inner valley. Water infiltrating from the Rio Grande initially moves vertically below the river, but, as flow spreads farther into the Rio Grande inner valley alluvial aquifer, flow becomes primarily horizontal. The slope of the water-table surface may be strongly controlled by the riverside drains and influenced by other more distal hydrologic boundary conditions, such as groundwater withdrawals by wells. Results from 35 slug tests performed in the Rio Grande inner valley alluvial aquifer during January and February 2009 indicate that hydraulic-conductivity values ranged from 5 feet per day to 160 feet per day with a median hydraulic-conductivity for all transects of 40 feet per day. Median annual horizontal hydraulic gradients in the Rio Grande inner valley alluvial aquifer ranged from 0.011 to 0.002. Groundwater fluxes through the alluvial aquifer calculated by using median slug-test results (qm slug ) and Darcy's law ranged from about 0.1 feet per day to about 0.7 feet per day. Groundwater fluxes calculated by using the Suzuki-Stallman method (qm heat ) ranged from 0.52 feet per day to 0.23 feet per day. Results from the Darcy's law and Suzuki-Stallman flux calculations were compared to discharge measured in riverside drains on both sides of the river north of the Montaño Bridge on February 26, 2009. Flow in the Corrales Riverside Drain increased by 1.4 cubic feet per second from mile 2 to mile 4, about 12 cubic feet per day per linear foot of drain. Flow in the Albuquerque Riverside Drain increased by 15 cubic feet per second between drain miles 0 and 3, about 82 cubic feet per day per linear foot of drain. The flux of water from the river to the aquifer was calculated to be 2.2 cubic feet per day per linear foot of river by using the median qm slug of 0.09 feet per day at Montaño transects west of the river. The total flux was calculated to be 6.0 cubic feet per day per linear foot of river by using the mean(qm heat of 0.24 feet per day for the Montaño transects west of the river. Assuming the Corrales Riverside Drain intercepted all of this flow, the qm slug or qm heat fluxes account for 18 to 50 percent, respectively, of the increase of flow in the drain. The flux of water from the river to the aquifer was calculated to be 15 cubic feet per day per linear foot of river by using the median qm slug of 0.30 feet per day at the Montaño transects east of the river. The flux of water from the river to the aquifer was calculated to be 17 cubic feet per day per linear foot of river by using the mean flux calculated from the Suzuki-Stallman method for the Montaño East transects of 0.34 feet per day. Assuming the Albuquerque Riverside Drain intercepted all this flow, the qm slug or (qm heat fluxes would only account for 18 to 21 percent, respectively, of the increase in flow in the drain. The comparison of these results with those of previous investigations suggests that calculated flux through the Rio Grande inner valley alluvial aquifer is strongly scale dependent and that the thickness of aquifer through which river water flows may be greater than indicated by the vertical temperature profiles.

Albuquerque, New Mexico↗

Hydrologic conditions in urban Miami-Dade County, Florida, and the effect of groundwater pumpage and increased sea level on canal leakage and regional groundwater flow

The extensive and highly managed surface-water system in southeastern Florida constructed during the 20th Century has allowed for the westward expansion of urban and agricultural activities in Miami-Dade County. In urban areas of the county, the surface-water system is used to (1) control urban flooding, (2) supply recharge to production well fields, and (3) control seawater intrusion. Previous studies in Miami-Dade County have determined that on a local scale, leakage from canals adjacent to well fields can supply a large percentage (46 to 78 percent) of the total groundwater pumpage from production well fields. Canals in the urban areas also receive seepage from the Biscayne aquifer that is derived from a combination of local rainfall and groundwater flow from Water Conservation Area 3 and Everglades National Park, which are west of urban areas of Miami-Dade County. To evaluate the effects of groundwater pumpage on canal leakage and regional groundwater flow, the U.S. Geological Survey (USGS) developed and calibrated a coupled surface-water/groundwater model of the urban areas of Miami-Dade County, Florida. The model was calibrated by using observation data collected from January 1997 through December 2004. The model calibration was verified using observation data collected from January 2005 through December 2010. A 1-year warmup period (January 1996 through December 1996) was added prior to the start of the calibration period to reduce the effects of inaccurate initial conditions on model calibration. The model is designed to simulate surface-water stage and discharge in the managed canal system and dynamic canal leakage to the Biscayne aquifer as well as seepage to the canal from the aquifer. The model was developed using USGS MODFLOW–NWT with the Surface-Water Routing (SWR1) Process to simulate surface-water stage, surface-water discharge, and surface-water/groundwater interaction and the Seawater Intrusion (SWI2) Package to simulate seawater intrusion, respectively. Automated parameter estimation software (PEST) and highly parameterized inversion techniques were used to calibrate the model to observed surface-water stage, surface-water discharge, net surface-water subbasin discharge, and groundwater level data from 1997 through 2004 by modifying hydraulic conductivity, specific storage coefficients, specific yield, evapotranspiration parameters, canal roughness coefficients (Manning’s n values), and canal leakance coefficients. Tikhonov regularization was used to produce parameter distributions that provide an acceptable fit between model outputs and observation data, while simultaneously minimizing deviations from preferred values based on field measurements and expert knowledge. Analytical and simulated water budgets for the period from 1996 through 2010 indicate that most of the water discharging through the salinity control structures is derived from within the urban parts of the study area and that, on average, the canals are draining the Biscayne aquifer. Simulated groundwater discharge from the urban areas to the coast is approximately 7 percent of the total surface-water inflow to Biscayne Bay and is consistent with previous estimates of fresh groundwater discharge to Biscayne Bay. Simulated groundwater budgets indicate that groundwater pumpage in some surface-water basins ranges from 13 to 27 percent of the sum of local sources of groundwater inflow. The largest percentage of groundwater pumpage to local sources of groundwater inflow occurs in the basins that have the highest pumping rates (C–2 and C–100 Basins). The ratio of groundwater pumpage to simulated local sources of groundwater inflow is less than values calculated in previous local-scale studies. The position of the freshwater-seawater interface at the base of the Biscayne aquifer did not change notably during the simulation period (1996–2010), consistent with the similar positions of the interface in 1984, 1995, and 2011 under similar hydrologic and groundwater pumping conditions. Landward movement of the freshwater-seawater interface above the base of the aquifer is more prone to occur during relatively dry years. The model was used to evaluate the effect of increased groundwater pumpage and (or) increased sea level on canal leakage, regional groundwater flow, and the position of the freshwater-seawater interface. Permitted groundwater pumping rates, which generally exceed historical groundwater pumping rates, were used for Miami-Dade County Water and Sewer Department groundwater pumping wells in the base-case future scenario. Base-case future and increased pumping scenario results suggest seawater intrusion may occur at the Miami-Springs well field if the Miami Springs, Hialeah, and Preston well fields are operated using current permitted groundwater pumping rates. Scenario simulations also show that, in general, the canal system limits the adverse effects of proposed groundwater pumpage increases on water-level changes and saltwater intrusion. Proposed increases (up to a 7 percent increase) in groundwater pumpage do not have a notable effect on movement of the freshwater-seawater interface. Increased groundwater pumpage increased lateral groundwater inflow into basins subject to additional groundwater pumpage; however, most (55 percent) of the additional groundwater extracted from pumping wells was supplied by changes in canal seepage and leakage in urban areas of the model. Increased sea level caused increased water-table elevations in urban areas and decreased hydraulic gradients across the system; the largest increases in water-table elevations occurred seaward of the salinity control structures. The extent of flood-prone areas and the percentage of time water-table elevations in flood-prone areas were less than 0.5 foot below land surface increased with increased sea level. Increased sea level also resulted in landward migration of the freshwater-seawater interface; the largest changes in the position of the interface occurred seaward of the salinity control structures except in parts of the model area that were inundated by increased sea level. Decreased water-table gradients reduced groundwater inflow, groundwater outflow, canal exchanges, surface-water inflow, and surface-water outflow through salinity control structures. Results for the scenario that evaluated the combination of increased groundwater pumpage and increased sea level did not differ substantially from the scenario that evaluated increased sea level alone. Groundwater inflow, groundwater outflow, and canal exchanges were reduced in urban areas of the study area as a result of decreased water-table gradients across the system, although reductions were less than those in the increased sea-level scenario. The decline in groundwater levels caused by increased groundwater pumpage was less under the increased sea-level scenario than under the increased groundwater-pumpage scenario. The largest reductions in surface-water outflow from the salinity control structures occurred with increased sea level and increased groundwater pumpage. The model was designed specifically to evaluate the effect of groundwater pumpage on canal leakage at the surface-water-basin scale and thus may not be appropriate for (1) predictions that are dependent on data not included in the calibration process (for example, subdaily simulation of high-intensity events and travel times) and (or) (2) hydrologic conditions that are substantially different from those during the calibration and verification periods. The reliability of the model is limited by the conceptual model of the surface-water and groundwater system, the spatial distribution of physical properties, the scale and discretization of the system, and specified boundary conditions. Some of the model limitations are manifested in model errors. Despite these limitations, however, the model represents the complexities of the interconnected surface-water and groundwater systems that affect how the systems respond to groundwater pumpage, sea-level rise, and other hydrologic stresses. The model also quantifies the relative effects of groundwater pumpage and sea-level rise on the surface-water and groundwater systems.

Florida↗

Simulation of the effects of rainfall and groundwater use on historical lake water levels, groundwater levels, and spring flows in central Florida

The urbanization of central Florida has progressed substantially in recent decades, and the total population in Lake, Orange, Osceola, Polk, and Seminole Counties more than quadrupled from 1960 to 2010. The Floridan aquifer system is the primary source of water for potable, industrial, and agricultural purposes in central Florida. Despite increases in groundwater withdrawals to meet the demand of population growth, recharge derived by infiltration of rainfall in the well-drained karst terrain of central Florida is the largest component of the long-term water balance of the Floridan aquifer system. To complement existing physics-based groundwater flow models, artificial neural networks and other data-mining techniques were used to simulate historical lake water level, groundwater level, and spring flow at sites throughout the area. Historical data were examined using descriptive statistics, cluster analysis, and other exploratory analysis techniques to assess their suitability for more intensive data-mining analysis. Linear trend analyses of meteorological data collected by the National Oceanic and Atmospheric Administration at 21 sites indicate 67 percent of sites exhibited upward trends in air temperature over at least a 45-year period of record, whereas 76 percent exhibited downward trends in rainfall over at least a 95-year period of record. Likewise, linear trend analyses of hydrologic response data, which have varied periods of record ranging in length from 10 to 79 years, indicate that water levels in lakes (307 sites) were about evenly split between upward and downward trends, whereas water levels in 69 percent of wells (out of 455 sites) and flows in 68 percent of springs (out of 19 sites) exhibited downward trends. Total groundwater use in the study area increased from about 250 million gallons per day (Mgal/d) in 1958 to about 590 Mgal/d in 1980 and remained relatively stable from 1981 to 2008, with a minimum of 559 Mgal/d in 1994 and a maximum of 773 Mgal/d in 2000. The change in groundwater-use trend in the early 1980s and the following period of relatively slight trend is attributable to the concomitant effects of increasing public-supply withdrawals and decreasing use of water by the phosphate industry and agriculture. On the basis of available historical data and exploratory analyses, empirical lake water-level, groundwater-level, and spring-flow models were developed for 22 lakes, 23 wells, and 6 springs. Input time series consisting of various frequencies and frequency-band components of daily rainfall (1942 to 2008) and monthly total groundwater use (1957 to 2008) resulted in hybrid signal-decomposition artificial neural network models. The final models explained much of the variability in observed hydrologic data, with 43 of the 51 sites having coefficients of determination exceeding 0.6, and the models matched the magnitude of the observed data reasonably well, such that models for 32 of the 51 sites had root-mean-square errors less than 10 percent of the measured range of the data. The Central Florida Artificial Neural Network Decision Support System was developed to integrate historical databases and the 102 site-specific artificial neural network models, model controls, and model output into a spreadsheet application with a graphical user interface that allows the user to simulate scenarios of interest. Overall, the data-mining analyses indicate that the Floridan aquifer system in central Florida is a highly conductive, dynamic, open system that is strongly influenced by external forcing. The most important external forcing appears to be rainfall, which explains much of the multiyear cyclic variability and long-term downward trends observed in lake water levels, groundwater levels, and spring flows. For most sites, groundwater use explains less of the observed variability in water levels and flows than rainfall. Relative groundwater-use impacts are greater during droughts, however, and long-term trends in water levels and flows were identified that are consistent with historical groundwater-use patterns. The sensitivity of the hydrologic system to rainfall is expected, owing to the well-drained karst terrain and relatively thin confinement of the Floridan aquifer system in much of central Florida. These characteristics facilitate the relatively rapid transmission of infiltrating water from rainfall to the water table and contribute to downward leakage of water to the Floridan aquifer system. The areally distributed nature of rainfall, as opposed to the site-specific nature of groundwater use, and the generally high transmissivity and low storativity properties of the semiconfined Floridan aquifer system contribute to the prevalence of water-level and flow patterns that mimic rainfall patterns. In general, the data-mining analyses demonstrate that the hydrologic system in central Florida is affected by groundwater use differently during wet periods, when little or no system storage is available (high water levels), compared to dry periods, when there is excess system storage (low water levels). Thus, by driving the overall behavior of the system, rainfall indirectly influences the degree to which groundwater use will effect persistent trends in water levels and flows, with groundwater-use impacts more prevalent during periods of low water levels and spring flows caused by low rainfall and less prevalent during periods of high water levels and spring flows caused by high rainfall. Differences in the magnitudes of rainfall and groundwater use during wet and dry periods also are important determinants of hydrologic response. An important implication of the data-mining analyses is that rainfall variability at subannual to multidecadal timescales must be considered in combination with groundwater use to provide robust system-response predictions that enhance sustainable resource management in an open karst aquifer system. The data-driven approach was limited, however, by the confounding effects of correlation between rainfall and groundwater use, the quality and completeness of the historical databases, and the spatial variations in groundwater use. The data-mining analyses indicate that available historical data when used alone do not contain sufficient information to definitively quantify the related individual effects of rainfall and groundwater use on hydrologic response. The knowledge gained from data-driven modeling and the results from physics-based modeling, when compared and used in combination, can yield a more comprehensive assessment and a more robust understanding of the hydrologic system than either of the approaches used separately.

Florida↗

Environmental tracer evidence for connection between shallow and bedrock aquifers and high intrinsic susceptibility to contamination of the conterminous U.S. glacial aquifer

Covering a large portion of the northern conterminous United States (1.87 x 106 km2), the glacial aquifer serves as the primary water supply for 39 million public and domestic water users. Mean groundwater age, groundwater age distribution, and susceptibility to land surface contamination, using a new metric (Susceptibility Index; SI) based on the full age distribution and less prone to bias than estimated mean age, is reported for 168 public and domestic wells across the aquifer. Comparison of groundwater age metrics between well networks of varying spatial scale suggest an extensive sample network of equally spaced, long screened interval wells can be used to characterize aquifer wide groundwater age. Estimated mean age ranges from 1 to 50,000 years and, according to the composite age distribution, approximately 63 percent of all sampled water recharged after 1950 (i.e., modern) and 18 percent of the sampled water was recharged greater than 10,000 years ago. The later finding strongly suggests a connection between the glacial aquifer and underlying bedrock aquifers. Statistical analysis of glacial aquifer hydrogeology and age metrics show groundwater ages are young (less than few 100 years) and more susceptible to land surface contamination (larger SI) in unconfined and shallow portions of the aquifer. Old groundwater (greater than 1000 years) is more often associated with thicker sequences of fine grain sediments and/or shallow bedrock. Calculated SI is shown to be more strongly related to the number of land surface contaminants detected than mean age or fraction modern. Statistical analysis of SI and hydrogeology indicates SI is largely dictated by well depth and confinement. This study demonstrates how sample network design can be used to characterize groundwater age of large aquifers with a limited number of samples and how interpretation of environmental tracers can be used to improve conceptual models of groundwater aquifers and identify groundwater susceptible to contamination.

Glacial aquifer↗

Water-level altitudes 2016 and water-level changes in the Chicot, Evangeline, and Jasper aquifers and compaction 1973–2015 in the Chicot and Evangeline aquifers, Houston-Galveston region, Texas

Most of the land-surface subsidence in the Houston-Galveston region, Texas, has occurred as a direct result of groundwater withdrawals for municipal supply, commercial and industrial use, and irrigation that depressured and dewatered the Chicot and Evangeline aquifers, thereby causing compaction of the aquifer sediments, mostly in the fine-grained silt and clay layers. This report, prepared by the U.S. Geological Survey 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, is one in an annual series of reports depicting water-level altitudes and water-level changes in the Chicot, Evangeline, and Jasper aquifers and measured cumulative compaction of subsurface sediments in the Chicot and Evangeline aquifers in the Houston-Galveston region. The report contains regional-scale maps depicting approximate 2016 water-level altitudes (represented by measurements made during December 2015–March 2016) for the Chicot, Evangeline, and Jasper aquifers; maps depicting 1-year (2015–16) water-level changes for each aquifer; maps depicting approximate contoured 5-year (2011–16) water-level changes for each aquifer; maps depicting approximate contoured long-term (1990–2016 and 1977–2016) water-level changes for the Chicot and Evangeline aquifers; a map depicting approximate contoured long-term (2000–16) water-level changes for the Jasper aquifer; a map depicting locations of borehole-extensometer sites; and graphs depicting measured long-term cumulative compaction of subsurface sediments at the extensometers during 1973–2015. Tables listing the water-level data used to construct each water-level map for each aquifer and the measured long-term cumulative compaction data for each extensometer site are included. Graphs depicting water-level measurement data also are included; these graphs can be used to approximate changes in effective stress caused by changes in groundwater withdrawal from the Chicot and Evangeline aquifers. In 2016, water-level-altitude contours for the Chicot aquifer ranged from 200 feet (ft) below the vertical datum (North American Vertical Datum of 1988; hereinafter, datum) in a localized area in northwestern Harris County to 200 ft above datum in west-central Montgomery County. Water-level changes during 2015–16 in the Chicot aquifer ranged from a 39-ft decline to a 26-ft rise. Contoured 5-year and long-term changes in water-level altitudes of the Chicot aquifer ranged from a 30-ft decline to a 20-ft rise (2011–16), from a 140-ft decline to a 160-ft rise (1990–2016), and from a 120-ft decline to a 200-ft rise (1977–2016). In 2016, water-level-altitude contours for the Evangeline aquifer ranged from 250 ft below datum in three separate areas in south-central Montgomery County and extending into north-central Harris County, in west-central Harris County, and in southwestern Harris County to 200 ft above datum in southeastern Grimes and northwestern Montgomery Counties. Water-level changes during 2015–16 in the Evangeline aquifer ranged from a 65-ft decline to a 61-ft rise. Contoured 5-year and long-term changes in water-level altitudes of the Evangeline aquifer ranged from a 60-ft decline to a 40-ft rise (2011–16), from a 160-ft decline to a 160-ft rise (1990–2016), and from a 320-ft decline to a 240-ft rise (1977–2016). In 2016, water-level-altitude contours for the Jasper aquifer ranged from 200 ft below datum in south-central Montgomery County extending into north-central Harris County to 250 ft above datum in northwestern Montgomery County and extending into eastern Grimes County and southwestern Walker County. Water-level changes during 2015–16 in the Jasper aquifer ranged from a 38-ft decline to a 51-ft rise. Contoured 5-year and long-term changes in water-level altitudes of the Jasper aquifer ranged from a 60-ft decline to a 40-ft rise (2011–16) and from a 220-ft decline to a 20-ft decline (2000–16). Compaction of subsurface sediments (mostly in the fine-grained silt and clay layers) in the Chicot and Evangeline aquifers was recorded continuously by using 13 extensometers at 11 sites that were either activated or installed between 1973 and 1980. During the period of record beginning in 1973 (or later depending on activation or installation date) and ending in December 2015, measured cumulative compaction at the 13 extensometers ranged from 0.095 ft at the Texas City-Moses Lake extensometer to 3.666 ft at the Addicks extensometer. From January through December 2015, the Northeast, Southwest, Addicks, Johnson Space Center, and Clear Lake (deep) extensometers recorded net decreases in land-surface elevation, but the Lake Houston, East End, Texas City-Moses Lake, Baytown C–1 (shallow), Baytown C–2 (deep), Seabrook, Clear Lake (shallow), and Pasadena extensometers recorded net increases in land-surface elevation. For the 11 extensometer sites during the selected years 1988, 1998, 2008, 2012, and 2015, the smallest effective stress (20.12 pounds per square inch [psi]) was estimated at the Texas City-Moses Lake extensometer site and was produced by a measured water level of 46.42 ft below land-surface datum (blsd) in January 2008. The corresponding net compaction during 2007 at this site was 0.001 ft. The largest effective stress (174.86 psi) was estimated at the Addicks extensometer site and was produced by a measured water level of 403.38 ft blsd in January 1998. The corresponding net compaction at the Addicks site was 0.067 ft in 1997. The 2011 drought caused water-level declines in the aquifers that were documented by the water-level-measurement data collected in January 2012. During the 2011 drought, the 13 extensometers recorded varying amounts of compaction that ranged from a net compaction value of 0.002 ft recorded by the Texas City-Moses Lake extensometer to a net compaction value of 0.192 ft recorded by the Pasadena extensometer. Water-level data for 1988, 1998, 2008, 2012, and 2015 and the corresponding net compaction values recorded by the extensometers for 1987, 1997, 2007, 2011, and 2014 were used to illustrate the cause and effect relations between changes in water level caused by groundwater withdrawals and resulting changes in effective stress. Changes in effective stress are related to changes in land-surface elevations caused by compaction of the fine-grained sediments composing the Chicot and Evangeline aquifers. The rate of compaction varies from site to site because of differences in rates of groundwater withdrawal in the areas adjacent to each extensometer site; differences among sites in the ratios of sand, silt, and clay and their corresponding compressibilities; and previously established preconsolidation heads. It is not appropriate, therefore, to extrapolate or infer a rate of compaction for an adjacent area on the basis of the rate of compaction recorded by proximal extensometers.

Texas↗

Estimation of groundwater use for a groundwater-flow model of the Lake Michigan Basin and adjacent areas, 1864-2005

The U.S. Geological Survey, at the request of Congress, is assessing the availability and use of the Nation's water resources to help characterize how much water is available now, how water availability is changing, and how much water can be expected to be available in the future. The Great Lakes Basin Pilot project of the U.S. Geological Survey national assessment of water availability and use focused on the Great Lakes Basin and included detailed studies of the processes governing water availability in the Great Lakes Basin. One of these studies included the development of a groundwater-flow model of the Lake Michigan Basin. This report describes the compilation and estimation of the groundwater withdrawals in those areas in Wisconsin, Michigan, Indiana, and Illinois that were needed for the Lake Michigan Basin study groundwater-flow model. These data were aggregated for 12 model time intervals spanning 1864 to 2005 and were summarized by model area, model subregion, category of water use, aquifer system, aquifer type, and hydrogeologic unit model layer. The types and availability of information on groundwater withdrawals vary considerably among states because water-use programs often differ in the types of data collected and in the methods and frequency of data collection. As a consequence, the methods used to estimate and verify the data also vary. Additionally, because of the different sources of data and different terminologies applied for the purposes of this report, the water-use data published in this report may differ from water-use data presented in other reports. These data represent only a partial estimate of groundwater use in each state because estimates were compiled only for areas in Wisconsin, Michigan, Indiana, and Illinois within the Lake Michigan Basin model area. Groundwater-withdrawal data were compiled for both nearfield and farfield model areas in Wisconsin and Illinois, whereas these data were compiled primarily for the nearfield model area in Michigan and Indiana. Overall water use for the selected areas in Wisconsin, Michigan, Indiana, and Illinois was less during early time intervals than during more recent intervals, with large increases beginning around the 1960s. Total estimated groundwater withdrawals for model input range from 18.01 million gallons per day (Mgal/d) for interval 1 (1864-1900) to 1,280.25 Mgal/d for interval 12 (2001-5). Withdrawals for the public-supply category make up the majority of the withdrawals in each of the four states. In Wisconsin and Michigan, the second largest withdrawals are for the irrigation category; in Indiana and Illinois, industrial withdrawals account for the second largest withdrawal amounts. The smallest withdrawals are for miscellaneous uses in Wisconsin and irrigation uses in Indiana and Illinois. Estimated groundwater withdrawals in the Southern Lower Peninsula of Michigan, Northeastern Illinois, and the farfield model area are generally larger than in the other model subregions. Withdrawals in Michigan and Indiana are predominantly from the Quaternary aquifer system, whereas withdrawals in Illinois are predominantly from the Cambrian-Ordovician aquifer systems. Withdrawals in Wisconsin are about equal from the Quaternary and Cambrian-Ordovician aquifer systems. Estimated groundwater withdrawals in Michigan and Indiana are predominantly from the unconfined unconsolidated aquifer type. Withdrawals in Illinois are largely from the deep confined bedrock aquifer type, although they decreased considerably in more recent time intervals. Wisconsin withdrawals are about equal from unconfined unconsolidated and deep confined bedrock aquifer types. Groundwater-withdrawal estimates in Wisconsin were compiled for the 47 easternmost counties within the boundary of the Lake Michigan Basin model, of which 32 counties, though not entirely contained, are at least partly within the Lake Michigan Basin. Overall, 6,457 withdrawal locations were estima

Scientific Investigations Report↗

Geospatial compilation of historical water-level changes in the Chicot and Evangeline aquifers 1977-2013 and Jasper aquifer 2000-13, 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 an annual series of reports that depict water-level changes in the Chicot, Evangeline, and Jasper aquifers of the Gulf Coast aquifer system in the Houston-Galveston region, Texas, from 1977 to 2013. Changes are determined from water-level measurements between December and March of each year from groundwater wells screened in one of the three aquifers. Existing published maps and unpublished geographic information system (GIS) datasets were compiled into a comprehensive geodatabase of all water-level-change maps produced as part of this multiagency effort. Annual water-level-change maps were georeferenced and digitized where existing GIS data were unavailable (1979–99). Existing GIS data available for 2000–13 were included in the geodatabase. The compilation contains 121 datasets showing water-level changes for each primary aquifer of the Gulf Coast aquifer system: 56 for the Chicot aquifer (1977; 1979–2013 and 1990; 1993–2013), 56 for the Evangeline aquifer (1977; 1979–2013 and 1990; 1993–2013), and 9 for the Jasper aquifer (2000; 2005–13).

Texas↗