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 1,099 records · Page 61Linked to original sources
USGS2008

Simulated Effects of Year 2030 Water-Use and Land-Use Changes on Streamflow near the Interstate-495 Corridor, Assabet and Upper Charles River Basins, Eastern Massachusetts

Continued population growth and land development for commercial, industrial, and residential uses have created concerns regarding the future supply of potable water and the quantity of ground water discharging to streams in the area of Interstate 495 in eastern Massachusetts. Two ground-water models developed in 2002-2004 for the Assabet and Upper Charles River Basins were used to simulate water supply and land-use scenarios relevant for the entire Interstate-495 corridor. Future population growth, water demands, and commercial and residential growth were projected for year 2030 by the Metropolitan Area Planning Council. To assess the effects of future development on subbasin streamflows, seven scenarios were simulated by using existing computer-based ground-water-flow models with the data projected for year 2030. The scenarios incorporate three categories of projected 2030 water- and land-use data: (1) 2030 water use, (2) 2030 land use, and (3) a combination of 2030 water use and 2030 land use. Hydrologic, land-use, and water-use data from 1997 through 2001 for the Assabet River Basin study and 1989 through 1998 for the Upper Charles River Basin study were used to represent current conditions - referred to as 'basecase' conditions - in each basin to which each 2030 scenario was compared. The effects of projected 2030 land- and water-use change on streamflows in the Assabet River Basin depended upon the time of year, the hydrologic position of the subbasin in the larger basin, and the relative areas of new commercial and residential development projected for a subbasin. Effects of water use and land use on streamflow were evaluated by comparing average monthly nonstorm streamflow (base flow) for March and September simulated by using the models. The greatest decreases in streamflow (up to 76 percent in one subbasin), compared to the basecase, occurred in September, when streamflows are naturally at their lowest level. By contrast, simulated March streamflows decreased less than 6.5 percent from basecase streamflows in all subbasins for all scenarios. The simulations showed similar effects in the Upper Charles River Basin, but increased water use contributed to decreased simulated streamflow in most subbasins. Simulated changes in March streamflows for 2030 in the Upper Charles River Basin were within +- 6 percent of the basecase for all scenarios and subbasins. Percentage decreases in simulated September streamflows for 2030 were greater than in March but less than the September decreases that resulted for some subbasins in the Assabet River Basin. Only two subbasins of the Upper Charles River Basin had projected decreases greater than 5 percent. In the Mill River subbasin, the decrease was 11 percent, and in the Mine Brook subbasin, 6.6 percent. Changes in water use and wastewater return flow generally were found to have the greatest effect in the summer months when streamflow and aquifer recharge rates are low and water use is high. September increases in main-stem streamflow of both basins were due mainly to increased discharge of treated effluent from wastewater-treatment facilities on the main-stem rivers. In the Assabet River Basin, wastewater-treatment-facility discharge became a smaller proportion of total streamflow with distance downstream. In contrast, wastewater-treatment facility discharge in the Upper Charles River Basin became a greater proportion of streamflow with distance downstream. The effects of sewer-line extension and low-impact development on streamflows in two different subbasins of the Assabet River Basin also were simulated. The result of extending sewer lines with a corresponding decrease in septic-system return flow caused September streamflows to decrease as much as 15 percent in the Fort Pond Brook subbasin. The effect of low-impact development was simulated in the Hop Brook subbasin in areas projected for commercial development. In this simulation, the greater the area where low-i

Scientific Investigations Report↗
USGS1997

Evaluation of ground-water flow and solute transport in the Lompoc area, Santa Barbara County, California

Ground-water quality in the Lompoc area, especially in the Lompoc plain, is only marginally acceptable for most uses. Demand for ground water has increased for municipal use since the late 1950's and has continued to be high for irrigation on the Lompoc plain, the principal agricultural area in the Santa Ynez River basin. As use has increased, the quality of ground water has deteriorated in some areas of the Lompoc plain. The dissolved-solids concentration in the main zone of the upper aquifer beneath most of the central and western plains has increased from less than 1,000 milligrams per liter in the 1940's to greater than 2,000 milligrams per liter in the 1960's. Dissolved- solids concentration have remained relatively constant since the 1960's. A three-dimensional finite-difference model was used to simulate ground-water flow in the Lompoc area and a two-dimensional finite-element model was used to simulate solute transport to gain a better understanding of the ground-water system and to evaluate the effects of proposed management plans for the ground-water basin. The aquifer system was simulated in the flow model as four horizontal layers. In the area of the Lompoc plain, the layers represent the shallow, middle, and main zones of the upper aquifer, and the lower aquifer. For the Lompoc upland and Lompoc terrace, the four layers represent the lower aquifer. The solute transport model was used to simulate dissolved-solids transport in the main zone of the upper aquifer beneath the Lompoc plain. The flow and solute-transport models were calibrated to transient conditions for 1941-88. A steady-state simulation was made to provide initial conditions for the transient-state simulation by using long-term average (1941-88) recharge rates. Model- simulated hydraulic heads generally were within 5 feet of measured heads in the main zone for transient conditions. Model-simulated dissolved- solids concentrations for the main zone generally differed less than 200milligrams per liter from concentrations in 1988. During 1941-88 about 1,096,000 acre-feet of water was pumped from the aquifer system. Average pumpage for this period (22,830 acre-feet per year) exceeded pumpage for the steady-state simulation by 16,590 acre-feet per year. The results of the transient simulation indicate that about 60 percent of this increase in pumpage was contributed by increased recharge, 28 percent by decreased natural discharge from the system (primarily discharge to the Santa Ynez River and transpiration), and 13 percent was withdrawn from storage. Total simulated downward leakage from the middle zone to the main zone in the central plain and upward leakage from the consolidated rocks to the main zone significantly increased in response to increased pumpage, which increased from about 6,240 to 30,870 acre-feet per year from 1941 to 1988. Average dissolved-solid concentration in the middle zone in 1987-88 ranged from 2,000 to 3,000 milligrams per liter beneath the northeastern plain and the dissolved-solids concentration of two samples from the consolidated rocks beneath the western plain averaged 4,300 milligrams per liter. Because the dissolved-solids concentration for the middle zone and the consolidated rocks is higher than the simulated steady-state dissolved-solids concentration of the main zone, the increase in the leakage from these two sources resulted in increased dissolved-solids concentration in the main zone during the transient period. The model results indicate that the main source of increased dissolved- solids concentration in the northeastern and central plains was downward leakage from the middle zone; whereas, upward leakage from the consolidated rocks was the main source of the increased dissolved-solids concentrations in the northwestern and western plains. The models were used to estimate changes in hydraulic head and in dissolved-solids concentration resulting from three proposed management alternatives: (1) average recharge

Water-Resources Investigations Report↗
USGS2010

Estimation of the effects of land use and groundwater withdrawals on streamflow for the Pomperaug River, Connecticut

A precipitation runoff model for the Pomperaug River watershed, Connecticut was developed to address issues of concern including the effect of development on streamflow and groundwater recharge, and the implications of water withdrawals on streamflow. The model was parameterized using a strategy that requires a minimum of calibration and optimization by establishing basic relations between the parameter value and physical characteristics of individual hydrologic response units (HRUs) that comprise the model. The strategy was devised so that the information needed can be obtained from Geographic Information System and other general databases for Connecticut. Simulation of groundwater recharge enabled evaluation of the temporal and spatial mapping of recharge variation across the watershed and the spatial effects of changes in land cover on base flow and surface runoff. The modeling indicated that over the course of a year, groundwater provides between 60 and 70 percent of flow in the Pomperaug River; the remainder is generated by more rapid flow through the shallow subsurface and runoff from impermeable surfaces and saturated ground. Groundwater is recharged primarily during periods of low evapotranspiration in the winter, spring, and fall. The largest amount of recharge occurs in the spring in response to snowmelt. During floods, the Weekeepeemee and Nonnewaug Rivers (tributaries that form the Pomperaug River) respond rapidly with little flood peak attenuation due to flood-plain storage. In the Pomperaug River, flood-plain storage is more important in attenuating floods; abandoned quarry ponds (O&G ponds) adjacent to the river provide substantial flood storage above specific river stages when flow from the river spills over the banks and fills the ponds. Discharge from the ponds also helps to sustain low flows in the Pomperaug River. Similarly, releases from the Bronson-Lockwood reservoir sustain flow in the Nonnewaug River and tend to offset the effect of groundwater withdrawals from a well field adjacent to the river during periods of natural low flow. The model indicated that under the current zoning, future development could reduce low flows by as much as 10 percent at the 99 percent exceedance level (99 percent of flows are greater than or equal to this flow), but would not substantially increase the highest flows. Simulation of projected and hypothetical development in the watershed shows, depending on how stormwater is managed, that between 10 and 20 percent effective impervious area in an HRU results in streamflow becoming dominated by the surface-runoff component. This shift from a groundwater-dominated system would likely result in substantial changes in water quality and instream habitat characteristics of the river. Base flow to streams in the Pomperaug River watershed is reduced by both increased impervious surface and increased groundwater withdrawals. For the watershed as a whole, increasing groundwater withdrawals have the potential for causing greater overall reductions in flow compared to increased development and impervious surfaces. Additionally, on the basis of groundwater-modeling simulations, the projected increase in development across the watershed and, to a lesser extent the increase in groundwater withdrawals, will increase the number of local losing reaches experiencing dry conditions and the duration of these dry periods. The location of the losing reaches tends to be in areas near the transition from the uplands to the valley bottoms that are filled with coarse glacial stratified deposits. The simulated increase in the duration and extent of localized dry stream reaches is most sensitive to local increase in impervious surface. Conversion of land from forest or developed land cover to pasture or agricultural land increases groundwater recharge and discharge to streams, while at the same time increasing overall streamflow (the opposite effect as increased impervious surface). These resu

Scientific Investigations Report↗
USGS1997

Hydrogeology and water quality of a surficial aquifer underlying an urban area, Manchester, Connecticut

The quality of water along flowpaths in a surficial aquifer system in Manchester, Connecticut, was studied during 1993-95 as part of the National Water Quality Assessment program. The flowpath study examined the relations among hydrogeology, land-use patterns, and the presence of contaminants in a surficial aquifer in an urban area, and evaluated ground water as a source of contamination to surface water. A two-dimensional, finite-difference groundwater- flow model was used to estimate travel distance, which ranged from about 50 to 11,000 feet, from the source areas to the sampled observation wells. Land use, land cover, and population density were determined in the source areas delineated by the ground-water-flow simulation. Source areas to the wells contained either high- or medium-density residential areas, and population density ranged from 629 to 8,895 people per square mile. Concentrations of selected inorganic constituents, including sodium, chloride, and nitrite plus nitrate nitrogen, were higher in the flowpath study wells than in wells in undeveloped areas with similar aquifer materials. One or more of 9 volatile organic compounds were detected at 12 of 14 wells. The three most commonly detected volatile organic compounds were chloroform, methyl-tert-butyl ether, and trichloroethene. Trichloroethene was detected at concentrations greater than the maximum contaminant level for drinking water (5 micrograms per liter) in samples from one well. Four pesticides, including dichloro diphenyl dichloroethylene, dieldrin, dichloroprop, and simazine were detected at low concentrations. Concentrations of sodium and chloride were higher in samples collected from wells screened in the top of the saturated zone than in samples collected from deeper zones. Volatile organic compounds and elevated concentrations of nitrite plus nitrate as nitrogen were detected at depths of as much as 60 feet below the water table, indicating that the effects of human activities on the ground-water quality extends to the bottom of the surficial aquifer. The age of ground water, as determined by tritium and 3helium concentrations, was 0.9 to 22.6 years. pH, alkalinity, and calcium were higher and concentrations of dissolved oxygen were lower in ground-water samples with ages of 10 years or more than in samples younger than 10 years. In addition, concentrations of sodium, chloride, and nitrite plus nitrate nitrogen were low in ground-water samples with ages of 10 years or more, indicating that concentrations of these compounds may be increasing with time or that the recharge areas to these wells may have had less intensive urban land use. Methyl-tert-butyl ether was detected only in wells with ground water ages of less than 11 years, which is consistent with the date of introduction of this compound as a gasoline additive in Connecticut. Analysis of additional samples collected for analysis of stable nitrogen isotopes indicated that the most likely source of elevated concentrations of nitrate nitrogen was lawn and garden fertilizers, but other sources, including wastewater effluents, soil organic nitrogen, and atmospheric deposition, may contribute to the total. Population density was positively correlated (at the 97 percent confidence level) to concentrations of nitrite plus nitrate as nitrogen. Water quality in the Hockanum River aquifer has been degraded by human activities, and, after discharge to surface water, affects the water quality in the Hockanum River. On an annual basis, ground-water discharge from the study area to the river (as measured at a downstream continuous-record gaging station) contributes about 5 percent of the annual load of nitrite plus nitrate nitrogen, but, during low flow, contributes 11 percent of the nitrite plus nitrate nitrogen, 32 percent of the calcium, and 16 percent of the chloride to the river.

Connecticut↗
USGS1998

Analytical methods, numerical modeling, and monitoring strategies for evaluating the effects of ground-water withdrawals on unconfined aquifers in the New Jersey Coastal Plain

Analytical and numerical solutions of ground-water withdrawals in the unconfined part of the Kirkwood-Cohansey aquifer system of the Coastal Plain of New Jersey were evaluated for their usefulness in predicting the area of influence of a pumped well and in determining hydraulic characteristics of an aquifer. Additionally, simulations of ground-water withdrawal using a finite-difference model provided information on the ways in which prudent well-location strategies can disperse the local effects of withdrawal over a larger part of an aquifer system. The design of a monitoring network that is sensitive to the ground-water hydraulics of streams and wetlands of the Coastal Plain of New Jersey also was considered for its utility in providing hydrologic data necessary to establish the baseline hydrologic conditions near wetlands and streams and in signaling when ground-water levels are being adversely affected by withdrawals elsewhere in the system. The application of methods based on the Theis analytical solution to ground-water flow in unconfined aquifers can lead to erroneous estimates of the size of the area of influence generated by ground-water withdrawals. Analysis oftime-drawdown data from an unconfined aquifer system are best evaluated by means of the Neuman solution, which accounts for the effects of gravity drainage; however, the pumped well must be far enough from streams so that ground water is not drawn from nearby streams. Time-drawdown data from a test well in Winslow County, N.J., were analyzed by means of the Neuman solution. Results indicate that the aquifer has a relatively high vertical to horizontal anisotropy of 1:198, and a specific yield of 0.028, an indication that the area of influence of a pumped well at the test site would be relatively large. Results from a finite-difference ground-water-flow model of the northeastern part of the Mullica River Basin near Chesilhurst, N.J., show that the area influenced by a long-term withdrawal is best estimated from a steady state ground-water-flow analysis that includes the effects of average areal recharge. Withdrawal simulations indicate an order-of-magnitude difference between the size of the area of influence generated from a 3-day (72 hour) withdrawal and the size of the area produced under steady-state conditions. An aquifer characterized by a low specific yield will cause the area of influence to extend farther away from the pumped well. The contributing area of flow to the pumped well includes areas on the water table that would, under natural conditions, be incorporated into the contributing areas of flow to streams. Ground water that is drawn to a pumped well is diverted from nearby streams; the withdrawal decreases the size of the contributing areas of flow to streams by an amount equal to the contributing area of flow to the well. Withdrawals made from a well close to a stream divert ground water that would, under natural conditions, flow to the stream. The diverted ground water causes the area of influence of the well to be smaller than it would if the well were far from the stream. Water-table declines caused by withdrawals near streams are, to some degree, mitigated by ground-water diversion from streams. However, the withdrawals can significantly reduce ground-water seepage to nearby streams, especially along stream reaches and wetlands close to the well. Alternatively, these effects can be dispersed over a large part of the aquifer if wells are located on surface-water divides. Measurements of seasonal water-level fluctuations in the Mullica River Basin indicate that the greatest fluctuations in water levels are found in upland areas, where the average fluctuation is 3.4 feet. Fluctuations in hydraulic head in the wetland areas averages 1.3 feet. The bimodal average of ranges in water levels show that upland areas are more sensitive to recharge than lowland areas. The pattern of yearly mean water levels fluctuates irregularly about a long-term mean value. Abnormally low or high yearly average values that are brought on by periods of drought or excess recharge are short lived; over time, hydrologic conditions shift back to average levels under natural conditions. Wetland areas in the New Jersey Coastal Plain are characterized by ground-water seepage into wide, shallow depressions. Periods of inundation are longest in the deepest part of the depression, whereas inundation of areas near the fringes of wetlands due to ground-water seepage is only seasonal. The seepage face in the fringe areas expand and contract in response to seasonal variation in water-table elevation and in response to precipitation. Values of the aquifer storage coefficient and transmissivity can, in some cases, be determined by use of hydraulic head or streamflow recession analysis as an alternative to aquifer testing. The recession curves developed from hydro graphs of Middle Branch and McDonalds Branch in the New Jersey Coastal Plain indicate that the aquifer near McDonalds Branch has about 2.6 times the storage capacity of the aquifer adjacent to Middle Branch; this finding is consistent with the relatively small ranges of water-level changes measured in McDonalds Branch compared to those measured in Middle Branch.

New Jersey↗
USGS2008

Simulation of Ground-Water Flow and Effects of Ground-Water Irrigation on Base Flow in the Elkhorn and Loup River Basins, Nebraska

Irrigated agriculture is vital to the livelihood of communities in the Elkhorn and Loup River Basins in Nebraska, and ground water is used to irrigate most of the cropland. Concerns about the sustainability of ground-water and surface-water resources have prompted State and regional agencies to evaluate the cumulative effects of ground-water irrigation in this area. To facilitate understanding of the effects of ground-water irrigation, a numerical computer model was developed to simulate ground-water flow and assess the effects of ground-water irrigation (including ground-water withdrawals, hereinafter referred to as pumpage, and enhanced recharge) on stream base flow. The study area covers approximately 30,800 square miles, and includes the Elkhorn River Basin upstream from Norfolk, Nebraska, and the Loup River Basin upstream from Columbus, Nebraska. The water-table aquifer consists of Quaternary-age sands and gravels and Tertiary-age silts, sands, and gravels. The simulation was constructed using one layer with 2-mile by 2-mile cell size. Simulations were constructed to represent the ground-water system before 1940 and from 1940 through 2005, and to simulate hypothetical conditions from 2006 through 2045 or 2055. The first simulation represents steady-state conditions of the system before anthropogenic effects, and then simulates the effects of early surface-water development activities and recharge of water leaking from canals during 1895 to 1940. The first simulation ends at 1940 because before that time, very little pumpage for irrigation occurred, but after that time it became increasingly commonplace. The pre-1940 simulation was calibrated against measured water levels and estimated long-term base flow, and the 1940 through 2005 simulation was calibrated against measured water-level changes and estimated long-term base flow. The calibrated 1940 through 2005 simulation was used as the basis for analyzing hypothetical scenarios to evaluate the effects of ground-water irrigation on stream base flow for 1940 through 2005 and for 2006 through 2045. Simulated base flows were compared for scenarios that alternately did or did not include a representation of the effects of ground-water irrigation. The difference between simulated base flows for the two scenarios represents the predicted effects of ground-water irrigation on base flow. Comparison of base flows between simulations with ground-water irrigation and no ground-water irrigation indicated that ground-water irrigation has cumulatively reduced streamflows from 1940 through 2005 by 888,000 acre-feet in the Elkhorn River Basin and by 2,273,000 acre-feet in the Loup River Basin. Generally, predicted cumulative effects of ground-water irrigation on base flow were 5 to 10 times larger from 2006 through 2045 than from 1940 through 2005, and were 7,678,000 acre-feet for the Elkhorn River Basin and 14,784,000 acre-feet for the Loup River Basin. The calibrated simulation also was used to estimate base-flow depletion as a percentage of pumping volumes for a 50-year future time period, because base-flow depletion percentages are used to guide the placement of management boundaries in Nebraska. Mapped results of the base-flow depletion analysis conducted for most of the interior of the study area indicated that pumpage of one additional theoretical well simulated for a future 50-year period generally would result in more than 80 percent depletion when it was located close to the stream, except in areas where depletion was partly offset by reduced ground-water discharge to evapotranspiration in wetland areas. In many areas, depletion for the 50-year future period composed greater than 10 percent of the pumped water volume for theoretical wells placed less than 7 or 8 miles from the stream, though considerable variations existed because of the heterogeneity of the natural system represented in the simulation. For a few streams, predicted future simulated base flows dec

Scientific Investigations Report↗
USGS1978

The hydrothermal system of Long Valley Caldera, California

Long Valley caldera, an elliptical depression covering 450 km 2 on the eastern front of the Sierra Nevada in east-central California, contains a hot-water convection system with numerous hot springs and measured and estimated aquifer temperatures at depths of 180?C to 280?C. In this study we have synthesized the results of previous geologic, geophysical, geochemical, and hydrologic investigations of the Long Valley area to develop a generalized conceptual and mathematical model which describes the gross features of heat and fluid flow in the hydrothermal system. Cenozoic volcanism in the Long Valley region began about 3.2 m.y. (million years) ago and has continued intermittently until the present time. The major event that resulted in the formation of the Long Valley caldera took place about 0.7 m.y. ago with the eruption of 600 km 3 or more of Bishop Tuff of Pleistocene age, a rhyolitic ash flow, and subsequent collapse of the roof of the magma chamber along one or more steeply inclined ring fractures. Subsequent intracaldera volcanism and uplift of the west-central part of the caldera floor formed a subcircular resurgent dome about 10 km in diameter surrounded by a moat containing rhyolitic, rhyodacitic, and basaltic rocks ranging in age from 0.5 to 0.05 m.y. On the basis of gravity and seismic studies, we estimate an aver- age thickness of fill of 2.4 km above the precaldera granitic and metamorphic basement rocks. A continuous layer of densely welded Bishop Tuff overlies the basement rocks, with an average thickness of 1.4 km; the fill above the welded Bishop Tuff consists of intercalated volcanic flows and tuffs and fluvial and lacustrine deposits. Assuming the average grain density of the fill is between 2.45 and 2.65 g/cm 3 , we calculate the average bulk porosity of the total fill as from 0.11 to 0.21. Comparison of published values of porosity of the welded Bishop Tuff exposed southeast of the caldera with calculated values indicates average bulk porosity for the welded tuff (including fracture porosity) from 0.05 to 0.10. Because of its continuity and depth and the likelihood of significant fracture permeability in the more competent rocks such as the welded tuff, our model of the hydrothermal system assumes that the Bishop Tuff provides the principal hot-water reservoir. However, because very little direct information exists from drill holes below 300 m, this assumption must be considered tentative. Long Valley caldera is drained by the Owens River and several tributaries which flow into Lake Crowley in the southeast end of the caldera. Streamflow and springflow measurements for water years 1964-74 indicate a total inflow to Lake Crowley of about 10,900 L/s. In contrast, the total discharge of hot water from the hydrothermal reservoir is about 300 L/s. For modeling purposes, the ground-water system is considered as comprising a shallow subsystem in the fill above the densely welded Bishop Tuff containing relatively cold ground water, and a deep subsystem or hydrothermal reservoir in the welded tuff containing relatively hot ground water. Hydrologic, isotopic, and thermal data indicate that recharge to the hydrothermal reservoir occurs in the upper Owens River drainage basin along the western periphery of the caldera. Temperature profiles in a 2.11- km-deep test well drilled by private industry in the southeastern part of the caldera suggest that an additional flux of relatively cool ground water recharges the deep subsystem around the northeast rim. Flow in the shallow ground-water subsystem is neglected in the model except in recharge areas and along Hot Creek gorge, where approximately 80 percent of the hot-water discharge from the hydrothermal reservoir moves upward along faults toward springs in the gorge. Heat-flow data from the Long Valley region indicate that the resurgent dome overlies a residual magma chamber more circular in plan than the original magma chamber that supplied the Bishop Tuff

California↗
USGS2002

Ground-water flow and numerical simulation of recharge from streamflow infiltration near Pine Nut Creek, Douglas County, Nevada

Ground-water flow and recharge from infiltration near Pine Nut Creek, east of Gardnerville, Nevada, were simulated using a single-layer numerical finite-difference model as part of a study made by the U.S. Geological Survey in cooperation with the Carson Water Subconservancy District. The model was calibrated to 190 water-level measurements made in 27 wells in December 2000, and in 9 wells from August 1999 through April 2001. The purpose of this study was to estimate reasonable limits for the approximate volume of water that may be stored by recharge through infiltration basins, and the rate at which recharged water would dissipate or move towards the valley floor. Measured water levels in the study area show that infiltration from the Allerman Canal and reservoir has created a water-table mound beneath them that decreases the hydraulic gradient east of the canal and increases the gradient west of the canal. North of Pine Nut Creek, the mound causes ground water to flow toward the northern end of the reservoir. South of Pine Nut Creek, relatively high water levels probably are maintained by the mound beneath the Allerman Canal and possibly by greater rates of recharge from the southeast. Water-level declines near Pine Nut Creek from August 1999 through April 2001 probably are caused by dissipation of recharge from infiltration of Pine Nut Creek streamflow in the springs of 1998 and 1999. Using the calibrated model, a simulation of recharge through a hypothetical infiltration basin covering 12.4 acres near Pine Nut Creek applied 700 acre-feet per year of recharge over a six-month period, for a total of 3,500 acre-feet after 5 consecutive years. This recharge requires a diversion rate of about 2 cubic feet per second and an infiltration rate of 0.3 foot per day. The simulations showed that recharge of 3,500 acre-feet caused water levels near the basin to rise over 70 feet, approaching land surface, indicating 3,500 acre-feet is the maximum that may be stored in a 5-year period, given the basin location and surface area used in the simulations. Greater amounts probably could be stored if separate infiltration basins were installed at different locations along the Pine Nut Creek alluvial fan, applying the recharge over a larger area. The water-table mound resulting from recharge extended 7,000 feet north, west, and south of the infiltration basin. After recharge ceased, water levels near the center of the mound declined rapidly to within 20 feet of initial levels after 2 years, and within 10 feet of initial levels after 7 years. The recharge mound dissipates laterally across the modeled area at decreasing rates over time. A water-level rise of 1 foot moved westward towards the valley floor 660 feet from peak conditions after 1 year, and averaged 550 feet, 440 feet, and 330 feet per year for the periods 1-4, 4-7, and 7-10 years, respectively, after recharge ceased. Simulations of subsequent pumping from hypothetical wells near the infiltration basin were made by applying pumping near the basin beginning 1 year after recharge of 3,500 acre-feet ceased. Pumping was applied over a 6-month period for 4 years from one well at 400 acre-feet per year, withdrawing 1,600 acre-feet or 45 percent of that recharged, and from two wells totaling 800 acre-feet per year, withdrawing 3,200 acre-feet or 90 percent of that recharged. Pumping of 1,600 acre-feet caused water-levels near the infiltration basin to decline only slightly below initial levels. Pumping of 3,200 acre-feet caused water-levels near the infiltration basin to decline a maximum of 30 feet below initial levels, with smaller declines extending laterally in all directions for 4,000 feet from the pumping wells. Water-level declines are a result of pumping at a rate sufficient to withdraw the majority of the water recharged through the infiltration basin. Although the declines may affect water levels in nearby domestic wells, the simulations show that water levels recover quickly after

Water-Resources Investigations Report↗
USGS1998

Ground-water hydrology and simulation of ground-water flow at Operable Unit 3 and surrounding region, U.S. Naval Air Station, Jacksonville, Florida

The Naval Air Station, Jacksonville (herein referred to as the Station), occupies 3,800 acres adjacent to the St. Johns River in Duval County, Florida. Operable Unit 3 (OU3) occupies 134 acres on the eastern side of the Station and has been used for industrial and commercial purposes since World War II. Ground water contaminated by chlorinated organic compounds has been detected in the surficial aquifer at OU3. The U.S. Navy and U.S. Geological Survey (USGS) conducted a cooperative hydrologic study to evaluate the potential for ground water discharge to the neighboring St. Johns River. A ground-water flow model, previously developed for the area, was recalibrated for use in this study. At the Station, the surficial aquifer is exposed at land surface and forms the uppermost permeable unit. The aquifer ranges in thickness from 30 to 100 feet and consists of unconsolidated silty sands interbedded with local beds of clay. The low-permeability clays of the Hawthorn Group form the base of the aquifer. The USGS previously conducted a ground-water investigation at the Station that included the development and calibration of a 1-layer regional ground-water flow model. For this investigation, the regional model was recalibrated using additional data collected after the original calibration. The recalibrated model was then used to establish the boundaries for a smaller subregional model roughly centered on OU3. Within the subregional model, the surficial aquifer is composed of distinct upper and intermediate layers. The upper layer extends from land surface to a depth of approximately 15 feet below sea level; the intermediate layer extends from the upper layer down to the top of the Hawthorn Group. In the northern and central parts of OU3, the upper and intermediate layers are separated by a low-permeability clay layer. Horizontal hydraulic conductivities in the upper layer, determined from aquifer tests, range from 0.19 to 3.8 feet per day. The horizontal hydraulic conductivity in the intermediate layer, determined from one aquifer test, is 20 feet per day. An extensive stormwater drainage system is present at OU3 and the surrounding area. Some of the stormwater drains have been documented to be draining ground water from the upper layer of the surficial aquifer, whereas other drains are only suspected to be draining ground water. The subregional model contained 78 rows and 148 columns of square model cells that were 100 feet on each side. Vertically, the surficial aquifer was divided into two layers; layer 1 represented the upper layer and layer 2 represented the intermediate layer. Steady-state ground-water flow conditions were assumed. The model was calibrated to head data collected on October 29 and 30, 1996. After calibration, the model matched all 67 measured heads to within the calibration criterion of 1 foot; and 48 of 67 simulated heads (72 percent) were within 0.5 foot. Model simulated recharge rates ranged from 0.4 inch per year in areas that were largely paved to 13.0 inches per year in irrigated areas. Simulated hydraulic conductivities in the upper layer at OU3 ranged from 0.5 foot per day in the north to 1.0 foot per day in the south. Simulated vertical leakance between the upper and intermediate layers ranged from 1.0x10-6 per day in an area with low-permeability clays to 4.3x10-2 per day in an area that had been dredged. Simulated transmissivities in the intermediate layer ranged from 25 feet squared per day in an area of low-permeability channel-fill deposits to a high of 1,200 feet squared per day in areas covering most of OU3. Simulated riverbed conductances ranged from 4 to 60 feet squared per day and simulated bottom conductances of leaking stormwater drains ranged from 5 to 20 feet squared per day. The direction and velocity of ground-water flow was determined using particle-tracking techniques. Ground-water flow in the upper layer was generally eastward toward the St. Johns River. However, leaking stormwater drains locally modified the flow system to create small areas with flow that was diverted to the drains. The flow velocities in the upper layer at OU3 were slow, averaging about 2 feet per year. The slow velocities were primarily the result of the low horizontal hydraulic conductivity and, secondarily, the result of the low recharge rate. The simulated rate at which ground water leaked into the stormwater drains was low, averaging about 0.0011 cubic feet per second per 100 feet of stormwater drainage conduit. Ground-water flow in the intermediate layer moved eastward toward and discharged into the St. Johns River. Flow velocities were significantly higher in this layer than in the upper layer. The velocity was about 35 and 12 feet per year in the northern and southern parts of OU3, respectively.

Florida↗
USGS2004

Hydrogeology, water quality, and ecology of Anderton Branch near the Quail Hollow Landfill, Bedford County, Tennessee, 1995-99

The Quail Hollow Landfill, located in southeastern Bedford County on the Highland Rim overlooking the Central Basin karst region of Tennessee, is constructed on the gravelly, clay-rich residuum of the Fort Payne Formation of Mississippian age. A conceptual hydrologic model of the landfill indicated that Anderton Branch was at risk of being affected by the landfill. Ground water flowing beneath the landfill mixes with percolating rainwater that has passed through the landfill and discharges to the surface from numerous weeps, seeps, and springs present in the area. Anderton Branch, adjacent to the landfill site on the north and east, receives most of the discharge from these weeps, seeps, and springs. Anderton Branch also receives water from the Powell Branch drainage basin to the west and south because of diverted flow of ground water through Harrison Spring Cave. The U.S. Geological Survey, in cooperation with the Bedford County Solid Waste Authority, conducted a study to evaluate the effect of the Quail Hollow Landfill on ground- and surface-water quality. During storm runoff, specific conductance was elevated, and cadmium, iron, manganese, lead, and nickel concentrations in Anderton Branch frequently exceeded maximum contaminant levels for drinking water for the State of Tennessee. High chloride inputs to Anderton Branch were detected at two locations?a barnyard straddling the stream and a tributary draining a pond that receives water directly from the landfill. The chloride inputs probably contribute to chloride load levels that are three times higher for Anderton Branch than for the control stream Anthony Branch. Although toxic volatile organic compounds were detected in water from monitoring wells at the landfill, no organic contaminants were detected in domestic water wells adjacent to the landfill or in Anderton Branch. Sons Spring, a karst spring near the landfill, has been affected by the landfill as indicated by an increase in chloride concentrations from 4 milligrams per liter in 1974 to 59 milligrams per liter in 1996. Analysis of water samples from Sons Spring detected concentrations of nickel that exceeded primary drinking-water standards and Tennessee Department of Environment and Conservation fish and aquatic life chronic standards. Trichloroethene, 1,1-dichloroethene, and 1,1-dichloroethane also were detected at Sons Spring. The presence of these chlorinated solvents imply the landfill origin of the contaminants in Sons Spring. Continuous monitoring at Sons Spring indicated a pattern of decreased specific conductance and lower contaminant concentrations after a storm. Contaminant concentrations increased with specific conductance to pre-storm levels after several days. The benthic macroinvertebrate community in Anderton Branch adjacent to the landfill was not different from the communities at control sites upstream and in Anthony Branch. Sons Spring, however, has low abundance and numbers of benthic macroinvertebrate taxa. Toxicity studies using Ceriodaphnia dubia indicated no toxicity in the base flow or storm water in Anderton Branch or in a tributary draining a pond that receives water from the landfill and Sons Spring; however, water collected from Sons Spring resulted in 100 percent mortality to all organisms within 48 hours. High concentrations of nickel were detected in crayfish tissue from control sites and Anderton Branch. Analysis of sediment samples also indicates nickel concentrations are high at control sites upstream of the landfill. Increased levels of the biomarker metallothionein detected in crayfish from Anderton Branch likely are not caused by nickel or cadmium because the levels present in the tissue are not correlated with metallothionein levels. Despite the high levels of certain metals in Anderton Branch during storm flow, the lack of toxicity and the health of the benthic community imply no detectable negative effect from the landfill to the stream. Sons Spring, howe

Scientific Investigations Report↗
USGS2009

Simulations of Groundwater Flow and Particle Tracking Analysis in the Area Contributing Recharge to a Public-Supply Well near Tampa, Florida, 2002-05

Shallow ground water in the north-central Tampa Bay region, Florida, is affected by elevated nitrate concentrations, the presence of volatile organic compounds, and pesticides as a result of groundwater development and intensive urban land use. The region relies primarily on groundwater for drinking-water supplies. Sustainability of groundwater quality for public supply requires monitoring and understanding of the mechanisms controlling the vulnerability of public-supply wells to contamination. A single public-supply well was selected for intensive study based on the need to evaluate the dominant processes affecting the vulnerability of public-supply wells in the Upper Floridan aquifer in the City of Temple Terrace near Tampa, Florida, and the presence of a variety of chemical constituents in water from the well. A network of 29 monitoring wells was installed, and water and sediment samples were collected within the area contributing recharge to the selected public-supply well to support a detailed analysis of physical and chemical conditions and processes affecting the water chemistry in the well. A three-dimensional, steady-state groundwater flow model was developed to evaluate the age of groundwater reaching the well and to test hypotheses on the vulnerability of the well to nonpoint source input of nitrate. Particle tracking data were used to calculate environmental tracer concentrations of tritium and sulfur hexafluoride and to calibrate traveltimes and compute flow paths and advective travel times in the model area. The traveltime of particles reaching the selected public-supply well ranged from less than 1 day to 127.0 years, with a median of 13.1 years; nearly 45 percent of the simulated particle ages were less than about 10 years. Nitrate concentrations, derived primarily from residential/commercial fertilizer use and atmospheric deposition, were highest (2.4 and 6.11 milligrams per liter as nitrogen, median and maximum, respectively) in shallow groundwater from the surficial aquifer system and lowest (less than the detection level of 0.06 milligram per liter) in the deeper Upper Floridan aquifer. Denitrification occurred near the interface of the surficial aquifer system and the underlying intermediate confining unit, within the intermediate confining unit, and within the Upper Floridan aquifer because of reducing conditions in this part of the flow system. However, simulations indicate that the rapid movement of water from the surficial aquifer system to the selected public-supply well through karst features (sinkholes) and conduit layers that bypass the denitrifying zones (short-circuits), coupled with high pumping rates, allow nitrate to reach the selected public-supply well in concentrations that resemble those of the overlying surficial aquifer system. Water from the surficial aquifer system with elevated concentrations of nitrate and low concentrations of some volatile organic compounds and pesticides is expected to continue moving into the selected public-supply well, because calculated flux-weighted concentrations indicate the proportion of young affected water contributing to the well is likely to remain relatively stable over time. The calculated nitrate concentration in the selected public-supply well indicates a lag of 1 to 10 years between peak concentrations of nonpoint source contaminants in recharge and appearance in the well.

Scientific Investigations Report↗
USGS2007

Simulated effects of projected 2010 withdrawals on ground-water flow and water levels in the New Jersey coastal plain – A task of the New Jersey Water Supply Plan, 2006 revision

A ground-water flow model previously developed as part of a Regional Aquifer System Analysis (RASA) of the New Jersey Coastal Plain was used to simulate ground-water flow in eight major confined aquifers to help evaluate ground-water resources in support of the New Jersey Department of Environmental Protection's revision of the New Jersey State Water Supply Plan. This model was calibrated to 1998 steady-state and transient conditions. Withdrawals at wells in operation in 1998 were varied in three scenarios to evaluate their effects on flow directions, water levels, and water budgets in the confined aquifers. The scenarios used to predict changes in pumpage from 1998 to 2010 were based on (1) a continuation of 1990-99 trends in water use, (2) public-supply withdrawals estimated from county population projections, and (3) restricted withdrawals in Water-Supply Critical Areas. Total withdrawals in these three scenarios were approximately 366, 362, and 355 million gallons per day, respectively. The results of these simulations are used by New Jersey water-management officials to help address water-supply concerns for the State. In the revision of the New Jersey State Water Supply Plan, the eight major confined aquifers of the New Jersey Coastal Plain and their outcrop areas are divided into 41 hydrologic budget areas (HBAs). Simulation results were used to assess the effects of changing ground-water withdrawals on water levels and the flow budgets in each budget area. Simulation results for each scenario were compared with 1998 (baseline) simulated water levels and flow budgets. The 41 hydrologic budget areas are in areas of large ground-water withdrawals, water-level declines, and (or) saltwater-intrusion potential. Their boundaries are based on various hydrologic, geohydrologic, and withdrawal conditions, such as aquifer extent, location of the 250-milligram-per-liter isochlor, aquifer outcrop area, and ground-water divides. The budget areas include primarily the onshore, freshwater portions of the aquifers. A budget analysis was done for each of the hydrologic budget areas for each scenario. Ground-water withdrawals, leakage to streams, net leakage to overlying and underlying aquifers, lateral flow to adjacent budget areas, and the flow direction at the 250-milligram-per-liter isochlor were evaluated. Although three different methods were applied to predict future pumping rates, the simulated water levels for scenarios 1 and 2 were generally within 2 feet of each other in most areas in the confined aquifers, but differences of more than 2 feet occurred locally. Differences in values of flow-budget components between scenarios 1 and 2 as a percentage change from 1998 values were generally within 2 percent in most hydrologic budget areas, but values of some budget components in some hydrologic budget areas differed by more than 2 percent. Simulated water levels recovered as much as 4 feet more in northeastern Camden and northwestern Burlington Counties in the Lower Potomac-Raritan-Magothy aquifer, and as much as 3 feet more in the same area in the Upper and Middle Potomac-Raritan-Magothy aquifers when pumpage restrictions were imposed in Critical Area 2 (scenario 3). In the Wenonah-Mount-Laurel aquifer, water levels declined continually in Monmouth County (HBA 8) downdip from the outcrop (in Critical Area 1) from 1988 to 2010 in all three scenarios, although most of the water levels farther downdip from this area in Critical Area 1 are still recovering because of mandated reductions in pumpage in the 1990s. In the Englishtown aquifer system, water levels declined continually in small areas in HBA 13 in central Monmouth County (in Critical Area 1) and in western Monmouth County downdip from the outcrop from 1988 to 2010 in all three scenarios, although most of the water levels farther downdip from this area are still recovering because of the mandated reductions in pumpage. In the Upper Potomac-Raritan-Magothy aquif

New Jersey↗
USGS2004

Geochemical characterization of ground-water flow in the Santa Fe Group aquifer system, Middle Rio Grande Basin, New Mexico

Chemical and isotopic data were obtained from ground water and surface water throughout the Middle Rio Grande Basin (MRGB), New Mexico, and supplemented with selected data from the U.S. Geological Survey (USGS) National Water Information System (NWIS) and City of Albuquerque water-quality database in an effort to refine the conceptual model of ground-water flow in the basin. The ground-water data collected as part of this study include major- and minor-element chemistry (30 elements), oxygen-18 and deuterium content of water, carbon-13 content and carbon-14 activity of dissolved inorganic carbon, sulfur-34 content of dissolved sulfate, tritium, and dissolved atmospheric gases including nitrogen, argon, helium, chlorofluorocarbons, and sulfur hexafluoride from 288 wells and springs in parts of the Santa Fe Group aquifer system. The surface-water data collected as part of this study include monthly measurements of major- and minor-element chemistry (30 elements), oxygen-18 and deuterium content of water, chlorofluorocarbons, and tritium content at 14 locations throughout the basin. Additional data include stable isotope analyses of precipitation and of ground water from City of Albuquerque production wells collected and archived from the early 1980?s, and other data on the chemical and isotopic composition of air, unsaturated zone air, plants, and carbonate minerals from throughout the basin. The data were used to identify 12 sources of water to the basin, map spatial and vertical extents of ground-water flow, map water chemistry in relation to hydrogeologic, stratigraphic, and structural properties of the basin, determine radiocarbon ages of ground water, and reconstruct paleo-environmental conditions in the basin over the past 30,000 years. The data indicate that concentrations of most elements and isotopes generally parallel the predominant north to south direction of ground-water flow. The radiocarbon ages of dissolved inorganic carbon in ground water range from modern (post-1950) to more than 30,000 years before present, and appear to be particularly well defined in the predominantly siliciclastic aquifer system. Major sources of water to the basin include (1) recharge from mountains along the north, east and southwest margins (median age 5,000-9,000 years); (2) seepage from the Rio Grande and Rio Puerco (median age 4,000-8,000 years), and from Abo and Tijeras Arroyos (median age 3,000-9,000 years); (3) inflow of saline water along the southwestern basin margin (median age 20,000 years); and (4) inflow along the northern basin margin that probably represents recharge from the Jemez Mountains during the last glacial period (median age 20,000 years). Water recharged from the Jemez Mountains during the last glacial period occurs at the water table in the central part of the basin and beneath younger recharge along the Rio Grande and the northern mountain front. In some parts of the basin, boundaries between hydrochemical zones appear to be near major faults that may affect ground-water flow. However, in other parts of the basin, such as along the east side of Albuquerque near the Sandia Fault zone, ground-water flow appears to be unaffected by major faults. Upward leakage of saline water occurs along some faults and can be a source of salinity and elevated arsenic concentrations in some ground water. A trough in the modern and predevelopment water table west of Albuquerque is centered along a zone of predominantly late Pleistocene age water through the center of the basin and is flanked and overlain along the trough boundary by water that infiltrated from the Rio Puerco on the west and the Rio Grande to the east. It is suggested that the groundwater trough is a relatively recent transient feature of the Santa Fe Group aquifer system. At Albuquerque, a distinct north-south boundary in deuterium content of ground water marks the division between recharge from the eastern mountain front and that from the Rio Grande.

New Mexico↗
USGS1997

Vulnerability of ground water to atrazine leaching in Kent County, Michigan

A steady-state model of pesticide leaching through the unsaturated zone was used with readily available hydrologic, lithologic, and pesticide characteristics to estimate the vulnerability of the near-surface aquifer to atrazine contamination from non-point sources in Kent County, Michigan. The modelcomputed fraction of atrazine remaining at the water table, RM , was used as the vulnerability criterion; time of travel to the water table also was computed. Model results indicate that the average fraction of atrazine remaining at the water table was 0.039 percent; the fraction ranged from 0 to 3.6 percent. Time of travel of atrazine from the soil surface to the water table averaged 17.7 years and ranged from 2.2 to 118 years. Three maps were generated to present three views of the same atrazine vulnerability characteristics using different metrics (nonlinear transformations of the computed fractions remaining). The metrics were chosen because of the highly (right) skewed distribution of computed fractions. The first metric, rm = RM λ (where λ was 0.0625), depicts a relatively uniform distribution of vulnerability across the county with localized areas of high and low vulnerability visible. The second metric, rm λ-0.5 , depicts about one-half the county at low vulnerability with discontinuous patterns of high vulnerability evident. In the third metric, rm λ-1.0 ( RM ), more than 95 percent of the county appears to have low vulnerability; small, distinct areas of high vulnerability are present. Aquifer vulnerability estimates in the RM metric were used with a steady-state, uniform atrazine application rate to compute a potential concentration of atrazine in leachate reaching the water table. The average estimated potential atrazine concentration in leachate at the water table was 0.16 μg/L (micrograms per liter) in the model area; estimated potential concentrations ranged from 0 to 26 μg/L. About 2 percent of the model area had estimated potential atrazine concentrations in leachate at the water table that exceeded the USEPA (U.S. Environmental Protection Agency) maximum contaminant level of 3 μg/L. Uncertainty analyses were used to assess effects of parameter uncertainty and spatial interpolation error on the variability of the estimated fractions of atrazine remaining at the water table. Results of Monte Carlo simulations indicate that parameter uncertainty is associated with a standard error of 0.0875 in the computed fractions (in the rm metric). Results of kriging analysis indicate that errors in spatial interpolation are associated with a standard error of 0.146 (in the rm metric). Thus, uncertainty in fractions remaining is primarily associated with spatial interpolation error, which can be reduced by increasing the density of points where the leaching model is applied. A sensitivity analysis indicated which of 13 hydrologic, lithologic, and pesticide characteristics were influential in determining fractions of atrazine remaining at the water table. Results indicate that fractions remaining are most sensitive to the unit changes in pesticide half life and in organic-carbon content in soils and unweathered rocks, and least sensitive to infiltration rates. The leaching model applied in this report provides an estimate of the vulnerability of the near-surface aquifer in Kent County to contamination by atrazine. The vulnerability estimate is related to water-quality criteria developed by the USEPA to help assess potential risks from atrazine to the near-surface aquifer. However, atrazine accounts for only 28 percent of the herbicide use in the county; additional potential for contamination exists from other pesticides and pesticide metabolites. Therefore, additional work is needed to develop a comprehensive understanding of the relative risks associated with specific pesticides. The modeling approach described in this report provides a technique for estimating relative vulnerabilities to specific pesticides and for helping to assess potential risks.

Michigan↗
USGS2002

Rainfall-runoff characteristics and effects of increased urban density on streamflow and infiltration in the eastern part of the San Jacinto River basin, Riverside County, California

To better understand the rainfall-runoff characteristics of the eastern part of the San Jacinto River Basin and to estimate the effects of increased urbanization on streamflow, channel infiltration, and land-surface infiltration, a long-term (1950?98) time series of monthly flows in and out of the channels and land surfaces were simulated using the Hydrologic Simulation Program- FORTRAN (HSPF) rainfall-runoff model. Channel and land-surface infiltration includes rainfall or runoff that infiltrates past the zone of evapotranspiration and may become ground-water recharge. The study area encompasses about 256 square miles of the San Jacinto River drainage basin in Riverside County, California. Daily streamflow (for periods with available data between 1950 and 1998), and daily rainfall and evaporation (1950?98) data; monthly reservoir storage data (1961?98); and estimated mean annual reservoir inflow data (for 1974 conditions) were used to calibrate the rainfall-runoff model. Measured and simulated mean annual streamflows for the San Jacinto River near San Jacinto streamflow-gaging station (North-South Fork subbasin) for 1950?91 and 1997?98 were 14,000 and 14,200 acre-feet, respectively, a difference of 1.4 percent. The standard error of the mean for measured and simulated annual streamflow in the North-South Fork subbasin was 3,520 and 3,160 acre-feet, respectively. Measured and simulated mean annual streamflows for the Bautista Creek streamflow-gaging station (Bautista Creek subbasin) for 1950?98 were 980 acre-feet and 991 acre-feet, respectively, a difference of 1.1 percent. The standard error of the mean for measured and simulated annual streamflow in the Bautista Creek subbasin was 299 and 217 acre-feet, respectively. Measured and simulated annual streamflows for the San Jacinto River above State Street near San Jacinto streamflow-gaging station (Poppet subbasin) for 1998 were 23,400 and 23,500 acre-feet, respectively, a difference of 0.4 percent. The simulated mean annual streamflow for the State Street gaging station at the outlet of the study basin and the simulated mean annual basin infiltration (combined infiltration from all the channels and land surfaces) were 8,720 and 41,600 acre-feet, respectively, for water years 1950-98. Simulated annual streamflow at the State Street gaging station ranged from 16.8 acre-feet in water year 1961 to 70,400 acre-feet in water year 1993, and simulated basin infiltration ranged from 2,770 acre-feet in water year 1961 to 149,000 acre-feet in water year 1983.The effects of increased urbanization on the hydrology of the study basin were evaluated by increasing the size of the effective impervious and non-effective impervious urban areas simulated in the calibrated rainfall-runoff model by 50 and 100 percent, respectively. The rainfall-runoff model simulated a long-term time series of monthly flows in and out of the channels and land surfaces using daily rainfall and potential evaporation data for water years 1950?98. Increasing the effective impervious and non-effective impervious urban areas by 100 percent resulted in a 5-percent increase in simulated mean annual streamflow at the State Street gaging station, and a 2.2-percent increase in simulated basin infiltration. Results of a frequency analysis of the simulated annual streamflow at the State Street gaging station showed that when effective impervious and non-effective impervious areas were increased 100 percent, simulated annual streamflow increased about 100 percent for low-flow conditions and was unchanged for high-flow conditions. The simulated increase in streamflow at the State Street gaging station potentially could infiltrate along the stream channel further downstream, outside of the model area.

Water-Resources Investigations Report↗
USGS2002

Ground-water flow simulation and chemical and isotopic mixing equation analysis to determine source contributions to the Missouri River alluvial aquifer in the vicinity of the Independence, Missouri, well field

The city of Independence, Missouri, operates a well field in the Missouri River alluvial aquifer. Steady-state ground-water flow simulation, particle tracking, and the use of chemical and isotopic composition of river water, ground water, and well-field pumpage in a two-component mixing equation were used to determine the source contributions of induced inflow from the Missouri River and recharge to ground water from precipitation in well-field pumpage. Steady-state flow-budget analysis for the simulation-defined zone of contribution to the Independence well field indicates that 86.7 percent of well-field pumpage is from induced inflow from the river, and 6.7 percent is from ground-water recharge from precipitation. The 6.6 percent of flow from outside the simulation-defined zone of contribution is a measure of the uncertainty of the estimation, and occurs because model cells are too large to uniquely define the actual zone of contribution. Flow-budget calculations indicate that the largest source of water to most wells is the Missouri River. Particle-tracking techniques indicate that the Missouri River supplies 82.3 percent of the water to the Independence well field, ground-water recharge from precipitation supplies 9.7 percent, and flow from outside defined zones of contribution supplies 8.0 percent. Particle tracking was used to determine the relative amounts of source water to total well-field pumpage as a function of traveltime from the source. Well-field pumpage that traveled 1 year or less from the source was 8.8 percent, with 0.6 percent from the Missouri River, none from precipitation, and 8.2 percent between starting cells. Well-field pumpage that traveled 2 years or less from the source was 10.3 percent, with 1.8 percent from the Missouri River, 0.2 percent from precipitation, and 8.3 percent between starting cells. Well-field pumpage that traveled 5 years or less from the source was 36.5 percent, with 27.1 percent from the Missouri River, 1.1 percent from precipitation, and 8.3 percent between starting cells. Well-field pumpage that traveled 10 years or less from the source was 42.7 percent, with 32.6 percent from the Missouri River, 1.8 percent from precipitation, and 8.3 percent between starting cells. Well-field pumpage that traveled 25 years or less from the source was 71.9 percent, with 58.9 percent from the Missouri River, 4.7 percent from precipitation, and 8.3 percent between starting cells. Results of chemical (calcium, sodium, iron, and fluoride) and isotopic (oxygen and hydrogen) analyses of water samples collected from the Missouri River, selected monitoring wells around the Independence well field, and combined well-field pumpage were used in a two component mixing equation to estimate the relative amount of Missouri River water in total well-field pumpage. The relative amounts of induced inflow from the Missouri River in well-field pumpage ranged from 49 percent for sodium to 80 percent for calcium, and sensitivities ranged from 0 percent for iron to plus or minus 35 percent for naturally occurring stable isotope ( 18 O). The average of all mixing equation results indicated that 61 percent of well-field pumpage was from induced inflow from the Missouri River. All methods used in the study indicate that more than one-half of the water in well-field pumpage was inflow from the Missouri River. River inflow estimates from ground-water simulation methods are larger and error values are smaller than those using chemical and isotopic data in the mixing equation, although substantial uncertainties exist for both estimation methods. Because of the complex hydrology of the aquifer near the Independence well field, the source estimates using particle tracking probably are the most reliable of the ground-water simulation methods. Mixing equation results are less reliable than those of the ground-water simulation for this study. However, more reliable results can be obtained from the mixing equation by increasing the number of samples and collecting samples for a longer period of time, and during different flow conditions. In the absence of a calibrated ground-water flow simulation, the mixing equation can provide a reasonable estimate of the sources of water to a well field at relatively low cost, if sources of error are clearly understood.

Water-Resources Investigations Report↗
USGS2005

Compilation of geologic, hydrologic, and ground-water flow modeling information for the Spokane Valley-Rathdrum Prairie aquifer, Spokane County, Washington, and Bonner and Kootenai Counties, Idaho

The U.S. Geological Survey, in cooperation with the Idaho Department of Water Resources and Washington Department of Ecology compiled and described geologic, hydrologic, and ground-water flow modeling information about the Spokane Valley-Rathdrum Prairie (SVRP) aquifer in northern Idaho and northeastern Washington. Descriptions of the hydrogeologic framework, water-budget components, ground- and surface-water interactions, computer flow models, and further data needs are provided. The SVRP aquifer, which covers about 370 square miles including the Rathdrum Prairie, Idaho and the Spokane valley and Hillyard Trough, Washington, was designated a Sole Source Aquifer by the U.S. Environmental Protection Agency in 1978. Continued growth, water management issues, and potential effects on water availability and water quality in the aquifer and in the Spokane and Little Spokane Rivers have illustrated the need to better understand and manage the region's water resources. The SVRP aquifer is composed of sand, gravel, cobbles, and boulders primarily deposited by a series of catastrophic glacial outburst floods from ancient Glacial Lake Missoula. The material deposited in this high-energy environment is coarser-grained than is typical for most basin-fill deposits, resulting in an unusually productive aquifer with well yields as high as 40,000 gallons per minute. In most places, the aquifer is bounded laterally by bedrock composed of granite, metasedimentary rocks, or basalt. The lower boundary of the aquifer is largely unknown except along the margins or in shallower parts of the aquifer where wells have penetrated its entire thickness and reached bedrock or silt and clay deposits. Based on surface geophysics, the thickness of the aquifer is about 500 ft near the Washington-Idaho state line, but more than 600 feet within the Rathdrum Prairie and more than 700 feet in the Hillyard trough based on drilling records. Depth to water in the aquifer is greatest in the northern Rathdrum Prairie (about 500 feet) and least near the city of Spokane along the Spokane River (less than about 50 feet). Ground-water flow is south from near the southern end of Lake Pend Oreille and Hoodoo Valley, through the Rathdrum Prairie, then west toward Spokane. In Spokane, the aquifer splits and water moves north through the Hillyard Trough as well as west through the Trinity Trough. From the Trinity Trough water flows north along the western arm of the aquifer. The aquifer's discharge area is along the Little Spokane River and near Long Lake, Washington. A compilation of estimates of water-budget components, including recharge (precipitation, irrigation, canal leakage, septic tank effluent, inflow from tributary basins, and flow from the Spokane River) and discharge (withdrawals from wells, flow to the Spokane and Little Spokane Rivers, evapotranspiration, and underflow to Long Lake) illustrates that these estimated values should be compared with caution due to several variables including the area and time period of interest as well as methods employed in making the estimates. Numerous studies have documented the dynamic ground-water and surface-water interaction between the SVRP aquifer and the Spokane and Little Spokane Rivers. Gains and losses vary throughout the year, as well as the locations of gains and losses. September 2004 streamflow measurements indicated that the upper reach of the Spokane River between Post Falls and downstream at Flora Road lost 321 cubic feet per second. A gain of 736 cubic feet per second was measured between the Flora Road site and downstream at Green Street Bridge. A loss of 124 cubic feet per second was measured for the reach between the Green Street Bridge and the Spokane River at Spokane gaging station. The river gained about 87 cubic feet per second between the Spokane River at Spokane gaging station and the TJ Meenach Bridge. Overall, the Spokane River gained about 284 cubic feet per second between the Post Falls,

Scientific Investigations Report↗
USGS2009

Geologic Map of MTM 35337, 40337, and 45337 Quadrangles, Deuteronilus Mensae Region of Mars

Deuteronilus Mensae, first defined as an albedo feature at lat 35.0 deg N., long 5.0 deg E., by U.S. Geological Survey (USGS) and International Astronomical Union (IAU) nomenclature, is a gradational zone along the dichotomy boundary in the northern mid-latitudes of Mars. The boundary in this location includes the transition from the rugged cratered highlands of Arabia Terra to the northern lowland plains of Acidalia Planitia. Within Deuteronilus Mensae, polygonal mesas are prominent along with features diagnostic of Martian fretted terrain, including lobate debris aprons, lineated valley fill, and concentric crater fill. Lobate debris aprons, as well as the valley and crater fill deposits, are geomorphic indicators of ground ice, and their concentration in Deuteronilus Mensae is of great interest because of their potential association with Martian climate change. The paucity of impact craters on the surfaces of debris aprons and the presence of ice-cemented mantle material imply young (for example, Amazonian) surface ages that are consistent with recent climate change in this region of Mars. North of Deuteronilus Mensae are the northern lowlands, a potential depositional sink that may have had large standing bodies of water or an ocean in the past. The northern lowlands have elevations that are several kilometers below the ancient cratered highlands with significantly younger surface ages. The morphologic and topographic characteristics of the Deuteronilus Mensae region record a diverse geologic history, including significant modification of the ancient highland plateau and resurfacing of low-lying regions. Previous studies of this region have interpreted a complex array of geologic processes, including eolian, fluvial and glacial activity, coastal erosion, marine deposition, mass wasting, tectonic faulting, effusive volcanism, and hydrovolcanism. The origin and age of the Martian crustal dichotomy boundary are fundamental questions that remain unresolved at the present time. Several scenarios for its formation, including single and multiple large impact events, have been proposed and debated in the literature. Endogenic processes whereby crust is thinned by internal mantle convection and tectonic processes have also been proposed. Planetary accretion models and isotopic data from Martian meteorites suggest that the crust formed very early in Martian history. Using populations of quasi-circular depressions extracted from the topography of Mars, other studies suggest that the age difference between the highlands and lowlands could be ~100 m.y.. Furthermore, understanding the origin and age of the dichotomy boundary has been made more complicated due to significant erosion and deposition that have modified the boundary and its adjacent regions. The resulting diversity of terrains and features is likely a combined result of ancient and recent events. Detailed geologic analyses of dichotomy boundary zones are important for understanding the spatial and temporal variations in highland evolution. This information, and comparisons to other highland regions, can help elucidate the scale of potential environmental changes. Previous geomorphic and geologic mapping investigations of the Deuteronilus Mensae region have been completed at local to global scales. The regional geology was first mapped by Lucchitta (1978) at 1:5,000,000 scale using Mariner 9 data. This study concluded that high crater flux early in Martian history formed overlapping craters and basins that were later filled by voluminous lava flows that buried the impacted surface, creating the highlands. After this period of heavy bombardment, fluvial erosion of the highlands formed the canyons and valleys, followed by dissection that created the small mesas and buttes, and later, formation of the steep escarpment marking the present-day northern highland margin. After valley dissection, mass wasting and eolian processes caused lateral retreat of mesas and buttes

Scientific Investigations Map↗