Geology Reports⌕ Search

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.

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

At least 559 records · Page 31Linked to original sources
USGS2011

Simulation of the effects of groundwater withdrawals on water-level altitudes in the Sparta aquifer in the Bayou Meto-Grand Prairie area of eastern Arkansas, 2007-37

A groundwater-flow model of the Mississippi embayment was used to evaluate changes in water-level altitudes before (scenario 1) and after (scenario 2) the addition of wells that simulate potential future pumping from the Sparta aquifer in the Bayou Meto-Grand Prairie area of eastern Arkansas for the 30-year period from 2007 through 2037. Water-level altitudes at six model cell locations from the two different scenarios were compared for the period 2007 through 2037. Potential future pumping wells were added to the Mississippi Embayment Regional Aquifer Study model at a rate of 13 wells per year within areas of potential future pumping. Change maps for the Bayou Meto-Grand Prairie area were constructed for each scenario and water-level hydrographs were constructed for each scenario for each of the six model cell locations. The additional pumping from wells in the Sparta aquifer created greater water-level declines in the Bayou Meto-Grand Prairie area. In scenario 1, simulated water-level altitude declines range from 20 to 40 feet from 2007 through 2037. In scenario 2, the cone of depression in Lonoke County is the deepest, with a maximum water-level decline of approximately 102 feet. Water-level altitude declines range from 40 to 50 feet over most of the remainder of the Bayou Meto-Grand Prairie area in scenario 2. Simulated water-level altitudes across the Bayou Meto-Grand Prairie area and at all six model cell locations indicate substantial declines when additional wells pumping from the Sparta aquifer are introduced into the model from 2007 through 2037.

Arkansas↗
USGS2012

Numerical simulation of groundwater movement and managed aquifer recharge from Sand Hollow Reservoir, Hurricane Bench area, Washington County, Utah

The Hurricane Bench area of Washington County, Utah, is a 70 square-mile area extending south from the Virgin River and encompassing Sand Hollow basin. Sand Hollow Reservoir, located on Hurricane Bench, was completed in March 2002 and is operated primarily as a managed aquifer recharge project by the Washington County Water Conservancy District. The reservoir is situated on a thick sequence of the Navajo Sandstone and Kayenta Formation. Total recharge to the underlying Navajo aquifer from the reservoir was about 86,000 acre-feet from 2002 to 2009. Natural recharge as infiltration of precipitation was approximately 2,100 acre-feet per year for the same period. Discharge occurs as seepage to the Virgin River, municipal and irrigation well withdrawals, and seepage to drains at the base of reservoir dams. Within the Hurricane Bench area, unconfined groundwater-flow conditions generally exist throughout the Navajo Sandstone. Navajo Sandstone hydraulic-conductivity values from regional aquifer testing range from 0.8 to 32 feet per day. The large variability in hydraulic conductivity is attributed to bedrock fractures that trend north-northeast across the study area. A numerical groundwater-flow model was developed to simulate groundwater movement in the Hurricane Bench area and to simulate the movement of managed aquifer recharge from Sand Hollow Reservoir through the groundwater system. The model was calibrated to combined steady- and transient-state conditions. The steady-state portion of the simulation was developed and calibrated by using hydrologic data that represented average conditions for 1975. The transient-state portion of the simulation was developed and calibrated by using hydrologic data collected from 1976 to 2009. Areally, the model grid was 98 rows by 76 columns with a variable cell size ranging from about 1.5 to 25 acres. Smaller cells were used to represent the reservoir to accurately simulate the reservoir bathymetry and nearby monitoring wells; larger cells were used in the northern and southern portions of the model where water-level data were limited. Vertically, the aquifer system was divided into 10 layers, which incorporated the Navajo Sandstone and Kayenta Formation. The model simulated recharge to the groundwater system as natural infiltration of precipitation and as infiltration of managed aquifer recharge from Sand Hollow Reservoir. Groundwater discharge was simulated as well withdrawals, shallow drains at the base of reservoir dams, and seepage to the Virgin River. During calibration, variables were adjusted within probable ranges to minimize differences among model-simulated and observed water levels, groundwater travel times, drain discharges, and monthly estimated reservoir recharge.

Utah↗
USGS2013

Potential depletion of surface water in the Colorado River and agricultural drains by groundwater pumping in the Parker-Palo Verde-Cibola area, Arizona and California

Water use along the lower Colorado River is allocated as “consumptive use,” which is defined to be the amount of water diverted from the river minus the amount that returns to the river. Diversions of water from the river include surface water in canals and water removed from the river by pumping wells in the aquifer connected to the river. A complication in accounting for water pumped by wells occurs if the pumping depletes water in drains and reduces measured return flow in those drains. In that case, consumptive use of water pumped by the wells is accounted for in the reduction of measured return flow. A method is needed to understand where groundwater pumping will deplete water in the river and where it will deplete water in drains. To provide a basis for future accounting for pumped groundwater in the Parker-Palo Verde-Cibola area, a superposition model was constructed. The model consists of three layers of finite-difference cells that cover most of the aquifer in the study area. The model was run repeatedly with each run having a pumping well in a different model cell. The source of pumped water that is depletion of the river, expressed as a fraction of the pumping rate, was computed for all active cells in model layer 1, and maps were constructed to understand where groundwater pumping depletes the river and where it depletes drains. The model results indicate that if one or more drains exist between a pumping well location and the river, nearly all of the depletion will be from drains, and little or no depletion will come from the Colorado River. Results also show that if a well pumps on a side of the river with no drains in the immediate area, depletion will come from the Colorado River. Finally, if a well pumps between the river and drains that parallel the river, a fraction of the pumping will come from the river and the rest will come from the drains. Model results presented in this report may be considered in development or refinement of strategies for accounting for groundwater pumping in the river aquifer connected to the Colorado River in the study area.

Arizona;California↗
USGS2015

Status and understanding of groundwater quality in the Cascade Range and Modoc Plateau study unit, 2010: California GAMA Priority Basin Project

Groundwater quality in the Cascade Range and Modoc Plateau study unit was investigated as part of the California State Water Resources Control Board’s Groundwater Ambient Monitoring and Assessment (GAMA) Program Priority Basin Project. The study was designed to provide a statistically unbiased assessment of untreated groundwater quality in the primary aquifer system. The depth of the primary aquifer system for the Cascade Range and Modoc Plateau study unit was delineated by the depths of the screened or open intervals of wells in the State of California’s database of public-supply wells. Two types of assessments were made: a status assessment that described the current quality of the groundwater resource, and an understanding assessment that made evaluations of relations between groundwater quality and potential explanatory factors representing characteristics of the primary aquifer system. The assessments characterize the quality of untreated groundwater, not the quality of treated drinking water delivered to consumers by water distributors. The status assessment was based on water-quality data collected in 2010 by the U.S. Geological Survey from 90 wells and springs (USGS-grid wells) and on water-quality data compiled from the State of California’s regulatory compliance database for samples collected from 240 public-supply wells between September 2007 and September 2010. To provide context, the water-quality data discussed in this report were compared to California and Federal drinking-water regulatory and non-regulatory benchmarks for treated drinking water. Groundwater quality is defined in terms of relative concentrations (RCs), which are calculated by dividing the concentration of a constituent in groundwater by the concentration of the benchmark for that constituent. The RCs for inorganic constituents (major ions, trace elements, nutrients, and radioactive constituents) were classified as “high” (the RC is greater than 1.0, indicating that the concentration is above the benchmark), “moderate” (the RC is from 1.0 to greater than 0.5), or “low” (the RC is less than or equal to 0.5). For organic constituents (volatile organic compounds and pesticides) and special-interest constituents (perchlorate), the boundary between moderate and low RCs was set at 0.1. All benchmarks used for organic constituents were health-based. For inorganic constituents, health-based and aesthetic-based benchmarks were used. Constituents without benchmarks were not considered in the status assessment . The primary metric used for quantifying regional-scale groundwater quality was the aquifer-scale proportion—the areal percentages of the primary aquifer system with high, moderate, and low RCs for a given constituent or class of constituents. The study unit was divided into six study areas on the basis of geologic differences (Eastside Sacramento Valley, Honey Lake Valley groundwater basin, Cascade Range and Modoc Plateau Low Use Basins, Quaternary Volcanic Areas, Shasta Valley and Mount Shasta Volcanic Area, and Tertiary Volcanic Areas), and each study area was divided into equal-area grid cells. Aquifer-scale proportions were calculated for individual constituents and constituent classes for each of the six study areas and for the study unit as a whole by using grid-based (one well per cell) and spatially weighted (many wells per cell) statistical methods. The status assessment showed that inorganic constituents were present at high and moderate RCs in greater proportions of the Cascade Range and Modoc Plateau study unit than were organic constituents. One or more inorganic constituents with health-based benchmarks were present at high RCs in 9.4 percent, and at moderate RCs in 14.7 percent of the primary aquifer system. Arsenic was present at high RCs in approximately 3 percent of the primary aquifer system; boron, molybdenum, uranium, and vanadium each were present at high RCs in approximately 2 percent of the primary aquifer system. One or more inorganic constituents with aesthetic-based benchmarks were present at high RCs in 15.1 percent of the primary aquifer system and at moderate RCs in 4.9 percent. Manganese, iron, and total dissolved solids were present at high RCs in approximately 12 percent, 5 percent, and 2 percent, respectively, of the primary aquifer system. Organic constituents were not detected at high or moderate RCs in the primary aquifer system, and one or more organic constituents were detected at low RCs in approximately 40 percent of the primary aquifer system. Two classes of organic constituents were detected in more than 10 percent of the primary aquifer system: trihalomethanes (chloroform only) and herbicides. The special interest constituent perchlorate was not detected at high RCs, but was detected at moderate RCs in approximately 2 percent of the primary aquifer system. The understanding assessment relied on statistical tests to evaluate relations between concentrations of constituents and values of potential explanatory factors representing geology, land use, well construction, hydrologic conditions, groundwater age, and geochemical conditions. The majority of the high and moderate RCs of arsenic, boron, molybdenum, uranium, and total dissolved solids were in samples from the Honey Lake Valley groundwater basin study area. Groundwater mixing with hydrothermal fluids present in the study area, evaporative concentration of groundwater in the Honey Lake playa, presence of uranium-bearing sediment derived from the adjacent Sierra Nevada, and release of arsenic and other trace elements from sediments under high pH and low dissolved oxygen conditions all appeared to contribute to these elevated concentrations. Thermal springs are in many parts of the Cascade Range and Modoc Plateau study unit and could account for locally elevated concentrations of arsenic, boron, molybdenum, and total dissolved solids in samples from the other study areas. Vanadium concentrations were greater in oxic samples than in anoxic samples, but were not correlated with pH, contrary to expectations from previous studies. Organic constituents were not detected at high or moderate RCs, and the occurrence of low organic constituents at low RCs ranged from 27 percent to 73 percent of the primary aquifers system in the six study areas. The Shasta Valley and Mount Shasta Volcanic study area had significantly greater occurrence of low RCs of herbicides compared to all of the other study areas, which could reflect the greater prevalence of modern groundwater in the Shasta Valley and Mount Shasta Volcanic study area and the presence of potential sources of herbicides, including applications to timberlands and roadside rights-of-way. The Eastside Sacramento Valley study area had the greatest occurrence of low concentrations of chloroform, and chloroform occurrence was most strongly associated with the combination of septic-tank density greater than two tanks per square kilometer and urban land use greater than 10 percent within a radius of 500 meters of the well. These conditions were most prevalent in the Eastside Sacramento Valley study area. The detection frequency of low concentrations of perchlorate was consistent with the probability of occurrence expected under natural conditions, except in the Eastside Sacramento Valley study area, where detection frequencies were much higher than expected and could not be explained by known anthropogenic sources of perchlorate.

California↗
USGS2018

Geochemical conditions and nitrogen transport in nearshore groundwater and the subterranean estuary at a Cape Cod embayment, East Falmouth, Massachusetts, 2013–14

Nitrogen transport and transformation were studied during 2013 to 2014 by the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, in a subterranean estuary beneath onshore locations on the Seacoast Shores peninsula, a residential area in Falmouth, Massachusetts, served by septic systems and cesspools, and adjacent offshore locations in the Eel River, a saltwater embayment connected to the ocean. The field investigation included installation and sampling of clusters of wells and temporary sampling points near a transect extending from about 35 meters (m) onshore to 18 m offshore. The fresh groundwater at the study site formed a lens about 11 m thick at the shoreline that was underlain by saline groundwater. Groundwater flow in the water-table aquifer was oriented northwestward toward the embayment. Nitrate concentrations in the fresh groundwater at a site about 35 m onshore increased in the downward direction from less than 500 micromoles per liter near the water table to about 1,700 micromoles per liter just above the freshwater/saltwater transition zone. Dissolved oxygen was largely absent in the onshore fresh groundwater. Distributions of salinity, dissolved oxygen, and nitrate at the shoreline and offshore generally were similar to those onshore; at some locations, however, shallow saline water was present above the freshwater, and there were scattered occurrences of elevated dissolved oxygen concentrations. Geochemical indicators of nitrate reduction, including concentrations of the reaction product nitrogen gas, stable isotope ratios of nitrate and nitrogen gas, and changes in alkalinity, provided evidence for nitrate reduction in two zones separated vertically by a zone 7–8 m thick with no evidence of nitrate reduction. The shallow nitrate-reduction zone was near the water table in fresh groundwater onshore, where nitrate reduction may be related to particular recharge conditions at nearby sources. The shallow nitrate-reduction zone also may be related to an interval of fine-grained sediments at about the same altitude (−1 to −6 m relative to the National Geodetic Vertical Datum of 1929), where flow is slower and reactive electron donors such as solid organic carbon, iron, or sulfide phases may be present to drive the reduction. The deep nitrate-reduction zone was near the freshwater/saltwater transition zone, where nitrate reduction may be related to mixing of freshwater containing nitrate and saltwater containing dissolved organic carbon and ammonium, or to fine-grained sediments near the transition zone. The maximum amount of nitrate converted to nitrogen gas was estimated to be less than or equal to 300 micromoles per liter in both nitrate-reduction zones. The presence of nitrate and low dissolved oxygen concentrations in the 7–8-meter-thick zone between the shallow and deep nitrate-reduction zones are conditions that could permit nitrate reduction. The absence of evidence of nitrate reduction in the high-nitrate zone may have resulted from the lack of reactive electron donors in that depth interval. The high-nitrate zone dissipated somewhat in the offshore direction, but the current study did not extend far enough to encompass the fresh groundwater discharge area or determine how much of the nitrate was removed prior to discharge. A shallow intertidal saltwater cell was formed during a spring tide by saltwater infiltration during tidal run-up on the beach. Nitrate reduction might have occurred if nitrate-containing fresh groundwater discharging to the estuary mixed with the saltwater containing dissolved organic carbon in this zone, but samples collected from the intertidal saltwater cell during this study were not analyzed for indicators of nitrate reduction. Elevated dissolved oxygen concentrations in fresh groundwater 9 m offshore may indicate that groundwater flow was partly oblique to the sampling transect or that groundwater from a regional flow system was converging under the river near the study area. Flow directions also may have been affected by aquifer heterogeneity such as the shallow fine-grained sediments onshore and at the bottom of the Eel River. Improved understanding of the fate of nitrate in this type of complex setting might be gained by including additional characterization of aquifer heterogeneity and groundwater flow and extending investigations of nitrate reduction to the shallow sediments in the intertidal saltwater cell and adjacent subtidal zone and to locations farther offshore beneath the estuary.

Massachusetts↗
USGS2018

Revised groundwater-flow model of the glacial aquifer system north of Aberdeen, South Dakota, through water year 2015

The city of Aberdeen, in northeastern South Dakota, requires an expanded and sustainable supply of water to meet current and future demands. Conceptual and numerical models of the glacial aquifer system in the area north of Aberdeen were developed by the U.S. Geological Survey in cooperation with the City of Aberdeen in 2012. The U.S. Geological Survey, in cooperation with the City of Aberdeen, completed a study to revise the original numerical groundwater-flow model using data through water year (WY) 2015 to aid the City of Aberdeen in their development of plans and strategies for a sustainable water supply and to increase understanding of the glacial aquifer system and groundwater-flow system near Aberdeen. The original model was revised to improve the fit between model-simulated values and observed (measured or estimated) data, provide greater insight into surface-water interactions, and improve the usefulness of the model for water-supply planning. The revised groundwater-flow model (hereafter referred to as the “revised model”) presented in this report supersedes the original model. The purpose of this report is to describe a revised groundwater-flow model including data collection, model calibration, and model results for the glacial aquifer system including the Elm, Middle James, and Deep James aquifers north of Aberdeen, South Dakota, using updated hydrologic data through WY 2015. The original numerical model was revised in several ways. The model was modified by adding four new layers, which included a surficial layer, two intervening confining layers, and a shale bedrock layer. The revised model provides an improved understanding of the groundwater-flow system in comparison to the original model. The principal aquifers of the model area include portions of the Elm, Middle James, and Deep James aquifers. The lithologic information used to define and describe the aquifers in the model area was unaltered; however, aquifer properties and boundary conditions were reviewed and updated using geological information reported by the South Dakota Department of Environmental and Natural Resources and information obtained from geophysical investigations for this study. The horizontal extent of the Elm, Middle James, and Deep James aquifers was unaltered from the original model. The thickness of the Deep James aquifer was modified based on interpretations from the geophysical investigations. In general, groundwater in the Elm aquifer flowed from northwest to southeast and locally towards rivers and streams. Similarly, in the Middle James and Deep James aquifers, groundwater also typically flowed southeast. The revisions made to the original model include use of the following MODFLOW stress packages: Recharge, Evapotranspiration, Time-Variant Specified Head, Wells, Drains, and Stream Flow Routing, all of which were updated from the original model except for the Stream Flow Routing Package, which replaced the River Package used in the original model. Model calibration is the process of estimating model parameters to minimize the differences, or residuals, between observed data and simulated values; therefore, Parameter ESTimation (PEST) software was used to optimize model input parameters by matching model-simulated values to observed data. Calibration parameters included horizontal hydraulic conductivity, vertical hydraulic conductivity, specific yield, specific storage, and vertical streambed conductance for stream and drain cells. Multipliers were used to calibrate the recharge and evapotranspiration stresses. Evapotranspiration extinction depth also was adjusted during model calibration. Comparisons to the original model are described to highlight the changes made in the revised model. In general, the revised model adequately simulates the natural system and compares favorably with observed hydrologic data. Simulated water levels were evaluated by comparing them to single water-level observations at selected well locations. The selected wells were the same wells used in the original model. The coefficient of determination value between simulated and observed water levels for the revised model was 0.89 and included simulated and observed values from October 1, 1974 (WY 1975), through September 30, 2015 (WY 2015). The coefficient of determination value for the original model was 0.94 and included simulated and observed values from October 1, 1974, through September 30, 2009. The difference may indicate that the original model could have been overfit to hydraulic head observations because base flow was not simulated. The additional data used in the revised model included some climatically wetter, more extreme periods, such as 2011, in which annual precipitation was 30.9 inches. Average annual precipitation for the original model timeframe, which included data from WYs 1975–2009, was 20.26 inches. Additional precipitation data for WYs 2010–15, included in the revised model timeframe, resulted in an average annual precipitation for WYs 1975–2015 in the model area of 20.6 inches. The larger variability in climate data coupled with the additional water-level data could explain the lower coefficient of determination for water levels in the revised model. The revised model was used to calculate various groundwater-budget components for steady-state and transient conditions for WYs 1975–2015. The time-variant specified-head cells in the revised model had the largest change when compared to the original steady-state model for inflows and outflows. Comparing the transient budget components between the original and the revised models indicated that inflow from recharge and time-variant specified-head cells had the greatest effect on groundwater inflows, and outflow from storage had the greatest effect on groundwater outflows. The simulated potentiometric contours from the revised model were compared with (1) the observed (interpreted) potentiometric surface (layer 2) and the hydraulic head values (layers 4 and 6) and (2) the simulated contours from the original model. The simulated hydraulic gradients and general direction of groundwater flow in the Elm aquifer in the revised model generally matched the observed potentiometric contours, the simulated potentiometric contours from the original model, and general flow directions interpreted to be perpendicular to the contours. Minor discrepancies between simulated potentiometric contours from the revised model and the observed potentiometric contours may be due to the lack of observed data in the model area. The revised model was designed to reduce the limitations of the original model. The revisions were validated by comparing the results of the original model with the revised model. A primary benefit of the revised model is the inclusion of the surficial deposits and the confining units as explicit layers in the model. The addition of the surficial layer was beneficial for three primary reasons: (1) more accurate representation of recharge from precipitation, (2) more accurate representation of groundwater evapotranspiration, and (3) more accurate representation of groundwater and surface-water interactions. The groundwater model is a numeric approximation of a complex physical hydrologic system, and the revised model data were interpolated in regions with sparse data. Additionally, model discretization included averaged and interpolated values for water use, withdrawal rates, and hydraulic conductivity. The revised model provides a useful estimate for hydraulic gradients, groundwater-flow directions, and aquifer response to groundwater withdrawals.

South Dakota↗
USGS2019

Potential interaction of groundwater and surface water including autonomous underwater vehicle reconnaissance at Nolin River Lake, Kentucky, 2016

The U.S. Geological Survey collaborated with the U.S. Army Corps of Engineers, Louisville District, on a synoptic study of water quality at Nolin River Lake during August 2016. The purpose of the study was to develop a better understanding of the potential for interaction between groundwater and surface water at Nolin River Lake, Kentucky. Groundwater can have properties that are measurably different from those in adjacent surface water, and inflows and outflows can be an important component of water quality and quantity. An improved understanding of potential interaction of groundwater and surface water at Nolin River Lake may be used to refine lake-management strategies. This study (1) compiled and interpreted existing information to characterize the hydrogeological setting and implications for potential interaction of groundwater and surface water in the Nolin River Lake watershed; (2) collected transects of onsite water-quality parameters using an autonomous underwater vehicle (AUV) in areas with potential for interaction of groundwater and surface water, including five sites on Nolin River Lake and one site on the Nolin River; and (3) collected discrete water-quality and phytoplankton community data at the same six sites. A review of existing hydrogeologic information did not indicate the presence of karst features adjacent to or beneath Nolin River Lake that would facilitate groundwater interaction with the reservoir. Observations leading to this conclusion include (1) limestone that is adjacent to the shoreline and perhaps beneath the lake, is overlain with siliciclastic rocks and fine-grained sediment that inhibits infiltration and development of karst features that encourage rapid groundwater flow; (2) the geologic deposits surrounding the reservoir are described as having limited or no potential for development of karst features, some exceptions may exist in tributary valleys; (3) very few karst features were mapped within 1 mile of the reservoir or in the area currently occupied by the reservoir; and (4) faults that intersect the reservoir but may not possess hydraulic properties that cause the faults to be conduits for groundwater flow. Groundwater interaction with reservoir tributaries is likely more common in areas of the watershed upstream from Nolin River Lake where karst hydrogeology is prevalent. Results of water-quality surveys using an AUV from August 15 to 19, 2016, did not identify areas of anomalous values that might indicate groundwater inflows through preferential flow zones. Spatial distributions of water-quality parameters were generally uniform within each constant-depth layer. The constant-depth layers were selected to be above, within, and below the thermocline and ranged from the water surface to 25 feet. Surveys near the bottom of the reservoir that might have been more sensitive to groundwater inflows were not done because presurvey data were not available to indicate locations of obstacles that could ensnare the AUV. Water-quality data collected with the AUV did identify water-quality anomalies where stream tributaries were discharging to the reservoir. The discrete water-quality samples indicated uniformity among the five reservoir sites. The riverine site that is immediately upstream from Nolin River Lake, however, had some unique water-quality characteristics relative to sites on the reservoir. The highest concentrations of nitrate plus nitrite as nitrogen (0.145 milligrams per liter [mg/L]), total phosphorous (0.07 mg/L), chlorophyll a (36.1 micrograms per liter), and pheophytin a (10.2 micrograms per liter) were measured at the Nolin River Lake riverine site (site 2NRR20034). The concentrations of nutrients and chlorophyll a at the riverine site did exceed the 25th percentile of median concentrations measured by the U.S. Environmental Protection Agency (EPA) at other lakes and reservoirs in EPA level IV ecoregion 71a. Concentrations of most nutrients and chlorophyll a at the five reservoir sites also exceeded the 25th percentile of median concentrations in EPA level IV ecoregion 72h. The exception was the concentrations of total phosphorus as phosphorus at the reservoir sites that were at or below the 25th percentile of median concentrations measured by EPA (0.03 mg/L). Concentrations of orthophosphate as phosphorus were less than the method detection limit of 0.004 mg/L at all sites. The phytoplankton community in Nolin River Lake was almost exclusively (greater than 90 percent of total phytoplankton abundance) cyanobacteria, also known as blue-green algae. A species of Cylindrospermopsis dominated the cyanobacterial community at the five reservoir sites, while Chroococcus microscopicus was most abundant at the riverine site. Cyanobacterial cell densities ranged from 10,000 to 198,067,460 cells per liter in five areas in the reservoir and from 4,800 to 73,751,253 cells per liter at the riverine site. Multiple potential sources of water to Nolin River Lake include direct precipitation, overland flow, interflow, groundwater, and surface water. Understanding the exact contribution of each of these components to the water budget at Nolin River Lake may help the U.S. Army Corps of Engineers manage the water quality, water quantity, and biological communities in the reservoir. Additional hydrogeologic and water-quality data that builds on the results of this study may refine the inferences of this study; for example, deeper AUV surveys that target the largest fault zones might further the understanding of the potential for groundwater flow through those features. A complete understanding of the reservoir hydrology, however, may require the use of scientific methods intended for water bodies as large as Nolin River Lake, such as aerial infrared photography and imagery; water mass, chemical, and isotopic balance studies; geophysical measurements; and numerical simulations.

Kentucky↗
USGS2015

Scoops3D: software to analyze 3D slope stability throughout a digital landscape

The computer program, Scoops3D, evaluates slope stability throughout a digital landscape represented by a digital elevation model (DEM). The program uses a three-dimensional (3D) method of columns approach to assess the stability of many (typically millions) potential landslides within a user-defined size range. For each potential landslide (or failure), Scoops3D assesses the stability of a rotational, spherical slip surface encompassing many DEM cells using a 3D version of either Bishop’s simplified method or the Ordinary (Fellenius) method of limit-equilibrium analysis. Scoops3D has several options for the user to systematically and efficiently search throughout an entire DEM, thereby incorporating the effects of complex surface topography. In a thorough search, each DEM cell is included in multiple potential failures, and Scoops3D records the lowest stability (factor of safety) for each DEM cell, as well as the size (volume or area) associated with each of these potential landslides. It also determines the least-stable potential failure for the entire DEM. The user has a variety of options for building a 3D domain, including layers or full 3D distributions of strength and pore-water pressures, simplistic earthquake loading, and unsaturated suction conditions. Results from Scoops3D can be readily incorporated into a geographic information system (GIS) or other visualization software. This manual includes information on the theoretical basis for the slope-stability analysis, requirements for constructing and searching a 3D domain, a detailed operational guide (including step-by-step instructions for using the graphical user interface [GUI] software, Scoops3D-i) and input/output file specifications, practical considerations for conducting an analysis, results of verification tests, and multiple examples illustrating the capabilities of Scoops3D. Easy-to-use software installation packages are available for the Windows or Macintosh operating systems; these packages install the compiled Scoops3D program, the GUI (Scoops3D-i), and associated documentation. Several Scoops3D examples, including all input and output files, are available as well. The source code is written in the Fortran 90 language and can be compiled to run on any computer operating system with an appropriate compiler.

Techniques and Methods↗
USGS2009

The Farm Process Version 2 (FMP2) for MODFLOW-2005 - Modifications and Upgrades to FMP1

The ability to dynamically simulate the integrated supply-and-demand components of irrigated agricultural is needed to thoroughly understand the interrelation between surface water and groundwater flow in areas where the water-use by vegetation is an important component of the water budget. To meet this need, the computer program Farm Process (FMP1) was updated and refined for use with the U.S. Geological Survey's MODFLOW-2005 groundwater-flow model, and is referred to as MF2005-FMP2. The updated program allows the simulation, analysis, and management of nearly all components of human and natural water use. MF2005-FMP2 represents a complete hydrologic model that fully links the movement and use of groundwater, surface water, and imported water for water consumption of irrigated agriculture, but also of urban use, and of natural vegetation. Supply and demand components of water use are analyzed under demand-driven and supply-constrained conditions. From large- to small-scale settings, the MF2005-FMP2 has the unique set of capabilities to simulate and analyze historical, present, and future conditions. MF2005-FMP2 facilitates the analysis of agricultural water use where little data is available for pumpage, land use, or agricultural information. The features presented in this new version of FMP2 along with the linkages to the Streamflow Routing (SFR), Multi-Node Well (MNW), and Unsaturated Zone Flow (UZF) Packages prevents mass loss to an open system and helps to account for 'all of the water everywhere and all of the time'. The first version, FMP1 for MODFLOW-2000, is limited to (a) transpiration uptake from unsaturated root zones, (b) on-farm efficiency defined solely by farm and not by crop type, (c) a simulation of water use and returnflows related only to irrigated agriculture and not also to non-irrigated vegetation, (d) a definition of consumptive use as potential crop evapotranspiration, (e) percolation being instantly recharged to the uppermost active aquifer, (f) automatic routing of returnflow from runoff either to reaches of tributary stream segments adjacent to a farm or to one reach nearest to the farm's lowest elevation, (g) farm-well pumping from cell locations regardless of whether an irrigation requirement from these cells exists or not, and (h) specified non-routed water transfers from an undefined source outside the model domain. All of these limitations are overcome in MF2005-FMP2. The new features include (a) simulation of transpiration uptake from variably saturated, fully saturated, or ponded root zones (for example, for crops like rice or riparian vegetation), (b) definition of on-farm efficiency not only by farm but also by crop, (c) simulation of water use and returnflow from non-irrigated vegetation (for example, rain-fed agriculture or native vegetation), (d) use of crop coefficients and reference evapotranspiration, (e) simulation of the delay between percolation from farms through the unsaturated zone and recharge into the uppermost active aquifer by linking FMP2 to the UZF Package, (f) an option to manually control the routing of returnflow from farm runoff to streams, (g) an option to limit pumping to wells located only in cells where an irrigation requirement exists, and (h) simulation of water transfers to farms from a series of well fields (for example, recovery well field of an aquifer-storage-and-recovery system, ASR). In addition to the output of an economic budget for each farm between irrigation demand and supply ('Farm Demand and Supply Budget' in FMP1), a new output option called 'Farm Budget' was created for FMP2, which allows the user to track all physical flows into and out of a water accounting unit at all times. Such a unit can represent individual farms, farming districts, natural areas, or urban areas. The example model demonstrates the application of MF2005-FMP2 with delayed recharge through an unsaturated zone, rejected infiltration in a riparian area, changes in de

Techniques and Methods↗
USGS2013

Documentation of the seawater intrusion (SWI2) package for MODFLOW

The SWI2 Package is the latest release of the Seawater Intrusion (SWI) Package for MODFLOW. The SWI2 Package allows three-dimensional vertically integrated variable-density groundwater flow and seawater intrusion in coastal multiaquifer systems to be simulated using MODFLOW-2005. Vertically integrated variable-density groundwater flow is based on the Dupuit approximation in which an aquifer is vertically discretized into zones of differing densities, separated from each other by defined surfaces representing interfaces or density isosurfaces. The numerical approach used in the SWI2 Package does not account for diffusion and dispersion and should not be used where these processes are important. The resulting differential equations are equivalent in form to the groundwater flow equation for uniform-density flow. The approach implemented in the SWI2 Package allows density effects to be incorporated into MODFLOW-2005 through the addition of pseudo-source terms to the groundwater flow equation without the need to solve a separate advective-dispersive transport equation. Vertical and horizontal movement of defined density surfaces is calculated separately using a combination of fluxes calculated through solution of the groundwater flow equation and a simple tip and toe tracking algorithm. Use of the SWI2 Package in MODFLOW-2005 only requires the addition of a single additional input file and modification of boundary heads to freshwater heads referenced to the top of the aquifer. Fluid density within model layers can be represented using zones of constant density (stratified flow) or continuously varying density (piecewise linear in the vertical direction) in the SWI2 Package. The main advantage of using the SWI2 Package instead of variable-density groundwater flow and dispersive solute transport codes, such as SEAWAT and SUTRA, is that fewer model cells are required for simulations using the SWI2 Package because every aquifer can be represented by a single layer of cells. This reduction in number of required model cells and the elimination of the need to solve the advective-dispersive transport equation results in substantial model run-time savings, which can be large for regional aquifers. The accuracy and use of the SWI2 Package is demonstrated through comparison with existing exact solutions and numerical solutions with SEAWAT. Results for an unconfined aquifer are also presented to demonstrate application of the SWI2 Package to a large-scale regional problem.

Techniques and Methods↗
USGS2000

Documentation of a computer program to simulate lake-aquifer interaction using the MODFLOW ground water flow model and the MOC3D solute-transport model

Heads and flow patterns in surficial aquifers can be strongly influenced by the presence of stationary surface-water bodies (lakes) that are in direct contact, vertically and laterally, with the aquifer. Conversely, lake stages can be significantly affected by the volume of water that seeps through the lakebed that separates the lake from the aquifer. For these reasons, a set of computer subroutines called the Lake Package (LAK3) was developed to represent lake/aquifer interaction in numerical simulations using the U.S. Geological Survey three-dimensional, finite-difference, modular ground-water flow model MODFLOW and the U.S. Geological Survey three-dimensional method-of-characteristics solute-transport model MOC3D. In the Lake Package described in this report, a lake is represented as a volume of space within the model grid which consists of inactive cells extending downward from the upper surface of the grid. Active model grid cells bordering this space, representing the adjacent aquifer, exchange water with the lake at a rate determined by the relative heads and by conductances that are based on grid cell dimensions, hydraulic conductivities of the aquifer material, and user-specified leakance distributions that represent the resistance to flow through the material of the lakebed. Parts of the lake may become "dry" as upper layers of the model are dewatered, with a concomitant reduction in lake surface area, and may subsequently rewet when aquifer heads rise. An empirical approximation has been encoded to simulate the rewetting of a lake that becomes completely dry. The variations of lake stages are determined by independent water budgets computed for each lake in the model grid. This lake budget process makes the package a simulator of the response of lake stage to hydraulic stresses applied to the aquifer. Implementation of a lake water budget requires input of parameters including those representing the rate of lake atmospheric recharge and evaporation, overland runoff, and the rate of any direct withdrawal from, or augmentation of, the lake volume. The lake/aquifer interaction may be simulated in both transient and steady-state flow conditions, and the user may specify that lake stages be computed explicitly, semi-implicitly, or fully-implicitly in transient simulations. The lakes, and all sources of water entering the lakes, may have solute concentrations associated with them for use in solute-transport simulations using MOC3D. The Stream Package of MODFLOW-2000 and MOC3D represents stream connections to lakes, either as inflows or outflows. Because lakes with irregular bathymetry can exist as separate pools of water at lower stages, that coalesce to become a single body of water at higher stages, logic was added to the Lake Package to allow the representation of this process as a user option. If this option is selected, a system of linked pools (sublakes) is identified in each time step and stages are equalized based on current relative sublake surface areas.

Water-Resources Investigations Report↗
USGS2003

Simulation of net infiltration and potential recharge using a distributed-parameter watershed model of the Death Valley region, Nevada and California

This report presents the development and application of the distributed-parameter watershed model, INFILv3, for estimating the temporal and spatial distribution of net infiltration and potential recharge in the Death Valley region, Nevada and California. The estimates of net infiltration quantify the downward drainage of water across the lower boundary of the root zone and are used to indicate potential recharge under variable climate conditions and drainage basin characteristics. Spatial variability in recharge in the Death Valley region likely is high owing to large differences in precipitation, potential evapotranspiration, bedrock permeability, soil thickness, vegetation characteristics, and contributions to recharge along active stream channels. The quantity and spatial distribution of recharge representing the effects of variable climatic conditions and drainage basin characteristics on recharge are needed to reduce uncertainty in modeling ground-water flow. The U.S. Geological Survey, in cooperation with the Department of Energy, developed a regional saturated-zone ground-water flow model of the Death Valley regional ground-water flow system to help evaluate the current hydrogeologic system and the potential effects of natural or human-induced changes. Although previous estimates of recharge have been made for most areas of the Death Valley region, including the area defined by the boundary of the Death Valley regional ground-water flow system, the uncertainty of these estimates is high, and the spatial and temporal variability of the recharge in these basins has not been quantified. To estimate the magnitude and distribution of potential recharge in response to variable climate and spatially varying drainage basin characteristics, the INFILv3 model uses a daily water-balance model of the root zone with a primarily deterministic representation of the processes controlling net infiltration and potential recharge. The daily water balance includes precipitation (as either rain or snow), snow accumulation, sublimation, snowmelt, infiltration into the root zone, evapotranspiration, drainage, water content change throughout the root-zone profile (represented as a 6-layered system), runoff (defined as excess rainfall and snowmelt) and surface water run-on (defined as runoff that is routed downstream), and net infiltration (simulated as drainage from the bottom root-zone layer). Potential evapotranspiration is simulated using an hourly solar radiation model to simulate daily net radiation, and daily evapotranspiration is simulated as an empirical function of root zone water content and potential evapotranspiration. The model uses daily climate records of precipitation and air temperature from a regionally distributed network of 132 climate stations and a spatially distributed representation of drainage basin characteristics defined by topography, geology, soils, and vegetation to simulate daily net infiltration at all locations, including stream channels with intermittent streamflow in response to runoff from rain and snowmelt. The temporal distribution of daily, monthly, and annual net infiltration can be used to evaluate the potential effect of future climatic conditions on potential recharge. The INFILv3 model inputs representing drainage basin characteristics were developed using a geographic information system (GIS) to define a set of spatially distributed input parameters uniquely assigned to each grid cell of the INFILv3 model grid. The model grid, which was defined by a digital elevation model (DEM) of the Death Valley region, consists of 1,252,418 model grid cells with a uniform grid cell dimension of 278.5 meters in the north-south and east-west directions. The elevation values from the DEM were used with monthly regression models developed from the daily climate data to estimate the spatial distribution of daily precipitation and air temperature. The elevation values were also used to simulate atmosp

Water-Resources Investigations Report↗
USGS2003

Conjunctive-use optimization model of the Mississippi River Valley alluvial aquifer of northeastern Arkansas

The Mississippi River Valley alluvial aquifer is a water-bearing assemblage of gravels and sands that underlies about 32,000 square miles of Missouri, Kentucky, Tennessee, Mississippi, Louisiana, and Arkansas. Because of the heavy demands placed on the aquifer, several large cones of depression over 100 feet deep have formed in the potentiometric surface, resulting in lower well yields and degraded water quality in some areas. A ground-water flow model of the alluvial aquifer was previously developed for an area covering 14,104 square miles, extending northeast from the Arkansas River into the northeast corner of Arkansas and parts of southeastern Missouri. The flow model showed that continued ground-water withdrawals at rates commensurate with those of 1997 could not be sustained indefinitely without causing water levels to decline below half the original saturated thickness of the aquifer. To develop estimates of withdrawal rates that could be sustained in compliance with the constraints of critical ground-water area designation, conjunctive-use optimization modeling was applied to the flow model of the alluvial aquifer in northeastern Arkansas. Ground-water withdrawal rates form the basis for estimates of sustainable yield from the alluvial aquifer and from rivers specified within the alluvial aquifer model. A management problem was formulated as one of maximizing the sustainable yield from all ground-water and surface-water withdrawal cells within limits imposed by plausible withdrawal rates, and within specified constraints involving hydraulic head and streamflow. Steady-state flow conditions were selected because the maximized withdrawals are intended to represent sustainable yield of the system (a rate that can be maintained indefinitely). Within the optimization model, 11 rivers are specified. Surface-water diversion rates that occurred in 2000 were subtracted from specified overland flow at the appropriate river cells. Included in these diversions were the planned diversions of 63,339,248 ft3/d for the Bayou Meto project area and 55,078,367 ft3/d for the Grand Prairie project area, which factor in an additional 30 and 40 percent transmission loss, respectively. Streamflow constraints were specified at all 1,165 river cells based on average 7-day minimum flows for 10 years. Sustainable yield for all rivers ranged from 0 (Current, Little Red, and Bayou Meto Rivers) to almost 5 billion cubic feet per day for the Arkansas River. Total sustainable yield from all rivers combined was 12.8 billion cubic feet per day, which represents a substantial source for supplementing ground water to meet the total water demand. Sustainable-yield estimates are affected by the allowable upper limit on withdrawals from wells specified in the optimization model. Ground-water withdrawal rates were allowed to vary as much as 200 percent of the withdrawal rate in 1997. As the overall upper limit on withdrawals is increased, the sustainable yield generally increases. Tests with the optimization model show that without limits on pumping, wells adjacent to sources of water would have optimized withdrawal rates that were orders of magnitude larger than rates corresponding to those of 1997. The sustainable yield from ground water for the entire study area while setting the maximum upper limit as the amount withdrawn in 1997 is 360 million cubic feet per day, which is only about 57 percent of the amount withdrawn in 1997 (635.6 million cubic feet per day). Optimal sustainable yields from within the Bayou Meto irrigation project area and within the Grand Prairie irrigation project area are 18.1 and 9.1 million cubic feet per day, respectively, assuming a maximum allowable withdrawal rate equal to 1997 rates. These values of sustainable yield represent 35 and 30 percent respectively of the amount pumped from these project areas in 1997. Unmet demand (defined as the difference between the optimized withdrawal rate or sustainable yield, and the anticipated demand) was calculated using different demand rates based on multiples of the 1997 withdrawal rate. Assuming that demand is the 1997 withdrawal rate, and that sustainable-yield estimates are those obtained using upper limits of withdrawal rates of 100-, 150-, and 200-percent of 1997 withdrawal rates, then the resulting unmet demand for the entire model area is 275.5, 190.9, and 110 million cubic feet per day, respectively. Whereas, if the demand is specified as 100-, 150-, and 200-percent of the 1997 withdrawal rate, and the sustainable-yield estimates remain the same, then the resulting unmet demand for the entire model area is 275.5, 508.8, and 745.8 million cubic feet per day, respectively. These unmet demands for ground water could be obtained from large sustainable surface-water withdrawals.

Arkansas↗
USGS1973

Geohydrology of the Cross-Florida Barge Canal area, with special reference to the Ocala vicinity

The Cross-Florida Barge Canal route commences at Palatka on the St. Johns River, about 75 miles upstream from the Atlantic Ocean, and extends 110 miles southwestward across Peninsular Florida into deep water in the Gulf of Mexico near Yankeetown. The canal will be equipped with five locks, each 600 feet long and 84 feet wide, and the channel will be a minimum of 12 feet deep and 150 feet wide. From near Ocala northeastward, the canal channel will replace much of the natural channel of the Oklawaha River, and will be excavated into beds of the so-called shallow sand aquifer of Miocene age and younger, which overlies limestone of the Floridan aquifer. Westward from Ocala, most of the canal will be excavated below the potentiometric surface into limestone and dolomite of the Floridan aquifer. Water levels of Rodman, Eureka, and Inglis Pools will be controlled by dams and spillways with the limited exchange of water between the pools and the aquifers. The water levels in the Summit Pools will fluctuate with the natural changes in the ground-water level of the Floridan aquifer, although the stage of the pool will be controlled partly by the stage held in the Eureka Pool. A dynamic inflow-outflow relationship will exist between the Summit Pool and the Floridan aquifer. The Floridan aquifer in the canal area is 1,000 to 1,200 feet thick and consists of limestone and dolomite of middle Eocene Miocene age, including from older to younger, the Lake City, Avon Park, and Ocala limestones plus permeable sandy, dolomitic limestone in the lower part of the Hawthorn Formation. It is possible that most of the flow to the two major springs in the area occurs in the upper 100 feet or so of the aquifer in the Ocala Limestone. The aquifer is underlain by the Oldsmar limestone of early Eocene age and is overlain by sand, clayey sand, clay and shell beds of Miocene through Holocene age, in thickness from a few feet to 300 feet. The permeable beds overlying the Floridan aquifer constitute the shallow aquifer, while the poorly permeable ones act as confining beds where the Floridan aquifer is under artesian conditions. A north-south line drawn separating the head of Silver Springs on the west from the Oklawaha River on the east marks the approximate western limit of a continuous blanket of materials of Miocene-Pliocene(?) age covering the rocks of the Floridan aquifer. East of the line, much of the aquifer is under artesian conditions, particularly in the Oklawaha River valley, although in some areas east of the valley, direct recharge through thick permeable Miocene-Pliocene(?) sands occurs. West of the line, only scattered remnants of a once continuous Miocene-Pliocene(?) cover remain. Lack of the cover is a result of erosion on the crest and flank of the Ocala Uplift, a broad northwest-southeast trending anticlinal upwarp, the axis of which is crossed by the canal route in the Dunnellon area. Over most of this area the Floridan aquifer is unconfined and receives direct recharge through a cover of a few tens of feet of sand and clayey sand of Quaternary age. Tensional stresses during the structural evolution of the Ocala Uplift produced an intersecting system of fractures and normal faults in rocks of the Floridan aquifer. The fractures and faults are important controls for orientation of solution channels and, therefore, for development of ground-water circulation patterns. When the system surface streams, which once drained the Barge Canal area, eroded the poorly permeable Miocene-Pliocene(?) cover from the flanks of the Ocala Uplift, surface runoff was reduced and precipitation began to directly infiltrate the underlying limestones. Now only principal streams remain, such as the Oklawaha and Withlacoochee Rivers and a few short tributaries, while one of the most highly developed subsurface drainage systems in the world has evolved in cavernous limestone of the Floridan aquifer. Two of the larger freshwater springs in the world now discharge from the Floridan aquifer in the canal area. Silver Springs near Ocala discharges an average 531 mgd (million gallons per day) down the 4-mile long Silver River, which flows on poorly permeable beds to the Oklawaha River. Rainbow Springs near Dunnellon discharges on average 468 mgd from numerous orifices in the bed of the 5-mile-long Rainbow River, which flows into the Withlacoochee River. The heads of the springs have migrated to their present positions partly because of a tendency of ground-water levels to decline as permeability in the aquifer is increased due to removal of limestone by solution, and because of mechanical erosion of the limestone in the vicinity of the spring heads. Also, points of principal spring discharges have shifted in the past due to changes in ground-water levels in response to changes in sea level. The subsurface drainage system is continuing to evolve today, as evidenced in part by frequent occurrence of new sinkholes and by the presence of significant amounts of calcium bicarbonate in the spring waters. Rodman Pool, at the east end of the canal, is separated from the Floridan aquifer by poorly permeable materials. The pool's operating water level will be only a few feet above the potentiometric surface at the downstream end, and at or slightly below the potentiometric surface at the upstream end. Little exchange of water between the Rodman Pool and the Floridan aquifer is expected. Eureka Pool, just upstream from Rodman Pool, will also be separated from the Floridan aquifer by poorly permeable beds. However, the stage of the pool will be about 15 feet higher than the natural potentiometric surface at the pool's downstream end, and some seepage into the Floridan aquifer is anticipated through faults and leaky parts of the poorly permeable beds, with a consequent rise in ground-water levels in areas adjacent to the lower end of the pool. Possibilities for particulate contamination of the aquifer will tend to be minimized because of the filtering capacity of the materials through which water must pass to reach the aquifer, although the natural filter will not preclude movement into the aquifer of contaminants which might become dissolved in the pool waters. No significant interchange of water between the pool and the aquifer is expected at the upstream end of the Eureka Pool. Present construction plans indicate an operating stage for Eureka Pool which will range between 38 and 40 feet above mean sea level, although it is possible to dredge the pool deep enough to permit a range in stage of 36 to 40 feet. A backwater effect extending up Silver River from the Eureka Pool is expected to regulate the stage at the head of Silver Springs between 39 and 44 feet above mean sea level if the pool ranges between 36 and 40 feet. If Eureka Pool ranges only between 38 and 40 feet, the range of stage at the head of the springs should be about 41 to 44 feet above mean sea level. From Inglis Lock west, the canal will have direct connection with the Gulf, and canal stage will fluctuate with the Gulf tide. Since the canal stage will be slightly lower than the adjacent ground-water levels along much of the reach, there will be some ground-water inflow to the canal. No significant changes in the existing ground-water regime are expected in the vicinity of Inglis Pool, the first step up in the canal east of the Gulf. Existing ground-water and surface-water levels in the area will not change appreciable, and the natural stage and flow of Rainbow Springs, which will flow by way of Rainbow River into Inglis Pool, should not be affected by canal operations. A possible adverse effect of the canal on the Inglis Pool area could result if sea water is locked up from the Gulf through Inglis Lock. However, the high step of 25 feet at the lock, flushing action of continuous flow from Inglis Pool to the lower reaches of the Withlacoochee River, and use of possible preventive locking procedures should minimize the problem. The potential for adverse effects on the ground-water regime is greatest in the area of the Summit Pool. Through most of the length of the pool, the canal channel will be excavated into limestone of the Floridan aquifer to depths of 12 to 27 feet below the potentiometric surface. Changes that will take place in the ground-water flow system in the Silver Springs drainage area, once the canal is completed, were estimated by flow-net analysis. Variation in aquifer transmissivity was determined by calculating transmissivity in 25 different flow cells surrounding Silver Springs. Transmissivity in the 25 cells averages about 15,600,000 gpd/ft [2,090,000 ft2/day (feet squared)], but transmissivity in the six cells through which the Summit Pool passes ranges from 9,000,000 to 44,000,000 gpd/ft (1,210,000 to 5,900,000 ft2/day). Transmissivity was used to compute static stage of the Summit Pool under given ground-water level conditions. Had the canal existed in May 1968 and had the stage of Eureka Pool been held at 36 feet at the time, the static stage in Summit Pool would have been about 42.1 feet above mean sea level. Thus, a conceptual model of the changes in the potentiometric surface wrought by the finished canal was drawn, and zones of ground-water inflow and outflow were delineated. Most outflow from the Summit Pool to the aquifer should be limited to one 4-mile-long zone along the north side of the pool, about 5 miles south of Silver Springs. It is estimated that a water volume equivalent to about 8 percent of the daily flow of Silver Springs will enter the Summit Pool each day from the southern one-third of the Silver Springs drainage area. A like amount will reenter the aquifer at the main zone of outflow and move toward Silver Springs at an estimated average velocity of about 200 feet per day, if something close to the natural static stage of the pool is maintained by return pumpage of the lockage losses. At a velocity of 200 feet per day, water from the Summit Pool would discharge at Silver Springs about 140 days later. However, any estimate of velocity in the highly cavernous limestone aquifer in the area should be used with caution, because difficult to measure changes in porosity and thickness of major zones of flow may cause large variations in velocity. If all lockage losses are returned to the Summit Pool by pumping from Eureka Pool, no net loss from the Silver Springs drainage area, except for some evaporation from the water surface in the canal and possible leakage around locks, will result from canal operations. The zone of outflow from the Summit Pool to the aquifer will be in a natural potentiometric trough, and the zone of inflow will be in a potentiometric ridge area. The equilibrium water level in the pool will tend to be about 1 foot higher than the altitude of the lowest level in the pre-canal potentiometric trough, and about 2 feet lower than the highest level on the pre-canal potentiometric ridge. West of the Silver Springs drainage area just east of Dunnellon Lock, in the area of a local potentiometric high, the water level in the Summit Pool is expected to be about 15 feet below the natural potentiometric surface. In most areas 2 to 3 miles away from the Summit Pool, effects of the canal on the natural potentiometric surface should be slight. The stage of the Summit Pool, judging from the 36-year record for ground-water level changes and the anticipated indirect effect of the controlled stage in Eureka Pool, should have a maximum of about 10.5 feet with a maximum water level of about 51.5 feet above mean sea level and a minimum of about 41.0 feet above mean sea level. Of particular importance in the Summit Pool is an implementation of well planned construction and operational procedures designed to minimize risks of ground-water contamination.

Florida↗
USGS1989

Hydrologic environments and water-quality characteristics at four landfills in Mecklenburg County, North Carolina, 1980-86

A water-quality study was conducted during 1980-86 at four landfills in Mecklenburg County, North Carolina. Each landfill has a three-layered hydrogeologic system typical of the Piedmont, consisting of (1) the regolith; (2) a transition zone; and (3) unweathered, fractured crystalline bedrock. As much as 7.6 inches per year of rainfall enters the ground-water system and has the potential to generate leachate within landfill cells. Ground water and leachate discharge to tributaries within the landfill sites or to streams adjacent to them. Water-quality samples were collected from 53 monitoring wells and 20 surface-water sites. Samples were analyzed for selected physical and biological characteristics, major inorganic ions, nutrients, trace elements, and organic compounds. Selected indicators of water quality, including specific conductance; hardness; and concentrations of chloride, manganese, dissolved solids, total organic carbon, and specific organic compounds were analyzed to determine the effects of each landfill on ground- and surface-water quality. Increases in concentrations of inorganic constituents above background levels were detected in ground water downgradient of the landfills. The increases were generally greatest in samples from wells in close proximity to the older landfill cells. In general, the increases in concentrations in downgradient wells were greater for calcium, magnesium, and chloride than for other major ions. Manganese exhibited the largest relative increase in concentration between upgradient and downgradient wells of any constituent, and manganese concentration data were effective in defining areas with extensive anaerobic biological activity. Differences between upgradient and downgradient concentrations of total organic carbon and specific organic compounds generally were not as apparent. The most frequently identified organic contaminants were the herbicides 2,4-D and 2,4,5-T. Chlorofluoromethanes were identified in three of four ground-water samples analyzed for volatile organic compounds. Landfills affected the water quality of several smaller streams but did not noticeably affect larger ones. Apparent effects on water quality were greatest at the oldest landfill, located on Statesville Road, where waste is in cells that are partly below the water table.

North Carolina↗
USGS1995

Radon-222 concentrations in ground water and soil gas on Indian reservations in Wisconsin

The weighted average radon-222 concentration of indoor air in homes located on Wisconsin Indian Reservations is 5.8 picocuries per liter, which exceeds the U.S. Environmental Protection Agency action limit of 4 picocuries per liter. Ground water is the principle source of drinking water on Wisconsin Indian Reservations and generally accounts for about 5 percent of the total indoor air radon-222 concentrations found in homes. To determine the distribution of radon-222, ground water from 29 private and community Wisconsin Indian Reservation wells and soil gas at a depth of about 3 feet below land surface adjacent to the wells were sampled. Sites with wells were distributed among the 11 Wisconsin Indian Reservations so that each Reservation contained at least 2 sites. The remaining seven sites were divided among the Reservation by acreage held by each tribe. Ground-water samples were collected using a syringe technique after the wells had been purged sufficiently to reach chemical stability. Samples were then sent by overnight mail to the U.S. Geological Survey, National Water Quality Lab for analysis by liquid scintillation counting. Soil gas was collected in Lucas cells, using a small-diameter soil probe and peristaltic pump. The cells were then analyzed by Lucas cell alpha scintillation counting with a portable radon detector. The highest radon-222 concentrations in ground water and soil gas are from sites with wells finished in the crystalline bedrock aquifer. Radon-222 concentrations in 29 ground-water samples collected from the sand and gravel and sedimentary and crystalline bedrock aquifers range from 260 to 22,000 picocuries per liter with a median concentration of 560 picocuries per liter. Only 2 of the 29 ground-water samples have radon-222 concentrations less than the U.S. Environmental Protection Agency proposed standard of 300 picocuries per liter. The highest radon-222 concentrations were found in ground water from wells in Shawano, Menominee, Forest, and Marathon counties. Radon-222 concentrations in soil gas range from 130 to 7,810 picocuries per liter with a median concentration of 560 picocuries per liter. For sites with wells finished in the sand and gravel aquifer, the coefficient of determination (R2) of the regression of concentration of radon-222 in ground water as a function of well depth is 0.003 and the significance level is 0.32, which indicates that there is not a statistically significant relation between radon-222 concentrations in ground water and well depth. The coefficient of determination of the regression of radon-222 in ground water and soil gas is 0.19 and the root mean square error of the regression line is 271 picocuries per liter. Even though the significance level (0.036) indicates a statistical relation, the root mean square error of the regression is so large that the regression equation would not give reliable predictions. Because of an inadequate number of samples, similar statistical analyses could not be performed for sites with wells finished in the crystalline and sedimentary bedrock aquifers.

Wisconsin↗
USGS1996

Description and field analysis of a coupled ground-water/surface-water flow model (MODFLOW/BRANCH) with modifications for structures and wetlands in southern Dade County, Florida

A coupled surface-water model (BRANCH) and ground-water model (MODFLOW) model were tested to simulate the interacting wetlands/surface-water/ ground-water system of southern Dade County. Several options created for the MODFLOW ground- ground-water model were used in representing this field situation. The primary option is the MODBRANCH interfacing software, which allows leakage to be accounted for between the MODFLOW ground-water model and the BRANCH dynamic model for simulation of flow in an interconnected network of open channels. A modification to an existing software routine, which is referred to as BCF2, allows cells in MODFLOW to rewet when dry--a requirement in representing the seasonal wetlands in Dade County. A companion to BCF2 is the modified evapotranspiration routine EVT2. The EVT2 routine changes the cells where evapotranspiration occurs, depending on which cells are wet. The Streamlink package represents direct connections between the canals and wetlands at locations where canals open directly into overland flow. Within the BRANCH model, the capability to represent the numerous hydraulic structures, gated spillways, gated culverts, and pumps was added. The application of these modifications to model surface-water/ground-water interactions in southern Dade County demonstrated the usefulness of the coupled MODFLOW/BRANCH model. Ground-water and surface-water flows are both simulated with dynamic models. Flow exchange between models, intermittent wetting and drying, evapotranspiration, and hydraulic structure operations are all represented appropriately. Comparison was made with a simulation using the RIV1 package instead of MODBRANCH to represent the canals. RIV1 represents the canals by user-defined stages, and computes leakage to the aquifer. Greater accuracy in reproducing measured ground- water heads was achieved with MODBRANCH, which also computes dynamic flow conditions in the canals, unlike RIV1. The surface-water integrated flow and transport two-dimensional model (SWIFT2D) was also applied to the southeastern coastal wetlands for comparison with the wetlands flow approximation made in MODFLOW. MODFLOW simulates the wetlands as a highly conductive upper layer of the aquifer, whereas SWIFT2D solves the hydrodynamic equations. Comparison in this limited test demonstrated no specific advantage for either method of representation. However, much additional testing on a wider variety of geometric and hydraulic situations, such as in areas with greater tidal or other dynamic forcing effects, is needed to make definite conclusions. A submodel of the existing southern Dade County model schematization was used to examine water-delivery alternatives proposed by the U.S. Army Corps of Engineers. For this application, the coupled MODFLOW/BRANCH model was used as a design tool. A new canal and several pumps to be tested to maintain lower water levels in a residential area (while water levels in the Everglades are raised) were added to the model schematization. The pumps were assumed to have infinite supply capacity in the model so that their maximum pumping rates during the simulation could be used to determine pump sizes.

Water-Resources Investigations Report↗
USGS1993

Numerical simulation of ground-water flow in the central part of the western San Joaquin Valley, California

The occurrence of selenium in agricultural drain water in the central part of the western San Joaquin Valley, California, has focused concern on strategies for managing shallow, saline ground water. To assess alternatives to agricultural drains, a three-dimensional, finite-difference numerical model of the regional groundwater flow system was developed. This report documents the mathematical approach used to model the flow system, the data base on which the model is based, and the methods used to calibrate the model. The 550-square-mile study area includes parts of the Panoche Creek alluvial fan and parts of the Little Panoche Creek and Cantua Creek alluvial fans. The model simulates transient flow in the semiconfined and confined zones above and below the Corcoran Clay Member of the Tulare Formation of Pleistocene age. The model incorporates areally distributed ground-water recharge, areally and vertically distributed pumping, regional-collector drains in the Wesdands Water District (operative from 1980 to 1985), on-farm drains in parts of the Panoche, Broadview, and Firebaugh Water Districts, and bare-soil evaporation (which occurs if the water table is within 7 feet of land surface). The model also incorporates texture-based estimates of hydraulic conductivity, where texture is defined as the fraction of coarse-grained deposits present in a given subsurface interval. The numerical model was developed using hydrologic data from 1972 to 1988. Most of the parameters incorporated into the model were evaluated independently of the model, including system geometry, the distribution of texture, the altitudes of the water table and potentiometric surface of the confined zone in 1972 (initial condition), the hydraulic conductivity of coarse-grained deposits derived from the Coast Ranges, the hydraulic conductivity of coarse-grained deposits derived from the Sierra Nevada, specific storage, recharge, pumping, and parameters needed to incorporate drains and bare-soil evaporation. Four parameters were calibration variables: the hydraulic conductivity of fine-grained deposits in the semiconfined zone, the hydraulic conductivity of the Corcoran Clay Member, specific yield, and the transmissivity of the confined zone. The model was calibrated in two phases. In the first phase, a steady-state model of the ground-water flow system in 1984 was used to constrain the relation between the hydraulic conductivity of fine-grained deposits in the semiconfined zone and the hydraulic conductivity of the Corcoran Clay Member, thus reducing the number of independent variables from four to three. In the second phase of calibration, the change in altitude of the water table from 1972 to 1984, the change in altitude of the potentiometric surface of the confined zone from 1972 to 1984, and the number of model cells subject to bare-soil evaporation from 1972 to 1988 were used to evaluate the remaining three variables. The calibrated model reproduces the average change in water-table altitude (1972-84) to within 0.4 foot (average measured change 11.5 feet) and the average change in confined zone head (1972- 84) to within 19 feet (average measured change 120 feet). The simulated time-series record of the total number of model cells subject to bare-soil evaporation (each cell is 1 mile square) is within the range of the measured data. The measured values are at a minimum in October and a maximum in July. The October values ranged from 103 in 1972 to 132 in 1984 (the drains were closed in 1985) to 151 in 1988. The July values ranged from 144 in 1973 to 198 in 1984, to 204 in 1988. The simulated values ranged from 103 in 1972 to 161 in 1984, to 208 in 1988.

Water Supply Paper↗