Search USGS⌕ Search

SEARCH · Search USGS

Results for “Ground Water Modeling Investigation Report”

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

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

At least 505 records · Page 28Linked to original sources

The computer model Hydrotherm, a three-dimensional finite-difference model to simulate ground-water flow and heat transport in the temperature range of 0 to 1,200 degrees C

Quantitative modeling of the deep parts of magmatic- hydrothermal systems has been limited by the lack of publicly available, documented computer models for multiphase, high-temperature flow. This report documents HYDROTHERM, a finite-difference model for three-dimensional, multiphase flow of pure water and heat over a temperature range of 0 to 1,200 degrees Celsius and a pressure range of 0.5 *0 10,000 bars.HYDROTHEP*M is an extension of multiphase geothermal models developed by the U.S. Geological Survey in the 1970's. It solves numerical approximations to mass- and energy-balance equations that are posed in terms of pressure and enthalpy. Major assumptions are that the rock matrix can be treated as a porous medium; that water and rock are in thermal equilibrium; and that capillary pressure is negligible.

Water-Resources Investigations Report↗

Approach for delineation of contributing areas and zones of transport to selected public-supply wells using a regional ground-water flow model, Palm Beach County, Florida

Rapid urban development and population growth in Palm Beach County, Florida, have been accompanied with the need for additional freshwater withdrawals from the surficial aquifer system. To maintain water quality, County officials protect capture areas and determine zones of transport of municipal supply wells. A multistep process was used to help automate the delineation of wellhead protection areas. A modular ground-water flow model (MODFLOW) Telescopic Mesh Refinement program (MODTMR) was used to construct an embedded flow model and combined with particle tracking to delineate zones of transport to supply wells; model output was coupled with a geographic information system. An embedded flow MODFLOW model was constructed using input and output file data from a preexisting three-dimensional, calibrated model of the surficial aquifer system. Three graphical user interfaces for use with the geographic information software, ArcView, were developed to enhance the telescopic mesh refinement process. These interfaces include AvMODTMR for use with MODTMR; AvHDRD to build MODFLOW river and drain input files from dynamically segmented linear (canals) data sets; and AvWELL Refiner, an interface designed to examine and convert well coverage spatial data layers to a MODFLOW Well package input file. MODPATH (the U.S. Geological Survey particle-tracking postprocessing program) and MODTOOLS (the set of U.S. Geological Survey computer programs to translate MODFLOW and MODPATH output to a geographic information system) were used to map zones of transport. A steady-state, five-layer model of the Boca Raton area was created using the telescopic mesh refinement process and calibrated to average conditions during January 1989 to June 1990. A sensitivity analysis of various model parameters indicates that the model is most sensitive to changes in recharge rates, hydraulic conductivity for layer 1, and leakance for layers 3 and 4 (Biscayne aquifer). Recharge (58 percent); river (canal) leakance (29 percent); and inflow through the northern, western, and southern prescribed flux model boundaries (10 percent) represent the major inflow components. Principal outflow components in the Boca Raton well field area include well discharge (56 percent), river (canal) leakance (27 percent), and water that discharges along the coast (10 percent). A particle-tracking analysis using MODPATH was conducted to better understand well-field ground-water flow patterns and time of travel. MODTOOLS was used to construct zones-of-transport spatial data for municipal supply wells. Porosity estimates were uniformly increased to study the effect of porosity on zones of transport. Where porosity was increased, the size of the zones of transport were shown to decrease.

Water-Resources Investigations Report↗

Documentation of revisions to the regional aquifer system analysis model of the New Jersey coastal plain

The model, which simulates flow in the New Jersey Coastal Plain sediments, developed for the U.S. Geological Survey Regional Aquifer System Analysis (RASA) program was revised. The RASA model was revised with (1) a rediscretization of the model parameters with a finer cell size, (2) a spatially variable recharge rate that is based on rates determined by recent studies and, (3) ground-water withdrawal data from 1981 to 1998. The RASA model framework, which subdivided the Coastal Plain sediments into 10 aquifers and 9 confining units, was preserved in the revised model. A transient model that simulates flow conditions from January 1, 1968 to December 31, 1998, was constructed using 21 stress periods. The model was calibrated by attempting to match the simulated results with (1) estimated base flow for five river basins, (2) measured water levels in long-term hydrographs for 28 selected observation wells, and (3) potentiometric surfaces in the model area for 1978, 1983, 1998, 1993, and 1998 conditions. The estimated and simulated base flow in the five river basins compare well. In general, the simulated water levels matched the interpreted potentiometric surfaces and the measured water levels of the hydrographs within 25 feet.

Water-Resources Investigations Report↗

Assessment of intrinsic bioremediation of jet fuel contamination in a shallow aquifer, Beaufort, South Carolina

Field and laboratory studies show that microorganisms indigenous to the ground-water system underlying Tank Farm C, Marine Corps Air Station Beaufort, S.C., degrade petroleum hydrocarbons under aerobic and anaerobic conditions. Under aerobic conditions, sediments from the shallow aquifer underlying the site mineralized radiolabeled (14C) toluene to 14CO2 with first-order rate constants of about -0.29 per day. Sediments incubated under anaerobic conditions mineralized radiolabeled toluene more slowly, with first-order rate constants of -0.001 per day. Although anaerobic rates of biodegradation are low, they are significant in the hydrologic and geochemical context of the site. Because of low hydraulic conductivities (1.9-9.1 feet per day) and low hydraulic gradients (about 0.004 feet per feet), ground water flows slowly (approximately 20 feet per year) at this site. Furthermore, aquifer sediments contain organic-rich peat that has a high sorptive capacity. Under these conditions, hydrocarbon contaminants have moved no further than 10 feet downgradient of the jet fuel free product. Digital solute-transport simulations, using the range of model parameters measured at the site, show that dissolved contaminants will be completely degraded before they are discharged from the aquifer into adjacent surface-water bodies. These results show that natural attenuation processes are containing the migration of soluble hydrocarbons, and that intrinsic bioremediation is a potentially effective remedial strategy at this site.

South Carolina↗

Preliminary estimates of spatially distributed net infiltration and recharge for the Death Valley region, Nevada-California

A three-dimensional ground-water flow model has been developed to evaluate the Death Valley regional flow system, which includes ground water beneath the Nevada Test Site. Estimates of spatially distributed net infiltration and recharge are needed to define upper boundary conditions. This study presents a preliminary application of a conceptual and numerical model of net infiltration. The model was developed in studies at Yucca Mountain, Nevada, which is located in the approximate center of the Death Valley ground-water flow system. The conceptual model describes the effects of precipitation, runoff, evapotranspiration, and redistribution of water in the shallow unsaturated zone on predicted rates of net infiltration; precipitation and soil depth are the two most significant variables. The conceptual model was tested using a preliminary numerical model based on energy- and water-balance calculations. Daily precipitation for 1980 through 1995, averaging 202 millimeters per year over the 39,556 square kilometers area of the ground-water flow model, was input to the numerical model to simulate net infiltration ranging from zero for a soil thickness greater than 6 meters to over 350 millimeters per year for thin soils at high elevations in the Spring Mountains overlying permeable bedrock. Estimated average net infiltration over the entire ground-water flow model domain is 7.8 millimeters per year.To evaluate the application of the net-infiltration model developed on a local scale at Yucca Mountain, to net-infiltration estimates representing the magnitude and distribution of recharge on a regional scale, the net-infiltration results were compared with recharge estimates obtained using empirical methods. Comparison of model results with previous estimates of basinwide recharge suggests that the net-infiltration estimates obtained using this model may overestimate recharge because of uncertainty in modeled precipitation, bedrock permeability, and soil properties for locations such as the Spring Mountains. Although this model is preliminary and uncalibrated, it provides a first approximation of the spatial distribution of net infiltration for the Death Valley region under current climatic conditions.

Water-Resources Investigations Report↗

Water resources of the Batavia Kill basin at Windham, Greene County, New York

The water resources of a 27.6-square-mile section of the Batavia Kill Basin near the village of Windham, N.Y., which has undergone substantial development, were evaluated. The evaluation entailed (1) estimation of the magnitude and distribution of several hydrologic components, including recharge, (2) measurement of discharge and chemical quality of the Batavia Kill and selected tributaries, (3) analysis of ground-water flow and chemistry, and (4) a conceptualization of the ground-water flow system. The region consists of deeply dissected, relatively flat-lying, clastic sedimentary sequences variably overlain by as much as 120 feet of glacial deposits. The types of bedrock fractures and their distribution in the Batavia Kill valley are consistent with valley stress-relief characteristics. Till predominates in the uplands, and stratified drift typically dominates within the valley of the Batavia Kill and the lower section of its largest tributary valley (Mitchell Hollow). Fractured bedrock is the most commonly used water source within the study area. The areas of highest yielding bedrock generally are with valleys, where the shallow fractures are saturated. Stratified-drift aquifers are also limited to the largest valleys; the greatest saturated thicknesses are in the Batavia Kill valley at Windham. A conceptual model of ground-water flow within the study areas suggests that the zones of most active flow are shallow fractured bedrock in upland areas and the shallow stratified drift in the largest valleys. The hydrogeologic system has been altered by development; major effects include (1) chemical alteration of natural ground-water and surface-water quality by point- and nonpoint-source contaminants, (2) hydraulic interconnection of other-wise isolated bedrock fractures by wellbores, and (3) drawdowns in wells within the Batavia Kill valley by pumping from the bedrock aquifer. Water resource development of the most promising unconsolidated aquifer beneath Windham may be precluded by the potential for contamination by leachate from an abandoned landfill, road-salt stockpiles, and domestic septic systems in the area.

New York↗

Hydrogeology, hydrologic budget, and water chemistry of the Medina Lake area, Texas

A three-phase study of the Medina Lake area in Texas was done to assess the hydrogeology and hydrology of Medina and Diversion Lakes combined (the lake system) and to determine what fraction of seepage losses from the lake system might enter the regional ground-water-flow system of the Edwards and (or) Trinity aquifers. Phase 1 consisted of revising the geologic framework for the Medina Lake area. Results of field mapping show that the upper member of the Glen Rose Limestone underlies Medina Lake and the intervening stream channel from the outflow of Medina Lake to the midpoint of Diversion Lake, where the Diversion Lake fault intersects Diversion Lake. A thin sequence of strata consisting primarily of the basal nodular and dolomitic members of the Kainer Formation of the Edwards Group, is present in the southern part of the study area. On the southern side of Medina Lake, the contact between the upper member of the Glen Rose Limestone and the basal nodular member is approximately 1,000 feet above mean sea level, and the contact between the basal nodular member and the dolomitic member is approximately 1,050 feet above mean sea level. The most porous and permeable part of the basal nodular member is about 1,045 feet above mean sea level. At these altitudes, Medina Lake is in hydrologic connection with rocks in the Edwards aquifer recharge zone, and Medina Lake appears to lose more water to the ground-water system along this bedding plane contact. Hydrologic budgets calculated during phase 2 for Medina Lake, Diversion Lake, and Medina/Diversion Lakes combined indicate that: (1) losses from Medina and Diversion Lakes can be quantified; (2) a portion of those losses are entering the Edwards aquifer; and (3) losses to the Trinity aquifer in the Medina Lake area are minimal and within the error of the hydrologic budgets. Hydrologic budgets based on streamflow, precipitation, evaporation, and change in lake storage were used to quantify losses (recharge) to the ground-water system from Medina Lake, Diversion Lake, and Medina/Diversion Lakes combined during October 1995–September 1996. Water losses from Medina Lake to the Edwards/Trinity aquifers ranged from -14.0 to 135 acre-feet per day; Diversion Lake ranged from -1.2 to 93.1 acre-feet per day; and Medina/Diversion Lakes combined ranged from 36.1 to 119 acre-feet per day. Monthly average recharge during December 1995–July 1996 was estimated using an alternative method developed during this study (current study method) and compared to monthly average recharge during December 1995–July 1996 estimated using the existing USGS method and the Trans-Texas method. Recharge to the Edwards aquifer estimated using the current study method was about 69 and 73 percent of the recharge estimated using the USGS and Trans-Texas methods, respectively. The USGS and Trans-Texas methods overestimated recharge from Medina Lake compared to the recharge estimated with the current study method when Medina Lake stage was between about 1,027 and 1,032 feet above mean sea level and underestimated recharge from Medina Lake when lake stage was between about 1,036 and 1,045 feet above mean sea level. The USGS and Trans-Texas methods underestimated recharge from Diversion Lake compared to the recharge estimated with the current study method when Diversion Lake stage was greater than 913 feet above mean sea level and overestimated recharge from Diversion Lake when lake stage was less than 913 feet above mean sea level. The water quality of Medina Lake and Medina River and in selected wells and springs in the Edwards and Trinity aquifers was characterized during phase 3 of the study. Environmental isotope analyses and geochemical modeling also were used to determine where water losses from the lake system might be entering the ground-water-flow system. Isotopic ratios of deuterium, oxygen, and strontium were analyzed in selected surface-water, lake-water, and ground-water samples to trace the isotopic “signature” of the lake water as it mixes with the ground water and to determine the fraction of lake water and ground water in selected Edwards aquifer wells. Isotopic data and geochemical modeling were used to show that lake water is moving into the Edwards aquifer in two fault blocks in the eastern Medina storage unit. One fault block is bounded on the north by the Vandenburg School fault and on the south by the Haby Crossing fault, and the second fault block is bounded on the north by the Diversion Lake fault and on the south by the Haby Crossing fault. In selected Edwards aquifer wells located southwest of Medina Lake and west of Diversion Lake, the proportion of lake water ranged from about 10 to 45 percent. Geochemical modeling using NETPATH confirms the degree of mixing between lake water and aquifer water shown by the isotopes.

Texas↗

A two-constituent solute-transport model for ground water having variable density

A numerical model has been developed to simulate solute transport and dispersion of either one or two constituents in groundwater where there is two-dimensional, density-dependent flow. The model is a modified version of the one documented by Konikow and Bredehoeft (1978), which uses finite-difference methods and the method of characteristics to solve the flow and transport equations. The model was tested on an idealized seawater intrusion problem for which an analytical solution has been developed. The results were nearly identical to those of other numerical models tested on the same problem. A description of the formats for the input data, a sample of input and output for a two-constituent example problem, and a listing of the Fortran program are presented. (Author 's abstract)

Water-Resources Investigations Report↗

Evaluation of availability of water from drift aquifers near the Pomme de Terre and Chippewa rivers, western Minnesota

Ground-water flow in the confined- and unconfined-drift aquifers near Appleton and Benson, Minnesota, was simulated with a three-dimensional finite-difference ground-water-flow model. Model results indicate that 98 percent of the total inflow to the modeled area is from precipitation. Of the total outflow, 38 percent is ground-water discharge to the Pom me de Terre and Chippewa Rivers, 36 percent is evapotranspiration, 17 percent is ground-water pumpage, and 8 percent is ground-water discharge to the Minnesota River. The model was used to simulate the effects of below-normal precipitation (drought) and hypothetical increases in ground-water development. Model results indicate that reduced recharge and increased pumping during a three-year extended drought probably would lower water levels 2 to 6 feet regionally in the surficial aquifer and in the Appleton and Benson-middle aquifers and as much as 11 feet near aquifer boundaries. Ground-water discharge to the Pomme de Terre and Chippewa Rivers in the modeled area probably would be reduced during the simulated drought by 15.2 and 7.4 cubic feet per second, respectively, compared to 1982 conditions. The addition of 30 hypothetical wells in the Benson-middle aquifer near Benson, pumping a total of 810 million gallons per year, resulted in water-level declines of as much as 1.3 and 2.7 feet in the surficial and Benson-middle aquifers, respectively. The addition of 28 hypothetical wells in the Appleton aquifer east and southeast of Appleton, pumping a total of 756 million gallons per year, lowered water levels in the surficial and Appleton confined aquifers as much as 5 feet.

Minnesota↗

Hydrogeology and simulation of ground-water flow at U.S. Marine Corps Air Station, Cherry Point, North Carolina, 1987-90

Geophysical and lithologic well-log data from 30 wells and chloride data, and water-level data from oil-test wells, supply wells, and observation wells were evaluated to define the hydrogeologic framework at the U.S. Marine Corps Air Station, Cherry Point, North Carolina. Elements of the hydrogeologic framework important to this study include six aquifers and their respective confining units. In descending order, these aquifers are the surficial, Yorktown, Pungo River, upper and lower Castle Hayne, and Beaufort. The upper and lower Castle Hayne and Beaufort aquifers and related confining units are relatively continuous throughout the study area. The surficial, Yorktown, Pungo River, and upper and lower Castle Hayne aquifers contain freshwater. The upper and lower Castle Hayne aquifers serve as the Air Station?s principal supply of freshwater. However, the lower Castle Hayne aquifer contains brackish water near its base and there is potential for upward movement of this water to supply wells completed in this aquifer. The potential for brackish-water encroachment is greatest if wells are screened too deep in the lower Castle Hayne aquifer or if pumping rates are too high. Lateral movement of brackish water into aquifers incised by estuarine streams is also possible if ground-water flow gradients toward these bodies are reversed by pumping. The potential for the reversed movement of water from the surficial aquifer downward to the water-supply aquifer is greatest in areas where clay confining units are missing. These missing clay units could indicate the presence of a paleochannel of the Neuse River. A quasi three-dimensional finite-difference ground-water flow model was constructed and calibrated to simulate conditions at and in the vicinity of the Air Station for the period of 1987-90. Comparisons of 94 observed and computed heads were made, and the average difference between them is -0.2 feet with a root mean square error of 5.7 feet. An analysis was made to evaluate the sensitivity of the model to the absence of the Yorktown and Pungo River confining units in a 1-square-mile area in the southern part of the Air Station. This analysis resulted in a maximum simulated head increase of 2 feet in one 0.11-square-mile model cell in the Pungo River aquifer.

North Carolina↗

Plan of study for the regional aquifer-system analysis of the Snake River plain, Idaho and eastern Oregon

The 15,600-square-mile Snake River Plain is largely in southern Idaho and includes one of the Nation 's major regional aquifers. A comprehensive investigation of the area 's ground-water resources will be made as part of the U.S. Geological Survey 's Regional Aquifer-System Analysis (RASA) program. Basaltic and sedimentary rocks in the Snake River Plain yield large quantities of water that are vital to the area 's agricultural economy. Basaltic rocks predominate in the eastern Snake River Plain and have especially high water-yielding capabilities. Surface water, largely from the Snake River, is extensively used for irrigation and is a major source of recharge to the ground-water system. Springs issuing from basaltic rocks that form the Snake River Canyon wall near Twin Falls are the major points of ground-water discharge. Increased use of ground water for irrigation is causing concern as to the effect of large-scale withdrawals on spring flow. Ground-water flow models will be used to improve understanding of the hydrologic system, and, if feasible, to aid in evaluating management alternatives. Ground-water quality will be defined and geochemical techniques used to determine the effects of water-rock reactions on water quality. Several reports are planned on different phases of the project, concluding with a summary report. (USGS)

Open-File Report↗

Simulated effects of ground-water withdrawals and artificial recharge on discharge to streams, springs, and riparian vegetation in the Sierra Vista Subwatershed of the Upper San Pedro Basin, southeastern Arizona

In the context of ground-water resources, “capture” or “streamflow depletion” refers to withdrawal-induced changes in inflow to or outflow from an aquifer. These concepts are helpful in understanding the effects of long-term development of ground-water resources. For the Upper San Pedro Basin in Arizona, USA and Sonora, Mexico, a recently developed ground-water flow model is available to help quantify capture of water from the river and riparian system. A common method of analysis is to compute curves of capture and aquifer-storage change for a range of time at select points of interest. This study, however, presents results of a method to show spatial distributions of total change in inflow and outflow from withdrawal or injection for select times of interest. The mapped areal distributions show the effect of a single well in terms of the ratio of the change in boundary flow rate to rate of withdrawal or injection by the well. To the extent that the system responds linearly to ground-water withdrawal or injection, fractional responses in the mapped distributions can be used to quantify response for any withdrawal or injection rate. Capture distributions calculated using the Upper San Pedro model include response to (1) withdrawal in the lower basin-fill aquifer for times of 10 and 50 years following the initiation of pumping from predevelopment conditions and (2) artificial recharge to the water table in the area underlain by the lower basin-fill aquifer after 10 and 50 years. The mapped distributions show that response to withdrawals and injections is greatest near the river/riparian system. Presence of clay layers in the vertical interval between withdrawal locations and the river/riparian system, however, can delay the response.

Arizona;Sonora↗

Development of a local meteoric water line for southeastern Idaho, western Wyoming, and south-central Montana

Linear-regression analysis was applied to stable hydrogen (H) and oxygen (O) isotope data in 72 snow-core and precipitation samples collected during 1999-2001 to determine the Local Meteoric Water Line (LMWL) for southeastern Idaho, western Wyoming, and south-central Montana. On the basis of (1) residuals from the regression model, (2) comparison of study-area deuterium-excess (d-excess) values with the global range of d-excess values, and (3) outlier analysis by means of Chauvenet's Criterion, values of four samples were excluded from final regression analysis of the dataset. Regression results for the 68 remaining samples yielded a LMWL defined by the equation ?H = 7.95 18O + 8.09 (r? = 0.98). This equation will be useful as a reference point for future studies in this area that use stable isotopes of H and O to determine sources of ground-water recharge, to determine water-mineral exchange, to evaluate surface-water and groundwater interaction, and to analyze many other geochemical and hydrologic problems.

Scientific Investigations Report↗

Hydrologic effects of storing liquified sewage sludge on strip-mined land, Fulton County, Illinois

The water table near four sewage storage basins in a strip-mined area of western Illinois, has risen about 10 feet since the basins were constructed in 1971. Two dimensional modeling of ground-water flow in the mine spoil indicates that the rise is caused by leakage from one storage basin. The hydrologic-parameter values producing the best fit between computed and observed head values are 0.000007 feet per second for the hydraulic conductivity of the mine spoil, 0.0000000004 feet per second for recharge from the leaking basin. The model indicates the volume of water leaking from the basin is 91,600 cubic yards per year. The principal components of the sewage sludge after the solids have been removed are alkalinity, nitrogen, phosphorous, and chloride. Components in higher concentrations near the basins were sodium, alkalinity, and chloride. Because the sodium and chloride concentrations in the sludge were too low to cause the higher concentrations in the ground water, the strip-mine spoil used in constructing the basins was considered to be the major source of these constituents. This spoil has been moved from its original location and unweathered surfaces exposed, which allowed dissolution of carbonate and chloride and release of sodium through cation exchange.

Illinois↗

Analysis of ground-water flow along a regional flow path of the Midwestern Basins and Arches aquifer system in Ohio

A cross-sectional analysis of ground-water flow in central-western and northwestern Ohio was done as part of the Midwestern Basins and Arches Regional Aquifer-System Analysis project. The Midwestern Basins and Arches aquifer system is composed of carbonate bedrock of Silurian and Devonian age and overlying glacial flow analysis of the Scioto and Blanchard rivers in the study area were used to describe patrems of ground-water flow, to evaluate stream-aquifer interaction, and to quantify recharge and discharge within the ground-water-flow system along a regional ground-water-flow path. The selected regional flow path begins at a regional topographic high in Logan County, Ohio, and ends in Sandusky Bay (Lake Erie), a regional topographic low. Recharge to the ground-water system along the selected regional flow path was estimated from hydrograph separation of streamflow and averaged 3.24 inches per year. Computer model simulations indicate that 84 percent of the water entering the ground-water system flows less than 5 miles from point of recharge to point of discharge and no deeper than the upper surficial aquifers. The distance and depth that ground water travels and traveltime from point of recharge to point of discharge is controlled largely by where ground water enters the flow system. Ground water entering the flow system in the vicinity of major surface- water divides generally travels further, deeper, and longer than ground water entering the flow system elsewhere along the regional flow path. Particle tracking simulations substantiate the concept that the 80-mile-long regional flow path is within a continuous ground-water basin. Estimated traveltimes for ground-water from the regional high to Sandusky Bay range from 22,000 to 40,700 years, given a range of porosities from 8 to 22 percent for the carbonate-rock aquifer.

Ohio↗

Estimating ground-water inflow to lakes in central Florida using the isotope mass-balance approach

The isotope mass-balance approach was used to estimate ground-water inflow to 81 lakes in the central highlands and coastal lowlands of central Florida. The study area is characterized by a subtropical climate and numerous lakes in a mantled karst terrain. Ground-water inflow was computed using both steady-state and transient formulations of the isotope mass-balance equation. More detailed data were collected from two study lakes, including climatic, hydrologic, and isotopic (hydrogen and oxygen isotope ratio) data. For one of these lakes (Lake Starr), ground-water inflow was independently computed from a water-budget study. Climatic and isotopic data collected from the two lakes were similar even though they were in different physiographic settings about 60 miles apart. Isotopic data from all of the study lakes plotted on an evaporation trend line, which had a very similar slope to the theoretical slope computed for Lake Starr. These similarities suggest that data collected from the detailed study lakes can be extrapolated to the rest of the study area. Ground-water inflow computed using the isotope mass-balance approach ranged from 0 to more than 260 inches per year (or 0 to more than 80 percent of total inflows). Steady-state and transient estimates of ground-water inflow were very similar. Computed ground-water inflow was most sensitive to uncertainty in variables used to calculate the isotopic composition of lake evaporate (isotopic compositions of lake water and atmospheric moisture and climatic variables). Transient results were particularly sensitive to changes in the isotopic composition of lake water. Uncertainty in ground-water inflow results is considerably less for lakes with higher ground-water inflow than for lakes with lower ground-water inflow. Because of these uncertainties, the isotope mass-balance approach is better used to distinguish whether ground-water inflow quantities fall within certain ranges of values, rather than for precise quantification. The lakes fit into three categories based on their range of ground-water inflow: low (less than 25 percent of total inflows), medium (25-50 percent of inflows), and high (greater than 50 percent of inflows). The majority of lakes in the coastal lowlands had low ground-water inflow, whereas the majority of lakes in the central highlands had medium to high ground-water inflow. Multiple linear regression models were used to predict ground-water inflow to lakes. These models help identify basin characteristics that are important in controlling ground-water inflow to Florida lakes. Significant explanatory variables include: ratio of basin area to lake surface area, depth to the Upper Floridan aquifer, maximum lake depth, and fraction of wetlands in the basin. Models were improved when lake water-quality data (nitrate, sodium, and iron concentrations) were included, illustrating the link between ground-water geochemistry and lake chemistry. Regression models that considered lakes within specific geographic areas were generally poorer than models for the entire study area. Regression results illustrate how more simplified models based on basin and lake characteristics can be used to estimate ground-water inflow. Although the uncertainty in the amount of ground-water inflow to individual lakes is high, the isotope mass-balance approach was useful in comparing the range of ground-water inflow for numerous Florida lakes. Results were also helpful in understanding differences in the geographic distribution of ground-water inflow between the coastal lowlands and central highlands. In order to use the isotope mass-balance approach to estimate inflow for multiple lakes, it is essential that all the lakes are sampled during the same time period and that detailed isotopic, hydrologic, and climatic data are collected over this same period of time. Isotopic data for Florida lakes can change over time, both seasonally and interannually, primarily because of differ

Water-Resources Investigations Report↗

Hydrogeology and simulation of groundwater flow in the Lucerne Valley groundwater basin, California

The Lucerne Valley is in the southwestern part of the Mojave Desert and is about 75 miles northeast of Los Angeles, California. The Lucerne Valley groundwater basin encompasses about 230 square miles and is separated from the Upper Mojave Valley groundwater basin by splays of the Helendale Fault. Since its settlement, groundwater has been the primary source of water for agricultural, industrial, municipal, and domestic uses. Groundwater withdrawal from pumping has exceeded the amount of water recharged to the basin, causing groundwater declines of more than 100 feet between 1917 and 2016 in the center of the basin. The continued withdrawal has resulted in an increase in pumping costs, reduced well efficiency, and land subsidence near Lucerne Lake. Although the volume of pumping has declined in recent years, there is concern that new agricultural growth and limits on imported water will continue to strain the sustainability of the groundwater system. To address these concerns, the U.S. Geological Survey entered into a cooperative agreement with the Mojave Water Agency to develop a better understanding of the Lucerne Valley hydrogeologic system and provide tools to help evaluate and manage the effects of future development in the Lucerne Valley. The objectives of this study were to (1) improve the understanding of the aquifer system, (2) improve the understanding of subsidence in the basin, and (3) incorporate the understanding into a groundwater-flow model that can be used to help manage the groundwater resources in the Lucerne Valley. The model developed for this study covers the period of 1942–2016 and can help evaluate various proposed water-management scenarios during different climatic and hydrologic conditions. The aquifer system consists of a shallow aquifer, a confining unit, and middle and lower aquifers. These layered water-bearing units were identified based on geologic units of the mostly unconsolidated sediments and hydrologic properties. These alluvial deposits consist of clay, silt, sand, and gravel; some places also contain clay and silty clay lacustrine deposits. Several faults act, at least in part, as barriers to groundwater flow on the eastern, southern, and western edges of the basin. Present-day natural recharge is primarily from the infiltration of runoff from the San Bernardino Mountains to the south; however, stable and radioactive isotopes show that groundwater from the middle of the Lucerne Valley was older than about 10,000 years and probably was recharged as infiltration from streams draining the mountains in the Mojave Desert to the north, which probably does not occur under present-day climatic conditions. The annual average natural recharge for 1942–2016, estimated by a Basin Characterization Model, was about 635 acre-feet per year; the average amount of treated wastewater effluent transferred to the Lucerne Valley for artificial recharge annually ranged from about 1,500 to 4,000 acre-feet per year during 1980–2016. Pumpage estimates for 1942–2016 ranged from about 3,000 acre-feet in 1942 to about 18,300 acre-feet in 1984. The total cumulative amount of groundwater removed from the basin by pumping between 1942 and 2016 was estimated to be about 700,000 acre-feet, which was about 10 times greater than the cumulative amount of recharge to the entire Lucerne Valley groundwater basin. Before groundwater development, the direction of groundwater flow was from the southern part of the basin northward to discharge areas near Lucerne Lake, where it discharged through springs along the Helendale Fault and by evapotranspiration. Since the early 1900s, groundwater-level declines have mostly eliminated the areas where natural discharge occurred and exceeded 100 feet in the middle of the basin between the early 1950s and mid-1990s, and as much as 25 feet near the margins from about the mid-1950s to 2000s. A decrease in the rate of pumping after the mid-1990s lessened the hydraulic stress on the middle and lower aquifers and enabled hydraulic heads in the middle of the basin to recover slightly as groundwater near the margins of the basin moved toward the pumping depression. Although trends in groundwater levels in the center of the basin have reversed since the mid-1990s, levels at the basin margins continue to decline as the movement of groundwater from the margins fills the pumping depression and gradually flattens the groundwater table throughout the basin. The long-term extraction of groundwater and associated dewatering of the fine-grained sediments present within the aquifer system has resulted in aquifer compaction and consequently land subsidence, primarily near Lucerne Lake. Analysis of interferometric synthetic aperture radar data shows that almost 11 inches of land subsidence has occurred south of Lucerne Lake between April 1992 and November 2009; less subsidence occurred elsewhere in the basin during this period. This differential land subsidence has caused fissures and cracks in the ground surface, which have buckled the pavement and undercut roads in several locations. The Lucerne Valley Hydrologic Model was developed using the finite-difference groundwater modeling software One Water Hydrologic Model to represent the hydrologic conditions and stresses during 1942–2016. The model has a uniform grid of approximately 92 acres per cell (2,000 feet by 2,000 feet) and has four layers representing the water-bearing units. The results from the calibrated model simulations indicated that groundwater pumpage exceeded recharge, resulting in an estimated net cumulative depletion of groundwater storage (discharge minus recharge) of about 465,000 acre-feet from 1942 to 2016. The model simulated as much as 7.5 feet (90 inches; 2,286 millimeters) of aquifer compaction, which indicates the extensive fine-grained deposits and measured subsidence near Lucerne Lake.

California↗

Effects of sanitary sewers on ground-water levels and streams in Nassau and Suffolk Counties, New York; part 3: Development and application of southern Nassau County model

By 1990, sanitary sewers in Nassau County Sewage Disposal Districts 2 and 3 and Suffolk County Southwest Sewer District will discharge to the ocean 140 cu ft of water per second that would otherwise be returned to the groundwater system through septic tanks and similar systems. To evaluate the effects of this loss on groundwater levels and streamflow, the U.S. Geological Survey developed a groundwater flow model that couples a fine-scale subregional model to a regional model of a larger scale. The regional model generates flux boundary conditions for the subregional model, and the subregional model provides detail in the area of concern. Results indicate that the water table will decline by as much as 90% from conditions in the early 1970's. This report is one of a three-part series describing the predicted hydrologic effects of sewers in southern Nassau and southwestern Suffolk Counties.

New York↗