SEARCH · Search USGS

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.

SourceResearch connectionFrom yearThrough year

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

At least 829 records · Page 46Linked to original sources

USGS2002

Simulated transport and biodegradation of chlorinated ethenes in a fractured dolomite aquifer near Niagara Falls, New York

Leakage of trichloroethene (TCE) from a neutralization pond at a former manufacturing facility near Niagara Falls, N.Y. during 1950-87 into the Guelph Formation of the Lockport Group, a fractured dolomite aquifer, created a plume of TCE and its metabolites that, by 1990, extended about 4,300 feet south of the facility. A smaller plume of dense, nonaqueous-phase liquids (DNAPL) probably serves as a continuing source of TCE. The presence of the TCE metabolites cis-1,2-dichloroethene (DCE), vinyl chloride (VC), and ethene in the plume, and the results of previous laboratory microcosm studies, indicate that the TCE is being degraded by naturally occurring microorganisms. Biodegradation rates of TCE and its metabolites were estimated through simulation with BIOMOC, a solute-transport model that represents multispecies reactions through Monod kinetics. A fracture zone in the Guelph Formation was represented as a porous medium containing an extensive, 3-foot thick layer with several interconnected fractures; this layer is bounded above and below by subhorizontal stratigraphic contacts. The Monod reaction constants were estimated through nonlinear regression to minimize the difference between computed concentrations of TCE and its metabolites, and the concentrations measured before and during 5 years of pump-and-treat remediation. Transport simulations indicated that, by April 1998, the chlorinated ethene plume had reached a dynamic equilibrium between the rate of TCE dissolution and the rate of removal through pumping and biodegradation. Biodegradation of chlorinated ethenes at the site can be simulated as first-order reactions because the concentrations are generally less than the half-saturation constants estimated for Monod kinetics (320 mg/L for TCE, 10 mg/L for DCE, and 1 mg/L for VC). Computed degradation rates are proportional to the estimated ground-water velocity, which could vary by more than an order magnitude at the site, as indicated by the estimated range of fracture porosity--3 to 0.3 percent. Half-lives corresponding to first-order rate constants estimated for the lower velocity (5 to 15 feet per day) ranged from 21 to 25 days for TCE, 170 to 230 days for DCE, and 18 to 23 days for VC. Chlorinated ethene concentrations of April 1998 were better reproduced when the TCE source was represented as a constant concentration than as a constant flux, because the latter predicted that the plume would dissipate after 5 years of pump-and-treat remediation. This result indicates that the rate of TCE dissolution is not limited by the mass of TCE in the DNAPL plume. Simulation of diffusion by the transport model MOC3D indicated that concentrations of these contaminants within the rock matrix surrounding the fracture zone were relatively unchanged after 5 years of pump-and-treat remediation. The principal sources of uncertainty in the prediction of biodegradation rates and of the fate of chlorinated ethenes at the site are the fracture porosity and DNAPL mass in the aquifer.

New York↗

USGS1998

Ground Water in the Southern Lihue Basin, Kauai, Hawaii

A multi-phased study of ground-water resources, including well drilling, aquifer tests, analysis of ground-water discharge, and numerical ground-water modeling, indicates that the rocks of the southern Lihue Basin, Kauai, have permeabilities that are much lower than in most other areas of ground-water development in the Hawaiian islands. The regional hydraulic conductivity of the Koloa Volcanics, which dominates fresh ground-water flow in the basin, is about 0.275 foot per day. The Waimea Canyon Basalt which surrounds the basin and underlies the Koloa Volcanics within the basin is intruded by dikes that reduce the bulk hydraulic conductivity of the rocks to about 1.11 feet per day. The low permeabilities result in steeper head gradients compared with other areas in the Hawaiian islands, and a higher proportion of ground-water discharging to streams than to the ocean. Water levels rise from near sea level at the coast to several hundreds of feet above sea level at the center of the basin a few miles inland. The high inland water levels are part of a completely saturated ground-water system. Because of the low regional hydraulic conductivity and high influx of water from recharge in the southern Lihue Basin, the rocks become saturated nearly to the surface and a variably saturated/unsaturated (perched) condition is not likely to exist. Streams incising the upper part of the aquifer drain ground water and keep the water levels just below the surface in most places. Streams thus play an important role in shaping the water table in the southern Lihue Basin. At least 62 percent of the ground water discharging from the aquifer in the southern Lihue Basin seeps to streams; the remainder seeps directly to the ocean or is withdrawn by wells.

Water-Resources Investigations Report↗

USGS2005

Reconnaissance of the hydrogeology of Ta'u, American Samoa

Analysis of existing data and information collected on a reconnaissance field visit supports a conceptual model of ground-water occurrence in Ta'u, American Samoa, in which a thin freshwater lens exists in a predominantly high-permeability aquifer that receives high rates of recharge. Because the freshwater lens is thin throughout most of the island, the productivity of wells, especially those near the coast where the lens is the thinnest, is likely to be limited by saltwater intrusion. The landfill in northwestern Ta'u is closer to the north coast of the island than to any of the existing or proposed well sites. Although this may indicate that ground water beneath the landfill would flow away from the existing and proposed well sites, this interpretation may change depending on the hydraulic properties of a fault and rift zone in the area. Of four plausible scenarios tested with a numerical ground-water flow model, only one scenario indicated that ground water from beneath the landfill would flow toward the existing and proposed well sites; the analysis does not, however, assess which of the four scenarios is most plausible. The analysis also does not consider the change in flow paths that will result from ground-water withdrawals, dispersion of contaminants during transport by ground water, other plausible hydrogeologic scenarios, transport of contaminants by surface-water flow, or that sources of contamination other than the landfill may exist. Accuracy of the hydrologic interpretations in this study is limited by the relatively sparse data available for Ta'u. Understanding water resources on Ta'u can be advanced by monitoring rainfall, stream-flow, evaporation, ground-water withdrawals, and water quality, and with accurate surveys of measuring point elevations for all wells and careful testing of well-performance. Assessing the potential for contaminants in the landfill to reach existing and proposed well sites can be improved with additional information on the landfill itself (history, construction, contents, water chemistry), surface-water flow directions, spatial distribution of ground-water levels, and the quality of water in nearby wells. Monitoring water levels and chemistry in one or more monitoring wells between the landfill and existing or proposed wells can provide a means to detect movement of contaminants before they reach production wells. Steps that can be implemented in the short term include analyzing water in the landfill and monitoring of water chemistry and water levels in all existing and new production wells. Placing future wells farther inland may mitigate saltwater intrusion problems, but the steep topography of Ta'u limits the feasibility of this approach. Alternative solutions include distributing ground-water withdrawal among several shallow-penetrating, low-yield wells.

Scientific Investigations Report↗

USGS1976

Finite difference model for aquifer simulation in two dimensions with results of numerical experiments

The model will simulate ground-water flow in an artesian aquifer, a water-table aquifer, or a combined artesian and water-table aquifer. The aquifer may be heterogeneous and anisotropic and have irregular boundaries. The source term in the flow equation may include well discharge, constant recharge, leakage from confining beds in which the effects of storage are considered, and evapotranspiration as a linear function of depth to water. The theoretical development includes presentation of the appropriate flow equations and derivation of the finite-difference approximations (written for a variable grid). The documentation emphasizes the numerical techniques that can be used for solving the simultaneous equations and describes the results of numerical experiments using these techniques. Of the three numerical techniques available in the model, the strongly implicit procedure, in general, requires less computer time and has fewer numerical difficulties than do the iterative alternating direction implicit procedure and line successive overrelaxation (which includes a two-dimensional correction procedure to accelerate convergence). The documentation includes a flow chart, program listing, an example simulation, and sections on designing an aquifer model and requirements for data input. It illustrates how model results can be presented on the line printer and pen plotters with a program that utilizes the graphical display software available from the Geological Survey Computer Center Division. In addition the model includes options for reading input data from a disk and writing intermediate results on a disk.

Techniques of Water-Resources Investigations↗

USGS1988

Water resources and effects of ground-water development in Pasco County, Florida

Ninety-nine percent of the 79.72 million gal/d of water used in Pasco County, Florida is groundwater from the Upper Floridan aquifer. In addition, 53.5 million gal/d is exported for use in Pinellas County. Chemical quality of the water generally falls within recommended limits for drinking water except near the coast. A groundwater flow model was used to predict drawdowns in water levels under five different development plans for west Pasco County. One of these plans was incorporated in an estimated groundwater development plan for 2035 for Pasco, Pinellas, and part of Hillsborough County. Reduction in evapotranspiration accounts for nearly all the water required for each development plan. Two plans showed less effect in drawdown and potential of saltwater intrusion than the other three. For the overall groundwater development plan for Pasco, Pinellas, and parts of Hillsborough Counties, the average potentiometric surface in 2035 would be from 8 ft higher to 20 ft lower than the average 1976-77 potentiometric surface. Reductions in head would increase the potential for infiltration of contaminants in areas where surficial materials are thin, sinkhole development in sinkhole prone areas, and upconing and lateral intrusion of saltwater. The potential for dewatering the surficial aquifer is great in the Cross Bar Ranch and Cypress Creek well-field areas. Evapotranspiration of groundwater would be reduced by 19% between 1976-77 and 2035. Exchange of water between the rivers and the aquifers will be reduced by 13%. Springflow will be reduced by 6%, and three springs will cease to flow. Model boundary inflow and outflow also will be reduced.

Florida↗

USGS1986

Ground-water flow in Melton Valley, Oak Ridge reservation, Roane County, Tennessee; preliminary model analysis

Shallow land burial of low-level radioactive waste has been practiced since 1951 in Melton Valley. Groundwater flow modeling was used to better understand the geohydrology of the valley, and to provide a foundation for future contaminant transport modeling. The three-dimensional, finite difference model simulates the aquifer as a two layer system that represents the regolith and bedrock. Transmissivities, which were adjusted during model calibration, range from 8 to 16 sq ft/day for the regolith, and from 0.2 to 1.5 sq ft/day for bedrock. An anisotropy ratio of 1:3 for strike-normal to strike-parallel transmissivity values, in conjunction with recharge rate = 6% of precipitation that is uniformly distributed over the model area, produces the best match between simulated and observed water levels. Simulated water levels generally compare well to observed or estimated 1978 groundwater conditions. Simulated water levels for the regolith for 39 of 69 comparison points are within +/- 10 ft of average 1978 levels. Simulated vertical flow components are in the observed direction for 9 of 11 comparison points. Preliminary simulations indicate that nearly all groundwater flow is within the regolith and discharges to either the Clinch River or the White Oak Creek-Melton Branch drainage systems. Less than 3% of the flow is between the regolith and bedrock, and < 1% of total groundwater flow discharges to the Clinch River through bedrock. Additional data needed to refine and further calibrate the model, include: (1) quantity and areal distribution of recharge; (2) water levels in the regolith near the model boundaries and beyond the Clinch River; (3) water levels and aquifer characteristics for bedrock; and (4) additional surface water data. (Author 's abstract)

Tennessee↗

USGS2004

Conceptualization and simulation of the Edwards aquifer, San Antonio region, Texas

A new numerical ground-water-flow model (Edwards aquifer model) that incorporates important components of the latest information and plausible conceptualization of the Edwards aquifer was developed. The model includes both the San Antonio and Barton Springs segments of the Edwards aquifer in the San Antonio region, Texas, and was calibrated for steady-state (1939–46) and transient (1947–2000) conditions, excluding Travis County. Transient simulations were conducted using monthly recharge and pumpage (withdrawal) data. The model incorporates conduits simulated as continuously connected (other than being separated in eastern Uvalde and southwestern Medina Counties), one-cell-wide (1,320 feet) zones with very large hydraulic-conductivity values (as much as 300,000 feet per day). The locations of the conduits were based on a number of factors, including major potentiometric-surface troughs in the aquifer, the presence of sinking streams, geochemical information, and geologic structures (for example, faults and grabens). The simulated directions of flow in the Edwards aquifer model are most strongly influenced by the presence of simulated conduits and barrier faults. The simulated flow in the Edwards aquifer is influenced by the locations of the simulated conduits, which tend to facilitate flow. The simulated subregional flow directions generally are toward the nearest conduit and subsequently along the conduits from the recharge zone into the confined zone and toward the major springs. Structures simulated in the Edwards aquifer model influencing groundwater flow that tend to restrict flow are barrier faults. The influence of simulated barrier faults on flow directions is most evident in northern Medina County. A water budget is an accounting of inflow to, outflow from, and storage change in the aquifer. For the Edwards aquifer model steady-state simulation, recharge (from seepage losses from streams and infiltration of rainfall) accounts for 93.5 percent of the sources of water to the Edwards aquifer, and inflow through the northern and northwestern model boundaries contributes 6.5 percent. The largest discharges are springflow (73.7 percent) and ground-water withdrawals by wells (25.7 percent). The principal source of water to the Edwards aquifer for the Edwards aquifer model transient simulation was recharge, constituting about 60 percent of the sources of water (excluding change in storage) to the Edwards aquifer during 1956, a drought period, and about 97 percent of the sources (excluding change in storage) during 1975, a period of above-normal rainfall and recharge. The principal discharges from the Edwards aquifer for the transient simulation were springflow and withdrawals by wells. During 1956, representing drought conditions, the change in storage (net water released from storage) was much greater than recharge, accounting for 75.9 percent of the total flow compared to 14.5 percent for recharge. Conversely, during 1975, representing above-normal rainfall and recharge conditions, recharge constituted 79.9 percent of the total flow, compared to 7.1 percent for the change in storage (net water added to storage). A series of sensitivity tests was made to ascertain how the model results were affected by variations greater than and less than the calibrated values of input data. Simulated hydraulic heads in the Edwards aquifer model were most sensitive to recharge, withdrawals, hydraulic conductivity of the conduit segments, and specific yield and were comparatively insensitive to spring-orifice conductance, northern boundary inflow, and specific storage. Simulated springflow in the Edwards aquifer model was most sensitive to recharge, withdrawals, hydraulic conductivity of the conduit segments, specific yield, and increases in northern boundary inflow and was comparatively insensitive to spring-orifice conductance and specific storage.

Texas↗

USGS1988

Simulation of flow in the Edwards Aquifer, San Antonio region, Texas, and refinement of storage and flow concepts

The Edwards aquifer is a complexly faulted, carbonate aquifer lying within the Balcones fault zone of south-central Texas. The aquifer consists of thin- to massive-bedded limestone and dolomite, most of which is in the form of mudstones and wackestones. Well-developed secondary porosity has formed in association with former erosional surfaces within the carbonate rocks, within dolomitized-burrowed tidal and evaporitic deposits, and along inclined fractures to produce an aquifer with transmissivities greater than 100 ft2/s. The aquifer is recharged mainly by streamflow losses in the outcrop area of the Edwards aquifer and is discharged by major springs located at considerable distances, as much as 150 mi, from the areas of recharge and by wells. Ground-water flow within the Edwards aquifer of the San Antonio region was simulated to investigate concepts relating to the storage and flow characteristics. The concepts of major interest were the effects of barrier faults on flow direction, water levels, springflow, and storage within the aquifer. A general-purpose, finite-difference model, modified to provide the capability of representing barrier faults, was used to simulate ground-water flow and storage in the aquifer. The approach in model development was to conduct a series of simulations beginning with a simple representation of the aquifer framework and then proceeding to subsequent representations of increasing complexity. The simulations investigated the effects of complex geologic structures and of significant changes in transmissivity, anisotropy, and storage coefficient. Initial values of transmissivity, anisotropy, and storage coefficient were estimated based on concepts developed in previous studies. Results of the simulations confirmed the original estimates of transmissivity values (greater than 100 square feet/s) in the confined zone of the aquifer between San Antonio and Comal Springs. A storage coefficient of 0.05 in the unconfined zone of the aquifer produced the best simulation of water levels and springflow. A major interpretation resulting from the simulations is that two essentially independent areas of regional flow were identified in the west and central part of the study area. Flows from the two areas converge at Comal Springs. The directions of computed flux vectors reflected the presence of major barrier faults, which locally deflect patterns of ground-water movement. The most noticeable deflection is the convergence of flow through a geologic structural opening, the Knippa gap, in eastern Uvalde County. A second significant interpretation is that ground-water flow in northeastern Bexar, Comal, and Hays Counties is diverted by barrier faults toward San Marcos Springs, a regional discharge point. Simulations showed that several barrier faults in the northwestern part of the San Antonio area had a significant effect on storage, water levels, and springflow within the Edwards aquifer.

Texas↗

USGS1986

Technique for predicting ground-water discharge to surface coal mines and resulting changes in head

Changes in seepage flux and head (groundwater level) from groundwater drainage into a surface coal mine can be predicted by a technique that considers drainage from the unsaturated zone. The user applies site-specific data to precalculated head and seepage-flux profiles. Groundwater flow through hypothetical aquifer cross sections was simulated using the U.S. Geological Survey finite-difference model, VS2D, which considers variably saturated two-dimensional flow. Conceptual models considered were (1) drainage to a first cut, and (2) drainage to multiple cuts, which includes drainage effects of an area surface mine. Dimensionless head and seepage flux profiles from 246 simulations are presented. Step-by-step instructions and examples are presented. Users are required to know aquifer characteristics and to estimate size and timing of the mine operation at a proposed site. Calculated groundwater drainage to the mine is from one excavated face only. First cut considers confined and unconfined aquifers of a wide range of permeabilities; multiple cuts considers unconfined aquifers of higher permeabilities only. The technique, developed for Illinois coal-mining regions that use area surface mining and evaluated with an actual field example, will be useful in assessing potential hydrologic impacts of mining. Application is limited to hydrogeologic settings and mine operations similar to those considered. Fracture flow, recharge, and leakage are nor considered. (USGS)

Water-Resources Investigations Report↗

USGS2003

Estimating Ground-Water Recharge from Precipitation on Whidbey and Camano Islands, Island County, Washington, Water Years 1998 and 1999

Ground-water recharge from precipitation to unconsolidated deposits on Whidbey and Camano Islands, Washington, was estimated for water years 1998-99 using a near-surface water-balance method and a chloride mass-balance method. A daily near-surface water-balance method, the Deep Percolation Model (DPM), was used to simulate water budgets for October 1, 1997 through September 30, 1999 (water years 1998-99) for six small drainage basins?four on Whidbey Island and two on Camano Island. Adjusted parameters from the DPM for each small basin were then used in island-wide DPM simulations. A spatial distribution of annual recharge was simulated for each island, with island averages of 5.71 inches per year for Whidbey Island and 5.98 inches per year for Camano Island. The spatial distribution of simulated annual recharge for each island reflects variations in precipitation amounts and the distribution of surficial materials. DPM results indicate that recharge generally is higher in areas underlain by coarse-grained deposits (outwash) than in areas underlain by fine-grained deposits (till). A chloride mass-balance method was used to estimate combined recharge to unconsolidated deposits on Whidbey and Camano Islands. The average combined recharge for Whidbey and Camano Islands estimated by this method was 2.00 inches per year. The range of chloride concentrations in ground-water samples from selected wells indicates that the average recharge to unconsolidated deposits ranges from 0.78 to 7.81 inches per year. Sources of chloride in ground water other than from the atmosphere would cause recharge estimated by the chloride mass-balance method to be less than the actual recharge, therefore, these estimates may represent lower limits.

Water-Resources Investigations Report↗

USGS2001

Simulation of flow in the upper North Coast Limestone Aquifer, Manati-Vega Baja area, Puerto Rico

A two-dimensional computer ground-water model was constructed of the Manati-Vega Baja area to improve the understanding of the unconfined upper aquifer within the North Coast Province of Puerto Rico. The modeled area covers approximately 79 square miles within the municipios of Manati and Vega Baja and small portions of Vega Alta and Barceloneta. Steady-state two-dimensional ground-water simulations were correlated to conditions prior to construction of the Laguna Tortuguero outlet channel in 1940 and calibrated to the observed potentiometric surface in March 1995. At the regional scale, the unconfined Upper North Coast Limestone aquifer is a diffuse ground-water flow system through the Aguada and Aymamon limestone units. The calibrated model input parameters for aquifer recharge varied from 2 inches per year in coastal areas to 18 inches per year in the upland areas south of Manati and Vega Baja. The calibrated transmissivity values ranged from less than 500 feet squared per day in the upland areas near the southern boundary to 70,000 feet squared per day in the areas west of Vega Baja. Increased ground-water withdrawals from 1.0 cubic foot per second for 1940 conditions to 26.3 cubic feet per second in 1995, has reduced the natural ground-water discharge to springs and wetland areas, and induced additional recharge from the rivers. The most important regional drainage feature is Laguna Tortuguero, which is the major ground-water discharge body for the upper aquifer, and has a drainage area of approximately 17 square miles. The discharge to the sea from Laguna Tortuguero through the outlet channel has been measured on a bi-monthly basis since 1974. The outflow represents a combination of ground- and surface-water discharge over the drainage area. Hydrologic conditions, prior to construction of the Laguna Tortuguero outlet channel in 1943, can be considered natural conditions with minimal ground-water pumpage (1.0 cubic foot per second), and heads in the lagoon were 2.4 feet higher. The model was calibrated to March 1995 conditions during a dry period of minimal aquifer recharge and relatively constant water levels in the upper aquifer. For the steady-state 1995 model simulation, however, ground-water pumpage had been increased to 26.3 cubic foot per second, due to increased demand for public water supply, the heads at 0.9 feet, and the outflow to the sea at Laguna Tortuguero had been lowered considerably. Simulated ground-water inflow for 1940 hydrologic conditions included 35.9 cubic feet per second from areal recharge, contributions from streamflow along the southern boundary of 1.6 cubic feet per second, and streamflow infiltration to the upper aquifer of 4.2 cubic feet per second. Simulated ground-water outflow for 1940 hydrologic conditions are discharge to springs of 17.4 cubic feet per second, total ground-water withdrawals of 1.0 cubic feet per second, and aquifer contribution to streamflow or wetland areas of 23.4 cubic feet per second. Simulated ground-water inflow for hydrologic conditions of March 1995 include d contributions from streamflow along the southern boundary of 1.6 cubic feet per second, areal recharge of 35.9 cubic feet per second, and streamflow infiltration to the upper aquifer of 11 cubic feet per second. Simulated ground-water outflow for hydrologic conditions of March 1995 are ground-water withdrawals of 26.3 cubic feet per second, discharge from springs of 7.3 cubic feet per second, and aquifer contribution to streamflow or wetland areas of 14 .9 cubic feet per second. The overall ground-water budget increased from 41.8 cubic feet per second for 1940 conditions to 48.6 cubic feet per second for the hydrologic conditions of March 1995. The increase in ground-water budget is a direct result of increased ground-water withdrawals, which induced greater streamflow infiltration. Simulated ground-water flux to Laguna Tortuguero for 1940 conditions was 11 cubic feet per second, which drop

Water-Resources Investigations Report↗

USGS1989

Simulation of steady-state ground water and spring flow in the upper Floridan aquifer of coastal Citrus and Hernando Counties, Florida

A digital groundwater flow model was developed to approximate steady-state predevelopment flow conditions in the Upper Floridan aquifer of coastal west-central Florida. The aquifer is the major source of water and natural spring flow in the area. The aquifer was simulated as a one-layer system with constant vertical recharge and discharge rates. Calibrated transmissivities ranged from 8,640 sq ft/day in the northern part of the area to nearly 13,000,000 sq ft/day near large springs. Calibrated inflows were about 2,708 cu ft/sec; of this, about 2,565 cu ft/sec discharged as natural spring flow and 137 cu ft/sec discharged as upward leakage along the coast. The model was used to show how the system might respond to large manmade stresses. Withdrawal of 116 cu ft/sec from a hypothetical regional well field resulted in potentiometric-surface drawdowns ranging from 0.1 to 1.7 ft and declines of generally less than 0.2 ft along the coast. Total spring flow decreased 5%, and the effect on individual springs varied from 0.1 to 8.0%. Withdrawal of 62 cu ft/sec from the 4-sq-mi node at each spring resulted in six of seven springs to the south of the Chassahowitzka River contributing 50% of their flow to pumpage. Springs located north of the Chassahowitzka River contributed as much as 18% of their flow to pumpage. (USGS)

Water-Resources Investigations Report↗

USGS2009

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↗

USGS1990

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↗

USGS1977

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-quality 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.

California↗

USGS1996

Hydrogeology and simulation of ground-water flow in the alluvial aquifer at Louisville, Kentucky

The alluvial aquifer at Louisville, Ky., lies in a valley eroded by glacial meltwater that was later partly filled with outwash sand and gravel deposits. The aquifer is primarily unconfined, and the direction of flow is from the adjacent limestone and shale valley wall toward the Ohio River and major pumping centers. Pumpage and water-level data indicate that the alluvial aquifer was in a steady-state condition in November 1962 and again in November 1983. Between these two dates, water-level data indicate a general rise in the water table. A two-dimensional finite-element ground-water-flow model of the alluvial aquifer was calibrated for both the steady-state and the transient-state period of 1962-83. The year 1962 represented a period in time when pumping was nearly three times that in 1983. The simulated steady-state water budget for 1962 indicated that of the total recharge to the aquifer of 5.19 million feet per day, 37.2 percent was flow from the river to pumped wells, 28.3 percent was recharge from rainfall, 19.7 percent was flow across the eastern valley wall, and 14.8 percent was upward flow from the bedrock. Discharge from the aquifer was to wells (68.9 percent) and to the Ohio River (31.1 percent). The simulated steady-state water budget for 1983 indicated that of the total recharge to the aquifer of 4.11 million feet per day, 42.6 percent was recharge from rainfall, 18.2 percent was flow across the eastern valley wall, 17.8 percent was flow from the river to pumped wells, 15.6 percent was upward flow from the bedrock, and 5.8 percent was flow from septic systems. The transient simulation resulted in an acceptable match between measured and simulated hydrographs. This gave additional confidence to the model calibration, choice of boundary conditions, and published values of specific yield. Both steady-state and transient-state models demonstrated that the main source of water needed to meet increased pumping requirements was induced flow from the Ohio River.

Water-Resources Investigations Report↗

USGS1999

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↗

USGS2003

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↗