USGS Science⌕ Search

SEARCH · USGS Science

Results for “Cell Reports Sustainability”

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.

38 records · Page 3Linked to original sources

Groundwater-flow model and effects of projected groundwater use in the Ozark Plateaus Aquifer System in the vicinity of Greene County, Missouri — 1907-2030

Recent and historical periods of rapid growth have increased the stress on the groundwater resources in the Ozark aquifer in the Greene County, Missouri area. Historical pumpage from the Ozark aquifer has caused a cone of depression beneath Springfield, Missouri. In an effort to ease its dependence on groundwater for supply, the city of Springfield built a pipeline in 1996 to bring water from Stockton Lake to the city. Rapid population growth in the area coupled with the expanding cone of depression raised concern about the sustainability of groundwater as a resource for future use. A groundwater-flow model was developed by the U.S. Geological Survey in cooperation with Greene County, Missouri, the U. S. Army Corps of Engineers, and the Missouri Department of Natural Resources to assess the effect that increased groundwater demand is having on the long-term availability of groundwater in and around Greene County, Missouri. Three hydrogeologic units were represented in the groundwater-flow model: the Springfield Plateau aquifer, the Ozark confining unit, and the Ozark aquifer. The Springfield Plateau aquifer is less than 350 feet thick in the model area and generally is a low yield aquifer suitable only for domestic use. The Ozark aquifer is composed of a more than 900-foot thick sequence of dolomite and sandstone in the model area and is the primary aquifer throughout most of southern Missouri. Wells open to the entire thickness of the Ozark aquifer typically yield 1,000 gallons per minute or more. Between the two aquifers is the Ozark confining unit composed of as much as 98 feet of shale and limestone. Karst features such as sinkholes, springs, caves, and losing streams are present in both aquifers, but the majority of these features occur in the Springfield Plateau aquifer. The solution-enlarged fracture and bedding plane conduits in the karst system, particularly in the Springfield Plateau aquifer, are capable of moving large quantities of groundwater through the aquifer in relatively short periods of time. Pumpage rates in the model area increased from 1,093,268 cubic feet per day in 1962 to 2,693,423 cubic feet per day in 1987 to 4,330,177 cubic feet per day in 2006. Annual precipitation ranged from 25.21 inches in 1953 to 62.45 inches in 1927 from 1915 to 2006 in the model area. Recharge to the model was calculated as 2.53 percent of the annual precipitation and was varied annually. Recharge was distributed over the model area based on land slope and was adjusted in the city limits of Springfield to account for the impervious surface. A groundwater model with annual stress periods from 1907 to 2030 was developed using a transient calibration period from 1987 to 2006 and a prediction period from 2007 to 2030 to simulate flow in the Springfield Plateau aquifer and the Ozark aquifer. For the model area of approximately 2,870 square miles, the model hydrogeologic units and hydraulic properties were discretized into 253 rows, 316 columns, and 3 layers with the layer boundaries crossing hydrogeologic unit boundaries in some areas. The horizontal cell spacing was 1,000 feet by 1,000 feet. The model was calibrated by minimizing the difference between simulated head and observed water levels and simulated and observed flows in rivers and springs. Population and the associated groundwater use were estimated for 12 communities and the unincorporated area of Greene County based on past growth. Each was analyzed individually, and a low and high annual rate of growth relative to the 2006 population was computed for each community or group. Low growth rates ranged from 0.215 percent per year in Springfield to 6.997 percent per year in Rogersville. Total growth from 2006 to 2030 at the low growth rate ranged from 5.2 percent in Springfield to 167.9 percent in Rogersville. High growth rates ranged from 0.236 percent per year in Springfield to 7.345 percent per year in Rogersville. Total growth from 2006 to 2030 at the high growth rate ranged from 5.7 percent in Springfield to 176.3 percent in Rogersville. Response of the flow system to selected hypothetical pumping stresses and recharge conditions was simulated using the calibrated model. Seven hypothetical scenarios were simulated from 2007 to 2030 to test the effects of various stresses on the head in the Ozark aquifer. Hypothetical scenario 1 continued the 2006 pumping rates without change to the end of 2030. Scenario 2 assumed a low population growth rate with a 4-year drought at the beginning of the prediction period. Scenario 3 assumed a low population growth rate with a 4-year drought at the end of the prediction period. Scenario 4 assumed a high population growth rate with a 4-year drought at the beginning of the prediction period. Scenario 5 assumed a high population growth rate with a 4-year drought at the end of the prediction period. Scenario 6 and 7 had one new industrial well installed within the city limits of Springfield and one new industrial well installed about 3.5 miles east of Rogersville. Scenario 6 assumed a low population growth rate and scenario 7 assumed a high population growth rate. Results were compared by examining differences in head at the end of the simulation period. All scenarios examined resulted in potentiometric-surface declines from 2006 levels. Results from scenario 1 indicated that even with no increase in pumping, the potentiometric surface in the Springfield area continued to decline. The maximum decline of approximately 62 feet from the 2006 potentiometric surface occurred in Springfield. The maximum decline from the 2006 potentiometric surface in scenarios 2 and 3 was approximately 203 feet and in scenarios 4 and 5 was approximately 207 feet. The drought occurring at the end of the simulation period tended to broaden the drawdown area relative to the drought at the beginning. Drought timing did not substantially affect the potentiometric surface in the Ozark aquifer except for where the Ozark aquifer was exposed. Although not a substantial difference, the high population growth rate scenarios tended to have larger declines than the low population growth rate scenarios. As in the previous scenarios, little difference was noted between the low and high growth rate in scenario 6 and 7. Scenarios 6 and 7 showed declines of more than 640 feet from the 2006 potentiometric surface at the new well located in Springfield. The drawdown at the new wells decreased relatively quickly with increased distance from the well. Simulated head in the nearby cities of Nixa, Ozark, and Republic was nearly the same for scenarios 2 through 7 and was lower than the head predicted for scenario 1. Results from scenarios 2 through 7 indicate that the potentiometric surface in 2030 near these cities could decline 100 feet or more from the 2006 levels. Because model layers 2 and 3, representing the Ozark confining unit and most of the thickness of the Ozark aquifer, were simulated as confined, drawdown in the wells in the area of the Ozark aquifer that is unconfined or becomes unconfined during the simulation period will likely be under predicted.

Missouri↗

Groundwater-flow and land-subsidence model of Antelope Valley, California

Antelope Valley, California, is a topographically closed basin in the western part of the Mojave Desert, about 50 miles northeast of Los Angeles. The Antelope Valley groundwater basin is about 940 square miles and is separated from the northern part of Antelope Valley by faults and low-lying hills. Prior to 1972, groundwater provided more than 90 percent of the total water supply in the valley; since 1972, it has provided between 50 and 90 percent. Most groundwater pumping in the valley occurs in the Antelope Valley groundwater basin, which includes the rapidly growing cities of Lancaster and Palmdale. Groundwater-level declines of more than 270 feet in some parts of the groundwater basin have resulted in an increase in pumping lifts, reduced well efficiency, and land subsidence of more than 6 feet in some areas. Future urban growth and limits on the supply of imported water may increase reliance on groundwater. In 2011, the Los Angeles County Superior Court of California ruled that the Antelope Valley groundwater basin is in overdraft—groundwater extractions are in excess of the Court-defined safe yield of the groundwater basin. The Court determined that the safe yield of the adjudicated area of the basin was 110,000 acre-feet per year (acre-ft/yr). Natural recharge is an important component of total groundwater recharge in Antelope Valley; however, the exact quantity and distribution of natural recharge, primarily in the form of mountain-front recharge, is uncertain, with total estimates ranging from 30,000 to 160,000 acre-ft/yr. Technical experts, retained by parties to the adjudication, used 60,000 acre-ft/yr to estimate the sustainable yield of the basin, and this value was used in this study. In order to better understand the uncertainty associated with natural recharge and to provide a tool to aid in groundwater management, a numerical model of groundwater flow and land subsidence in the Antelope Valley groundwater basin was developed using old and new geohydrologic information. The groundwater-flow system consists of three aquifers: the upper, middle, and lower aquifers. The three aquifers, which were identified on the basis of the hydrologic properties, age, and depth of the unconsolidated deposits, consist of gravel, sand, silt, and clay alluvial deposits and clay and silty clay lacustrine deposits. Prior to groundwater development in the valley, recharge was primarily the infiltration of runoff from the surrounding mountains. Groundwater flowed from the recharge areas to discharge areas around the playas where it discharged from the aquifer system as either evapotranspiration or from springs. Partial barriers to horizontal groundwater flow, such as faults, have been identified in the groundwater basin. Water-level declines owing to groundwater development have eliminated the natural sources of discharge, and pumping for agricultural and urban uses have become the primary source of discharge from the groundwater system. Infiltration of return flow from agricultural irrigation has become an important source of recharge to the aquifer system. The groundwater-flow model of the basin was discretized horizontally into a grid of 130 rows and 118 columns of square cells 1 kilometer (0.621 mile) on a side, and vertically into four layers representing the upper (two layers), middle (one layer), and lower (one layer) aquifers. Faults that were thought to act as horizontal-flow barriers were simulated in the model. The model was calibrated to simulate steady-state conditions, represented by 1915 water levels and transient-state conditions during 1915–95, by using water-level and subsidence data. Initial estimates of the aquifer-system properties and stresses were obtained from a previously published numerical model of the Antelope Valley groundwater basin; estimates also were obtained from recently collected hydrologic data and from results of simulations of groundwater-flow and land-subsidence models of the Edwards Air Force Base area. Some of these initial estimates were modified during model calibration. Groundwater pumpage for agriculture was estimated on the basis of irrigated crop acreage and crop consumptive-use data. Pumpage for public supply, which is metered, was compiled and entered into a database used for this study. Estimated annual agricultural pumpage peaked at 395,000 acre-feet (acre-ft) in 1951 and then declined because of declining agricultural production. Recharge from irrigation return flows was assumed to be 30 percent of agricultural pumpage; delays associated with return flow moving through the unsaturated zone were also simulated. The annual quantity of mountain-front recharge initially was based on estimates from previous studies. The model was calibrated using the PEST software suite; prior information from the area was incorporated through the use of Tikhonov regularization. During model calibration, the estimated mountain-front recharge was reduced from the previous estimate of 30,300 acre-ft/yr to 29,150 acre-ft/yr. Results of the simulations using the calibrated model indicate that simulated groundwater pumpage exceeded recharge in most years, resulting in an estimated cumulative depletion in groundwater storage of 8,700,000 acre-ft during the transient-simulation period (1915–2005). About 15,000,000 acre-ft of cumulative groundwater pumpage was simulated during the transient-simulation period (1915–2005), reaching a maximum rate of about 400,000 acre-ft/yr in 1951. Groundwater pumpage resulted in simulated hydraulic heads declining by more than 150 feet (ft) compared to 1915 conditions in agricultural areas. The decline in hydraulic head in the groundwater basin is the result of this depletion of groundwater storage. In turn, the simulated decline in hydraulic head in the groundwater basin has resulted in the decrease in natural discharge from the basin and has caused compaction of aquitards, resulting in land subsidence. The areal distribution of total simulated land subsidence for 2005, after about 90 years of groundwater development, indicates that land subsidence occurred throughout almost the entire Lancaster subbasin, with a maximum of about 9.4 ft in the central and eastern parts of the subbasin. An important objective of this study was to systematically address the uncertainty in estimates of natural recharge and related aquifer parameters by using the groundwater-flow and land-subsidence model with observational data and expert knowledge. After the model was calibrated to the observations and a reasonable parameter set obtained, the parameter null space—parameter values that do not appreciably affect the model calibration but may have importance for prediction—was identified. The effect of parameter uncertainty on the estimation of mountain-front recharge was addressed using the Null-Space Monte Carlo method. The Pareto trade-off method of visualizing uncertainty was also used to portray the reasonableness of larger natural-recharge rates. Results indicate that the total mountain-front recharge likely ranges between 28,000 and 44,000 acre-ft/yr, which is appreciably less than published estimates of 60,000 acre-ft/yr. Additionally, expected errors associated with agricultural pumpage estimates used in this study were found to have relatively little effect on the estimates of mountain-front recharge, reflecting the difficulty in increasing recharge through manipulation of other components of the water budget. The calibrated model was used to simulate the response of the aquifer to potential future pumping scenarios: (1) no change in the distribution of pumpage, or status quo; (2) redistribution of pumpage; and (3) artificial recharge. All three of these scenarios specify a total pumpage throughout the Antelope Valley of 110,000 acre-ft/yr according to the safe yield value ruled by the Los Angeles County Superior Court of California. This reduction in groundwater pumpage is assumed uniform throughout the basin, based on a 10-percent reduction of the total pumpage in 2005 to achieve the 110,000 acre-ft/yr level. The calibrated Antelope Valley groundwater-flow and land-subsidence model was used to simulate the hydrologic effects of the three groundwater-management scenarios during a 50-year period by using the reduced, temporally constant, pumpage distribution. Results from the first scenario indicated that the total drawdown observed since predevelopment would continue, with values exceeding 325 ft near Palmdale; consequently, land subsidence would also continue, with additional subsidence (since 2005) exceeding 3 ft in the central part of the Lancaster subbasin. The second scenario evaluated redistributing pumpage from areas in the Lancaster subbasin where simulated hydraulic-head declines were the greatest to areas where declines were smallest. Neither a formal optimization algorithm nor water-rights allocations were considered when redistributing the pumpage. Results indicated that hydraulic heads near Palmdale, where the pumpage was reduced, would recover by about 200 ft compared to 2005 conditions, with only 30 ft of additional drawdown in the northwestern part of the Lancaster subbasin, where the pumpage was increased. The magnitude of the simulated additional land subsidence decreased slightly compared to the first, status quo, scenario but land subsidence continued to be simulated throughout most of the northern part of the Lancaster subbasin. The third scenario consisted of two artificial-recharge simulations along the Upper Amargosa Creek channel and at a site located north of Antelope Buttes. Results indicate that applying artificial recharge at these sites would yield continued drawdowns and associated land subsidence. However, the magnitudes of drawdown and subsidence would be smaller than those simulated in the status quo scenario, indicating that artificial-recharge operations in the Antelope Valley could be expected to reduce the magnitude and extent of continued water-level declines and associated land subsidence.

California↗