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 685 records · Page 38Linked to original sources

Conceptual Model of Hydrologic and Thermal Conditions of the Eastbank Aquifer System near Rocky Reach Dam, Douglas County, Washington

The Lower and Combined Aquifers of the Eastbank Aquifer system, located in a river-terrace deposit along the Columbia River near Rocky Reach Dam, Washington, are primarily recharged by the Columbia River and provide water to the Eastbank Hatchery and the regional water system servicing the cities of Wenatchee, East Wenatchee, and parts of unincorporated Chelan and Douglas Counties. In 2006, mean annual pumpage from the aquifers by the hatchery and regional water system was about 43 and 16 cubic feet per second, respectively. Reportedly, temperatures of ground water pumped by the hatchery have been increasing, thereby making water potentially too warm for salmonid fish production. An evaluation of hourly ground-water and river temperatures from January 1991 through August 2007 indicates increasing interannual trends in temperatures in most of the Lower and Combined Aquifers from 1999 through 2006 that correspond to increasing trends in the annual mean and annual maximum river temperatures during the same period of 0.07 and 0.17?C per year, respectively. There were no trends in the annual minimum river temperatures from 1999 through 2006, and there were no trends in the annual minimum, mean, and maximum river temperatures from 1991 through 1998 and from 1991 through 2007. Increases in river temperatures from 1999 through 2006 are within the natural variability of the river temperatures. Most of the Lower and Combined Aquifers reached thermal equilibrium?defined by constant time lags between changes in river temperatures and subsequent changes in ground-water temperatures?during 1991?98. The only exceptions are the Combined Aquifer north of the well field of the regional water system, which had not reached thermal equilibrium by 2006, and the Lower Aquifer west of the well fields of the hatchery and the regional water system, which reached thermal equilibrium prior to 1991. Because most of the Lower and Combined Aquifers were in thermal equilibrium from 1999 through 2006 and seasonal pumpage patterns were relatively stable, reported trends of increasing temperatures of water pumped by the hatchery well field are most likely explained by increasing trends in river temperatures. Most of the water pumped by the hatchery well field recharges in an area west to southwest of the well field about 2 months prior to the time it is pumped from the aquifer. The northern extent of the hatchery well field may pump some colder water from a bedrock depression to the north and west of the well field. The conceptual model of hydrologic and thermal conditions is supported by analyses of historical water temperatures, water-level data collected on July 18, 2007, and dissolved-constituent and bacterial concentrations in samples collected on August 20?22, 2007.

Scientific Investigations Report↗

The Cerrillos Uplift, the La Bajada Constriction, and hydrogeologic framework of the Santo Domingo Basin, Rio Grande Rift, New Mexico

The geologic, geophysical, and hydrogeologic properties of the La Bajada constriction and Santo Domingo Basin, northern New Mexico, result from tectonic and volcanic processes of the late Tertiary and Quaternary Rio Grande rift. An integrated geologic and geophysical assessment in the La Bajada constriction allows development of a geologic framework that can provide input for regional ground-water flow models. These models then can provide better estimates of future water supplies in a region that largely subsists on aquifers in Rio Grande rift basins. The combination of surface geologic investigations (stratigraphic and structural studies; chapters A, B, C, and E), airborne geophysics (aeromagnetic and time-domain electromagnetic surveys; chapters D and F), ground geophysical measurements (gravity and magnetotelluric surveys; chapters D and F), and data from the few wells in the area (chapter G) provides new constraints on the hydrogeologic framework of this area. Summary results of our investigations are synthesized in chapter G. Through-going aquifers consisting of ancestral Rio Grande axial-river sand and gravel and of coarse western-piedmont gravel form the predominant ground-water pathways through the partly buried structural trough defining the La Bajada constriction between Española and Santo Domingo Basins. Thick, clay-rich Cretaceous marine shales of low hydraulic conductivity form a pervasive regional confining unit within the Cerrillos uplift on the southeast flank of the constriction. Numerous, dominantly north-northwest-striking, intrabasin faults that project part way across the La Bajada constriction create a matrix of laterally and vertically variable hydrogeologic compartments that locally partition and deflect ground-water flow parallel to faults.

New Mexico↗

Reconnaissance Assessment of the Potential for Roadside Dry Wells to Affect Water Quality on the Island of Hawai'i

The County of Hawai'i Department of Public Works (DPW) uses dry wells to dispose of stormwater runoff from roads. Recently, concern has been raised that water entering the dry wells may transport contaminants to groundwater and affect the quality of receiving waters. The DPW operates 2,052 dry wells. Compiling an inventory of these dry wells and sorting it on the basis of presence or absence of urbanization in the drainage area, distance between the bottom of the dry well and the water table, and proximity to receiving waters helps identify the dry wells having greatest potential to affect the quality of receiving waters so that future studies or mitigation efforts can focus on a smaller number of dry wells. The drainage areas of some DPW dry wells encompass urbanized areas, which could be a source of contaminants. Some dry wells penetrate close to or through the water table, eliminating or substantially reducing opportunities for contaminant attenuation between the ground surface and water table. Dry wells that have drainage areas that encompass urbanization, penetrate to near the water table, and are near the coast have the highest potential to affect the quality of coastal waters (this study did not consider specific sections of coastline that may be of greater concern than others). Some DPW dry wells, including a few that have drainage areas that encompass urbanization, lie within the areas contributing recharge (ACR) to drinking-water wells. Numerical groundwater modeling studies by previous investigators indicate that water infiltrating those dry wells could eventually be pumped at drinking-water wells. Dry wells that have a high potential for affecting coastal receiving waters or drinking-water wells can be the focus of studies to further understand the effect of the dry wells on the quality of receiving waters. Possible study approaches include sampling for contaminants at the dry well and receiving water, injecting and monitoring the movement of tracers, and numerical modeling. To fully assess whether dry wells actually pose a significant contamination threat to receiving waters, results from modeling or monitoring must be compared to limits for contaminant concentration at receiving waters. These limits are usually established by the agencies tasked with protecting those waters.

Scientific Investigations Report↗

Simulation of Ground-Water Flow in the Irwin Basin Aquifer System, Fort Irwin National Training Center, California

Ground-water pumping in the Irwin Basin at Fort Irwin National Training Center, California resulted in water-level declines of about 30 feet from 1941 to 1996. Since 1992, artificial recharge from wastewater-effluent infiltration and irrigation-return flow has stabilized water levels, but there is concern that future water demands associated with expansion of the base may cause a resumption of water-level declines. To address these concerns, a ground-water flow model of the Irwin Basin was developed to help better understand the aquifer system, assess the long-term availability and quality of ground water, and evaluate ground-water conditions owing to current pumping and to plan for future water needs at the base. Historical data show that ground-water-level declines in the Irwin Basin between 1941 and 1996, caused the formation of a pumping depression near the pumped wells, and that recharge from the wastewater-treatment facility and disposal area caused the formation of a recharge mound. There have been two periods of water-level recovery in the Irwin Basin since the development of ground water in this basin; these periods coincide with a period of decreased pumpage from the basin and a period of increased recharge of water imported from the Bicycle Basin beginning in 1967 and from the Langford Basin beginning in 1992. Since 1992, artificial recharge has exceeded pumpage in the Irwin Basin and has stabilized water-level declines. A two-layer ground-water flow model was developed to help better understand the aquifer system, assess the long-term availability and quality of ground water, and evaluate ground-water conditions owing to current pumping and to plan for future water needs at the base. Boundary conditions, hydraulic conductivity, altitude of the bottom of the layers, vertical conductance, storage coefficient, recharge, and discharge were determined using existing geohydrologic data. Rates and distribution of recharge and discharge were determined from existing data and estimated when unavailable. Results of predictive simulations indicate that in 50 years, if artificial recharge continues to exceed pumpage in Irwin Basin, water levels could rise as much as 65 feet beneath the pumping depression, and as much as 10 feet in the wastewater-treatment facility and disposal area. Particle-tracking simulations were used to determine the pathlines and the traveltimes of water high in dissolved solids into the main pumping area. The pathlines of particles from two areas with high dissolved-solids concentrations show that in 50 years water from these areas almost reaches the nearest pumped well.

California↗

Spatial and temporal variations in streamflow, dissolved solids, nutrients, and suspended sediment in the Rio Grande Valley study unit, Colorado, New Mexico, and Texas, 1993–95

Streamflow and water quality vary spatially and temporally in the Rio Grande from Del Norte, Colorado, to El Paso, Texas. The variations in streamflow and in concentrations of selected waterquality constituents—dissolved solids, dissolved nitrite plus nitrate as nitrogen, total phosphorus, and suspended sediment—are described in this report. A multivariate linear regression model, ESTIMATOR2000, was used to estimate loads for selected constituents. Streamflow decreases in the downstream direction throughout most of the basin because outflows (due to agricultural use, leakage to ground water, and evapotranspiration) are greater than inflows. Streamflow increases between Rio Grande above the mouth of Trinchera Creek, near Lasauses, Colorado, to Rio Grande at Otowi Bridge, near San Ildefonso, New Mexico, because ground-water and tributary inflow are greater than outflow. Concentrations of dissolved solids, dissolved nitrite plus nitrate, total phosphorus, and suspended sediment generally increase in the downstream direction. Concentrations of dissolved solids, dissolved nitrite plus nitrate, and total phosphorus decrease between Rio Grande above the mouth of Trinchera Creek, near Lasauses, Colorado, and Rio Grande at Otowi Bridge, near San Ildefonso, New Mexico, because of dilution by tributary inflow. Concentrations of dissolved nitrite plus nitrate, total phosphorus, and suspended sediment decrease between Rio Grande Floodway at San Marcial, New Mexico, and Rio Grande below Leasburg Dam, near Leasburg, New Mexico, because of reservoir effects (nutrient uptake and settling of sediment). Several instances of decreasing streamflow and increasing loads indicate the presence of inflows with large constituent concentrations (relative to those of the Rio Grande immediately upstream from that inflow); this occurs (1) between Rio Grande near Del Norte, Colorado, and Rio Grande above the mouth of Trinchera Creek, near Lasauses, Colorado, for dissolved solids, (2) between Rio Grande at Otowi Bridge, near San Ildefonso, New Mexico, and Rio Grande Floodway at San Marcial, NewMexico, for all constituents, and (3) between Rio Grande below Leasburg Dam, near Leasburg, New Mexico, and Rio Grande at El Paso, Texas, for all constituents. Streamflow increases along every reach of the Rio Grande between the streamflow-gaging station Rio Grande above the mouth of Trinchera Creek, near Lasauses, Colorado, and the station Rio Grande at Otowi Bridge, near San Ildefonso, NewMexico. These increases in streamflow result in increases in the loads of dissolved solids, total phosphorus, and suspended sediment regardless of changes in concentrations.

Water-Resources Investigations Report↗

Hydrology and water quality of Whitewater and Rice lakes in southeastern Wisconsin, 1990-91

The hydrology and water quality of Whitewater and Rice Lakes was studied by the U.S. Geological Survey during November 15, 1990November 14, 1991, in cooperation with the Whitewater-Rice Lakes Management District, Walworth County, Wis. Whitewater and Rice Lakes are small, shallow lakes; surface areas are 697 and 162 acres and mean depths are 8.4 feet and 5.8 feet, respectively. Although both lakes have surface outlets, water levels were below the dam crests during the study, and no water left the lake through the outlets. The drainage basin of Whitewater Lake is 10.9 square miles and that of Rice Lake is 11.8 square miles; but, because of large amounts of depressional areas, only 1.4 square miles and 0.2 square mile, respectively, contribute surface runoff to the lakes. Whitewater Lake is an artificial lake created in 1947 by the damming of three smaller lakes. Rice Lake is an artificial lake created in 1954 by the damming of Whitewater Creek, which drains Whitewater Lake. Maintaining the lake levels at the elevations of their dam crests has been difficult since the lakes were created. For most years, water levels were below the lakes' dam crests. Ground water, precipitation, and evaporation are important components in the hydrologic budgets of the lakes. For Whitewater Lake, ground water was the dominant source of water, accounting for 57 percent of the inflow budget; precipitation accounted for 26 percent. Ground water also dominated the outflow, accounting for 81 percent of the outflow budget. The remaining 19 percent of the outflow budget was evaporation. For Rice Lake, precipitation was the dominant source of water, accounting for 88 percent of the inflow budget; ground water accounted for 8 percent. Evaporation dominated the outflow budget, at 70 percent, whereas ground water accounted for 30 percent. The external phosphorus budget for Whitewater Lake showed that shoreline drainage was the largest source of phosphorus to the lake42 percent of the total input of 558 pounds. Other sources of phosphorus were septic systems, 19 percent of the total; precipitation, 18 percent; a spring inlet at base flow, 13 percent; and ground water, 8 percent. The external phosphorus budget for Rice Lake showed that shoreline drainage also was the largest source of phosphorus to the lake-59 percent of the total input of 63 pounds; other sources were precipitation, 38 percent of the total; and ground water, 3 percent. Application of Vollenweider's phosphorus loading model fairly accurately predicted the lakes' spring turnover phosphorus concentrations and suggested that the external loading of phosphorus would result in mesotrophic to eutrophic conditions for Whitewater Lake and mesotrophic conditions for Rice Lake. Dillon and Rigler's model further suggested additional phosphorus from internal recycling was required to result in the high chlorophyll-a concentrations experienced in both systems during summer. Internal recycling of phosphorus in addition to external loading seems to also cause waterquality problems in both lakes. The amount of phosphorus recycled from the lake sediments was estimated from a mass-balance approach for April 1-November 14, 1991. For Whitewater Lake, the internal load of 582 pounds was slightly greater than the annual external load of 558 pounds. For Rice Lake, the internal load of 295 pounds far exceeded the annual external load of 63 pounds.

Wisconsin↗

Simulated changes in the sources of ground water for public-supply wells, ponds, streams, and coastal areas on western Cape Cod, Massachusetts

The U.S. Geological Survey, in cooperation with the Air Force Center for Environmental Excellence, developed a three-dimensional ground-water-flow model to simulate sources of water to public-supply wells, ponds, streams, and coastal areas for average and changing pumping and recharge conditions on western Cape Cod, Massachusetts. Understanding the changing sources of water that flows to these discharge areas under changing hydrologic conditions and stresses is needed to help water-resource managers protect the quality of the drinking water and the water resources of western Cape Cod. The interconnection of hydrologic features through the ground-water-flow system on Cape Cod complicates the determination of the sources of water for changing pumping and recharge conditions. When pumping and recharge conditions change, the numerical model can help evaluate changes in the recharge areas to pumped wells, ponds, streams, and coastal areas. The steady-state recharge area and time-of-travel zones to a pumping well represent the source of water and time of travel from the water table to the well once all of the water discharging to the well entered the aquifer after pumping began. Initially, the water that discharges to a pumped well over a specified period consists of a blend of waters that entered the aquifer over time. When pumping begins, a deep well may capture water that entered the aquifer many years earlier at a great distance from the well; therefore, determining the source of water discharging to a well for specified pumping times may be more important for protecting public-water supplies than determining only the steady-state recharge area and time-of-travel zones to the well.

Massachusetts↗

Ground-water-withdrawal and water-level data used to simulate regional flow in the major coastal plain aquifers of New Jersey

The report documents groundwater withdrawal and level data for the major aquifers of the New Jersey Coastal Plain. These data have been used in the Regional Aquifer System Analysis flow model of 10 major aquifers in the New Jersey Coastal Plain. Groundwater withdrawal data from more than 1,400 wells are tabulated by year for 1956 through 1980. Total yearly withdrawals for 171 public supply and industrial purveyors for 1918 through 1955 are documented. Best estimates are given for percentages of total withdrawals by aquifer for the 1918 through 1955 data. Groundwater withdrawal data is also tabulated yearly from 1918 through 1980 by county for each aquifer. Long-term hydrographs for 89 Coastal Plain wells provide information on water level fluctuations over time. Potentiometric maps showing earliest recorded water levels in Coastal Plain aquifers delineate prepumping conditions. Groundwater withdrawals from the Coastal Plain have increased from less than 50 million gal/day in 1918 to more than 350 million gal/day in 1980. In response to these increased withdrawals, groundwater levels have been steadily declining. In some areas where wells once flowed freely, water levels are now 100 to 200 ft below land surface. (USGS)

Water-Resources Investigations Report↗

Geohydrology and simulation of ground-water flow in the aquifer system near Calvert City, Kentucky

The U.S. Geological Survey, in cooperation with the Kentucky Natural Resources and Environmental Protection Cabinet, constructed a two-dimensional, steady-state ground-water-flow model to estimate hydraulic properties, contributing areas to discharge boundaries, and the average linear velocity at selected locations in an aquifer system near Calvert City, Ky. Nonlinear regression was used to estimate values of model parameters and the reliability of the parameter estimates. The regression minimizes the weighted difference between observed and calculated hydraulic heads and rates of flow. The calibrated model generally was better than alternative models considered, and although adding transmissive faults in the bedrock produced a slightly better model, fault transmissivity was not estimated reliably. The average transmissivity of the aquifer was 20,000 feet squared per day. Recharge to two outcrop areas, the McNairy Formation of Cretaceous age and the alluvium of Quaternary age, were 0.00269 feet per day (11.8 inches per year) and 0.000484 feet per day (2.1 inches per year), respectively. Contributing areas to wells at the Calvert City Water Company in 1992 did not include the Calvert City Industrial Complex. Since completing the fieldwork for this study in 1992, the Calvert City Water Company discontinued use of their wells and began withdrawing water from new wells that were located 4.5 miles east-southeast of the previous location; the contributing area moved farther from the industrial complex. The extent of the alluvium contributing water to wells was limited by the overlying lacustrine deposits. The average linear ground-water velocity at the industrial complex ranged from 0.90 feet per day to 4.47 feet per day with a mean of 1.98 feet per day.

Water-Resources Investigations Report↗

Geophysical Data from Spring Valley to Delamar Valley, East-Central Nevada

Cenozoic basins in eastern Nevada and western Utah constitute major ground-water recharge areas in the eastern part of the Great Basin and these were investigated to characterize the geologic framework of the region. Prior to these investigations, regional gravity coverage was variable over the region, adequate in some areas and very sparse in others. Cooperative studies described herein have established 1,447 new gravity stations in the region, providing a detailed description of density variations in the middle to upper crust. All previously available gravity data for the study area were evaluated to determine their reliability, prior to combining with our recent results and calculating an up-to-date isostatic residual gravity map of the area. A gravity inversion method was used to calculate depths to pre-Cenozoic basement rock and estimates of maximum alluvial/volcanic fill in the major valleys of the study area. The enhanced gravity coverage and the incorporation of lithologic information from several deep oil and gas wells yields a much improved view of subsurface shapes of these basins and provides insights useful for the development of hydrogeologic models for the region.

Open-File Report↗

Mineral and Vegetation Maps of the Bodie Hills, Sweetwater Mountains, and Wassuk Range, California/Nevada, Generated from ASTER Satellite Data

Multispectral remote sensing data acquired by the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) were analyzed to identify and map minerals, vegetation groups, and volatiles (water and snow) in support of geologic studies of the Bodie Hills, Sweetwater Mountains, and Wassuk Range, California/Nevada. Digital mineral and vegetation mapping results are presented in both portable document format (PDF) and ERDAS Imagine format (.img). The ERDAS-format files are suitable for integration with other geospatial data in Geographic Information Systems (GIS) such as ArcGIS. The ERDAS files showing occurrence of 1) iron-bearing minerals, vegetation, and water, and 2) clay, sulfate, mica, carbonate, Mg-OH, and hydrous quartz minerals have been attributed according to identified material, so that the material detected in a pixel can be queried with the interactive attribute identification tools of GIS and image processing software packages (for example, the Identify Tool of ArcMap and the Inquire Cursor Tool of ERDAS Imagine). All raster data have been orthorectified to the Universal Transverse Mercator (UTM) projection using a projective transform with ground-control points selected from orthorectified Landsat Thematic Mapper data and a digital elevation model from the U.S. Geological Survey (USGS) National Elevation Dataset (1/3 arc second, 10 m resolution). Metadata compliant with Federal Geographic Data Committee (FGDC) standards for all ERDAS-format files have been included, and contain important information regarding geographic coordinate systems, attributes, and cross-references. Documentation regarding spectral analysis methodologies employed to make the maps is included in these cross-references.

Scientific Investigations Map↗

Salmonids, stream temperatures, and solar loading--modeling the shade provided to the Klamath River by vegetation and geomorphology

The U.S. Geological Survey is studying approaches to characterize the thermal regulation of water and the dynamics of cold water refugia. High temperatures have physiological impacts on anadromous fish species. Factors affecting the presence, variability, and quality of thermal refugia are known, such as riverine and watershed processes, hyporheic flows, deep pools and bathymetric factors, thermal stratification of reservoirs, and other broader climatic considerations. This research develops a conceptual model and methodological techniques to quantify the change in solar insolation load to the Klamath River caused by riparian and floodplain vegetation, the morphology of the river, and the orientation and topographic characteristics of its watersheds. Using multiple scales of input data from digital elevation models and airborne light detection and ranging (LiDAR) derivatives, different analysis methods yielded three different model results. These models are correlated with thermal infrared imagery for ground-truth information at the focal confluence with the Scott River. Results from nonparametric correlation tests, geostatistical cross-covariograms, and cross-correlograms indicate that statistical relationships between the insolation models and the thermal infrared imagery exist and are significant. Furthermore, the use of geostatistics provides insights to the spatial structure of the relationships that would not be apparent otherwise. To incorporate a more complete representation of the temperature dynamics in the river system, other variables including the factors mentioned above, and their influence on solar loading, are discussed. With similar datasets, these methods could be applied to any river in the United States—especially those listed as temperature impaired under Section 303(d) of the Clean Water Act—or international riverine systems. Considering the importance of thermal refugia for aquatic species, these methods can help investigate opportunities for riparian restoration, identify problematic reaches unlikely to provide good habitat, and simulate changes to solar loading estimates from alternative landscape configurations.

California↗

Simulation of ground-water flow at Anchorage, Alaska, 1955-83

The groundwater system at Anchorage, Alaska was analyzed by using a two-layer three-dimensional mathematical model. By use of existing data, both nonpumping and pumping steady-state conditions and transient conditions were simulated. Under steady-state conditions, calculated directions of groundwater flow were similar to observed flow patterns, and calculated stream discharges generally were within 10% of observed values. However, in many parts of the modeled area computed hydraulic head values were more than 20 ft higher or lower than observed values. Hydraulic conductivity and transmissivity are the most sensitive hydraulic parameters under steady-state conditions. Under steady-state conditions, a pumping rate of 18.8 Mgal/d lowers heads in the confined aquifer by as much as 30 ft, but reduces streamflow by less than 5%. Transient conditions show that drawndowns due to withdrawals by production wells follow similar patterns of nearby observation wells. On the basis of analytical techniques, the confining layer does not appear to hold significant quantities of water. (USGS)

Water-Resources Investigations Report↗

Hydrogeologic framework and simulation of ground-water flow and travel time in the shallow aquifer system in the area of Naval Support Activity Memphis, Millington, Tennessee

Naval Support Activity (NSA) Memphis is a Department of the Navy facility located at the City of Millington, Tennessee, about 5 miles north of Memphis. Contaminants have been detected in surface-water, sediment, and ground-water samples collected at the facility. As part of the Installation Restoration Program, the Navy is considering remedial-action options to prevent or lessen the effect of ground-water contamination at the facility and to control the movement and discharge of contaminants. A numerical model of the ground-water-flow system in the area of NSA Memphis was constructed and calibrated so that quantifiable estimates could be made of ground-water-flow rates, direction, and time-of-travel. The sediments beneath NSA Memphis, to a depth of about 200 feet, form a shallow aquifer system. From youngest to oldest, the stratigraphic units that form the shallow aquifer system are alluvium, loess, fluvial deposits, and the Cockfield and Cook Mountain Formations. The shallow aquifer system is organized into five hydrogeologic units: (1) a confining unit composed of the relatively low permeability sediments of the upper alluvium and the loess; (2) the A1 aquifer comprising sand and gravel of the lower alluvium and the fluvial deposits, and sand lenses in the upper part of the preserved section of the Cockfield Formation; (3) a confining unit composed of clay and silt within the upper part of the Cockfield Formation; (4) the Cockfield aquifer comprising sand lenses within the lower part of the preserved section of the Cockfield Formation; and (5) a confining unit formed by low permeability sediments of the Cook Mountain Formation that composes the upper confining unit for the Memphis aquifer. Thicknesses of individual units vary considerably across the facility. Structural and depositional features that affect the occurrence of ground water in the shallow aquifer system include faulting, an erosional scarp, and 'windows' in the confining units. Underlying the shallow aquifer system is the Memphis aquifer, the primary source of water for NSA Memphis and the City of Memphis, Tennessee. Analyses of sediment cores, aquifer and well specific-capacity tests, and numerical modeling were used to estimate the hydraulic characteristics of units of the shallow aquifer system. The vertical hydraulic conductivity of core samples of the alluvium-loess confining unit ranged from about 8.5 x 10-5 to 1.6 x 10-2 feet per day, and the total porosity of the samples ranged from about 35 to 48 percent. The results of the aquifer test were used to estimate a horizontal hydraulic conductivity of about 5 feet per day for the alluvial-fluvial deposits aquifer. The total porosity of core samples of the alluvial-fluvial deposits aquifer ranged from about 22 to 39 percent. The vertical hydraulic conductivity of core samples of the Cockfield confining unit ranged from about 4.5 x 10-5 to 2.5 x 10-3 feet per day, and the total porosity ranged from about 41 to 55 percent. Well specific-capacity tests indicate that the horizontal hydraulic conductivity of sand units that compose the Cockfield aquifer range from about 0.5 to 3 feet per day. The vertical hydraulic conductivity of core samples of the Cook Mountain confining unit ranged from about 5.0 x 10-6 to 9.9 x 10-4 feet per day. Total porosity of core samples of the Cook Mountain confining unit ranged from about 30 to 42 percent. Ground-water flow and time-of-travel in the shallow aquifer system were simulated using the MODFLOW finite-difference model and the -particle-tracking program MODPATH. A three-layer, steady-state model of the shallow aquifer system was constructed and calibrated to the potentiometric surface of the A1 aquifer. Results of numerical modeling support the proposed conceptual hydrogeologic model of the shallow aquifer system. Ground-water time-of-travel in the A1 aquifer was simulated using an assumed effective porosity of 25 percent. Typical ground-water-flow velocities were on the order of 15 to 25 feet per year in the layer representing the Al aquifer in the model. The average residence time of particles seeded in this layer was about 800 years. Ground-water travel times were simulated at three sites within the Al aquifer: (1) the former N-6 hangar location, (2) the "grassy" area near Solid Waste Management Unit 7, and (3) at Solid Waste Management Unit 2. Results indicate close agreement between the particle-tracking simulations and the measured extent of contaminant plumes at the former N-6 hangar area and the grassy area near Solid Waste Management Unit 7. Based on the results of particle-tracking analyses of ground-water flow and the estimated locations of contaminant plumes at these two sites, the potential for contaminants to reach the Memphis aquifer in the next 100 years is negligible. However, particle-tracking analysis of ground-water flow at Solid Waste Management Unit 2 suggests that the time-of-travel of contaminants to Big Creek Drainage Canal could be less than 30 years.

Tennessee↗

Simulation of climate-change effects on streamflow, lake water budgets, and stream temperature using GSFLOW and SNTEMP, Trout Lake Watershed, Wisconsin

Although groundwater and surface water are considered a single resource, historically hydrologic simulations have not accounted for feedback loops between the groundwater system and other hydrologic processes. These feedbacks include timing and rates of evapotranspiration, surface runoff, soil-zone flow, and interactions with the groundwater system. Simulations that iteratively couple the surface-water and groundwater systems, however, are characterized by long run times and calibration challenges. In this study, calibrated, uncoupled transient surface-water and steady-state groundwater models were used to construct one coupled transient groundwater/surface-water model for the Trout Lake Watershed in north-central Wisconsin, USA. The computer code GSFLOW (Ground-water/Surface-water FLOW) was used to simulate the coupled hydrologic system; a surface-water model represented hydrologic processes in the atmosphere, at land surface, and within the soil-zone, and a groundwater-flow model represented the unsaturated zone, saturated zone, stream, and lake budgets. The coupled GSFLOW model was calibrated by using heads, streamflows, lake levels, actual evapotranspiration rates, solar radiation, and snowpack measurements collected during water years 1998–2007; calibration was performed by using advanced features present in the PEST parameter estimation software suite. Simulated streamflows from the calibrated GSFLOW model and other basin characteristics were used as input to the one-dimensional SNTEMP (Stream-Network TEMPerature) model to simulate daily stream temperature in selected tributaries in the watershed. The temperature model was calibrated to high-resolution stream temperature time-series data measured in 2002. The calibrated GSFLOW and SNTEMP models were then used to simulate effects of potential climate change for the period extending to the year 2100. An ensemble of climate models and emission scenarios was evaluated. Downscaled climate drivers for the period 2010–2100 showed increases in maximum and minimum temperature over the scenario period. Scenarios of future precipitation did not show a monotonic trend like temperature. Uncertainty in the climate drivers increased over time for both temperature and precipitation. Separate calibration of the uncoupled groundwater and surface-water models did not provide a representative initial parameter set for coupled model calibration. A sequentially linked calibration, in which the uncoupled models were linked by means of utility software, provided a starting parameter set suitable for coupled model calibration. Even with sequentially linked calibration, however, transmissivity of the lower part of the aquifer required further adjustment during coupled model calibration to attain reasonable parameter values for evaporation rates off a small seepage lake (a lake with no appreciable surface-water outlets) with a long history of study. The resulting coupled model was well calibrated to most types of observed time-series data used for calibration. Daily stream temperatures measured during 2002 were successfully simulated with SNTEMP; the model fit was acceptable for a range of groundwater inflow rates into the streams. Forecasts of potential climate change scenarios showed growing season length increasing by weeks, and both potential and actual evapotranspiration rates increasing appreciably, in response to increasing air temperature. Simulated actual evapotranspiration rates increased less than simulated potential evapotranspiration rates as a result of water limitation in the root zone during the summer high-evapotranspiration period. The hydrologic-system response to climate change was characterized by a reduction in the importance of the snow-melt pulse and an increase in the importance of fall and winter groundwater recharge. The less dynamic hydrologic regime is likely to result in drier soil conditions in rainfed wetlands and uplands, in contrast to less drying in groundwater-fed systems. Seepage lakes showed larger forecast stage declines related to climate change than did drainage lakes (lakes with outlet streams). Seepage lakes higher in the watershed (nearer to groundwater divides) had less groundwater inflow and thus had larger forecast declines in lake stage; however, ground-water inflow to seepage lakes in general tended to increase as a fraction of the lake budgets with lake-stage decline because inward hydraulic gradients increased. Drainage lakes were characterized by less simulated stage decline as reductions in outlet streamflow of set losses to other water flows. Net groundwater inflow tended to decrease in drainage lakes over the scenario period. Simulated stream temperatures increased appreciably with climate change. The estimated increase in annual average temperature ranged from approximately 1 to 2 degrees Celsius by 2100 in the stream characterized by a high groundwater inflow rate and 2 to 3 degrees Celsius in the stream with a lower rate. The climate drivers used for the climate-change scenarios had appreciable variation between the General Circulation Model and emission scenario selected; this uncertainty was reflected in hydrologic flow and temperature model results. Thus, as with all forecasts of this type, the results are best considered to approximate potential outcomes of climate change.

Wisconsin↗

Estimates of hydraulic properties from a one-dimensional numerical model of vertical aquifer-system deformation, Lorenzi site, Las Vegas, Nevada

Land subsidence related to aquifer-system compaction and ground-water withdrawals has been occurring in Las Vegas Valley, Nevada, since the 1930's, and by the late 1980's some areas in the valley had subsided more than 5 feet. Since the late 1980's, seasonal artificial-recharge programs have lessened the effects of summertime pumping on aquifer-system compaction, but the long-term trend of compaction continues in places. Since 1994, the U.S. Geological Survey has continuously monitored water-level changes in three piezometers and vertical aquifer-system deformation with a borehole extensometer at the Lorenzi site in Las Vegas, Nevada. A one-dimensional, numerical, ground-water flow model of the aquifer system below the Lorenzi site was developed for the period 1901-2000, to estimate aquitard vertical hydraulic conductivity, aquitard inelastic skeletal specific storage, and aquitard and aquifer elastic skeletal specific storage. Aquifer water-level data were used in the model as the aquifer-system stresses that controlled simulated vertical aquifer-system deformation. Nonlinear-regression methods were used to calibrate the model, utilizing estimated and measured aquifer-system deformation data to minimize a weighted least-squares objective function, and estimate optimal property values. Model results indicate that at the Lorenzi site, aquitard vertical hydraulic conductivity is 3 x 10-6 feet per day, aquitard inelastic skeletal specific storage is 4 x 10-5 per foot, aquitard elastic skeletal specific storage is 5 x 10-6 per foot, and aquifer elastic skeletal specific storage is 3 x 10-7 per foot. Regression statistics indicate that the model and data provided sufficient information to estimate the target properties, the model adequately simulated observed data, and the estimated property values are accurate and unique.

Water-Resources Investigations Report↗

Surface-water hydrology of the Western New York Nuclear Service Center Cattaraugus County, New York

Precipitation data were collected from October 1980 through September 1983 from three recording gages at the Western New York Nuclear Service Center, and surface water data were collected at three continuous-record gaging stations and one partial-record gage on streams that drain a 0.7 sq km part of the site. Seepage from springs was measured periodically during the study. The data were used to identify runoff characteristics at the waste burial ground and the reprocessing plant area, 400 meters to the north. Preliminary water budgets for April 1982 through March 1983 were calculated to aid in the development of groundwater flow models to the two areas. Nearly 80% of the measured runoff from the burial ground area was storm runoff; the remaining 20% was base flow. In contrast, only 30% of the runoff leaving the reprocessing plant area was storm runoff, and 70% was base flow. This difference is attributed to soil composition. The burial ground soil consists of clayey silty till that limits infiltration and causes most precipitation to flow to local channels as direct runoff. In contrast, the reprocessing plant area is overlain by alluvial sand and gravel that allows rapid infiltration of precipitation and subsequent steady discharge from the water table to nearby stream channels and seepage faces. Measured total annual runoff and estimated evapotranspiration from the reprocessing plant area exceeded the precipitation by 35%, which suggests that the groundwater basin is larger than the surface water basin. The additional outflow probably includes underflow from bedrock upgradient from the plant, water leakage from plant facilities, and groundwater flow from adjacent basins. (Author 's abstract)

New York↗

Simulation of freshwater-saltwater interfaces in the Brooklyn-Queens aquifer system, Long Island, New York

The seaward limit of the fresh ground-water system underlying Kings and Queens Counties on Long Island, N.Y., is at the freshwater-saltwater transition zone. This zone has been conceptualized in transient-state, three-dimensional models of the aquifer system as a sharp interface between freshwater and saltwater, and represented as a stationary, zero lateral-flow boundary. In this study, a pair of two-dimensional, four-layer ground-water flow models representing a generalized vertical section in Kings County and one in adjacent Queens County were developed to evaluate the validity of the boundary condition used in three-dimensional models of the aquifer system. The two-dimensional simulations used a model code that can simulate the movement of a sharp interface in response to transient stress. Sensitivity of interface movement to four factors was analyzed; these were (1) the method of simulating vertical leakage between freshwater and saltwater; (2) recharge at the normal rate, at 50-percent of the normal rate, and at zero for a prolonged (3-year) period; (3) high, medium, and low pumping rates; and (4) pumping from a hypothetical cluster of wells at two locations. Results indicate that the response of the interfaces to the magnitude and duration of pumping and the location of the hypothetical wells is probably sufficiently slow that the interfaces in three-dimensional models can reasonably be approximated as stationary, zero-lateral- flow boundaries.

New York↗