USGS ScienceSearch

Geology topics

J.H. Ellis

Publications and source records attributed to J.H. Ellis.

4 recordsLinked to original sources

Streamflow, base flow, and precipitation trends and simulated effects of Rush Springs aquifer groundwater withdrawals on base flows upgradient from Fort Cobb Reservoir, western Oklahoma

To better understand the relation between groundwater use in the Rush Springs aquifer and inflows to the Fort Cobb Reservoir, the U.S. Geological Survey, in cooperation with the Bureau of Reclamation, used a previously published numerical groundwater-flow model and historical streamflow records to evaluate four scenarios to investigate how changing groundwater withdrawals could affect base flows in streams that flow into Fort Cobb Reservoir. These scenarios consisted of observing simulated base-flow response by (1) scaling the 20-year equal-proportionate-share groundwater-withdrawal rate by various percentages over a 50-year period; (2) scaling the historical groundwater-withdrawal rates by various percentages across the entire Rush Springs aquifer; (3) scaling the historical groundwater-withdrawal rates within various subareas (zones) of the Fort Cobb Reservoir surface watershed; and (4) simulating a base-flow-depletion scenario. Cobb, Lake, and Willow Creeks are the major streams upgradient from the Fort Cobb Reservoir (listed from highest to lowest mean annual base flow). The results of scenarios 1 and 2 indicated that Willow Creek is the most susceptible to drying, but Cobb Creek was the most likely to have reduced base flow. Scenarios 3 and 4 indicated that groundwater withdrawals affect Cobb Creek base flows over a broader watershed area compared to Lake and Willow Creeks. In scenario 4, Cobb Creek base-flow depletion was higher across a larger area than Lake Creek and Willow Creek. Groundwater withdrawals in the Cobb Creek watershed tended to affect total inflows into Fort Cobb Reservoir more than other areas in the extent of the Rush Springs aquifer.

Oklahoma

Hydrogeology, land-surface subsidence, and documentation of the Gulf Coast Land Subsidence and Groundwater-Flow (GULF) model, southeast Texas, 1897–2018

Executive Summary As a part of the Texas Water Development Board groundwater availability modeling program, the U.S. Geological Survey developed the Gulf Coast Land Subsidence and Groundwater-Flow model (hereinafter, the “GULF model”) and ensemble to simulate groundwater flow and land-surface subsidence in the northern part of the Gulf Coast aquifer system (the study area) in Texas from predevelopment (1897) through 2018. Since the publication of a previous groundwater model for the greater Houston area in 2012, there have been changes to the distribution of groundwater withdrawals and advances in modeling tools. To reflect these changes and to simulate more recent conditions, the GULF model was developed in cooperation with the Harris-Galveston and Fort Bend Subsidence Districts to provide an updated Groundwater Availability Model. Since the early 1900s, most of the groundwater withdrawals in the study area have been from three of the hydrogeologic units that compose the Gulf Coast aquifer system—the Chicot, Evangeline, and Jasper aquifers and, more recently, from the Catahoula confining unit. Withdrawals from these hydrogeologic units are used for municipal supply, commercial and industrial use, and irrigation purposes. Withdrawals of large quantities of groundwater in the greater Houston area have caused widespread groundwater-level declines in the Chicot, Evangeline, and Jasper aquifers of more than 300 feet (ft). Early development of the aquifer system, which began before 1900, resulted in nearly 50 percent of the eventual historical groundwater-level minimums having been reached as early as 1946 in some areas. These groundwater-level declines led to more than 9 ft of land-surface subsidence—historically in central and southeastern Harris County and Galveston County, but more recently in northern, northwestern, and western Harris County, Montgomery County, and northern Fort Bend County—from depressurization and compaction of clay and silt layers interbedded in the aquifer sediments. In a generalized conceptual model of the Gulf Coast aquifer system, water enters the groundwater system in topographically high outcrops of the hydrogeologic units in the northwestern part of the aquifer system. Groundwater that does not discharge to streams flows to intermediate and deep zones of the aquifer system southeastward of the outcrop areas where it is discharged by wells and by upward leakage in topographically low areas near the coast. The uppermost parts of the aquifer system, which include outcrop areas, are under water-table (unconfined) conditions where the groundwater is not confined under pressure. As depth increases in the aquifer system and interbedded clay and silt layers accumulate, water-table conditions evolve into confined conditions where the groundwater is under pressure. Groundwater flow and land-surface subsidence in the GULF model and ensemble were simulated by using MODFLOW 6 with the Skeletal Storage, Compaction, and Subsidence package. The model consists of six layers, one for each of the five hydrogeologic units in the northern part of the Gulf Coast aquifer system and a surficial top layer that includes part of each hydrogeologic unit. Transient groundwater flow was simulated during 1897–2018 by using a combination of multiyear, annual, and monthly stress periods. An initial steady-state stress period was configured to represent predevelopment mean annual inflows and outflows. The subsidence package used in the GULF model and ensemble uses a head-based subsidence formulation that simulates the delayed drainage response from clay and silt sediment to changes in groundwater levels. The GULF model and ensemble were history matched to groundwater-level observations at selected wells, land-surface subsidence at benchmarks, aquifer compaction at borehole extensometers, and vertical displacement from Global Positioning System stations. A Bayesian framework was used to represent uncertainty in modeled parameters and simulated outputs of interest. History matching and uncertainty quantification were performed by using a Monte Carlo approach enabled through iterative ensemble smoother software to produce an ensemble of models fit to historical data. The iterative ensemble smoother substantially reduced the computational demand of parameter estimation by approximating the first-order relation between model inputs and outputs, thereby allowing 183,207 adjustable parameters to be used for history matching at a relatively low computational and time cost. The history-matched parameter values are within the ranges of previously published values and agree with the current understanding of the spatial and temporal patterns of parameter uncertainty for the Gulf Coast aquifer system. A good agreement between the observed (or estimated) and simulated groundwater levels, land-surface subsidence, compaction, and vertical displacement was obtained across the modeled area based on qualitative and quantitative comparisons. Ensemble mean annual groundwater-flow rates to the Chicot, Evangeline, Jasper aquifers and Catahoula confining unit were 0.0–0.49 inch (in.), 0.09–0.33 in., 0.01–0.07 in., and 0.01–0.05 in., respectively. GULF model mean annual groundwater-flow rates to the Chicot, Evangeline, and Jasper aquifers and Catahoula confining unit were 0.31 in., 0.19 in., 0.03 in., and 0.03 in., respectively. The GULF-model-simulated recharge to the outcrop area was the largest inflow (75 percent), and recharge to other areas was 25 percent of the model inflow. The simulated outflows included (1) net surface-water/groundwater exchange with study area streams (50 percent), (2) groundwater use (49 percent), and (3) net surface-water/groundwater exchange with the Gulf of Mexico (1 percent). The sum of the simulated values of the outflows (1,041,973 acre-feet per year [acre-ft/yr]) and the elastic expansion of the fine-grained sediment and numerical solver error (339 acre-ft/yr) minus the inflows (654,172 acre-ft/yr) represents the reduction of storage from the Gulf Coast aquifer system (388,140 acre-ft/yr). Most of the storage depletion is caused by the long-term groundwater-level declines that have resulted primarily in inelastic compaction. The GULF model was used to estimate Jasper aquifer compaction at selected benchmarks in Montgomery County and northern Harris County, which are the primary locations of Jasper aquifer groundwater use. Simulated Jasper aquifer compaction in northern Harris County was between 0.2 and 0.5 ft, or between about 5 and 16 percent of simulated subsidence at the benchmark locations. Simulated Jasper aquifer compaction in Montgomery County was between 0.8 and 1.2 ft, or between about 33 and 57 percent of simulated subsidence at the benchmark locations.

Texas

Evaluating the effects of downscaled climate projections on groundwater storage and simulated base-flow contribution to the North Fork Red River and Lake Altus, southwest Oklahoma (USA)

Potential effects of projected climate variability on base flow and groundwater storage in the North Fork Red River aquifer, Oklahoma (USA), were estimated using downscaled climate model data coupled with a numerical groundwater-flow model. The North Fork Red River aquifer discharges groundwater to the North Fork Red River, which provides inflow to Lake Altus. To approximate future conditions, Coupled Model Intercomparison Project Phase 5 climate data were downscaled to the watershed and a time-series of scaling factors were developed and interpolated for three climate scenarios (central tendency, warmer and drier, and less warm and wetter) representing future climate conditions for the period 2045–2074. These scaling factors were then applied to a soil-water-balance model to produce groundwater recharge and evapotranspiration estimates. A MODFLOW groundwater-flow model of the North Fork Red River aquifer used the scaled recharge and evapotranspiration data to estimate changes in base flow and water-surface elevation of Lake Altus. Compared to a baseline scenario, the mean percent change in annual base flow during 2045–2074 was −10.8 and −15.9% for the central tendency and warmer/drier scenarios, respectively; the mean percent change in annual base flow for the less-warm/wetter scenario was +15.7%. The mean annual percent change in groundwater storage for the central tendency, warmer/drier, and less-warm/wetter climate scenarios and the baseline are −2.7, −3.2, and +3.0%, respectively. The range of outcomes from the climate scenarios may be influenced by variability in the downscaled climate data for precipitation more than for temperature.

Oklahoma

Simulation of groundwater flow and analysis of projected water use for the Rush Springs aquifer, western Oklahoma

The U.S. Geological Survey, in cooperation with the Bureau of Reclamation and the Oklahoma Water Resources Board, (1) quantified the groundwater resources of the Rush Springs aquifer in western Oklahoma by developing a numerical groundwater-flow model, (2) evaluated the effects of estimated equal-proportionate-share (EPS) pumping rates on aquifer storage and streamflow for time periods of 20, 40, and 50 years into the future, (3) assessed the uncertainty in the EPS scenario results, and (4) evaluated the effects of (a) projected groundwater-use rates extended 50 years into the future and (b) sustained hypothetical drought conditions over a 10-year period on stream base flow and groundwater in storage. The Rush Springs aquifer is an important source of water for municipal and irrigation use by many communities and agricultural users in the study area. The study area is composed of about 4,970 square miles (3,181,003 acres) of Rush Springs aquifer bedrock deposits located in 14 counties. The study area also includes the alluvium and terrace deposits of the Canadian and Washita Rivers, as well as alluvium along the Little Washita River, Deer Creek, and a number of smaller tributaries of the Washita River that overlie the bedrock. A numerical groundwater-flow model of the Rush Springs aquifer was constructed by using MODFLOW with the Newton solver. Groundwater flow was simulated for January 1979–December 2015 by using monthly stress periods, and an initial steady-state stress period was configured to represent mean annual inflows and outflows. The model was calibrated to groundwater-level observations at selected wells, monthly base flow at nine streamgages, stream seepage as estimated for the conceptual water budget, and Fort Cobb Reservoir stage. The EPS scenarios for the Rush Springs aquifer were run for periods of 20, 40, and 50 years. The 20-, 40-, and 50-year EPS pumping rates under normal recharge conditions were 0.82, 0.49, and 0.43 acre-foot per acre per year, respectively. Given the 2,954,545-acre aquifer area used for the EPS scenarios, the 20-year rate corresponds to an annual yield of about 2,422,727 acre-feet per year. Groundwater storage at the end of the 20-year EPS scenario was about 13,321,000 acre-feet, or about 31,516,437 acre-feet (70 percent) less than the starting EPS scenario storage. This decrease in storage was equivalent to a mean groundwater-level decline of about 152 feet. Water availability under the EPS pumping rate was primarily from the western area of the model. Saturation was sustained though the entire EPS scenario where the aquifer was sufficiently thick or a shallow hydraulic gradient was present. Fort Cobb Reservoir stage was below the dead-pool stage after about 5 years of 20-year EPS pumping. An uncertainty analysis was conducted to assess the uncertainty in the EPS scenario results. An ensemble of 400 random sets of possible parameter values was performed for the uncertainty analysis by using a multivariate normal distribution centered on the calibrated parameter values. The parameter bounds for the uncertainty analysis were determined by using the posterior covariance matrix, which allows for the incorporation of knowledge gained during the calibration process as well as observation uncertainty and the correlation between estimated parameters. The uncertainty results indicate a 95-percent confidence interval for the 20-year EPS pumping rate between 0.73 and 0.95 acre-foot per acre per year. Projected 50-year pumping scenarios were used to simulate the effects of selected well withdrawal rates on groundwater storage of the Rush Springs aquifer. The effects of well withdrawals were evaluated by comparing changes in groundwater storage between four 50-year scenarios using (1) no groundwater use, (2) mean groundwater use for the study period (1979–2015), (3) increasing groundwater use, and (4) groundwater use at the 2015 rate. The increasing-use scenario assumed a 38-percent increase in pumping over 50 years on the basis of 2010–60 demand projections for western Oklahoma. Simulated groundwater storage changes ranged between an increase of 6.3 percent for the scenario with no groundwater use, and 0.9 percent for the scenario with 2015 groundwater-use rates. For the Fort Cobb Reservoir surface watershed, simulated groundwater storage changes ranged between an increase of 23.6 percent for the scenario with no groundwater use and a decrease of 4.0 percent for the increasing groundwater-use scenario. Groundwater-level changes were generally greater in areas with a large concentration of groundwater wells and groundwater use such as the Fort Cobb Reservoir surface watershed. A hypothetical 10-year drought scenario was used to simulate the effects of a prolonged period of reduced recharge on the Rush Springs aquifer groundwater storage and Fort Cobb Reservoir stage and storage. Drought effects were quantified by comparing the results of the drought scenario to those of the calibrated numerical model. To simulate the hypothetical drought, recharge in the calibrated numerical model was reduced by 50 percent during the simulated drought period (1983–1992), and upstream inflows to the Canadian and Washita Rivers and associated tributaries were reduced by 37 percent. Groundwater storage at the end of the hypothetical drought period in December 1992 was about 42,983,000 acre-feet, or about 3,525,000 acre-feet (7.6 percent) less than the groundwater storage of the calibrated numerical model. This change in groundwater storage is equivalent to a mean groundwater-level decline of 15.8 feet. Simulated mean base-flow declines at the Canadian and Washita River streamgages were between 39 and 59 percent during the drought period. The minimum stage in Fort Cobb Reservoir at the end of the hypothetical drought period was 1,311 feet, indicating a storage capacity of only 10 percent of active conservation pool storage. The Fort Cobb Reservoir storage declines mostly resulted from reduced base flows in Cobb, Lake, and Willow Creeks upstream from the reservoir.

Oklahoma