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 703 records · Page 39Linked to original sources
USGS1995

Numerical Simulation of Regional Changes in Ground-Water Levels and in the Freshwater-Saltwater Interface Induced by Increased Pumpage at Barbers Point Shaft, Oahu, Hawaii

The effect on the regional ground-water system of southern Oahu from increased pumpage at Barbers Point shaft was estimated by a numerical ground-water model developed for the Oahu Regional Aquifer Systems Analysis (RASA) study. The RASA model was updated by revising pumping and ground-water recharge data. Pumpage data used in the new simulations were based on the allocated pumping rates for 1995 as set by the State Commission on Water Resource Management. On the basis of numerical simulation, Barbers Point shaft can sustain a withdrawal rate of 4.34 million gallons per day without adversely affecting wells in the Waianae aquifer. From results of numerical simulations, it is estimated that, as a result of increasing pumpage in Barbers Point shaft by 2 million gallons per day above the 1995-allocated rate of 2.337 million gallons per day, regional declines in ground-water levels will be about 0.4 to 0.7 feet throughout the Waianae aquifer and about 0.8 ft at the shaft. The corresponding rise of the freshwater-saltwater interface, as a result of declines in ground-water levels, is estimated to be about 20 to 30 feet. Numerical simulation also indicates that changes in ground-water levels greater than about 0.1 feet do not extend across either the Waianae-Koolau unconformity or the south Schofield barrier. The model-estimated position of the freshwater-saltwater interface, as a result of additional pumpage, ranges from 500 to 860 feet below sea level in the southern and northern parts of the aquifer, respectively, and about 540 feet below sea level at the shaft. On the basis of an estimate of the thickness of the transition-zone, the freshwater lens would remain about 240 feet thick below the shaft. In addition, the estimated declines in ground-water levels throughout the aquifer are small compared with the thickness of the freshwater lens and these declines would not be expected to affect the yields of other wells in terms of quantity. Chloride concentrations in the water pumped at Barbers Point shaft were about 240 milligrams per liter in 1992. The estimated background chloride concentration is 200 to 220 milligrams per liter because of low rainfall and the contamination of recharge water from natural salt accumulation in the soil. A reduction in irrigation through 1995 is expected to reduce recharge to the aquifer from irrigation-return water and chloride concentrations associated with the irrigation water throughout the Waianae aquifer. As a result of these combined effects, chloride concentrations of water pumped from the Barbers Point shaft will likely decrease, although the length of time required for this lowering is unknown.

Water-Resources Investigations Report↗
USGS2005

Snowmelt discharge characteristics Sierra Nevada, California

Alpine snow is an important water resource in California and the western U.S. Three major features of alpine snowmelt are the spring pulse (the first surge in snowmelt-driven river discharge in spring), maximum snowmelt discharge, and base flow (low river discharge supported by groundwater in fall). A long term data set of hydrologic measurements at 24 gage locations in 20 watersheds in the Sierra Nevada was investigated to relate patterns of snowmelt with stream discharge In wet years, the daily variations in snowmelt discharge at all the gage locations in the Sierra Nevada correlate strongly with the centrally located Merced River at Happy Isles, Yosemite National Park (i.e., in 1983, the mean of the 23 correlations was R= 0.93 + 0.09) ; in dry years, however, this correlation breaks down (i.e., in year 1977, R=0.72 + 0.24). A general trend towards earlier snowmelt was found and modeled using correlations with the timing of the spring pulse and the river discharge center of mass. For the 24 river and creek gage locations in this study, the spring pulse appeared to be a more sensitive measure of early snowmelt than the center of mass. The amplitude of maximum daily snowmelt discharge correlates strongly with initial snow water equivalent. Geologic factors, base rock permeability and soil-to-bedrock ratio, influence snowmelt flow pathways. Although both surface and ground water flows and water levels increase in wet years compared to dry years, the increase was greater for surface water in a watershed with relatively impermeable base rock than for surface water in a watershed with highly permeable base rock The relation was the opposite for base flow (ground water). The increase was greater for groundwater in a watershed with permeable rock compared to ground water in a watershed with impermeable rock. A similar, but weaker, surface/groundwater partitioning was observed in relatively impermeable granitic watersheds with differing soil-to-bedrock ratios. The increase in surface flow was greater in a watershed with a low, compared to a high, soil-to-bedrock ratio; whereas the increase in ground water flow was greater in a watershed with a high, compared to a low, soil-to-bedrock ratio. Transects that include long-term observations of shallow well-water depth and chemistry would complement traditional hydroclimate data and provide a more complete understanding of hydrologic controls of snowmelt.

California↗
USGS1984

Simulated ground-water flow in the Potomac aquifers, New Castle County, Delaware

Flow in three aquifers of the Potomac Formation in New Castle County, Delaware, was simulated. Results of the calibrated flow model show the maximum transmissivity of the lower aquifer is between 1,000 and 1,500 square feet per day; the middle aquifer, between 3,000 and 3,500 square feet per day; and the upper aquifer, between 5,000 and 6,000 square feet per day. Vertical leakance values for the three confining beds are similar and range 1 x 10 to the (1-8) power per day in downdip areas and .01 per day in subcrop areas. A storage coefficient of 0.00056 was used for each aquifer. A specific storage value of 0.000006 per foot was used for each confining bed. Estimated storage coefficients were not changed during calibration. Changes in water levels resulting from five scenarios of future pumpage were evaluated. Results of these simulations show: (1) the reduction in pumpage at Amoco causing a 120 foot local recovery of heads; (2) water levels generally will change less than 25 feet if there is no change in pumpage; (3) proposed pumpage in western Delaware may cause drawdowns below the tops of the aquifers; and (4) expected increases in pumpage would cause the greatest head decline, about 40 feet, in the lower aquifer at Getty. (USGS)

Water-Resources Investigations Report↗
USGS1997

Geochemistry and microbiology of iron-related well-screen encrustation and aquifer biofouling in Suffolk County, Long Island, New York

Iron-related well-screen encrustation and aquifer biofouling has decreased the specific capacity of several production wells in Suffolk County, N.Y., and has forced the Suffolk County Water Authority to adopt a costly well-reconditioning and replacement program. The specific-capacity declines are the result of the precipitation of iron oxyhydroxides and the growth of iron bacteria on the well screens and in the pore spaces of the surrounding formation. Mineralogic and chemical analyses indicate that the inorganic part of the encrusting material consists primarily of amorphous ferric hydroxide (Fe(OH)3 ); minor components of the material include goethite (FeOOH), hematite (Fe2 O 3 ), and quartz (SiO 2 ). The weight percent of ferric hydroxide in the material ranged from 32.3 to 98.6 percent and averaged 64.3 percent. Equilibrium modeling indicated that during pumping the well waters were supersaturated with respect to goethite, hematite, magnetite, and quartz and were under-saturated with respect to ferric hydroxide. Theoretical Eh values computed for the ferrous/ferric-iron redox couple and the oxygen/water redox couple averaged 390 millivolts and 810 millivolts, respectively, indicating that the waters were in a state of redox disequilibrium. The disequilibrium condition arises from the mixing of ground water with a low dissolved-oxygen concentration with oxygenated ground water during operation of the well. The low pH of the ground water contributes to the disequilibrium condition by slowing the rate of iron oxidation after the introduction of oxygen. Chemical and mineralogical data indicate that most of the encrusting material in the wells was deposited while the wells were shut down, probably in response to the use of treated water of higher pH to keep pump turbines wet while the wells were not in operation; the increased pH of water in the static water column increases the rate of ferrous-iron oxidation and causes the well water to become increasingly saturated with respect to ferric hydroxide. The median half-time of oxidation in samples of untreated ground water (pH 4-5) was 4.19 days, whereas the average half-time of oxidation in treated water (pH 7-8) was 11.9 minutes Equilibrium modeling indicated that treated waters generally were supersaturated with respect to ferric hydroxide, whereas untreated well waters were not. Field and laboratory data indicate that iron bacteria play an important role in the encrustation and biofouling process in Suffolk County. Filamentous iron bacteria were common in the affected wells. The most common species was Gallionella ferruginea, an effective biofouling agent that prefers water with low, but detectable, dissolved-oxygen concentrations and high dissolved-iron concentrations; this species was more common in biofilm samples from the Magothy aquifer than in those from the upper glacial aquifer. Iron bacteria also were found in sediment cores from several locations in the aquifer and in drilling water. Lignite could act as a carbon source for heterotrophic iron bacteria, which could accelerate the formation of iron-bacteria biofilms in wells screened in some parts of the Magothy aquifer. Iron-bacteria biofilms alter the chemistry of well water by removing iron, manganese, and sulfate from solution and by increasing the pH. Sulfur-reducing bacteria and iron-sulfide mineral phases were observed in some samples of encrusting material, indicating that these bacteria could contribute to well-screen encrustation in some geochemical environments.

Water-Resources Investigations Report↗
USGS1989

Simulation analysis of the ground-water system in Mesozoic rocks in the Four Corners area, Utah, Colorado, Arizona, and New Mexico

The steady-state groundwater system in Mesozoic rocks in the Four Corners area, Utah, Colorado, Arizona, and New Mexico, was simulated with a finite-difference digital-computer model to improve the understanding of the system. The simulated area is 4 ,100 sq mi, and it includes three aquifers. The Entrada-Navajo aquifer includes the Wingate, Navajo, and Entrada Sandstones. The Morrison aquifer includes the sandstone units of the Morrison Formation. The Dakota aquifer includes the Burro Canyon Formation and Dakota Sandstone. The simulation of the groundwater system had a mean error (error is absolute value of residual) of 70 ft for the Entrada-Navajo aquifer, 67 ft for the Morrison aquifer and 79 ft for the Dakota aquifer. The hydraulic conductivity used in the simulation ranged from 0.38 to 0.47 ft/day. Simulated inflow to the groundwater system was 30,000 acre-ft/yr. 48% of the inflow is from infiltration of precipitation within the simulated area, and 42% is from infiltration in 145 sq mi of mountain areas adjacent to the simulated area. Simulations indicated that some vertical inflow of water is needed between the Entrada-Navajo and Morrison aquifers to develop a reasonable representation of the system.

Arizona, Colorado, New Mexico, Utah↗
USGS1990

Ground-water resources of Honey Lake Valley, Lassen County, California, and Washoe County, Nevada

Honey Lake Valley is a 2,200 sq-mi, topographically closed basin about 35 miles northwest of Reno, Nevada. Unconsolidated basin-fill deposits on the valley floor and fractured volcanic rocks in northern and eastern uplands are the principal aquifers. In the study area, about 130,000 acre- ft of water recharges the aquifer system annually, about 40% by direct infiltration of precipitation and about 60% by infiltration of streamflow and irrigation water. Balancing this is an equal amount of groundwater discharge, of which about 65% evaporates from the water table or is transpired by phreatophytes, about 30 % is withdrawn from wells, and about 5% leaves the basin as subsurface outflow to the east. Results of a groundwater flow model of the eastern part of the basin, where withdrawals for public supply have been proposed, indicate that if 15,000 acre-ft of water were withdrawn annually, a new equilibrium would eventually be established by a reduction of about 60% in both evapotranspiration and subsurface outflow to the east. Hydrologic effects would be minimal at the western boundary of the flow-model area. Within the modeled area, the increased withdrawals cause an increase in the simulated net flow of groundwater eastward across the California-Nevada State line from about 670 acre-ft/yr to about 2,300 acre-ft/yr. (USGS)

California, Nevada↗
USGS2002

Probability distributions of hydraulic conductivity for the hydrogeologic units of the Death Valley regional ground-water flow system, Nevada and California

The use of geologic information such as lithology and rock properties is important to constrain conceptual and numerical hydrogeologic models. This geologic information is difficult to apply explicitly to numerical modeling and analyses because it tends to be qualitative rather than quantitative. This study uses a compilation of hydraulic-conductivity measurements to derive estimates of the probability distributions for several hydrogeologic units within the Death Valley regional ground-water flow system, a geologically and hydrologically complex region underlain by basin-fill sediments, volcanic, intrusive, sedimentary, and metamorphic rocks. Probability distributions of hydraulic conductivity for general rock types have been studied previously; however, this study provides more detailed definition of hydrogeologic units based on lithostratigraphy, lithology, alteration, and fracturing and compares the probability distributions to the aquifer test data. Results suggest that these probability distributions can be used for studies involving, for example, numerical flow modeling, recharge, evapotranspiration, and rainfall runoff. These probability distributions can be used for such studies involving the hydrogeologic units in the region, as well as for similar rock types elsewhere. Within the study area, fracturing appears to have the greatest influence on the hydraulic conductivity of carbonate bedrock hydrogeologic units. Similar to earlier studies, we find that alteration and welding in the Tertiary volcanic rocks greatly influence hydraulic conductivity. As alteration increases, hydraulic conductivity tends to decrease. Increasing degrees of welding appears to increase hydraulic conductivity because welding increases the brittleness of the volcanic rocks, thus increasing the amount of fracturing.

Water-Resources Investigations Report↗
USGS2008

Modeling Water-Surface Elevations and Virtual Shorelines for the Colorado River in Grand Canyon, Arizona

Using widely-available software intended for modeling rivers, a new one-dimensional hydraulic model was developed for the Colorado River through Grand Canyon from Lees Ferry to Diamond Creek. Solving one-dimensional equations of energy and continuity, the model predicts stage for a known steady-state discharge at specific locations, or cross sections, along the river corridor. This model uses 2,680 cross sections built with high-resolution digital topography of ground locations away from the river flowing at a discharge of 227 m3/s; synthetic bathymetry was created for topography submerged below the 227 m3/s water surface. The synthetic bathymetry was created by adjusting the water depth at each cross section up or down until the model?s predicted water-surface elevation closely matched a known water surface. This approach is unorthodox and offers a technique to construct one-dimensional hydraulic models of bedrock-controlled rivers where bathymetric data have not been collected. An analysis of this modeling approach shows that while effective in enabling a useful model, the synthetic bathymetry can differ from the actual bathymetry. The known water-surface profile was measured using elevation data collected in 2000 and 2002, and the model can simulate discharges up to 5,900 m3/s. In addition to the hydraulic model, GIS-based techniques were used to estimate virtual shorelines and construct inundation maps. The error of the hydraulic model in predicting stage is within 0.4 m for discharges less than 1,300 m3/s. Between 1,300-2,500 m3/s, the model accuracy is about 1.0 m, and for discharges between 2,500-5,900 m3/s, the model accuracy is on the order of 1.5 m. In the absence of large floods on the flow-regulated Colorado River in Grand Canyon, the new hydraulic model and the accompanying inundation maps are a useful resource for researchers interested in water depths, shorelines, and stage-discharge curves for flows within the river corridor with 2002 topographic conditions.

Arizona↗
USGS1996

Hydrogeology of, and ground-water flow in, a valley-fill and carbonate-rock aquifer system near Long Valley in the New Jersey Highlands

The hydrogeology of and ground-water flow in a valley-fill and carbonate-rock aquifer system were evaluated by using numerical-modeling techniques and geochemical interpretations to address concerns about the adequacy of the aquifer system to meet increasing demand for water. The study was conducted during 1987-90 by the U.S. Geological Survey, in cooperation with the New Jersey Department of Environmental Protection and Energy. The effects of recent and anticipated ground-water withdrawals on water levels, stream base flows, and water budgets were estimated. Simulation results indicate that recent withdrawals of 4.7 million gallons per day have resulted in water-level declines of up to 35 feet. Under conditions of increases in withdrawals of 121 percent, water levels would decline up to an additional 28 feet. The magnitude of predicted average base-flow depletion, when compared with historic low flows, indicates that projected increases in withdrawals may substantially deplete seasonal low flow of Drakes Brook and South Branch Raritan River. Results of a water-budget analysis indicate that the sources of water to additional supply wells would include leakage from the overlying valley-fill aquifer and induced leakage of surface water into the aquifer system. Results of water-quality analyses indicate that human activities are affecting the quality of the ground water. With the exception of an elevated iron concentration in water from one well, concentrations of inorganic constituents in water from 75 wells did not exceed New Jersey primary or secondary drinking-water regulations. Volatile organic compounds were detected in water from several wells; in two samples, concentrations of specific compounds exceeded drinking-water regulations.

Water-Resources Investigations Report↗
USGS2004

Estimation of a water budget for 1972-2000 for the Grasslands Area, central part of the Western San Joaquin Valley, California

Equitable implementation of regulations restricting discharges from agricultural drains into the San Joaquin River requires a greater understanding of the influence of extreme precipitation events on the ground-water flow system. As part of a larger investigation, this study estimated ground-water recharge and ground-water pumpage, two important components of the water budget in the Grasslands drainage area in the central part of the western San Joaquin Valley, California, for the water years 1972 through 2000. These estimates will be used as inputs to a numerical simulation model of the regional ground-water flow system in the continuing investigation. Crop-acreage and surface-water delivery data were compiled for 14 water districts and 6 other areas comprising approximately 97 percent of the 600-square-mile study area. Little ground-water pumpage data exists for the study area. A climate-based approach was employed to estimate annual water-table recharge flux and ground-water pumpage for 11 water-budget areas. Ground-water pumpage was estimated from the residual irrigation demand after crop consumption of surface water. Estimated recharge flux to the water table for the entire study area averaged 0.8 ft/yr, and estimated ground-water pumpage per unit area for the entire study area averaged 0.5 ft/yr. Increased discharges from agricultural drains in the late 1990s may have been due partly to 4 years of high recharge from precipitation over the 6-year period from 1993 to 1998. Knowledge of the ratio of annual crop water demand to annual potential evapotranspiration, expressed as an aggregate crop coefficient, Kd, will facilitate estimation of annual water-budget components in future studies. Annual aggregate crop coefficients, calculated each year for the entire study area, were nearly constant at 0.59 from 1983 to 2000, and reasonably constant at 0.53 prior to 1983. The overall trend suggests continuous reductions in recharge from irrigation over time. This reduction is most likely due to gradual improvements in irrigation management. The recharge and pumpage estimates are both sensitive to the consumption-distribution ratio, and the pumpage estimate is sensitive to the cropped acreage.

Scientific Investigations Report↗
USGS2003

Raster dataset showing the probability of elevated concentrations of nitrate in ground water in Colorado, hydrogeomorphic regions and fertilizer use estimates included

This dataset is one of eight datasets produced by this study. Four of the datasets predict the probability of detecting atrazine and(or) desethyl-atrazine (a breakdown product of atrazine) in ground water in Colorado; the other four predict the probability of detecting elevated concentrations of nitrate in ground water in Colorado. The four datasets that predict the probability of atrazine and(or) desethyl-atrazine (atrazine/DEA) are differentiated by whether or not they incorporated atrazine use and whether or not they incorporated hydrogeomorphic regions. The four datasets that predict the probability of elevated concentrations of nitrate are differentiated by whether or not they incorporated fertilizer use and whether or not they incorporated hydrogeomorphic regions. Each of the eight datasets has its own unique strengths and weaknesses. The user is cautioned to read Rupert (2003, Probability of detecting atrazine/desethyl-atrazine and elevated concentrations of nitrate in ground water in Colorado: U.S. Geological Survey Water-Resources Investigations Report 02-4269, 35 p., http://water.usgs.gov/pubs/wri/wri02-4269/) to determine if he(she) is using the most appropriate dataset for his(her) particular needs. This dataset specifically predicts the probability of detecting elevated concentrations of nitrate in ground water in Colorado with hydrogeomorphic regions and fertilizer use included. The following text was extracted from Rupert (2003). Draft Federal regulations may require that each State develop a State Pesticide Management Plan for the herbicides atrazine, alachlor, metolachlor, and simazine. Maps were developed that the State of Colorado could use to predict the probability of detecting atrazine/DEA in ground water in Colorado. These maps can be incorporated into the State Pesticide Management Plan and can help provide a sound hydrogeologic basis for atrazine management in Colorado. Maps showing the probability of detecting elevated nitrite plus nitrate as nitrogen (nitrate) concentrations in ground water in Colorado also were developed because nitrate is a contaminant of concern in many areas of Colorado. Maps showing the probability of detecting atrazine/DEA at or greater than concentrations of 0.1 microgram per liter and nitrate concentrations in ground water greater than 5 milligrams per liter were developed as follows: (1) Ground-water quality data were overlaid with anthropogenic and hydrogeologic data by using a geographic information system (GIS) to produce a dataset in which each well had corresponding data on atrazine use, fertilizer use, geology, hydrogeomorphic regions, land cover, precipitation, soils, and well construction. These data then were downloaded to a statistical software package for analysis by logistic regression. (2) Relations were observed between ground-water quality and the percentage of land-cover categories within circular regions (buffers) around wells. Several buffer sizes were evaluated; the buffer size that provided the strongest relation was selected for use in the logistic regression models. (3) Relations between concentrations of atrazine/DEA and nitrate in ground water and atrazine use, fertilizer use, geology, hydrogeomorphic regions, land cover, precipitation, soils, and well-construction data were evaluated, and several preliminary multivariate models with various combinations of independent variables were constructed. (4) The multivariate models that best predicted the presence of atrazine/DEA and elevated concentrations of nitrate in ground water were selected. (5) The accuracy of the multivariate models was confirmed by validating the models with an independent set of ground-water quality data. (6) The multivariate models were entered into a geographic information system and the probability GRIDS were constructed.

Open-File Report↗
USGS2010

Marine benthic habitat mapping of Muir Inlet, Glacier Bay National Park and Preserve, Alaska, with an evaluation of the Coastal and Marine Ecological Classification Standard III

Seafloor geology and potential benthic habitats were mapped in Muir Inlet, Glacier Bay National Park and Preserve, Alaska, using multibeam sonar, ground-truth information, and geological interpretations. Muir Inlet is a recently deglaciated fjord that is under the influence of glacial and paraglacial marine processes. High glacially derived sediment and meltwater fluxes, slope instabilities, and variable bathymetry result in a highly dynamic estuarine environment and benthic ecosystem. We characterize the fjord seafloor and potential benthic habitats using the Coastal and Marine Ecological Classification Standard (CMECS) recently developed by the National Oceanic and Atmospheric Administration (NOAA) and NatureServe. Substrates within Muir Inlet are dominated by mud, derived from the high glacial debris flux. Water-column characteristics are derived from a combination of conductivity temperature depth (CTD) measurements and circulation-model results. We also present modern glaciomarine sediment accumulation data from quantitative differential bathymetry. These data show Muir Inlet is divided into two contrasting environments: a dynamic upper fjord and a relatively static lower fjord. The accompanying maps represent the first publicly available high-resolution bathymetric surveys of Muir Inlet. The results of these analyses serve as a test of the CMECS and as a baseline for continued mapping and correlations among seafloor substrate, benthic habitats, and glaciomarine processes.

Scientific Investigations Map↗
USGS1991

Documentation of a computer program to simulate aquifer-system compaction using the modular finite-difference ground-water flow model

Removal of ground water by pumping from aquifers may result in compaction of compressible fine-grained beds that are within or adjacent to the aquifers. Compaction of the sediments and resulting land subsidence may be permanent if the head declines result in vertical stresses beyond the previous maximum stress. The process of permanent compaction is not routinely included in simulations of ground-water flow. To simulate storage changes from both elastic and inelastic compaction, a computer program was written for use with the U.S. Geological Survey modular finite-difference ground- water flow model. The new program, the Interbed-Storage Package, is designed to be incorporated into this model. In the Interbed-Storage Package, elastic compaction or expansion is assumed to be proportional to change in head. The constant of proportionality is the product of the skeletal component of elastic specific storage and the thickness of the sediments. Similarly, inelastic compaction is assumed to be proportional to decline in head. The constant of proportionality is the product of the skeletal component of inelastic specific storage and the thickness of the sediments. Storage changes are incorporated into the ground-water flow model by adding an additional term to the right-hand side of the flow equation. Within a model time step, the package appropriately apportions storage changes between elastic and inelastic components on the basis of the relation of simulated head to the previous minimum (preconsolidation) head. Two tests were performed to verify that the package works correctly. The first test compared model-calculated storage and compaction changes to hand-calculated values for a three-dimensional simulation. Model and hand-calculated values were essentially equal. The second test was performed to compare the results of the Interbed-Storage Package with results of the one-dimensional Helm compaction model. This test problem simulated compaction in doubly draining confining beds stressed by head changes in adjacent aquifers. The Interbed-Storage Package and the Helm model computed essentially equal values of compaction. Documentation of the Interbed-Storage Package includes data input instructions, flow charts, narratives, and listings for each of the five modules included in the package. The documentation also includes an appendix describing input instructions and a listing of a computer program for time-variant specified-head boundaries. That package was developed to reduce the amount of data input and output associated with one of the Interbed-Storage Package test problems.

Techniques of Water-Resources Investigations↗
USGS1985

Simulated effects of increased recharge on the ground-water flow system of Yucca Mountain and vicinity, Nevada-California

A study was performed to assess the potential effects of changes in future climatic conditions on the groundwater system in the vicinity of Yucca Mountain, the site of a potential mined geologic repository for high-level nuclear wastes. These changes probably would result in greater rates of precipitation and, consequently, greater rates of recharge. The study was performed by simulating the groundwater system, using a two-dimensional, finite-element, groundwater flow model. The simulated position of the water table rose as much as 130 meters near the U.S. Department of Energy 's preferred repository area at Yucca Mountain for a simulation involving a 100-percent increase in precipitation compared to modern-day conditions. Despite the water table rise, no flooding of the potential repository would occur at its current proposed location. According to the simulation, springs would discharge south and west of Timber Mountain, along Fortymile Canyon, in the Amargosa Desert near Lathrop Wells and Franklin Lake playa, and near Furnace Creek Ranch in Death Valley, where they presently discharge. Simulated directions of groundwater flow paths near the potential repository area generally would be the same for the baseline (modern-day climate) and the increased-recharge simulations, but the magnitude of flow would increase by 2 to 4 times that of the baseline-simulation flow. (USGS)

California, Nevada↗
USGS2001

Simulation of ground-water flow in the Mojave River basin, California

The proximity of the Mojave River ground-water basin to the highly urbanized Los Angeles region has led to rapid growth in population and, consequently, to an increase in the demand for water. The Mojave River, the primary source of surface water for the region, normally is dry-except for a small stretch of perennial flow and periods of flow after intense storms. Thus, the region relies almost entirely on ground water to meet its agricultural and municipal needs. Ground-water withdrawal since the late 1800's has resulted in discharge, primarily from pumping wells, that exceeds natural recharge. To better understand the relation between the regional and the floodplain aquifer systems and to develop a management tool that could be used to estimate the effects that future stresses may have on the ground-water system, a numerical ground-water flow model of the Mojave River ground-water basin was developed, in part, on the basis of a previously developed analog model. The ground-water flow model has two horizontal layers; the top layer (layer 1) corresponds to the floodplain aquifer and the bottom layer (layer 2) corresponds to the regional aquifer. There are 161 rows and 200 columns with a horizontal grid spacing of 2,000 by 2,000 feet. Two stress periods (wet and dry) per year are used where the duration of each stress period is a function of the occurrence, quantity of discharge, and length of stormflow from the headwaters each year. A steady-state model provided initial conditions for the transient-state simulation. The model was calibrated to transient-state conditions (1931-94) using a trial-and-error approach. The transient-state simulation results are in good agreement with measured data. Under transient-state conditions, the simulated floodplain aquifer and regional aquifer hydrographs matched the general trends observed for the measured water levels. The simulated streamflow hydrographs matched wet stress period average flow rates and times of no flow at the Barstow and Afton Canyon gages. Steady-state particle-tracking was used to estimate travel times for mountain-front and streamflow recharge. The simulated travel times for mountain-front recharge to reach the area west of Victorville were about 5,000 to 6,000 years; this result is in reasonable agreement with published results. Steady-state particle-tracking results for streamflow recharge indicate that in most subareas along the river, the particles quickly leave and reenter the river. The complaint that resulted in the adjudication of the Mojave River ground-water basin alleged that the cumulative water production upstream of the city of Barstow had overdrafted the ground-water basin. In order to ascertain the effect of pumping on ground-water and surface-water relations along the Mojave River, two pumping simulations were compared with the 1931-90 transient-state simulation (base case). The first simulation assumed 1931-90 pumping in the upper region (Este, Oeste, Alto, and Transition zone model subareas) but with no pumping in the remainder of the basin, and the second assumed 1931-90 pumping in the lower region (Centro, Harper Lake, Baja, Coyote Lake, and Afton Canyon model subareas) but with no pumping in remainder of the basin. In the upper region, assuming pumping only in the upper region, there was no change in storage, recharge from the Mojave River, ground-water discharge to the Mojave River, or evapotranspiration when compared with the base case. In the lower region, assuming pumping only in the upper region, there was storage accretion, decreased recharge from the Mojave River, increased ground-water discharge to the Mojave River, and increased evapotranspiration when compared with the base case. In the upper region, assuming pumping only in the lower region, there was storage accretion, decreased recharge from the Mojave River, increased ground-water discharge to the Mojave River, and increased evapotranspiration when compared with the base case. In the lower region, assuming pumping only in the lower region, there was less storage depletion, increased recharge from the Mojave River, increased ground-water discharge to the Mojave River, and increased evapotranspiration when compared with the base case. Overall, pumping in the lower region does not negatively affect the upper region; however, pumping in the upper region negatively affects the lower region by decreasing recharge from the Mojave River. Streamflow, pumpage, and water-level data from calendar years 1995-99 were used to validate the calibrated ground-water flow model, that is, to test that the ground-water flow model will duplicate measured data for a noncalibration period without modification of the model parameters. In general, the simulated results are in good agreement with the measured data, and the simulated hydrographs for wells in the floodplain and regional aquifers follow the measured water-level trends. Simulated streamflow data for the 1995-99 wet and dry stress periods at the Lower Narrows, Barstow, and Afton Canyon were compared with the measured data for average streamflow for the same periods; in general, the model reflects 1995-99 streamflow conditions. The simulation results also indicate that the streambed conductance values calibrated to the 1931-94 conditions reasonably simulate the 1995-99 conditions and therefore can be used for predictive purposes. To visualize the magnitude, spatial distribution, and timing of water-level changes in the basin through time, simulated hydraulic heads for 1932-99 were compared with simulated hydraulic heads for 1931. Greater than average annual inflows to the Mojave River from the headwaters during the late 1930's and throughout much of the 1940's resulted in simulated hydraulic heads that were higher than the 1931 hydraulic heads along the Mojave River in most model subareas. Parts of the Baja and Harper Lake model subareas had declines in the simulated hydraulic head because of the increase in agricultural pumpage. By 1960, the simulated hydraulic heads were lower than the simulated hydraulic heads for 1931 in all model subareas of the floodplain and the regional aquifers because of pumpage. After 1960, the size and the magnitude of the areas of the regional aquifer for which simulated hydraulic heads were lower than those for 1931 continued to increase until the end of the simulation (1999). Along the Mojave River, hydraulic heads fluctuated in the floodplain aquifer in response to recharge during years with large inflows with little apparent effect on the simulated hydraulic heads in the regional aquifer. Three water-management alternatives were evaluated to determine their effect on ground-water resources using the calibrated ground-water flow model. The water-management alternatives consider the artificial recharge of imported water allocated to the Mojave Water Agency (MWA): the first assumes that zero percent of the MWA allocation is available (alternative 1), the second assumes that 50 percent of the MWA allocation is available (alternative 2), and the third assumes that 100 percent of the MWA allocation is available (alternative 3). Each of the three water-management alternatives were evaluated for a 20-year drought. Streamflow conditions were simulated using the 20-year drought of 1945-64 with associated calibrated stream parameters. Management alternative 1 results in a reduction in ground-water recharge from the Mojave River compared with average recharge for 1995-99; this reduction is reflected in simulated hydraulic-head declines between 1999 and 2019 of as much as 45 feet. Management alternatives 2 and 3 result in no change in recharge from the Mojave River for management alternative 2 and a small increase for management alternative 3 when compared with recharge for management alternative 1. The artificial recharge of imported water causes increases in simulated hydraulic head for both management alternatives at each of the artificial-recharge sites. Some of the increases are related to water that recharges into areas of low transmissivity which implies that the recharge operations may benefit from being distributed over a larger area.

California↗
USGS1982

A digital-computer model of the Big Sioux aquifer in Minnehaha County, South Dakota

A finite-difference digital model was used to simulate steady-state conditions of the Big Sioux aquifer in Minnehaha County. Average water levels and average base flow discharge (4.9 cu ft/s) of the Big Sioux River were based on data from 1970 through 1979. The computer model was calibrated for transient conditions by simulating monthly historic conditions for 1976. During 1976, pumpage was offset mostly by surface-water recharge to the aquifer from January through June and ground-water discharge from storage from July through December. Measured drawdowns during 1976 generally were less than 2 feet except in the Sioux Falls city well field where drawdowns were as much as 15 feet. The model was used to study the effects of increased withdrawals under three hypothetical hydrologic situations. One hypothetical situation consisted of using 1976 pumping rates, recharge, and evapotranspiration but the Big Sioux River dry. The pumping rate after 16 months was decreased by 40 percent from the actual pumping rate for that month in order to complete the monthly simulation without the storage being depleted at a nodal area. The second hypothetical situation consisted of a pumpage rate of 44.4 cubic feet per second from 60 wells spaced throughout the aquifer under historic 1976 hydrologic conditions. The results were that the aquifer could supply the additional withdrawal. The third hypothetical situation used the same hydrologic conditions as the second except that recharge was zero and the Big Sioux River was dry downstream from row 54. After 18 monthly simulations, the pumping rate was decreased by 44 percent to prevent pumping wells from depleting the aquifer, and, at that rate, 63 percent of the water being pumped was being replaced by water from the river. (USGS)

Water-Resources Investigations Report↗
USGS1987

Movement and fate of crude-oil in contaminants in the subsurface environment at Bemidji, Minnesota: Chapter C in U.S. Geological Survey program on toxic waste--ground-water contamination: Proceedings of the Third technical meeting, Pensacola, Florida, March 23-27, 1987

On August 20, 1979, a pipeline break in a remote area near Bemidji, Minn. (fig. C-l), resulted in the release of 1.5x10 5 L (liters) of crude oil. Although about 1.1x10 5 L were removed from the site as part of the cleanup, some crude oil infiltrated the ground and percolated to the water table. The spill occurred in the recharge area of a local flow system that discharges to a small lake 300 m (meters) downgradient (Hult, 1984). The aquifer is a pitted and dissected outwash plain underlain at a depth of about 20 m by low-permeability till. Crude oil is floating on the water table about 8 m below land surface and has migrated about 20 m as a separate fluid phase. Soluble petroleum derivatives have dissolved in and are moving with ground water. Volatile constituents are migrating through the unsaturated zone by diffusion. The abstracts presented in this chapter provide an overview of ongoing efforts to combine the results of interdependent, interdisciplinary research into the comprehensive understanding of the physical, chemical, and biological processes that will be needed to develop predictive models of contaminant mobilization, transport, and fate. Franzi investigates the relationship between depositional and post-depositional processes and the heterogeneity and anisotropy of the aquifer. Post-depositional subsidence and collapse caused by melting of stagnant ice has created significant geomorphic and hydrogeologic structures. Miller describes a preliminary ground-water flow and chemical-transport model used to assess the effect of these hydrologic discontinuities, estimate aquifer properties, and to guide continuing field work. Abstracts by Siegel, Bennett, and Berndt describe studies of dissolved inorganic constituents and parameters done to characterize the geochemical environments in and around the site. An anoxic zone, probably plume shaped, extends downgradient of the oil pool and into an oxic zone that surrounds the contamination. Carbonate minerals are actively dissolving in the oxic spray zone upgradient of the oil pool where oil is being mineralized, whereas quartz is dissolving underneath the oil by organic-acid complexation. Selective leaching experiments show that fractionation of metals is also occurring where the aquifer matrix is contaminated. Morphological analysis of quartz sand grains clearly show unusual rock-water interactions are occurring in the anaerobic contaminant zone. Preliminary analysis indicates the possibility of an organic/quartz interaction that is significantly increasing the mobility of silica through the system. Predictions of the evolution and ultimate geometry of contaminant plumes resulting from spills require quantitative descriptions of the rate of mass transfer from the organic fluid to ground water. Pfannkuch presents laboratory and field work that describe how the the rate of oil dissolution, and therefore the strength of the contaminant source, is controlled by fluctuations in ground-water velocity and water-table fluctuations.

Minnesota↗
USGS1996

A three-dimensional method-of-characteristics solute-transport model (MOC3D)

This report presents a model, MOC3D, that simulates three-dimensional solute transport in flowing ground water. The model computes changes in concentration of a single dissolved chemical constituent over time that are caused by advective transport, hydrodynamic dispersion (including both mechanical dispersion and diffusion), mixing (or dilution) from fluid sources, and mathematically simple chemical reactions (including linear sorption, which is represented by a retardation factor, and decay). The transport model is integrated with MODFLOW, a three-dimensional ground-water flow model that uses implicit finite-difference methods to solve the transient flow equation. MOC3D uses the method of characteristics to solve the transport equation on the basis of the hydraulic gradients computed with MODFLOW for a given time step. This implementation of the method of characteristics uses particle tracking to represent advective transport and explicit finite-difference methods to calculate the effects of other processes. However, the explicit procedure has several stability criteria that may limit the size of time increments for solving the transport equation; these are automatically determined by the program. For improved efficiency, the user can apply MOC3D to a subgrid of the primary MODFLOW grid that is used to solve the flow equation. However, the transport subgrid must have uniform grid spacing along rows and columns. The report includes a description of the theoretical basis of the model, a detailed description of input requirements and output options, and the results of model testing and evaluation. The model was evaluated for several problems for which exact analytical solutions are available and by benchmarking against other numerical codes for selected complex problems for which no exact solutions are available. These test results indicate that the model is very accurate for a wide range of conditions and yields minimal numerical dispersion for advection-dominated problems. Mass-balance errors are generally less than 10 percent, and tend to decrease and stabilize with time.

Water-Resources Investigations Report↗