USGS Science⌕ Search

SEARCH · USGS Science

Results for “Ground Water Modeling Investigation Report”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

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

At least 919 records · Page 51Linked to original sources
USGS2005

Simulated ground-water flow for a pond-dominated aquifer system near Great Sandy Bottom Pond, Pembroke, Massachusetts

A ground-water flow simulation for a 66.4-square-mile area around Great Sandy Bottom (GSB) Pond (105 acres) near Pembroke, Massachusetts, was developed for use by local and State water managers to assess the yields for public water supply of local ponds and wells for average climatic and drought conditions and the effects of water withdrawals on nearby water levels and streamflows. Wetlands and ponds cover about 30 percent of the study area and the aquifer system is dominated by interactions between ground water and the ponds. The three largest surface-water bodies in the study area are Silver Lake (640 acres), Monponsett Pond (590 acres), and Oldham Pond (236 acres). The study area is drained by tributaries of the Taunton River to the southwest, the South and North Rivers to the northeast, and the Jones River to the southeast. In 2002, 10.8 million gallons per day of water was exported from ponds and 3.5 million gallons per day from wells was used locally for public supply. A transient ground-water-flow model with 69 monthly stress periods spanning the period from January 1998 through September 2003 was calibrated to stage at GSB Pond and nearby Silver Lake and streamflow and water levels collected from September 2002 through September 2003. The calibrated model was used to assess hydrologic responses to a variety of water-use and climatic conditions. Simulation of predevelopment (no pumping or export) average monthly (1949-2002) water-level conditions caused the GSB Pond level to increase by 6.3 feet from the results of a simulation using average 2002 pumping for all wells, withdrawals, and exports. Most of this decline can be attributed to pumping, withdrawals, and exports of water from sites away from GSB Pond. The effects of increasing the export rate from GSB Pond by 1.25 and 1.5 times the 2002 rate were a lowering of pond levels by a maximum of 1.6 and 2.8 feet, respectively. Simulated results for two different drought conditions, one mild drought similar to that of 1979-82 and a more severe drought similar to that of 1963-66, but with current (2002) pumping, were compared to results for average monthly recharge conditions (1949-2002). Simulated mild drought conditions showed a reduction of GSB Pond level of about 1.3 feet and a lower streamflow of about 1.7 percent in the nearby stream. Simulated severe drought conditions reduced the pond level at GSB Pond by almost 7 feet and lowered streamflow by about 37 percent. Varying cranberry-irrigation practices had little effect on simulated GSB Pond water levels, but may be important in other ponds. The model was most sensitive to changes in areal recharge. An increase and decrease of 22 percent in recharge produced changes in the GSB Pond water level of +1.4 feet and -2.4 feet, respectively. The accuracy of simulation results was best in the central portion of the study area in the immediate location of GSB Pond. The model was developed with the study-area boundary far enough away from the GSB Pond area that the boundary would have minimal effect on the water levels in GSB Pond, nearby ponds, and the underlying aquifer system. The model is best suited for use by local and State water managers to assess the effects of different withdrawal scenarios for wells and ponds near GSB Pond and for general delineation of areas contributing recharge to wells and ponds in the vicinity of GSB Pond. The model in its current form may not be well suited to detailed analyses of water budgets and flow patterns for parts of the study area farther from GSB Pond without further investigation, calibration, and data collection.

Scientific Investigations Report↗
USGS1985

Projected ground-water development, ground-water levels, and stream-aquifer leakage in the South Fork Solomon River Valley between Webster Reservoir and Waconda Lake, north-central Kansas, 1979-2020

A two-dimensional finite difference computer model was used to project changes in the potentiometric surface, saturated thickness, and stream aquifer leakage in an alluvial aquifer resulting from four instances of projected groundwater development. The alluvial aquifer occurs in the South Fork Solomon River valley between Webster Reservoir and Waconda Lake in north-central Kansas. In the first two projections, pumpage for irrigation was held constant at 1978 rates throughout the projection period (1979-2020). In the second two projections, the 1978 pumpage was progressively increased each yr through 2020. In the second and fourth projections, surface water diversions in the Osborne Irrigation Canal were decreased by 50 %. For the third and fourth projections, each grid-block in the modeled area was classified initially as one of six types according to whether it represented irrigable or nonirrigable land, to its saturated thickness, to its location inside or outside the canal-river area, and to its pumping rate. The projected base-flow rates (leakage from the aquifer to the river) were lower during the irrigation season (June, July, and August) than during the other months of the yr because of the decline in hydraulic head produced by groundwater pumpage. Stream depletion, calculated as a decrease below the average (1970-78) estimated winter base-flow rate of 16.5 cu ft/sec, varied inversely with base flow. For the first two projections, a constant annual cycle of well pumpage and recharge was used throughout the projection period. Aquifer leakage to the river was nearly constant by the mid-to-late 1990's, implying that flow conditions had attained a stabilized annual cycle. The third and fourth projections never attained an annual stabilized cycle because the irrigation pumpage rate was increased each year. By the early 1980's, the hydraulic head had fallen below river stage, reversing the hydraulic gradient at the stream-aquifer interface and resulting in net leakage from the river to the aquifer during the summer months. By the early 1990 's, the projected potentiometric surface of the aquifer was lower than the river stage even during the winter and spring months. (Author 's abstract)

Water-Resources Investigations Report↗
USGS2001

Simulated ground-water flow and water quality of the Mississippi River alluvium near Burlington, Iowa, 1999

The City of Burlington, Iowa, obtains some of its public water supply by withdrawing ground water from the Mississippi River alluvium, an alluvial aquifer adjacent to the Mississippi River. The U.S. Geological Survey, in cooperation with the City of Burlington, conducted a hydrologic study of the Mississippi River alluvium near Burlington in 1999 to improve understanding of the flow system, evaluate the effects of hypothetical pumping scenarios on the flow system, and evaluate selected water-quality constituents in parts of the alluvium. A steady-state, ground-water flow model was constructed for a 7-square-mile area of the alluvium using October 1999 hydrologic conditions to help conceptualize the flow system, identify sources of water to the alluvium, and assess potential effects from additional hypothetical ground-water withdrawals from the lower alluvium. The model was discretized into a 70-row by 68-column grid using cells measuring 200 feet by 200 feet. Three model layers were used to represent flow in the upper part of the alluvium, lower part of the alluvium, and bedrock. The primary sources of ground water to the alluvium were subsurface flow from areas of the alluvium adjacent to the modeled area, recharge from precipitation, subsurface flow from Flint River streamchannel deposits adjacent to the alluvium, and river leakage. The primary components of outflow from the flow system were river leakage, municipal ground-water withdrawals (pumpage), and leakage to drainage ditches. Three hypothetical pumping scenarios were used to assess the potential effects of increased ground-water withdrawals from the lower part of the alluvium: (1) pumping a second existing municipal well at a rate of 0.5 million gallons per day, (2) pumping a hypothetical well completed in an area between the city water-treatment facility and Flint River at a rate of 1.0 million gallons per day, and (3) pumping a hypothetical well completed in an area south of the Flint River at a rate of 1.0 million gallons per day. Maximum additional simulated drawdown in the upper alluvium ranged from less than 3 feet (for scenario 1) to about 9 feet (for scenario 3). Maximum additional simulated drawdown in the lower alluvium ranged from about 12 feet (for scenario 1) to about 34 feet (for scenario 3). Water budgets for each scenario indicated future additional withdrawals from the flow system near Burlington’s existing municipal wells would significantly increase the amount of river leakage into the flow system. Water samples collected from the alluvium indicated ground water can be classified as a calcium-magnesium-bicarbonate type. Reducing conditions likely occur in some localized areas of the alluvium, as suggested by relatively large concentrations of dissolved iron (4,390 micrograms per liter) and manganese (2, 430 micrograms per liter) in some ground-water samples. Nitrite plus nitrate was detected at concentrations greater than or equal to 8 milligrams per liter in three samples collected from observation wells completed in close proximity to cropland; the nitrite plus nitrate concentration in one groundwater sample exceeded the U.S. Environmental Protection Agency Maximum Contaminant Level for nitrate in drinking water (10 milligrams per liter as N). Triazine herbicides (atrazine, cyanazine, propazine, simazine, and selected degradation products) and chloroacetanilide herbicides (acetochlor, alachlor, and metolachlor) were detected in some water samples. A greater number of herbicide compounds were detected in surface-water samples than in ground-water samples. Herbicide concentrations typically were at least an order of magnitude greater in surfacewater samples than in ground-water samples. The Maximum Contaminant Level for alachlor (2 micrograms per liter) was exceeded in a sample from Dry Branch Creek at Tama Road and for atrazine (3 micrograms per liter) was exceeded in samples collected from Dry Branch Creek at Tama Road and the county drainage ditch at Tama Road.

Iowa↗
USGS1996

Hydrologic and chemical interaction of the Arkansas River and the Equus Beds aquifer between Hutchinson and Wichita, south-central Kansas

Large chloride concentrations in Arkansas River water have the potential to degrade water quality in the adjacent Equus beds aquifer between Hutchinson and Wichita, Kansas. The aquifer is an important source of water for municipal, industrial, agricultural, and domestic uses. A three-dimensional, finite-difference, ground-water flow-model program (MODFLOW) was used with data from past studies and data collected during 1988-91 to simulate aquifer and stream conditions during the late 1930's, during 1940-89, and during 1990-2019. Results of ground-water flow-model simulations indicated that declining water levels in the Equus beds aquifer since the 1940's have caused base flow in the Arkansas and Little Arkansas Rivers to decrease. In 1940, the Arkansas and Little Arkansas Rivers had simulated net base-flow gains within the model area of about 21 and about 67 ft 3 /s (cubic feet per second), respectively. By the end of 1989, the Arkansas River had a simulated net base-flow loss of about 52 ft 3 /s, and the Little Arkansas River had a net base-flow gain of about 27 ft 3 /s. Simulations for 1990-2019 showed that the water-level changes in a selected model cell located in the central part of the Wichita well field could range from -0.2 to -78 feet. Waterlevel changes in a selected model cell located near the Arkansas River could range from +1.3 to -1.2 feet. In model simulations where only pumpage varied, net base-flow loss from the Arkansas River to the aquifer ranged from about 59 ft 3 /s (no increase in pumpage since 1989) to 117 ft 3 /s (a 3-percent per year increase in pumpage since 1989) by 2019. Assuming a chloride concentration of 630 milligrams per liter, the median concentration in Arkansas River water collected during 1988-91, the quantity of chloride discharged from the Arkansas River to the aquifer was estimated to have increased from about 21 tons per day in 1940 to about 100 tons per day in 1989. By 2019, chloride discharge was indicated to range from about 110 tons per day (associated with no increase in pumpage since 1989) to 200 tons per day (associated with a 3-percent per year increase in pumpage since 1989). A particle-tracking program (MODPATH), which used the results from the flow model, was used to simulate the distribution in the aquifer of chloride from the river during the same time periods. Particle-tracking simulations show that, during 1940-89, the simulated distribution of particles representing chloride from the Arkansas River expanded from relatively narrow bands near the river to a wider distribution within the aquifer and the Wichita well field. Particle-tracking simulations indicate that chloride discharge from the Arkansas River may have reached the edge of the Wichita well field as early as 1963.

Kansas↗
USGS1998

Simulated response to pumping stress in the Sparta aquifer of southeastern Arkansas and north-central Louisiana, 1998-2027

The Sparta aquifer in southeastern Arkansas and north-central Louisiana is a major water resource for municipal, industrial, and agricultural uses. In recent years, the demand for water in some areas has resulted in withdrawals from the Sparta that significantly exceed recharge to the aquifer. Considerable drawdown has occurred in the potentiometric surface, and water users and managers alike have begun to question the ability of the aquifer to supply water for the long term. Large cones of depression are centered beneath the Grand Prairie area and the cities of Pine Bluff and El Dorado in Arkansas, and Monroe in Louisiana. Water levels in the aquifer have declined at rates greater than 1 foot per year for more than a decade in much of southern Arkansas and northern Louisiana and are now below the top of the formation in parts of Union and Columbia Counties, Arkansas, and in several areas of Louisiana. Problems related to over draft in the Sparta could result in increased drilling and pumping costs, loss of yield, salt-water intrusion, and decrease in water quality in areas of large drawdown. The effects of current ground-water withdrawals and potential future withdrawals on water availability are major concerns of water managers and users as well as the general public in the two States. The Sparta model-a regional scale, digital ground-water flow model-was first calibrated in the mid-1980's. The model was updated and reverified using 1995-97 data. Visual inspection of the observed (1996-97) and simulated potentiometric surfaces, statistical analysis of the error for the original calibration and current reverification, and comparison of observed versus simulated hydro graphs indicates that the model is simulating conditions in the aquifer within acceptable error, and the quality of current (1998) model results is similar to the original model calibration results. When stressed with current withdrawal volumes and distributions, the model is able to simulate currently observed heads effectively as heads were simulated in the original calibration period. Five pumping scenarios were simulated over a 30-year period based on (1) current pumping rates, (2) current rates of change in pumping, (3) decreased pumping in selected areas, (4) increased pumping in selected areas, and (5) redistribution and increase of pumping in selected areas. Model results show that although continued pumping at current rates will result in relatively minor declines in water levels (scenario 1 above), continued pumping at currently observed rates of change will result in drastic declines across large areas of focused withdrawals (scenario 2). Under the first scenario-in which current pumping rates are input to the model for the 30-year simulation period-water levels in the middle of the cones of depression centered on El Dorado and Monroe decrease less than 10 feet. In the second scenario-in which the current rate of change in pumpage is applied to the model-substantial declines occur in the proximity of most major pumpage centers. During the 1998-2027 model period, predicted water levels decline from 307 feet below sea level to 438 feet below sea level near El Dorado, from 58 feet below sea level to 277 feet below sea level near Pine Bluff, but only by about 25 feet-from 202 feet below sea level to 225 feet below sea level near Monroe. In the third scenario-in which minimum predicted water use figures supplied by selected facilities in Arkansas and decreased pumping estimates for Louisiana are applied to the model-simulated water levels are substantially higher at cones of depression around the major pumping centers of Monroe and El Dorado as compared to initial (1997) values. During the 1998-2027 model period, predicted water levels near Monroe increase from 202 feet below sea level to 133 feet below sea level; water levels near El Dorado increase from 307 feet below sea level to 123 feet below sea level. For the fourth scenario-in which maxi mum pr

Water-Resources Investigations Report↗
USGS1990

Regional evaluation of hydrologic factors and effects of pumping, St Peter-Jordan aquifer, Iowa

The St. Peter-Jordan aquifer includes the Cambrian Jordan Sandstone and the overlying Ordovician Prairie du Chien Group and St. Peter Sandstone. The aquifer is present throughout Iowa and is confined beneath other aquifers in much of the State. Information on the aquifer available from drillers and contractors, provided estimates of aquifer transmissivity values ranging from about 500 to about 3,000 square feet per day. The largest transmissivity values are for dolomite and dolomite-cemented sandstone, indicating that permeability in much of the aquifer is due to secondary fractures. The aquifer is vertically bounded by an upper leaky confining unit with a vertical hydraulic conductivity of 10 -10 feet per second. The aquifer was simulated using a two-layer finite-difference ground-water flow model. The upper layer simulated a source bed in aquifers composed of Silurian and Devonian rocks overlying the St. Peter-Jordan aquifer. The lower layer simulated flow in the St. Peter-Jordan aquifer. Lateral boundaries assigned in the model include constant heads in northeastern Iowa, where the aquifer is in contact with the Mississippi River or is unconfined, and no-flow boundaries in western and northwestern Iowa, where the rocks are insufficiently permeable to form an aquifer. The aquifer boundaries to the north, east, and south of Iowa were determined by geohydrologic conditions and the relation of the St. Peter-Jordan aquifer with the lateral extent of adjacent aquifers. An assumption that the largest part of recharge to the aquifer is from outcrop areas in northeastern Iowa and from Minnesota is not supported by the results of this study. Vertical leakage from overlying rocks accounted for most of the recharge to the aquifer in northwestern Iowa. Discharge is mostly through lateral boundaries and to rivers. Pumping has caused changes in the flow system that include regional declines in the potentiometric surface of the aquifer. Simulation indicates that pumping through 1980 increased net vertical leakage into the aquifer to about double the predevelopment rate. Discharge across lateral boundaries has been substantially reduced or reversed by pumping. Aquifer storage provided about one-third of the water required to supply pumping in the 1970's. Simulation of future conditions, assuming no increase in pumping rates, indicates that the rate of decline in water levels will decrease by the year 2020. As equilibrium with pumping is approached in 2020, 75 percent of the pumpage will be balanced by vertical leakage, eight percent by water released from aquifer storage, and 17 percent by increases in boundary recharge or decreases in boundary discharge. Future pumping at an increasing rate of about 10 percent per decade of the average pumping rate in 1975 will require about one and one-half times the vertical leakage of the 1971-1980 period and about fivetimes the net inflow from lateral boundaries; however, the rate of water released from aquifer storage will be about half the 1970's rate. Under these conditions, the head in the aquifer will continue to decline at an almost constant rate until 2020.

Illinois, Indiana, Iowa, Minnesota, Missouri, Wiis↗
USGS2009

Modeling Flood Plain Hydrology and Forest Productivity of Congaree Swamp, South Carolina

An ecological field and modeling study was conducted to examine the flood relations of backswamp forests and park trails of the flood plain portion of Congaree National Park, S.C. Continuous water level gages were distributed across the length and width of the flood plain portion - referred to as 'Congaree Swamp' - to facilitate understanding of the lag and peak flood coupling with stage of the Congaree River. A severe and prolonged drought at study start in 2001 extended into late 2002 before backswamp zones circulated floodwaters. Water levels were monitored at 10 gaging stations over a 4-year period from 2002 to 2006. Historical water level stage and discharge data from the Congaree River were digitized from published sources and U.S. Geological Survey (USGS) archives to obtain long-term daily averages for an upstream gage at Columbia, S.C., dating back to 1892. Elevation of ground surface was surveyed for all park trails, water level gages, and additional circuits of roads and boundaries. Rectified elevation data were interpolated into a digital elevation model of the park trail system. Regression models were applied to establish time lags and stage relations between gages at Columbia, S.C., and gages in the upper, middle, and lower reaches of the river and backswamp within the park. Flood relations among backswamp gages exhibited different retention and recession behavior between flood plain reaches with greater hydroperiod in the lower reach than those in the upper and middle reaches of the Congaree Swamp. A flood plain inundation model was developed from gage relations to predict critical river stages and potential inundation of hiking trails on a real-time basis and to forecast the 24-hour flood In addition, tree-ring analysis was used to evaluate the effects of flood events and flooding history on forest resources at Congaree National Park. Tree cores were collected from populations of loblolly pine (Pinus taeda), baldcypress (Taxodium distichum), water tupelo (Nyssa aquatica), green ash (Fraxinus pennslyvanica), laurel oak (Quercus laurifolia), swamp chestnut oak (Quercus michauxii), and sycamore (Plantanus occidentalis) within Congaree Swamp in highand low-elevation sites characteristic of shorter and longer flood duration and related to upriver flood controls and dam operation. Ring counts and dating indicated that all loblolly pine trees and nearly all baldcypress collections in this study are postsettlement recruits and old-growth cohorts, dating from 100 to 300 years in age. Most hardwood species and trees cored for age analysis were less than 100 years old, demonstrating robust growth and high site quality. Growth chronologies of loblolly pine and baldcypress exhibited positive and negative inflections over the last century that corresponded with climate history and residual effects of Hurricane Hugo in 1989. Stemwood production on average was less for trees and species on sites with longer flood retention and hydroperiod affected more by groundwater seepage and site elevation than river floods. Water level data provided evidence that stream regulation and operations of the Saluda Dam (post-1934) have actually increased the average daily water stage in the Congaree River. There was no difference in tree growth response by species or hydrogeomorphic setting to predam and postdam flood conditions and river stage. Climate-growth analysis showed that long-term growth variation is controlled more by spring/ summer temperatures in loblolly pine and by spring/summer precipitation in baldcypress than flooding history.

South Carolina↗
USGS2001

Ground-water recharge and flowpaths near the edge of the Decorah-Platteville-Glenwood confining unit, Rochester, Minnesota

The primary source of ground water for the city of Rochester, Olmsted County, southeastern Minnesota is the St. Peter-Prairie du Chien-Jordan aquifer. Based on results of a previous U.S. Geological Survey investigation in the Rochester area, relatively high rates of areal recharge to the St. Peter-Prairie du Chien-Jordan aquifer occur along the edge of the overlying Decorah-Platteville-Glenwood confining unit. The primary source of water to the zone of increased recharge along the edge of the confining unit is the upper carbonate aquifer. Ground-water recharge rates to the St. Peter-Prairie du Chien-Jordan aquifer during all of 1998 ranged from 1.9 to 25.5 in./yr (inches per year). Recharge rates were greatest near the edge of the Decorah-Platteville-Glenwood confining unit and least where the confining unit is thick and is overlain by the upper carbonate aquifer (mean of 2.0 in./yr). Recharge rates downslope from the edge of the confining unit were greatest to the St. Peter on the slope entering the main South Fork Zumbro River Valley. Results of groundwater age dating using chlorofluorocarbons (CFCs) indicated recharge dates ranging from (1) the mid-1950’s to the early 1990’s for the St. Peter, (2) the late 1960’s to approximately 1990 for the Prairie du Chien, and (3) the early to mid-1950’s for the Jordan. Cross-sectional model simulations indicated that most of the areal recharge entering the aquifer system through the upper carbonate aquifer discharges from springs and seeps. Of the 2.28 ft3/s (cubic feet per second) of areal recharge that enters the upper carbonate aquifer, 2.23 ft3/s is discharged from the aquifer by springs and seeps. Results indicate that areal recharge to the upper Prairie du Chien moves primarily westward and discharges to Bear Creek. Areal recharge rates derived from hydrograph analysis, CFC age-dating, and cross-sectional model analysis were much greater to the St. Peter downslope from the edge of the Decorah-Platteville-Glenwood confining unit (25.5, 35.3, and 23.75 in./yr, respectively) than occurs to the Prairie du Chien in any hydrogeologic setting. The model-simulated discharge from springs and seeps in the lower part of the upper carbonate aquifer represents a potential source of water of 33 in./yr to the St. Peter unit, similar to the estimated areal recharge rates derived from hydrograph analysis and CFC age-dating. The water withdrawn by pumped wells or discharged to Bear Creek is derived predominantly from areal recharge near the edge of the Decorah-Platteville-Glenwood confining unit (0.47 ft 3 /s), rather than from water that has leaked downward through the Decorah unit (0.03 ft 3 /s). Model simulated discharge through springs and seeps in the lower part of the upper carbonate aquifer (0.21 ft 3 /s) represents a potential source of water to the St. Peter-Prairie du Chien-Jordan aquifer.

Minnesota↗
USGS1999

Health of native riparian vegetation and its relation to hydrologic conditions along the Mojave River, southern California

The health of native riparian vegetation and its relation to hydrologic conditions were studied along the Mojave River mainly during the growing seasons of 1997 and 1998. The study concentrated on cottonwood-willow woodlands (predominantly Populus fremontii and Salix gooddingii ) and mesquite bosques (predominantly Prosopis glandulosa ). Tree-growth characteristics were measured at 16 cottonwood-willow woodland sites and at 3 mesquite bosque sites. Density of live and dead trees, tree diameter and height, canopy density, live-crown volume, leaf-water potential, leaf-area index, mortality, and reproduction were measured or noted at each site. The sites included healthy and reproducing woodlands and bosques, stressed woodlands and bosques with no reproduction, and woodlands and bosques with high mortality. Tree roots were studied at seven sites to determine the vertical distribution of the root system and their relation to the water table at healthy, stressed, and high-mortality cottonwood-willow woodlands. In the six trenches that were dug for this study in May 1997, no cottonwood roots were observed that reached the water table. The root systems of healthy trees typically ended 1 to 2 feet above the water table. At sites with high mortality, the main root mass was commonly 7 to 8 feet above the water table. Water-table depth was monitored at each of the study sites. In addition, volumetric soil moisture and soil-water potential were monitored at varying depths at three cottonwood-willow woodland study sites and at two mesquite bosque sites. Ground, soil, river, lake, and plant (xylem sap) water were analyzed for concentrations of stable hydrogen and oxygen isotopes to determine the source of water used by the trees. On the basis of the root-distribution, soil- and leaf-water potential, and isotope data, it was concluded that cottonwood, willow, and mesquite trees mainly rely on ground water for their perennial sustained supply of water. The trees mainly utilize ground water that has moved upward from the water table into the capillary fringe and into unsaturated soil nearer to land surface. Most precipitation (average is 4 to 6 inches per year) is lost by evaporation and by transpiration of shallow-rooted xeric plants, and very little reaches the root zone of trees along the Mojave River. Water-table depth had no strong correlation to many individual tree-growth characteristics, such as density, diameter, height, and live-crown volume. However, leaf-area index (corrected for stem area) of both healthy and stressed cottonwood-willow woodlands had a highly significant statistical relation to water-table depth, and a curvilinear regression model was defined. As in cottonwood-willow woodlands, leaf-area index of mesquite bosques also decreased with increased water-table depth. However, because of the small number of sites, no significant statistical relation could be defined for mesquite bosques. Because it can be accurately measured repeatedly at the same locations, leaf-area index (corrected for stem area) is recommended as the primary growth characteristic that should be monitored. Future vegetation changes along the Mojave River can be quantified using the sites established for this study. Mortality was as high as 39 percent in healthy cottonwood-willow woodlands, but mortality of 50 to 100 percent was common where water-table depth was greater than about 7 feet or in areas where permanent water-table declines greater than about 5 feet had occurred. At a healthy mesquite bosque where the water-table depth ranged from about 8 to 11 feet, mortality was about 20 percent. Where the water table had been lowered an additional 10 to 25 feet by pumping, mortality of the mesquite was extremely high (80 to 99 percent). On the basis of observations of plant reproduction, it was concluded that established cottonwood-willow woodlands probably will reproduce, mainly by root sprouting of mature trees, if the water-table depth during spring and early summer is less than about 5 feet and if there is sufficient sunlight. Likewise, reproduction in established mesquite bosques will probably occur by root sprouting if the water-table depth during spring and early summer is less than about 8 feet. Germination of cottonwood and willow seeds is common during spring in dry channels following winter stormflows in the Mojave River. However, most of the seedlings die before their roots reach a perennial water supply near the water table, and established stands of trees resulting from seed germination are rare.

California↗
USGS1996

Hydrogeology of the area near the J4 test cell, Arnold Air Force Base, Tennessee

The U.S. Air Force operates a major aerospace systems testing facility at Arnold Engineering Development Center (AEDC) in Coffee County, Tennessee. Dewatering operations at one of the test facilities, the J4 test cell, has affected the local ground-water hydrology. The J4 test cell is approximately 100 feet in diameter, extends approximately 250 feet below land surface, and penetrates several aquifers. Ground water is pumped continuously from around the test cell to keep the cell structurally intact. Because of the test cell's depth, dewatering has depressed water levels in the aquifers surrounding the site. The depressions that have developed exhibit anisotropy that is controlled by zones of high permeability in the aquifers. Additionally, contaminants - predominately volatile organic compounds - are present in the ground-water discharge from the test cell and in ground water at several other Installation Restoration Program (IRP) sites within the AEDC facility. The dewatering activities at J4 are drawing these contaminants from the nearby sites. The effects of dewatering at the J4 test cell were investigated by studying the lithologic and hydraulic characteristics of the aquifers, investigating the anisotropy and zones of secondary permeability using geophysical techniques, mapping the potentiometric surfaces of the underlying aquifers, and developing a conceptual model of the ground-water-flow system local to the test cell. Contour maps of the potentiometric surfaces in the shallow, Manchester, and Fort Payne aquifers (collectively, part of the Highland Rim aquifer system) show anisotropic water-level depressions centered on the J4 test cell. This anisotropy is the result of features of high permeability such as chert-gravel zones in the regolith and fractures, joints, and bedding planes in the bedrock. The presence of these features of high permeability in the Manchester aquifer results in complex flow patterns in the Highland Rim aquifers near the J4 test cell. The occurrence, distribution, and orientation of these features has a great effect on ground-water flow to the J4 test cell. The depression caused by dewatering extends out horizontally through the aquifers along the most permeable pathways. Since the aquifers above the Chattanooga Shale are not separated by distinct confining units, areas in adjacent aquifers above and below these zones of high permeability in the Manchester aquifer are also dewatered. Conditions in all Highland Rim aquifers approximate steady-state equilibrium because ground-water withdrawal at the test cell has been continuous since the late 1960's. The average ground-water discharge from the dewatering system at the J4 test cell was 105 gallons per minute, for 1992-95. The ground-water capture areas in each aquifer extend into all or parts of landfill #2 and leaching pit #2 (IRP site 1), the main testing area (IRP site 7), and the old fire training area (IRP site 10). IRP sites 8 and 12 are outside the ground-water capture areas. Of the 35 sampled wells in the J4 area, 10 produced water samples containing chlorinated organic compounds such as 1,2-dichloroethane (1,2-DCA), 1,1-dichloroethylene (1,1-DCE), and trichloroethylene (TCE) in concentrations which exceeded the Tennessee Department of Environment and Conservation Maximum Contaminant Levels (MCL's) for public water-supply systems. The highest concentrations were detected in samples from well AEDC-274 with 45 micrograms per liter (mg/L) 1,2-DCA, 320 mg/L 1,1-DCE, and 1,200 mg/L TCE. These compounds are synthetic and do not occur naturally in the environment. A sample of the ground-water discharge from the J4 test cell also contained concentrations of these compounds that exceed MCL's. Chlorinated organic compounds, including 1,2-DCA; 1,1-DCE; and TCE also have been detected at IRP sites 1, 7, 8, nd 10. The six dewatering wells surrounding the J4 test cell penetrate the Chattanooga Shale and are open to the Highland Rim aquifer system, there

Water-Resources Investigations Report↗
USGS2006

Sources of water to wells in updip areas of the Wenonah-Mount Laurel aquifer, Gloucester and Camden Counties, New Jersey

Since 1996, when the New Jersey Department of Environmental Protection (NJDEP) restricted ground-water withdrawals from the Potomac-Raritan-Magothy aquifer system in the southern New Jersey Coastal Plain as a result of excessive drawdown, Coastal Plain communities have been interested in developing alternate sources of water supply for their residents. The use of ground water from areas near the updip parts of the overlying confined aquifers where withdrawals are not restricted is being considered to meet the demand for drinking water. Concerns have arisen, however, regarding the potential effects of increased withdrawals from these areas on ground-water flow to streams and wetlands as well as to the deeper, confined parts of the aquifers. Therefore, the U.S. Geological Survey, in cooperation with the NJDEP, conducted a study to investigate the sources of water to currently inactive wells in the updip part of the Wenonah-Mount Laurel aquifer in Gloucester and Camden Counties, New Jersey. Of particular interest is whether the primary source of the increased withdrawals is likely to be the aquifer outcrop or the downdip, confined part of the aquifer. The outcrop of the Wenonah-Mount Laurel aquifer covers nearly 8 mi2 (square miles), or about 46 percent of Deptford Township's 17.56-mi2 area. The Deptford Township Municipal Utilities Authority owns six currently (2005) inactive wells in the Wenonah-Mount Laurel aquifer at the southeastern boundary of Deptford Township, 1.25 mi (miles) from the outcrop. For the purposes of this study, an existing ground-water-flow model of the New Jersey Coastal Plain aquifers was used to simulate ground-water-flow conditions in Gloucester and Camden Counties in 1998. Two alternative withdrawal scenarios were superimposed on the results of the 1998 simulation. In the first (the 'full-allocation' scenario), full-allocation withdrawal rates established by the NJDEP were applied to 45 existing wells in the Deptford Township area. In the second (the 'additional-withdrawal' scenario), the full-allocation scenario was modified by adding an additional withdrawal of 1.62 million gallons per day from the six inactive Deptford Township withdrawal wells. Simulated drawdown for the full-allocation scenario is zero to near zero in Deptford Township. Changes are greatest downdip from Deptford Township, where a broad area of 5- to 10-ft (feet) drawdowns is simulated; maximum drawdown at the center of the cone of depression is 20 ft. Water levels declined as much as 10 ft around individual wells whose current withdrawals are only a small percentage of their allotted allocation. Simulated drawdown for the additional-withdrawal scenario exceeds 40 ft and is centered around the six inactive Deptford Township withdrawal wells. The area in which the simulated drawdown is 5 ft extends approximately 3.75 mi downdip from the wells and 2 mi updip, into the outcrop. Water budgets based on the simulation results for the full-allocation and additional-withdrawal scenarios were calculated and compared, with particular focus on a 75-mi2 area in and around Deptford Township that includes the outcrop of the Wenonah-Mount Laurel aquifer and part of the area downdip from the outcrop (budget zone 2). The comparison of the two water budgets for zone 2 shows that 46 percent of the withdrawals from the six inactive Deptford Township wells would result from reduced stream base flow in the outcrop of the Wenonah-Mount Laurel aquifer and 35 percent would result from increased downward flow from the overlying Vincentown aquifer. Four percent would result from increased flow from the downdip areas of the Wenonah-Mount Laurel aquifer, 5 percent would result from decreased flow to the downdip areas of the Wenonah-Mount Laurel aquifer, and 5 percent would result from decreased flow to the underlying Englishtown aquifer system. The remaining 4 percent was attributed to decreased upward flow to the overlying Vincentown aquifer.

Scientific Investigations Report↗
USGS2002

Estimates of ground-water recharge from precipitation to glacial-deposit and bedrock aquifers on Lopez, San Juan, Orcas, and Shaw islands, San Juan County, Washington

An important source of fresh water on Lopez, San Juan, Orcas, and Shaw Islands in San Juan County off the northwestern coast of Washington is glacial-deposit and bedrock aquifers. Two methods were used to estimate recharge from precipitation to the water tables on the islands. A daily near-surface water-balance method, the Deep Percolation Model (DPM), was used to simulate water budgets for the period October 1, 1996, through September 30, 1998 (water years 1997-98) for six small drainage basins?three on Lopez Island and one each on San Juan, Orcas, and Shaw Islands. The calibrated soil and subsoil parameters from the DPM for each small basin were then used in island-wide applications of the DPM where the direct runoff component (which is not available on an island-wide basis) was simulated, rather than input, and calibration was not required. A spatial distribution of annual recharge was simulated for each island, with island averages of: Lopez Island, 2.49 inches per year; San Juan Island, 1.99 inches per year; Orcas Island, 1.46 inches per year; and Shaw Island, 1.44 inches per year.A chloride mass-balance method that requires measurements of atmospheric chloride deposition, precipitation, streamflow, and chloride concentrations in ground water was used to estimate recharge to the glacial-deposit aquifers of Lopez Island. Only average recharge could be estimated using this method rather than area-specific recharge. Average recharge for Lopez Island estimated by this method was only 0.63 inch per year. The range of chloride concentrations in ground-water samples from selected wells indicates that the average recharge in areas of glacial deposits is between 0.29 and 1.95 inches per year. Recharge simulated using the DPM for two drainage basins on Lopez Island overlain by glacial deposits are 2.76 and 2.64 inches per year. Sources of chloride in ground water other than from the atmosphere would cause the recharge estimated by the chloride mass-balance method to be less than the actual recharge, therefore these estimates may represent lower limits which are, at least, consistent with the higher simulated recharge from the DPM. The average island-wide recharge is most closely related to the amount of area overlain by glacial deposits. Thus, even though Lopez Island receives the least precipitation, it has the most recharge per square mile because it proportionally has the largest area overlain by glacial deposits. Recharge simulated by the DPM for areas of shallow to outcropping bedrock generally were less than 1.5 inches per year, but recharge simulated in areas of glacial deposits ranged from less than 0.5 to 3 inches per year, with recharge as high as 9 inches per year in some small areas.

Water-Resources Investigations Report↗
USGS2005

Simulation of ground-water flow in the basin-fill aquifer of the Tularosa Basin, south-central New Mexico, predevelopment through 2040

The hydrology of the basin-fill aquifer in the Tularosa Basin was evaluated through construction and calibration of steady-state and transient three-dimensional ground-water-flow simulations. Simulations were made using the U.S. Geological Survey finite-difference modular ground-water-flow computer software MODFLOW-96. The transient simulation covered 1948-2040. Both steady-state and transient simulations were calibrated by matching simulation output to available ground-water-level measurements. The root-mean-square error of the steady-state calibration in the well-calibrated area of the ground-water-flow simulation was 6.3 meters, and root-mean-square errors of individual transient-calibration points ranged from 0.8 to 17.0 meters. The areal distribution of water-level measurements used in the steady-state and transient calibrations restricts the well-calibrated area of the model to the eastern side of the Tularosa Basin. Water levels in the La Luz Creek subbasin area were underestimated by both the steady-state and transient models, suggesting that the hydrology of this area is not well represented in the model. About 143,000 cubic meters per day of recharge is estimated to enter the basin-fill aquifer from subbasins that rim the Tularosa Basin. The estimated recharge is about 4-5 percent of total precipitation in most subbasins. Approximately 88 percent of total recharge left the basin-fill aquifer as evapotranspiration under predevelopment conditions. Water levels were simulated for 1948, 1995, and 2040 under scenarios of zero and maximum return flows. Estimated return flows from municipalities were calculated on the basis of data in the Tularosa Basin Regional Water Plan for 2000-2040. Agricultural return flows were estimated primarily on the basis of ground-water-withdrawal, ground-water-depletion, surface-water-withdrawal, and surface-water-depletion data for the Tularosa Basin. The ground-water-flow simulation was sensitive to the return-flow scenario in the agricultural area near Tularosa and decreasingly sensitive to the south. Declines in simulated water levels near Tularosa between 1948 and 1995 were as large as 30 meters under the zero return-flow scenario and 15 meters under the maximum return-flow scenario. Declines in simulated water levels between 1995 and 2040 were as large as 25 meters under the zero return-flow scenario and 15 meters under the maximum return-flow scenario. Comparison of water levels measured near Tularosa in 1991 and water levels simulated under the maximum return-flow scenario for 1991 suggests that declines in simulated water levels near Tularosa may be overestimated under the zero return-flow scenario. Declines in simulated water levels near the City of Alamogordo well field between 1948 and 1995 were as large as 15 meters under the zero return-flow scenario and 10 meters under the maximum return-flow scenario. Simulated declines in water levels between 1995 and 2040 were nearly 15 meters under both return-flow scenarios assuming that all projected increases in withdrawal came from existing City of Alamogordo public-supply wells and all withdrawal from the wells came from the basin-fill aquifer. Declines in simulated water levels near the Holloman Air Force Base well fields between 1948 and 1995 and between 1995 and 2040 were less than 5 meters under both the zero and maximum return-flow scenarios. In 1995 under the zero return-flow scenario, an estimated 56,000 cubic meters of water per day was removed from aquifer storage. Of the approximately 199,000 cubic meters of water per day that left the aquifer under 1995 conditions, 40 percent left the basin-fill aquifer as ground-water withdrawal, 51 percent as evapotranspiration, 7 percent by interbasin ground-water flow into the Hueco Bolson, and 2 percent by flow into creeks and springs. Generalized directions of ground-water flow were simulated for 1948, 1995, and 2040 for much of the eastern part of the Tularosa Basin. Localized

Scientific Investigations Report↗
USGS2001

Simulation of a long-term aquifer test conducted near the Rio Grande, Albuquerque, New Mexico

A long-term aquifer test was conducted near the Rio Grande in Albuquerque during January and February 1995 using 22 wells and piezometers at nine sites, with the City of Albuquerque Griegos 1 production well as the pumped well. Griegos 1 discharge averaged about 2,330 gallons per minute for 54.4 days. A three-dimensional finite-difference ground-water-flow model was used to estimate aquifer properties in the vicinity of the Griegos well field and the amount of infiltration induced into the aquifer system from the Rio Grande and riverside drains as a result of pumping during the test. The model was initially calibrated by trial-and-error adjustments of the aquifer properties. The model was recalibrated using a nonlinear least-squares regression technique. The aquifer system in the area includes the middle Tertiary to Quaternary Santa Fe Group and post-Santa Fe Group valley- and basin-fill deposits of the Albuquerque Basin. The Rio Grande and adjacent riverside drains are in hydraulic connection with the aquifer system. The hydraulic-conductivity values of the upper part of the Santa Fe Group resulting from the model calibrated by trial and error varied by zone in the model and ranged from 12 to 33 feet per day. The hydraulic conductivity of the inner-valley alluvium was 45 feet per day. The vertical to horizontal anisotropy ratio was 1:140. Specific storage was 4 x 10-6 per foot of aquifer thickness, and specific yield was 0.15 (dimensionless). The sum of squared errors between the observed and simulated drawdowns was 130 feet squared. Not all aquifer properties could be estimated using nonlinear regression because of model insensitivity to some aquifer properties at observation locations. Hydraulic conductivity of the inner-valley alluvium, middle part of the Santa Fe Group, and riverbed and riverside-drain bed and specific yield had low sensitivity values and therefore could not be estimated. Of the properties estimated, hydraulic conductivity of the upper part of the Santa Fe Group was estimated to be 12 feet per day, the vertical to horizontal anisotropy ratio was estimated to be 1:82, and specific storage was estimated to be 1.2 x 10-6 per foot of aquifer thickness. The overall sum of squared errors between the observed and simulated drawdowns was 87 feet squared, a significant improvement over the model calibrated by trial and error. At the end of aquifer-test pumping, induced infiltration from the Rio Grande and riverside drains was simulated to be 13 percent of the total amount of water pumped. The remainder was water removed from aquifer storage. After pumping stopped, induced infiltration continued to replenish aquifer storage. Simulations estimated that 5 years after pumping began (about 4.85 years after pumping stopped), 58 to 72 percent of the total amount of water pumped was replenished by induced infiltration from the Rio Grande surface-water system.

Water-Resources Investigations Report↗
USGS2018

Hydrogeologic controls and geochemical indicators of groundwater movement in the Niles Cone and southern East Bay Plain groundwater subbasins, Alameda County, California

Beginning in the 1970s, Alameda County Water District began infiltrating imported water through ponds in repurposed gravel quarries at the Quarry Lakes Regional Park, in the Niles Cone groundwater subbasin, to recharge groundwater and to minimize intrusion of saline, San Francisco Bay water into freshwater aquifers. Hydraulic connection between distinct aquifers underlying Quarry Lakes allows water to recharge the upper aquifer system to depths of 400 feet below land surface, and the Deep aquifer to depths of more than 650 feet. Previous studies of the Niles Cone and southern East Bay Plain groundwater subbasins suggested that these two subbasins may be hydraulically connected. Characterization of storage capacities and hydraulic properties of the complex aquifers and the structural and stratigraphic controls on groundwater movement aids in optimal storage and recovery of recharged water and provides information on the ability of aquifers shared by different water management agencies to fulfill competing storage and extraction demands. The movement of recharge water through the Niles Cone groundwater subbasin from Quarry Lakes and the possible hydraulic connection between the Niles Cone and the southern East Bay Plain groundwater subbasins were investigated using interferometric synthetic aperture radar (InSAR), water-chemistry, and isotopic data, including tritium/helium-3, helium-4, and carbon-14 age-dating techniques. InSAR data collected during refilling of the Quarry Lakes recharge ponds show corresponding ground-surface displacement. Maximum uplift was about 0.8 inches, reasonable for elastic expansion of sedimentary materials experiencing an increase in hydraulic head that resulted from pond refilling. Sodium concentrations increase while calcium and magnesium concentrations in groundwater decrease along groundwater flowpaths from the Niles Cone groundwater subbasin through the Deep aquifer to the northwest toward the southern East Bay Plain groundwater subbasin. Residual effects of pre-1970s intrusion of saline water from San Francisco Bay, including high chloride concentrations in groundwater, are evident in parts of the Niles Cone subbasin. Noble gas recharge temperatures indicate two primary recharge sources (Quarry Lakes and Alameda Creek) in the Niles Cone groundwater subbasin. Although recharge at Quarry Lakes affects hydraulic heads as far as the transition zone between the Niles Cone and East Bay Plain groundwater subbasins (about 5 miles), the effect of recharged water on water quality is only apparent in wells near (less than 2 miles) recharge sources. Groundwater chemistry from upper aquifer system wells near Quarry Lakes showed an evaporated signal (less negative oxygen and hydrogen isotopic values) relative to surrounding groundwater and a tritium concentration (2 tritium units) consistent with recently recharged water from a surface-water impoundment. Uncorrected carbon-14 activities measured in water sampled from wells in the Niles Cone groundwater subbasin range from 16 to 100 percent modern carbon (pmC). The geochemical reaction modeling software NETPATH was used to interpret carbon-14 ages along a flowpath from Quarry Lakes toward the East Bay Plain groundwater subbasin. Model results indicate that changes in groundwater chemistry are controlled by cation exchange on clay minerals and weathering of primary silicate minerals. Old groundwater (lower carbon-14 activities) is characterized by high dissolved silica and pH. Interpreted carbon-14 ages ranged from 830 to more than 7,000 years before present and are less than helium-4 ages that range from 2,000 to greater than 11,000 years before present. The average horizontal groundwater velocity along the studied flowpath, as calculated using interpreted carbon-14 ages, through the Deep aquifer of the Niles Cone groundwater subbasin is between 3 and 12 feet per year. The groundwater velocity decreases near the boundary of the transition zone to the southern East Bay Plain groundwater subbasin to about 0.5 feet per year. These changes may result from water recharged from different sources converging in flowpaths north of the transition zone, or a boundary to flow between the Niles Cone and southern East Bay Plain groundwater subbasins, likely owing to changes in lithology caused by depositional patterns.

California↗
USGS1998

Areas contributing ground water to the Peconic Estuary, and ground-water budgets for the north and south forks and Shelter Island, eastern Suffolk County, New York

The Peconic Estuary, at the eastern end of Long Island, has been plagued by a recurrent algal bloom, locally referred to as ?Brown Tide,? that has caused the severe decline of local marine resources. Although the factors that trigger Brown Tide blooms remain uncertain, groundwater discharge has previously been shown to affect surface-water quality in the western part of the estuary. A U.S. Geological Survey groundwater- flow model of the main body of Long Island indicates that a total of about 7.5 x 106 ft3/d (cubic feet per day) of freshwater discharges to the western part of the estuary, but the model does not include the ground-water flow systems on the North and South Forks and Shelter Island, which contribute significant amounts of freshwater to the central and eastern parts of the estuary. The need for information on freshwater discharge to the entire estuary prompted the U.S. Geological Survey to evaluate ground-water discharge from the North and South Forks and Shelter Island. Source areas that contribute ground water to the Peconic Estuary were delineated, and groundwater budgets for these areas were developed, to evaluate the distribution and magnitude of ground-water discharge to the central and eastern parts of the estuary. Contributing-area boundaries that were delineated coincide with the hydraulic boundaries of the fresh ground-water-flow systems of the North and South Forks and Shelter Island; these boundaries are of two types? external (saltwater bodies) and internal (groundwater divides). Hydrologic components that were evaluated include recharge from precipitation, public-supply withdrawal and return flow, and agricultural withdrawal. Values for each of these components were calculated or estimated for the individual freshwater flow subsystems that form each ground-water-budget area, then summed to obtain the total discharge of fresh ground water to tidewater. Ground-water discharge to the Peconic Estuary is about 3.8 x 106 ft3/d from the North Fork, 11 x 106 ft3/d from the South Fork, and 1.7 x 106 ft3/d from Shelter Island. The total contribution to the estuary from these areas is about 16 x 106 ft3/d?roughly twice the total contribution from the main body of Long Island. In contrast to the freshwater contribution from the main body of Long Island, which is concentrated near the head of the estuary, the contributions from the North and South Forks and Shelter Island are distributed along the east-west length of the estuary. Changes in water-table altitude and the resulting changes in total discharge to the Peconic Estuary were estimated from the relative changes in annual mean water level at observation wells. The 1985-95 interval included 7 years (1985-88, 1991- 92, 1995) of generally below-average water-table altitudes that presumably caused similar decreases in ground-water discharge to the estuary; intense Brown Tide blooms coincided with six of these years (1985-88, 1991, 1995), and localized blooms coincided with the remaining year (1992). Watertable altitudes in the remaining 4 years of the 1985-95 interval (1989-90, 1993-94) were nearly average or above average, and presumably produced comparably near-average or increased amounts of ground-water discharge to the estuary; none of these years saw any widespread Brown Tide blooms. Fluctuations in the amounts of ground-water discharge to the estuary appear to affect the occurrence of Brown Tide blooms, although the factors that trigger the blooms have not been determined.

Water-Resources Investigations Report↗
USGS2002

Hydrogeology and simulation of ground-water flow and land-surface subsidence in the Chicot and Evangeline aquifers, Houston area, Texas

In November 1997, the U.S. Geological Survey, in cooperation with the City of Houston Utilities Planning Section and the City of Houston Department of Public Works & Engineering, began an investigation of the Chicot and Evangeline aquifers in the greater Houston area in Texas to better understand the hydrology, flow, and associated land-surface subsidence. The principal part of the investigation was a numerical finite-difference model (MODFLOW) developed to simulate ground-water flow and land-surface subsidence in an 18,100-square-mile area encompassing greater Houston. The focus of the study was Harris and Galveston Counties, but other counties were included to achieve the appropriate boundary conditions. The model was vertically discretized into three 103-row by 109-column layers resulting in a total of 33,681 grid cells. Layer 1 represents the water table using a specified head, layer 2 represents the Chicot aquifer, and layer 3 represents the Evangeline aquifer. Simulations were made under transient conditions for 31 ground-water-withdrawal (stress) periods spanning 1891–1996. The years 1977 and 1996 were chosen as potentiometric-surface calibration periods for the model. Simulated and measured potentiometric surfaces of the Chicot and Evangeline aquifers for 1977 match closely. Waterlevel measurements indicate that by 1977, large ground-water withdrawals in east-central and southeastern areas of Harris County had caused the potentiometric surfaces to decline as much as 250 feet below sea level in the Chicot aquifer and as much as 350 feet below sea level in the Evangeline aquifer. Simulated and measured potentiometric surfaces of the Chicot and Evangeline aquifers for 1996 also match closely. The large potentiometric-surface decline in 1977 in the southeastern Houston area showed significant recovery by 1996. The 1996 centers of potentiometric-surface decline are located much farther northwest. Potentiometric-surface declines of more than 200 feet below sea level in the Chicot aquifer and more than 350 feet below sea level in the Evangeline aquifer were measured in observation wells and simulated in the flow model. Simulation of land-surface subsidence and water released from storage in the clay layers was accomplished using the Interbed-Storage Package of the MODFLOW model. Land-surface subsidence was calibrated by comparing simulated long-term (1891–1995) and short-term (1978–95) land-surface subsidence with published maps of land-surface subsidence for about the same period until acceptable matches were achieved. Simulated 1996 Chicot aquifer flow rates indicate that a net flow of 562.5 cubic feet per second enters the Chicot aquifer in the outcrop area, and a net flow of 459.5 cubic feet per second passes through the Chicot aquifer into the Evangeline aquifer. The remaining 103.0 cubic feet per second of flow is withdrawn as pumpage, with a shortfall of about 84.9 cubic feet per second supplied to the wells from storage in sands and clays. Water simulated from storage in clays in the Chicot aquifer is about 19 percent of the total water withdrawn from the aquifer. Simulated 1996 Evangeline aquifer flow rates indicate that a net flow of 14.8 cubic feet per second enters the Evangeline aquifer in the outcrop area, and a net flow of 459.5 cubic feet per second passes through the Chicot aquifer into the Evangeline aquifer for a total inflow of 474.3 cubic feet per second. A greater amount, 528.6 cubic feet per second, is withdrawn by wells; the shortfall of about 54.8 cubic feet per second is supplied from storage in sands and clays. Water simulated from storage in clays in the Evangeline aquifer is about 10 percent of the total water withdrawn from the aquifer.

Texas↗
USGS2006

Simulation of saltwater movement in the Upper Floridan aquifer in the Savannah, Georgia-Hilton Head Island, South Carolina, area, predevelopment-2004, and projected movement for 2000 pumping conditions

A digital model was developed to simulate ground-water flow and solute transport for the Upper Floridan aquifer in the Savannah, Georgia–Hilton Head Island, South Carolina, area. The model was used to (1) simulate trends of saltwater intrusion from predevelopment to the present day (1885–2004), (2) project these trends from the present day into the future, and (3) evaluate the relative influence of different assumptions regarding initial and boundary conditions and physical properties. The model is based on a regional, single-density ground-water flow model of coastal Georgia and adjacent parts of South Carolina and Florida. Variable-density ground-water flow and solute transport were simulated using the U.S. Geological Survey finite-element, variable-density solute-transport simulator SUTRA, 1885–2004. The model comprises seven layers: the surficial aquifer system, the Brunswick aquifer system, the Upper Floridan aquifer, the Lower Floridan aquifer, and the intervening confining units. The model was calibrated to September 1998 water levels, for single-density freshwater conditions, then refined using variable density and chloride concentration to give a reasonable match to the trend in the chloride distribution in the Upper Floridan aquifer inferred from field measurements of specific conductance made during 2000, 2002, 2003, and 2004. The model was modified to simulate solute transport by allowing saltwater to enter the system through localized areas near the northern end of Hilton Head Island, at Pinckney Island, and near the Colleton River, and was calibrated to match chloride concentrations inferred from field measurements of specific conductance. This simulation is called the "Base Case." Water-level residuals ranged from –5.3 to 23.4 feet for September 1998 conditions and single-density freshwater conditions. When chloride transport was simulated, water-level residuals ranged from –12.5 to 23.3 feet. The simulated chloride distribution captures the general trends in the field data. Chloride transport is sensitive to the permeabilities assigned to the confining units in the source areas and the porosity assigned to the Upper Floridan aquifer. Results of the study indicate that if present-day (year 2000) pumping conditions are maintained, plumes of saltwater in the Upper Floridan aquifer will continue to expand and move toward Savannah and across Hilton Head Island; the rate of movement of the 250-mg/L (milligram per liter) isochlor toward Savannah is between 144 feet per year and 190 feet per year and that the 250-mg/L isochlor could reach the pumping center at Savannah in 800 years; if effective porosities are lower than those used in the model, as is likely, higher rates of solute transport would result; and plumes may have occurred along the northern shore of Hilton Head Island before substantial development began in the mid-1960s, and lesser amounts of intrusion may have already occurred prior to the onset of pumping during 1885. Model limitations include uncertainty in (1) field data, (2) the conceptual model, (3) the physical properties and representation of the hydrogeologic framework, and (4) uncertainty in the boundary and initial conditions. Results of simulations projected far into the future must be interpreted with caution because they are based on an assumed future pumping distribution and fixed boundary conditions, and because these conditions may differ substantially from those for which the model is calibrated.

South Carolina↗