USGS Science⌕ Search

Geology topics

Groundwater & aquifers

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

At least 469 records · Page 26Linked to original sources

Aquifer geochemistry at potential aquifer storage and recovery sites in coastal plain aquifers in the New York city area, USA

The effects of injecting oxic water from the New York city (NYC) drinking-water supply and distribution system into a nearby anoxic coastal plain aquifer for later recovery during periods of water shortage (aquifer storage and recovery, or ASR) were simulated by a 3-dimensional, reactive-solute transport model. The Cretaceous aquifer system in the NYC area of New York and New Jersey, USA contains pyrite, goethite, locally occurring siderite, lignite, and locally varying amounts of dissolved Fe and salinity. Sediment from cores drilled on Staten Island and western Long Island had high extractable concentrations of Fe, Mn, and acid volatile sulfides (AVS) plus chromium-reducible sulfides (CRS) and low concentrations of As, Pb, Cd, Cr, Cu and U. Similarly, water samples from the Lloyd aquifer (Cretaceous) in western Long Island generally contained high concentrations of Fe and Mn and low concentrations of other trace elements such as As, Pb, Cd, Cr, Cu and U, all of which were below US Environmental Protection Agency (USEPA) and NY maximum contaminant levels (MCLs). In such aquifer settings, ASR operations can be complicated by the oxidative dissolution of pyrite, low pH, and high concentrations of dissolved Fe in extracted water.The simulated injection of buffered, oxic city water into a hypothetical ASR well increased the hydraulic head at the well, displaced the ambient groundwater, and formed a spheroid of injected water with lower concentrations of Fe, Mn and major ions in water surrounding the ASR well, than in ambient water. Both the dissolved O2 concentrations and the pH of water near the well generally increased in magnitude during the simulated 5-a injection phase. The resultant oxidation of Fe2+ and attendant precipitation of goethite during injection provided a substrate for sorption of dissolved Fe during the 8-a extraction phase. The baseline scenario with a low (0.001M) concentration of pyrite in aquifer sediments, indicated that nearly 190% more water with acceptably low concentrations of dissolved Fe could be extracted than was injected. Scenarios with larger amounts of pyrite in aquifer sediments generally resulted in less goethite precipitation, increased acidity, and increased concentrations of dissolved Fe in extracted water. In these pyritic scenarios, the lower amounts of goethite precipitated and the lower pH during the extraction phase resulted in decreased sorption of Fe2+ and a decreased amount of extractable water with acceptably low concentrations of dissolved Fe (5.4??10-6M). A linear decrease in recovery efficiency with respect to dissolved Fe concentrations is caused by pyrite dissolution and the associated depletion of dissolved O2 (DO) and increase in acidity. Simulations with more than 0.0037M of pyrite, which is the maximum amount dissolved in the baseline scenario, had just over a 50% recovery efficiency. The precipitation of ferric hydroxide minerals (goethite) at the well screen, and a possible associated decrease in specific capacity of the ASR well, was not apparent during the extraction phase of ASR simulations, but the model does not incorporate the microbial effects and biofouling associated with ferric hydroxide precipitation.The host groundwater chemistry in calcite-poor Cretaceous aquifers of the NYC area consists of low alkalinity and moderate to low pH. The dissolution of goethite in scenarios with unbuffered injectate indicates that corrosion of the well could occur if the injectate is not buffered. Simulations with buffered injectate resulted in greater precipitation of goethite, and lower concentrations of dissolved Fe, in the extracted water. Dissolved Fe concentrations in extracted water were highest in simulations of aquifers (1) in which pyrite and siderite in the aquifer were in equilibrium, and (2) in coastal areas affected by saltwater intrusion, where high dissolved-cation concentrations provide a greater exchange of Fe2+ (FeX2). Results indicate that ASR in pyrite-beari

Applied Geochemistry↗

Influence of groundwater recharge and well characteristics on dissolved arsenic concentrations in southeastern Michigan groundwater

Arsenic concentrations exceeding 10 ??g/l, the United States maximum contaminant level and the World Health Organization guideline value, are frequently reported in groundwater from bedrock and unconsolidated aquifers of southeastern Michigan. Although arsenic-bearing minerals (including arsenian pyrite and oxide/hydroxide phases) have been identified in Marshall Sandstone bedrock of the Mississippian aquifer system and in tills of the unconsolidated aquifer system, mechanisms responsible for arsenic mobilization and subsequent transport in groundwater are equivocal. Recent evidence has begun to suggest that groundwater recharge and characteristics of well construction may affect arsenic mobilization and transport. Therefore, we investigated the relationship between dissolved arsenic concentrations, reported groundwater recharge rates, well construction characteristics, and geology in unconsolidated and bedrock aquifers. Results of multiple linear regression analyses indicate that arsenic contamination is more prevalent in bedrock wells that are cased in proximity to the bedrock-unconsolidated interface; no other factors were associated with arsenic contamination in water drawn from bedrock or unconsolidated aquifers. Conditions appropriate for arsenic mobilization may be found along the bedrock-unconsolidated interface, including changes in reduction/oxidation potential and enhanced biogeochemical activity because of differences between geologic strata. These results are valuable for understanding arsenic mobilization and guiding well construction practices in southeastern Michigan, and may also provide insights for other regions faced with groundwater arsenic contamination. ?? Springer-Verlag 2008.

Environmental Geochemistry and Health↗

Processes affecting geochemistry and contaminant movement in the middle Claiborne aquifer of the Mississippi embayment aquifer system

Groundwater chemistry and tracer-based age data were used to assess contaminant movement and geochemical processes in the middle Claiborne aquifer (MCA) of the Mississippi embayment aquifer system. Water samples were collected from 30 drinking-water wells (mostly domestic and public supply) and analyzed for nutrients, major ions, pesticides, volatile organic compounds (VOCs), and transient age tracers (chlorofluorocarbons, tritium and helium-3, and sulfur hexafluoride). Redox conditions are highly variable throughout the MCA. However, mostly oxic groundwater with low dissolved solids is more vulnerable to nitrate contamination in the outcrop areas east of the Mississippi River in Mississippi and west Tennessee than in mostly anoxic groundwater in downgradient areas in western parts of the study area. Groundwater in the outcrop area was relatively young (apparent age of less than 40 years) with significantly ( p < 0.05) higher dissolved oxygen and nitrate&ndash;N concentrations and higher detections of pesticides and VOCs compared to water samples from wells in downgradient areas. Oxygen reduction and denitrification rates were low compared to other aquifers in the United States (zero order rate constants for oxygen reduction and denitrification were 4.7 and 5&ndash;10 &mu;mol/L/year, respectively). Elevated concentrations of nitrate&ndash;N, and detections of pesticides and VOCs in some deep public supply wells (>50 m depth) indicated contaminant movement from shallow parts of the aquifer into deeper oxic zones. Given the persistence of nitrate in young oxic groundwater that was recharged several decades ago, and the lack of a confining unit, the downward movement of young contaminated water may result in higher nitrate concentrations over time in deeper parts of the aquifer containing older oxic water.

Claiborne Aquifer;Mississippi Embayment↗

Simulation of groundwater flow, 1895–2010, and effects of additional groundwater withdrawals on future stream base flow in the Elkhorn and Loup River Basins, central Nebraska—Phase three

The U.S. Geological Survey, in cooperation with the Lewis and Clark, Lower Elkhorn, Lower Loup, Lower Platte North, Lower Niobrara, Middle Niobrara, Upper Elkhorn, and the Upper Loup Natural Resources Districts, designed a study to refine the spatial and temporal discretization of a previously modeled area. This updated study focused on a 30,000-square-mile area of the High Plains aquifer and constructed regional groundwater-flow models to evaluate the effects of groundwater withdrawal on stream base flow in the Elkhorn and Loup River Basins, Nebraska. The model was calibrated to match groundwater-level and base-flow data from the stream-aquifer system from pre-1940 through 2010 (including predevelopment [pre-1895], early development [1895–1940], and historical development [1940 through 2010] conditions) using an automated parameter-estimation method. The calibrated model then was used to simulate hypothetical development conditions (2011 through 2060). Predicted changes to stream base flow based on simulated changes to groundwater withdrawal will aid in developing strategies for management of hydrologically connected water supplies. Additional wells were simulated throughout the model domain and pumped for 50 years to assess the effect of wells on aquifer depletions, including stream base flow. The percentage of withdrawal for each well after 50 years, which was compensated by aquifer reductions to stream base flow, storage, or evapotranspiration, was computed and mapped. These depletions are influenced by aquifer properties, time, and distance from the well. Stream base-flow depletion results showed that the closer the added well was to a stream, the greatest the effect on the stream base flow. Areas of stream base-flow depletion percentages greater than 80 percent were generally within 1 mile (mi) from the stream. The distance increased to 6 mi near the confluence of the Dismal and Middle Loup Rivers, and the North Loup and Calamus Rivers. The percentage of stream base-flow depletion decreased as the distance from the stream increased. Areas more than 10 mi from the stream generally had a stream base-flow depletion of 10 percent or less. Evapotranspiration depletion was largest in areas closest to streams, specifically in the Elkhorn River watershed. It was also larger in areas of interdunal wetlands within the Sand Hills. Evapotranspiration depletion was negligible in areas greater than 5 mi from a stream, with the exception of interdunal areas in Cherry, Grant, and Arthur Counties. The storage depletion percentage increased as the distance from a stream increased. Storage depletion was largest in areas between streams. Areas experiencing the smallest amount of storage depletion were adjacent to streams. Calibrated model outputs and streamflow depletion analysis are publicly available online. Accuracy of the simulations is affected by input data limitations, system simplifications, assumptions, and resources available at the time of the simulation construction and calibration. Most of the important limitations relate either to data used as simulation inputs or to data used to estimate simulation inputs. Development of the regional simulations focused on generalized hydrogeologic characteristics within the study area and did not attempt to describe variations important to local-scale conditions. These simulations are most appropriate for analyzing groundwater-management scenarios for large areas and during long periods and are not suitable for analysis of small areas or short periods.

Nebraska↗

Groundwater – The disregarded component in lake water and nutrient budgets. Part 2: effects of groundwater on nutrients

Lacustrine groundwater discharge (LGD) transports nutrients from a catchment to a lake, which may fuel eutrophication, one of the major threats to our fresh waters. Unfortunately, LGD has often been disregarded in lake nutrient studies. Most measurement techniques are based on separate determinations of volume and nutrient concentration of LGD: Loads are calculated by multiplying seepage volumes by concentrations of exfiltrating water. Typically low phosphorus (P) concentrations of pristine groundwater often are increased due to anthropogenic sources such as fertilizer, manure or sewage. Mineralization of naturally present organic matter might also increase groundwater P. Reducing redox conditions favour P transport through the aquifer to the reactive aquifer-lake interface. In some cases, large decreases of P concentrations may occur at the interface, for example, due to increased oxygen availability, while in other cases, there is nearly no decrease in P. The high reactivity of the interface complicates quantification of groundwater-borne P loads to the lake, making difficult clear differentiation of internal and external P loads to surface water. Anthropogenic sources of nitrogen (N) in groundwater are similar to those of phosphate. However, the environmental fate of N differs fundamentally from P because N occurs in several different redox states, each with different mobility. While nitrate behaves essentially conservatively in most oxic aquifers, ammonium's mobility is similar to that of phosphate. Nitrate may be transformed to gaseous N 2 in reducing conditions and permanently removed from the system. Biogeochemical turnover of N is common at the reactive aquifer-lake interface. Nutrient loads from LGD were compiled from the literature. Groundwater-borne P loads vary from 0.74 to 2900 mg PO 4 -P m −2 year −1 ; for N, these loads vary from 0.001 to 640 g m −2 year −1 . Even small amounts of seepage can carry large nutrient loads due to often high nutrient concentrations in groundwater. Large spatial heterogeneity, uncertain areal extent of the interface and difficult accessibility make every determination of LGD a challenge. However, determinations of LGD are essential to effective lake management.

Hydrological Processes↗

Groundwater quality and groundwater levels in Dougherty County, Georgia, April 2020 through January 2023

The Upper Floridan aquifer is the uppermost reliable groundwater source in southwest Georgia. The aquifer lies on top of the Claiborne, Clayton, and Cretaceous aquifers, all of which exhibited water-level declines in the 1960s and 1970s. The U.S. Geological Survey has been working cooperatively with Albany Utilities to monitor groundwater quality and availability in these aquifers since 1977. Flow direction in the Upper Floridan aquifer is to the south and toward the Flint River. During the past 3 years, water levels varied above and below period-of-record median values. Water levels in the Upper Floridan aquifer were primarily above or at median levels during 2020 and 2021 and at or below median levels during 2022. Water levels in the Claiborne aquifer were above median levels, whereas water levels in the Clayton aquifer were at or below median levels, and in the Cretaceous aquifer system were close to median levels. During January 2021, eight wells were sampled for major ions, including nitrate plus nitrite as nitrogen (N). Nitrate plus nitrite as N concentrations ranged from 2.3 to 10.5 milligrams per liter (mg/L). During December 2021, seven wells were sampled for major ions, including nitrate plus nitrite as N. Nitrate plus nitrite as N concentrations ranged from 3.9 to 9.9 mg/L. During November 2022, eight wells were sampled for major ions, including nitrate plus nitrite as N. Nitrate plus nitrite as N concentrations ranged from 3.9 to 10.0 mg/L. Two wells were also sampled for per- and polyfluoroalkyl substances during November 2022.

Georgia↗

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

Most of the subsidence in the Houston-Galveston region has occurred as a direct result of groundwater withdrawals for municipal supply, industrial use, and irrigation that depressured and dewatered the Chicot and Evangeline aquifers causing compaction of the clay layers of the aquifer sediments. This report, prepared by the U.S. Geological Survey, in cooperation with the Harris-Galveston Subsidence District, City of Houston, Fort Bend Subsidence District, and Lone Star Groundwater Conservation District, is one in an annual series of reports depicting water-level altitudes and water-level changes in the Chicot, Evangeline, and Jasper aquifers and compaction in the Chicot and Evangeline aquifers in the Houston-Galveston region. The report contains maps showing 2010 water-level altitudes for the Chicot, Evangeline, and Jasper aquifers, respectively; maps showing 1-year (2009-10) water-level-altitude changes for each aquifer; maps showing 5-year (2005-10) water-level-altitude changes for each aquifer; maps showing long-term (1990-2010 and 1977-2010) water-level-altitude changes for the Chicot and Evangeline aquifers; a map showing long-term (2000-10) water-level-altitude change for the Jasper aquifer; a map showing locations of borehole extensometer sites; and graphs showing measured compaction of subsurface material at the extensometers from 1973, or later, through 2009. Tables listing the data used to construct each aquifer-data map and the compaction graphs are included. Water levels in the Chicot, Evangeline, and Jasper aquifers were measured during December 2009-March 2010. In 2010, water-level-altitude contours for the Chicot aquifer ranged from 200 feet below National Geodetic Vertical Datum of 1929 or North American Vertical Datum of 1988 (hereinafter, datum) in a small area in southwestern Harris County to 200 feet above datum in central to southwestern Montgomery County. Water-level-altitude changes in the Chicot aquifer ranged from a 49-foot decline to a 67-foot rise (2009-10), from a 25-foot decline to a 35-foot rise (2005-10), from a 40-foot decline to an 80-foot rise (1990-2010), and from a 140-foot decline to a 200-foot rise (1977-2010). In 2010, water-level-altitude contours for the Evangeline aquifer ranged from 300 feet below datum in north-central Harris County to 200 feet above datum at the boundary of Waller, Montgomery, and Grimes Counties. Water-level-altitude changes in the Evangeline aquifer ranged from a 58-foot decline to a 69-foot rise (2009-10), from an 80-foot decline to an 80-foot rise (2005-10), from a 200-foot decline to a 220-foot rise (1990-2010), and from a 320-foot decline to a 220-foot rise (1977-2010). In 2010, water-level-altitude contours for the Jasper aquifer ranged from 200 feet below datum in south-central Montgomery County to 250 feet above datum in eastern-central Grimes County. Water-level-altitude changes in the Jasper aquifer ranged from a 39-foot decline to a 39-foot rise (2009-10), from a 110-foot decline to no change (2005-10), and from a 180-foot decline to no change (2000-10). Compaction of subsurface materials (mostly in the clay layers) composing the Chicot and Evangeline aquifers was recorded continuously at 13 borehole extensometers at 11 sites. For the period of record beginning in 1973, or later, and ending in December 2009, cumulative clay compaction data measured by 12 extensometers ranged from 0.088 foot at the Texas City-Moses Lake site to 3.559 foot at the Addicks site. The rate of compaction varies from site to site because of differences in groundwater withdrawals near each site and differences among sites in the clay-to-sand ratio in the subsurface materials. Therefore, it is not possible to extrapolate or infer a rate of clay compaction for an area based on the rate of compaction measured at a nearby extensometer.

Texas↗

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

Most of the subsidence in the Houston–Galveston region has occurred as a direct result of groundwater withdrawals for municipal supply, industrial use, and irrigation that depressured and dewatered the Chicot and Evangeline aquifers causing compaction of the clay layers of the aquifer sediments. This report, prepared by the U.S. Geological Survey, in cooperation with the Harris–Galveston Subsidence District, City of Houston, Fort Bend Subsidence District, and Lone Star Groundwater Conservation District, is one in an annual series of reports depicting water-level altitudes and water-level changes in the Chicot, Evangeline, and Jasper aquifers and compaction in the Chicot and Evangeline aquifers in the Houston–Galveston region. The report contains maps showing 2011 water-level altitudes for the Chicot, Evangeline, and Jasper aquifers; maps showing 1-year (2010–11) water-level-altitude changes for each aquifer; maps showing 5-year (2006–11) water-level-altitude changes for each aquifer; maps showing long-term (1990–2011 and 1977–2011) water-level-altitude changes for the Chicot and Evangeline aquifers; a map showing long-term (2000–11) water-level-altitude change for the Jasper aquifer; a map showing locations of borehole extensometer sites; and graphs showing measured compaction of subsurface material at the extensometers from 1973, or later, through 2010. Tables listing the data used to construct each aquifer-data map and the compaction graphs are included. Water levels in the Chicot, Evangeline, and Jasper aquifers were measured during December 2010–February 2011. In 2011, water-level-altitude contours for the Chicot aquifer ranged from 200 feet below North American Vertical Datum of 1988 (hereinafter, datum) in a small area in southwestern Harris County to 200 feet above datum in central to southwestern Montgomery County. Water-level-altitude changes in the Chicot aquifer ranged from a 40-foot decline to a 33-foot rise (2010–11), from a 10-foot decline to an 80-foot rise (2006–11), from a 140-foot decline to a 100-foot rise (1990–2011), and from a 120-foot decline to a 200-foot rise (1977–2011). In 2011, water-level-altitude contours for the Evangeline aquifer ranged from 300 feet below datum in north-central Harris County to 200 feet above datum at the boundary of Waller, Montgomery, and Grimes Counties. Water-level-altitude changes in the Evangeline aquifer ranged from a 43-foot decline to a 73-foot rise (2010–11), from a 40-foot decline to a 160-foot rise (2006–11), from a 200-foot decline to a 240-foot rise (1990–2011), and from a 340-foot decline to a 260-foot rise (1977–2011). In 2011, water-level-altitude contours for the Jasper aquifer ranged from 200 feet below datum in south-central Montgomery County to 250 feet above datum in east-central Grimes County. Water-level-altitude changes in the Jasper aquifer ranged from a 45-foot decline to a 29-foot rise (2010–11), from a 90-foot decline to a 10-foot rise (2006–11), and from a 190-foot decline to no change (2000–11). Compaction of subsurface materials (mostly in the clay layers) composing the Chicot and Evangeline aquifers was recorded continuously at 13 borehole extensometers at 11 sites. For the period of record beginning in 1973, or later, and ending in December 2010, cumulative clay compaction data measured by 12 extensometers ranged from 0.100 foot at the Texas City–Moses Lake site to 3.544 foot at the Addicks site. The rate of compaction varies from site to site because of differences in groundwater withdrawals near each site and differences among sites in the clay-to-sand ratio in the subsurface materials. Therefore, it is not possible to extrapolate or infer a rate of clay compaction for an area based on the rate of compaction measured at a nearby extensometer.

Texas↗

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

The U.S. Geological Survey, in cooperation with the Arkansas Natural Resources Commission and the Arkansas Geological Survey, has monitored water levels in the Sparta Sand of Claiborne Group and Memphis Sand of Claiborne Group (herein referred to as “the Sparta Sand” and “the Memphis Sand,” respectively) since the 1920s. Groundwater withdrawals have increased while water levels have declined since monitoring was initiated. Herein, aquifers in the Sparta Sand and Memphis Sand will be referred to as “the Sparta-Memphis aquifer” throughout Arkansas. During the spring of 2011, 291 water levels were measured in wells completed in the Sparta-Memphis aquifer and used to produce a regional potentiometric-surface map. During the summer of 2011, groundwater-quality samples were collected and measured from 61 wells for specific conductance, pH, and temperature. In the northern half of Arkansas, the regional direction of groundwater flow in the Sparta-Memphis aquifer is generally to the south-southeast and flows east and south in the southern half of Arkansas. The groundwater in the southern half of Arkansas flows away from the outcrop area except where affected by large depressions in the potentiometric surface. The highest and lowest water-level altitudes measured in the Sparta-Memphis aquifer were 326 feet above and 120 feet below National Geodetic Vertical Datum of 1929 (NGVD 29), respectively. Five depressions are located in the following counties: Arkansas, Cleveland, Jefferson, Lincoln, and Prairie; Union; Cross, Poinsett, St. Francis, and Woodruff; Columbia; and Bradley. Two large depressions, centered in Jefferson and Union Counties, are the result of large withdrawals for industrial, irrigation, or public supply. The depression centered in Jefferson County has expanded in recent years into Arkansas and Prairie Counties as a result of large withdrawals for irrigation and public supply. The lowest water-level altitude measured in this depression is approximately 20 feet (ft) higher in 2011 than in 2009. The area enclosed within the 40-ft contour on the 2011 potentiometric-surface map has decreased in area, shifting north in Lincoln County and west in Arkansas County when compared with the 2009 potentiometric-surface map. The depression in Union County is roughly circular within the -60-ft contour. The lowest water-level altitude measurement was 157 ft below NGVD 29 in 2009, with a 37-ft rise to 120 ft below NGVD 29 in 2011. The depression in Union County has diminished and encloses a smaller area than in recent years. In 1993, the -60-ft contour enclosed 632 square miles (mi 2 ). In 2011, the -60-ft contour enclosed 375 mi 2 , a decrease of 41 percent from 1993. The lowest water-level altitude measurement during 2011 in the center of the depression in Union County represents a rise of 79 ft since 2003. The area enclosed by the lowest altitude contour, 120 ft below NGVD 29, on the 2011 potentiometric-surface map is less than 10 percent of the area enclosed by that same contour on the 2009 potentiometric-surface map. A broad depression in western Poinsett and Cross Counties was first shown in the 1995 potentiometric-surface map. In 2011, the lowest water-level altitude measurement in this depression, 129 ft above NGVD 29, is 2 ft lower than in 2009. The 140-ft contour has extended southwest into northwestern St. Francis and east-central Woodruff Counties in 2011. In Columbia County in 2011, the area of the depression has decreased, with water levels rising about 1 ft since 2005 in the well with the lowest water-level altitude measurement. The depression in Bradley County in 2011 has decreased in area compared to 2007. A water-level difference map was constructed using the difference between water-level measurements made during 2007 and 2011 at 247 wells. The differences in water level between 2007 and 2011 ranged from -17.3 to 45.4 ft, with a mean of 4.1 ft. Water levels generally declined in the northern half of the study area and generally increased in the southern half of the study area. Areas with a general decline in water levels include Lonoke and western Prairie Counties; northern Arkansas County; Miller County; and Craighead, Poinsett, Cross, and Woodruff Counties. Areas with a general rise in water levels include Lafayette, Columbia, Union, Calhoun, and Bradley Counties; Grant, Jefferson, southern Arkansas, Lincoln, Drew, and Desha Counties; and Phillips County. Hydrographs from 183 wells with a minimum of 25 years of water-level measurements were constructed. During the period 1987–2011, county mean annual water levels generally declined. Mean annual declines were between 0.5 foot per year (ft/yr) and 0.0 ft/yr in Ashley, Chicot, Crittenden, Drew, Grant, Jefferson, Lafayette, Mississippi, Monroe, Ouachita, Phillips, Pulaski, St. Francis, and Woodruff Counties. Mean annual declines were between 1.0 ft/yr and 0.5 ft/yr in Bradley, Calhoun, Cleveland, Craighead, Cross, Desha, Lonoke, Miller, Poinsett, and Prairie Counties. Mean annual declines were between 1.5 ft/yr and 1.0 ft/yr in Arkansas, Lee, and Lincoln Counties. The county mean annual water level rose in Columbia, Dallas, and Union Counties about 0.3 ft/yr, 0.1 ft/yr, and 1.2 ft/yr, respectively. Water samples were collected in the summer of 2011 from 61 wells completed in the Sparta-Memphis aquifer and measured onsite for specific conductance, temperature, and pH. Although there is a regional increase in specific conductance to the east and south, anomalous increases occur in some parts of the study area. Specific conductance ranged from 35 microsiemens per centimeter (μS/cm) in Ouachita County to 1,380 μS/cm in Monroe County. Relatively large specific conductance values (greater than 700 mS/cm) occur in samples from wells in Arkansas, Ashley, Clay, Monroe, Phillips, and Union Counties.

Arkansas↗

Simulation of groundwater flow and interaction of groundwater and surface water on the Lac du Flambeau Reservation, Wisconsin

The Lac du Flambeau Band of Lake Superior Chippewa and Indian Health Service are interested in improving the understanding of groundwater flow and groundwater/surface-water interaction on the Lac du Flambeau Reservation (Reservation) in southwest Vilas County and southeast Iron County, Wisconsin, with particular interest in an understanding of the potential for contamination of groundwater supply wells and the fate of wastewater that is infiltrated from treatment lagoons on the Reservation. This report describes the construction, calibration, and application of a regional groundwater flow model used to simulate the shallow groundwater flow system of the Reservation and water-quality results for groundwater and surface-water samples collected near a system of waste-water-treatment lagoons. Groundwater flows through a permeable glacial aquifer that ranges in thickness from 60 to more than 200 feet (ft). Seepage and drainage lakes are common in the area and influence groundwater flow patterns on the Reservation. A two-dimensional, steady-state analytic element groundwater flow model was constructed using the program GFLOW. The model was calibrated by matching target water levels and stream base flows through the use of the parameter-estimation program, PEST. Simulated results illustrate that groundwater flow within most of the Reservation is toward the Bear River and the chain of lakes that feed the Bear River. Results of analyses of groundwater and surface-water samples collected downgradient from the wastewater infiltration lagoons show elevated levels of ammonia and dissolved phosphorus. In addition, wastewater indicator chemicals detected in three downgradient wells and a small downgradient stream indicate that infiltrated wastewater is moving southwest of the lagoons toward Moss Lake. Potential effects of extended wet and dry periods (within historical ranges) were evaluated by adjusting precipitation and groundwater recharge in the model and comparing the resulting simulated lake stage and water budgets to stages and water budgets from the calibrated model. Simulated lake water budgets and water level changes illustrate the importance of understanding the position of a lake within the hydrologic system (headwater or downstream), the type of lake (surface-water drainage or seepage lake), and the role of groundwater in dampening the effects of large-scale changes in weather patterns on lake levels. Areas contributing recharge to drinking-water supply wells on the Reservation were delineated using forward particle tracking from the water table to the well. Monte Carlo uncertainty analyses were used to produce maps showing the probability of groundwater capture for areas around each well nest. At the Main Pumphouse site near the Village of Lac du Flambeau, most of the area contributing recharge to the wells occurs downgradient from a large wetland between the wells and the wastewater infiltration lagoons. Nonetheless, a small potential for the wells to capture infiltrated wastewater is apparent when considering uncertainty in the model parameter values. At the West Pumphouse wells south of Flambeau Lake, most of the area contributing recharge is between the wells and Tippecanoe Lake. The extent of infiltrated wastewater from two infiltration lagoons was tracked using the groundwater flow model and Monte Carlo uncertainty analyses. Wastewater infiltrated from the lagoons flows predominantly south toward Moss Lake as it integrates with the regional groundwater flow system. The wastewater-plume-extent simulations support the area-contributing-recharge simulations, indicating that there is a possibility, albeit at low probability, that some wastewater could be captured by water-supply wells. Comparison of simulated water-table contours indicate that the lagoons may mound the water table approximately 4 ft, with diminishing levels of mounding outward from the lagoons. Four scenarios, representing potential alternatives for wastewater management, were simulated (at current discharge rates) to evaluate the potential extent of wastewater in the aquifer and discharge to surface-water bodies associated with each management scenario. Wastewater simulated to infiltrate through a hypothetical diffuser below a wetland south of the current lagoons appears to discharge to the overlying wetland and would likely discharge to Moss Lake as overland flow. Wastewater simulated to discharge to a small lake (Mindy Lake) between Moss and Fence Lakes appears to spread radically over a large area between the lakes. Wastewater simulated to discharge to lagoons south and northeast of the current lagoons also appears to spread radially, but the areas of the aquifer with the highest probability of encountering waste-water contamination would likely be between the lagoons and the nearest lake, where the wastewater would eventually discharge. Probability results for the wastewater-plume-extent scenarios are sensitive to the number of mathematical water particles used to represent infiltrating wastewater and the level of detail in the synthetic grid used for the probability analysis. Thus, probability results from wastewater-plume-extent simulations are qualitative only; however, it is expected that illustrations of relatively high or low probability will be useful as a general guide for decision making. Management problems requiring quantitative estimates of probability are best re-cast into problems evaluating the area that contributes recharge to the location of interest, which is not dependent upon the number of simulated particles or the resolution of a synthetic grid.

Wisconsin↗

Groundwater flow conceptualization of the Pahute Mesa–Oasis Valley Groundwater Basin, Nevada—A synthesis of geologic, hydrologic, hydraulic-property, and tritium data

This report provides a groundwater-flow conceptualization that integrates geologic, hydrologic, hydraulic-property, and radionuclide data in the Pahute Mesa–Oasis Valley (PMOV) groundwater basin, southern Nevada. Groundwater flow in the PMOV basin is of interest because 82 underground nuclear tests were detonated, most near or below the water table. A potentiometric map and nine sets of hydrostratigraphic and hydrologic cross sections supplement the conceptualization. Potentiometric contours indicate that groundwater in the PMOV basin generally flows south-southwest and discharges at Oasis Valley. Groundwater encounters an alternating sequence of low- and high-transmissivity rocks, referred to as dams and pools, respectively, as it moves from east to west across eastern Pahute Mesa. Flow from all Pahute Mesa nuclear tests is to Oasis Valley and is well-constrained by water-level data. Flow converges along a corridor of high transmissivity between Pahute Mesa and Oasis Valley. The location of the lateral PMOV basin boundary is well defined, and this boundary, with a few minor exceptions, represents a no-flow boundary. Some boundary uncertainty exists in the northeastern part of the basin, but potential flow-rate estimates across the northeastern boundary resulting from this uncertainty are small relative to the basin groundwater budget. Recharge in the PMOV basin is derived from episodic pulses of modern water and the diffuse percolation of old water (greater than 1,000 years). Episodic recharge is a minor recharge component observed as a rise in groundwater levels that occurs 3 months to 1 year following a wet winter. Minor amounts of episodic recharge through an unsaturated zone in excess of 1,000 feet (ft) requires preferential flow through faults and fractures. The dominant recharge component is slow, steady, diffuse percolation of old water through the unsaturated zone. A large component of old water recharging the groundwater system is consistent with observations of isotopically light deuterium and oxygen 18 compositions in water from wells on Pahute Mesa and central Oasis Valley. About half the recharge in the PMOV basin is derived from the eastern Pahute Mesa area. The remaining recharge is derived primarily from other highland areas including Timber Mountain, Belted and Kawich Ranges, and Black Mountain. The PMOV groundwater system is nearly steady state, where recharge is balanced by the 5,900 acre-feet per year of natural discharge at Oasis Valley. This assumption is reasonable because the basin is dominated by steady-state conditions, where long-term changes in groundwater storage are minimal. Total groundwater withdrawals from 1963 to 2018 have amounted to less than 10 percent of annual groundwater discharge and less than 0.2 percent of the basin’s groundwater storage. Therefore, present-day (2020) conditions are considered representative of predevelopment (pre-1950) conditions in nearly all areas of the basin. The lower PMOV basin boundary is defined at 4,000 ft below the water table to encompass all underground nuclear tests and tritium plumes. This boundary defines the lower boundary of radionuclide migration. However, nearly all flow and tritium transport occur in the upper 1,600 ft of the saturated zone because a transmissivity-with-depth relation indicates that greater than 90 percent of the transmissivity contributing to groundwater flow occurs within 1,600 ft of the water table. Rocks at deeper depths have low transmissivity because argillic and mineralized alterations plug the fractures. Volcanic rocks form the primary aquifers and confining units in the PMOV basin. Volcanic hydrogeologic units (HGUs) and hydrostratigraphic units (HSUs) have transmissivity distributions that span up to eight orders of magnitude with considerable overlap between distributions. Despite the large overlap between units, mean transmissivities of aquifers are one-to-two orders of magnitude greater than the confining units. However, all volcanic-rock HGUs and HSUs are composite units, meaning that they can function spatially as either an aquifer or confining unit .

Nevada↗

Use of chlorine-36 to determine regional-scale aquifer dispersivity, eastern Snake River Plain aquifer, Idaho/USA

Chlorine-36 (36Cl) derived from processed nuclear waste that was disposed at the US Department of Energy's Idaho National Engineering and Environmental Laboratory (INEEL) through a deep injection well in 1958, was detected 24-28 yr later in groundwater monitoring wells approximately 26 km downgradient from the source. Groundwater samples covering the period 1966-1995 were selected from the US Geological Survey's archived-sample library at the INEEL and analyzed for 36Cl by accelerator mass spectrometry (AMS). The smaller 36Cl peak concentrations in water from the far-field monitoring wells relative to the input suggest that aquifer dispersivity may be large. However, the sharpness of the 1958 disposal peak of 36Cl is matched by the measured 36Cl concentrations in water from these wells. This implies that a small aquifer dispersivity may be attributed to preferential groundwater flowpaths. Assuming that tracer arrival times at monitoring wells are controlled by preferential flow, a 1-D system-response model was used to estimate dispersivity by comparing the shape of predicted 36Cl-concentration curves to the shape of 36Cl-concentration curves measured in water from these observation wells. The comparisons suggest that a 1-D dispersivity of 5 m provides the best fit to the tracer data. Previous work using a 2-D equivalent porous-media model concluded that longitudinal dispersivity (equivalent to 1-D dispersivity in our model) was 90 m (Ackerman, 1991). A 90 m dispersivity value eliminates the 1958 disposal peak in our model output curves. The implications of the arrival of 36Cl at downgradient monitoring wells are important for three reasons: (1) the arrival times and associated 36Cl concentrations provide quantitative constraints on residence times, velocities, and dispersivities in the aquifer; (2) they help to refine our working hypotheses of groundwater flow in this aquifer and (3) they may suggest a means of estimating the distribution of preferential flowpaths in the aquifer. ?? 2000 Elsevier Science B.V. All rights reserved.

Nuclear Instruments and Methods in Physics Researc↗

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

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

Florida↗

Statistical analysis of lake levels and field study of groundwater and surface-water exchanges in the northeast Twin Cities Metropolitan Area, Minnesota, 2002 through 2015: Chapter A of Water levels and groundwater and surface-water exchanges in lakes of the northeast Twin Cities Metropolitan Area, Minnesota, 2002 through 2015

Water levels declined from 2003 to 2011 in many lakes in Ramsey and Washington Counties in the northeast Twin Cities Metropolitan Area, Minnesota; however, water levels in other northeast Twin Cities Metropolitan Area lakes increased during the same period. Groundwater and surface-water exchanges can be important in determining lake levels where these exchanges are an important component of the water budget of a lake. An understanding of groundwater and surface-water exchanges in the northeast Twin Cities Metropolitan Area has been limited by the lack of hydrologic data. The U.S. Geological Survey, in cooperation with the Metropolitan Council and Minnesota Department of Health, completed a field and statistical study assessing lake-water levels and regional and local groundwater and surface-water exchanges near northeast Twin Cities Metropolitan Area lakes. This report documents the analysis of collected hydrologic, water-quality, and geophysical data; and existing hydrologic and geologic data to (1) assess the effect of physical setting and climate on lake-level fluctuations of selected lakes, (2) estimate potential percentages of surface-water contributions to well water across the northeast Twin Cities Metropolitan Area, (3) estimate general ages for waters extracted from the wells, and (4) assess groundwater inflow to lakes and lake-water outflow to aquifers downgradient from White Bear Lake. Statistical analyses of lake levels during short-term (2002–10) and long-term (1925–2014) periods were completed to help understand lake-level changes across the northeast Twin Cities Metropolitan Area. Comparison of 2002–10 lake levels to several landscape and geologic characteristics explained variability in lake-level changes for 96 northeast Twin Cities Metropolitan Area lakes. Application of several statistical methods determined that (1) closed-basin lakes (without an active outlet) had larger lake-level declines than flow-through lakes with an outlet; (2) closed-basin lake-level changes reflected groundwater-level changes in the Quaternary, Prairie du Chien, and Jordan aquifers; (3) the installation of outlet-control structures, such as culverts and weirs, resulted in smaller multiyear lake-level changes than lakes without outlet-control structures; (4) water levels in lakes primarily overlying Superior Lobe deposits were significantly more variable than lakes primarily overlying Des Moines Lobe deposits; (5) lake-level declines were larger with increasing mean lake-level elevation; and (6) the frequency of some of these characteristics varies by landscape position. Flow-through lakes and lakes with outlet-control structures were more common in watersheds with more than 50 percent urban development compared to watersheds with less than 50 percent urban development. A comparison of two 35-year periods during 1925–2014 revealed that variability of annual mean lake levels in flow-through lakes increased when annual precipitation totals were more variable, whereas variability of annual mean lake levels in closed-basin lakes had the opposite pattern, being more variable when annual precipitation totals were less variable. Oxygen-18/oxygen-16 and hydrogen-2/hydrogen-1 ratios for water samples from 40 wells indicated the well water was a mixture of surface water and groundwater in 31 wells, whereas ratios from water sampled from 9 other wells indicated that water from these wells receive no surface-water contribution. Of the 31 wells with a mixture of surface water and groundwater, 11 were downgradient from White Bear Lake, likely receiving water from deeper parts of the lake. Age dating of water samples from wells indicated that the age of water in the Prairie du Chien and Jordan aquifers can vary widely across the northeast Twin Cities Metropolitan Area. Estimated ages of recharge for 9 of the 40 wells sampled for chlorofluorocarbon concentrations ranged widely from the early 1940s to mid-1970s. The wide range in estimated ages of recharge may have resulted from the wide range in the open-interval lengths and depths for the wells. Results from stable isotope analyses of water samples, lake-sediment coring, continuous seismic-reflection profiling, and water-level and flow monitoring indicated that there is groundwater inflow from nearshore sites and lake-water outflow from deep-water sites in White Bear Lake. Continuous seismic-reflection profiling indicated that deep sections of White Bear, Pleasant, Turtle, and Big Marine Lakes have few trapped gases and little organic material, which indicates where groundwater and lake-water exchanges are more likely. Water-level differences between White Bear Lake and piezometer and seepage measurements in deep waters of the lake indicate that groundwater and lake-water exchange is happening in deep waters, predominantly downgradient from the lake and into the lake sediment. Seepage fluxes measured in the nearshore sites of White Bear Lake generally were higher than seepage fluxes measured in the deep-water sites, which indicates that groundwater-inflow rates at most of the nearshore sites are higher than lake-water outflow from the deep-water sites.

Minnesota↗

Sources of high-chloride water and managed aquifer recharge in an alluvial aquifer in California, USA

As a result of pumping in excess of recharge, water levels in alluvial aquifers within the Eastern San Joaquin Groundwater Subbasin, 130 km east of San Francisco (California, USA), declined below sea level in the early 1950s and have remained so to the present. Chloride concentrations in some wells increased during that time and exceeded the US Environmental Protection Agency&rsquo;s secondary maximum contaminant level of 250 mg/L, resulting in removal of some wells from service. Sources of high-chloride water include irrigation return in 16 % of sampled wells and water from delta sediments and deeper groundwater in 50 % of sampled wells. Chloride concentrations resulting from irrigation return commonly did not exceed 100 mg/L, although nitrate concentrations were as high as 25 mg/L as nitrogen. Chloride concentrations ranged from less than 100&ndash;2,050 mg/L in wells affected by water from delta sediments and deeper groundwater. Sequential electromagnetic logs show movement of high-chloride water from delta sediments to pumping wells through permeable interconnected aquifer layers. &delta;D and &delta;18O data show most groundwater originated as recharge along the front of the Sierra Nevada, but tritium and carbon-14 data suggest recharge rates in this area are low and have decreased over recent geologic time. Managed aquifer recharge at two sites show differences in water-level responses to recharge and in the physical movement of recharged water with depth related to subsurface geology. Well-bore flow logs also show rapid movement of water from recharge sites through permeable interconnected aquifer layers to pumping wells.

California↗

Groundwater quality in the Monterey Bay and Salinas Valley groundwater basins, California

The Monterey-Salinas study unit is nearly 1,000 square miles and consists of the Santa Cruz Purisima Formation Highlands, Felton Area, Scotts Valley, Soquel Valley, West Santa Cruz Terrace, Salinas Valley, Pajaro Valley, and Carmel Valley groundwater basins (California Department of Water Resources, 2003; Kulongski and Belitz, 2011). These basins were grouped into four study areas based primarily on geography. Groundwater basins in the north were grouped into the Santa Cruz study area, and those to the south were grouped into the Monterey Bay, the Salinas Valley, and the Paso Robles study areas (Kulongoski and others, 2007). The study unit has warm, dry summers and cool, moist winters. Average annual rainfall ranges from 31 inches in Santa Cruz in the north to 13 inches in Paso Robles in the south. The study areas are drained by several rivers and their principal tributaries: the Salinas, Pajaro, and Carmel Rivers, and San Lorenzo Creek. The Salinas Valley is a large intermontane valley that extends southeastward from Monterey Bay to Paso Robles. It has been filled, up to a thickness of 2,000 feet, with Tertiary and Quaternary marine and terrestrial sediments that overlie granitic basement. The Miocene-age Monterey Formation and Pliocene- to Pleistocene-age Paso Robles Formation, and Pleistocene to Holocene-age alluvium contain freshwater used for supply. The primary aquifers in the study unit are defined as those parts of the aquifers corresponding to the perforated intervals of wells listed in the California Department of Public Health database. Public-supply wells are typically drilled to depths of 200 to 650 feet, consist of solid casing from the land surface to depths of about 175 to 500 feet, and are perforated below the solid casing. Water quality in the primary aquifers may differ from that in the shallower and deeper parts of the aquifer system. Groundwater movement is generally from the southern part of the Salinas Valley north towards the Monterey Bay. Land use in the study unit is about 44 percent (%) natural (mostly grassland and forests), 43% agricultural, and 13% urban. The primary agricultural uses are row crops, pasture, hay, and vineyards. The largest urban areas are the cities of Santa Cruz, Watsonville, Monterey, Salinas, King City, and Paso Robles. Recharge to the groundwater system is primarily from stream-channel infiltration from the major rivers and their tributaries, and from infiltration of water from precipitation and irrigation. The primary sources of discharge are water pumped for irrigation and municipal supply, evaporation, and discharge to streams.

Fact Sheet↗

Simulated effects of Lower Floridan aquifer pumping on the Upper Floridan aquifer at Barbour Pointe, Chatham County, Georgia

Steady-state simulations using a revised regional groundwater-flow model based on MODFLOW were run to assess the potential long-term effects on the Upper Floridan aquifer (UFA) of pumping the Lower Floridan aquifer (LFA) at well 36Q398, located at Barbour Pointe in coastal Georgia near Savannah. Simulated pumping of well 36Q398 at a rate of 750 gallons per minute (gal/min; or 1.08 million gallons per day [Mgal/d]) indicated a maximum drawdown of about 2.19 feet (ft) in the UFA directly above the pumped well and at least 1 ft of drawdown within a nearly 190-square-mile area (scenario A). Induced vertical leakage from the UFA provided about 98 percent of the water to the pumped well. Simulated pumping of well 36Q398 caused increased downward leakage in all layers above the LFA, decreased upward leakage in all layers above the LFA, increased inflow to and decreased outflow from lateral specified-head boundaries in the UFA and LFA, and an increase in the volume of induced inflow from the general-head boundary representing outcrop units. Water budgets for scenario A indicated that changes in inflows and outflows through general-head boundaries would compose about 45 percent of the simulated pumpage from well 36Q398, with the remaining 55 percent of the pumped water derived from flow across lateral specified-head boundaries. Additional steady-state simulations were run to evaluate a pumping rate in the UFA of 240 gal/min (0.346 Mgal/d), which would produce an equivalent maximum drawdown in the UFA as pumping from well 36Q398 in the LFA at a rate of 750 gal/min (called the “drawdown offset”; scenario B). Simulated pumping in the UFA for the drawdown offset produced about 2.18 ft of drawdown, comparable to 2.19 ft of drawdown in the UFA simulated in scenario A. Water budgets for scenario B also provided favorable comparisons with scenario A, indicating that 42 percent of the drawdown-offset pumpage (0.346 Mgal/d) in the UFA originates as increased inflow and decreased outflow across general-head boundaries from overlying units in the surficial and Brunswick aquifer systems and that the remaining simulated pumpage originates as flow across general- and specified-head boundaries within the UFA and LFA. The revised model was evaluated for sensitivity by first altering horizontal and vertical hydraulic conductivity in the Lower Floridan semiconfining unit and then adjusting horizontal and vertical hydraulic conductivity in the LFA to match the 35.6 ft of drawdown at pumping well 36Q398. These adjustments also affected the maximum simulated drawdown in the UFA and the equivalent offset pumping in the UFA that would produce the same amount of drawdown. The maximum drawdown in the UFA ranged from 1.82 to 2.57 ft and the equivalent offset pumping in the UFA ranged from 199 to 278 gal/min. The revised model reasonably depicts changes in groundwater levels resulting from pumping the LFA at Barbour Pointe at a rate of 750 gal/min. Results are limited, however, by the same model assumptions and design as the original model, and placement of boundaries and type of boundary used exert the greatest control on overall groundwater flow and interaquifer leakage in the system. Simulation results have improved regional characterization of the Floridan aquifer system, which could be used by State officials in evaluating requests for groundwater withdrawal from the LFA.

Georgia↗

Preliminary delineation and description of the regional aquifers of Tennessee: The east Tennessee aquifer system

The east Tennessee aquifer system occurs in the Valley and Ridge and the Blue Ridge provinces of Tennessee. These areas are underlain by rocks of Precambrian to Pennsylvanian age which have been structurally deformed and faulted during the Appalachian orogeny. Groundwater in the Valley and Ridge occurs primarily in solution openings in carbonate rocks and in fractures in sandstone and shale. Fractures in the crystalline rocks store and transmit most of the groundwater in the Blue Ridge province. The east Tennessee aquifer system is important as a source of rural and municipal drinking water. Within 300 ft of land surface, groundwater generally contains less than 500 milligrams/L dissolved solids. At greater depths, fractures and solution openings are smaller and fewer in number. There are very little data to define groundwater occurrence at depths greater than about 300 ft. Groundwater flow may be restricted and the dissolved solids concentrations in the groundwater may reach thousands or even ten thousands of milligrams/L.

Tennessee↗