Search USGSSearch

SEARCH · Search USGS

Results for “Texas Water Development Board 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.

74 records · Page 5Linked to original sources

Hydrogeology, water use, and simulation of flow in the High Plains aquifer in northwestern Oklahoma, southeastern Colorado, southwestern Kansas, northeastern New Mexico, and northwestern Texas

The U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, began a three-year study of the High Plains aquifer in northwestern Oklahoma in 1996. The primary purpose of this study was to develop a ground-water flow model to provide the Water Board with the information it needs to manage the quantity of water withdrawn from the aquifer. The study area consists of about 7,100 square miles in Oklahoma and about 20,800 square miles in adjacent states to provide appropriate hydrologic boundaries for the flow model. The High Plains aquifer includes all sediments from the base of the Ogallala Formation to the potentiometric surface. The saturated thickness in Oklahoma ranges from more than 400 feet to less than 50 feet. Natural recharge to the aquifer from precipitation occurs throughout the area but is extremely variable. Dryland agricultural practices appear to enhance recharge from precipitation, and part of the water pumped for irrigation also recharges the aquifer. Natural discharge occurs as discharge to streams, evapotranspiration where the depth to water is shallow, and diffuse ground-water flow across the eastern boundary. Artificial discharge occurs as discharge to wells. Irrigation accounted for 96 percent of all use of water from the High Plains aquifer in the Oklahoma portion of the study area in 1992 and 93 percent in 1997. Total estimated water use in 1992 for the Oklahoma portion of the study area was 396,000 acre-feet and was about 3.2 million acre-feet for the entire study area. Since development of the aquifer, water levels have declined more than 100 feet in small areas of Texas County, Oklahoma, and more than 50 feet in areas of Cimarron County. Only a small area of Beaver County had declines of more than 10 feet, and Ellis County had rises of more than 10 feet. A flow model constructed using the MODFLOW computer code had 21,073 active cells in one layer and had a 6,000- foot grid in both the north-south and east-west directions. The model was used to simulate the period before major development of the aquifer and the period of development. The model was calibrated using observed conditions available as of 1998. The predevelopment-period model integrated data or estimates on the base of aquifer, hydraulic conductivity, streambed and drain conductances, and recharge from precipitation to calculate the predevelopment altitude of the water table, discharge to the rivers and streams, and other discharges. Hydraulic conductivity, recharge, and streambed conductance were varied during calibration so that the model produced a reasonable representation of the observed water table altitude and the estimated discharge to streams. Hydraulic conductivity was reduced in the area of salt dissolution in underlying Permianage rocks. Recharge from precipitation was estimated to be 4.0 percent of precipitation in greater recharge zones and 0.37 percent in lesser recharge zones. Within Oklahoma, the mean difference between water levels simulated by the model and measured water levels at 86 observation points is -2.8 feet, the mean absolute difference is 44.1 feet, and the root mean square difference is 52.0 feet. The simulated discharge is much larger than the estimated discharge for the Beaver River, is somewhat larger for Cimarron River and Wolf Creek, and is about the same for Crooked Creek. The development-period model added specific yield, pumpage, and recharge due to irrigation and dryland cultivation to simulate the period 1946 through 1997. During calibration, estimated specific yield was reduced by 15 percent in Oklahoma east of the Cimarron-Texas County line. Simulated recharge due to irrigation ranges from 24 percent for the 1940s and 1950s to 2 percent for the 1990s. Estimated recharge due to dryland cultivation is about 3.9 percent of precipitation. The mean difference between the simulated and observed waterlevel changes from predevelopment to 1998 at 162 observation points in Oklahoma is less than 0.01 foot, the mean absolute difference is 13.1 feet, and the root mean square difference is 17.9 feet. The model simulates 7.8 cubic feet per second discharge to the Beaver River above Optima Reservoir at the end of 1997, whereas the river was actually dry to this point by this time. The model simulates a decrease in discharge to the Cimarron River and Wolf Creek that appears to be reasonable. The sensitivity of the predevelopment-period model to changes in recharge and hydraulic conductivity was tested. Simulated water levels are sensitive to changes in both recharge and hydraulic conductivity. Simulated discharge to streams is sensitive to recharge but is insensitive to hydraulic conductivity. Recharge and hydraulic conductivity are closely related with respect to water levels but are not related with respect to discharge to streams. The sensitivity of the development-period model to changes in specific yield and recharge due to dryland cultivation was tested. The model appears more sensitive to specific yield, but the two inputs were not varied over the same range in percent change. The calibrated development-period model was used to simulate water-level changes from 1998 to 2020 using mean 1996-97 pumpage. The largest simulated water-level changes in Oklahoma occur in Texas County where water levels are simulated to decline 25 to 50 additional feet over a large area. Water levels also are simulated to decline 10 to 25 additional feet in two large areas of Cimarron County, two areas in Beaver County, and one area in Ellis County, all in Oklahoma. Water levels are simulated to decline more than 100 additional feet in several areas in Kansas and were simulated to decline 50 to 100 additional feet in several areas in Texas.

Colorado, Kansas, New Mexico, Oklahoma, Texas

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