Evaluation of ground-water flow and land-surface subsidence caused by hypothetical withdrawals in the northern part of the Gulf Coast Aquifer system, Texas

Mark C. Kasmarek; Brian D. Reece; Natalie A. Houston

doi:10.3133/sir20055024

Home/Archive/USGS

USGS · sir20055024

Evaluation of ground-water flow and land-surface subsidence caused by hypothetical withdrawals in the northern part of the Gulf Coast Aquifer system, Texas

Mark C. Kasmarek·Brian D. Reece·Natalie A. Houston

View original source ↗Read source document ↗Export metadata ↓

Abstract

During 2003–04 the U.S. Geological Survey, in cooperation with the Texas Water Development Board (TWDB) and the Harris-Galveston Coastal Subsidence District (HGCSD), used the previously developed Northern Gulf Coast Ground-Water Availability Modeling (NGC GAM) model to evaluate the effects of hypothetical projected withdrawals on ground-water flow in the northern part of the Gulf Coast aquifer system and land-surface subsidence in the NGC GAM model area of Texas. The Gulf Coast aquifer system comprises, from the surface, the Chicot and Evangeline aquifers, the Burkeville confining unit, the Jasper aquifer, and the Catahoula confining unit. Two withdrawal scenarios were simulated. The first scenario comprises historical withdrawals from the aquifer system for 1891–2000 and hypothetical projected withdrawals for 2001–50 compiled by the TWDB (TWDB scenario). The projected withdrawals compiled by the TWDB are based on ground-water demands estimated by regional water planning groups. The second scenario is a “merge” of the TWDB scenario with an alternate set of projected withdrawals from the Chicot and Evangeline aquifers in the Houston metropolitan area for 1995–2030 provided by the HGCSD (HGCSD scenario). Under the TWDB scenario withdrawals from the entire system are projected to be about the same in 2050 as in 2000. The simulated potentiometric surfaces of the Chicot aquifer for 2010, 2020, 2030, 2040, and 2050 show relatively little change in configuration from the simulated 2000 potentiometric surface (maximum water-level depths in southern Harris County 150–200 feet below NGVD 29). The simulated decadal potentiometric surfaces of the Evangeline aquifer show the most change between 2000 and 2010. The area of water levels 250– 400 feet below NGVD 29 in western Harris County in 2000 shifts southeastward to southern Harris County, and water levels recover to 200–250 feet below NGVD 29 by 2010. Water levels in southern Harris County recover to 150–200 feet below NGVD 29 by 2020 and remain in that range through 2050. A relatively small cone of depression in southern Montgomery County that did not appear in the 2000 surface develops and enlarges during the projected period, with a maximum depth of 250–300 feet below NGVD 29 in 2030, 2040, and 2050. The simulated decadal potentiometric surfaces of the Jasper aquifer each have a major cone of depression centered in southern Montgomery County that was minimally developed in 2000 but reaches depths of 550–650 feet below NGVD 29 in the 2020, 2030, 2040, and 2050 surfaces. Under the TWDB scenario the percentage of withdrawals supplied by net recharge increases from 75 percent in 2000 to 87 percent in 2050, and the percentage of withdrawals supplied by storage decreases from 25 percent in 2000 to 13 percent in 2050. Under the HGCSD scenario, withdrawals from the Chicot and Evangeline aquifers increase about 74 percent during 1995–2030; Jasper aquifer withdrawals are unchanged from those of the TWDB scenario. For the 2010, 2020, and 2030 potentiometric surfaces of the Chicot and Evangeline aquifers, the substantially greater withdrawals of the HGCSD scenario relative to those of the TWDB scenario result in progressively deeper cones of depression than those in the potentiometric surfaces associated with the TWDB scenario—for the Chicot aquifer in southern Harris County, 400–450 feet below NGVD 29 in 2030; for the Evangeline aquifer in southern Montgomery County, 700–750 feet below NGVD 29 in 2030. Although Jasper aquifer withdrawals are the same for both scenarios, the major cone of depression centered in southern Montgomery County in the 2030 potentiometric surface is 50 feet deeper at its center (600–700 feet below NGVD 29) than the cone in the 2030 surface under the TWDB scenario. Under the HGCSD scenario, the percentage of withdrawals supplied by net recharge decreases from 72 percent in 1995 to 57 percent in 2030, and the percentage of withdrawals supplied by storage increases from 28 percent in 2000 to 43 percent in 2030. About 85 percent of the increase supplied by storage is from the compaction of clay. Land-surface subsidence in the major area of subsidence centered in Harris and Galveston Counties during 2000–50 that results from simulating the TWDB withdrawal scenario expands slightly to the west and increases in places. The maximum change occurs in the Conroe area where subsidence increases from about 4 to about 13 feet during the projected period. Land-surface subsidence in the major area of subsidence during 1995–2030 that results from simulating the HGCSD withdrawal scenario increases substantially. For example, in east-central Harris County maximum subsidence increases from about 10–11 feet in 1995 to 22 feet in 2030. The hypothetical projected withdrawal scenarios are estimates of future withdrawals and might not represent actual future withdrawals. The simplifying assumptions that the downdip limit of freshwater flow in each hydrogeologic unit is a stable, sharp interface across which no flow occurs and that the base of the system is a no-flow boundary become less realistic and thus increase the uncertainty in results as drawdowns increase. The presence of uncertainty dictates that the results of the predictive simulations described in this report be used with caution in any decision-making process.

Explore related subjects

Scientific Investigations Report Report

Source-reported bounding extent: 28.110748760633534° to 31.89621446335144° latitude; -97.734375° to -92.867431640625° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Mark C. Kasmarek, Brian D. Reece, Natalie A. Houston. 2005. Evaluation of ground-water flow and land-surface subsidence caused by hypothetical withdrawals in the northern part of the Gulf Coast Aquifer system, Texas. https://doi.org/10.3133/sir20055024

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

Evaluating the number of realizations (train-test splits) necessary to ensure robust estimates of conventional hydrothermal resource favorability using Monte Carlo cross-validation, the Great Basin, Western USA

Supervised machine learning approaches used to produce favorability maps in resource assessments commonly partition an input dataset into training and testing subsets (that is, through train-test splits). However, if a single split is not representative of underlying resource conditions, the resulting favorability map may be biased, particularly with sparse or heterogeneous datasets. To mitigate the potential for unrepresentative favorability maps, recent works by U.S. Geological Survey researchers have used Monte Carlo cross-validation to implement a large number ( n) of randomized train-test splits to produce n realizations. The realizations are then aggregated using cell-wise median predictions to produce the resulting final favorability map, recognizing that the median will remove extreme values that might have been estimated from a comparatively few unrepresentative realizations. Theoretically, the median prediction at every cell converges to a single estimate as n approaches infinity, and the question addressed herein is: How many realizations are necessary to get a robust estimate (that is, an estimate unaffected by additional realizations)? Therefore, convergence metrics are described to identify the number of train-test splits and realizations needed to make a reliable map using the data and supervised machine learning workflow for conventional hydrothermal resource assessment of the Great Basin of the Western United States. For this study, robust prediction occurs with n approximately equal to 600 realizations, and smaller ensembles (for example, n is approximately equal to 200 or even 100) may be suitable during preliminary or exploratory analyses. Although the resulting threshold numbers will vary depending on the data and modeling approach used, the convergence metrics proposed and applied herein are applicable to other supervised machine learning studies that use Monte Carlo cross-validation.

Scientific Investigations Report↗

USGS2026

Groundwater budget for the Portneuf River Basin, southeastern Idaho

The U.S. Geological Survey, in cooperation with the Idaho Department of Water Resources, developed the first groundwater budget for the upper Portneuf River Basin in southeastern Idaho, a tributary of the Snake River. A groundwater budget for the Portneuf River Basin supports water resource and administrative management decisions within the basin and for downstream users of the eastern Snake River Plain aquifer. The analysis quantifies inflows, outflows, and residuals for an aquifer system defined by the Idaho Department of Water Resources under three recent precipitation conditions: an average year (2016), a dry year (2021), and a wet year (2023). The budget was calculated for the basin and for two subbasins delineated by the U.S. Geological Survey streamgage Portneuf River at Topaz, Idaho (13073000): the above Topaz subbasin, which includes the Portneuf-Gem valley, and the below Topaz subbasin, which includes Marsh Valley. Aquifer inflows include mountain-front recharge; areal recharge; recharge from applied water, including recharge from irrigation, canal seepage, and effluent returns; groundwater flow into the basin; and reservoir seepage. Outflows include groundwater pumping for agricultural, municipal, and domestic uses; subsurface outflow through the Portneuf gap; and discharge to streams as base flow. Estimated annual inflows for the full basin ranged from 116,880 acre-feet in 2021 to 154,880 acre-feet in 2023, and outflows ranged from 132,190 acre-feet in 2016 to 158,490 acre-feet in 2023. Basinwide residuals, the change in storage plus uncertainty, were −2,410 acre-feet in 2016, −18,410 acre-feet in 2021, and −3,610 acre-feet in 2023, and each subbasin had a mixture of negative and positive residuals. These basinwide results, interpreted as change in storage plus water-budget uncertainty, indicate minor net losses that increase with dry conditions, consistent with well depth-to-water measurements indicating gradual multi-decadal declines. Patterns in subbasin residuals are driven by irrigation practices and surface-water management; for example, the above Topaz subbasin relies more on groundwater for irrigation, resulting in more negative residuals during dry years. In contrast, the below Topaz subbasin depends more on surface water, producing positive residuals during dry years because of irrigation losses and canal seepage. The analyses for both subbasins indicate decades-long groundwater declines, which agree with the overall negative residuals calculated. The largest sources of uncertainty include climate inputs (precipitation and evapotranspiration) and poorly constrained components such as groundwater inflow and tributary stream infiltration. These findings provide insight into how climatic variables and irrigation practices affect groundwater inflows and outflows within the Portneuf River Basin.

Idaho↗

USGS2026

World minerals outlook—Antimony, copper, graphite, niobium, potash, silicon, silver, and tungsten through 2030

Demand for mineral commodities continues to increase owing to rapid technological advancements and global economic growth. To support this increase in demand, information about current and near-term supplies of these mineral commodities is essential. Such data help evaluate the health and stability of global supply chains and protect against potential supply risks. Current production and production capacity data coupled with information regarding near-term expected changes (such as facility closures, planned expansions, or the completion of new mineral production projects) can provide insight into expected future market developments and expected available supply of mineral commodities. This report within the “World Minerals Outlook” series analyzes data for antimony, copper, natural graphite, niobium, potash, silicon metal, silver, and tungsten—eight mineral commodities listed in the 2025 List of Critical Minerals. World production data from 2019 through 2025, as well as estimated world capacity data from 2023 through 2025 were analyzed. Using these data, near-term (2026–2030) expected world capacities were estimated. For antimony, copper, niobium, silicon metal, silver, and tungsten, modest increases in production capacity of 4 to 10 percent were expected between 2025 and 2030, whereas capacities for natural graphite and potash were expected to increase by about 18 to 20 percent.

Scientific Investigations Report↗