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 739 records · Page 41Linked to original sources

Assessment of Eutrophication in the Lower Yakima River Basin, Washington, 2004-07

In response to concerns that excessive plant growth in the lower Yakima River in south-central Washington was degrading water quality and affecting recreational use, the U.S. Geological Survey and the South Yakima Conservation District conducted an assessment of eutrophication in the lower 116 miles of the river during the 2004-07 irrigation seasons (March - October). The lower Yakima River was divided into three distinct reaches based on geomorphology, habitat, aquatic plant and water-quality conditions. The Zillah reach extended from the upstream edge of the study area at river mile (RM) 116 to RM 72, and had abundant periphyton growth and sparse macrophyte growth, the lowest nutrient concentrations, and moderately severe summer dissolved oxygen and pH conditions in 2005. The Mabton reach extended from RM 72 to RM 47, and had sparse periphyton and macrophyte growth, the highest nutrient conditions, but the least severe summer dissolved oxygen and pH conditions in 2005. The Kiona reach extended from RM 47 to RM 4, and had abundant macrophyte and epiphytic algae growth, relatively high nutrient concentrations, and the most severe summer dissolved oxygen and pH conditions in 2005. Nutrient concentrations in the lower Yakima River were high enough at certain times and locations during the irrigation seasons during 2004-07 to support the abundant growth of periphytic algae and macrophytes. The metabolism associated with this aquatic plant growth caused large daily fluctuations in dissolved oxygen concentrations and pH levels that exceeded the Washington State water-quality standards for these parameters between July and September during all 4 years, but also during other months when streamflow was unusually low. The daily minimum dissolved oxygen concentration was strongly and negatively related to the preceding day's maximum water temperature - information that could prove useful if a dissolved oxygen predictive model is developed for the lower Yakima River. Periphytic algal growth generally was not nutrient-limited and frequently reached nuisance levels in the Zillah reach, where some surface-water nutrient concentrations were below the reference concentrations suggested by the U.S. Environmental Protection Agency. Although lowering nutrient concentrations in this reach might limit periphytic algal growth enough to improve dissolved oxygen and pH conditions, ground water inflow at some locations might still provide an adequate supply of nutrients for periphytic algal growth. Macrophyte growth in the Kiona reach was dominated by water stargrass (Heteranthera dubia), was far greater compared to the other two reaches, varied greatly between years, and was negatively related to greater spring runoff due to lower light availability. Lowering nutrient concentrations in the Kiona reach might not impact the level of macrophyte growth because macrophytes with extensive root systems such as water stargrass can get nutrients from river sediment. In addition, the results from this study did not indicate any nutrient uptake by the macrophytes from the water column (nutrient uptake from the sediment was not examined). Creating the prolonged turbid and deep conditions during spring necessary to suppress macrophyte growth in this reach would not be possible in years with low streamflow. In addition, because of the relatively stable substrate present in much of this reach, the macrophyte root systems would likely not be disturbed under all but the most extremely high streamflows that occur in the lower Yakima River.

Scientific Investigations Report↗

A method to estimate canal leakage to the Biscayne Aquifer, Dade County, Florida

The leakage characteristics of channels that partially penetrate the Biscayne aquifer and have reduced bed permeability were studied. Leakage characteristics were described in terms of a reach transmissivity-defined as the volume flow rate out of the channel per unit length of the channel per unit drawdown, where drawdown is defined as the difference in altitude between the water surface in the canal and the water table in the adjacent aquifer. A theoretical expression was developed to relate the reach transmissivity to the transmissivity of the formation, mean channel width, distance of drawdown measurement from the channel centerline, ratio of drawdowns on both sides of the channel, and local reach transmissivity associated with reduced bed permeability. This theoretical expression was verified using a fine-scale numerical model, which gave accurate results when drawdowns were measured beyond 10 aquifer depths from the side of the channel. Using the theoretical formulation, it is shown that the reach transmissivity employed in regional ground-water models, which are based on average drawdowns within a cell, depends on the size of the cell as well as the transmissivity of the formation, channel width, and local reach transmissivity due to reduced bed permeability. The theoretical reach transmissivity function was compared with field measurements at L-31N Canal and Snapper Creek Extension Canal in Dade County, Florida. Analyses of the data for both canals showed good agreement between the estimated and measured reach transmissivities. At L- 31N Canal, field measurements indicated that the local reach transmissivity was relatively uniform over a 2-mile reach of the channel (averaging 630 cubic feet per second per mile per foot), and the formation transmissivity was 1.8 x106 feet squared per day. At Snapper Creek Extension Canal, an approximate analysis was necessary due to the inability of the acoustic velocity meter to measure very low water velocities in the channel. Assuming an aquifer transmissivity of 1 x 106 feet squared per day, drawdown measurements indicated that the local reach transmissivity was about 400 cubic feet per second per mile per foot. The theoretical relation, combined with the local reach transmissivity and formation transmissivity, was sufficient to predict the leakage out of L-31N Canal and Snapper Creek Extension Canal for any drawdown scenario.

Florida↗

Areal studies aid protection of ground-water quality in Illinois, Indiana, and Wisconsin

In 1991, the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, initiated studies designed to characterize the ground-water quality and hydrogeology in northern Illinois, and southern and eastern Wisconsin (with a focus on the north-central Illinois cities of Belvidere and Rockford, and the Calumet region of northeastern Illinois and northwestern Indiana). These areas are considered especially susceptible to ground-water contamination because of the high density of industrial and waste-disposal sites and the shallow depth to the unconsolidated sand and gravel aquifers and the fractured, carbonate bedrock aquifers that underlie the areas. The data and conceptual models of ground-water flow and contaminant distribution and movement developed as part of the studies have allowed Federal, State, and local agencies to better manage, protect, and restore the water supplies of the areas. Water-quality, hydrologic, geologic, and geophysical data collected as part of these areal studies indicate that industrial contaminants are present locally in the aquifers underlying the areas. Most of the contaminants, particularly those at concentrations that exceeded regulatory water-quality levels, were detected in the sand and gravel aquifers near industrial or waste-disposal sites. In water from water-supply wells, the contaminants that were present generally were at concentrations below regulatory levels. The organic compounds detected most frequently at concentrations near or above regulatory levels varied by area. Trichloroethene, tetrachloroethene, and 1,1,1-trichloroethane (volatile chlorinated compounds) were most prevalent in north-central Illinois; benzene (a petroleum-related compound) was most prevalent in the Calumet region. Differences in the type of organic compounds that were detected in each area likely reflect differences in the types of industrial sites that predominate in the areas. Nickel and aluminum were the trace metals detected most frequently at concentrations above regulatory levels in both areas. Contaminants in the shallow sand and gravel aquifers and carbonate aquifers appear to have moved with ground water discharging to local lakes, streams, and wetlands. Ground-water flow and possibly contaminant movement is concentrated in the weathered surface zones and in deeper fractures of the carbonate aquifers underlying both areas.

Illinois, Indiana, Wisconsin↗

Potential for saltwater intrusion into the lower Tamiami aquifer near Bonita Springs, southwestern Florida

A study was conducted to examine the potential for saltwater intrusion into the lower Tamiami aquifer beneath Bonita Springs in southwestern Florida. Field data were collected, and constant- and variable-density ground-water flow simulations were performed that: (1) spatially quantified modern and seasonal stresses, (2) identified potential mechanisms of saltwater intrusion, and (3) estimated the potential extent of saltwater intrusion for the area of concern. MODFLOW and the inverse modeling routine UCODE were used to spatially quantify modern and seasonal stresses by calibrating a constant-density ground-water flow model to field data collected in 1996. The model was calibrated by assuming hydraulic conductivity parameters were accurate and by estimating unmonitored ground-water pumpage and potential evapotranspiration with UCODE. Uncertainty in these estimated parameters was quantified with 95-percent confidence intervals. These confidence intervals indicate more uncertainty (or less reliability) in the estimates of unmonitored ground-water pumpage than estimates of pan-evaporation multipliers, because of the nature and distribution of observations used during calibration. Comparison of simulated water levels, streamflows, and net recharge with field data suggests the model is a good representation of field conditions. Potential mechanisms of saltwater intrusion into the lower Tamiami aquifer include: (1) lateral inland movement of the freshwater-saltwater interface from the southwestern coast of Florida; (2) upward leakage from deeper saline water-bearing zones through natural upwelling and upconing, both of which could occur as diffuse upward flow through semiconfining layers, conduit flow through karst features, or pipe flow through leaky artesian wells; (3) downward leakage of saltwater from surface-water channels; and (4) movement of unflushed pockets of relict seawater. Of the many potential mechanisms of saltwater intrusion, field data and variable-density ground-water flow simulations suggest that upconing is of utmost concern, and lateral encroachment is of second-most concern. This interpretation is uncertain, however, because the predominance of saltwater intrusion through leaky artesian wells with connection to deeper, more saline, and higher pressure aquifers was difficult to establish. Effective management of ground-water resources in southwestern Florida requires an understanding of the potential extent of saltwater intrusion in the lower Tamiami aquifer near Bonita Springs. Variable-density, ground-water flow simulations suggest that when saltwater is at dynamic equilibrium with 1996 seasonal stresses, the extent of saltwater intrusion is about 100 square kilometers areally and 70,000 hectare-meters volumetrically. The volumetric extent of saltwater intrusion was most sensitive to changes in recharge, ground-water pumpage, sea level, salinity of the Gulf of Mexico, and the potentiometric surface of the sandstone aquifer, respectively.

Florida↗

Water levels and ground-water discharge, regional aquifer system of the midwestern Basins and Arches Region, in parts of Indiana, Ohio, Illinois, and Michigan

Aquifers in Quaternary glacial deposits and the underlying Silurian and Devonian carbonate bedrock in parts of Indiana, Ohio, Illinois, and Michigan compose the regional aquifer system under investigation as part of the Midwestern Basins and Arches Regional Aquifer System Analysis (Midwestern Basins and Arches—RASA) project of the U.S. Geological Survey (USGS). The Midwestern Basins and Arches—RASA is part of a USGS program to assess the regional hydrology, geology, and water quality of the Nation's most important aquifers (Sun, 1986). An objective specific to the Midwestern Basins and Arches—RASA project is to conceptualize and describe regional ground-water flow in the glacial-deposit and carbonate-bedrock aquifer system, including regional recharge and discharge areas and regional relations between surface and ground water (Bugliosi, 1990). Water-level and ground-water discharge data were collected and (or) analyzed to help meet the above objective. Specifically, data from the USGS Ground-Water Site Inventory (GWSI) data base were used to determine relations between land-surface altitude and water levels in glacial-deposit aquifers. Water levels in the carbonate-bedrock aquifer were synoptically measured during July 1990, and the data were used to construct a potentiometric surface map of the aquifer. Regional hydraulic gradients and general directions of regional flow in the carbonate-bedrock aquifer can be inferred from this map. Steady-state groundwater discharge to streams that drain the area underlain by the glacial-deposit and carbonate bedrock aquifer system was estimated from base-flow daily values computed from streamflow records. Water-level and ground-water-discharge data collectively form the sample information necessary to develop calibration targets for calibration of a ground-water-flow model (Anderson and Woessner, 1992). Such a ground-water-flow model of the glacial-deposit and carbonate-bedrock aquifer system was constructed to help conceptualize and describe regional ground-water flow in the aquifer system. The model was calibrated to the water-level and ground-water-discharge data presented in this atlas.

Hydrologic Atlas↗

Ground-water availability and water quality at Southbury and Woodbury, Connecticut

Increases in population and commercial and industrial development during the past 20 years have increased the demand for water in the Towns of Southbury and Woodbury, Connecticut. The stratified-drift aquifer, underlying much of the Pomperaug River valley, is the most practical source for additional large supplies. The yield of the aquifer was evaluated with a two-dimensional, digital flow model. The model was constructed with hydrologic data from previous studies, and test boring logs , seismic profiles, water-level measurements, and other information collected during the present study. Simulations made with the calibrated model indicate that, with no pumpage, groundwater levels in the aquifer will fall about 4.6 ft below average during low-recharge (least-favorable) periods, and rise about 0.6 ft above average during high-recharge (most-favorable) periods. Simulated withdrawals from 10 wells indicate that from 5.0 to 8.8 million gallons/day are available as total recharge rates range from 21.4 to 36.1 inches/year. If these pumpages were consumed or exported from the basin, estimated average flow reductions of the Pomperaug River would range from 7.7 to 12.9 cu ft/sec. The quality of the water from the stratified-drift aquifer is generally excellent in most areas and meets State drinking-water standards. Chemical analyses of groundwater from 11 wells in the Middle Quarter area of Woodbury indicate that organohalide compounds are present. A maximum trichloroethane concentration of 260 micrograms/L has been reported and groundwater in the area is presently being monitored for organohalides. The water meets standards established by the State. Surface water samples collected at 7 sites in the study area meet the Connecticut drinking water standards for all constituents except coliform bacteria. Complete conventional treatment of surface water from some wells will be required to meet State drinking water standards relative to coliforms. (USGS)

Water-Resources Investigations Report↗

A science plan for a comprehensive assessment of water supply in the region underlain by fractured rock in Maryland

The fractured rock region of Maryland, which includes land areas north and west of the Interstate 95 corridor, is the source of water supply for approximately 4.4 million Marylanders, or approximately 76 percent of the State's population. Whereas hundreds of thousands of residents rely on wells (both domestic and community), millions rely on surface-water sources. In this region, land use, geology, topography, water withdrawals, impoundments, and other factors affect water-flow characteristics. The unconfined groundwater systems are closely interconnected with rivers and streams, and are affected by seasonal and climatic variations. During droughts, groundwater levels drop, thereby decreasing well yields, and in some cases, wells have gone dry. Low ground-water levels contribute to reduced streamflows, which in turn, can lead to reduced habitat for aquatic life. Increased demand, over-allocation, population growth, and climate change can affect the future sustainability of water supplies in the region of Maryland underlain by fractured rock. In response to recommendations of the 2008 Advisory Committee on the Management and Protection of the State's Water Resources report, the Maryland Department of the Environment's Water Supply Program, the Maryland Geological Survey, the Maryland Department of Natural Resources, Monitoring and Non-Tidal Assessment (MANTA) Division, and the U.S. Geological Survey have developed a science plan for a comprehensive assessment that will provide new scientific information, new data analysis, and new tools for the State to better manage water resources in the fractured rock region of Maryland. The science plan lays out five goals for the comprehensive assessment: (1) develop tools for the improved management and investigation of groundwater and surface-water resources; (2) characterize factors affecting reliable yields of individual groundwater and surface-water supplies; (3) investigate impacts on nearby water withdrawal users caused by groundwater and surface-water withdrawals; (4) assess the role of streamflow and water withdrawals on the ecological integrity of streams; and (5) improve understanding of the distribution of water-quality conditions in fractured rock aquifers. To accomplish these goals, accurate data collection, review, and analysis are needed, including the study of "Research Watersheds" that can provide detailed information about the potential effects that climate change and water withdrawals may have on groundwater, streamflow, and aquatic life. The assessment planning started in 2009 and is being conducted with close interagency coordination. A Fractured Rock Aquifer Information System is currently (2012) undergoing initial development. Other major tasks that will be performed include the development of work plans for each science goal, the estimation of daily streamflow at ungaged streams, and the design and implementation of Research Watersheds. Finally, scenarios will be modeled to evaluate current water allocation permitting methodologies, investigate effects on nearby water withdrawal users caused by groundwater and surface-water withdrawals, and assess the potential impacts of climate change on water resources. Desktop and Web-based tools will be developed in order to meet the diverse research needs of the assessment. These tools, including the Fractured Rock Aquifer Information System will be continuously improved during the assessment to store relevant groundwater and surface-water data in spatially referenced databases, estimate streamflows, locate higher-yielding wells, estimate the impacts of withdrawals on nearby users, and assess the cumulative impacts of withdrawals on the aquatic resource. Tools will be developed to serve the needs of many audiences, including water resource managers, water suppliers, planners, policymakers, and other scientific investigators.

Maryl↗

A hydrogeologic approach to identify land uses that overlie ground-water flow paths, Broward County, Florida

A hydrogeologic approach that integrates the use of hydrogeologic and spatial tools aids in the identification of land uses that overlie ground- water flow paths and permits a better understanding of ground-water flow systems. A mathematical model was used to simulate the ground-water flow system in Broward County, particle-tracking software was used to determine flow paths leading to the monitor wells in Broward County, and a Geographic Information System was used to identify which land uses overlie the flow paths. A procedure using a geographic information system to evaluate the output from a ground-water flow model has been documented. The ground-water flow model was used to represent steady-state conditions during selected wet- and dry-season months, and an advective flow particle- tracking program was used to simulate the direction of ground-water flow in the aquifer system. Digital spatial data layers were created from the particle pathlines that lead to the vicinity of the open interval of selected wells in the Broward County ground-water quality monitoring network. Buffer zone data layers were created, surrounding the particle pathlines to represent the area of contribution to the water sampled from the monitor wells. Spatial data layers, combined with a land-use data layer, were used to identify the land uses that overlie the ground-water flow paths leading to the monitor wells. The simulation analysis was performed on five Broward County wells with different hydraulic parameters to determine the source of ground-water stress, determine selected particle pathlines, and identify land use in buffer zones in the vicinity of the wells. The flow paths that lead to the grid cells containing wells G-2355, G-2373, and G-2373A did not vary between the wet- and dry-season conditions. Changes in the area of contribution for wells G-2345X and G-2369 were attributed to variations in rainfall patterns, well-field pumpage, and surface-water management practices. Additionally, using a different open interval at a site, such as for wells G-2373 and G-2373A, can result in a very different area that overlies the flow path leading to the monitor well.

Water-Resources Investigations Report↗

Development of a geodatabase and conceptual model of the hydrogeologic units beneath air force plant 4 and Naval Air Station-Joint Reserve Base Carswell Field, Fort Worth, Texas

Air Force Plant 4 and adjacent Naval Air Station-Joint Reserve Base Carswell Field at Fort Worth, Texas, constitute a government-owned, contractor-operated facility that has been in operation since 1942. Contaminants from AFP4, primarily volatile organic compounds and metals, have entered the ground-water-flow system through leakage from waste-disposal sites and from manufacturing processes. The U.S. Geological Survey developed a comprehensive geodatabase of temporal and spatial environmental information associated with the hydrogeologic units (alluvial aquifer, Goodland-Walnut confining unit, and Paluxy aquifer) beneath the facility and a three-dimensional conceptual model of the hydrogeologic units integrally linked to the geodatabase. The geodatabase design uses a thematic layer approach to create layers of feature data using a geographic information system. The various features are separated into relational tables in the geodatabase on the basis of how they interact and correspond to one another. Using the geodatabase, geographic data at the site are manipulated to produce maps, allow interactive queries, and perform spatial analyses. The conceptual model for the study area comprises computer-generated, three-dimensional block diagrams of the hydrogeologic units. The conceptual model provides a platform for visualization of hydrogeologic-unit sections and surfaces and for subsurface environmental analyses. The conceptual model is based on three structural surfaces and two thickness configurations of the study area. The three structural surfaces depict the altitudes of the tops of the three hydrogeologic units. The two thickness configurations are those of the alluvial aquifer and the Goodland-Walnut confining unit. The surface of the alluvial aquifer was created using a U.S. Geological Survey 10-meter digital elevation model. The 2,130 point altitudes of the top of the Goodland-Walnut unit were compiled from lithologic logs from existing wells, available soil-boring logs, and previous studies. Data from 120 wells, primarily from existing reports, were used to create a map of the approximate altitude of the Paluxy aquifer.

Texas↗

Ground-water contamination at Wurtsmith Air Force Base, Michigan

A sand and gravel aquifer of glacial origin underlies Wurtsmith Air Force Base in northeastern lower Michigan. The aquifer overlies a thick clay layer at an average depth of 65 feet. The water table is about 10 feet below land surface in the western part of the Base and about 25 feet below land surface in the eastern part. A ground-water divide cuts diagonally across the Base from northwest to southeast. South of the divide, ground water flows to the Au Sable River; north of the divide, it flows to Van Etten Creek and Van Etten Lake. Mathematical models were used to aid in calculating rates of groundwater flow. Rates range from about 0.8 feet per day in the eastern part of the Base to about 0.3 feet per day in the western part. Models also were used as an aid in making decisions regarding purging of contaminated water from the aquifer. In 1977, trichloroethylene was detected in the Air Force Base water-supply system. It had leaked from a buried storage tank near Building 43 in the southeastern part of the Base and moved northeastward under the influence of the natural ground-water gradient and the pumping of Base water-supply wells. In the most highly contaminated part of the plume, concentrations are greater than 1,000 micrograms per liter. Current purge pumping is removing some of the trichloroethylene, and seems to have arrested its eastward movement. Pumping of additional purge wells could increase the rate of removal. Trichloroethylene has also been detected in ground water in the vicinity of the Base alert apron, where a plume from an unknown source extends northeastward off Base. A smaller, less well-defined area of contamination also occurs just north of the larger plume. Trichloroethylene, identified near the waste-treatment plant, seepage lagoons, and the northern landfill area, is related to activities and operations in these areas. Dichloroethylene and trichloroethylene occur in significant quantities westward of Building 43, upgradient from the major spill site. Benzene, indicative of ground-water contamination by a fuel substance, occurs in an area northeast of the bulk-fuel storage area. Analysis of a variety of chemical, physical, and biologic characteristics of water on the Base indicate that there is a measurable affect on ground-water quality from landfills, the seepage lagoon, and the waste-treatment plant.

Michigan↗

Ground-Water Resource Assessment in the Rio Grande de Manati Alluvial Plain, Rio Arriba Saliente Area, Puerto Rico

The alluvial aquifer within a 160-acre area of the Rio Grande de Manati alluvial plain was investigated to evaluate its potential as a water-supply source for the Barrios Rio Arriba Saliente and Pugnado Afuera, municipio of Manati, Puerto Rico. Analysis of well boring samples and the results of electric resistivity surveys indicate that the average thickness of the unconsolidated alluvial deposits in the study area is about 100 to 110 feet. The alluvium is a mixture of sand and gravel, which generally has a porosity of 0.2 to 0.35. Short-duration pump tests in small-diameter piezometers indicate that the alluvial aquifer has a hydraulic conductivity of about 200 feet per day and a transmissivity of about 7,900 feet squared per day. Analyses of water levels in piezometers, combined with stage measurements at a series of surveyed reference points along the Rio Grande de Manati channel, indicate that the water-table gradient in the alluvial aquifer is about 0.001, and that ground-water flow is generally from south to north, in the general direction of river flow. The water-table data indicate that the Rio Grande de Manati is the principal source of ground-water recharge to the alluvial aquifer in the study area. Because base flow for the Rio Grande de Manati is usually greater than 44 cubic feet per second, a continuous withdrawal rate of 0.5 to 1.0 cubic foot per second (225 to 450 gallons per minute) from a production well is possible. Chemical analysis of a ground-water sample indicates that the alluvial aquifer water meets U.S. Environmental Protection Agency secondary standards for selected constituents. Bacteriological analysis of ground-water samples indicates that the ground water contains little or no fecal coliform or fecal streptococcus bacteria. Although long-term data from upstream of the study area indicate high levels of fecal coliform and fecal streptococcus prior to 1996, bacteriological analyses of Rio Grande de Manati water samples obtained during the present study indicate that fecal coliform and fecal streptococcus concentrations are within the standards for surface water intended for use (or with the potential for use) as a raw source of public water supply in Puerto Rico. If a production well were constructed in the study area, it would be located close to the river channel (within 500 to 800 feet). Pumping from the porous and permeable alluvial aquifer close to the river channel could substantially enhance recharge from the Rio Grande de Manati channel to the aquifer. Enhanced recharge could shorten travel times for ground water in the aquifer, which might not allow sufficient time to attenuate bacteria and viruses. Travel times for bacteria moving from the river channel to a hypothetical production well were estimated using the numerical transport model MODFLOW/MT3DMS with an uncalibrated model of the alluvial aquifer. The model assumes a well pumping at 1 cubic foot per second. The transport of particles from the river to the well is most sensitive to the porosity of the aquifer and the pumping rate of the well. Sensitivity analysis indicates that a decrease in pumpage will increase the time of travel for particles to move from the river to the pumping well. The model indicates that the leading edge of a plume would reach the production well in about 40 days assuming a porosity of 0.20, 60 days assuming a porosity of 0.275, and about 70 days assuming a porosity of 0.35. If the well were moved 50 feet further from the river, the leading edge of the plume would reach the well in about 50 days assuming a porosity of 0.20 and about 70 days assuming a porosity of 0.275. These estimates are considered worse case estimates because no decay rate was included in the simulation, and because the hypothetical well was located in the center of the alluvial plain rather than further eastward, away from the river channel.

Water-Resources Investigations Report↗

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

The Mississippi River Valley alluvial aquifer is a water-bearing assemblage of gravels and sands that underlies about 32,000 square miles of Missouri, Kentucky, Tennessee, Mississippi, Louisiana, and Arkansas. Because of the heavy demands placed on the aquifer, several large cones of depression have formed in the potentiometric surface, resulting in lower well yields and degraded water quality in some areas. A ground-water flow model of the alluvial aquifer was previously developed for an area covering 3,826 square miles, extending south from the Arkansas River into the southeastern corner of Arkansas, parts of northeastern Louisiana, and western Mississippi. The flow-model results indicated that continued ground-water withdrawals at rates commensurate with those of 1997 could not be sustained indefinitely without causing water levels to decline below half the original saturated thickness of the aquifer. Conjunctive-use optimization modeling was applied to the flow model of the alluvial aquifer to develop withdrawal rates that could be sustained relative to the constraints of critical ground-water area designation. These withdrawal rates form the basis for estimates of sustainable yield from the alluvial aquifer and from rivers specified within the alluvial aquifer model. A management problem was formulated as one of maximizing the sustainable yield from all ground-water and surface-water withdrawal cells within limits imposed by plausible withdrawal rates, and within specified constraints involving hydraulic head and streamflow. Steady-state conditions were selected because the maximized withdrawals are intended to represent sustainable yield of the system (a rate that can be maintained indefinitely).One point along the Arkansas River and one point along Bayou Bartholomew were specified for obtaining surface-water sustainable-yield estimates within the optimization model. Streamflow constraints were specified at two river cells based on average 7-day low flows with 10-year recurrence intervals. Sustainable-yield estimates were affected by the allowable upper limit on withdrawals from wells specified in the optimization model. Withdrawal rates were allowed to increase to 200 percent of the withdrawal rate in 1997. As the overall upper limit is increased, the sustainable yield generally increases. Tests with the optimization model show that without limits on pumping, wells adjacent to sources of water, such as large rivers, would have optimal withdrawal rates that were orders of magnitude larger than rates corresponding to those of 1997. Specifying an upper withdrawal limit of 100 percent of the 1997 withdrawal rate, the sustainable yield from ground water for the entire study area is 70.3 million cubic feet per day, which is about 96 percent of the amount withdrawn in 1997 (73.5 million cubic feet per day). If the upper withdrawal limit is increased to 150 percent of the 1997 withdrawal rate, the sustainable yield from ground water for the entire study area is 80.6 million cubic feet per day, which is about 110 percent of the amount withdrawn in 1997. If the upper withdrawal limit is increased to 200 percent of the 1997 withdrawal rate, the sustainable yield from ground water for the entire study area is 110.2 million cubic feet per day, which is about 150 percent of the amount withdrawn in 1997. Total sustainable yield from the Arkansas River and Bayou Bartholomew is about 4,900 million cubic feet per day, or about 6,700 percent of the amount of ground-water withdrawn in 1997. The large, sustainable yields from surface water represent a potential source of water that could supplement ground water and meet the total water demand. Unmet demand (defined as the difference between the optimized withdrawal rate or sustainable yield, and the anticipated demand) was calculated using different demand rates based on multiples of the 1997-withdrawal rate. Assuming that demand is the 1997 withdrawal rate, and that sustainable-

Water-Resources Investigations Report↗

Recalibration of a ground-water flow model of the Mississippi River Valley alluvial aquifer of northeastern Arkansas, 1918-1998, with simulations of water levels caused by projected ground-water withdrawals through 2049

A digital model of the Mississippi River Valley alluvial aquifer in eastern Arkansas was used to simulate ground-water flow for the period from 1918 to 2049. The model results were used to evaluate effects on water levels caused by demand for ground water from the alluvial aquifer, which has increased steadily for the last 40 years. The model results showed that water currently (1998) is being withdrawn from the aquifer at rates greater than what can be sustained for the long term. The saturated thickness of the alluvial aquifer has been reduced in some areas resulting in dry wells, degraded water quality, decreased water availability, increased pumping costs, and lower well yields. The model simulated the aquifer from a line just north of the Arkansas-Missouri border to south of the Arkansas River and on the east from the Mississippi River westward to the less permeable geologic units of Paleozoic age. The model consists of 2 layers, a grid of 184 rows by 156 columns, and comprises 14,118 active cells each measuring 1 mile on a side. It simulates time periods from 1918 to 1998 along with further time periods to 2049 testing different pumping scenarios. Model flux boundary conditions were specified for rivers, general head boundaries along parts of the western side of the model and parts of Crowleys Ridge, and a specified head boundary across the aquifer further north in Missouri. Model calibration was conducted for observed water levels for the years 1972, 1982, 1992, and 1998. The average absolute residual was 4.69 feet and the root-mean square error was 6.04 feet for the hydraulic head observations for 1998. Hydraulic-conductivity values obtained during the calibration process were 230 feet per day for the upper layer and ranged from 230 to 730 feet per day for the lower layer with the maximum mean for the combined aquifer of 480 feet per day. Specific yield values were 0.30 throughout the model and specific storage values were 0.000001 inverse-feet throughout the model. Areally specified recharge rates ranged from 0 to about 30 inches and total recharge increased from 1972 to 1998 by a factor of about four. Water levels caused by projected ground-water withdrawals were simulated using the calibrated model. Simulations represented a period of 50 years into the future in three scenarios with either unchanged pumpage, pumpage increased by historic trends, or pumpage increased by historic trends except in two areas of the Grand Prairie. If pumping remains at 1997 rates, this produces extreme water-level declines (areas where model cells have gone dry or where the water level in the aquifer is equal to or less than the original saturated thickness, assuming confined conditions in the aquifer everywhere in the formation in predevelopment times) in the aquifer in two areas of the aquifer (one in the Grand Prairie area between the Arkansas and White Rivers and the other west of Crowleys Ridge along the Cache River) with about 400 square miles going dry. Increasing the pumping rates to that which would be projected using historic data led to increased extreme water-level declines in both areas with about 1,300 square miles going dry. Declines in both scenarios generally occurred most rapidly between 2009 and 2019. Reducing the pumping rates to 90 percent of that used for projected historic rates in areas between the Arkansas and White Rivers relating to two diversion projects of the U.S. Army Corps of Engineers and other agencies did little to decrease the extreme water-level declines. However, these pumpage reductions are small (amounting to about 16 percent of the reductions that could result from implementation of these diversion projects).

Water-Resources Investigations Report↗

Simulation of ground-water flow and movement of the freshwater-saltwater interface in the New Jersey coastal plain

The confined aquifers of the New Jersey Coastal Plain are sands that range in thickness from 50 to 600 feet and are separated by confining units. The confining units are composed of silts and clays that range in thickness from 500 to 1,000 feet. The aquifers are recharged by precipitation on their outcrop areas. This water then flows laterally downdip and vertically to the deeper confined aquifers. The confined aquifers ultimately discharge to the Raritan and Delaware Bays and to the Atlantic Ocean. In 1988, ground-water withdrawals from confined and unconfined New Jersey Coastal Plain aquifers were approximately 345 million gallons per day, more than 75 percent of which was pumped from the confined aquifers. These withdrawals have created large cones of depression in several Coastal Plain aquifers near populated areas, particularly in Camden and Monmouth Counties. The continued decline of water levels in confined aquifers can cause saltwater intrusion, reduce stream discharge near the outcrop areas, and threaten the quality of the ground-water supply. SHARP, a quasi-three-dimensional finite-difference computer model that can simulate freshwater and saltwater flow, was used to simulate the ground-water flow system in the New Jersey Coastal Plain, including the location and movement of the freshwater-saltwater interface in nine aquifers and eight intervening confining units. The freshwater-saltwater interface is defined as the hypothetical line seaward of which the chloride concentration is equal to or greater than 10,000 milligrams per liter. Model simulations were used to estimate the location and movement of the freshwater-saltwater interface resulting from (1) eustatic sea-level changes over the past 84,000 years, (2) ground-water withdrawals from 1896 through 1988, (3) and future ground-water withdrawals from 1988 to 2040 from Coastal Plain aquifers. Simultion results showed that the location and movement of the freshwater-saltwater interface are more dependent on the historical sea level than on the stresses imposed on the flow system by ground-water withdrawals from the Coastal Plain aquifers from 1896 to 1988. Results of a predictive simulation in which pumpage from existing wells was increased by 30 percent indicate that additional withdrawals from each of the eight confined aquifers in the Coastal Plain would broaden and deepen the existing cones of depression and result in significant drawdowns from the 1988 potentiometric surfaces. Drawdowns of 30 feet were simulated at the center of the cone of depression in the Upper, Middle, and Lower Potomac-Raritan-Magothy aquifers in Camden and Ocean Counties. Simulated drawdowns exceeded 80 feet at the center of the cone of depression in the Wenonah-Mount Laurel and Englishtown aquifers in Monmouth County. Drawdowns of 30 feet were simulated in the lower Kirkwood-Cohansey and confined Kirkwood aquifers in Cape May County. Simulation results showed that the increase in ground-water withdrawals would result in only minimal movement of the freshwater-saltwater interface by 2040, despite large drawdowns.

New Jersey↗

Preliminary evaluation of the importance of existing hydraulic-head observation locations to advective-transport predictions, Death Valley regional flow system, California and Nevada

When a model is calibrated by nonlinear regression, calculated diagnostic statistics and measures of uncertainty provide a wealth of information about many aspects of the system. This report presents a method of ranking the likely importance of existing observation locations using measures of prediction uncertainty. It is suggested that continued monitoring is warranted at more important locations, and unwarranted or less warranted at less important locations. The report develops the methodology and then demonstrates it using the hydraulic-head observation locations of a three-layer model of the Death Valley regional flow system. The predictions of interest are subsurface transport from beneath Yucca Mountain and 14 Underground Test Areas. The advective component of transport is considered because it is the component most affected by the system dynamics represented by the scale model being used. The problem is addressed using the capabilities of the U.S. Geological Survey computer program MODFLOW-2000, with its ADVective-Travel Observation (ADV) Package, and an additional computer program developed for this work. The methods presented in this report are used in three ways. (1) The ratings for individual observations are obtained by manipulating the measures of prediction uncertainty, and do not involve recalibrating the model. In this analysis, observation locations are each omitted individually and the resulting increase in uncertainty in the predictions is calculated. The uncertainty is quantified as standard deviations on the simulated advective transport. The increase in uncertainty is quantified as the percent increase in the standard deviations caused by omitting the one observation location from the calculation of standard deviations. In general, observation locations associated with larger increases are rated as more important. (2) Ratings for largely geographically based groups are obtained using a straightforward extension of the method used for individual observation locations. This analysis is needed where observations are clustered to determine whether the area is important to the predictions of interest. (3) Finally, the method is used to evaluate omitting a set of 100 observation locations. The locations were selected because they had low individual ratings and were not one of the few locations at which hydraulic heads from deep in the system were measured. The major results of the three analyses, when applied to the three-layer DVRFS ground-water flow system, are described in the following paragraphs. The discussion is labeled using the numbers 1 to 3 to clearly relate it to the three ways the method is used, as listed above. (1) The individual observation location analysis indicates that three observation locations are most important. They are located in Emigrant Valley, Oasis Valley, and Beatty. Of importance is that these and other observations shown to be important by this analysis are far from the travel paths considered. This displays the importance of the regional setting within which the transport occurs, the importance of including some sites throughout the area in the monitoring network, and the importance of including sites in these areas in particular. The method considered in this report indicates that the 19 observation locations that reflect hydraulic heads deeper in the system (in model layers 1, 2, and 3) are not very important. This appears to be because the locations of these observations are in the vicinity of shallow observation locations that also generally are rated as low importance, and because the model layers are hydraulically well connected vertically. The value of deep observations to testing conceptual models, however, is stressed. As a result, the deep observations are rated higher than is consistent with the results of the analysis presented, and none of these observations are omitted in the scenario discussed under (3) below. (2) The geographic grouping of th

California, Nevada↗

Projected effects of proposed chloride-control projects on shallow ground water– Preliminary results for the Wichita River basin, Texas

The U.S. Army Corps of Engineers' plan to control the natural chloride pollution in the Wichita River basin includes the construction of Truscott Brine Lake on a tributary of the North Wichita River. In connection with the proposed brine lake, the U.S. Geological Survey was requested to: (1) Define the existing ground-water conditions in the shallow fresh-water system of the project area; and (2) project the post-construction effects of the proposed lake on the fresh-water aquifer, especially in relation to hydraulic-head changes but also with respect to possible changes in the chemical quality of the ground water. The fresh-water aquifer in the project area is a shallow water-table system with relatively fresh water (generally a calcium sulfate type) that contains approximately 500 to 5,000 milligrams per liter of dissolved solids. The aquifer consists of Permian rocks with very small values of hydraulic conductivity and overlies a brine system (sodium chloride type) that is even less permeable. A thin transition zone separates the fresh-water aquifer and the brine system. Small quantities of infiltration from precipitation throughout the area's watershed constitute the recharge to the aquifer. Discharge from the aquifer consists of the small base flow along creeks; well discharge is negligible. Two-dimensional mathematical computer models were developed for aquifer simulation of: (1) Steady-state conditions in a fresh-water system and (2) transient conditions in a brine- fresh-water system where the density effects of the brine are considered. The main results 'of projecting the effects of the proposed Truscott Brine Lake on the fresh-water aquifer are: (1) Hydraulic head rises of 5 to 40 feet would be confined to areas near the proposed dam and along the lake shoreline, and (2) migration of salt water downstream from the dam generally would be limited to less than 1 mile and apparently would not reach equilibrium during the 100-year duration of the project. The modeling efforts did not include possible effects related to hydrodynamic dispersion in the brine- fresh-water system. Possible changes in the hydraulic conductivity of the aquifer, due to physical and chemical interactions in the brine and fresh-water environments, also were not considered.

Texas↗

Hydrogeologic and geochemical characteristics of the Ogallala and White River aquifers, Cheyenne, Wyoming

The Ogallala aquifer and the underlying White River aquifer are important ground-water resources of public and private drinking water in the Cheyenne, Wyoming area. In 1997, as part of a cooperative project between the Cheyenne Board of Public Utilities and the U.S. Geological Survey, a well was installed to develop information for those two aquifers. Information provided for the Ogallala aquifer included core descriptions, geophysical logs, water levels, aquifer transmissivity, water quality, isotopic analysis, and geochemical modeling. Information for the White River aquifer was limited to core descriptions and geophysical logs.Evaluation of the core obtained from the drill hole showed the sediments to be primarily sands, silts, and clays. The thickness of the Ogallala Formation at the well site was estimated to be 246 feet. Water levels and precipitation from October 1, 1998 to September 30, 1999 indicated that water levels responded to precipitation. During that time, water levels ranged from 6,002.41 feet to 6,004.27 feet above mean sea level. The transmissivity was estimated to be 1.1 feet squared/day. The temperature differences between the municipal water and the Ogallala aquifer water were examined in relation to selected hydraulic conductivities and it was found that the warmer municipal water would slightly increase the hydraulic conductivity if the water were injected into the Ogallala aquifer.The water quality of a sample from the Ogallala 1 well indicated the predominant major ions were calcium, magnesium, and bicarbonate. Isotopic analyses of hydrogen-2, tritium, chlorine-36, carbon-14, and carbon-13 indicated the water was a mixture of pre- and post-1953 recharge. A simple geochemical mixing model indicated there was the potential for dissolution of anhydrite, calcite, gypsum, and dolomite and precipitation of goethite, hematite, pyrolusite, and amorphous ferric hydroxide if municipal and Ogallala aquifer waters were mixed.

Water-Resources Investigations Report↗

Application of digital profile modeling techniques to ground-water solute transport at Barstow, California

This study investigated the use of a two-dimensional profile-oriented water-quality model for the simulation of head and water-quality changes through the saturated thickness of an aquifer. The profile model is able to simulate confined or unconfined aquifers with nonhomogeneous anisotropic hydraulic conductivity, nonhomogeneous specific storage and porosity, and nonuniform saturated thickness. An aquifer may be simulated under either steady or nonsteady flow conditions provided that the ground-water flow path along which the longitudinal axis of the model is oriented does not move in the aquifer during the simulation time period. The profile model parameters are more difficult to quantify than are the corresponding parameters for an areal-oriented water-fluality model. However, the sensitivity of the profile model to the parameters may be such that the normal error of parameter estimation will not preclude obtaining acceptable model results. Although the profile model has the advantage of being able to simulate vertical flow and water-quality changes in a single- or multiple-aquifer system, the types of problems to which it can be applied is limited by the requirements that (1) the ground-water flow path remain oriented along the longitudinal axis of the model and (2) any subsequent hydrologic factors to be evaluated using the model must be located along the land-surface trace of the model. Simulation of hypothetical ground-water management practices indicates that the profile model is applicable to problem-oriented studies and can provide quantitative results applicable to a variety of management practices. In particular, simulations of the movement and dissolved-solids concentration of a zone of degraded ground-water quality near Barstow, Calif., indicate that halting subsurface disposal of treated sewage effluent in conjunction with pumping a line of fully penetrating wells would be an effective means of controlling the movement of degraded ground water.

Water Supply Paper↗