USGS Science⌕ Search

SEARCH · USGS Science

Results for “Groundwater”

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 217 records · Page 12Linked to original sources
USGS2010

Groundwater-quality data in the South Coast Range-Coastal study unit, 2008: Results from the California GAMA Program

Groundwater quality in the approximately 766-square-mile South Coast Range–Coastal (SCRC) study unit was investigated from May to December 2008, as part of the Priority Basins Project of the Groundwater Ambient Monitoring and Assessment (GAMA) Program. The GAMA Priority Basins Project was developed in response to legislative mandates (Supplemental Report of the 1999 Budget Act 1999-00 Fiscal Year; and, the Groundwater Quality Monitoring Act of 2001 [Sections 10780-10782.3 of the California Water Code, Assembly Bill 599]) to assess and monitor the quality of groundwater in California, and is being conducted by the U.S. Geological Survey (USGS) in cooperation with the California State Water Resources Control Board (SWRCB). The SCRC study unit was the 25th study unit to be sampled as part of the GAMA Priority Basins Project. The SCRC study unit was designed to provide a spatially unbiased assessment of untreated groundwater quality in the primary aquifer systems and to facilitate statistically consistent comparisons of untreated groundwater quality throughout California. The primary aquifer systems (hereinafter referred to as primary aquifers) were defined as that part of the aquifer corresponding to the perforation interval of wells listed in the California Department of Public Health (CDPH) database for the SCRC study unit. The quality of groundwater in shallow or deep water-bearing zones may differ from the quality of groundwater in the primary aquifers; shallow groundwater may be more vulnerable to surficial contamination. In the SCRC study unit, groundwater samples were collected from 70 wells in two study areas (Basins and Uplands) in Santa Barbara and San Luis Obispo Counties. Fifty-five of the wells were selected using a spatially distributed, randomized grid-based method to provide statistical representation of the study unit (grid wells), and 15 wells were selected to aid in evaluation of specific water-quality issues (understanding wells). In addition to the 70 wells sampled, 3 surface-water samples were collected in streams near 2 of the sampled wells in order to better comprehend the interaction between groundwater and surface water in the area. The groundwater samples were analyzed for organic constituents (volatile organic compounds [VOC], pesticides and pesticide degradates, polar pesticides and metabolites, and pharmaceutical compounds), constituents of special interest (perchlorate, N-nitrosodimethylamine [NDMA], and 1,2,3-TCP), naturally occurring inorganic constituents (trace elements, nutrients, dissolved organic carbon [DOC], major and minor ions, silica, total dissolved solids [TDS], and alkalinity), and radioactive constituents (gross alpha and gross beta radioactivity). Naturally occurring isotopes (stable isotopes of hydrogen and oxygen in water, stable isotopes of nitrogen and oxygen in dissolved nitrate, stable isotopes of sulfur in dissolved sulfate, stable isotopes of carbon in dissolved inorganic carbon, activities of tritium, and carbon-14 abundance), and dissolved gases (including noble gases) also were measured to help identify the sources and ages of the sampled groundwater. In total, 298 constituents and field water-quality indicators were investigated. Three types of quality-control samples (blanks, replicates, and matrix-spikes) were collected at approximately 3 to 12 percent of the wells in the SCRC study unit, and the results for these samples were used to evaluate the quality of the data for the groundwater samples. Field blanks rarely contained detectable concentrations of any constituent, suggesting that contamination from sample collection procedures was not a significant source of bias in the data for the groundwater samples. Differences between replicate samples generally were less than 10 percent relative and/or standard deviation, indicating acceptable analytical reproducibility. Matrix-spike recoveries were within the acceptable range (70 to 130 percent) for approximately 84 percent of the compounds. This study did not attempt to evaluate the quality of drinking water delivered to consumers; after withdrawal from the ground, untreated groundwater typically is treated, disinfected, and/or blended with other waters to maintain water quality. Regulatory thresholds apply to water that is served to the consumer, not to untreated groundwater. However, to provide some context for the results, concentrations of constituents measured in the untreated groundwater were compared with regulatory and non-regulatory health-based thresholds established by the U.S. Environmental Protection Agency (USEPA) and CDPH, and to non-regulatory thresholds established for aesthetic concerns by CDPH. Comparisons between data collected for this study and thresholds for drinking water are for illustrative purposes only and are not indicative of compliance or noncompliance with those thresholds. Most organic and inorganic constituents that were detected in groundwater samples from the 55 grid wells in the SCRC study unit were detected at concentrations less than drinking-water thresholds. In addition, all detections of organic constituents in SCRC grid well samples were less than health-based thresholds. In total, VOCs were detected in 33 percent of the 55 grid wells sampled and pesticides and pesticide degradates were detected in 27 percent of grid wells sampled in the SCRC study unit. In the Basins study area, VOCs and pesticides and pesticide degradates were detected in approximately 33 percent of the 39 grid wells. In the Uplands study area, VOCs were detected in approximately 31 percent and pesticides and pesticide degradates were detected in approximately 13 percent of the 16 grid wells. Trace elements and minor ions were sampled for at 32 grid wells and nutrients at 33 grid wells in the SCRC study unit, and most detections were less than health-based thresholds. Exceptions in the Basins study area include one detection of arsenic greater than the USEPA maximum contaminant level (MCL-US) of 10 µg/L and three detections of nitrite plus nitrate, as nitrogen (NO2-+NO3-) greater than the MCL-US of 10 mg/L. Exceptions in the Uplands study area include two detections of arsenic greater than the MCL-US and eight detections of molybdenum greater than the USEPA lifetime health advisory level (HAL-US) of 40 µg/L. All detections of major and minor ions and gross alpha and gross beta radioactivity from the SCRC grid wells were less than health-based thresholds. Results for trace elements, major ions, and TDS with non-enforceable thresholds set for aesthetic concerns from 16 Basins study area grid wells showed that iron concentrations greater than the CDPH secondary maximum contaminant level (SMCL-CA) of 300 µg/L were detected in grid wells. Manganese concentrations greater than the SMCL-CA of 50 µg/L were detected in six grid wells. Chloride concentrations greater than the recommended SMCL-CA threshold of 250 mg/L were detected in one grid well. Sulfate concentrations greater than the recommended SMCL-CA threshold of 250 mg/L were measured in 12 grid wells and 3 of these wells also were greater than the upper SMCL-CA threshold of 500 mg/L. TDS concentrations greater than the SMCL-CA recommended threshold of 500 mg/L were measured in 14 of the 16 Basins study area grid wells and concentrations in 5 of these wells also were greater than the SMCL-CA upper threshold of 1,000 mg/L. In the Uplands study area, iron concentrations greater than the SMCL-CA were detected in 2 of 16 grid wells and manganese concentrations greater than the SMCL-CA were detected in 3 grid wells. TDS and sulfate concentrations greater than the recommended SMCL-CA thresholds were detected in 11 and 2 grid wells, respectively, but none of these concentrations were greater than the SMCL-CA upper thresholds.

California↗
USGS2019

Groundwater availability in the Ozark Plateaus aquifer system

Executive Summary The study described in this report, initiated by the U.S. Geological Survey in 2014, was designed to evaluate fresh groundwater resources within the Ozark Plateaus, central United States, as an area within a broader national assessment of groundwater availability. The goals of the Ozark study were to evaluate historical effects of human activities on water levels and groundwater availability, quantify groundwater resources now and under probable future pumping and climate conditions, and evaluate existing monitoring networks for their value in making better predictions of future groundwater resources. Previous studies include simulation of local-scale groundwater flow under varying temporal scales, or simulation of the regional system under steady-state conditions. While these studies are useful, particularly for the problem for which they were designed, there is a need to look at the larger regional system under transient conditions to fully evaluate the water resource over time. This study focused on multiple spatial and temporal scales to examine changes in groundwater pumping, storage, and water-level declines. The regional scale provides a broad view of the sources and demands on the system with time. The study area covers approximately 68,000 square miles in the central United States in parts of Missouri, Arkansas, Kansas, and Oklahoma and encompasses the Ozark Plateaus Physiographic Province (Ozark Plateaus), including the Salem Plateau, Springfield Plateau, and Boston Mountains. Groundwater is withdrawn from the Ozark Plateaus aquifer system (Ozark system) for public supply and for domestic, agriculture (including irrigation and aquaculture), livestock, and non-agricultural use (including industrial, thermoelectric power generation, mining, and commercial). The Ozark system provides an important drinking-water supply for people living in the Ozark Plateaus because public supply and domestic use combined constitute the largest groundwater use. Precipitation is the ultimate source of freshwater to the Ozark system; most rainfall occurs during April, May, and June, and precipitation increases generally from north to south across the study area. Groundwater use currently accounts for only 10 percent of the total water use in the areas overlying the Ozark system, but provides a critical drinking-water resource because public supply and domestic groundwater withdrawals are largely from groundwater resources. The 380 million gallons per day of groundwater withdrawn from the Ozark system in 2010 accounts for approximately 2 percent of recharge. Although groundwater use represents a small component of the hydrologic budget, because of low storage in aquifer units, cones of depression with steep water-level gradients can develop quickly around pumping centers. The amount of water entering and leaving the aquifer system from 1900 to about 1965 was relatively constant at a rate of about 13 billion gallons per day (Bgal/d). Much of this inflow of water is discharged through streams in the system to balance the hydrologic budget. Changes in storage over time (from outflows to inflows) reflect the large variability in recharge: if recharge decreases, water levels will decrease, resulting in less groundwater discharge to streams and more water released from aquifer storage. Conversely, when recharge increases, water levels increase, more groundwater discharges to streams, and aquifer storage is replenished. Although pumping generally increased from 1900 to 2016, it does not appear to correlate with the change in storage over the same time period. Regionally, simulated change in groundwater storage corresponds with changes in recharge, more so than with increases in pumping. Average recharge was 11.6 Bgal/d for the period 1900 to 2016. Recharge was generally above average from predevelopment to 1965, followed by a period of below-average recharge from 1965 to about 1980. Recharge remained consistently above average from 1980 to about 1988, after which there was a period of average or below-average recharge, reflected by a decline through the mid-2000s. The implications and potential effects of increased pumping and long-term climate change on the Ozark Plateaus hydrologic system and groundwater availability are a concern for communities and resource managers in the area. Pumping varies from year to year, but is generally expected to moderately increase with population, industrial, and agricultural needs. Most climate models predict warmer minimum and maximum air temperatures by midcentury in the Ozark Plateaus area, especially from midspring through early fall. Three scenarios were developed to simulate possible future conditions from 2016 to 2060 and assess the potential effects on the hydrologic system and availability of water resources. For each scenario, changes in water levels and hydrologic budget components were evaluated from predevelopment (1900) to present (2016) and 45 years into the future (2060). The baseline scenario represents an extension of the average (1996 to 2016) seasonal pumping and recharge values. The pumping scenario is an extension of the average (1996 to 2016) seasonal recharge values with increases in pumping following the historical trend for the period 2016–2060 of up to 120 percent of the 1996 to 2016 average seasonal pumping values. The general circulation model (GCM) scenario is an extension of the average (1996 to 2016) seasonal pumping values and variable recharge based on seasonal averages of soil water storage from a water-balance model using temperature and precipitation from multiple GCMs. The general patterns of water-level decline are similar for each scenario. The areas of water-level decline in southwest Missouri and northeast Oklahoma are only marginally different by 2060 from those of 2009. In one area south of Springfield, Mo., water-level declines are less in the baseline and GCM scenarios than in 2009. This may be the result of a transition from groundwater use to surface-water supplies for a larger percentage of the demand in the area. For all three scenarios, forecasted pumping, recharge, and aquifer properties play an important role in determining the uncertainty of water-level forecasts at 94 real-time observation wells. Simulated aquifer properties in the productive middle and lower Ozark aquifers and the St. Francois confining unit of the Ozark system contribute most to predictive uncertainty in water levels at approximately 35 percent of the real-time observation wells. Out of the 94 real-time observation wells, 82 are developed in the lower Ozark aquifer.

Arkansas, Kansas, Missouri, Oklahoma↗
USGS2012

Assessing potential effects of changes in water use with a numerical groundwater-flow model of Carson Valley, Douglas County, Nevada, and Alpine County, California

Rapid growth and development within Carson Valley in Douglas County, Nevada, and Alpine County, California, has caused concern over the continued availability of groundwater, and whether the increased municipal demand could either impact the availability of water or result in decreased flow in the Carson River. Annual pumpage of groundwater has increased from less than 10,000 acre feet per year (acre-ft/yr) in the 1970s to about 31,000 acre-ft/yr in 2004, with most of the water used in agriculture. Municipal use of groundwater totaled about 10,000 acre-feet in 2000. In comparison, average streamflow entering the valley from 1940 to 2006 was 344,100 acre-ft/yr, while average flow exiting the valley was 297,400 acre-ft/yr. Carson Valley is underlain by semi-consolidated Tertiary sediments that are exposed on the eastern side and dip westward. Quaternary fluvial and alluvial deposits overlie the Tertiary sediments in the center and western side of the valley. The hydrology of Carson Valley is dominated by the Carson River, which supplies irrigation water for about 39,000 acres of farmland and maintains the water table less than 5 feet (ft) beneath much of the valley floor. Perennial and ephemeral watersheds drain the Carson Range and the Pine Nut Mountains, and mountain-front recharge to the groundwater system from these watersheds is estimated to average 36,000 acre-ft/yr. Groundwater in Carson Valley flows toward the Carson River and north toward the outlet of the Carson Valley. An upward hydraulic gradient exists over much of the valley, and artesian wells flow at land surface in some areas. Water levels declined as much as 15 ft since 1980 in some areas on the eastern side of the valley. Median estimated transmissivities of Quaternary alluvial-fan and fluvial sediments, and Tertiary sediments are 316; 3,120; and 110 feet squared per day (ft 2 /d), respectively, with larger transmissivity values in the central part of the valley and smaller values near the valley margins. A groundwater-flow model of Quaternary and Tertiary sediments in Carson Valley was developed using MODFLOW and calibrated to simulate historical conditions from water years 1971 through 2005. The 35-year transient simulation represented quarterly changes in precipitation, streamflow, pumping and irrigation. Inflows to the groundwater system simulated in the model include mountain-front recharge from watersheds in the Carson Range and Pine Nut Mountains, valley recharge from precipitation and land application of wastewater, agricultural recharge from irrigation, and septic-tank discharge. Outflows from the groundwater system simulated in the model include evapotranspiration from the water table and groundwater withdrawals for municipal, domestic, irrigation and other water supplies. The exchange of water between groundwater, the Carson River, and the irrigation system was represented with a version of the Streamflow Routing (SFR) package that was modified to apply diversions from the irrigation network to irrigated areas as recharge. The groundwater-flow model was calibrated through nonlinear regression with UCODE to measured water levels and streamflow to estimate values of hydraulic conductivity, recharge and streambed hydraulic-conductivity that were represented by 18 optimized parameters. The aquifer system was simulated as confined to facilitate numerical convergence, and the hydraulic conductivity of the top active model layers that intersect the water table was multiplied by a factor to account for partial saturation. Storage values representative of specific yield were specified in parts of model layers where unconfined conditions are assumed to occur. The median transmissivity ( T ) values (11,000 and 800 ft 2 /d for the fluvial and alluvial-fan sediments, respectively) are both within the third quartile of T values estimated from specific-capacity data, but T values for Tertiary sediments are larger than the third quartile estimated from specific-capacity data. The estimated vertical anisotropy for the Quaternary fluvial sediments (9,000) is comparable to the value estimated for a previous model of Carson Valley. The estimated total volume of mountain-front recharge is equivalent to a previous estimate from the Precipitation-Runoff Modeling System (PRMS) watershed models, but less recharge is estimated for the Carson Range and more recharge is estimated for the Pine Nut Mountains than the previous estimate. Simulated flow paths indicate that groundwater flows faster through the center of Carson Valley and slower through the lower hydraulic-conductivity Tertiary sediments to the east. Shallow flow in the center of the valley is towards drainage channels, but deeper flow is generally directed toward the basin outlet to the north. The aquifer system is in a dynamic equilibrium with large inflows from storage in dry years and large outflows to storage in wet years. Pumping has historically been less than 10 percent of outflows from the groundwater system, and agricultural recharge has been less than 10 percent of inflows to the groundwater system. Three principal sources of uncertainty that affect model results are: (1) the hydraulic characteristics of the Tertiary sediments on the eastern side of the basin, (2) the composition of sediments beneath the alluvial fans and (3) the extent of the confining unit represented within fluvial sediments in the center of the basin. The groundwater-flow model was used in five 55-year predictive simulations to evaluate the long-term effects of different water-use scenarios on water-budget components, groundwater levels, and streamflow in the Carson River. The predictive simulations represented water years 2006 through 2060 using quarterly stress periods with boundary conditions that varied cyclically to represent the transition from wet to dry conditions observed from water years 1995 through 2004. The five scenarios included a base scenario with 2005 pumping rates held constant throughout the simulation period and four other scenarios using: (1) pumping rates increased by 70 percent, including an additional 1,340 domestic wells, (2A) pumping rates more than doubled with municipal pumping increased by a factor of four over the base scenario, (2B) pumping rates of 2A with 2,040 fewer domestic wells, and (3) pumping rates of 2A with 3,700 acres removed from irrigation. The 55-year predictive simulations indicate that increasing groundwater withdrawals under the scenarios considered would result in as much as 40 ft and 60 ft of water-table decline on the west and east sides of Carson Valley, respectively. The water table in the central part of the valley would remain essentially unchanged, but water-level declines of as much as 30 ft are predicted for the deeper, confined aquifer. The increased withdrawals would reduce the volume of groundwater storage and decrease the mean downstream flow in the Carson River by as much as 16,500 acre-ft/yr. If, in addition, 3,700 acres were removed from irrigation, the reduction in mean downstream flow in the Carson River would be only 6,500 acre-ft/yr. The actual amount of flow reduction is uncertain because of potential changes in irrigation practices that may not be accounted for in the model. The projections of the predictive simulations are sensitive to rates of mountain-front recharge specified for the Carson Range and the Pine Nut Mountains. The model provides a tool that can be used to aid water managers and planners in making informed decisions. A prudent management approach would include continued monitoring of water levels on both the east and west sides of Carson Valley to either verify the predictions of the groundwater-flow model or to provide additional data for recalibration of the model if the predictions prove inaccurate.

California;Nevada↗
USGS2013

Status and understanding of groundwater quality in the South Coast Range-Coastal study unit, 2008: California GAMA Priority Basin Project

Groundwater quality in the South Coast Range–Coastal (SCRC) study unit was investigated from May through November 2008 as part of the Priority Basin Project of the Groundwater Ambient Monitoring and Assessment (GAMA) Program. The study unit is located in the Southern Coast Range hydrologic province and includes parts of Santa Barbara and San Luis Obispo Counties. The GAMA Priority Basin Project is conducted by the U.S. Geological Survey (USGS) in collaboration with the California State Water Resources Control Board and the Lawrence Livermore National Laboratory. The GAMA Priority Basin Project was designed to provide a statistically unbiased, spatially distributed assessment of untreated groundwater quality within the primary aquifer system. The primary aquifer system is defined as that part of the aquifer corresponding to the perforation interval of wells listed in the California Department of Public Health (CDPH) database for the SCRC study unit. The assessments for the SCRC study unit were based on water-quality and ancillary data collected in 2008 by the USGS from 55 wells on a spatially distributed grid, and water-quality data from the CDPH database. Two types of assessments were made: (1) status, assessment of the current quality of the groundwater resource, and (2) understanding, identification of the natural and human factors affecting groundwater quality. Water-quality and ancillary data were collected from an additional 15 wells for the understanding assessment. The assessments characterize untreated groundwater quality, not the quality of treated drinking water delivered to consumers by water purveyors. The first component of this study, the status assessment of groundwater quality, used data from samples analyzed for anthropogenic constituents such as volatile organic compounds (VOCs) and pesticides, as well as naturally occurring inorganic constituents such as major ions and trace elements. Although the status assessment applies to untreated groundwater, Federal and California regulatory and non-regulatory water-quality benchmarks that apply to drinking water are used to provide context for the results. Relative-concentrations (sample concentration divided by benchmark concentration) were used for evaluating groundwater. A relative-concentration greater than (>) 1.0 indicates a concentration greater than the benchmark and is classified as high. Inorganic constituents are classified as moderate if relative-concentrations are >0.5 and less than or equal to (≤) 1.0, or low if relative-concentrations are ≤0.5. For organic constituents, the boundary between moderate and low relative-concentrations was set at 0.1. Aquifer-scale proportion was used in the status assessment as the primary metric for evaluating regional-scale groundwater quality. High aquifer-scale proportion is defined as the areal percentage of the primary aquifer system with a high relative-concentration for a particular constituent or class of constituents. Moderate and low aquifer-scale proportions were defined as the areal percentage of the primary aquifer system with moderate and low relative-concentrations, respectively. Two statistical approaches—grid-based and spatially weighted—were used to evaluate aquifer-scale proportions for individual constituents and classes of constituents. Grid-based and spatially weighted estimates were comparable for the study (within 90 percent confidence intervals). For inorganic constituents with human-health benchmarks, relative-concentrations were high for at least one constituent for 33 percent of the primary aquifer system in the SCRC study unit. Arsenic, molybdenum, and nitrate were the primary inorganic constituents with human-health benchmarks that were detected at high relative-concentrations. Inorganic constituents with aesthetic benchmarks, referred to as secondary maximum contaminant levels (SMCLs), had high relative-concentrations for 35 percent of the primary aquifer system. Iron, manganese, total dissolved solids (TDS), and sulfate were the inorganic constituents with SMCLs detected at high relative-concentrations. In contrast to inorganic constituents, organic constituents with human-health benchmarks were not detected at high relative-concentrations in the primary aquifer system in the SCRC study unit. Of the 205 organic constituents analyzed, 21 were detected—13 with human-health benchmarks. Perchloroethene (PCE) was the only VOC detected at moderate relative-concentrations. PCE, dichlorodifluoromethane (CFC-12), and chloroform were detected in more than 10 percent of the primary aquifer system. Of the two special-interest constituents, one was detected; perchlorate, which has a human-health benchmark, was detected at moderate relative-concentrations in 29 percent of the primary aquifer system and had a detection frequency of 60 percent in the SCRC study unit. The second component of this study, the understanding assessment, identified the natural and human factors that may have affected groundwater quality in the SCRC study unit by evaluating statistical correlations between water-quality constituents and potential explanatory factors. The potential explanatory factors evaluated were land use, septic tank density, well depth and depth to top-of-perforations, groundwater age, density and distance to the nearest formerly leaking underground fuel tank (LUFT), pH, and dissolved oxygen (DO) concentration. Results of the statistical evaluations were used to explain the occurrence and distribution of constituents in the study unit. DO was the primary explanatory factor influencing the concentrations of many inorganic constituents. Arsenic, iron, and manganese concentrations increased as DO concentrations decreased, consistent with patterns expected as a result of reductive dissolution of iron and (or) manganese oxides in aquifer sediments. Molybdenum concentrations increased in anoxic conditions and in oxic conditions with high pH, reflecting two mechanisms for the mobilization of molybdenum—reductive dissolution and pH-dependent desorption under oxic conditions from aquifer sediments. Nitrate concentrations decreased as DO concentrations decreased which would be consistent with degradation of nitrate under anoxic conditions (denitrification). It also is possible that nitrate concentrations decreased in relation to increasing depth and groundwater age and not as a result of denitrification. Groundwater age was another explanatory factor frequently correlated to several inorganic constituents. Iron and manganese concentrations were higher in pre-modern (water recharged before 1952) or mixed-age groundwater. This correlation is one indication that iron and manganese are from natural sources. Nitrate, TDS, and sulfate concentrations were higher in modern groundwater (water recharged since 1952) and may indicate that human activities increase concentrations of nitrate, TDS, and sulfate. Land use was a third explanatory factor frequently correlated with inorganic constituents. Nitrate, TDS, and sulfate concentrations were higher in agricultural land-use areas than in natural land-use areas, indicating that increased concentrations may be a result of agricultural practices. Organic constituents usually were detected at low relative-concentrations; therefore, statistical analyses of relations to explanatory factors usually were done for classes of constituents (for example, pesticides or solvents) as well as for selected constituents. The number of VOCs detected in a well was not correlated to any of the explanatory factors evaluated. The number of pesticide and solvent detections and PCE and CFC-12 concentrations were higher in modern groundwater than in pre-modern groundwater. PCE and CFC-12 also were positively correlated to the density of LUFTs. PCE was negatively correlated to natural land use. Chloroform concentrations were positively correlated to the density of septic systems. Perchlorate concentrations were greater in agricultural areas than in urban or natural areas. Correlation of perchlorate with DO may indicate that perchlorate biodegradation under anoxic conditions may occur. Anthropogenic sources have contributed perchlorate to groundwater in the SCRC study unit, although low levels of perchlorate may occur naturally.

California↗
USGS2014

Simulation of groundwater flow in the Edwards-Trinity and related aquifers in the Pecos County region, Texas

The Edwards-Trinity aquifer is a vital groundwater resource for agricultural, industrial, and public supply uses in the Pecos County region of western Texas. The U.S. Geological Survey completed a comprehensive, integrated analysis of available hydrogeologic data to develop a numerical groundwater-flow model of the Edwards-Trinity and related aquifers in the study area in parts of Brewster, Jeff Davis, Pecos, and Reeves Counties. The active model area covers about 3,400 square miles of the Pecos County region of Texas west of the Pecos River, and its boundaries were defined to include the saturated areas of the Edwards-Trinity aquifer. The model is a five-layer representation of the Pecos Valley, Edwards-Trinity, Dockum, and Rustler aquifers. The Pecos Valley aquifer is referred to as the alluvial layer, and the Edwards-Trinity aquifer is divided into layers representing the Edwards part of the Edwards-Trinity aquifer and the Trinity part of the Edwards-Trinity aquifer, respectively. The calibration period of the simulation extends from 1940 to 2010. Simulated hydraulic heads generally were in good agreement with observed values; 1,684 out of 2,860 (59 percent) of the simulated values were within 25 feet of the observed value. The average root mean square error value of hydraulic head for the Edwards-Trinity aquifer was 34.2 feet, which was approximately 4 percent of the average total observed change in groundwater-level altitude (groundwater level). Simulated spring flow representing Comanche Springs exhibits a pattern similar to observed spring flow. Independent geochemical modeling corroborates results of simulated groundwater flow that indicates groundwater in the Edwards-Trinity aquifer in the Leon-Belding and Fort Stockton areas is a mixture of recharge from the Barilla and Davis Mountains and groundwater that has upwelled from the Rustler aquifer. The model was used to simulate groundwater-level altitudes resulting from prolonged pumping to evaluate sustainability of current and projected water-use demands. Each of three scenarios utilized a continuation of the calibrated model. Scenario 1 extended recent (2008) irrigation and nonirrigation pumping values for a 30-year period from 2010 to 2040. Projected groundwater-level changes in and around the Fort Stockton area under scenario 1 change little from current conditions, indicating that the groundwater system is near equilibrium with respect to recent (2008) pumping stress. Projected groundwater-level declines in the eastern part of the model area ranging from 5.0 to 15.0 feet are likely the result of nonequilibrium conditions associated with recent increases in pumping after a prolonged water-level recovery period of little or no pumping. Projected groundwater-level declines (from 15.0 to 31.0 feet) occurred in localized areas by the end of scenario 1 in the Leon-Belding area. Scenario 2 evaluated the effects of extended recent (2008) pumping rates as assigned in scenario 1 with year-round maximum permitted pumping rates in the Belding area. Results of scenario 2 are similar in water-level decline and extent as those of scenario 1. The extent of the projected groundwater-level decline in the range from 5.0 to 15.0 feet in the Leon-Belding irrigation area expanded slightly (about a 2-percent increase) from that of scenario 1. Maximum projected groundwater-level declines in the Leon-Belding irrigation area were approximately 31.3 feet in small isolated areas. Scenario 3 evaluated the effects of periodic increases in pumping rates over the 30-year extended period. Results of scenario 3 are similar to those of scenario 2 in terms of the areas of groundwater-level decline; however, the maximum projected groundwater-level decline increased to approximately 34.5 feet in the Leon-Belding area, and the extent of the decline was larger in area (about a 17-percent increase) than that of scenario 2. Additionally, the area of projected groundwater-level declines in the eastern part of the model area increased from that of scenario 2—two individual areas of decline coalesced into one larger area. The localized nature of the projected groundwater-level declines is a reflection of the high degree of fractured control on storage and hydraulic conductivity in the Edwards-Trinity aquifer. Additionally, the finding that simulated spring flow is highly dependent on the transient nature of hydraulic heads in the underlying aquifer indicates the importance of adequately understanding and characterizing the entire groundwater system.

Texas↗
USGS2015

Numerical simulation of groundwater flow, resource optimization, and potential effects of prolonged drought for the Citizen Potawatomi Nation Tribal Jurisdictional Area, central Oklahoma

A hydrogeological study including two numerical groundwater-flow models was completed for the Citizen Potawatomi Nation Tribal Jurisdictional Area of central Oklahoma. One numerical groundwater-flow model, the Citizen Potawatomi Nation model, encompassed the jurisdictional area and was based on the results of a regional-scale hydrogeological study and numerical groundwater flow model of the Central Oklahoma aquifer, which had a geographic extent that included the Citizen Potawatomi Nation Tribal Jurisdictional Area. The Citizen Potawatomi Nation numerical groundwater-flow model included alluvial aquifers not in the original model and improved calibration using automated parameter-estimation techniques. The Citizen Potawatomi Nation numerical groundwater-flow model was used to analyze the groundwater-flow system and the effects of drought on the volume of groundwater in storage and streamflow in the North Canadian River. A more detailed, local-scale inset model was constructed from the Citizen Potawatomi Nation model to estimate available groundwater resources for two Citizen Potawatomi Nation economic development zones near the North Canadian River, the geothermal supply area and the Iron Horse Industrial Park. Groundwater pumping rates at potential well locations were optimized using the most recent version of the U.S. Geological Survey Groundwater-Management Process for MODFLOW. The objectives of optimization were to determine if a total pumping rate of 500 gallons per minute could be pumped from 5 wells at the geothermal supply area and to maximize discharge from 16 wells at the Iron Horse Industrial Park without exceeding specified head drawdown constraints at the pumping wells and thus prevent groundwater depletion. The inset model was used to estimate North Canadian River streamflow depletion caused by optimized pumping at the Iron Horse Industrial Park because water quality was a concern, and the river may have degraded water quality compared to water in other parts of the alluvial aquifer. The fate of streamflow that infiltrates into groundwater because of pumping was not directly determined, but it was assumed that this water could end up in the well discharge, and was considered to be a maximum proportion of well discharge derived from the North Canadian River. The total optimized continuous pumping rate from five managed wells at the geothermal supply area was 638 gallons per minute, which exceeded the target pumping rate of 500 gallons per minute. The total continuous pumping rate from 16 wells at the Iron Horse Industrial Park was 1,472 gallons per minute, which induced stream infiltration of approximately 4.1 gallons per minute (approximately 0.3 percent of the total well discharge) from the North Canadian River. To estimate the effects of drought on water resources in the Citizen Potawatomi Nation Tribal Jurisdictional Area, a hypothetical 10-year drought during which precipitation would decrease by 50 percent was simulated by decreasing model groundwater recharge by the same proportion for the period 1990–2000 of the transient model. The effects of the drought were estimated by calculating the change in the volume of groundwater storage and groundwater flow to streams at the end of the drought period, and the change in simulated streamflow in the North Canadian River at the streamflow-gaging station at Shawnee, Okla., during and after the drought. The hypothetical decrease in recharge during the simulated drought caused groundwater in storage over the entire model in the study area to decrease by 361,500 acre-feet (14,100 acre-feet in the North Canadian River alluvial aquifer and 347,400 acre-feet in the Central Oklahoma aquifer), or approximately 0.2 percent of the total groundwater in storage over the drought period. This small percentage of groundwater loss showed that the Central Oklahoma aquifer as a bedrock aquifer has relatively low rates of recharge from the surface relative to the approximate storage. The budget for base flow to the North Canadian River indicated that the change in groundwater flow to the North Canadian River decreased during the 10-year drought by 386,500 acre-feet, or 37 percent. In all other parts of the Citizen Potawatomi Nation Tribal Jurisdictional Area, base flow decreased by 292,000 acre-feet, or 28 percent. Streamflow in the North Canadian River at the streamflow-gaging station at Shawnee, Okla., decreased during the hypothetical drought by as much as 28 percent, and the mean change in streamflow decreased as much as 16 percent. Streamflow at the Shawnee streamflow-gaging station did not recover to nondrought conditions until about 3 years after the simulated drought ended, during the relatively wet year of 2007.

Oklahoma↗
USGS2019

Groundwater quality and hydrology with emphasis on selenium mobilization and transport in the Lower Gunnison River Basin, Colorado, 2012–16

Dissolved selenium is a contaminant of concern in the lower Gunnison River Basin, Colorado. Selenium is naturally present in the Cretaceous Mancos Shale and is leached to groundwater and surface water by irrigation. The groundwater on the east side of the Uncompahgre River in Delta and Montrose Counties is one of the primary sources of selenium concentration and load to surface water in the lower Gunnison River Basin. Although little information about the contribution of groundwater to surface water has been historically available, groundwater has often been implicated as an appreciable source of selenium to surface water. From 2012 to 2016, the U.S. Geological Survey, in cooperation with the Bureau of Reclamation, the Colorado Water Conservation Board, and the Gunnison Basin Selenium Management Program, established a 30-well groundwater-monitoring network on irrigated land to characterize the hydrology and groundwater quality of the shallow groundwater system on the east side of the Uncompahgre River in the lower Gunnison River Basin. The installation of the 30-well network and the data collected allowed for the development of a conceptual model of selenium mobilization and transport in the shallow groundwater system. Monitoring wells were completed in surficial deposits and in weathered Mancos Shale, which generally exhibited unconfined and confined conditions, respectively. Groundwater-quality monitoring provides information on the distribution of selenium and the geochemical processes controlling selenium concentrations in shallow groundwater. Monitoring wells were sampled between August 2013 and March 2015 to understand groundwater quality, seasonality, sources of recharge, and groundwater age. Concentrations of dissolved selenium ranged from below the limit of detection to 4,100 micrograms per liter (µg/L), with a median concentration of 14 µg/L. Concentrations showed a high degree of spatial variability and no seasonal difference. Similarly, no seasonal pattern was observed in specific conductance values of groundwater despite the considerably lower specific conductance value of irrigation water. Reduction-oxidation processes are important controls on selenium mobility. Nitrate derived from geologic material was a primary control on reduction-oxidation conditions in groundwater and inhibited selenium reduction to less mobile forms. Nitrate was reduced by denitrification in groundwater, but it was not reduced to the extent necessary to allow for selenium reduction. Groundwater ages were determined for groundwater samples from eight wells and ranged from 6 to 20 years old. Isotopic data indicate groundwater was recharged by irrigation water; no information collected supported an older, deeper source of recharge to the shallow groundwater system. Data on water level in all wells showed response to irrigation practices, but the response was delayed in some wells, which may be an indication of distance from recharge source.

Colorado↗
USGS2021

Characterization of groundwater quality and discharge with emphasis on selenium in an irrigated agricultural drainage near Delta, Colorado, 2017–19

Selenium is a water-quality constituent of concern for aquatic ecosystems in the lower Gunnison River Basin. Selenium is derived from bedrock of the Mancos Shale and is mobilized and transported to groundwater and surface water by application of irrigation water. Although it is recognized that groundwater contributes an appreciable amount of selenium to surface water, few studies have addressed interactions between the two. The U.S. Geological Survey in cooperation with the Colorado Water Conservation Board conducted a study during 2017–19 to characterize the quality and quantity of groundwater discharging to an agricultural drainage near Delta, Colorado, locally known as Sunflower Drain. Water quality in the study area is characterized by high dissolved solids with elevated concentrations of selenium and nitrate resulting from dissolution of soluble salts in the Mancos Shale. Selenium concentrations have decreased by 50 percent since the early 2000s, possibly in response to irrigation system improvements. Stable water isotopes indicate streamflow is dominated by canal water during the irrigation season (April to October) and, during the nonirrigation season (November to March), is dominated by groundwater that has undergone some degree of evaporation. Pesticide and pharmaceutical compounds were infrequently detected, and results indicate they were derived from sources outside the study area such that they do not appear to be useful as tracers of groundwater sources. Stable isotopes of nitrate indicate that nitrate originates from the Mancos Shale, and the isotopic composition is enriched by denitrification in the groundwater system. Using a mass-balance approach, estimated groundwater discharge rates to Sunflower Drain ranged from 0.15 to 0.27 cubic feet per second per mile with one losing reach identified. Selenium, sulfate, and nitrate concentrations in groundwater estimated by mass-balance calculations were similar to concentrations measured in the Poly 17 observation well, located in a largely irrigated area in east tributary. One tributary reach had higher concentrations of selenium, sulfate, and nitrate likely reflecting localized inputs of more concentrated groundwater, similar to the concentrations in the Poly 7 observation well, which is downgradient from a residential area in the west tributary. Three pilot studies were conducted, including fiber optic distributed temperature sensing to detect groundwater discharge zones in the stream channel, a passive seismic technique to estimate depth to bedrock, and use of radon-222 as a geochemical tracer of groundwater discharge. All three techniques show promise as additional approaches for investigating groundwater discharge surface-water systems in irrigated drainage areas on Mancos Shale. The factors that affect groundwater movement mainly include when and where irrigation water is transported and applied, and the distribution of bedrock of the Mancos Shale and overlying alluvial deposits. The average groundwater recharge rate for the study area was estimated at 8.1 inches per year, based on mass balance calculations from synoptic survey data. Along the western tributary of Sunflower Drain, there was evidence that spills from the East Canal may recharge the groundwater aquifer adjacent to the stream channel. Groundwater movement to the stream channel may be controlled by the topography of the alluvial/bedrock interface or focused along human-made features, such as tile drains and ditches constructed around irrigated fields. On larger scales, the top of bedrock was also important, creating a topographic constriction that caused a zone of groundwater discharge. The groundwater system is complex, and further study could better define the system, possibly through application of a groundwater flow model and more extensive studies using some of the exploratory methods evaluated in this study.

Colorado↗
USGS2024

Simulated mean monthly groundwater-transported nitrogen loads in watersheds on the north shore of Long Island Sound, 1993–2022

Elevated nitrogen loads are pervasive in the Long Island Sound, an estuary that receives freshwater and nutrients from both surface-water and groundwater discharge. Surface-water nitrogen loads to the Long Island Sound are relatively well characterized, but less is known about groundwater-transported nitrogen loads. Prior work on the northern shore of Long Island Sound (Connecticut and areas of New York and Rhode Island) suggested that groundwater travel times are relatively short (median less than 2 years) and that decade-long nutrient legacies are not widespread. Because the travel times are short, groundwater flow and nutrient loads likely vary substantially between months. In the current study, the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency’s Long Island Sound Study and the Connecticut Department of Energy and Environmental Protection, developed a set of models to better characterize spatial and temporal patterns of groundwater-transported nitrogen loading from atmospheric deposition, septic systems, and fertilizers within the study area. The models provide an estimate, with uncertainty, of groundwater-transported nitrogen loads in the study area, filling a key gap in the nitrogen budget for Long Island Sound. The models also highlight the spatial and temporal variation in nitrogen loading throughout the study area. The modeling workflow involved four models. (1) A soil-water-balance model was developed by using the Soil-Water-Balance software to simulate groundwater recharge across the study area for water years 2005 through 2022. The simulated mean monthly recharge from the soil-water-balance model was used as input into a groundwater-flow model. (2) The groundwater-flow model was developed by using the MODFLOW 6 software and data for water years 1993 through 2022 and simulates average monthly hydrologic conditions. The groundwater-flow model was calibrated by using the Iterative Ensemble Smoother method within the PEST++ software. The Iterative Ensemble Smoother method generates an ensemble of sets of parameter values, with each set producing reasonable simulated hydrologic parameter values. (3) An ensemble of MODPATH particle-tracking simulations were run to generate particle flow paths and travel times, with each simulation using a different set of the flow model parameters. (4) A nitrogen load model uses the MODPATH simulation outputs to track nitrogen from the land surface through multiple attenuation zones until it discharges into fresh or saline surface water. As with the groundwater-flow model, the nitrogen model simulated average monthly groundwater-transported nitrogen loads for water years 1993 through 2022. One novel aspect of the nitrogen load model is that the nitrogen attenuation parameters were calibrated to observed nitrogen loads. Across the ensemble of simulated nitrogen loads, the median study-area-wide monthly simulated nitrogen loads from the aquifer to Long Island Sound throughout the year ranged from 900 to 18,600 kilograms of nitrogen per day, with a median load of 5,100 kilograms of nitrogen per day. The simulated loads were based on average monthly conditions for water years 1993 through 2022. Loads were highest during the winter and early spring and lowest during the late summer. However, simulated travel times for groundwater and nitrogen loads discharged to Long Island Sound during summer were longer than travel times for groundwater and loads discharged during the winter, indicating that, on average, groundwater discharged during summer traveled along different, and longer, flow paths, than groundwater discharged during winter. This indicates that summer loads would respond more slowly to changes in nitrogen inputs at the water table than winter loads. Over the entire study area, approximately 15 percent of the simulated load is from atmospheric deposition sources, 30 to 40 percent is from fertilizer, and 50 to 60 percent is from septic systems. The final analysis of the study involved simulating the change in groundwater-transported nitrogen load in response to upgrading septic systems or reducing fertilizing inputs to areas of turf grass. Both management interventions reduced the groundwater-transported nitrogen load, and reductions were greater in areas with greater loads from septic systems or turf-grass fertilizers. The delay between management actions and substantial reductions in groundwater-transported nitrogen loads varied seasonally; loads during the late summer months remained elevated longer than the winter loads.

Connecticut, Rhode Island↗
USGS2025

Quality of groundwater used for domestic supply in the Gilroy-Hollister basin and surrounding areas, California, 2022

More than 2 million Californians rely on groundwater from domestic wells for drinking-water supply. This report summarizes a 2022 California Groundwater Ambient Monitoring and Assessment Priority Basin Project (GAMA-PBP) water-quality survey of 33 domestic and small-system drinking-water supply wells in the Gilroy-Hollister Valley groundwater basin and the surrounding areas, where more than 20,000 residents are estimated to utilize privately owned domestic wells. The study area includes the Llagas subbasin in the north, the North San Benito subbasin in the south, and the surrounding uplands. The study was focused on groundwater resources used for domestic drinking-water supply, which are mostly drawn from shallower parts of aquifer systems rather than those of groundwater resources used for public drinking-water supply in the same area. This assessment characterized the quality of ambient groundwater in the aquifer before filtration or treatment, rather than the quality of drinking water delivered to the tap. To provide context, the measured concentrations of constituents in groundwater were compared to Federal and California State regulatory and non-regulatory benchmarks for drinking-water quality. A grid-based method was used to estimate the areal proportions of groundwater resources used for domestic drinking wells that have water-quality constituents present at high concentrations (above the benchmark), moderate concentrations (between one-half of the benchmark and the benchmark for inorganic constituents, or between one-tenth of the benchmark and the benchmark for organic constituents), and low concentrations (less than one-half or one-tenth the benchmark for inorganic and organic constituents, respectively). This method provides statistically representative results at the study-area scale and permits comparisons to other GAMA-PBP study areas. In the study area, inorganic constituents in groundwater were greater than regulatory benchmarks (U.S. Environmental Protection Agency [EPA] or State of California maximum contaminant levels [MCLs]) for public drinking-water quality in 24 percent of domestic groundwater resources. The inorganic constituents present at concentrations greater than MCLs for drinking water were nitrate (as nitrogen), barium, chromium, and selenium. Total dissolved solids (TDS) or manganese were present at concentrations greater than the secondary maximum contaminant levels (SMCLs) that the State of California uses as aesthetic-based benchmarks in 48 percent of domestic groundwater resources. No volatile organic compounds or pesticide constituents were present at concentrations greater than regulatory benchmarks. Total coliform bacteria and enterococci were detected in 4 percent of domestic groundwater resources. Per- and polyfluoroalkyl substances (PFAS) were detected in 19 percent of domestic groundwater resources, and 10 percent had concentrations greater than recently enacted (April 2024) EPA MCLs. Physical and chemical factors from natural and anthropogenic sources that could affect the groundwater quality were evaluated using results from statistical testing of associations between constituent concentrations and potential explanatory variables. In this study, relevant physical factors include well construction characteristics, groundwater age, site proximity to groundwater recharge or discharge zones, and potential sources of contamination. Relevant chemical factors include the initial chemistry of the recharge water, the mineralogy of the aquifer sediments, and the subsequent shifts in chemistry as biologic and geologic reactions alter groundwater in the subsurface. Nitrate concentrations were correlated to agricultural land use, distance from the boundary of the Gilroy-Hollister Valley groundwater basin, and the proportion of modern (post-1950s) water captured by the well. Denitrification under anoxic redox conditions can mitigate some nitrate derived from fertilizer application. Total dissolved solids primarily were derived from water-rock interactions with soils and aquifer materials in the study area, but there were high concentrations where agricultural practices contributed additional TDS. Mineralogy of aquifer sediments and rocks also affect barium, selenium, boron, and chromium concentrations in the Gilroy-Hollister Valley groundwater basin. PFAS were positively correlated with urban land use and the proportion of modern water captured by the well.

California↗
USGS2026

Streamflow, base-flow, and precipitation trends and simulated effects of groundwater withdrawals from the North Fork Red River aquifer on base flows upgradient from Lake Altus, western Oklahoma, 1980–2022

The U.S. Geological Survey, in cooperation with the Bureau of Reclamation, used five scenarios created from a previously published numerical groundwater-flow model (1980–2013) and historical streamflow records (1980–2022) to investigate the relation between groundwater withdrawals from the North Fork Red River aquifer and inflows to Lake Altus from the North Fork Red River in western Oklahoma. The five scenarios were (1) a scaled-equal-proportionate-share (EPS) groundwater-withdrawal scenario, (2) a study-area-scaled-reported groundwater-withdrawal scenario, (3) a zonal-scaled-reported groundwater-withdrawal scenario, (4) a historical drought-threshold scenario, and (5) a base-flow and evapotranspiration depletion scenario. For the scaled-EPS groundwater-withdrawal scenario, EPS groundwater withdrawals were often much higher than reported groundwater withdrawals and greatly decreased base flows for most scale factors. For the study-area-scaled-reported groundwater-withdrawal scenario, base flows were reduced more but by smaller percentages during wet periods than during dry periods when scaling simulated reported groundwater withdrawals. For the zonal-scaled-reported groundwater-withdrawal scenario, scaling simulated reported groundwater withdrawals within selected zones with more groundwater withdrawals did not always affect base flows more than scaling reported groundwater withdrawals within zones with less groundwater withdrawals. For the historical drought-threshold scenario, curtailing groundwater withdrawals at the drought thresholds increased annual base flows to Lake Altus by about 1,169 to 3,665 acre-feet. For the base-flow and evapotranspiration depletion scenario, the distance between a groundwater well and a stream was a major factor affecting base flow to the North Fork Red River when increasing groundwater withdrawals; however, spatially variable hydrologic properties and saturated-zone evapotranspiration could also affect the relation between base flows and groundwater withdrawals.

Oklahoma↗
USGS2017

Groundwater flow model for the Little Plover River basin in Wisconsin’s Central Sands

The Little Plover River is a groundwater-fed stream in the sand plains region of central Wisconsin. In this region, sandy sediment deposited during or soon after the last glaciation forms an important unconfined sand and gravel aquifer. This aquifer supplies water for numerous high-capacity irrigation, municipal, and industrial wells that support a thriving agricultural industry. In recent years the addition of many new wells, combined with observed diminished flows in the Little Plover and other nearby rivers, has raised concerns about the impacts of the wells on groundwater levels and on water levels and flows in nearby lakes, streams, and wetlands. Diverse stakeholder groups, including well operators, Growers, environmentalists, local land owners, and regulatory and government officials have sought a better understanding of the local groundwater-surface water system and have a shared desire to balance the water needs of the he liagricultural, industrial, and urban users with the maintenance and protection of groundwater-dependent natural resources. To help address these issues, the Wisconsin Department of Natural Resources requested that the Wisconsin Geological and Natural History Survey and U.S. Geological Survey cooperatively develop a groundwater flow model that could be used to demonstrate the relationships among groundwater, surface water, and well withdrawals and also be a tool for testing and evaluating alternative water management strategies for the central sands region. Because of an abundance of previous studies, data availability, local interest, and existing regulatory constraints the model focuses on the Little Plover River watershed, but the modeling methodology developed during this study can apply to much of the larger central sands of Wisconsin. The Little Plover River groundwater flow model simulates three-dimensional groundwater movement in and around the Little Plover River basin under steady-state and transient conditions. This model explicitly includes all high-capacity wells in the model domain and simulates seasonal variations in recharge and well pumping. The model represents the Little Plover River, and other significant streams and drainage ditches in the model domain, as fully connected to the groundwater system, computes stream base flow resulting from groundwater discharge, and routes the flow along the stream channel. A separate soil-water-balance (SWB) model was used to develop groundwater recharge arrays as input for the groundwater flow model. The SWB model uses topography, soils, land use, and climatic data to estimate recharge as deep drainage from the soil zone. The SWB model explicitly includes recharge originating as irrigation water, and computes irrigation using techniques similar to those used by local irrigation operators. The groundwater flow model uses the U.S. Geological Survey’s MODFLOW modeling code which is freely available, widely accepted, and commonly used by the groundwater community. The groundwater flow model and the SWB model use identical high-resolution numerical grids having model cells 100 feet on a side, with physical properties assigned to each grid cell. This grid allows accurate geographic placement of wells, streams, and other model features. The 3-dimensional grid has three layers; layers 1 and 2 represent the sand and gravel aquifer and layer 3 represents the underlying sandstone. The distribution of material properties in the model (hydraulic conductivity, aquifer thickness, etc.) comes from previous published geologic studies of the region, updated by calibration to recent streamflow and groundwater level data. The SWB model operates on a daily time step. The groundwater flow model was calibrated to monthly stress periods with time steps ranging from 1 to 16 days. More detailed time discretization is possible. The groundwater model was calibrated to water-level and streamflow data collected during 2013 and 2014 by adjusting model parameters (primarily hydraulic conductivity, storage, and recharge) until the model produced a conditionally optimal fit between field observations and model output, subject to consistency with previously published geologic studies. Calibration was performed under both steady and transient conditions, and used a sophisticated parameter-estimation procedure (PEST) for the calibration process and to identify important model parameters. For the Little Plover River, the two most important parameters are the global recharge multiplier and the hydraulic conductivity of the stream bed. The calibrated model produces water-level and mass-balance results that are consistent with field observations and previous studies of the area. The completed model is a powerful tool for testing and demonstrating alternative water-management scenarios. Example model applications described in this report include simulating how the cumulative impacts of pumping and land-use change have affected average baseflow in the Little Plover River. Depletion-potential mapping represents a method for predicting which wells and well locations have the greatest impact on nearby surface-water resources. The completed model is publicly available, along with a companion user’s guide to assist with its operation, at http://wgnhs.org/littleplover- river-groundwater-model.

Wisconsin↗
USGS1996

Nitrate in groundwater and water sources used by riparian trees in an agricultural watershed: A chemical and isotopic investigation in southern Minnesota

This study evaluates processes that affect nitrate concentrations in groundwater beneath riparian zones in an agricultural watershed. Nitrate pathways in the upper 2 m of groundwater were investigated beneath wooded and grass-shrub riparian zones next to cultivated fields. Because trees can be important components of the overall nitrate pathway in wooded riparian zones, water sources used by riparian trees and possible effects of trees on nitrate concentrations in groundwater were also investigated. Average nitrate concentrations in shallow groundwater beneath the cultivated fields were 5.5 mg/L upgradient of the wooded riparian zone and 3.5 mg/L upgradient of the grass-shrub zone. Shallow groundwater beneath the fields passed through the riparian zones and discharged into streams that had average nitrate concentrations of 8.5 mg/L (as N). Lateral variations of δD values in groundwater showed that mixing among different water sources occurred beneath the riparian zones. In the wooded riparian zone, nitrate concentrations in shallow groundwater were diluted by upwelling, nitrate-poor, deep groundwater. Upwelling deep groundwater contained ammonium with a δ 15 N of 5‰ that upon nitrification and mixing with nitrate in shallow groundwater caused nitrate δ 15 N values in shallow groundwater to decrease by as much as 19.5‰. Stream water penetrated laterally beneath the wooded riparian zone as far as 19 m from the stream's edge and beneath the grass-shrub zone as far as 27 m from the stream's edge. Nitrate concentrations in shallow groundwater immediately upgradient of where it mixed with stream water averaged 0.4 mg/L in the wooded riparian zone and 0.8 mg/L near the grass-shrub riparian zone. Nitrate concentrations increased toward the streams because of mixing with nitrate-rich stream water. Because nitrate concentrations were larger in stream water than shallow groundwater, concentrated nitrate in the streams cannot have come from shallow groundwater at these sites. Water sources of riparian trees were identified by comparing δD values of sap water, soil water, groundwater, and stream water. Soil water was the main water source for trees in the outer 4 to 6 m of one part of the wooded riparian zone and outer 10 m of another part. Groundwater was a significant water source for trees closer to the streams where the water table was less than about 2.1 to 2.7 m below the surface. No evidence was found in the nitrate concentration profiles that trees close to the streams that took up groundwater through their roots also took up nitrate from groundwater. The lack of such evidence is attributed to the nitrate concentration profiles being insufficiently sensitive indicators of nitrate removal by trees.

Minnesota↗
USGS2010

Groundwater-quality data for the Sierra Nevada study unit, 2008: Results from the California GAMA program

Groundwater quality in the approximately 25,500-square-mile Sierra Nevada study unit was investigated in June through October 2008, as part of the Priority Basin Project of the Groundwater Ambient Monitoring and Assessment (GAMA) Program. The GAMA Priority Basin Project is being conducted by the U.S. Geological Survey (USGS) in cooperation with the California State Water Resources Control Board (SWRCB). The Sierra Nevada study was designed to provide statistically robust assessments of untreated groundwater quality within the primary aquifer systems in the study unit, and to facilitate statistically consistent comparisons of groundwater quality throughout California. The primary aquifer systems (hereinafter, primary aquifers) are defined by the depth of the screened or open intervals of the wells listed in the California Department of Public Health (CDPH) database of wells used for public and community drinking-water supplies. The quality of groundwater in shallower or deeper water-bearing zones may differ from that in the primary aquifers; shallow groundwater may be more vulnerable to contamination from the surface. In the Sierra Nevada study unit, groundwater samples were collected from 84 wells (and springs) in Lassen, Plumas, Butte, Sierra, Yuba, Nevada, Placer, El Dorado, Amador, Alpine, Calaveras, Tuolumne, Madera, Mariposa, Fresno, Inyo, Tulare, and Kern Counties. The wells were selected on two overlapping networks by using a spatially-distributed, randomized, grid-based approach. The primary grid-well network consisted of 30 wells, one well per grid cell in the study unit, and was designed to provide statistical representation of groundwater quality throughout the entire study unit. The lithologic grid-well network is a secondary grid that consisted of the wells in the primary grid-well network plus 53 additional wells and was designed to provide statistical representation of groundwater quality in each of the four major lithologic units in the Sierra Nevada study unit: granitic, metamorphic, sedimentary, and volcanic rocks. One natural spring that is not used for drinking water was sampled for comparison with a nearby primary grid well in the same cell. Groundwater samples were analyzed for organic constituents (volatile organic compounds [VOC], pesticides and pesticide degradates, and pharmaceutical compounds), constituents of special interest (N-nitrosodimethylamine [NDMA] and perchlorate), naturally occurring inorganic constituents (nutrients, major ions, total dissolved solids, and trace elements), and radioactive constituents (radium isotopes, radon-222, gross alpha and gross beta particle activities, and uranium isotopes). Naturally occurring isotopes and geochemical tracers (stable isotopes of hydrogen and oxygen in water, stable isotopes of carbon, carbon-14, strontium isotopes, and tritium), and dissolved noble gases also were measured to help identify the sources and ages of the sampled groundwater. Three types of quality-control samples (blanks, replicates, and samples for matrix spikes) each were collected at approximately 10 percent of the wells sampled for each analysis, and the results for these samples were used to evaluate the quality of the data for the groundwater samples. Field blanks rarely contained detectable concentrations of any constituent, suggesting that contamination from sample collection, handling, and analytical procedures was not a significant source of bias in the data for the groundwater samples. Differences between replicate samples were within acceptable ranges, with few exceptions. Matrix-spike recoveries were within acceptable ranges for most compounds. This study did not attempt to evaluate the quality of water delivered to consumers; after withdrawal from the ground, groundwater typically is treated, disinfected, or blended with other waters to maintain water quality. Regulatory benchmarks apply to finished drinking water that is served to the consumer, not to untreated groundwater. However, to provide some context for the results, concentrations of constituents measured in the groundwater were compared with regulatory and nonregulatory health-based benchmarks established by the U.S. Environmental Protection Agency (USEPA) and CDPH and with nonregulatory aesthetic and technical benchmarks established by CDPH. Comparisons between data collected for this study and drinking-water benchmarks are for illustrative purposes only and do not indicate compliance or noncompliance with regulatory benchmarks. All organic constituents and most inorganic constituents that were detected in groundwater samples from the 30 primary grid wells in the Sierra Nevada study unit were detected at concentrations less than drinking-water benchmarks. Of the 150 organic and special-interest constituents analyzed, 21 were detected in groundwater samples; all concentrations were less than regulatory and nonregulatory health-based benchmarks, and most were less than 1/10th of benchmark levels. One or more organic constituents were detected in 37 percent of the primary grid wells, and perchlorate was detected in 27 percent of the primary grid wells. Most samples analyzed for inorganic and radioactive constituents had concentrations or activities less than regulatory and nonregulatory health-based benchmarks. Nutrients were not detected at concentrations greater than health-based benchmarks. Six of the 30 primary grid wells (20 percent) and 7 of the 53 lithologic grid wells had concentrations of or activities for one or two constituents that were greater than the benchmarks. Constituents present in one or more samples at concentrations or activities greater than health-based benchmarks were arsenic (5 wells, MCL-US), gross alpha particle activity (4 wells, MCL-US), boron (2 wells, NL-CA), fluoride (1 well, MCL-CA), and selenium (1 well, MCL-US). Two of the wells that had high gross alpha particle activities had uranium concentrations (MCL-CA) and uranium activities (MCL-CA) greater than the benchmark levels. Four of the 29 samples analyzed had activities of radon-222 greater than the proposed alternative MCL-US. Most samples analyzed for inorganic constituents that had nonregulatory, aesthetic-based benchmarks (SMCLs) had concentrations less than these benchmarks. Total dissolved solids concentrations were less than the upper SMCL-CA in all 83 primary and lithologic grid well samples, and TDS concentrations were less than the recommended SMCL-CA in 79 of these samples. Manganese concentrations were greater than the SMCL-CA in 2 of the 30 primary grid wells (7 percent) and in 6 of the 53 lithologic grid wells, and iron concentrations were greater than the SMCL-CA in the same 2 primary grid wells and in 5 of the same lithologic grid wells.

California↗
USGS2013

Groundwater-quality data in the Cascade Range and Modoc Plateau study unit, 2010-Results from the California GAMA Program

Groundwater quality in the 39,000-square-kilometer Cascade Range and Modoc Plateau (CAMP) study unit was investigated by the U.S. Geological Survey (USGS) from July through October 2010, as part of the California State Water Resources Control Board (SWRCB) Groundwater Ambient Monitoring and Assessment (GAMA) Program’s Priority Basin Project (PBP). The GAMA PBP was developed in response to the California Groundwater Quality Monitoring Act of 2001 and is being conducted in collaboration with the SWRCB and Lawrence Livermore National Laboratory (LLNL). The CAMP study unit is the thirty-second study unit to be sampled as part of the GAMA PBP. The GAMA CAMP study was designed to provide a spatially unbiased assessment of untreated-groundwater quality in the primary aquifer system and to facilitate statistically consistent comparisons of untreated-groundwater quality throughout California. The primary aquifer system is defined as that part of the aquifer corresponding to the open or screened intervals of wells listed in the California Department of Public Health (CDPH) database for the CAMP study unit. The quality of groundwater in shallow or deep water-bearing zones may differ from the quality of groundwater in the primary aquifer system; shallow groundwater may be more vulnerable to surficial contamination. In the CAMP study unit, groundwater samples were collected from 90 wells and springs in 6 study areas (Sacramento Valley Eastside, Honey Lake Valley, Cascade Range and Modoc Plateau Low Use Basins, Shasta Valley and Mount Shasta Volcanic Area, Quaternary Volcanic Areas, and Tertiary Volcanic Areas) in Butte, Lassen, Modoc, Plumas, Shasta, Siskiyou, and Tehama Counties. Wells and springs were selected by using a spatially distributed, randomized grid-based method to provide statistical representation of the study unit (grid wells). Groundwater samples were analyzed for field water-quality indicators, organic constituents, perchlorate, inorganic constituents, radioactive constituents, and microbial indicators. Naturally occurring isotopes and dissolved noble gases also were measured to provide a dataset that will be used to help interpret the sources and ages of the sampled groundwater in subsequent reports. In total, 221 constituents were investigated for this study. Three types of quality-control samples (blanks, replicates, and matrix spikes) were collected at approximately 10 percent of the wells in the CAMP study unit, and the results for these samples were used to evaluate the quality of the data for the groundwater samples. Blanks rarely contained detectable concentrations of any constituent, suggesting that contamination from sample collection procedures was not a significant source of bias in the data for the groundwater samples. Replicate samples generally were within the limits of acceptable analytical reproducibility. Matrix-spike recoveries were within the acceptable range (70 to 130 percent) for approximately 90 percent of the compounds. This study did not attempt to evaluate the quality of water delivered to consumers; after withdrawal from the ground, untreated groundwater typically is treated, disinfected, and (or) blended with other waters to maintain water quality. Regulatory benchmarks apply to water that is served to the consumer, not to untreated groundwater. However, to provide some context for the results, concentrations of constituents measured in the untreated groundwater were compared with regulatory and non-regulatory health-based benchmarks established by the U.S. Environmental Protection Agency (USEPA) and CDPH, and to non-regulatory benchmarks established for aesthetic concerns by CDPH. Comparisons between data collected for this study and benchmarks for drinking water are for illustrative purposes only and are not indicative of compliance or non-compliance with those benchmarks. All organic constituents and most inorganic constituents that were detected in groundwater samples from the 90 grid wells in the CAMP study unit were detected at concentrations less than drinking-water benchmarks. Of the 148 organic constituents analyzed, 27 were detected in groundwater samples; concentrations of all detected constituents were less than regulatory and nonregulatory health-based benchmarks, and all were less than 1/10 of benchmark levels. One or more organic constituents were detected in 52 percent of the grid wells in the CAMP study unit: VOCs were detected in 30 percent, and pesticides and pesticide degradates were detected in 31 percent. Trace elements, major ions, nutrients, and radioactive constituents were sampled for at 90 grid wells in the CAMP study unit, and most detected concentrations were less than health-based benchmarks. Exceptions include three detections of arsenic greater than the USEPA maximum contaminant level (MCL-US) of 10 micrograms per liter (µg/L), two detections of boron greater than the CDPH notification level (NL-CA) of 1,000 µg/L, two detections of molybdenum greater than the USEPA lifetime health advisory level (HAL-US) of 40 µg/L, two detections of vanadium greater than the CDPH notification level (NL-CA) of 50 µg/L, one detection of nitrate, as nitrogen, greater than the MCL-US of 10 milligrams per liter (mg/L), two detections of uranium greater than the MCL-US of 30 µg/L and the MCL-CA of 20 picocuries per liter (pCi/L), one detection of radon-222 greater than the proposed MCL-US of 4,000 pCi/L, and two detections of gross alpha particle activity greater than the MCL-US of 15 pCi/L. Results for inorganic constituents with non-regulatory benchmarks set for aesthetic concerns showed that iron concentrations greater than the CDPH secondary maximum contaminant level (SMCL-CA) of 300 µg/L were detected in four grid wells. Manganese concentrations greater than the SMCL-CA of 50 µg/L were detected in nine grid wells. Chloride and TDS were detected at concentrations greater than the upper SMCL-CA benchmarks of 500 mg/L and 1,000 mg/L, respectively, in one grid well. Microbial indicators (total coliform and Escherichia coli [E. coli]) were detected in 11 percent of the 83 grid wells sampled for these analyses in the CAMP study unit. The presence of total coliform was detected in nine grid wells, and the presence of E. coli was detected in one of these same grid wells.

California↗
USGS2012

Preliminary investigation of the effects of sea-level rise on groundwater levels in New Haven, Connecticut

Global sea level rose about 0.56 feet (ft) (170 millimeters (mm)) during the 20th century. Since the 1960s, sea level has risen at Bridgeport, Connecticut, about 0.38 ft (115 mm), at a rate of 0.008 ft (2.56 mm + or - 0.58 mm) per year. With regional subsidence, and with predicted global climate change, sea level is expected to continue to rise along the northeast coast of the United States through the 21st century. Increasing sea levels will cause groundwater levels in coastal areas to rise in order to adjust to the new conditions. Some regional climate models predict wetter climate in the northeastern United States under some scenarios. Scenarios for the resulting higher groundwater levels have the potential to inundate underground infrastructure in lowlying coastal cities. New Haven is a coastal city in Connecticut surrounded and bisected by tidally affected waters. Monitoring of water levels in wells in New Haven from August 2009 to July 2010 indicates the complex effects of urban influence on groundwater levels. The response of groundwater levels to recharge and season varied considerably from well to well. Groundwater temperatures varied seasonally, but were warmer than what was typical for Connecticut, and they seem to reflect the influence of the urban setting, including the effects of conduits for underground utilities. Specific conductance was elevated in many of the wells, indicating the influence of urban activities or seawater in Long Island Sound. A preliminary steady-state model of groundwater flow for part of New Haven was constructed using MODFLOW to simulate current groundwater levels (2009-2010) and future groundwater levels based on scenarios with a rise of 3 ft (0.91 meters (m)) in sea level, which is predicted for the end of the 21st century. An additional simulation was run assuming a 3-ft rise in sea level combined with a 12-percent increase in groundwater recharge. The model was constructed from existing hydrogeologic information for the New Haven area and from new information on groundwater levels collected during October 2009-June 2010. For the scenario with a 3-ft rise in sea level and no increase in recharge, simulated groundwater levels near the coast rose 3 ft; this increased water level tapered off toward a discharge area at the only nontidal stream in the study area. Simulated stream discharge increased at the nontidal stream because of the increased gradient. Although groundwater levels rose, the simulated difference between the groundwater levels in the aquifer and the increased sea level declined, indicating that the depth to the interface between freshwater and saltwater may possibly decline. Simulated water levels were affected by rise in sea level even in areas where the water table was at 17-24 ft (5.2-7.3 m) above current (2011) sea level. For the scenario with increased recharge, simulated groundwater levels were as much as an additional foot higher at some locations in the study area. The results of this preliminary investigation indicate that groundwater levels in coastal areas can be expected to rise and may rise higher if groundwater recharge also increases. This finding has implications for the disposal of stormwater through infiltration, a low-impact development practice designed to improve water quality and reduce overland peak discharge. Other implications include increased risk of basement flooding and increased groundwater seepage into underground sewer pipes and utility corridors in some areas. These implications will present engineering challenges to New Haven and Yale University. The preliminary model developed for this study can be the starting point for further simulation of future alternative scenarios for sea-level rise and recharge. Further simulations could identify those areas of New Haven where infrastructure may be at greatest risk from rising levels of groundwater. The simulations described in this report have limitations due to the preliminary scope of the work. Approaches to improve simulations include but are not limited to incorporating: * The variable density of seawater into the model in order to understand the current and future location of the interface between freshwater and saltwater; * Collection of additional data in order to better resolve temporal and spatial patterns in water levels in the aquifer; * Improved estimates of recharge through direct and indirect measurements of freshwater discharge from the study area; and * Transient simulations for greater understanding of the amount of time required for water levels and the position of the interface between freshwater and saltwater to adjust to changes in sea level and recharge.

Connecticut↗
USGS2011

Conceptual model and numerical simulation of the groundwater-flow system of Bainbridge Island, Washington

Groundwater is the sole source of drinking water for the population of Bainbridge Island. Increased use of groundwater supplies on Bainbridge Island as the population has grown over time has created concern about the quantity of water available and whether saltwater intrusion will occur as groundwater usage increases. A groundwater-flow model was developed to aid in the understanding of the groundwater system and the effects of groundwater development alternatives on the water resources of Bainbridge Island. Bainbridge Island is underlain by unconsolidated deposits of glacial and nonglacial origin. The surficial geologic units and the deposits at depth were differentiated into aquifers and confining units on the basis of areal extent and general water-bearing characteristics. Eleven principal hydrogeologic units are recognized in the study area and form the basis of the groundwater-flow model. A transient variable-density groundwater-flow model of Bainbridge Island and the surrounding area was developed to simulate current (2008) groundwater conditions. The model was calibrated to water levels measured during 2007 and 2008 using parameter estimation (PEST) to minimize the weighted differences or residuals between simulated and measured hydraulic head. The calibrated model was used to make some general observations of the groundwater system in 2008. Total flow through the groundwater system was about 31,000 acre-ft/ yr. The recharge to the groundwater system was from precipitation and septic-system returns. Groundwater flow to Bainbridge Island accounted for about 1,000 acre-ft/ yr or slightly more than 5 percent of the recharge amounts. Groundwater discharge was predominately to streams, lakes, springs, and seepage faces (16,000 acre-ft/yr) and directly to marine waters (10,000 acre-ft/yr). Total groundwater withdrawals in 2008 were slightly more than 6 percent (2,000 acre-ft/yr) of the total flow. The calibrated model was used to simulate predevelopment conditions, during which no groundwater pumping or secondary recharge occurred and currently developed land was covered by conifer forests. Simulated water levels in the uppermost aquifer generally were slightly higher at the end of 2008 than under predevelopment conditions, likely due to increased recharge from septic returns and reduced evapotranspiration losses due to conversion of land cover from forests to current conditions. Simulated changes in water levels for the extensively used sea-level aquifer were variable, although areas with declines between zero and 10 feet were common and generally can be traced to withdrawals from public-supply drinking wells. Simulated water-level declines in the deep (Fletcher Bay) aquifer between predevelopment and 2008 conditions ranged from about 10 feet in the northeast to about 25 feet on the western edge of the Island. These declines are related to groundwater withdrawals for public-supply purposes. The calibrated model also was used to simulate the possible effects of increased groundwater pumping and changes to recharge due to changes in land use and climactic conditions between 2008 and 2035 under minimal, expected, and maximum impact conditions. Drawdowns generally were small for most of the Island (less than 10 ft) for the minimal and expected impact scenarios, and were larger for the maximum impact scenario. No saltwater intrusion was evident in any scenario by the year 2035. The direction of flow in the deep Fletcher Bay aquifer was simulated to reverse direction from its predevelopment west to east direction to an east to west direction under the maximum impact scenario.

Scientific Investigations Report↗
USGS2012

Flooded area and plant zonation in isolated wetlands in well fields in the Northern Tampa Bay Region, Florida, following reductions in groundwater-withdrawal rates

The extent and duration of the flooded area were compared in two reference wetlands and nine wetlands in well fields in the northern Tampa Bay region, Florida, to determine whether reductions in well-field groundwater-withdrawal rates resulted in increases in wetland flooded area. Flooded area, expressed as a percentage of the total wetland area, was used to provide a quantitative and comparable line of evidence for describing the hydrologic conditions in isolated wetlands of different sizes and locations. Flooded-area frequencies were quantified for periods with different groundwater-withdrawal rates that bracket reductions in well-field groundwater withdrawals. Four-year pre-reduction and post-reduction periods were applied to wetlands in Cypress Creek and Cross Bar Ranch well fields, whereas 3-year periods were applied to wetlands in Starkey well field. The reduced groundwater-withdrawal rates in Cypress Creek and Cross Bar Ranch well fields were 30 and 24 percent less than their pre-reduction rates, respectively. The reduced groundwater-withdrawal rate in the Starkey well field was 64 percent less. Total rainfall amounts were similar (differed by 1 percent or less) in the respective pre- and post-reduction periods, which minimized the effect that rainfall variability had on the analysis. Flooded-area patterns at the reference wetlands, which were unaffected by groundwater withdrawals, were similar during pre- and post-reduction periods, indicating that short-term rainfall variability within those periods did not affect the longer-term patterns of flooded-area extent and duration. One well-field wetland (W-33) experienced an extent and duration of flooded area similar to that observed at the reference wetlands. About 61–100 percent of W-33 was flooded 41 percent of the time during the pre-reduction period and 45 percent of the time in the post-reduction period. The amount of time the wetland was dry decreased from 40 percent in the pre-reduction period to 26 percent in the post-reduction period. The median elevation of the potentiometric surface of the Upper Floridan aquifer increased beneath this wetland by about 4 feet after reductions in groundwater-withdrawal rates. Four well-field wetlands (W-17, W-56, Starkey N, and Starkey 108) had substantial increases in the extent and duration of the flooded area after reductions in groundwater-withdrawal rates. These four wetlands were dry for 25–45 percent less time during the post-reduction period, when the pre- and post-reduction periods were compared. Up to 20 percent of W-56 was flooded more than three times as long after reductions in groundwater-withdrawal rates. All parts of W-17 were flooded for as much as 10 percent of the time in the post-reduction period. Parts of Starkey N and Starkey 108 were flooded for more than twice as much time after reductions in groundwater-withdrawal rates. The median elevation of the potentiometric surface of the Upper Floridan aquifer was about 4–8 feet higher beneath W-17 and W-56 after reductions in groundwater-withdrawal rates, whereas the median elevation increased beneath Starkey N and Starkey 108 by about 4 feet after reductions in groundwater-withdrawal rates. Four other well-field wetlands (W-41, Q-1, Starkey D, and Starkey E) were mostly dry before reductions in groundwater-withdrawal rates and remained mostly dry after the reductions. W-41 was dry 23 percent less time in the post-reduction period, but most of the increase in flooded area was confined to less than 20 percent of the total wetland area. Q-1 was dry for only 12 percent less time in the post-reduction period. The median elevation of the potentiometric surface of the Upper Floridan aquifer increased beneath W-41 by about 5 feet and beneath Q-1 by about 2 feet after reductions in groundwater-withdrawal rates. The extent and duration of the flooded area was unchanged at Starkey D when the post-reduction period was compared to the pre-reduction period. At Starkey E the extent of the flooded area decreased slightly during the post-reduction period. Even though groundwater-withdrawal rates at Starkey well field decreased in the post-reduction period, the median elevation of the potentiometric surface of the Upper Floridan aquifer did not increase beneath Starkey D and Starkey E after reductions in groundwater-withdrawal rates from this well field. Factors such as the high permeability of sediments beneath the wetlands, subsidence, or sinkholes could contribute to continued downward leakage from these four wetlands and the lack of recovery of wetland water levels. Plant zonation in the two reference wetlands and the nine well-field wetlands was described using data collected by the Southwest Florida Water Management District and Tampa Bay Water, a regional utility, in their Wetland Assessment Procedure (WAP). A scoring system was used to describe the distribution of trees, woody shrubs, and groundcover in zones at three depths along a transect line through each wetland. The locations of the three zones were identified on contoured wetland bathymetry maps and were discussed in relation to areas of the wetland bottom that flooded for different periods of time during the study. Higher scores are characteristic of a greater extent and duration of wetland flooded area. WAP scores and weighted average scores for wetland vegetation were generally consistent with the results of the flooded area analysis. The WAP scores and weighted average scores were higher overall and did not decline with time at four wetlands in well fields (W-33, W-56, Starkey N, and Starkey 108) during the years following reductions in groundwater-withdrawal rates. These four wetlands also had increases in the extent and duration of the flooded area during the post-reduction period. Scores for trees were more consistent than scores for shrubs and groundcover. WAP scores remained relatively low or generally declined at five well-field wetlands (Q-1, W-17, W-41, Starkey D, and Starkey E) during the years following reductions in groundwater-withdrawal rates, and weighted average scores either declined over time or remained low. These five wetlands either did not have an increase in the extent and duration of the flooded area, or if there was an increase, it was small.

Florida↗