Geology Reports⌕ Search

SEARCH · Geology Reports

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 1,117 records · Page 62Linked to original sources

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↗

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 then 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 highly 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.

Florida↗

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↗

Simulation of runoff and water quality for 1990 and 2008 land-use conditions in the Reedy Creek watershed, east-central Florida

Hydrologic and water-quality data have been collected within the 177-square-mile Reedy Creek, Florida, watershed, beginning as early as 1939, but the data have not been used to evaluate relations among land use, hydrology, and water quality. A model of the Reedy Creek watershed was developed and applied to the period January 1990 to December 1995 to provide a computational foundation for evaluating the effects of future land-use changes on hydrology and water quality in the watershed. The Hydrological Simulation Program-Fortran (HSPF) model was used to simulate hydrology and water quality of runoff for pervious land areas, impervious land areas, and stream reaches. Six land-use types were used to characterize the hydrology and water quality of pervious and impervious land areas in the Reedy Creek watershed: agriculture, rangeland, forest, wetlands, rapid infiltration basins, and urban areas. Hydrologic routing and water-quality reactions were simulated to characterize hydrologic and water-quality processes and the movement of runoff and its constituents through the main stream channels and their tributaries. Because of the complexity of the stream system within the Reedy Creek Improvement District (RCID) (hydraulic structures, retention ponds) and the anticipated difficulty of modeling the system, an approach of calibrating the model parameters for a subset of the gaged watersheds and confirming the usefulness of the parameters by simulating the remainder of the gaged sites was selected for this study. Two sub-watersheds (Whittenhorse Creek and Davenport Creek) were selected for calibration because both have similar land use to watersheds within the RCID (with the exception of urban areas). Given the lack of available rainfall data, the hydrologic calibration of the Whittenhorse Creek and Davenport Creek sub-watersheds was considered acceptable (for monthly data, correlation coefficients, 0.86 and 0.88, and coefficients of model-fit efficiency, 0.72 and 0.74, respectively). The hydrologic model was tested by applying the parameter sets developed for Whittenhorse Creek and Davenport Creek to other land areas within the Reedy Creek watershed, and by comparing the simulated results to observed data sets for Reedy Creek near Vineland, Bonnet Creek near Vineland, and Reedy Creek near Loughman. The hydrologic model confirmation for Reedy Creek near Vineland (correlation coefficient, 0.91, and coefficient of model fit efficiency, 0.78, for monthly flows) was acceptable. Flows for Bonnet Creek near Vineland were substantially under simulated. Consideration of the ground-water contribution to Bonnet Creek could improve the water balance simulation for Bonnet Creek near Vineland. On longer time scales (monthly or over the 72-month simulation period), simulated discharges for Reedy Creek near Loughman agreed well with observed data (correlation coefficient, 0.88). For monthly flows the coefficient of model-fit efficiency was 0.77. On a shorter time scale (less than a month), however, storm volumes were greatly over simulated and low flows (less than 8 cubic feet per second) were greatly under simulated. A primary reason for the poor results at low flows is the diversion of an unknown amount of water from the RCID at the Bonnet Creek near Kissimmee site. Selection of water-quality constituents for simulation was based primarily on the availability of water-quality data. Dissolved oxygen, nitrogen, and phosphorus species were simulated. Representation of nutrient cycling in HSPF also required simulation of biochemical oxygen demand and phytoplankton populations. The correlation coefficient for simulated and observed daily mean dissolved oxygen concentration values at Reedy Creek near Vineland was 0.633. Simulated time series of total phosphorus, phosphate, ammonia nitrogen, and nitrate nitrogen generally agreed well with periodically observed values for the Whittenhorse Creek and Davenport Creek sites. Simulated water-quality constituents at the Bonnet Creek and Reedy Creek near Vineland sites varied as to how well the values agreed with periodically observed constituent concentrations. Simulated water-quality constituent concentrations for the Reedy Creek near Loughman site generally agreed well with observed constituent concentrations. Simulation of a future land-use scenario for the Reedy Creek watershed was based on the hydrologic and water-quality simulations, projected 2008 land use within the RCID, and assuming no change in existing land use for other areas within the Reedy Creek watershed but external to the RCID. The percentages of forest and urban-impervious land use showed the most change between existing and future land use; forest areas decreased by 50 percent and urban-impervious areas increased by 300 percent. Simulated values of mean total phosphorus, phosphate, ammonia nitrogen, and nitrate nitrogen concentrations for existing and future land-use simulations were within 0.01 milligrams per liter of each other. The simulated maximum daily load increased an average of 10 percent for all constituents. Maximum daily nitrate nitrogen load increased about 17 percent, the greatest increase of all daily constituent loads. Duration curves of daily total phosphorus, phosphate, ammonia nitrogen, and nitrate nitrogen load indicated an increase in the likelihood of exceeding a given load throughout the range of daily constituent loads at Reedy Creek near Loughman.

Florida↗

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↗

Flow-system analysis of the Madison and Minnelusa aquifers in the Rapid City area, South Dakota — Conceptual model

The conceptual model of the Madison and Minnelusa aquifers in the Rapid City area synthesizes the physical geography, hydraulic properties, and ground-water flow components of these important aquifers. The Madison hydrogeologic unit includes the karstic Madison aquifer, which is defined as the upper, more permeable 100 to 200 ft of the Madison Limestone, and the Madison confining unit, which consists of the lower, less permeable part of the Madison Limestone and the Englewood Formation. Overlying the Madison hydrogeologic unit is the Minnelusa hydrogeologic unit, which includes the Minnelusa aquifer in the upper, more permeable 200 to 300 ft and the Minnelusa confining unit in the lower, less permeable part. The Madison and Minnelusa hydrogeologic units outcrop in the study area on the eastern flank of the Black Hills where recharge occurs from streamflow losses and areal recharge. The conceptual model describes streamflow recharge, areal recharge, ground-water flow, storage in aquifers and confining units, unsaturated areas, leakage between aquifers, discharge from artesian springs, and regional outflow. Effective transmissivities estimated for the Madison aquifer range from 500 to 20,000 ft 2 /d and for the Minnelusa aquifer from 500 to 10,000 ft 2 /d. Localized anisotropic transmissivity in the Madison aquifer has tensor ratios as high as 45:1. Vertical hydraulic conductivities for the Minnelusa confining unit determined from aquifer tests range from 1.3x10 -3 to 3.0x10 -1 ft/d. The confined storage coefficient of the Madison and Minnelusa hydrogeologic units was estimated as 3x10 -4 ft/d. Specific yield was estimated as 0.09 for the Madison and Minnelusa aquifers and 0.03 for the Madison and Minnelusa confining units. Potentiometric surfaces for the Madison and Minnelusa aquifers have a general easterly gradient of about 70 ft/mi with local variations. Temporal change in hydraulic head in the Madison and Minnelusa aquifers ranged from about 5 to 95 ft in water years 1988-97. The unconfined areas were estimated at about 53 and 36 mi 2 for the Madison and Minnelusa hydrogeologic units, respectively, in contrast to an aquifer analysis area of 629 mi 2 . Dye-tracer tests, stable isotopes, and hydrogeologic features were analyzed conjunctively to estimate generalized ground-water flowpaths in the Madison aquifer and their influences on the Minnelusa aquifer. The western Rapid City area between Boxelder Creek and Spring Creek was characterized as having undergone extensive tectonic activity, greater brecciation in the Minnelusa Formation, large transmissivities, generally upward hydraulic gradients from the Madison aquifer to the Minnelusa aquifer, many karst springs, and converging flowpaths. Water-budget analysis included: (1) a dry-period budget for declining water levels; October 1, 1987, to March 31, 1993; (2) a wet-period budget for rising water levels, April 1, 1993, to September 30, 1997; and (3) a full 10-year period budget for water years 1988-97. By simultaneously balancing these water budgets, initial estimates of recharge, discharge, change in storage, and hydraulic properties were refined. Inflow rates for the 10-year budget included streamflow recharge of about 45 ft 3 /s or 61 percent of the total budget and areal recharge of 22 ft 3 /s or 30 percent. Streamflow recharge to the Madison hydrogeologic unit was about 86 percent of the total streamflow recharge. Outflow for the 10-year budget included springflow of 31 ft 3 /s or 42 percent of the total budget, water use of about 10 ft 3 /s or 14 percent, and regional outflow of 22 ft 3 /s or 30 percent. Ground-water storage increased 9 ft 3 /s during the 10-year period, and net ground-water movement from the Madison to Minnelusa hydrogeologic unit was about 8 ft 3 /s.

South Dakota↗

Terrestrial Ecosystems - Topographic Moisture Potential of the Conterminous United States

As part of an effort to map terrestrial ecosystems, the U.S. Geological Survey has generated topographic moisture potential classes to be used in creating maps depicting standardized, terrestrial ecosystem models for the conterminous United States, using an ecosystems classification developed by NatureServe. A biophysical stratification approach, developed for South America and now being implemented globally, was used to model the ecosystem distributions. Substrate moisture regimes strongly influence the differentiation and distribution of terrestrial ecosystems, and therefore topographic moisture potential is one of the key input layers in this biophysical stratification. The method used to produce these topographic moisture potential classes was based on the derivation of ground moisture potential using a combination of computed topographic characteristics (CTI, slope, and aspect) and mapped National Wetland Inventory (NWI) boundaries. This method does not use climate or soil attributes to calculate relative topographic moisture potential since these characteristics are incorporated into the ecosystem model though other input layers. All of the topographic data used for this assessment were derived from the USGS 30-meter National Elevation Dataset (NED ) including the National Compound Topographic Index (CTI). The CTI index is a topographically derived measure of slope for a raster cell and the contributing area from upstream raster cells, and thus expresses potential for water flow to a point. In other words CTI data are 'a quantification of the position of a site in the local landscape', where the lowest values indicate ridges and the highest values indicate stream channels, lakes and ponds. These CTI values were compared to independent estimates of water accumulation by obtaining geospatial data from a number of sample locations representing two types of NWI boundaries: freshwater emergent wetlands and freshwater forested/shrub wetlands. Where these shorelines (the interface between the NWI wetlands and adjacent land) occurred, the CTI values were extracted and a histogram of their statistical distributions was calculated. Based on an evaluation of these histograms, CTI thresholds were developed to separate periodically saturated or flooded land, mesic uplands (moderately moist), and uplands. After the range of CTI values for these three different substrate moisture regimes was determined, the CTI values were grouped into three initial topographic moisture potential classes. As a final step in the generation of this national data layer, the uplands classification was subdivided into either very dry uplands or dry uplands. Very dry uplands were defined as uplands with relatively steep, south-facing slopes, and identification of this class was based on the slope and aspect datasets derived from the NED. The remaining uplands that did not meet these additional criteria were simply re-classified as dry uplands. The final National Topographic Moisture Potential dataset for the conterminous United States contains four classes: periodically saturated or flooded land (CTI = 18.5), mesic uplands (12 =< CTI < 18.5), dry uplands (CTI < 12), and very dry uplands (CTI < 12, Slope > 24 degrees and 91 degrees =< Aspect =< 314 degrees). This map shows a smoothed and generalized image of the four topographic moisture potential classes. Additional information about this map and any of the data developed for the ecosystems modeling of the conterminous United States is available online at http://rmgsc.cr.usgs.gov/ecosystems/.

Scientific Investigations Map↗

Simulation of the effects of nearby quarrying operations on ground-water flow at the South Well Field, Franklin County, Ohio

The City of Columbus, Ohio, operates a municipal well field in southern Franklin County that is adjacent to a sand and gravel mining operation. Mining operations have the potential to alter ground-water flowpaths and change the sources of water to pumped wells. Previous ground-water-flow modeling of the area has shown that water pumped from the supply wells is derived from infiltration from nearby rivers and surrounding bedrock. Some of that water flows through existing quarries. Because water quality differs among these sources and is affected by the path along which water flows to the wells, five flow conditions were simulated to evaluate the influence of different mining scenarios on sources of water as related to the size and shape of contributing recharge areas (CRAs) to wells. The first simulation was based on a revision of an existing model by Schalk (1996). The second and third simulations included one in which a 20-foot layer of undisturbed aquifer material within the quarry above the bedrock is left intact, and another in which the 20-foot layer is removed. The fourth and fifth simulations included one in which the 20-foot layer of undisturbed aquifer material is left above the bedrock and the quarry is backfilled with fine- grained sand and silt (a byproduct of the mining operations), and another in which the 20-foot layer is removed before the quarry is backfilled with the fine-grained sand and silt. The results of the five model simulations indicate that the overall volumetric budgets among models change only slightly in response to changing conditions at the quarry. The most significant change is noted in the amount of water that the aquifers gained from constant head and river leakage. This change is due to the way the quarries were simulated and lower heads in the aquifers compared to those in simulations made with earlier models. Previously published model simulations showed that the 5-year CRAs did not extend into the area of the newest sand and gravel quarry. In this study, however, the size and shape of the CRAs differ in response to the addition of two supply wells, and the 5-year CRA for one of these wells intersects the sand and gravel quarry. Particle-tracking analysis indicates that the proportions of water from the principal water sources -- rivers, other quarries (excluding the quarry of interest), and bedrock -- stayed relatively constant among the model simulations. The number of particles that originate in the quarry of interest increased from about 1 percent of the total particles withdrawn through all wells when the 20-foot layer of undisturbed aquifer material was left intact, to 2 percent when the 20-foot layer was removed. In simulations in which the quarry of interest was filled with fine-grained sand and silt, the percentage of particles that travel through or originate in that quarry decreased to less than 1 percent of the total number of particles. Thus, the mining activities at the quarry of interest, combined with increased pumping at the two supply wells have only a small potential to change the overall water quality of water withdrawn from supply wells at the South Well Field.

Water-Resources Investigations Report↗

Simulation of the Ground-Water Flow System in 1992, and Simulated Effects of Projected Ground-Water Withdrawals in 2020 in the New Jersey Coastal Plain

In 1992, ground-water withdrawals from the unconfined and confined aquifers in the New Jersey Coastal Plain totaled about 300 million gallons per day, and about 70 percent (200 million galllons per day) of this water was pumped from confined aquifers. The withdrawals have created large cones of depression in several Coastal Plain aquifers near populated areas, particularly in Camden and Ocean Counties. The continued decline of water levels in confined aquifers could cause saltwater intrusion, reduction of stream discharge near the outcrop areas of these aquifers, and depletion of the ground-water supply. Because of this, withdrawals from wells located within these critical areas have been reduced in the Potomac-Raritan-Magothy aquifer system, the Englishtown aquifer system, and the Wenonah-Mount Laurel aquifer. A computer-based model that simulates freshwater and saltwater flow was used to simulate transient ground-water flow conditions and the location of the freshwater-saltwater interface during 1989-92 in the New Jersey Coastal Plain. This simulation was used as the baseline for comparison of water levels and flow budgets. Four hypothetical withdrawal scenarios were simulated in which ground-water withdrawals were either increased or decreased. In scenario 1, withdrawals from wells located within critical area 2 in the Potomac-Raritan-Magothy aquifer system were reduced by amounts ranging from 0 to 35 percent of withdrawals prior to 1992. Critical area 2 is mainly located in Camden County, and most of Burlington and Gloucester Counties. With the reductions, water levels recovered about 30 feet in the regional cone of depression centered in Camden County in the Upper Potomac-Raritan-Magothy aquifer and by 20 ft in the Lower and Middle Potomac-Raritan-Magothy aquifers. In scenarios 2 to 4, withdrawals projected for 2020 were input to the model. In scenario 2, withdrawal restrictions within the critical areas were imposed in the Potomac-Raritan-Magothy aquifer system, the Englishtown aquifer system, and the Wenonah-Mount Laurel aquifer, but withdrawals were increased outside the critical areas to the projected 2020 demand. With withdrawals restrictions in critical areas, water levels recovered about 20 feet at the center of the regional cone of depression in the Upper Potomac-Raritan Magothy aquifer. Water levels recovered by about 20 feet at the center of a regional cone of depression in the Englishtown aquifer system in Ocean County, and by about 20 feet in the Wenonah-Mount Laurel aquifer in the same area. In scenario 3, withdrawals were increased to the projected 2020 demand inside and outside the critical areas. As a result, water levels declined as much as 20 feet at the center of a regional cone of depression in the Englishtown aquifer system in Ocean County, and as much as 10 feet in the Wenonah-Mounty Laurel aquifer near this area. The Englishtown aquifer system and the Wenonah-Mount Laurel aquifer are particularly sensitive to increases and decreases in withdrawals because in certain areas the transmissivities of these aquifers are lower than the transmissivities of other confined aquifers of the New Jersey Coastal Plain, and because these aquifers are hydraulically connected. Simulated water levels declined by as much as 10 ft at the center of the regional cone of depression in Atlantic County. In scenario 4, withdrawal amounts were equal to that in scenario 2, except an additional 13.2 million gallons per day was withdrawn from hypothetical wells located outside the critical areas in the Upper Potomac-Raritan-Magothy aquifer, Englishtown aquifer system, and the Wenonah-Mount Laurel aquifer. The additional withdrawals resulted in increased leakage from overlying aquifers to the Wenonah-Mount Laurel aquifer and subsequently to the Englishtown aquifer system.

Water-Resources Investigations Report↗

Hydrogeology and steady-state numerical simulation of groundwater flow in the Lost Creek Designated Ground Water Basin, Weld, Adams, and Arapahoe Counties, Colorado

The Lost Creek Designated Ground Water Basin (Lost Creek basin) is an important alluvial aquifer for irrigation, public supply, and domestic water uses in northeastern Colorado. Beginning in 2005, the U.S. Geological Survey, in cooperation with the Lost Creek Ground Water Management District and the Colorado Water Conservation Board, collected hydrologic data and constructed a steady-state numerical groundwater flow model of the Lost Creek basin. The model builds upon the work of previous investigators to provide an updated tool for simulating the potential effects of various hydrologic stresses on groundwater flow and evaluating possible aquifer-management strategies. As part of model development, the thickness and extent of regolith sediments in the basin were mapped, and data were collected concerning aquifer recharge beneath native grassland, nonirrigated agricultural fields, irrigated agricultural fields, and ephemeral stream channels. The thickness and extent of regolith in the Lost Creek basin indicate the presence of a 2- to 7-mile-wide buried paleovalley that extends along the Lost Creek basin from south to north, where it joins the alluvial valley of the South Platte River valley. Regolith that fills the paleovalley is as much as about 190 ft thick. Average annual recharge from infiltration of precipitation on native grassland and nonirrigated agricultural fields was estimated by using the chloride mass-balance method to range from 0.1 to 0.6 inch, which represents about 1-4 percent of long-term average precipitation. Average annual recharge from infiltration of ephemeral streamflow was estimated by using apparent downward velocities of chloride peaks to range from 5.7 to 8.2 inches. Average annual recharge beneath irrigated agricultural fields was estimated by using passive-wick lysimeters and a water-balance approach to range from 0 to 11.3 inches, depending on irrigation method, soil type, crop type, and the net quantity of irrigation water applied. Estimated average annual recharge beneath irrigated agricultural fields represents about 0-43 percent of net irrigation. The U.S. Geological Survey modular groundwater modeling program, MODFLOW-2000, was used to develop a steady-state groundwater flow model of the Lost Creek basin. Groundwater in the basin is simulated generally to flow from the basin margins toward the center of the basin and northward along the paleovalley. The largest source of inflow to the model occurs from recharge beneath flood- and sprinkler-irrigated agricultural fields (14,510 acre-feet per year [acre-ft/yr]), which represents 39.7 percent of total simulated inflow. Other substantial sources of inflow to the model are recharge from precipitation and stream-channel infiltration in nonirrigated areas (13,810 acre-ft/yr) seepage from Olds Reservoir (4,280 acre-ft/yr), and subsurface inflow from ditches and irrigated fields outside the model domain (2,490 acre-ft/yr), which contribute 37.7, 11.7, and 6.8 percent, respectively, of total inflow. The largest outflow from the model occurs from irrigation well withdrawals (26,760 acre-ft/yr), which represent 73.2 percent of total outflow. Groundwater discharge (6,640 acre-ft/yr) at the downgradient end of the Lost Creek basin represents 18.2 percent of total outflow, and evapotranspiration (3,140 acre-ft/yr) represents about 8.6 percent of total outflow.

Scientific Investigations Report↗

A precipitation-runoff model for the analysis of the effects of water withdrawals and land-use change on streamflow in the Usquepaug–Queen River Basin, Rhode Island

The 36.1-square-mile Usquepaug–Queen River Basin in south-central Rhode Island is an important water resource. Streamflow records indicate that withdrawals may have diminished flows enough to affect aquatic habitat. Concern over the effect of withdrawals on streamflow and aquatic habitat prompted the development of a Hydrologic Simulation Program–FORTRAN (HSPF) model to evaluate the water-management alternatives and land-use change in the basin. Climate, streamflow, and water-use data were collected to support the model development. A logistic-regression equation was developed for long-term simulations to predict the likelihood of irrigation, the primary water use in the basin, from antecedent potential evapotranspiration and precipitation for generating irrigation demands. The HSPF model represented the basin by 13 pervious-area and 2 impervious-area land-use segments and 20 stream reaches. The model was calibrated to the period January 1, 2000 to September 30, 2001, at three continuous streamflow-gaging stations that monitor flow from 10, 54, and 100 percent of the basin drainage area. Hydrographs and flow-duration curves of observed and simulated discharges, along with statistics compiled for various model-fit metrics, indicate a satisfactory model performance. The calibrated HSPF model was modified to evaluate streamflow (1) under no withdrawals to streamflow under current (2000–01) withdrawal conditions under long-term (1960–2001) climatic conditions, (2) under withdrawals by the former Ladd School water-supply wells, and (3) under fully developed land use. The effects of converting from direct-stream withdrawals to ground-water withdrawals were evaluated outside of the HSPF model by use of the STRMDEPL program, which calculates the time delayed response of ground-water withdrawals on streamflow depletion. Simulated effects of current withdrawals relative to no withdrawals indicate about a 20-percent decrease in the lowest mean daily streamflows at the basin outlet, but withdrawals have little effect on flows that are exceeded less than about 90 percent of the time. Tests of alternative model structures to evaluate model uncertainty indicate that the lowest mean daily flows ranged between 3 and 5 cubic feet per second (ft3/s) without withdrawals and 2.2 to 4 ft3/s with withdrawals. Changes in the minimum daily streamflows are more pronounced, however; at the upstream streamflow-gaging station, a minimum daily flow of 0.2 ft3/s was sustained without withdrawals, but simulations with withdrawals indicate that the reach would stop flowing part of a day about 5 percent of the time. The effect on streamflow of potential ground-water withdrawals of 0.20, 0.90, and 1.78 million gallons per day (Mgal/d) at the former Ladd School near the central part of the basin were evaluated. The lowest daily mean flows in model reach 3, the main stem of the Queen River closest to the pumped wells, decreased by about 50 percent for withdrawals of 0.20 Mgal/d (from about 0.4 to 0.2 ft3/s) in comparison to current withdrawals. Reach 3 would occasionally stop flowing during part of the day at the 0.20-Mgal/d withdrawal rate because of diurnal fluctuation in streamflow. The higher withdrawal rates (0.90 and 1.78 Mgal/d) would cause reach 3 to stop flowing about 10 to 20 percent of the time, but the effects of pumping rapidly diminished downstream because of tributary inflows. Simulation results indicate little change in the annual 1-, 7-, and 30-day low flows at the 0.20 Mgal/d pumping rate, but at the 1.78 Mgal/d pumping rate, reach 3 stopped flowing for nearly a 7-day period every year and for a 30-day period about every other year. At the 0.90 Mgal/d pumping rate, reach 3 stopped flowing about every other year for a 7-day period and about once every 5 years for a 30-day period. Land-use change was simulated by converting model hydrologic-response units (HRUs) representing undeveloped areas to HRUs representing developed areas on the basis of development suitability and town zoning. About 55 percent of the basin is suitable for development; this area would accommodate about 4,300 new low-density residential homes under current zoning. Increases in storm volume and peak flows, and decreases in base flow, typically associated with urbanization, were not evident in buildout simulations because the effective impervious area was assumed to increase by only 2 percent. Under fully developed conditions, withdrawals from self-supply wells were estimated to reach 1.2 Mgal/d. Potential increases in water withdrawals for a fully developed basin have only a minor impact on the main stem streamflow, but the effects of urbanization could be more pronounced in localized areas where development is concentrated. Streamflow-depletion rates were calculated for varying distances of a pumped irrigation well from a stream. For the irrigation rates and aquifer conditions tested, streamflow depletion, relative to the pumping rate, decreases rapidly as the pumped well was moved away from the stream. Streamflow depletion, relative to the peak withdrawal rate, decreased by about 60, 80, and 90 percent by locating the pumped well 500, 1,000, and 1,500 feet from the stream, respectively.

Rhode Island↗

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↗

Influence of Locally Derived Recharge on the Water Quality and Temperature of Springs in Hot Springs National Park, Arkansas

The hot springs of Hot Springs National Park consist of a mixture of water from two recharge components: a primary hot-water component and a secondary cold-water component. Widespread distribution of fractures enables mixing of the hot- and cold-water components of flow near the discharge area for the springs. Urbanization in the area near the hot springs of Hot Springs National Park has increased the potential for degradation of the quality of surface-water runoff and locally derived ground-water recharge to the hot springs. Previous studies by the U.S. Geological Survey have indicated that water from some cold-water springs and wells in the vicinity of Hot Springs, Arkansas, showed evidence of contamination and that water from locally derived cold-water recharge might contribute 25 percent of the total flow to the hot springs after storms. Water samples were collected during base-flow conditions at nine hot springs and two cold-water springs in September 2000. Nine hot springs and one cold-water spring were resampled in October 2001 after a storm that resulted in a measurable decrease in water temperature in selected hot springs. Water samples were analyzed for a variety of dissolved chemical constituents (nutrients, major ions, trace elements, pesticides, semivolatile compounds, isotopes, and radiochemicals), physical properties, field measurements, and bacteria. Comparison of analyses of samples collected during base-flow conditions from the springs in 2000 and during a storm event in 2001 with the results from earlier studies dating back to the late 1800's indicates that little change in major, minor, and trace constituent chemistry has occurred and that the water continues to be of excellent quality. Water-quality data show distinguishable differences in water chemistry of the springs during base-flow and stormflow conditions, indicating changing input of cold-water recharge relative to hot-water recharge. Silica, total dissolved solids, strontium, barium, and sulfate show statistically significant differences between the median values of base-flow and stormflow samples. While variations in these constituents do not degrade water quality, the differences do provide evidence of variability in the factors controlling water quality of the hot springs and show that water quality is influenced by the locally derived, cold-water component of flow to the springs. Water temperature was measured continuously (3-minute intervals) between August 2000 and October 2002 at four hot springs. Continuous water-temperature data at the springs provide no indication of persistent long-term change in water temperature through time. Short time-scale water-temperature decreases occur in response to mixing of hot-springs water with locally derived recharge after storm events; the magnitude of these decreases varied inversely with the amount of rainfall. Maximum decreases in water temperature for specific storms had a non-linear relation with the amount of precipitation measured for the events. Response time for water temperature to begin decreasing from baseline temperature as a result of storm recharge was highly variable. Some springs began decreasing from baseline temperature as quickly as 1 hour after the beginning of a storm; one spring had an 8-hour minimum response time to show a storm-related temperature decrease. Water-quality, water-temperature, isotopic, and radiochemical data provide multiple lines of evidence supporting the importance of the contribution of cold-water recharge to hot springs. All the springs sampled indicated some measure of influence from local recharge. Binary mixing models using silica and total dissolved solids indicate that cold-water recharge from stormflow contributes an estimated 10 to 31 percent of the flow of hot springs. Models using water temperature indicate that cold-water recharge from stormflow contributes an estimated 1 to 35 percent of the flow of the various hot springs. Alth

Scientific Investigations Report↗