USGS Science⌕ Search

SEARCH · USGS Science

Results for “Ground Water Modeling Investigation Report”

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 937 records · Page 52Linked to original sources

Vertical gradients in water chemistry and age in the Northern High Plains Aquifer, Nebraska, 2003

The northern High Plains aquifer is the primary source of water used for domestic, industrial, and irrigation purposes in parts of Colorado, Kansas, Nebraska, South Dakota, and Wyoming. Despite the aquifer’s importance to the regional economy, fundamental ground-water characteristics, such as vertical gradients in water chemistry and age, remain poorly defined. As part of the U.S. Geological Survey’s National Water-Quality Assessment Program, water samples from nested, short-screen monitoring wells installed in the northern High Plains aquifer were analyzed for major ions, nutrients, trace elements, dissolved organic carbon, pesticides, stable and radioactive isotopes, dissolved gases, and other parameters to evaluate vertical gradients in water chemistry and age in the aquifer. Chemical data and tritium and radiocarbon ages show that water in the aquifer was chemically and temporally stratified in the study area, with a relatively thin zone of recently recharged water (less than 50 years) near the water table overlying a thicker zone of older water (1,800 to 15,600 radiocarbon years). In areas where irrigated agriculture was an important land use, the recently recharged ground water was characterized by elevated concentrations of major ions and nitrate and the detection of pesticide compounds. Below the zone of agricultural influence, major-ion concentrations exhibited small increases with depth and distance along flow paths because of rock/water interactions. The concentration increases were accounted for primarily by dissolved calcium, sodium, bicarbonate, sulfate, and silica. In general, the chemistry of ground water throughout the aquifer was of high quality. None of the approximately 90 chemical constituents analyzed in each sample exceeded primary drinking-water standards. Mass-balance models indicate that changes in groundwater chemistry along flow paths in the aquifer can be accounted for by small amounts of feldspar and calcite dissolution; goethite and clay-mineral precipitation; organic-carbon and pyrite oxidation; oxygen reduction and denitrification; and cation exchange. Mixing with surface water affected the chemistry of ground water in alluvial sediments of the Platte River Valley. Radiocarbon ages in the aquifer, adjusted for carbon mass transfers, ranged from 1,800 to 15,600 14C years before present. These results have important implications with respect to development of ground-water resources in the Sand Hills. Most of the water in the aquifer predates modern anthropogenic activity so excessive removal of water by pumping is not likely to be replenished by natural recharge in a meaningful timeframe. Vertical gradients in ground-water age were used to estimate long-term average recharge rates in the aquifer. In most areas, the recharge rates ranged from 0.02 to 0.05 foot per year. The recharge rate was 0.2 foot per year in one part of the aquifer characterized by large downward hydraulic gradients. Nitrite plus nitrate concentrations at the water table were 0.13 to 3.13 milligrams per liter as nitrogen, and concentrations substantially decreased with depth in the aquifer. Dissolved-gas and nitrogen-isotope data indicate that denitrification in the aquifer removed 0 to 97 percent (average = 50 percent) of the nitrate originally present in recharge. The average amount of nitrate removed by denitrification in the aquifer north of the Platte River (Sand Hills) was substantially greater than the amount removed south of the river (66 as opposed to 0 percent), and the extent of nitrate removal appears to be related to the presence of thick deposits of sediment on top of the Ogallala Group in the Sand Hills that contained electron donors, such as organic carbon and pyrite, to support denitrification. Apparent rates of dissolved-oxygen reduction and denitrification were estimated on the basis of decreases in dissolved-oxygen concentrations and increases in concentrations of excess nitrogen gas and ground-water ages along flow paths from the water table to deeper wells. Median rates of dissolved-oxygen reduction and denitrification south of the Platte River were at least 10 times smaller than the median rates north of the river in the Sand Hills. The relatively large denitrification rates in the Sand Hills indicate that the aquifer in that area may have a greater capacity to attenuate nitrate contamination than the aquifer south of the river, depending on rates of ground-water movement in the two areas. Small denitrification rates south of the river indicate that nitrate contamination in that part of the aquifer would likely persist for a longer period of time.

Nebraska↗

Geophysical Methods for Investigating Ground-Water Recharge

While numerical modeling has revolutionized our understanding of basin-scale hydrologic processes, such models rely almost exclusively on traditional measurements?rainfall, streamflow, and water-table elevations?for calibration and testing. Model calibration provides initial estimates of ground-water recharge. Calibrated models are important yet crude tools for addressing questions about the spatial and temporal distribution of recharge. An inverse approach to recharge estimation is taken of necessity, due to inherent difficulties in making direct measurements of flow across the water table. Difficulties arise because recharging fluxes are typically small, even in humid regions, and because the location of the water table changes with time. Deep water tables in arid and semiarid regions make recharge monitoring especially difficult. Nevertheless, recharge monitoring must advance in order to improve assessments of ground-water recharge. Improved characterization of basin-scale recharge is critical for informed water-resources management. Difficulties in directly measuring recharge have prompted many efforts to develop indirect methods. The mass-balance approach of estimating recharge as the residual of generally much larger terms has persisted despite the use of increasing complex and finely gridded large-scale hydrologic models. Geophysical data pertaining to recharge rates, timing, and patterns have the potential to substantially improve modeling efforts by providing information on boundary conditions, by constraining model inputs, by testing simplifying assumptions, and by identifying the spatial and temporal resolutions needed to predict recharge to a specified tolerance in space and in time. Moreover, under certain conditions, geophysical measurements can yield direct estimates of recharge rates or changes in water storage, largely eliminating the need for indirect measures of recharge. This appendix presents an overview of physically based, geophysical methods that are currently available or under development for recharge monitoring. The material is written primarily for hydrogeologists. Uses of geophysical methods for improving recharge monitoring are explored through brief discussions and case studies. The intent is to indicate how geophysical methods can be used effectively in studying recharge processes and quantifying recharge. As such, the material constructs a framework for matching the strengths of individual geophysical methods with the manners in which they can be applied for hydrologic analyses. The appendix is organized in three sections. First, the key hydrologic parameters necessary to determine the rate, timing, and patterns of recharge are identified. Second, the basic operating principals of the relevant geophysical methods are discussed. Methods are grouped by the physical property that they measure directly. Each measured property is related to one or more of the key hydrologic properties for recharge monitoring. Third, the emerging conceptual framework for applying geophysics to recharge monitoring is presented. Examples of the application of selected geophysical methods to recharge monitoring are presented in nine case studies. These studies illustrate hydrogeophysical applications under a wide range of conditions and measurement scales, which vary from tenths of a meter to hundreds of meters. The case studies include practice-proven as well as emerging applications of geophysical methods to recharge monitoring.

Professional Paper↗

Chemical Characteristics, Water Sources and Pathways, and Age Distribution of Ground Water in the Contributing Recharge Area of a Public-Supply Well near Tampa, Florida, 2002-05

In 2001, the National Water-Quality Assessment (NAWQA) Program of the U.S. Geological Survey began a series of studies on the transport of anthropogenic and natural contaminants (TANC) to public-supply wells. The main goal of the TANC program was to better understand the source, transport, and receptor factors that control contaminant movement to public-supply wells in representative aquifers of the United States. Studies were first conducted at regional scales at four of the eight TANC study areas during 2002-03 and at small (local) scales during 2003-05 in California, Nebraska, Connecticut, and Florida. In the Temple Terrace study area near Tampa, Florida, multiple chemical indicators and geochemical and ground-water flow modeling techniques were used to assess the vulnerability of a public-supply well in the karstic Upper Floridan aquifer to contamination from anthropogenic and naturally occurring contaminants. During 2003-05, water samples were collected from the public-supply well and 13 surrounding monitoring wells that all tap the Upper Floridan aquifer, and from 15 monitoring wells in the overlying surficial aquifer system and the intermediate confining unit that are located within the modeled ground-water contributing recharge area of the public-supply well. Six volatile organic compounds and four pesticides were detected in trace concentrations (well below drinking-water standards) in water from the public-supply well, which had an open interval from 36 to 53 meters below land surface. These contaminants were detected more frequently in water samples from monitoring wells in the overlying clastic surficial aquifer system than in water from monitoring wells in the Upper Floridan aquifer in the study area. Likewise, nitrate-N concentrations in the public-supply well (0.72-1.4 milligrams per liter) were more similar to median concentrations in the oxic surficial aquifer system (2.1 milligrams per liter) than to median nitrate-N concentrations in the anoxic Upper Floridan aquifer (0.06 milligram per liter) under sulfate-reducing conditions. High concentrations of radon-222 and uranium in the public-supply well compared to those in monitoring wells in the Upper Floridan aquifer appear to originate from water moving downward through sands and discontinuous clay lenses that overlie the aquifer. Water samples also were collected from three overlapping depth intervals (38-53, 43-53, and 49-53 meters below land surface) in the public-supply well. The 49- to 53-meter interval was identified as a high-flow zone during geophysical logging of the wellbore. Water samples were collected from these depth intervals at a low pumping rate by placing a low-capacity submersible pump (less than 0.02 cubic meter per minute) at the top of each interval. To represent higher pumping conditions, a large-capacity portable submersible pump (1.6 cubic meters per minute) was placed near the top of the open interval; water-chemistry samples were collected using the low-capacity submersible pump. The 49- to 53-meter depth interval had distinctly different chemistry than the other two sampled intervals. Higher concentrations of nitrate-N, atrazine, radon, trichloromethane (chloroform), and arsenic (and high arsenic (V)/arsenic (III) ratios); lower concentrations of dissolved solids, strontium, iron, manganese, and lower nitrogen and sulfur isotope ratios were found in this highly transmissive zone in the limestone than in water from the two other depth intervals. Movement of water likely occurs from the overlying sands and clays of the oxic surficial aquifer system and intermediate confining unit (that contains high radon-222 and nitrate-N concentrations) into the anoxic Upper Floridan aquifer (that contains low radon-222 and nitrate-N concentrations). Differences in arsenic concentrations in water from the various depth intervals in the public-supply well (3.2-19.0 micrograms per liter) were related to pumping conditions. The high arsenic

Scientific Investigations Report↗

Hydrogeology and water quality of glacial-drift aquifers in the Bemidji-Bagley area, Beltrami, Clearwater, Cass, and Hubbard Counties, Minnesota

Unconfined and the upper confined aquifers in glacial drift are the primary sources of water in a 1,600 square-mile area including parts of Beltrami, Cass, Clearwater, and Hubbard Counties, Minnesota. The unconfineddrift aquifer consists of coarse sand and gravel in the center of the study area. The total area underlain by the unconfined-drift aquifer is approximately 550 square miles. The unconfined aquifer ranges in thickness from 0 to 130 feet, and is greater than 20 feet thick over an area of 280 square miles. On the basis of scant data, the transmissivity of the unconfined aquifer ranges from less than 70 feet squared per day in the south and west to greater than 8,900 feet squared per day in an area west of Bemidji. Well yields from 10 to 300 gallons per minute are possible in some areas. The unconfined and upper confined-drift aquifers are separated by a fine-grained confining unit of till or lake deposits. The thickness of the upper confined-drift aquifer ranges from 0 to 60 feet in the Bemidji area. On the basis of specific-capacity and aquifer-thickness data, and results of model simulations, the transmissivity of the upper confined-drift aquifer ranges from less than 100 feet squared per day in the south and west parts of the aquifer to about 12,800 feet squared per day in the area around Bemidji. Well yields of 10 to 2,100 gallons per minute are possible in some areas. The direction of ground-water flow in both unconfined and upper confineddrift aquifers is toward the Mississippi and Clearwater Rivers. These rivers are the major discharge points for both aquifers. Ground-water divides, which separate the ground-water flow systems that discharge to these rivers, are approximately coincidental with surface-water divides between the rivers. Water from both aquifers generally is of the calcium bicarbonate type and is very hard, averaging 309 and 267 milligrams per liter as CaCO 3 from confined and unconfined-drift aquifers, respectively. Water from both aquifers generally is suitable for drinking, crop irrigation, and most other uses. Concentrations of ammonia, boron, chromium, iron, manganese, and phenols, however, locally exceed recommended limits for drinking water (Minnesota Pollution Control Agency, 1988). Longer residence time and leakage through glacial till is believed to cause higher concentrations of common inorganic constituents in water from confined-drift aquifers than concentrations in water from the unconfined-drift aquifer. Elevated concentrations of nutrients, chloride, and phenols in the unconfined-drift aquifer may be related to land-use practices. Statistical comparisons of common chemical constituents in water from wells completed in the unconfined-drift aquifer in several land-use areas suggest that concentrations of many constituents and physical properties are generally greater for wells in areas of commercial and residential land-use than for wells in areas of agriculture or forest land-use. These constituents include ammonia plus organic nitrogen, phosphorus, calcium, sodium, potassium, chloride, sulfate, silica, dissolved solids, and specific conductances. The mean values of ammonia nitrogen, magnesium, and fluoride are generally greater for wells in commercial land-use type areas than for wells in forested and agricultural land-use type areas. The mean concentration of nitrogen (N0 2 + NO 3 , dissolved) is generally greater for wells in residential land-use type areas than for wells in forested and agricultural land-use type areas. The Kruskil-Wallis test, a nonparametric that for 12 of the 21 constituents sampled in groups in the unconfined-drift aquifer, a of these constituents and land use was found statistical technique, indicated common in all land-use type relation between the concentration to be statistically significant.

Minnesota↗

Water Quality, Hydrology, and Response to Changes in Phosphorus Loading of Nagawicka Lake, a Calcareous Lake in Waukesha County, Wisconsin

Nagawicka Lake is a 986-acre, usually mesotrophic, calcareous lake in southeastern Wisconsin. Because of concern over potential water-quality degradation of the lake associated with further development in its watershed, a study was conducted by the U.S. Geological Survey from 2002 to 2006 to describe the water quality and hydrology of the lake; quantify sources of phosphorus, including those associated with urban development; and determine the effects of past and future changes in phosphorus loading on the water quality of the lake. All major water and phosphorus sources were measured directly, and minor sources were estimated to construct detailed water and phosphorus budgets for the lake. The Bark River, near-lake surface inflow, precipitation, and ground water contributed 74, 8, 12, and 6 percent of the inflow, respectively. Water leaves the lake primarily through the Bark River outlet (88 percent) or by evaporation (11 percent). The water quality of Nagawicka Lake has improved dramatically since 1980 as a result of decreasing the historical loading of phosphorus to the lake. Total input of phosphorus to the lake was about 3,000 pounds in monitoring year (MY) 2003 and 6,700 pounds in MY 2004. The largest source of phosphorus entering the lake was the Bark River, which delivered about 56 percent of the total phosphorus input, compared with about 74 percent of the total water input. The next largest contributions were from the urbanized near-lake drainage area, which disproportionately accounted for 37 percent of the total phosphorus input but only about 5 percent of the total water input. Simulations with water-quality models within the Wisconsin Lakes Modeling Suite (WiLMS) indicated the response of Nagawicka Lake to 10 phosphorus-loading scenarios. These scenarios included historical (1970s) and current (base) years (MY 2003-04) for which lake water quality and loading were known, six scenarios with percentage increases or decreases in phosphorus loading from controllable sources relative to the base years 2003-04, and two scenarios corresponding to specific management actions. Because of the lake's calcareous character, the average simulated summer concentration of total phosphorus for Nagawicka Lake was about 2 times that measured in the lake. The models likely over-predict because they do not account for coprecipitation of phosphorus and dissolved organic matter with calcite, negligible release of phosphorus from the deep sediments, and external phosphorus loading with abnormally high amounts of nonavailable phosphorus. After adjusting the simulated results for the overestimation of the models, a 50-percent reduction in phosphorus loading resulted in an average predicted phosphorus concentration of 0.008 milligrams per liter (mg/L) (a decrease of 46 percent). With a 50-percent increase in phosphorus loading, the average predicted concentration was 0.020 mg/L (an increase of 45 percent). With the changes in land use under the assumed future full development conditions, the average summer total phosphorus concentration should remain similar to that measured in MY 2003-04 (approximately 0.014 mg/L). However, if stormwater and nonpoint controls are added to achieve a 50-percent reduction in loading from the urbanized near-lake drainage area, the average summer total phosphorus concentration should decrease from the present conditions (MY 2003-04) to 0.011 mg/L. Slightly more than a 25-percent reduction in phosphorus loading from that measured in MY 2003-04 would be required for the lake to be classified as oligotrophic.

Scientific Investigations Report↗

Assessment of the spatial extent and height of flooding in Lake Champlain during May 2011, using satellite remote sensing and ground-based information

Landsat 5 and moderate resolution imaging spectro-radiometer satellite imagery were used to map the area of inundation of Lake Champlain, which forms part of the border between New York and Vermont, during May 2011. During this month, the lake’s water levels were record high values not observed in the previous 150 years. Lake inundation area determined from the satellite imagery is correlated with lake stage measured at three U.S. Geological Survey lake level gages to provide estimates of lake area at different lake levels (stage/area rating) and also compared with the levels of the high-water marks (HWMs) located on the Vermont side of the lake. The rating developed from the imagery shows a somewhat different relation than a similar stage/area rating developed from a medium-resolution digital elevation model (DEM) of the region. According to the rating derived from the imagery, the lake surface area during the peak lake level increased by about 17 percent above the average or “normal” lake level. By using a comparable rating developed from the DEM, the increase above average is estimated to be about 12 percent. The northern part of the lake (north of Burlington) showed the largest amount of flooding. Based on intersecting the inundation maps with the medium-resolution DEM, lake levels were not uniform around the lake. This is also evident from the lake level gage measurements and HWMs. The gage data indicate differences up to 0.5 feet between the northern and southern end of the lake. Additionally, the gage data show day-to-day and intradaily variation of the same range (0.5 foot). The high-water mark observations show differences up to 2 feet around the lake, with the highest level generally along the south- and west-facing shorelines. The data suggest that during most of May 2011, water levels were slightly higher and less variable in the northern part of the lake. These phenomena may be caused by wind effects as well as proximity to major river inputs to the lake. The inundation areas generated from the imagery generally coincide with flood mapping as estimated by the Federal Emergency Management Agency (FEMA) and shown on its digital flood insurance rate maps. Where areas in the flood inundation map derived from the imagery and the FEMA estimated flooded areas differ substantially, this difference may be due to differences between the flood magnitude at the time of the image and the assumed flood condition used for the FEMA modeling and mapping, wind/storage effects not accounted for by the FEMA modeling, and the resolution of the image compared to the DEM used in the FEMA mapping.

Lake Champlain↗

A simulation-optimization model for water-resources management, Santa Barbara, California

In times of drought, the local water supplies of the city of Santa Barbara, California, are insufficient to satisfy water demand. In response, the city has built a seawater desalination plant and gained access to imported water in 1997. Of primary concern to the city is delivering water from the various sources at a minimum cost while satisfying water demand and controlling seawater intrusion that might result from the overpumping of ground water. A simulation-optimization model has been developed for the optimal management of Santa Barbara?s water resources. The objective is to minimize the cost of water supply while satisfying various physical and institutional constraints such as meeting water demand, maintaining minimum hydraulic heads at selected sites, and not exceeding water-delivery or pumping capacities. The model is formulated as a linear programming problem with monthly management periods and a total planning horizon of 5 years. The decision variables are water deliveries from surface water (Gibraltar Reservoir, Cachuma Reservoir, Cachuma Reservoir cumulative annual carryover, Mission Tunnel, State Water Project, and desalinated seawater) and ground water (13 production wells). The state variables are hydraulic heads. Basic assumptions for all simulations are that (1) the cost of water varies with source but is fixed over time, and (2) only existing or planned city wells are considered; that is, the construction of new wells is not allowed. The drought of 1947?51 is Santa Barbara?s worst drought on record, and simulated surface-water supplies for this period were used as a basis for testing optimal management of current water resources under drought conditions. Assumptions that were made for this base case include a head constraint equal to sea level at the coastal nodes; Cachuma Reservoir carryover of 3,000 acre-feet per year, with a maximum carryover of 8,277 acre-feet; a maximum annual demand of 15,000 acre-feet; and average monthly capacities for the Cachuma and the Gibraltar Reservoirs. The base-case results indicate that water demands can be met, with little water required from the most expensive water source (desalinated seawater), at a total cost of $5.56 million over the 5-year planning horizon. The simulation model has drains, which operate as nonlinear functions of heads and could affect the model solutions. However, numerical tests show that the drains have little effect on the optimal solution. Sensitivity analyses on the base case yield the following results: If allowable Cachuma Reservoir carryover is decreased by about 50 percent, then costs increase by about 14 percent; if the peak demand is decreased by 7 percent, then costs will decrease by about 14 percent; if the head constraints are loosened to -30 feet, then the costs decrease by about 18 percent; if the heads are constrained such that a zero hydraulic gradient condition occurs at the ocean boundary, then the optimization problem does not have a solution; if the capacity of the desalination plant is constrained to zero acre-feet, then the cost increases by about 2 percent; and if the carryover of State Water Project water is implemented, then the cost decreases by about 0.5 percent. Four additional monthly diversion distribution scenarios for the reservoirs were tested: average monthly Cachuma Reservoir deliveries with the actual (scenario 1) and proposed (scenario 2) monthly distributions of Gibraltar Reservoir water, and variable monthly Cachuma Reservoir deliveries with the actual (scenario 3) and proposed (scenario 4) monthly distributions of Gibraltar Reservoir water. Scenario 1 resulted in a total cost of about $7.55 million, scenario 2 resulted in a total cost of about $5.07 million, and scenarios 3 and 4 resulted in a total cost of about $4.53 million. Sensitivities of the scenarios 1 and 2 to desalination-plant capacity and State Water Project water carryover were tested. The scenario 1 sensitivity analysis indicated that incorpo

Water-Resources Investigations Report↗

Hydrogeologic characterization of the Hualapai Plateau on the western Hualapai Indian Reservation, northwestern Arizona

This study was developed to assess if groundwater from the western Hualapai Plateau could be used to supply developments in the Grand Canyon West area of the Hualapai Indian Reservation and to collect hydrogeologic data for future use in a numerical groundwater model for the reservation. Ground-based geophysical surveys; existing well, spring, and other hydrogeologic information from previous studies; and new well and spring data collected for this study were used to provide a better understanding of the hydrogeology of the western Hualapai Plateau. Surface geophysical data provided information on the depth and geologic structure of lower Paleozoic rock units and Proterozoic crystalline and metamorphic rocks that underlie the western Hualapai Plateau. The surface geophysical data and discharge information from springs were used to select a site to drill and develop the U.S. Geological Survey Hualapai Test Well. The Hualapai Test Well was drilled to understand the geophysical properties of geologic formations at depth. These data were used to verify the results of surface geophysical data and to evaluate if sufficient water was present in the Hualapai Test Well for potential groundwater development. The Hualapai Test Well was drilled to a depth of 2,468 feet and bottomed in Proterozoic granite. Water was expected in the lower part of the Muav Limestone, but water was not observed until the Tapeats Sandstone at a depth of 2,400 feet. The Tapeats Sandstone was determined to be confined with a hydrostatic head of over 900 feet. A 48-hour pumping test was conducted to determine aquifer properties. Low specific capacity indicated that although groundwater is present in the Tapeats Sandstone, well yields are likely to be small. A water-quality sample indicated the sample had a calcium, magnesium-bicarbonate water type with a total dissolved-solids concentration of 371 milligrams per liter. Alpha radioactivity of the sample, 18.3 picocuries per liter, exceeded the U.S. Environmental Protection Agency maximum contaminant level of 15 picocuries per liter for drinking water. Concentrations of iron and manganese in the water sample also exceeded the U.S. Environmental Protection Agency secondary maximum contaminant levels for drinking water. An inventory of wells and springs provided insight into the occurrence of groundwater on the western Hualapai Plateau. Data from 56 springs on and adjacent to the western Hualapai Plateau were compiled for this study, and new data were collected at 31 springs. Discharge from springs visited for this study ranged from dry to about 345 gallons per minute. The temporal data from springs, where repeat measurements were available, indicated that spring flow is highly variable and likely related to seasonal and annual precipitation. Water levels from 36 wells on and adjacent to the western Hualapai Plateau were compiled for this study, and new water levels were collected at 5 wells. The spring and well data in conjunction with the Hualapai Test Well results indicated that on the western Hualapai Plateau, bedrock aquifers have limited discrete flow paths that make extensive groundwater development unlikely.

Arizona↗

Hydrogeology, digital solute-transport simulation, and geochemistry of the Lower Cretaceous aquifer system near Baltimore, Maryland

This study was made to develop information on the hydrogeology and ground-water geochemistry of the Patuxent and Patapsco aquifers (Lower Cretaceous) near Baltimore, Maryland. This information is needed to evaluate the availability and chemical quality of water from these aquifers. The Patuxent aquifer unconformably overlies Lower Paleozoic and Precambrian basement rocks and consists primarily of medium- to coarse-grained quartz sand. Discontinuous lenses of gravel and silty clay are commonly interbedded with the sand-sized material. The Patuxent aquifer in this area attains a thickness of 250 feet and transmissivities range from 2,000 to 8,000 feet squared per day. The Patuxent is the most productive source of ground water in the Baltimore area. In 1982, approximately 11 million gallons of water per day was produced from this unit. Several cones of depression, ranging from 30 to 50 feet below sea level, have developed in response to this pumping stress. The Arundel Formation conformably overlies the Patuxent aquifer. The Arundel is composed predominantly of clay and ranges from 0 to 150 feet thick. The Arundel exhibits very low vertical hydraulic conductivities that are on the order of 10 -9 to 10 -11 feet per second. This unit acts as the upper confining bed of the Patuxent aquifer in much of the project area. The Patapsco aquifer unconformably overlies the Arundel Formation and is a medium- to fine-grained quartz sand. The Patapsco functions as a water-table aquifer in much of the project area. Although the Patapsco has been heavily pumped in the past, pumpage from that aquifer in Baltimore was negligible in 1982. Brackish-water contamination of the Patuxent and Patapsco aquifers has been a major water-quality problem since the early 1900's. The Patuxent aquifer presently (1982) contains a circular plume of brackish-water contamination about 5 miles in diameter. This plume is centered on the Harbor district and has enlarged measurably since 1945. The Patapsco aquifer has a smaller zone of brackish-water contamination that has decreased in size since 1945. Borehole data demonstrate that the Arundel Formation has been breached by Pleistocene river channels near the Harbor district. These erosional channels provide a conduit for brackish water to intrude into the Patuxent aquifer. A two-dimensional areal solute-transport model of the Patuxent aquifer was constructed. This model was designed to estimate the future movement of the brackish-water plume based on alternative scenarios of aquifer use. Model simulations suggest that the plume will remain relatively immobile if 1982 pumping patterns continue into the foreseeable future. However, increased pumpage in the Marley Neck peninsula could draw the plume to the southeast and increase contamination of the Fairfield area. The water quality of the Patuxent aquifer is extremely variable. Because of this variability, it is useful to divide the aquifer into three water-quality zones: Zone 1 -- This zone corresponds to the plume of brackish-water contamination. Zone 2 -- This zone exhibits evidence of urbanization-related contamination such as measurable concentrations of organic chemicals and elevated concentrations of trace elements and total organic carbon. Zone 3 -- Water composition in this zone is controlled exclusively by naturally occurring chemical processes. These processes are dominated by reactions involving dissolved iron. Near the outcrop area, oxidation of pyrite and lignite consumes dissolved oxygen and produces ledges of iron hydroxide-cemented sandstones and conglomerates. The predominant dissolved iron species in oxic water is Fe(OH) 2 + . Downgradient, the water becomes anoxic and sulfate reduction becomes an important process. The predominant dissolved iron species in anoxic water is Fe 2+ .

Maryland↗

Development of regression equations to estimate flow durations and low-flow-frequency statistics in New Hampshire streams

Regression equations and basin-characteristic digital datasets were developed to help water-resource managers estimate surface-water resources during periods of low flow in New Hampshire. The regression equations were developed to estimate statistics for the seasonal and annual low-flow-frequency and seasonal period-of-record and period-of-record flow durations. Because streamflow is maintained by ground-water discharge during periods of low flow, these equations also will aid in the assessment of ground-water availability. Ultimately, the equations and datasets developed herein can be combined with data on water withdrawals, discharges, and interbasin transfers in a geographic information system (GIS) to allow assessments of water use and water availability in any drainage basin in the State of New Hampshire. Regression equations developed in this study provide estimates of the seasonal (spring, summer, fall, and winter) and annual 7-day 2-year (7Q2) and 7-day 10-year (7Q10) low-flow-frequency values, as well as seasonal period-of-record and period-of-record flow durations (60-, 70-, 80-, 90-, 95-, and 98-percent exceedences) for ungaged reaches of unregulated New Hampshire streams. Regression equations were developed using seasonal and annual low-flow statistics from 58 to 60 continuous-record stream-gaging stations in New Hampshire and nearby areas in neighboring states, and measurements of various characteristics of the drainage basins that contribute flow to those stations. The estimating equations for the seasonal and annual 7Q2 and 7Q10 values were developed using generalized-least-squares (GLS) regression analyses. The GLS equations developed for these flow statistics gave average prediction errors that ranged from 11 to 61 percent. The estimating equations for flow-duration exceedence frequency values were developed using ordinary-least-squares (OLS) regression analysis. The OLS equations developed for these flow statistics gave average prediction errors ranging from 14 to 79 percent. A total of 93 measurable drainage-basin characteristics were selected as possible predictor variables. Of these 93 variables, the following 10 were determined to be statistically significant predictors for at least one of the dependent variables: drainage area, average basin slope, maximum basin elevation, average summer gage precipitation for 1961-90, average spring gage precipitation for 1961-90, average mean annual basin temperature for 1961-90, average mean summer basin temperature for 1961-90, average winter basin-centroid precipitation for 1961-90, percent of the basin that is coniferous, and percent of the basin that is mixed coniferous and deciduous. These 10 basin characteristics were selected because they were statistically significant based on several statistical parameters that evaluated which combination of characteristics contributed the most to the predictive accuracy of the regression-equation models. A GIS is required to measure the values of the predictor variables for the equations developed in this study.

New Hampshire↗

Numerical simulation of ground-water withdrawals in the southern Lihue Basin, Kauai, Hawaii

Numerical simulations indicate that ground-water withdrawals from the Hanamaulu and Puhi areas of the southern Lihue Basin will result in a decline in water levels and reductions in base flows of streams near proposed new water-supply wells. Most of the changes will be attained within 10 to 20 years of the start of pumping. Except for areas such as Puhi and Kilohana, the freshwater lens in most inland areas of the southern Lihue Basin is thick and model simulations indicate that changes in water level and the position of the freshwater- saltwater interface in response to pumping will be small relative to the present thickness of the freshwater lens. Effects of the proposed withdrawals on streamflow depend on withdrawal rate and proximity of the wells to streams. Placing pumped wells away from streams with low base flow and toward streams with high base flow can reduce the relative effect on individual streams. Simulation of the 0.42-million-gallon-per-day increase in withdrawal projected for 2000 indicates that the resulting changes in water levels and interface position, relative to conditions prior to the withdrawal increase, will be small, and that stream base flow will be reduced by less than 10 percent. Simulation of the 0.83-million-gallon-per-day withdrawal projected for 2010 indicates further thinning of the freshwater lens in the Puhi area, where the lens already may be thin, as well as base-flow reduction in Nawiliwili Stream. Simulation of an alternative distribution of the 0.83-million-gallon-per-day withdrawal indicates that the effects can be reduced by shifting most of the new withdrawal to the Hanamaulu area where the freshwater lens is thicker and stream base flows are greater. Simulation of the 1.16-million-gallon-per-day increase in withdrawal projected for 2020 indicates that if withdrawal is distributed only among Hana-maulu wells 1, 3, and 4, and Puhi well 5A, further thinning of the already-thin freshwater lens in the Puhi area would occur. Such a distribution would also exceed the maximum draft recommended by the water-systems standards used in Hawaii. Another simulation in which part of the 1.16 million gallons per day was distributed among three additional hypothetical wells in the Hanamaulu area showed that the pumping effects could be shifted from the Puhi area to the Hanamaulu area, where the freshwater lens is thicker, but that base flow in Hanamaulu Stream may decrease by as much as 16 percent.

Hawaii↗

Simulation of subsurface storage and recovery of treated effluent injected in a saline aquifer, St. Petersburg, Florida

The potential for subsurface storage and recovery of treated effluent into the uppermost producing zone (zone A) of the Upper Floridan aquifer in St. Petersburg, Florida, is being studied by the U.S. Geological Survey, in cooperation with the city of St. Petersburg and the Southwest Florida Water Management District. A measure of the success of this practice is the recovery efficiency, or the quantity of water relative to the quantity injected, that can be recovered before the water that is withdrawn fails to meet water-quality standards. The feasibility of this practice will depend upon the ability of the injected zone to receive, store, and discharge the injected fluid. A cylindrical model of ground-water flow and solute transport, incorporating available data on aquifer properties and water quality, was developed to determine the relation of recovery efficiency to various aquifer and fluid properties that could prevail in the study area. The reference case for testing was a base model considered representative of the saline aquifer underlying St. Petersburg. Parameter variations in the tests represent possible variations in aquifer conditions in the area. The model also was used to study the effect of various cyclic injection and withdrawal schemes on the recovery efficiency of the well and aquifer system. A base simulation assuming 15 days of injection of effluent at a rate of 1.0 million gallons per day and 15 days of withdrawal at a rate of 1.0 million gallons per day was used as reference to compare changes in various hydraulic and chemical parameters on recovery efficiency. A recovery efficiency of 20 percent was estimated for the base simulation. For practical ranges of hydraulic and fluid properties that could prevail in the study area, the model analysis indicates that (1) the greater the density contrast between injected and resident formation water, the lower the recovery efficiency, (2) recovery efficiency decreases significantly as dispersion increases, (3) high formation permeability favors low recovery efficiencies, and (4) porosity and anisotropy have little effect on recovery efficiencies. In several hypothetical tests, the recovery efficiency fluctuated between about 4 and 76 percent. The sensitivity of recovery efficiency to variations in the rate and duration of injection (0.25, 0.50, 1.0, and 2.0 million gallons per day) and withdrawal cycles (60, 180, and 365 days) was determined. For a given operational scheme, recovery efficiency increased as the injection and withdrawal rate is increased. Model results indicate that recovery efficiencies of between about 23 and 37 percent can be obtained for different subsurface storage and recovery schemes. Five successive injection, storage, and recovery cycles can increase the recovery efficiency to about 46 to 62 percent. There is a larger rate of increase at smaller rates than at larger rates. Over the range of variables studied, recovery efficiency improved with successive cycles, increasing rapidly during initial cycles tyhen more slowly at later cycles. The operation of a single well used for subsurface storage and recovery appears to be technically feasible under moderately favorable conditions; however, the recovery efficiency is higly dependent upon local physical and operational parameters. A combination of hydraulic, chemical, and operational parameters that minimize dispersion and buoyancy flow, maximizes recovery efficiency. Recovery efficiency was optimal where resident formation water density and permeabilities were relatively similar and low.

Water-Resources Investigations Report↗

Quality characteristics of ground water in the Ozark aquifer of northwestern Arkansas, southeastern Kansas, southwestern Missouri and northeastern Oklahoma, 2006-07

Because of water quantity and quality concerns within the Ozark aquifer, the State of Kansas in 2004 issued a moratorium on most new appropriations from the aquifer until results were made available from a cooperative study between the U.S. Geological Survey and the Kansas Water Office. The purposes of the study were to develop a regional ground-water flow model and a water-quality assessment of the Ozark aquifer in northwestern Arkansas, southeastern Kansas, southwestern Missouri, and northeastern Oklahoma (study area). In 2006 and 2007, water-quality samples were collected from 40 water-supply wells completed in the Ozark aquifer and spatially distributed throughout the study area. Samples were analyzed for physical properties, dissolved solids and major ions, nutrients, trace elements, and selected isotopes. This report presents the results of the water-quality assessment part of the cooperative study. Water-quality characteristics were evaluated relative to U.S. Environmental Protection Agency drinking-water standards. Secondary Drinking-Water Regulations were exceeded for dissolved solids (11 wells), sulfate and chloride (2 wells each), fluoride (3 wells), iron (4 wells), and manganese (2 wells). Maximum Contaminant Levels were exceeded for turbidity (3 wells) and fluoride (1 well). The Maximum Contaminant Level Goal for lead (0 milligrams per liter) was exceeded in water from 12 wells. Analyses of isotopes in water from wells along two 60-mile long ground-water flow paths indicated that water in the Ozark aquifer was at least 60 years old but the upper age limit is uncertain. The source of recharge water for the wells along the flow paths appeared to be of meteoric origin because of isotopic similarity to the established Global Meteoric Water Line and a global precipitation relation. Additionally, analysis of hydrogen-3 (3H) and carbon-14 (14C) indicated that there was possible leakage of younger ground water into the lower part of the Ozark aquifer. This may be caused by cracks or fissures in the confining unit that separates the upper and lower parts of the aquifer, poorly constructed or abandoned wells, or historic mining activities. Analyses of major ions in water from wells along the flow paths indicated a transition from freshwater in the east to saline water in the west. Generally, ground water along flow paths evolved from a calcium magnesium bicarbonate type to a sodium calcium bicarbonate or a sodium calcium chloride bicarbonate type as water moved from recharge areas in Missouri into Kansas. Much of this evolution occurred within the last 20 to 25 miles of the flow paths along a water-quality transition zone near the Kansas-Missouri State line and west. The water quality of the Kansas part of the Ozark aquifer is degraded compared to the Missouri part. Geophysical and well-bore flow information and depth-dependent water-quality samples were collected from a large-capacity (1,900-2,300 gallons per minute) municipal-supply well to evaluate vertical ground-water flow accretion and variability in water-quality characteristics at different levels. Although the 1,050-foot deep supply well had 500 feet of borehole open to the Ozark aquifer, 77 percent of ground-water flow entering the borehole came from two 20-foot thick rock layers above the 1,000-foot level. For the most part, water-quality characteristics changed little from the deepest sample to the well-head sample, and upwelling of saline water from deeper geologic formations below the well was not evident. However, more saline water may be present below the bottom of the well.

Arkansas, Kansas, Missouri, Oklahoma↗

Simulated response of the Sparta Aquifer to outcrop area recharge augmentation, southeastern Arkansas

Recharge augmentation by construction of infiltration impoundments is a potential means of increasing aquifer water levels and aquifer yield that is under consideration for the Sparta aquifer in southeastern Arkansas. The aquifer is a major water resource for municipal, industrial, and agricultural uses, and approximately 287 million gallons per day was pumped from the aquifer in Arkansas in 1995; this is double the amount pumped in 1975. Historically, the Sparta aquifer has provided abundant water of high quality. In recent years, however, the demand for water in some areas has resulted in withdrawals from the Sparta that significantly exceed recharge to the aquifer, and considerable declines have occurred in the potentiometric surface. To better manage the Sparta aquifer, water users in Arkansas are evaluating and implementing a variety of management practices and assessing alternative, surface-water sources to reduce stress upon the Sparta aquifer. One approach to managing and maximizing use of the Sparta aquifer is augmenting recharge to the aquifer by construction of infiltration lakes or canals within the recharge area. The basic concept of augmented recharge is simply to increase the amount of water being introduced into the aquifer so that more water will be available for use. Ground-water flow model simulations were conducted to assess the effectiveness of constructing lakes or canals to augment recharge. Results show that construction of five new lakes in the Sparta recharge area upgradient from major pumping centers or construction of a series of canals along the length of the recharge area yield notable benefit to aquifer conditions when compared with simulations entailing no augmentation of recharge. Augmentation of recharge in the Sparta aquifer with emplacement of lakes provides slight increase to aquifer water levels. The presence of the lakes increased simulated aquifer water levels 0.5 foot or more across a broad area comprising all or a substantial part of 19 counties after the 30-year simulation period. Substantial increases of 5 feet or greater are limited to a smaller area proximal to the lakes. Increases of 5 feet or more are seen in El Dorado, Pine Bluff, and Stuttgart. The positive effect of the lakes on aquifer water levels is rapidly realized after emplacement of the lakes. For example, in the El Dorado area more than 3 feet of a total of 8 feet of water-level increase is seen in the first 5 years of the simulation; in the Pine Bluff area 9 feet of a total of 16 feet of increase occurs within 5 years. Sustainable yield from the aquifer could be expected to be increased within the zone of influence of the lakes. Augmentation of recharge in the Sparta aquifer with emplacement of canals provides considerable increase of aquifer water levels. The zone of influence in the aquifer with canal-augmented recharge extends from the recharge area eastward to the Mississippi River. Aquifer water levels exhibit an increase of 5 feet or more across a broad area comprising all or a substantial part of 15 counties. Increases of 20 feet or more are seen in El Dorado, Pine Bluff, and Stuttgart. The amount of water moving into the aquifer is substantially increased under this scenario, and the amount of water removed from storage is decreased, thereby, increasing aquifer conditions considerably. Sustainable yield from the aquifer could be expected to be greater within the zone of influence of the canals as compared to either the scenario without recharge augmentation or recharge augmentation with lakes. The effect of the canal on aquifer water levels is rapidly realized after emplacement of the canals. For example, in the El Dorado area, 22 feet of a total of 30 feet of increase is seen in the first 5 years of the simulation; in the Pine Bluff area, 15 feet of a total of 24 feet of increase occurs within 5 years. As constructed, the model simulations imply that any lakes or canals constructed would maintain exce

Arkansas↗

Secondary Mineral Deposits and Evidence of Past Seismicity and Heating of the Proposed Repository Horizon at Yucca Mountain, Nevada

The Drift Degradation Analysis (DDA) (BSC, 2003) for the proposed high-level radioactive waste repository at Yucca Mountain, Nevada, describes model simulations of the effects of pre- and post-closure seismicity and waste-induced heating on emplacement drifts. Based on probabilistic seismic hazard analyses of the intensity and frequency of future seismic events in the region (CRWMS M&O, 1998), the DDA concludes that future seismicity will lead to substantial damage to emplacement drifts, particularly those in the lithophysal tuffs, where some simulations predict complete collapse of the drift walls. Secondary mineral studies conducted by the U.S. Geological Survey since 1995 indicate that secondary calcite and silica have been deposited in some fractures and lithophysal cavities in the unsaturated zone (UZ) at Yucca Mountain during at least the past 10 million years (m.y.), and probably since the tuffs cooled to less than 100?C. Tuff fragments, likely generated by past seismic activity, have commonly been incorporated into the secondary mineral depositional sequences. Preliminary observations indicate that seismic activity has generated few, if any, tuff fragments during the last 2 to 4 m.y., which may be inconsistent with the predictions of drift-wall collapse described in the DDA. Whether or not seismicity-induced tuff fragmentation occurring at centimeter to decimeter scales in the fracture and cavity openings relates directly to failure of tuff walls in the 5.5-m-diameter waste emplacement drifts, the deposits do provide a potential record of the spatial and temporal distribution of tuff fragments in the UZ. In addition, the preservation of weakly attached coatings and (or) delicate, upright blades of calcite in the secondary mineral deposits provides an upper limit for ground motion during the late stage of deposition that might be used as input to future DDA simulations. Finally, bleaching and alteration at a few of the secondary mineral sites indicate that they were subjected to heated gases at approximately the temperatures expected from waste emplacement. These deposits provide at least limited textural and mineralogic analogs for waste-induced, high-humidity thermal alteration of emplacement drift wall rocks.

Water-Resources Investigations Report↗

Hydrogeology of Two Areas of the Tug Hill Glacial-Drift Aquifer, Oswego County, New York

Two water-production systems, one for the Village of Pulaski and the other for the Villages of Sandy Creek and Lacona in Oswego County, New York, withdraw water from the Tug Hill glacial-drift aquifer, a regional sand and gravel aquifer along the western flank of the Tug Hill Plateau, and provide the sole source of water for these villages. As a result of concerns about contamination of the aquifer, two studies were conducted during 2001 to 2004, one for each water-production system, to refine the understanding of ground-water flow surrounding these water-production systems. Also, these studies were conducted to determine the cause of the discrepancy between ground-water ages estimated from previously constructed numerical ground-water-flow models for the Pulaski and Sandy Creek/Lacona well fields and the apparent ground-water ages determined using concentrations of tritium and chlorofluorocarbons. The Village of Pulaski withdrew 650,000 gallons per day in 2000 from four shallow, large-diameter, dug wells finished in glaciolacustrine deposits consisting of sand with some gravelly lenses 3 miles east of the village. Four 2-inch diameter test wells were installed upgradient from each production well, hydraulic heads were measured, and water samples collected and analyzed for physical properties, inorganic constituents, nutrients, bacteria, tritium, dissolved gases, and chlorofluorocarbons. Recharge to the Tug Hill glacial-drift aquifer is from precipitation directly over the aquifer and from upland sources in the eastern part of the recharge area, including (1) unchannelized runoff from till and bedrock hills east of the aquifer, (2) seepage to the aquifer from streams that drain the Tug Hill Plateau, (3) ground-water inflow from the till and bedrock on the adjoining Tug Hill Plateau. Water-quality data collected from four piezometers near the production wells in November 2003 indicated that the water is a calcium-bicarbonate type with iron concentrations that slightly exceeded the U.S. Environmental Protection Agency?s Secondary Maximum Contaminant Level in three of the four samples. The relatively small concentrations of major ions and nutrients in the samples indicate that there is no contamination from septic-tank effluent and dissolved road salt in the ground water at the Village of Pulaski well field. Three of the four samples were analyzed for total coliform bacteria and Escherichia coli (E. coli), and only total coliform bacteria were detected in all three samples. E. coli was not detected in any samples. The Villages of Sandy Creek and Lacona use about 270,000 gallons of water per day for public consumption?alternating withdrawals from northern and southern well fields located in glaciolacustrine beach, glaciofluvial outwash, and alluvial inwash deposits consisting mostly of silty sand and gravel. Four test wells were drilled, hydraulic heads were measured, and water samples collected between 2001 and 2003 and analyzed for physical properties, inorganic constituents, nutrients, bacteria, tritium, dissolved gases, and chlorofluorocarbons. The aquifer in the Sandy Creek/Lacona area is highly susceptible to contamination because the aquifer (1) is unconfined, (2) is highly transmissive, (3) is thin (10 to 25 feet) and narrow (about 0.5 miles wide), and (4) has relatively short ground-water flowpaths from recharge to discharge areas (including wells). Additionally, drainage ditches east of the southern well field intercept westward-flowing ground water, which then is routed to an area just upgradient from the production wells where some of the water seeps back into the aquifer and is captured by these wells, effectively shortcircuiting the ground-water-flow system. Water-quality samples collected from three wells in the Sandy Creek/Lacona southern well field indicate that the ground water is a sodium-bicarbonate/sulfate type and of good quality; however, in one water sample, the sodium concentration of 6

Scientific Investigations Report↗

Evapotranspiration Rates of Riparian Forests, Platte River, Nebraska, 2002-06

Evapotranspiration (ET) in riparian areas is a poorly understood component of the regional water balance in the Platte River Basin, where competing demands have resulted in water shortages in the ground-water/surface-water system. From April 2002 through March 2006, the U.S. Geological Survey, Nebraska Platte River Cooperative Hydrology Study Group, and Central Platte Natural Resources District conducted a micrometeorological study of water and energy balances at two sites in central Nebraska near Odessa and Gothenburg to improve understanding of ET rates and factors affecting them in Platte River riparian forests. A secondary objective of the study was to constrain estimates of ground-water use by riparian vegetation to satisfy ET consumptive demands, a useful input to regional ground-water flow models. Both study sites are located on large islands within the Platte River characterized by a cottonwood-dominated forest canopy on primarily sandy alluvium. Although both sites are typical of riparian forests along the Platte River in Nebraska, the Odessa understory is dominated by deciduous shrubs, whereas the Gothenburg understory is dominated by eastern redcedars. Additionally, seasonal ground-water levels fluctuated more at Odessa than at Gothenburg. The study period of April 2002 through March 2006 encompassed precipitation conditions ranging from dry to wet. This study characterized the components of the water balance in the riparian zone of each site. ET was evaluated from eddy-covariance sensors installed on towers above the forest canopy at a height of 26.1 meters. Precipitation was measured both above and below the forest canopy. A series of sensors measured soil-moisture availability within the unsaturated zone in two different vertical profiles at each site. Changes in ground-water altitude were evaluated from piezometers. The areal footprint represented in the water balance extended up to 800 meters from each tower. During the study, ET was less variable than precipitation. Annual ET fluctuated about 7 percent from the 4-year mean, ranging from about 514 to 586 millimeters per year (551 on average) at the Odessa site and 535 to 616 millimeters per year (575 on average) at the Gothenburg site. Conversely, annual precipitation fluctuated by about 35 percent from the 4-year mean, ranging from 429 to 844 millimeters per year at Odessa and 359 to 791 millimeters per year at Gothenburg. Of this precipitation, 14 to 15 percent was intercepted by the forest canopy before it could infiltrate into the soil. For the 4-year period, annual ground-water recharge from the riparian measurement zone averaged 76 and 13 millimeters at Odessa and Gothenburg, respectively, to satisfy the water balance at each site. This indicates that, from an annual perspective, ground-water reductions caused by ET may be minimal. This effect varied somewhat and primarily was affected by fluctuations in precipitation. Ground-water discharge occurred during the driest study year (2002), whereas ground-water recharge occurred from 2003 to 2005. These results do not exclude ground water as an important source of water to riparian vegetation - especially to phreatophytes that have the capability of directly using water from the saturated zone - during periods of high ET in the summer, particularly during periods of lower than normal precipitation. However, the calculations indicate that, on an annual (or longer) net-flux basis, ground-water use by riparian forests is likely to be balanced by periods of recharge from excess precipitation at other times of the year. In contrast to more arid settings, where scientific literature indicates that ground water may supply a large fraction of the water used for ET by riparian vegetation, precipitation along the Platte River of Nebraska was great enough - and generally greater than ET - that most or all of the annual ET demand was satisfied by available precipitation. Crop coefficients developed for 15-

Scientific Investigations Report↗

Simulation of the water-table altitude in the Biscayne Aquifer, southern Dade County, Florida, water years 1945-89

A digital model of the flow system in the highly permeable surficial Biscayne aquifer of southern Dade County, Florida, was constructed for the purposes of better understanding processes that influence the flow system and of supporting the construction of a subregional model of the transport of brackish water from a flowing artesian well. Problems that needed resolution in this endeavor included the development of methods to represent the influence of flowing surface water in seasonally inundated wetlands and the influence of a network of controlled canals developed in stages during the simulation time period (water years 1945-89). An additional problem was the general lack of natural aquifer boundaries near the boundaries of the study area. The model construction was based on a conceptual description of the Biscayne aquifer developed from the results of previous U.S. Geological Survey investigations. Modifications were made to an existing three-dimensional finite-difference simulator of ground-water flow to enable an upper layer of the grid to represent seasonally occurring overland sheetflow in a series of transient simulations of water levels from 1945 to 1989. A rewetting procedure was developed for the simulator that permitted resaturation of cells in this layer when the wet season recurred. An "equivalent hydraulic conductivity" coefficient was assigned to the overland flow layer that was analogous, subject to various approximations, to the use of the Manning equation. The surficial semiconfining peat and marl layers, levees, canals, and control structures were also represented as part of the model grid with the appropriate choices of hydraulic coefficient values. For most of the Biscayne aquifer grid cells, the value assigned to hydraulic conductivity for model calibration was 30,000 feet per day and the value assigned to porosity was 20 percent. Boundary conditions were specified near data sites having long-term records of surface-water stages or water-table altitudes, and modifications to the simulator permitted the specification of time-varying pressures at boundary grid cells. Rainfall data from a station in Homestead generally were used as an areally uniform rainfall specification throughout the modeled region. Maximum evapotranspiration rates ranged seasonally from a minimum of 0.08 inch per day in January to a maximum of 0.21 inch per day between June and October. Shallow-root and deep-root zone depths for the evapotranspiration calculation were 3 and 20 feet in the coastal ridge and were 0.10 and 5 feet in the glades regions where peat and marl covers occurred. Results of sensitivity analyses indicated that the simulations of stages and water levels were relatively unresponsive to 50 percent changes in aquifer hydraulic conductivity, porosity, and the equivalent hydraulic conductivity of overland flow. However, 20 percent changes in rainfall and maximum evapotranspiration rates produced significantly different water levels, as did interchange of coastal ridge and glades deep-root zone (extinction) depths. Water levels were simulated very well at most measurement sites. Sensitivity analyses illustrated the significant influence of the uncontrolled agricultural drainage canals on pre-1968 regional water levels and the further influence of Black Creek Canal in draining a region of high water after 1961. Other analyses indicated that the flood-control system of 1968-82 lowered peak water levels in the affected region by as much as 1.5 feet in the wet summers of 1968, 1969, and 1981, and that Levee 67 Extended channeled flows from the S-12 spillway structures and raised overland flow stages in Shark River Slough. Hypothetical scenarios of well-field pumping in the vicinity of Levee 31N indicated that the pumping induced a significant amount of recharge from the adjacent borrow canal, the degree of which depended on the distance between the canal and the well field. The computed ratio of evapotranspiration to rainfall recharge ranged from 88 to 94 percent during water years 1945-82. The ratio increased to about 97.9 percent during water years 1983-89, possibly because of changing water-management practices and deficient rainfall.

Florida↗