SEARCH · Geology Reports

Results for “Cells”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 829 records · Page 46Linked to original sources

USGS2013

Investigations of groundwater system and simulation of regional groundwater flow for North Penn Area 7 Superfund site, Montgomery County, Pennsylvania

Groundwater in the vicinity of several industrial facilities in Upper Gwynedd Township and vicinity, Montgomery County, in southeast Pennsylvania has been shown to be contaminated with volatile organic compounds (VOCs), the most common of which is the solvent trichloroethylene (TCE). The 2-square-mile area was placed on the National Priorities List as the North Penn Area 7 Superfund site by the U.S. Environmental Protection Agency (USEPA) in 1989. The U.S. Geological Survey (USGS) conducted geophysical logging, aquifer testing, and water-level monitoring, and measured streamflows in and near North Penn Area 7 from fall 2000 through fall 2006 in a technical assistance study for the USEPA to develop an understanding of the hydrogeologic framework in the area as part of the USEPA Remedial Investigation. In addition, the USGS developed a groundwater-flow computer model based on the hydrogeologic framework to simulate regional groundwater flow and to estimate directions of groundwater flow and pathways of groundwater contaminants. The study area is underlain by Triassic- and Jurassic-age sandstones and shales of the Lockatong Formation and Brunswick Group in the Mesozoic Newark Basin. Regionally, these rocks strike northeast and dip to the northwest. The sequence of rocks form a fractured-sedimentary-rock aquifer that acts as a set of confined to partially confined layers of differing permeabilities. Depth to competent bedrock typically is less than 20 ft below land surface. The aquifer layers are recharged locally by precipitation and discharge locally to streams. The general configuration of the potentiometric surface in the aquifer is similar to topography, except in areas affected by pumping. The headwaters of Wissahickon Creek are nearby, and the stream flows southwest, parallel to strike, to bisect North Penn Area 7. Groundwater is pumped in the vicinity of North Penn Area 7 for industrial use, public supply, and residential supply. Results of field investigations by USGS at the site and results from other studies support, and are consistent with, a conceptual model of a layered leaky aquifer where the dip of the beds has a strong control on hydraulic connections in the groundwater system. Connections within and (or) parallel to bedding tend to be greater than across bedding. Transmissivities of aquifer intervals isolated by packers ranged over three orders of magnitude [from about 2.8 to 2,290 square feet per day (ft 2 /d) or 0.26 to 213 square meters per day (m 2 /d)], did not appear to differ much by mapped geologic unit, but showed some relation to depth being relatively smaller in the shallowest and deepest intervals (0 to 50 ft and more than 250 ft below land surface, respectively) compared to the intermediate depth intervals (50 to 250 ft below land surface) tested. Transmissivities estimated from multiple-observation well aquifer tests ranged from about 700 to 2,300 ft 2 /d (65 to 214 m 2 /d). Results of chemical analyses of water from isolated intervals or monitoring wells open to short sections of the aquifer show vertical differences in concentrations; chloride and silica concentrations generally were greater in shallow intervals than in deeper intervals. Chloride concentrations greater than 100 milligrams per liter (mg/L), combined with distinctive chloride/bromide ratios, indicate a different source of chloride in the western part of North Penn Area 7 than elsewhere in the site. Groundwater flow at a regional scale under steady-state conditions was simulated by use of a numerical model (MODFLOW-2000) for North Penn Area 7 with different layers representing saprolite/highly weathered rock near the surface and unweathered competent bedrock. The sedimentary formations that underlie the study area were modeled using dipping model layers for intermediate and deep zones of unweathered, fractured rock. Horizontal cell model size was 100 meters (m) by 100 meters (328 ft by 328 ft), and model layer thickness ranged from 6 m (19.7 ft) representing shallow weathered rock and saprolite up to 200 m (656 ft) representing deeper dipping bedrock. The model did not include detailed structure to account for local-scale differences in hydraulic properties, with the result that local-scale groundwater flow may not be well simulated. Additional detailed multi-well aquifer tests would be needed to establish the extent of interconnection between intervals at the local scale to address remediation of contamination at each source area. This regional groundwater-flow model was calibrated against measured groundwater levels (1996, 2000, and 2005) and base flow estimated from selected streamflow measurements by use of nonlinear-regression parameter-estimation algorithms to determine hydraulic conductivity and anisotropy of hydraulic conductivity, streambed hydraulic conductivity, and recharge during calibration periods. Results of the simulation using the calibrated regional model indicate that the aquifer appears to be anisotropic where hydraulic conductivity is greatest parallel to the orientation of bedding of the formations underlying the area and least in the cross-bed direction. The maximum hydraulic conductivity is aligned with the average regional strike of the formations, which is “subhorizontal” in the model because the altitudes of the beds and model cells vary in the strike, as well as dip, direction. Estimated subhorizontal hydraulic conductivities (in strike direction parallel to dipping beds) range from 0.001 to 1.67 meters per day (0.0032 to 5.5 feet per day). The ratio of minimum (dip direction) to maximum (strike direction) subhorizontal hydraulic conductivity ranges from 1/3.1 to 1/8.6, and the ratio of vertical to horizontal hydraulic conductivity ranges from 1/1 to 1/478. However, limited available field data precluded rigorous calibration of vertical anisotropy in the model. Estimated recharge rates corresponding to calibration periods in 1996, 2000, and 2005 are 150, 109, and 124 millimeters per year (5.9, 4.3, and 4.9 inches per year), respectively. The calibrated groundwater-flow model was used to simulate groundwater flow under steady-state conditions during periods of relatively high withdrawals (pumpage) (1990) and relatively low withdrawals (2000 and 2005). Groundwater-flow paths originating from recharge areas near known areas of soil contamination (sources) were simulated. Pumped industrial and production wells captured more groundwater from several of these sources during 1990 than after 1990 when pumping declined or ceased and greater amounts of contaminated groundwater moved away from North Penn Area 7 Superfund site to surrounding areas. Uncertainty in simulated groundwater-flow paths from contaminant sources and contributing areas, resulting from uncertainty in estimated hydraulic properties of the model, was illustrated through Monte Carlo simulations. The effect of uncertainty in the vertical anisotropy was not included in the Monte Carlo simulations. Contributing areas indicating the general configuration of groundwater flow towards production well MG-202 (L-22) in the study area also were simulated for the different time periods; as simulated, the flow paths do not pass through any identified contaminant source in North Penn Area 7. However, contributing areas to wells, such as MG-202, located near many pumped wells are particularly complex and, in some cases, include areas that contribute flow to streams that subsequently recharge the aquifer through stream loss. In these cases, water-quality constituents, including contaminants that are present in surface water may be drawn into the aquifer to nearby pumped wells. Results of a simulated shutdown of well MG-202 under steady-state 2005 conditions showed that the area contributing recharge for nearby production well MG-76 (L-17), when MG-202 is not pumping, shifts downstream and is similar to the area contributing recharge for MG-202 when both wells are pumping. Concentrations of constituents in groundwater samples collected in fall 2005 or spring 2006 were compared to simulated groundwater-flow paths for the year 2005 to provide a qualitative assessment of model results. The observed spatial distribution of selected constituents, including TCE, CFC-11, and CFC-113 in groundwater in 2005 and the chloride/bromide mass ratios in 2006, generally were consistent with the model results of the simulated 2005 groundwater-flow paths at North Penn Area 7, indicating the presence of several separate sources of contaminants within North Penn Area 7.

Pennsylvania↗

USGS2013

Conceptual and numerical models of groundwater flow in the Ogallala aquifer in Gregory and Tripp Counties, South Dakota, water years 1985--2009

The Ogallala aquifer is an important water resource for the Rosebud Sioux Tribe in Gregory and Tripp Counties in south-central South Dakota and is used for irrigation, public supply, domestic, and stock water supplies. To better understand groundwater flow in the Ogallala aquifer, conceptual and numerical models of groundwater flow were developed for the aquifer. A conceptual model of the Ogallala aquifer was used to analyze groundwater flow and develop a numerical model to simulate groundwater flow in the aquifer. The MODFLOW–NWT model was used to simulate transient groundwater conditions for water years 1985–2009. The model was calibrated using statistical parameter estimation techniques. Potential future scenarios were simulated using the input parameters from the calibrated model for simulations of potential future drought and future increased pumping. Transient simulations were completed with the numerical model. A 200-year transient initialization period was used to establish starting conditions for the subsequent 25-year simulation of water years 1985–2009. The 25-year simulation was discretized into three seasonal stress periods per year and used to simulate transient conditions. A single-layer model was used to simulate flow and mass balance in the Ogallala aquifer with a grid of 133 rows and 282 columns and a uniform spacing of 500 meters (1,640 feet). Regional inflow and outflow were simulated along the western and southern boundaries using specified-head cells. All other boundaries were simulated using no-flow cells. Recharge to the aquifer occurs through precipitation on the outcrop area. Model calibration was accomplished using the Parameter Estimation (PEST) program that adjusted individual model input parameters and assessed the difference between estimated and model-simulated values of hydraulic head and base flow. This program was designed to estimate parameter values that are statistically the most likely set of values to result in the smallest differences between simulated and observed values, within a given set of constraints. The potentiometric surface of the aquifer calculated during the 200-year initialization period established initial conditions for the transient simulation. Water levels for 38 observation wells were used to calibrate the 25-year simulation. Simulated hydraulic heads for the transient simulation were within plus or minus 20 feet of observed values for 95 percent of observation wells, and the mean absolute difference was 5.1 feet. Calibrated hydraulic conductivity ranged from 0.9 to 227 feet per day (ft/d). The annual recharge rates for the transient simulation (water years 1985–2009) ranged from 0.60 to 6.96 inches, with a mean of 3.68 inches for the Ogallala aquifer. This represents a mean recharge rate of 280.5 ft 3 /s for the model area. Discharge from the aquifer occurs through evapotranspiration, discharge to streams through river leakage and flow from springs and seeps, and well withdrawals. Water is withdrawn from wells for irrigation, public supply, domestic, and stock uses. Simulated mean discharge rates for water years 1985–2009 were about 185 cubic feet per second (ft 3 /s) for evapotranspiration, 66.7 ft 3 /s for discharge to streams, and 5.48 ft 3 /s for well withdrawals. Simulated annual evapotranspiration rates ranged from about 128 to 254 ft 3 /s, and outflow to streams ranged from 52.2 to 79.9 ft 3 /s. A sensitivity analysis was used to examine the response of the calibrated model to changes in model parameters for horizontal hydraulic conductivity, recharge, evapotranspiration, and spring and riverbed conductance. The model was most sensitive to recharge and maximum potential evapotranspiration and least sensitive to riverbed and spring conductances. Two potential future scenarios were simulated: a potential drought scenario and a potential increased pumping scenario. To simulate a potential drought scenario, a synthetic drought record was created, the mean of which was equal to 60 percent of the mean estimated recharge rate for the 25-year simulation period. Compared with the results of the calibrated model (non-drought simulation), the simulation representing a potential drought scenario resulted in water-level decreases of as much as 30 feet for the Ogallala aquifer. To simulate the effects of potential future increases in pumping, well withdrawal rates were increased by 50 percent from those estimated for the 25-year simulation period. Compared with the results of the calibrated model, the simulation representing an increased pumping scenario resulted in water-level decreases of as much as 26 feet for the Ogallala aquifer. Groundwater budgets for the potential future scenario simulations were compared with the transient simulation representing water years 1985–2009. The simulation representing a potential drought scenario resulted in lower aquifer recharge from precipitation and decreased discharge from streams, springs, seeps, and evapotranspiration. The simulation representing a potential increased pumping scenario was similar to results from the transient simulation, with a slight increase in well withdrawals and a slight decrease in discharge from river leakage and evapotranspiration. This numerical model is suitable as a tool that could be used to better understand the flow system of the Ogallala aquifer, to approximate hydraulic heads in the aquifer, and to estimate discharge to rivers, springs, and seeps in the study area. The model also is useful to help assess the response of the aquifer to additional stresses, including potential drought conditions and increased well withdrawals.

South Dakota↗

USGS2013

Changing patterns in the use, recycling, and material substitution of mercury in the United States

Environmental concerns have led to numerous regulations that have dramatically decreased the reported production and use of mercury in the United States since the 1980s. Government legislation and subsequent industry actions have led to increased collection of mercury-containing materials and the recovery of mercury through recycling. Mercury emissions have been reduced and effective alternatives to mercury products have been developed for many applications. This study updates and quantifies the changes in demand, supply, use, and material flow for mercury in various sectors in the United States that have taken place since 1996. Nearly all primary mercury produced in the United States is derived as a byproduct of processing of gold and silver ore in Nevada. Since 2001, annual production of mercury from gold and silver mining in Nevada has decreased by 22 percent overall because ore from greater depths containing low grade mercury is recovered, and mercury emissions from this source have decreased by 95 percent as a result of increased regulation and improved collection and suppression technology. The distribution of consumption of mercury in the United States has changed as a result of regulation (elimination of large-scale mercury use in the paint and battery sectors), reduction by consumers (decommissioning of mercury-cell chloralkali manufacturing capacity), and technological advances (improvements in dental, lighting, and wiring sectors). Mercury use in the chloralkali sector, the leading end-use sector in the United States in 1996, has declined by 98 percent from 136 metric tons (t) in 1996 to about 0.3 t in 2010 because of increased processing and recycling efficiencies and plant closures or conversion to other technologies. As plants were closed, mercury recovered from the infrastructure of decommissioned plants has been exported, making the United States a net exporter of mercury, even though no mercury has been produced as the primary product from mines in the United States since 1992. In 1996, the three leading end-use sectors for mercury in the United States were chloralkali manufacturing (accounting for 38 percent of consumption), electrical and electronic instrumentation (13 percent of consumption), and instruments and measuring devices (11 percent of consumption). In 2010, the three leading end-use sectors were dental amalgam (accounting for between 35 and 57 percent of consumption), electrical and electronic instrumentation (29 percent of consumption), and batteries (8 percent of consumption). Mercury use in lighting is increasing because incandescent lights are being phased out in favor of mercury-containing compact fluorescent bulbs, but the demand for mercury per unit produced is small. Dental amalgam constituted the largest amount of mercury in use in the United States. One study reported about 290 t of mercury in dental amalgam was estimated to be contained in human mouths, an estimated 30 t of mercury amalgam was treated as waste, 28.5 t of mercury amalgam was released to the environment, 6 t of amalgam was recycled, and 3.5 t was treated and stored in landfills in 2009. Mercury contained in products recovered by State, municipal, or industry collection activities is recycled, but the estimated overall recycling rate is less than 10 percent. Increasingly, the U.S. mercury recycling industry has been processing a significant amount of mercury-containing material derived from foreign gold mining operations or decommissioned mercury-cell chloralkali plants. Regulation of mercury export and storage is expected to result in surplus mercury inventories in the United States. The Mercury Export Ban Act of 2008 limits elemental mercury exports for unregulated uses such as artisanal gold mining after January 1, 2013, and requires development of adequate long-term storage facilities in the United States for elemental mercury. During the past 4 years, producers and recyclers of elemental mercury have been exporting large quantities of mercury in anticipation of this regulation, but the U.S. inventory of mercury in 2010 was estimated to have exceeded 7,000 t from Government stockpiles and industry stocks. Costs attributed to long-term storage may affect the competitiveness of mercury recycling.

Scientific Investigations Report↗

USGS2013

Revised shallow and deep water-level and storage-volume changes in the Equus Beds Aquifer near Wichita, Kansas, predevelopment to 1993

Beginning in the 1940s, the Wichita well field was developed in the Equus Beds aquifer in southwestern Harvey County and northwestern Sedgwick County to supply water to the city of Wichita. The decline of water levels in the aquifer was noted soon after the development of the Wichita well field began. Development of irrigation wells began in the 1960s. City and agricultural withdrawals led to substantial water-level declines. Water-level declines enhanced movement of brines from past oil and gas activities near Burrton, Kansas and enhanced movement of natural saline water from the Arkansas River into the well field area. Large chloride concentrations may limit use or require the treatment of water from the well field for irrigation or public supply. In 1993, the city of Wichita adopted the Integrated Local Water Supply Program (ILWSP) to ensure an adequate water supply for the city through 2050 and as part of its effort to effectively manage the part of the Equus Beds aquifer it uses. ILWSP uses several strategies to do this including the Equus Beds Aquifer Storage and Recovery (ASR) project. The purpose of the ASR project is to store water in the aquifer for later recovery and to help protect the aquifer from encroachment of a known oilfield brine plume near Burrton and saline water from the Arkansas River. As part of Wichita’s ASR permits, Wichita is prohibited from artificially recharging water into the aquifer in a Basin Storage area (BSA) grid cell if water levels in that cell are above the January 1940 water levels or are less than 10 feet below land surface. The map previously used for this purpose did not provide an accurate representation of the shallow water table. The revised predevelopment water-level altitude map of the shallow part of the aquifer is presented in this report. The city of Wichita’s ASR permits specify that the January 1993 water-level altitudes will be used as a lower baseline for regulating the withdrawal of artificial rechage credits from the Equus Beds aquifer by the city of Wichita. The 1993 water levels correspond to the lowest recorded levels and largest storage declines since 1940. Revised and new water-level maps of shallow and deep layers were developed to better represent the general condition of the aquifer. Only static water levels were used to better represent the general condition of the aquifer and comply with Wichita’s ASR permits. To ensure adequate data density, the January 1993 period was expanded to October 1992 through February 1993. Static 1993 water levels from the deep aquifer layer of the Equus Beds aquifer possibly could be used as the lower baseline for regulatory purposes. Previously, maps of water-level changes used to estimate the storage-volume changes included a combination of static (unaffected by pumping or nearby pumping) and stressed (affected by pumping or nearby pumping) water levels from wells. Some of these wells were open to the shallow aquifer layer and some were open to the deep aquifer layer of the Equus Beds aquifer. In this report, only static water levels in the shallow aquifer layer were used to determine storage-volume changes. The effects on average water-level and storage-volume change from the use of mixed, stressed water levels and a specific yield of 0.20 were compared to the use of static water levels in the shallow aquifer and a specific yield of 0.15. This comparison indicates that the change in specific yield causes storage-volume changes to decrease about 25 percent, whereas the use of static water levels in the shallow aquifer layer causes an increase of less than 4 percent. Use of a specific yield of 0.15 will result in substantial decreases in the amount of storage-volume change compared to those reported previously that were calculated using a specific yield of 0.20. Based on these revised water-level maps and computations, the overall decline and change in storage from predevelopment to 1993 represented a loss in storage of about 6 percent (-202,000 acre-feet) of the overall storage volume within the newly defined study area.

Kansas↗

USGS2014

Documentation of a groundwater flow model (SJRRPGW) for the San Joaquin River Restoration Program study area, California

To better understand the potential effects of restoration flows on existing drainage problems, anticipated as a result of the San Joaquin River Restoration Program (SJRRP), the U.S. Geological Survey (USGS), in cooperation with the U.S. Bureau of Reclamation (Reclamation), developed a groundwater flow model (SJRRPGW) of the SJRRP study area that is within 5 miles of the San Joaquin River and adjacent bypass system from Friant Dam to the Merced River. The primary goal of the SJRRP is to reestablish the natural ecology of the river to a degree that restores salmon and other fish populations. Increased flows in the river, particularly during the spring salmon run, are a key component of the restoration effort. A potential consequence of these increased river flows is the exacerbation of existing irrigation drainage problems along a section of the river between Mendota and the confluence with the Merced River. Historically, this reach typically was underlain by a water table within 10 feet of the land surface, thus requiring careful irrigation management and (or) artificial drainage to maintain crop health. The SJRRPGW is designed to meet the short-term needs of the SJRRP; future versions of the model may incorporate potential enhancements, several of which are identified in this report. The SJRRPGW was constructed using the USGS groundwater flow model MODFLOW and was built on the framework of the USGS Central Valley Hydrologic Model (CVHM) within which the SJRRPGW model domain is embedded. The Farm Process (FMP2) was used to simulate the supply and demand components of irrigated agriculture. The Streamflow-Routing Package (SFR2) was used to simulate the streams and bypasses and their interaction with the aquifer system. The 1,300-square mile study area was subdivided into 0.25-mile by 0.25-mile cells. The sediment texture of the aquifer system, which was used to distribute hydraulic properties by model cell, was refined from that used in the CVHM to better represent the natural heterogeneity of aquifer-system materials within the model domain. In addition, the stream properties were updated from the CVHM to better simulate stream-aquifer interactions, and water-budget subregions were refined to better simulate agricultural water supply and demand. External boundary conditions were derived from the CVHM. The SJRRPGW was calibrated for April 1961 to September 2003 by using groundwater-level observations from 133 wells and streamflow observations from 19 streamgages. The model was calibrated using public-domain parameter estimation software (PEST) in a semi-automated manner. The simulated groundwater-level elevations and trends (including seasonal fluctuations) and surface-water flow magnitudes and trends reasonably matched observed data. The calibrated model is planned to be used to assess the potential effects of restoration flows on agricultural lands and the relative capabilities of proposed SJRRP actions to reduce these effects.

California↗

USGS2014

Modeling uncertainty in coal resource assessments, with an application to a central area of the Gillette coal field, Wyoming

Standards for the public disclosure of mineral resources and reserves do not require the use of any specific methodology when it comes to estimating the reliability of the resources. Unbeknownst to most intended recipients of resource appraisals, such freedom commonly results in subjective opinions or estimations based on suboptimal approaches, such as use of distance methods. This report presents the results of a study of the third of three coal deposits in which drilling density has been increased one order of magnitude in three stages. Applying geostatistical simulation, the densest dataset was used to check the results obtained by modeling the sparser drillings. We have come up with two summary displays of results based on the same simulations, which individually and combined provide a better assessment of uncertainty than traditional qualitative resource classifications: (a) a display of cell 90 percent confidence interval versus cumulative cell tonnage, and (b) a histogram of total resources. The first graph allows classification of data into any number of bins with dividers to be decided by the assessor on the basis of a discriminating variable that is statistically accepted as a measure of uncertainty, thereby improving the quality and flexibility of the modeling. The second display expands the scope of the modeling by providing a quantitative measure of uncertainty for total tonnage, which is a fundamental concern for stockholders, geologists, and decision makers. Our approach allows us to correctly model uncertainty issues not possible to predict with distance methods, such as (a) different levels of uncertainty for individual beds with the same pattern and density of drill holes, (b) different local degrees of reduction of uncertainty with drilling densification reflecting fluctuation in the complexity of the geology, (c) average reduction in uncertainty at a disproportionately lesser rate than the reduction in area per drill hole, (d) the proportional effect of higher uncertainty in areas of higher tonnages, despite a regular drilling pattern, (e) the possibility of a local increase in uncertainty despite drilling densification to reflect a more complex geology as the deposit is known in more detail, and (f) for exactly the same drilling pattern, tonnage per individual beds with different uncertainty than the aggregated tonnage. These results should be considered realistic improvements over distance methods used for quantitative classification of uncertainty in coal resource, such as U.S. Geological Survey Circular 891. 1 The approach should be a welcome addition to the toolkit of Competent Persons preparing public disclosures according to international mineral codes such as those promoted by the Combined Reserves International Reporting Standards Committee 2 and the Joint Ore Reserve Committee. 3 1 Wood, G.H., Jr., Kehn, T.M., Carter, M.D., and Culbertson, W.C., 1983, Coal resources classification system of the U.S. Geological Survey: U.S. Geological Survey Circular 891, 65 p. 2 CRIRSCO (Combined Reserves International Reporting Standards Committee), 2013, International reporting template for the public reporting of exploration results, mineral resources and mineral reserves: Accessed February 2014 at http://www.crirsco.com/crirsco_template_may2013.pdf. 3 JORC (Joint Ore Reserves Committee), 2012, Australasian code for reporting of exploration results, mineral resources and ore reserves: Accessed September 2014 at http://www.jorc.org/docs/jorc_code2012.pdf.

Wyoming↗

USGS2015

Potential groundwater recharge for the State of Minnesota using the Soil-Water-Balance model, 1996-2010

Groundwater recharge is one of the most difficult components of a water budget to ascertain, yet is an important boundary condition necessary for the quantification of water resources. In Minnesota, improved estimates of recharge are necessary because approximately 75 percent of drinking water and 90 percent of agricultural irrigation water in Minnesota are supplied from groundwater. The water that is withdrawn must be supplied by some combination of (1) increased recharge, (2) decreased discharge to streams, lakes, and other surface-water bodies, and (3) removal of water that was stored in the system. Recent pressure on groundwater resources has highlighted the need to provide more accurate recharge estimates for various tools that can assess the sustainability of long-term water use. As part of this effort, the U.S. Geological Survey, in cooperation with the Minnesota Pollution Control Agency, used the Soil-Water-Balance model to calculate gridded estimates of potential groundwater recharge across Minnesota for 1996‒2010 at a 1-kilometer (0.621-mile) resolution. The potential groundwater recharge estimates calculated for Minnesota from the Soil-Water Balance model included gridded values (1-kilometer resolution) of annual mean estimates (that is, the means for individual years from 1996 through 2010) and mean annual estimates (that is, the mean for the 15-year period 1996−2010). The Soil-Water-Balance model uses a modified Thornthwaite-Mather soil-water-balance approach, with components of the soil-water balance calculated on a daily basis. A key advantage of this approach includes the use of commonly available geographic information system data layers that incorporate land cover, soil properties, and daily meteorological data to produce temporally and spatially variable gridded estimates of potential recharge. The Soil-Water-Balance model was calibrated by using a combination of parameter estimation techniques, making manual adjustments of model parameters, and using parameter values from previously published Soil-Water-Balance models. Each calibration simulation compared the potential recharge estimate from the model against base-flow estimates derived from three separate hydrograph separation techniques. A total of 35 Minnesota watersheds were selected for the model calibration. Meteorological data necessary for the model included daily precipitation, minimum daily temperature, and maximum daily temperature. All of the meteorological data were provided by the Daymet dataset, which included daily continuous surfaces of key climatological data. Land-cover data were provided by the 2001 and 2006 National Land Cover Database: the 2001 classification was used from 1994 through 2003, and the 2006 classification was used from 2004 through 2010. Soil data used in the model included hydrologic soils group and the available soil-water capacity. These soil data were obtained from the Natural Resources Conservation Service Soil Survey Geographic (SSURGO) database and the State Soil Geographic (STATSGO) database. The statewide mean annual potential recharge rate from 1996–2010 was 4.9 inches per year. Potential recharge estimates increased from west to east across Minnesota. The mean annual potential recharge estimates across Minnesota at a 1-km resolution for the overall simulation period (1996–2010) ranged from less than 0.1 to 17.8 inches per year. Some of the lowest potential recharge rates for the simulation period were in the Red River of the North Basin of northwestern Minnesota, and generally were between 1.0 and 1.5 inches per year. The highest potential recharge rates were in northeastern Minnesota and the Anoka Sand Plain in central Minnesota. Eighty-eight percent of the potential recharge rates (by grid cell) were between 2 and 8 inches per year from 1996–2010. Only about 3 percent of all the potential recharge estimates (by grid cell) were less than 2 inches per year, and 9 percent of estimates were greater than 8 inches per year. On an annual basis, however, potential recharge rates were as high as 27.2 inches per year. The highest annual mean recharge estimate across the State was for 2010, and the lowest mean recharge estimate was for 2003. Although precipitation variability partially explained the annual differences in potential recharge estimates, precipitation alone did not account for these differences, and other factors such as antecedent moisture conditions likely were important. Also, because precipitation gradients across the State can vary from year to year, the dominant land-cover class and hydrologic soil group combinations for a particular region had a large effect on the resulting potential recharge value. During 1996–2010, April had the greatest monthly mean potential recharge compared to all other months, accounting for a mean of 30 percent of annual potential recharge in this single month.

Minnesota↗

USGS2015

Median nitrate concentrations in groundwater in the New Jersey Highlands Region estimated using regression models and land-surface characteristics

Nitrate-concentration data are used in conjunction with land-use and land-cover data to estimate median nitrate concentrations in groundwater underlying the New Jersey (NJ) Highlands Region. Sources of data on nitrate in 19,670 groundwater samples are from the U.S. Geological Survey (USGS) National Water Information System (NWIS) and the NJ Private Well Testing Act (PWTA). In a study conducted by the USGS, in cooperation with the New Jersey Department of Environmental Protection, logistic regression was used to relate measured nitrate concentrations to five explanatory variables (percent urban and agricultural land use, septic-system density, total length of streams, and number of known contaminated sites) quantified in 610-meter-square grid cells. A method for calculating the median concentrations of nitrate from a series of logistic regression models was developed. Two calibration and two validation procedures showed that the logistic-regression-based method can estimate groundwater-nitrate concentrations in the Highlands Region accurately to within 0.1 milligram per liter as nitrogen (mg/L as N). Limitations of the logistic-regression-based method include the inability to select a logistic model with exactly 0.5 probability of exceeding the threshold value and lack of an algorithm to directly calculate the median value. Quantile regression was evaluated as a suitable alternative and was slightly less accurate than the logistic-regression method in estimating median groundwater nitrate concentrations in the Highlands Region. Multiple-linear regression with log-transformed nitrate-concentration data and the same five explanatory values was less accurate than either logistic or quantile regression in estimating median nitrate concentrations. On the basis of 4,516 2000 x 2000 foot grid cells that contain wells with data stored in NWIS and the PWTA database, the estimated median nitrate concentration for the entire Highlands Region is about 1.25 mg/L as N, and estimated median concentrations range from about 1.05 to 1.78 mg/L as N among 11 smaller administratively defined areas within the Highlands Region that vary in percentages of urban land use, agricultural land use, and septic-system density. The Kaplan-Meier method of estimating summary statistics from left-censored data was applied in order to include nondetects (left-censored data) in median nitrate-concentration calculations. Median concentrations also were determined using three alternative methods of handling nondetects. Treatment of the 23 percent of samples that were nondetects had little effect on estimated median nitrate concentrations because method detection limits were mostly less than median values.

New Jersey↗

USGS2016

The source of groundwater and solutes to Many Devils Wash at a former uranium mill site in Shiprock, New Mexico

The Shiprock Disposal Site is the location of the former Navajo Mill (Mill), a uranium ore-processing facility, located on a terrace overlooking the San Juan River in the town of Shiprock, New Mexico. Following the closure of the Mill, all tailings and associated materials were encapsulated in a disposal cell built on top of the former Mill and tailings piles. The milling operations, conducted at the site from 1954 to 1968, created radioactive tailings and process-related wastes that are now found in the groundwater. Elevated concentrations of constituents of concern—ammonium, manganese, nitrate, selenium, strontium, sulfate, and uranium—have also been measured in groundwater seeps in the nearby Many Devils Wash arroyo, leading to the inference that these constituents originated from the Mill. These constituents have also been reported in groundwater that is associated with Mancos Shale, the bedrock that underlies the site. The objective of this report is to increase understanding of the source of water and solutes to the groundwater beneath Many Devils Wash and to establish the background concentrations for groundwater that is in contact with the Mancos Shale at the site. This report presents evidence on three working hypotheses: (1) the water and solutes in Many Devils Wash originated from the operations at the former Mill, (2) groundwater in deep aquifers is upwelling under artesian pressure to recharge the shallow groundwater beneath Many Devils Wash, and (3) the groundwater beneath Many Devils Wash originates as precipitation that infiltrates into the shallow aquifer system and discharges to Many Devils Wash in a series of springs on the east side of the wash. The solute concentrations in the shallow groundwater of Many Devils Wash would result from the interaction of the water and the Mancos Shale if the source of water was upwelling from deep aquifers or precipitation. In order to compare the groundwater from various wells to groundwater that has been affected by Mill activities, a classification system was developed to determine which wells were most likely to have been affected. Affects to groundwater by the Mill were determined by using the reported uranium alpha activity ratios measured in groundwater samples, along with the concentration of the uranium and the location of the wells relative to the Mill. Activity ratios of 1.2 or less were determined to be the most reliable indicator of Mill-affected groundwater. Wells with samples that had a reported activity ratio of 1.2 or less were classified as Mill affected. To compare groundwater with background water-quality, data from groundwater seeps and springs in the Upper Eagle Nest Arroyo and Salt Creek Wash, located north of the San Juan River, are also presented and analyzed. Based on groundwater elevations and tritium concentrations measured in wells located between the disposal cell and Many Devils Wash, Mill water is not likely to reach Many Devils Wash. The tritium concentrations also indicate that groundwater from the Mill has not substantially affected Many Devils Wash in the past. Upwelling from deep aquifers was also determined to be an unlikely source, primarily by comparing the composition of the stable isotopes of water in the shallow groundwater with those reported in groundwater samples from the deeper aquifers. The stable-isotope compositions of the shallow groundwater around the site are enriched relative to the San Juan River and local meteoric lines, which suggests that most of the shallow groundwater has been influenced by evaporation and therefore was recharged at the surface. Several observations indicate that focused recharge is the likely source of groundwater in the area of Many Devils Wash. The visible erosional features in Many Devils Wash provide evidence of piping and groundwater sapping, and the distribution and type of vegetation in Many Devils Wash suggest that the focused recharge of precipitation is occurring. The estimated recharge from precipitation was calculated to be 0.0008 inches per year (in/yr) by using the mass-balance approach from reported seep discharge and 0.0011 in/yr using the chloride mass-balance approach. A conceptual model of groundwater quality beneath Many Devils Wash is presented to explain the source of solutes in the groundwater beneath Many Devils Wash. The major-ion concentrations and geochemical evolution in the groundwater beneath Many Devils Wash and across the study area support the conceptual model that the underlying Mancos Shale is the source of solutes. Differences in the major-ion composition between groundwater samples collected around the site, result from the degree of weathering to the Mancos Shale. The cation distribution appears to be an indicator of effects from the Mill, with samples from the Mill-affected wells largely having a calcium/magnesium-sulfate composition that resembles the reported compositions of more weathered shale; however, that composition could change if the Mill-processed water flowed into areas where the Mancos Shale was less weathered. On the basis of the widespread presence of uranium in the Mancos Shale and the distribution of aqueous uranium in the analog sites and other sites in the region, it appears likely that uranium in the groundwater of Many Devils Wash is naturally sourced from the Mancos Shale.

New Mexico↗

USGS2016

Simulated effects of groundwater withdrawals from aquifers in Ocean County and vicinity, New Jersey

Rapid population growth since the 1930s in Ocean County and vicinity, New Jersey, has placed increasing demands upon the area’s freshwater resources. To examine effects of groundwater withdrawals, a three-dimensional groundwater-flow model was developed to simulate the groundwater-flow systems of five area aquifers: the unconfined Kirkwood-Cohansey aquifer system and Vincentown aquifer, and three confined aquifers— the Rio Grande water-bearing zone, the Atlantic City 800-foot sand, and the Piney Point aquifer. The influence of withdrawals is evaluated by using transient groundwater-flow model simulations that incorporate three withdrawal schemes. These are (1) no-withdrawal conditions; (2) 2000–03 withdrawal conditions, using reported monthly withdrawals at all production wells from January 2000 through December 2003; and (3) maximum-allocation withdrawal conditions using the maximum withdrawal allowed by New Jersey Department of Environmental Protection permits at each well. Particle tracking analysis, using results from model simulations, delineated particle flow paths from production wells to the point of recharge, and estimated particle travel times. Compared with no-withdrawal conditions, 2000–03 withdrawal conditions reduced the amount of groundwater flow out of the Kirkwood-Cohansey aquifer system into streams, increased the net flow of water into other layers, reduced net flow into or out of storage, and reduced flow from the Kirkwood-Cohansey aquifer system to constant head cells. Freshwater discharging to the Barnegat Bay-Little Egg Harbor estuary from streams and groundwater is essential to maintaining the ecology of the bay. Examination of selected stress periods indicates that simulated base flow in streams flowing into the Barnegat Bay-Little Egg Harbor estuary is reduced by as much as 49 cubic feet per second for 2000 to 2003 withdrawal conditions when compared with no-withdrawal conditions. In the three confined aquifers, water levels during periods of low recharge and high withdrawals, and high recharge and low withdrawals, were examined to determine seasonal effects on the confined flow systems. The simulated potentiometric surface of the Rio Grande water-bearing zone and the Atlantic City 800-foot sand during selected stress periods indicates substantial declines from no-withdrawal conditions to 2000–03 conditions as a result of groundwater withdrawals. Cones of depression in Toms River Township, Seaside Heights and Seaside Park Boroughs, and Barnegat Light Borough developed in the potentiometric surface of the Piney Point aquifer in response to withdrawals. Maximum-allocation withdrawals decreased flow out of the Kirkwood-Cohansey aquifer system to constant head cells, increased flow out of the aquifer system to adjacent and lower layers, and reduced groundwater discharge to streams when compared with 2000–03 withdrawal conditions. Increases in withdrawals from the Rio Grande water-bearing zone, the Atlantic City 800-foot sand, and the Piney Point aquifer result in an increase in simulated net groundwater flow into these aquifers. Base-flow reduction from 2000–03 conditions to maximum-allocation conditions of 25 to 29 cubic feet per second in all streams draining to the Barnegat Bay-Little Egg Harbor also is indicated. Potentiometric surfaces of the Rio Grande water-bearing zone, Atlantic City 800-foot sand, and the Piney Point aquifer during two stress periods of simulated maximum-allocation withdrawal conditions indicated the expansion of several cones of depression developed during 2000–03 withdrawals. Simulation of average 2000–03 withdrawal conditions indicated the extent to which the groundwater-flow system is susceptible to potential saltwater intrusion into near-shore wells. Travel time from recharge to discharge location ranged from 11 to approximately 50,700 years in near-shore Kirkwood-Cohansey aquifer system wells. Those in Seaside Heights Borough, in Island Beach State Park (Berkeley Township), and in Ship Bottom Borough have particle travel times from 140 to 12,000 years and flow paths that originated under Barnegat Bay or the Atlantic Ocean from the simulation of average maximum-allocation withdrawal conditions. Travel time along flow paths to wells screened in the Rio Grande water-bearing zone and the Atlantic City 800-foot sand from recharge to discharge point ranged from nearly 530 years to greater than 3.73 million years from the simulation of average 2000–03 withdrawal conditions. Particle tracking indicated that most wells screened in these aquifers derived a large part of their recharge from the Oswego River Basin, with a small portion of flow originating either beneath Barnegat Bay or to the east beneath the Atlantic Ocean. Travel time along flow paths that start beneath either Barnegat Bay or the Atlantic Ocean ranged from 2,300 to approximately 134,000 years from the simulation of average maximum-allocation withdrawal conditions."

New Jersey↗

USGS2016

Numerical simulation of the groundwater-flow system of the Kitsap Peninsula, west-central Washington

A groundwater-flow model was developed to improve understanding of water resources on the Kitsap Peninsula. The Kitsap Peninsula is in the Puget Sound lowland of west-central Washington, is bounded by Puget Sound on the east and by Hood Canal on the west, and covers an area of about 575 square miles. The peninsula encompasses all of Kitsap County, Mason County north of Hood Canal, and part of Pierce County west of Puget Sound. The peninsula is surrounded by saltwater, and the hydrologic setting is similar to that of an island. The study area is underlain by a thick sequence of unconsolidated glacial and interglacial deposits that overlie sedimentary and volcanic bedrock units that crop out in the central part of the study area. Twelve hydrogeologic units consisting of aquifers, confining units, and an underlying bedrock unit form the basis of the groundwater-flow model. Groundwater flow on the Kitsap Peninsula was simulated using the groundwater-flow model, MODFLOW‑NWT. The finite difference model grid comprises 536 rows, 362 columns, and 14 layers. Each model cell has a horizontal dimension of 500 by 500 feet, and the model contains a total of 1,227,772 active cells. Groundwater flow was simulated for transient conditions. Transient conditions were simulated for January 1985–December 2012 using annual stress periods for 1985–2004 and monthly stress periods for 2005–2012. During model calibration, variables were adjusted within probable ranges to minimize differences between measured and simulated groundwater levels and stream baseflows. As calibrated to transient conditions, the model has a standard deviation for heads and flows of 47.04 feet and 2.46 cubic feet per second, respectively. Simulated inflow to the model area for the 2005–2012 period from precipitation and secondary recharge was 585,323 acre-feet per year (acre-ft/yr) (93 percent of total simulated inflow ignoring changes in storage), and simulated inflow from stream and lake leakage was 43,905 acre-ft/yr (7 percent of total simulated inflow). Simulated outflow from the model primarily was through discharge to streams, lakes, springs, seeps, and Puget Sound (594,595 acre-ft/yr; 95 percent of total simulated outflow excluding changes in storage) and through withdrawals from wells (30,761 acre-ft/yr; 5 percent of total simulated outflow excluding changes in storage). Six scenarios were formulated with input from project stakeholders and were simulated using the calibrated model to provide representative examples of how the model could be used to evaluate the effects on water levels and stream baseflows of potential changes in groundwater withdrawals, in consumptive use, and in recharge. These included simulations of a steady-state system, no-pumping and return flows, 15-percent increase in current withdrawals in all wells, 80-percent decrease in outdoor water to simulate effects of conservation efforts, 15-percent decrease in recharge from precipitation to simulate a drought, and particle tracking to determine flow paths. Changes in water-level altitudes and baseflow amounts vary depending on the stress applied to the system in these various scenarios. Reducing recharge by 15 percent between 2005 and 2012 had the largest effect, with water-level altitudes declining throughout the model domain and baseflow amounts decreasing by as much as 18 percent compared to baseline conditions. Changes in pumping volumes had a smaller effect on the model. Removing all pumping and resulting return flows caused increased water-level altitudes in many areas and increased baseflow amounts of between 1 and 3 percent.

Washington↗

USGS2018

Construction and calibration of a groundwater-flow model to assess groundwater availability in the uppermost principal aquifer systems of the Williston Basin, United States and Canada

The U.S. Geological Survey developed a groundwater-flow model for the uppermost principal aquifer systems in the Williston Basin in parts of Montana, North Dakota, and South Dakota in the United States and parts of Manitoba and Saskatchewan in Canada as part of a detailed assessment of the groundwater availability in the area. The assessment was done because of the potential for increased demands and stresses on groundwater associated with large-scale energy development in the area. As part of this assessment, a three-dimensional groundwater-flow model was developed as a tool that can be used to simulate how the groundwater-flow system responds to changes in hydrologic stresses at a regional scale. The three-dimensional groundwater-flow model was developed using the U.S. Geological Survey’s numerical finite-difference groundwater model with the Newton-Rhapson solver, MODFLOW–NWT, to represent the glacial, lower Tertiary, and Upper Cretaceous aquifer systems for steady-state (mean) hydrological conditions for 1981‒2005 and for transient (temporally varying) conditions using a combination of a steady-state period for pre-1960 and transient periods for 1961‒2005. The numerical model framework was constructed based on existing and interpreted hydrogeologic and geospatial data and consisted of eight layers. Two layers were used to represent the glacial aquifer system in the model; layer 1 represented the upper one-half and layer 2 represented the lower one-half of the glacial aquifer system. Three layers were used to represent the lower Tertiary aquifer system in the model; layer 3 represented the upper Fort Union aquifer, layer 4 represented the middle Fort Union hydrogeologic unit, and layer 5 represented the lower Fort Union aquifer. Three layers were used to represent the Upper Cretaceous aquifer system in the model; layer 6 represented the upper Hell Creek hydrogeologic unit, layer 7 represented the lower Hell Creek aquifer, and layer 8 represented the Fox Hills aquifer. The numerical model was constructed using a uniform grid with square cells that are about 1 mile (1,600 meters) on each side with a total of about 657,000 active cells. Model calibration was completed by linking Parameter ESTimation (PEST) software with MODFLOW–NWT. The PEST software uses statistical parameter estimation techniques to identify an optimum set of input parameters by adjusting individual model input parameters and assessing the differences, or residuals, between observed (measured or estimated) data and simulated values. Steady-state model calibration consisted of attempting to match mean simulated values to measured or estimated values of (1) hydraulic head, (2) hydraulic head differences between model layers, (3) stream infiltration, and (4) discharge to streams. Calibration of the transient model consisted of attempting to match simulated and measured temporally distributed values of hydraulic head changes, stream base flow, and groundwater discharge to artesian flowing wells. Hydraulic properties estimated through model calibration included hydraulic conductivity, vertical hydraulic conductivity, aquifer storage, and riverbed hydraulic conductivity in addition to groundwater recharge and well skin. The ability of the numerical model to accurately simulate groundwater flow in the Williston Basin was assessed primarily by its ability to match calibration targets for hydraulic head, stream base flow, and flowing well discharge. The steady-state model also was used to assess the simulated potentiometric surfaces in the upper Fort Union aquifer, the lower Fort Union aquifer, and the Fox Hills aquifer. Additionally, a previously estimated regional groundwater-flow budget was compared with the simulated steady-state groundwater-flow budget for the Williston Basin. The simulated potentiometric surfaces typically compared well with the estimated potentiometric surfaces based on measured hydraulic head data and indicated localized groundwater-flow gradients that were topographically controlled in outcrop areas and more generalized regional gradients where the aquifers were confined. The differences between the measured and simulated (residuals) hydraulic head values for 11,109 wells were assessed, which indicated that the steady-state model generally underestimated hydraulic head in the model area. This underestimation is indicated by a positive mean residual of 11.2 feet for all model layers. Layer 7, which represents the lower Hell Creek aquifer, is the only layer for which the steady-state model overestimated hydraulic head. Simulated groundwater-level changes for the transient model matched within plus or minus 2.5 feet of the measured values for more than 60 percent of all measurements and to within plus or minus 17.5 feet for 95 percent of all measurements; however, the transient model underestimated groundwater-level changes for all model layers. A comparison between simulated and estimated base flows for the steady-state and transient models indicated that both models overestimated base flow in streams and underestimated annual fluctuations in base flow. The estimated and simulated groundwater budgets indicate the model area received a substantial amount of recharge from precipitation and stream infiltration. The steady-state model indicated that reservoir seepage was a larger component of recharge in the Williston Basin than was previously estimated. Irrigation recharge and groundwater inflow from outside the Williston Basin accounted for a relatively small part of total groundwater recharge when compared with recharge from precipitation, stream infiltration, and reservoir seepage. Most of the estimated and simulated groundwater discharge in the Williston Basin was to streams and reservoirs. Simulated groundwater withdrawal, discharge to reservoirs, and groundwater outflow in the Williston Basin accounted for a smaller part of total groundwater discharge. The transient model was used to simulate discharge to 571 flowing artesian wells within the model area. Of the 571 established flowing artesian wells simulated by the model, 271 wells did not flow at any time during the simulation because hydraulic head was always below the land-surface altitude. As hydraulic head declined throughout the simulation, 68 of these wells responded by ceasing to flow by the end of 2005. Total mean simulated discharge for the 571 flowing artesian wells was 55.1 cubic feet per second (ft 3 /s), and the mean simulated flowing well discharge for individual wells was 0.118 ft 3 /s. Simulated discharge to individual flowing artesian wells increased from 0.039 to 0.177 ft 3 /s between 1961 and 1975 and decreased to 0.102 ft 3 /s by 2005. The mean residual for 34 flowing wells with measured discharge was 0.014 ft 3 /s, which indicates the transient model overestimated discharge to flowing artesian wells in the model area. Model limitations arise from aspects of the conceptual model and from simplifications inherent in the construction and calibration of a regional-scale numerical groundwater-flow model. Simplifying assumptions in defining hydraulic parameters in space and hydrologic stresses and time-varying observational data in time can limit the capabilities of this tool to simulate how the groundwater-flow system responds to changes in hydrologic stresses, particularly at the local scale; nevertheless, the steady-state model adequately simulated flow in the uppermost principal aquifer systems in the Williston Basin based on the comparison between the simulated and estimated groundwater-flow budget, the comparison between simulated and estimated potentiometric surfaces, and the results of the calibration process.

Montana, North Dakota, South Dakota, Wyoming↗

USGS2018

Understanding the effect of salinity tolerance on cyanobacteria associated with a harmful algal bloom in Lake Okeechobee, Florida

In an effort to simulate the survival of cyanobacteria as they are transported from Lake Okeechobee to the estuarine habitats that receive waters from the lake, a bioassay encompassing a range of salinities was performed. An overall decline in cyanobacteria health in salinity treatments greater than 18 practical salinity units (psu) was indicated by loss of cell membrane integrity based on SYTOX® Green staining, but this loss varied by the kind of cyanobacteria present. Microcystis aeruginosa was tolerant of salinities up to 18 psu; however, higher salinities caused leaking of microcystin from the cells. Dolichospermum circinale, another common bloomformer in this system, did not tolerate salinities greater than 7.5 psu. Stimulation of mucilage production was observed and is likely a mechanism used by both species to protect organism viability. At 7.5 psu, microcystin increased relative to chlorophyll-a, providing some evidence of biosynthesis when M. aeruginosa is exposed to this salinity. This study indicates that as freshwater cyanobacteria are transported to brackish and marine waters, there will be a loss of membrane integrity which will lead to the release of cellular microcystin into the surrounding waterbody. Additional research would be needed to determine the exact effect of salinity on this relationship.

Florida↗

USGS2019

Spatial distribution of nutrients, chloride, and suspended sediment concentrations and loads determined by using different sampling methods in a cross section of the Trenton Channel of the Detroit River, Michigan, November 2014–November 2015

The Detroit River separates the United States and Canada as it flows from Lake St. Clair to Lake Erie. The Trenton Channel is a 13-kilometer-long branch of the Detroit River that flows to the west of Grosse Ile before rejoining the Detroit River near its mouth, just before the Detroit River flows into Lake Erie. The U.S. Environmental Protection Agency has listed both the Trenton Channel and Detroit River as Areas of Concern because of a list of Beneficial Use Impairments such as interrupted drinking-water services, loss of aquatic life, and reduced recreational use. Phosphorus loading from tributaries such as the Trenton Channel is one of the primary drivers of eutrophication in Lake Erie. The complex flow patterns and variable distribution of chemical constituents in the Trenton Channel make it difficult to accurately characterize the concentrations and loads of nutrients and other constituents conveyed through the channel to Lake Erie. In order to better understand the Trenton Channel’s contributions of nutrients (total phosphorus, orthophosphate, total nitrogen, and ammonia), chloride, and suspended sediment to Lake Erie and evaluate differences in results obtained by using different sample methodologies, the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency and Environment Canada, completed 12 sampling campaigns on the Trenton Channel in Detroit, Michigan, from November 2014 through November 2015. Acoustic Doppler current profiler (ADCP) techniques were used to characterize the distribution of velocity components within a cross section corresponding to a transect of the Trenton Channel at U.S. Geological Survey station 041686401 Trenton Channel of Detroit River at Grosse Ile, Mich. Three methods of collecting water-quality data at the same transect of the Trenton Channel were used: multiple-vertical depth-integrated (MVDI), fixed-point, and discrete samples. Horizontal and vertical variations in concentrations of nutrients, chloride, and suspended sediment were analyzed from discrete samples to better understand distributions of these constituents throughout the channel. Constituent loads were calculated by using individual sample concentrations and ADCP measurements for discharge made on the same day that the water-quality samples were collected. Constituent loads calculated from MVDI and fixed-point sampling methods were compared. The relation between MVDI and fixed-point samples helped quantify the differences between the sampling methods. Linear regression equations depicting the relation between concentrations measured by using MVDI and fixed-point samples were prepared. ADCP data indicates that velocities throughout the sampled transect remain uniform except for one location around 200 meters from the west bank of the channel. Secondary flow vectors suggest the presence of counter-rotating helical flow cells, and these helical flow cells could affect the mixing of constituents in transport by preventing cross-channel mixing. Flow discharges throughout the sampling campaign showed small variations, although lower flow rates were observed in the early winter months than in the summer months. Discrete sampling methods results displayed both heterogeneity throughout the channel horizontally, representing limited horizontal mixing in the channel, and displayed homogeneity throughout vertical transects, indicating mixing vertically. Comparisons between MVDI and fixed-point methods found consistently higher concentrations were measured in MVDI samples compared to concentrations measured in fixed-point samples. To correct for this bias between MVDI and fixed-point sample results, simple linear-regression equations were developed for all major constituents to help estimate constituent concentrations from fixed-point samples equivalent to those measured by using MVDI sampling techniques. Instantaneous constituent loads were developed by using velocity and discharge data obtained from ADCPs and constituent concentrations obtained from MVDI and fixed-point samples.

Michigan↗

USGS2019

ModelMuse Version 4: A graphical user interface for MODFLOW 6

ModelMuse, a graphical user interface for groundwater-modeling software, was modified to support MODFLOW 6. ModelMuse works with two types of spatial discretization in MODFLOW 6: structured grids (DIS) and discretization by vertices (DISV). Quadtree refinement is used to generate a DISV model from a structured-grid model. The locations and weights for ghost nodes used to improve DISV model accuracy are computed automatically by ModelMuse using a new algorithm. ModelMuse does not support other types of DISV grids and unstructured grids. ModelMuse supports options in MODFLOW 6 that designate individual cells as confined or convertible and remove inactive cells associated with discontinuous layers, thereby reducing the computational burden. ModelMuse can specify fully three-dimensional (3D), spatially variable anisotropy in hydraulic conductivity. Although MODFLOW 6 does not support the parameters supported by MODFLOW–2005, ModelMuse provides backward compatibility by allowing ModelMuse parameters to specify scale-factor variables in MODFLOW 6 time-series files within packages that support time-series files. ModelMuse can automatically convert the data for many of the packages from other MODFLOW models to the new data for these packages in MODFLOW 6. Some packages, such as the Streamflow-Routing (SFR) package, changed significantly enough that only a partial conversion is possible. Head and flow observations in older models are also converted to observation locations in the MODFLOW 6 Observation utility. ModelMuse accommodates the ability of MODFLOW 6 to store specific discharge components by allowing the user to visualize the components of a specific discharge on the model grid. ModelMuse supports the versions of MODPATH and ZONEBUDGET compatible with MODFLOW 6.

Scientific Investigations Report↗

USGS2021

Hydrogeology and model-simulated groundwater availability in the Salt Fork Red River aquifer, southwestern Oklahoma, 1980–2015

The 1973 Oklahoma Water Law (82 OK Stat § 82-1020.5) requires that the Oklahoma Water Resources Board (OWRB) conduct hydrologic investigations of the State’s groundwater basins to support a determination of the maximum annual yield for each groundwater basin (hereinafter referred to as an “aquifer”). The maximum annual yield allocated per acre of land is known as the equal-proportionate-share (EPS) pumping rate. At present (2021), the OWRB has not yet established a maximum annual yield and EPS pumping rate for the Salt Fork Red River aquifer. To provide updated information to the OWRB that could support evaluation and determination of an appropriate maximum annual yield, the U.S. Geological Survey (USGS), in cooperation with the OWRB, conducted a hydrologic investigation and evaluated the effects of potential groundwater withdrawals on groundwater availability in the Salt Fork Red River aquifer. The Salt Fork Red River aquifer in Greer, Harmon, and Jackson Counties of southwestern Oklahoma is composed of about 274.5 square miles of alluvium and terrace deposits associated with the Salt Fork Red River. The mean annual recharge rate to the Salt Fork Red River aquifer for the period 1980–2015 was estimated to be about 2.94 inches per year, or 10.0 percent of the mean annual precipitation for the same period (29.4 inches per year). This 1980–2015 mean annual recharge rate is equivalent to a mean annual recharge rate of about 38,000 acre-feet per year (acre-ft/yr) for the Salt Fork Red River aquifer excluding about 19,764 acres comprising the Mulberry Creek and Horse Creek terraces. The mean annual recharge rates upgradient and downgradient from USGS streamgage 07300500 Salt Fork Red River at Mangum, Okla. (hereinafter referred to as the “Mangum gage”), apportioned by aquifer area (41.5 and 58.5 percent, respectively), were about 16,000 and 22,000 acre-ft/yr, respectively. Mean annual groundwater use for the study period (1980–2015) was 3,532.7 acre-ft/yr; about 77 percent of that groundwater use was for irrigation, and about 23 percent was for public supply. Most groundwater use for irrigation was associated with wells in the Martha terrace. A hydrogeologic framework was developed for the Salt Fork Red River aquifer and included a definition of the aquifer extent and potentiometric surface, as well as a description of the textural and hydraulic properties of aquifer materials. The hydrogeologic framework was used in the construction of the numerical groundwater-flow model of the Salt Fork Red River aquifer described in this report. A conceptual model for the Salt Fork Red River aquifer that reasonably represents the groundwater-flow system was developed to constrain the construction and calibration of the numerical model. The conceptual-model water budget estimated mean annual inflows to, and outflows from, the Salt Fork Red River aquifer for the period 1980–2015 and included a subaccounting of mean annual inflows and outflows for the portions of the aquifer that were upgradient and downgradient from the Mangum gage. The numerical groundwater-flow model of the Salt Fork Red River aquifer was constructed by using MODFLOW-2005 with the Newton formulation solver. The model of the Salt Fork Red River aquifer was spatially discretized into 1,050 rows, 1,125 columns, about 170,000 active cells measuring 200 by 200 feet (ft), and a single convertible layer. The model was temporally discretized into 432 monthly transient stress periods (each with two time steps to improve model stability). An initial steady-state stress period represented mean annual inflows to, and outflows from, the aquifer and produced a solution that was used as the initial condition for subsequent transient stress periods as well as some groundwater-availability scenarios. The model was calibrated to water-table-altitude observations at selected wells and base-flow observations at selected streamgages. The simulated saturated thickness of the Salt Fork Red River aquifer was determined by subtracting the altitude of the aquifer base from the simulated water-table altitude at the end of the numerical-model period (2015). The simulated saturated thickness was more than 75 ft in a paleochannel in the Dodson terrace near the Texas border. The mean aquifer thickness (sum of saturated and unsaturated) was 49.62 ft, and the mean saturated thickness was 28.55 ft. A simulated mean transmissivity of 1,024 feet squared per day was computed from the calibrated hydraulic conductivity and saturated thickness of each cell. The simulated available water in storage at the end of the numerical-model period (2015) was 526,117 acre-feet (acre-ft); about 42 percent of that total was available upgradient from the Mangum gage, and about 58 percent of that total was available downgradient from the Mangum gage (including the Mangum terrace). Three types of groundwater-availability scenarios were run using the calibrated numerical model. These scenarios were used to (1) estimate the EPS pumping rate that ensures a minimum 20-, 40-, and 50-year life of the aquifer, (2) quantify the potential effects of projected well withdrawals on groundwater storage over a 50-year period, and (3) simulate the potential effects of a hypothetical 10-year drought on base flow and groundwater storage. The 20-, 40-, and 50-year EPS pumping rates under normal recharge conditions were about 0.51, 0.48, and 0.48 acre-foot per acre per year, respectively. Given the 155,929-acre modeled aquifer area, these rates correspond to annual yields of about 78,800, 74,900, and 74,700 acre-ft/yr, respectively. For the 20-year EPS scenario, decreasing and increasing recharge by 10 percent resulted in a 6-percent change in the EPS pumping rate in both cases; for the 40- and 50-year EPS scenarios, decreasing and increasing recharge by 10 percent resulted in a 7-percent change in the EPS pumping rate in both cases. Projected 50-year pumping scenarios were used to simulate the effects of selected well withdrawal rates on groundwater storage of the Salt Fork Red River aquifer and base flows in the Salt Fork Red River. The effects of well withdrawals were evaluated by quantifying differences in groundwater storage and base flow in four 50-year scenarios, which applied (1) no groundwater pumping, (2) mean pumping rates for the study period (1980–2015), (3) 2015 pumping rates, and (4) increasing demand pumping rates at simulated wells. The increasing demand pumping rates assumed a cumulative 20.4-percent increase in pumping over 50 years based on 2010–60 demand projections for southwestern Oklahoma. Groundwater storage after 50 years with no pumping was 535,000 acre-ft, or 8,900 acre-ft (1.7 percent) greater than the initial groundwater storage; this groundwater storage increase is equivalent to a mean water-table-altitude increase of 0.48 ft. Groundwater storage after 50 years of pumping at the mean rate for the study period (1980–2015) was 519,900 acre-ft, or 6,200 acre-ft (1.2 percent) less than the initial groundwater storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.34 ft. Groundwater storage at the end of the 50-year period with 2015 pumping rates was 513,100 acre-ft, or 13,000 acre-ft (2.5 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.71 ft. Groundwater storage at the end of the 50-year period with increasing demand pumping rates was 509,700 acre-ft, or 16,500 acre-ft (3.1 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.89 ft. A hypothetical 10-year drought scenario was used to simulate the effects of a prolonged period of reduced recharge on groundwater storage. The period January 1983–December 1992 was chosen as the simulated drought period. Drought effects were quantified by comparing the results of the drought scenario to those of the calibrated numerical model (no drought) at the end of the simulated drought period (1992). To simulate the hypothetical drought, recharge in the calibrated numerical model was reduced by 50 percent during the simulated drought period (1983–92). Upstream inflows from the Salt Fork Red River, Turkey Creek, and Bitter Creek were reduced by 75 percent. Groundwater storage at the end of the drought period (1992) was 479,200 acre-ft, or 53,200 acre-ft (10.0 percent) less than the groundwater storage of the calibrated numerical model at the end of the drought period. This decrease in groundwater storage is equivalent to a mean water-table-altitude decline of 2.9 ft. At the end of the 10-year hypothetical drought period, simulated base flows at the Mangum gage and USGS streamgage 07301110 Salt Fork Red River near Elmer, Okla., had decreased by about 80 and 70 percent, respectively.

Oklahoma↗

USGS2021

Simulation of groundwater budgets and travel times for watersheds on the north shore of Long Island Sound, with implications for nitrogen-transport studies

Aquatic systems in and around the Long Island Sound (LIS) provide a variety of ecological and economic benefits, but in some areas of the LIS, aquatic ecosystems have become degraded by excess nitrogen. A substantial fraction of the nitrogen inputs to the LIS are transported through the groundwater-flow system. Because groundwater travel times in surficial aquifers can exceed 100 years, multiyear lags are introduced between inputs at the water table in recharge areas and discharge to inland or coastal receiving waters. The U.S. Geological Survey, in cooperation with the Connecticut Department of Energy and Environmental Protection and the U.S. Environmental Protection Agency’s Long Island Sound Study, developed a steady-state groundwater model of the watersheds draining from the northern shore of the LIS for the purpose of calculating groundwater budgets and travel times to coastal waters. The model was developed by using the MODFLOW–NWT software and existing spatial data on aquifers, river networks, land-surface altitudes, land cover, groundwater recharge, and water use. Coastal waters were delineated on the basis of the National Wetland Inventory; all non-coastal waters were collectively termed “inland waters.” A coarse-resolution model was calibrated by using the PEST++ software, long-term records of water levels in 65 wells, stream altitudes from 477 streams, base-flow records for 14 streamgages that are relatively unaffected by withdrawals, and error metrics based on incorrectly simulated flooding and incorrectly simulated dry streams. The calibrated values were used in a fine-resolution model in which the mean absolute residuals were 4.5 meters for groundwater levels, 1.3 meters for stream altitudes, and 7,200 cubic meters per day (2.9 cubic feet per second) for base flow. About 89 percent of the terrestrial cells were correctly simulated with the water table below land surface, and nearly 90 percent of the cells representing streams were correctly simulated as having the water table above the stream bottom. Together, these metrics suggest that this model is robust for simulating regional-scale groundwater patterns. Simulated groundwater budgets were compiled for the entire study area, for each HUC12 (Hydrologic Unit Code no. 12) watershed and its adjacent coastal waters, if applicable, within the study area, and for 14 coastal-embayment watersheds. Most groundwater (90.6 percent of inflows) discharged to inland waters, with smaller fractions to coastal waters (7.0 percent) and well withdrawals (2.4 percent). When computed for HUC12 watersheds with coastal discharge, the portions of groundwater discharging to coastal waters ranged from 0.02 to 66 percent of groundwater outflows, with a median of 13 percent. Within priority-embayment watersheds, the portions of groundwater discharging to coastal waters ranged from 2 to 56 percent, with a median of 15 percent. Groundwater travel times also were simulated for the entire study area, for each HUC12 watershed and its adjacent coastal waters, if applicable, within the study area and for 14 priority coastal embayments. Within the entire study area, the median groundwater travel time was 1.9 years, with an interquartile range of 0.1 to 5.9 years. Sensitivity analysis of groundwater travel times within a subbasin in the study area indicates that the travel times are a function of the grid resolution, with coarser grids resulting in shorter median travel times. Travel times for groundwater discharging to coastal waters were similar to travel times for groundwater discharging to inland waters, with a median of 1.9 years. Median travel times for the HUC12 watersheds ranged from 0.9 to 53.5 years, with a median of 1.8 years. Among HUC12 watersheds that include coastal areas, travel times for groundwater discharging to coastal waters ranged from less than 1 to 61.6 years, with a median of 2.8 years. The HUC12 watersheds with the longest simulated travel times were in the urban area near New York City where the model performance is less accurate. Median travel times for groundwater discharging to coastal waters within the priority-embayment watersheds ranged from less than 1 to 18.6 years, with a median of 2.3 years. A more focused analysis was conducted for the Niantic River watershed to demonstrate the applicability of the regional model to local-scale nitrogen-transport analyses by using nitrogen-input and -attenuation rates from literature sources. Nitrogen inputs were estimated by using land-cover-based loading factors, and attenuation was estimated by using attenuation factors based on geologic zones and soil properties. Based on this analysis, groundwater transports an estimated 22,000 kilograms of nitrogen per year (2.9 kilograms of nitrogen per hectare per year) to streams, rivers, and coastal waters within the Niantic River watershed. Approximately 36 percent of discharging nitrogen is from atmospheric-deposition sources, 38 percent is from fertilizers, and 26 percent is from septic systems. Most of the groundwater-transported nitrogen (88 percent) discharges first to streams and rivers, with only 12 percent discharging directly to coastal waters. Travel times for groundwater-transported nitrogen ranged from less than 1 day to more than 100 years, with a median of 1.6 years.

Connecticut, New York, Rhode Island↗

USGS2024

Groundwater model of the Harney Basin, southeastern Oregon

Groundwater development, mainly for large-scale irrigation, has increased substantially in the Harney Basin of southeastern Oregon since 2010. Concurrently, some areas of the basin experienced groundwater-level declines of more than 100 feet, and some shallow wells have gone dry. The Oregon Water Resources Department has limited new groundwater development in the basin until an improved understanding of the groundwater-flow system is available. The groundwater resources report by Gingerich and others (2022, U.S. Geological Survey Scientific Investigations Report 2021–5103, https://doi.org/10.3133/sir20215103 ) provides that understanding. This report describes the development of a numerical groundwater-flow model that can be used as a tool to help improve that understanding. The Harney Basin Groundwater Model was developed using the finite-difference groundwater-modeling software U.S. Geological Survey modular finite-difference groundwater-flow model (MODFLOW 6) and associated Python pre- and post-processing routines. The groundwater model encompasses the entire 5,240-square-mile Harney Basin and adjacent areas and is calibrated to the hydrologic conditions from 1930 to 2018. The model has a uniform grid consisting of 78,064 nearly square cells, each covering 2,005 by 2,007 feet (about 92 acres) and has 10 layers (780,640 total cells) representing the vertical distribution of hydrogeologic units. The results from the calibrated model simulations indicate that groundwater pumpage exceeded recharge since about the mid-1980s, resulting in an estimated net cumulative depletion of groundwater storage (discharge minus recharge) of about 840,000 acre-feet and also indicated declines in groundwater evapotranspiration and spring and stream discharge. Model simulations show as much as 100 feet of groundwater-level decline in some areas and more than 40 feet of decline in widespread areas in recent decades. Model simulations are consistent with field observations of groundwater levels through time.

Oregon↗