Geology topics

Natalie A. Houston

Publications and source records attributed to Natalie A. Houston.

At least 19 recordsLinked to original sources

USGS2021

Estimates of water use associated with continuous oil and gas development in the Permian Basin, Texas and New Mexico, 2010–19, with comparisons to the Williston Basin, North Dakota and Montana

The Permian Basin, in west Texas and southeastern New Mexico is one of the largest conventional oil and gas reservoirs in the United States and is becoming one of the world’s largest continuous oil and gas (COG) reservoirs. Advances in technology have enabled oil and gas to be extracted from reservoirs that historically were developed using conventional, or vertical, well drilling techniques. Conventional oil and gas reservoirs have discrete deposits that are well defined and are typically trapped by an overlying geologic formation or caprock, whereas COG reservoirs contain deposits that are distributed evenly throughout the geologic formation, typically have much lower permeability (the capacity of a porous rock to transmit a fluid) than the conventional deposits, and require specialized horizontal extraction techniques. The methods to extract the oil from the two different reservoirs require differing amounts of water, and the horizontal extraction methods typically require substantially more water. In 2015, the U.S. Geological Survey started a topical study to quantify water used during COG development. The Permian Basin, which contains both types of reservoirs (continuous and conventional), was the second basin in the United States in the U.S. Geological Survey’s topical study to quantify water used during COG development.

Montana, New Mexico, North Dakota, Texas↗

USGS2021

Estimates of water use associated with continuous oil and gas development in the Permian Basin, Texas and New Mexico, 2010–19

In 2015, the U.S. Geological Survey started a topical study to quantify water use in areas of continuous oil and gas (COG) development. The first phase of the study was completed in 2019 and analyzed the Williston Basin. The second phase of the study analyzed the Permian Basin using the same techniques and approaches used for the Williston Basin analysis. The Permian Basin was selected for the second phase of water-use analysis for the following reasons: (1) the basin has the largest undiscovered technically recoverable oil and gas resource in the United States, (2) the basin has a continuous resource in tight shale that primarily produces oil, and (3) the basin is within the contiguous United States. This study used data from 60 counties in Texas and New Mexico with spatial coverage based on the Permian Basin extent defined by the U.S. Energy Information Administration, a representation of the geologically defined Permian Basin. Data from several sources were used in the analysis of direct, indirect, and ancillary water use associated with COG development in the Permian Basin and are available in an associated data release. Hydraulic fracturing water-use data were used to determine the start of the recent (before 2019) COG development boom in oil production in the Permian Basin in the same way that the data were used for the Williston Basin study. Water-use data were aggregated by county and year, which were the sampling units used in the analysis. The water-use analysis of the Permian Basin contained three elements: (1) estimates of water use, in million gallons, by county and year; (2) coefficients of water use from regression models, in million gallons per developed oil and gas well; and (3) performance (based on goodness-of-fit metrics) of the regression models in estimating the observed water use. Coefficients from the linear and quantile regression models of direct, indirect, and ancillary water use in the Permian Basin were produced as aggregate values for the counties and years. The mean estimate of direct water use had a 95-percent confidence interval of 4.13–5.45 million gallons (Mgal) per developed oil and gas well. The coefficient from the linear regression model of indirect water use was 0.111 Mgal per well, with a 95-percent confidence interval of 0.104–0.117 Mgal per well. The mean estimate of ancillary water use in the Permian Basin was 1.09 Mgal per well, with a 95-percent confidence interval of 1.05–1.13 Mgal per well. Model performance was evaluated with goodness-of-fit metrics including coefficient of determination ( R 2 ), root mean square error, and the ratio of root mean square error to standard deviation of observations computed from leave-one-out cross validation of the linear and quantile regression models of direct, indirect, and ancillary water use. Model performance for direct water use was acceptable, with an R 2 value of 0.91. The model performance of indirect water use was acceptable, with an R 2 value of 0.89. Values of R 2 for the ancillary water-use categories were at least 0.89. Annual mean estimates for hydraulic fracturing, cementing, drilling, indirect, and ancillary water use per well for the years 2010–17 were comparable between the Permian and Williston Basins. Hydraulic fracturing water use increased similarly from 2010 to 2015 in the Permian Basin and the Williston Basin, increasing from 0.6 Mgal per well in 2010 to 5.4 Mgal per well in 2015 in the Permian Basin and from 1.4 Mgal per well in 2010 to 4.7 Mgal per well in 2015 in the Williston Basin. By design, the Permian water-use assessment is a simplification of a complex and continually developing system and therefore has uncertainty and limitations in the interpretation of results. Despite the modeling limitations, the results summarized in the report, when compared to other studies, compare well with water-use estimations. The favorable comparison highlights the transferability of the water-use methodology to other areas of COG development.

New Mexico, Texas↗

USGS2021

Hydrogeologic framework and groundwater characterization in selected alluvial basins in the upper Rio Grande basin, Colorado, New Mexico, and Texas, United States, and Chihuahua, Mexico, 1980 to 2015

Increasing demand for the limited water resources of the United States continues to put pressure on resource management agencies to balance the competing needs of ecosystem health with municipal, agricultural, and other uses. To meet these needs, the U.S. Geological Survey conducted a multiyear study to evaluate water resources in the upper Rio Grande Basin in the southwestern United States. The upper Rio Grande Basin extends from south-central Colorado, through New Mexico, into west Texas near Fort Quitman, including parts of Chihuahua, Mexico. The upper Rio Grande Basin consists of a sequence of alluvial basins that formed in the Rio Grande rift approximately 30 million years ago. This report describes the hydrogeology of the upper Rio Grande Basin and how the groundwater resources in the basin have changed from 1980 to 2015. The hydrogeologic framework includes the horizontal delineation of the alluvial basins within the upper Rio Grande Basin from the headwaters in Colorado to Fort Quitman, Texas, including part of Mexico. Groundwater-level measurements from existing State and Federal data were used to construct groundwater-level altitude and groundwater-level change maps. Of the 2,699 wells with groundwater-level data used in this study, 1,055 wells had data for only a single 5-year period, 703 wells had data for 50 percent or more of the 35 years of the study, and only 57 wells have 5-year groundwater-level data for the entire study period. The median decline in water levels in the upper Rio Grande Basin was 0.13 foot (ft) per 5-year period, and declines were measured in 53 percent of the 703 wells that contained data for 50 percent or more of the study period. Rates of groundwater-level decline greater than 1 ft per 5-year period were measured in 17 percent of the wells, greater than 2 ft per 5-year period, in 3 percent of the wells, and greater than 3 ft per 5-year period, in 1 percent of the wells. Overall, groundwater levels rose in 6 percent of the 703 wells that contained data for 50 percent or more of the study period, and in 4 percent of the wells, groundwater levels rose by 1 ft or more per 5-year period. Groundwater-level changes in wells with consecutive 5-year measurement periods exhibited the most variability in the Española, Middle Rio Grande, and Mesilla/Conejos-Médanos alluvial basins. The largest declines in groundwater-level altitudes in individual wells were observed in the Española alluvial basin during 1995–2000, in the Palomas alluvial basin during 2010–2015, and in the Jornada del Muerto alluvial basin during 2005–10. The largest rises in groundwater-level altitudes in individual wells were observed in the Española alluvial basin during 2005–10, in the Middle Rio Grande alluvial basin during 1995–2000, and in the Mesilla/Conejos-Médanos alluvial basin during 1980–85. Changes in groundwater storage throughout the study period varied by alluvial basin, likely based largely on changes in groundwater withdrawals because of increased demands during drier periods and population growth. All alluvial basins except the Tularosa-Hueco alluvial basin were evaluated for changes in groundwater storage from 1980 to 2015. Extremely limited data availability in 2010–15 for the Tularosa-Hueco alluvial basin led to this 5-year period being dropped from the groundwater-level change map and storage analysis for this basin. In the San Luis Valley in southern Colorado, efforts to reverse groundwater depletion in the unconfined aquifer recovered approximately 250,000 acre-feet in storage between late 2013 and early 2018, following the implementation of a “pay-to-pump” groundwater program. However, severe drought that persists in the upper Rio Grande Basin, particularly in southern Colorado, has undone some of the conservation efforts. Within the Española alluvial basin, groundwater storage varied because municipal demand increased the demand on groundwater resources and conservation efforts were implemented. A groundwater-flow model evaluated for the Española alluvial basin indicated declines in groundwater storage from 1947 through 1982. Groundwater storage decreased in the Española alluvial basin in 1980–85, 1985–90, 1990–95, 1995–2000, and 2005–10 and increased in 2000–05 and 2010–15 leading to groundwater storage in 2015 about even with that in 1985. Based on gridded groundwater-level altitudes, groundwater storage decreased in the Middle Rio Grande Basin from 1980 to 2015, except for during the 1980–85, 2000–05, and 2010–15 periods with an overall cumulative storage decrease from 1980 to 2015. Groundwater-flow models evaluated for the Middle Rio Grande alluvial basin showed groundwater storage in the Middle Rio Grande alluvial basin has been reduced since the mid-1950s through the end of the study period except for a brief recovery (reduction in storage outflow) in the mid-1980s. Simulated groundwater storage has also decreased in parts of the Palomas and Mesilla/Conejos-Médanos alluvial basins, and the northern part of the Conejos-Médanos alluvial basin starting in 1995 (excluding 2005 and 2007) and in the Tularosa-Hueco alluvial basin from the early 1940s to the end of the study period. Groundwater storage increased in the Mesilla/Conejos-Médanos alluvial basin during 1980–85 and slightly during 1990–95 and then decreased in the other 5-year periods. Groundwater storage in the Tularosa-Hueco alluvial basin increased from 1985 to 1990, but otherwise decreased, leading to an overall net groundwater-level decline in this part of the basin from 1980 to 2010.

Colorado, New Mexico, Texas, Chihuahua↗

USGS2021

Risk-based decision-support groundwater modeling for the lower San Antonio River Basin, Texas, USA

A numerical surface-water/groundwater model was developed for the lower San Antonio River Basin to evaluate the responses of low base flows and groundwater levels within the basin under conditions of reduced recharge and increased groundwater withdrawals. Batch data assimilation through history matching used a simulation of historical conditions (2006-2013); this process included history-matching to groundwater levels and base-flow estimates at several gages, and was completed in a high-dimensional (highly parameterized) framework. The model was developed in an uncertainty framework such that parameters, observations, and scenarios of interest are envisioned stochastically as distributions of potential values. Results indicate that groundwater contributions to surface water during periods of low flow may be reduced from 6% to 25% with a corresponding 25% reduction in recharge and a 25% increase in groundwater pumping over an 8-year planning period. Furthermore, results indicate groundwater-level reductions in some hydrostratigraphic units are more likely than in other hydrostratigraphic units over an 8-year period under drought conditions with the higher groundwater withdrawal scenario.

Texas↗

USGS2019

Pecos River Basin salinity assessment, Santa Rosa Lake, New Mexico, to the confluence of the Pecos River and the Rio Grande, Texas, 2015

The elevated salinity of the Pecos River throughout much of its length is of paramount concern to water users and water managers. Dissolved-solids concentrations in the Pecos River exceed 3,000 milligrams per liter in many of its reaches in the study area, from Santa Rosa Lake, New Mexico, to the confluence of the Pecos River with the Rio Grande, Texas. The salinity of the Pecos River increases downstream and affects the availability of useable water in the Pecos River Basin. In this report, “salinity” and “dissolved-solids concentration” are considered synonymous; both terms are used to refer to the total ionic concentration of dissolved minerals in water. The sources of salinity in the Pecos River Basin are natural (geologic) and anthropogenic, including but not limited to groundwater discharge, springs, and irrigation return flows. Previous studies in the Pecos River Basin were project specific and designed to address salinity issues in specific parts of the basin; therefore, in 2015, the U.S. Geological Survey in cooperation with the U.S. Army Corps of Engineers, New Mexico Interstate Stream Commission, Texas Commission on Environmental Quality, and Texas Water Development Board assessed the major sources of salinity throughout the extent of the basin where elevated salinity in the Pecos River is well documented (that is, in the drainage area of the Pecos River from Santa Rosa Lake to the confluence of the Pecos River and the Rio Grande). The goal was to gain a better understanding of how specific areas might be contributing to the elevated salinity in the Pecos River and how salinity of the Pecos River has changed over time. This assessment includes a literature review and compilation of previously published salinity-related data, which guided the collection of additional water-quality samples and streamflow gain-loss measurements. Differences in water quality of surface-water and groundwater samples, streamflow measurements, and geophysical data were assessed to gain new insights regarding sources of salinity in the Pecos River Basin and a more detailed assessment of potential areas of elevated salinity in the basin. The datasets compiled for this assessment are available in a companion data release. The literature review identified several potential sources of salinity inputs to the Pecos River in New Mexico and Texas. In New Mexico, sources of salinity inputs included sinkhole springs discharging into El Rito Creek, the Bitter Lake National Wildlife Refuge inflow to the Pecos River, inflow from the Rio Hondo, including the main channel and a restored channel at the Bitter Lake National Wildlife Refuge referred to as the “Rio Hondo spring channel,” the outflow from Lea Lake at Bottomless Lakes State Park, and the Malaga Bend region of the Pecos River. In Texas, sources of salinity inputs included Salt Creek downstream from Red Bluff Reservoir and the area near the Horsehead Crossing ford on the Pecos River. The compilation of historical water-quality data revealed a lack of consistent sampling of the same constituents at the same sites along the main stem of the Pecos River, which results in data gaps that hinder the ability to effectively analyze long-term changes in water quality that may help with the understanding of how salinity in the Pecos River has changed over time and identifying the sources of salinity in the Pecos River Basin. To help fill these data gaps, water-quality and streamflow data were collected in the study area in February 2015 by the U.S. Geological Survey. Historical water-quality data and newly collected data from February 2015 were evaluated for selected major-ion concentrations, dissolved-solids concentrations, and deuterium, oxygen, and strontium isotopes. Analysis of the data indicated several areas of increasing salinity in the Pecos River. Most notable increases were in two subreaches of the river, between Acme, N. Mex., and Artesia, N. Mex., and between Orla, Tex., and Grandfalls, Tex. Increasing sodium and chloride concentrations from Acme to Artesia coincided with changes in isotopic ratios within the Pecos River Basin. Changes in isotopic ratios in this reach indicate a likely inflow from an isotopically different source of water compared to the water in the main stem of the Pecos River, such as groundwater inflow, inflow from surface-water features distinct from the main stem of the Pecos River, or both. In the subreach between Orla and Grandfalls, an increase in dissolved-solids concentrations was observed along with a shift in isotope values, indicating that neither evaporative processes in Red Bluff Reservoir nor inflow from Salt Creek likely solely influences the salinity of the Pecos River in this subreach. The highest dissolved-solids concentrations in the Pecos River Basin were measured downstream from Grandfalls, where dissolved-solids concentrations are greater than 16,000 milligrams per liter near Iraan, Tex. Changes in isotopic values (deuterium, oxygen, and strontium) indicate mixing of different waters at several areas along the main stem of the Pecos River. The spatial distribution of the areas of interest from the literature review and the water-quality data are available in the companion data release.

Texas, New Mexico↗

USGS2015

Water-level altitudes 2015 and water-level changes in the Chicot, Evangeline, and Jasper aquifers and compaction 1973-2014 in the Chicot and Evangeline aquifers, Houston-Galveston region, Texas

Most of the land-surface subsidence in the Houston-Galveston region, Texas, has occurred as a direct result of groundwater withdrawals for municipal supply, commercial and industrial use, and irrigation that depressured and dewatered the Chicot and Evangeline aquifers, thereby causing compaction of the aquifer sediments, mostly in the fine-grained silt and clay layers. This report, prepared by the U.S. Geological Survey in cooperation with the Harris-Galveston Subsidence District, City of Houston, Fort Bend Subsidence District, Lone Star Groundwater Conservation District, and Brazoria County Groundwater Conservation District, is one in an annual series of reports depicting water-level altitudes and water-level changes in the Chicot, Evangeline, and Jasper aquifers and measured cumulative compaction of subsurface sediments in the Chicot and Evangeline aquifers in the Houston-Galveston region. The report contains regional-scale maps depicting approximate 2015 water-level altitudes (represented by measurements made during December 2014–March 2015) for the Chicot, Evangeline, and Jasper aquifers; maps depicting 1-year (2014–15) water-level changes for each aquifer; maps depicting approximate contoured 5-year (2010–15) water-level changes for each aquifer; maps depicting approximate contoured long-term (1990–2015 and 1977–2015) water-level changes for the Chicot and Evangeline aquifers; a map depicting approximate contoured long-term (2000–15) water-level changes for the Jasper aquifer; a map depicting locations of borehole-extensometer sites; and graphs depicting measured cumulative compaction of subsurface sediments at the borehole extensometers during 1973–2014. Three tables listing the water-level data used to construct each water-level map for each aquifer and a table listing the measured cumulative compaction data for each extensometer site and graphs are included. In 2015, water-level-altitude contours for the Chicot aquifer ranged from 175 feet (ft) below the vertical datum (the National Geodetic Vertical Datum of 1929 or the North American Vertical Datum of 1988; hereinafter, datum) in a localized area in northwestern Harris County to 200 ft above datum in northern and western Montgomery County. Water-level changes for 2014–15 in the Chicot aquifer ranged from a 24-ft decline to a 31-ft rise. Contoured 5-year and long-term water-level changes in the Chicot aquifer ranged from a 40-ft decline to a 40-ft rise (2010–15), from a 100-ft decline to a 100-ft rise (1990–2015), and from a 100-ft decline to a 200-ft rise (1977–2015). In 2015, water-level-altitude contours for the Evangeline aquifer ranged from 250 ft below datum in a localized area extending from south-central Montgomery County into north-central Harris County and in an additional area located in central Harris County to 200 ft above datum in southeastern Grimes and northwestern Montgomery Counties. Water-level changes for 2014–15 in the Evangeline aquifer ranged from a 66-ft decline to a 78-ft rise. Contoured 5-year and long-term water-level changes in the Evangeline aquifer ranged from a 60-ft decline to an 80-ft rise (2010–15), from a 200-ft decline to a 240-ft rise (1990–2015), and from a 320-ft decline to a 240-ft rise (1977–2015). In 2015, water-level-altitude contours for the Jasper aquifer ranged from 200 ft below datum in south-central Montgomery County that extends into north-central Harris County to 250 ft above datum in northwestern Montgomery County. Water-level changes for 2014–15 in the Jasper aquifer ranged from a 17-ft decline to a 35-ft rise. Contoured 5-year and long-term water-level changes in the Jasper aquifer ranged from a 60-ft decline to four small, localized areas of 10-ft rises (2010–15) and from a 220-ft decline to no change (2000–15). Compaction of subsurface sediments (mostly in the fine-grained silt and clay layers) composing the Chicot and Evangeline aquifers was recorded continuously by using analog technology at the 13 borehole extensometers at 11 sites that were either activated or installed between 1973 and 1980. For the period of record beginning in 1973 (or later depending on activation or installation date) and ending in December 2014, measured cumulative compaction at the 13 extensometers ranged from 0.101 ft at the Texas City-Moses Lake extensometer to 3.668 ft at the Addicks extensometer. During 2014, a total of 10 of the 13 extensometers recorded a slight net decrease of land-surface elevation; the extensometers at the Lake Houston and Clear Lake (shallow) sites recorded slight net increases of land-surface elevation, and the extensometer at the Texas City-Moses Lake site recorded no change in elevation. The rate of compaction varies from site to site because of differences in rates of groundwater withdrawal in the areas adjacent to each extensometer site and differences among sites in the ratios of sand, silt, and clay and compressibilities of the subsurface sediments. It is not appropriate, therefore, to extrapolate or infer a rate of compaction for an adjacent area on the basis of the rate of compaction measured at nearby extensometers.

Texas↗

USGS2015

Updated numerical model with uncertainty assessment of 1950-56 drought conditions on brackish-water movement within the Edwards aquifer, San Antonio, Texas

In 2010, the U.S. Geological Survey, in cooperation with the San Antonio Water System, began a study to assess the brackish-water movement within the Edwards aquifer (more specifically the potential for brackish-water encroachment into wells near the interface between the freshwater and brackish-water transition zones, referred to in this report as the transition-zone interface) and effects on spring discharge at Comal and San Marcos Springs under drought conditions using a numerical model. The quantitative targets of this study are to predict the effects of higher-than-average groundwater withdrawals from wells and drought-of-record rainfall conditions of 1950–56 on (1) dissolved-solids concentration changes at production wells near the transition-zone interface, (2) total spring discharge at Comal and San Marcos Springs, and (3) the groundwater head (head) at Bexar County index well J-17. The predictions of interest, and the parameters implemented into the model, were evaluated to quantify their uncertainty so the results of the predictions could be presented in terms of a 95-percent credible interval. The model area covers the San Antonio and Barton Springs segments of the Edwards aquifer; the history-matching effort was focused on the San Antonio segment. A previously developed diffuse-flow model of the Edwards aquifer, which forms the basis for the model in this assessment, is primarily based on a conceptualization in which flow in the aquifer is predominately through a network of numerous small fractures and openings. Primary updates to this model include an extension of the active area downdip, a conversion to an 8-layer SEAWAT variable-density flow and transport model to simulate dissolved-solids concentration effects on water density, history matching to 1999–2009 conditions, and parameter estimation in a highly parameterized context using automated methods in PEST (a model-independent Parameter ESTimation code). In addition to the best-fit parameter values derived from history matching, the uncertainty of model parameters was also estimated by using linear uncertainty analysis. Comparison of “prior” (before history matching) and “posterior” (after history matching) variances of parameters indicate that the information within the observation dataset used for history matching informs many parameters. The concentration threshold parameters were well-informed by the observation dataset as their posterior distributions were much narrower than their prior distributions. The transition-zone scaling parameters of hydraulic conductivity, effective porosity, and specific storage were all informed by the observation dataset, as evidenced by the difference between the prior and posterior variances. Saline-zone scaling parameters, alternatively, were not informed by the observation dataset for effective porosity and specific storage. Resulting posterior drier-month, wetter-month, and annual recharge multiplier parameter variances are important to understanding how well recharge is estimated and implemented within the model. The shifts of the posterior distributions left and right indicate that there were zones where less or more water was needed in the model. The widths of the distributions were not decreased substantially, indicating that many of the best-fit recharge parameters are not statistically different from the initial values specified in the history-matching effort. Recharge from rainfall is the driving force behind groundwater flow and heads in the aquifer; therefore, an increase in understanding of this process would benefit model development by potentially decreasing the uncertainty of this parameter. The history-matching effort was most helpful in informing the parameters in the model that control discharge at springs, namely, the spring orifice (drain) altitude and drain conductance parameters for each spring. The uncertainty assessment of the predictive model (a hypothetical recurrence of 1950–56 drought conditions and higher-than-average groundwater withdrawals from wells) provided insights into the potential effects of these conditions on dissolved-solids concentration changes at production wells near the transition-zone interface, discharges at Comal and San Marcos Springs, and heads at Bexar County index well J-17. Results at the 25 production wells near the transition-zone interface indicate that the uncertainty of model input parameters based on expert knowledge yielded an upper bound of the 95-percent credible interval of dissolved-solids concentrations that exceeds the secondary drinking water standards of 1,000 milligrams per liter (mg/L) of the Texas Commission on Environmental Quality (TCEQ) for many wells. However, the history-matching process provided key information to inform prediction-sensitive model parameters and therefore, contributed to a substantial decrease of the upper bound of the 95-percent credible interval to below the secondary drinking water standards. Reductions in dissolved-solids concentration changes were on the order of 400 mg/L to 1,300 mg/L. The reduction in uncertainty in regards to this prediction implies that this prediction of dissolved-solids concentration change can be made with some certainty using this current model and that those parameters that control this prediction are informed by the observation dataset. Even though predictive uncertainty was reduced for this prediction, dissolved-solids concentration changes were still greater than zero, indicating a minimal increase in concentration at these 25 production wells during the 7-year simulation period is likely. However, this minimal concentration increase indicates a small potential for movement of the brackish-water transition zone near these wells during the 7-year simulation period of drought-ofrecord (1950–56) rainfall conditions with higher-than-average groundwater withdrawals by wells. Predictive results of total spring discharge during the 7-year period, as well as head predictions at Bexar County index well J-17, were much different than the dissolved-solids concentration change results at the production wells. These upper bounds are an order of magnitude larger than the actual prediction which implies that (1) the predictions of total spring discharge at Comal and San Marcos Springs and head at Bexar County index well J-17 made with this model are not reliable, and (2) parameters that control these predictions are not informed well by the observation dataset during historymatching, even though the history-matching process yielded parameters to reproduce spring discharges and heads at these locations during the history-matching period. Furthermore, because spring discharges at these two springs and heads at Bexar County index well J-17 represent more of a cumulative effect of upstream conditions over a larger distance (and longer time), many more parameters (with their own uncertainties) are potentially controlling these predictions than the prediction of dissolved-solids concentration change at the prediction wells, and therefore contributing to a large posterior uncertainty.

Texas↗

USGS2014

Simulation of groundwater flow in the Edwards-Trinity and related aquifers in the Pecos County region, Texas

The Edwards-Trinity aquifer is a vital groundwater resource for agricultural, industrial, and public supply uses in the Pecos County region of western Texas. The U.S. Geological Survey completed a comprehensive, integrated analysis of available hydrogeologic data to develop a numerical groundwater-flow model of the Edwards-Trinity and related aquifers in the study area in parts of Brewster, Jeff Davis, Pecos, and Reeves Counties. The active model area covers about 3,400 square miles of the Pecos County region of Texas west of the Pecos River, and its boundaries were defined to include the saturated areas of the Edwards-Trinity aquifer. The model is a five-layer representation of the Pecos Valley, Edwards-Trinity, Dockum, and Rustler aquifers. The Pecos Valley aquifer is referred to as the alluvial layer, and the Edwards-Trinity aquifer is divided into layers representing the Edwards part of the Edwards-Trinity aquifer and the Trinity part of the Edwards-Trinity aquifer, respectively. The calibration period of the simulation extends from 1940 to 2010. Simulated hydraulic heads generally were in good agreement with observed values; 1,684 out of 2,860 (59 percent) of the simulated values were within 25 feet of the observed value. The average root mean square error value of hydraulic head for the Edwards-Trinity aquifer was 34.2 feet, which was approximately 4 percent of the average total observed change in groundwater-level altitude (groundwater level). Simulated spring flow representing Comanche Springs exhibits a pattern similar to observed spring flow. Independent geochemical modeling corroborates results of simulated groundwater flow that indicates groundwater in the Edwards-Trinity aquifer in the Leon-Belding and Fort Stockton areas is a mixture of recharge from the Barilla and Davis Mountains and groundwater that has upwelled from the Rustler aquifer. The model was used to simulate groundwater-level altitudes resulting from prolonged pumping to evaluate sustainability of current and projected water-use demands. Each of three scenarios utilized a continuation of the calibrated model. Scenario 1 extended recent (2008) irrigation and nonirrigation pumping values for a 30-year period from 2010 to 2040. Projected groundwater-level changes in and around the Fort Stockton area under scenario 1 change little from current conditions, indicating that the groundwater system is near equilibrium with respect to recent (2008) pumping stress. Projected groundwater-level declines in the eastern part of the model area ranging from 5.0 to 15.0 feet are likely the result of nonequilibrium conditions associated with recent increases in pumping after a prolonged water-level recovery period of little or no pumping. Projected groundwater-level declines (from 15.0 to 31.0 feet) occurred in localized areas by the end of scenario 1 in the Leon-Belding area. Scenario 2 evaluated the effects of extended recent (2008) pumping rates as assigned in scenario 1 with year-round maximum permitted pumping rates in the Belding area. Results of scenario 2 are similar in water-level decline and extent as those of scenario 1. The extent of the projected groundwater-level decline in the range from 5.0 to 15.0 feet in the Leon-Belding irrigation area expanded slightly (about a 2-percent increase) from that of scenario 1. Maximum projected groundwater-level declines in the Leon-Belding irrigation area were approximately 31.3 feet in small isolated areas. Scenario 3 evaluated the effects of periodic increases in pumping rates over the 30-year extended period. Results of scenario 3 are similar to those of scenario 2 in terms of the areas of groundwater-level decline; however, the maximum projected groundwater-level decline increased to approximately 34.5 feet in the Leon-Belding area, and the extent of the decline was larger in area (about a 17-percent increase) than that of scenario 2. Additionally, the area of projected groundwater-level declines in the eastern part of the model area increased from that of scenario 2—two individual areas of decline coalesced into one larger area. The localized nature of the projected groundwater-level declines is a reflection of the high degree of fractured control on storage and hydraulic conductivity in the Edwards-Trinity aquifer. Additionally, the finding that simulated spring flow is highly dependent on the transient nature of hydraulic heads in the underlying aquifer indicates the importance of adequately understanding and characterizing the entire groundwater system.

Texas↗

USGS2013

Geodatabase compilation of hydrogeologic, remote sensing, and water-budget-component data for the High Plains aquifer, 2011

The High Plains aquifer underlies almost 112 million acres in the central United States. It is one of the largest aquifers in the Nation in terms of annual groundwater withdrawals and provides drinking water for 2.3 million people. The High Plains aquifer has gained national and international attention as a highly stressed groundwater supply primarily because it has been appreciably depleted in some areas. The U.S. Geological Survey has an active program to monitor the changes in groundwater levels for the High Plains aquifer and has documented substantial water-level changes since predevelopment: the High Plains Groundwater Availability Study is part of a series of regional groundwater availability studies conducted to evaluate the availability and sustainability of major aquifers across the Nation. The goals of the regional groundwater studies are to quantify current groundwater resources in an aquifer system, evaluate how these resources have changed over time, and provide tools to better understand a systems response to future demands and environmental stresses. The purpose of this report is to present selected data developed and synthesized for the High Plains aquifer as part of the High Plains Groundwater Availability Study. The High Plains Groundwater Availability Study includes the development of a water-budget-component analysis for the High Plains completed in 2011 and development of a groundwater-flow model for the northern High Plains aquifer. Both of these tasks require large amounts of data about the High Plains aquifer. Data pertaining to the High Plains aquifer were collected, synthesized, and then organized into digital data containers called geodatabases. There are 8 geodatabases, 1 file geodatabase and 7 personal geodatabases, that have been grouped in three categories: hydrogeologic data, remote sensing data, and water-budget-component data. The hydrogeologic data pertaining to the northern High Plains aquifer is included in three separate geodatabases: (1) base data from a groundwater-flow model; (2) hydrogeology and hydraulic properties data; and (3) groundwater-flow model data to be used as calibration targets. The remote sensing data for this study were developed by the U. S. Geological Survey Earth Resources Observation and Science Center and include historical and predicted land-use/land-cover data and actual evapotranspiration data by using remotely sensed temperature data. The water-budget-component data contains selected raster data from maps in the “Selected Approaches to Estimate Water-Budget Components of the High Plains, 1940 Through 1949 and 2000 Through 2009” report completed in 2011 ( http://pubs.usgs.gov/sir/2011/5183/ ). Federal Geographic Data Committee compliant metadata were created for each spatial and tabular data layer in the geodatabases.

Colorado, Kansas, Nebraska, New Mexico, Oklahoma, ↗

USGS2013

A conceptual hydrogeologic model for the hydrogeologic framework, geochemistry, and groundwater-flow system of the Edwards-Trinity and related aquifers in the Pecos County region, Texas

The Edwards-Trinity aquifer is a vital groundwater resource for agricultural, industrial, and municipal uses in the Trans-Pecos region of west Texas. A conceptual model of the hydrogeologic framework, geochemistry, and groundwater-flow system in the 4,700 square-mile study area was developed by the U.S. Geological Survey (USGS) in cooperation with the Middle Pecos Groundwater Conservation District, Pecos County, City of Fort Stockton, Brewster County, and Pecos County Water Control and Improvement District No. 1. The model was developed to gain a better understanding of the groundwater system and to establish a scientific foundation for resource-management decisions. Data and information were collected or obtained from various sources to develop the model. Lithologic information obtained from well reports and geophysical data were used to describe the hydrostratigraphy and structural features of the groundwater system, and aquifer-test data were used to estimate aquifer hydraulic properties. Groundwater-quality data were used to evaluate groundwater-flow paths, water and rock interaction, aquifer interaction, and the mixing of water from different sources. Groundwater-level data also were used to evaluate aquifer interaction as well as to develop a potentiometric-surface map, delineate regional groundwater divides, and describe regional groundwater-flow paths. Several previous studies have been done to compile or collect physical and chemical data, describe the hydrogeologic processes, and develop conceptual and numerical groundwater-flow models of the Edwards-Trinity aquifer in the Trans-Pecos region. Documented methods were used to compile and collect groundwater, surface-water, geochemical, geophysical, and geologic information that subsequently were used to develop this conceptual model.

Texas↗

USGS2012

Data collection and compilation for a geodatabase of groundwater, surface-water, water-quality, geophysical, and geologic data, Pecos County Region, Texas, 1930-2011

The U.S. Geological Survey, in cooperation with Middle Pecos Groundwater Conservation District, Pecos County, City of Fort Stockton, Brewster County, and Pecos County Water Control and Improvement District No. 1, compiled groundwater, surface-water, water-quality, geophysical, and geologic data for site locations in the Pecos County region, Texas, and developed a geodatabase to facilitate use of this information. Data were compiled for an approximately 4,700 square mile area of the Pecos County region, Texas. The geodatabase contains data from 8,242 sampling locations; it was designed to organize and store field-collected geochemical and geophysical data, as well as digital database resources from the U.S. Geological Survey, Middle Pecos Groundwater Conservation District, Texas Water Development Board, Texas Commission on Environmental Quality,and numerous other State and local databases. The geodatabase combines these disparate database resources into a simple data model. Site locations are geospatially enabled and stored in a geodatabase feature class for cartographic visualization and spatial analysis within a Geographic Information System. The sampling locations are related to hydrogeologic information through the use of geodatabase relationship classes. The geodatabase relationship classes provide the ability to perform complex spatial and data-driven queries to explore data stored in the geodatabase.

Texas↗

USGS2012

A conceptual model of the hydrogeologic framework, geochemistry, and groundwater-flow system of the Edwards-Trinity and related aquifers in the Pecos County region, Texas

A conceptual model of the hydrogeologic framework, geochemistry, and groundwater-flow system of the Edwards-Trinity and related aquifers, which include the Pecos Valley, Igneous, Dockum, Rustler, and Capitan Reef aquifers, was developed as the second phase of a groundwater availability study in the Pecos County region in west Texas. The first phase of the study was to collect and compile groundwater, surface-water, water-quality, geophysical, and geologic data in the area. The third phase of the study involves a numerical groundwater-flow model of the Edwards-Trinity aquifer in order to simulate groundwater conditions based on various groundwater-withdrawal scenarios. Resource managers plan to use the results of the study to establish management strategies for the groundwater system. The hydrogeologic framework is composed of the hydrostratigraphy, structural features, and hydraulic properties of the groundwater system. Well and geophysical logs were interpreted to define the top and base surfaces of the Edwards-Trinity aquifer units. Elevations of the top and base of the Edwards-Trinity aquifer generally decrease from the southwestern part of the study area to the northeast. The thicknesses of the Edwards-Trinity aquifer units were calculated using the interpolated top and base surfaces of the hydrostratigraphic units. Some of the thinnest sections of the aquifer were in the eastern part of the study area and some of the thickest sections were in the Pecos, Monument Draw, and Belding-Coyanosa trough areas. Normal-fault zones, which formed as growth and collapse features as sediments were deposited along the margins of more resistant rocks and as overlying sediments collapsed into the voids created by the dissolution of Permian-age evaporite deposits, were delineated based on the interpretation of hydrostratigraphic cross sections. The lowest aquifer transmissivity values were measured in the eastern part of the study area; the highest transmissivity values were measured in a faulted area of the Monument Draw trough. Hydraulic conductivity values generally exhibited the same trends as the transmissivity values. Groundwater-quality data and groundwater-level data were used in context with the hydrogeologic framework to assess the chemical characteristics of water from different sources, regional groundwater-flow paths, recharge sources, the mixing of water from different sources, and discharge in the study area. Groundwater-level altitudes generally decrease from southwest to northeast and regional groundwater flow is from areas of recharge south and west to the north and northeast. Four principal sources of recharge to the Edwards-Trinity aquifer were identified: (1) regional flow that originated as recharge northwest of the study area, (2) runoff from the Barilla, Davis, and Glass Mountains, (3) return flow from irrigation, and (4) upwelling from deeper aquifers. Results indicated Edwards-Trinity aquifer water in the study area was dominated by mineralized, regional groundwater flow that most likely recharged during the cooler, wetter climates of the Pleistocene with variable contributions of recent, local recharge. Groundwater generally flows into the down-dip extent of the Edwards-Trinity aquifer where it discharges into overlying or underlying aquifer units, discharges from springs, discharges to the Pecos River, follows a regional flow path east out of the study area, or is withdrawn by groundwater wells. Structural features such as mountains, troughs, and faults play a substantial role in the distribution of recharge, local and regional groundwater flow, spring discharge, and aquifer interaction.

Texas↗

USGS2011

Sources and Delivery of Nutrients to the Northwestern Gulf of Mexico from Streams in the South-Central United States

SPAtially Referenced Regressions On Watershed attributes (SPARROW) models were developed to estimate nutrient inputs [total nitrogen (TN) and total phosphorus (TP)] to the northwestern part of the Gulf of Mexico from streams in the South‐Central United States (U.S.). This area included drainages of the Lower Mississippi, Arkansas‐White‐Red, and Texas‐Gulf hydrologic regions. The models were standardized to reflect nutrient sources and stream conditions during 2002. Model predictions of nutrient loads (mass per time) and yields (mass per area per time) generally were greatest in streams in the eastern part of the region and along reaches near the Texas and Louisiana shoreline. The Mississippi River and Atchafalaya River watersheds, which drain nearly two‐thirds of the conterminous U.S., delivered the largest nutrient loads to the Gulf of Mexico, as expected. However, the three largest delivered TN yields were from the Trinity River/Galveston Bay, Calcasieu River, and Aransas River watersheds, while the three largest delivered TP yields were from the Calcasieu River, Mermentau River, and Trinity River/Galveston Bay watersheds. Model output indicated that the three largest sources of nitrogen from the region were atmospheric deposition (42%), commercial fertilizer (20%), and livestock manure (unconfined, 17%). The three largest sources of phosphorus were commercial fertilizer (28%), urban runoff (23%), and livestock manure (confined and unconfined, 23%).

Texas, Oklahoma, Arkansas, Louisiana, Missouri, Ka↗

USGS2011

Hydrogeology, chemical characteristics, and water sources and pathways in the zone of contribution of a public-supply well in San Antonio, Texas

In 2001, the National Water-Quality Assessment (NAWQA) Program of the U.S. Geological Survey initiated a series of studies on the transport of anthropogenic and natural contaminants (TANC) to public-supply wells (PSWs). The main goal of the TANC project was to better understand the source, transport, and receptor factors that control contaminant movement to PSWs in representative aquifers of the United States. Regional- and local-scale study areas were selected from within existing NAWQA study units, including the south-central Texas Edwards aquifer. The local-scale TANC study area, nested within the regional-scale NAWQA study area, is representative of the regional Edwards aquifer. The PSW selected for study is within a well field of six production wells. Although a single PSW was initially selected, because of constraints of well-field operation, samples were collected from different wells within the well field for different components of the study. Data collected from all of the well-field wells were considered comparable because of similar well construction, hydrogeology, and geochemistry. An additional 38 PSWs (mostly completed in the confined part of the aquifer) were sampled throughout the regional aquifer to characterize water quality. Two monitoring well clusters, with wells completed at different depths, were installed to the east and west of the well field (the Zarzamora and Timberhill monitoring well clusters, respectively). One of the monitoring wells was completed in the overburden to evaluate potential hydrologic connectivity with the Edwards aquifer. Geophysical and flowmeter logs were collected from one of the well-field PSWs to determine zones of contribution to the wellbore. These contributing zones, associated with different hydrogeologic units, were used to select monitoring well completion depths and groundwater sample collection depths for depth-dependent sampling. Depth-dependent samples were collected from the PSW from three different depths and under three different pumping conditions. Additionally, selected monitoring wells and one of the well-field PSWs were sampled several times in response to a rainfall and recharge event to assess short-term (event-scale) temporal variations in water quality. For comparison purposes, groundwater samples were categorized as being from regional aquifer PSWs, from the well field (wellhead samples), from the monitoring wells (excluding the overburden well), from the overburden well, from the PSW depth-dependent sampling, and from temporal sampling. Groundwater samples were analyzed for inorganic, organic, isotopic, and age-dating tracers to characterize geochemical conditions in the aquifer and provide understanding of the mechanisms of mobilization and movement of selected constituents from source areas to a PSW. Sources, tracers, and conditions used to assess water quality and processes affecting the PSW and the aquifer system included (1) carbonate host rock composition; (2) physicochemical constituents; (3) major and trace element concentrations; (4) saturation indices with respect to minerals in aquifer rocks; (5) elemental ratios, such as magnesium to calcium ratios, that are indicative of water-rock interaction processes; (6) oxidation-reduction conditions; (7) nutrient concentrations, in particular nitrate concentrations; (8) the isotopic composition of nitrate, which can point to specific nitrate sources; (9) strontium isotopes; (10) stable isotopes of hydrogen and oxygen; (11) organic contaminant concentrations, including pesticides and volatile organic compounds; (12) age tracers, apparent-age distribution, and dissolved gas data used in age interpretations; (13) depth-dependent water chemistry collected from the PSW under different pumping conditions to assess zones of contribution; and (14) temporal variability in groundwater composition from the PSW and selected monitoring wells in response to an aquifer recharge event. Geochemical results indicate that the well-field and monitoring well samples were largely representative of groundwater in the regional confined aquifer. Constituents of concern in the Edwards aquifer for the long-term sustainability of the groundwater resource include the nutrient nitrate and anthropogenic organic contaminants. Nitrate concentrations (as nitrogen) for regional aquifer PSWs had a median value of 1.9 milligrams per liter, which is similar to previously reported values for the regional aquifer. Nitrate-isotope compositions for groundwater samples collected from the well-field PSWs and monitoring wells had a narrow range, with values indicative of natural soil organic values. A comparison with historical nitrate-isotope values, however, suggests that a component of nitrate in groundwater from biogenic sources might have increased over the last 30 years. Several organic contaminants (the pesticide atrazine, its degradate deethylatrazine, trichloromethane (chloroform; a drinking-water disinfection byproduct), and the solvent tetrachloroethene (PCE)) were widely distributed throughout the regional aquifer and in the local-scale TANC study area at low concentrations (less than 1 microgram per liter). Higher concentrations of PCE were detected in samples from the well-field PSWs and Zarzamora monitoring wells relative to the regional aquifer PSWs. The urban environment is a likely source of contaminants to the aquifer, and these results indicate that one or more local urban sources might be supplying PCE to the Zarzamora monitoring wells and the well-field wells. Samples from the well field also had high concentrations of chloroform relative to the monitoring wells and regional aquifer PSWs. For samples from the regional aquifer PSWs, the most frequently detected organic contaminants generally decreased in concentration with increasing well depth. Deeper wells might intercept longer regional flow paths with higher fractions of older water or water recharged in rural recharge areas in the western part of the aquifer that have been less affected by anthropogenic contaminants. A scenario of hypothetical contaminant loading was evaluated by using results from groundwater-flow-model particle tracking to assess the response of the aquifer to potential contamination. Results indicate that the aquifer responds quickly (less than 1 year to several years) to contaminant loading; however, it takes a relatively long time (decades) for concentrations to reach peak values. The aquifer also responds quickly (less than 1 year to several years) to the removal of contaminant loading; however, it also takes a relatively long time (decades) to reach near background concentrations. Interpretation of geochemical age tracers in this well-mixed karst system was complicated by contamination of a majority of measured tracers and complexities of extensive mixing. Age-tracer results generally indicated that groundwater samples were composed of young, recently recharged water with piston-flow model ages ranging from less than 1 to 41 years, with a median of 17 years. Although a piston-flow model is typically not valid for karst aquifers, the model ages provide a basis for comparing relative ages of different samples and a reference point for more complex hydrogeologic models for apparent-age interpretations. Young groundwater ages are consistent with particle-tracking results from hydrogeologic modeling for the local-scale TANC study area. Age-tracer results compared poorly with other geochemical indicators of groundwater residence time and anthropogenic effects on water quality, indicating that hydrogeologic conceptual models used in groundwater age interpretations might not adequately account for mixing in this karst system. Groundwater samples collected from the well field under a variety of pumping conditions were relatively homogeneous and well mixed for numerous geochemical constituents (with the notable exception of age tracers). Groundwater contributions to the PSW were dominated by well-mixed, relatively homogeneous groundwater, typical of the regional confined aquifer. Zones of preferential flow were determined for the PSW, but groundwater samples from different stratigraphic units were not geochemically distinct. Variations in chemical constituents in response to a rainfall and aquifer recharge event occurred but were relatively minor in the PSW and monitoring wells. This observation is consistent with the hypothesis that the response to individual recharge events in the confined aquifer, unless intersecting conduit flow paths, might be attenuated by mixing processes along regional flow paths. Results of this study are consistent with the existing conceptual understanding of aquifer processes in this karst system and are useful for water-resource development and management practices.

Texas↗

USGS2011

Simulations of groundwater flow and particle-tracking analysis in the zone of contribution to a public-supply well in San Antonio, Texas

In 2006, a public-supply well in San Antonio, Texas, was selected for intensive study to assess the vulnerability of public-supply wells in the Edwards aquifer to contamination by a variety of compounds. A local-scale, steady-state, three-dimensional numerical groundwater-flow model was developed and used in this study to evaluate the movement of water and solutes from recharge areas to the selected public-supply well. Particle tracking was used to compute flow paths and advective traveltimes throughout the model area and to delineate the areas contributing recharge and zone of contribution for the selected public-supply well. The local-scale model grid has a finer vertical discretization than do previous regional Edwards aquifer models and incorporates refined parameter zones corresponding with multiple (10) hydrogeologic units representing the Edwards aquifer. In the Edwards aquifer, high matrix porosity and permeability likely are overshadowed by high permeability developed in structurally influenced karstic conduit systems that transmit water into, through, and out of the aquifer system. The complexity of the aquifer system in the local-scale study area is further increased by numerous faults with varying vertical displacements. The extensive faulting results in the juxtaposition of hydrogeologic units with differing hydraulic properties and has appreciable effects on groundwater flow in the Edwards aquifer. The local-scale model simulations use the MODFLOW Hydrogeologic-Unit Flow Package and include two hydrogeologic units with high hydraulic conductivities (one or more orders of magnitude higher than for the other simulated hydrogeologic units) that are intended to simulate fast flow paths attributable to karst features. The two “conduit” hydrogeologic units of the Edwards aquifer represent the lower 8 meters of the leached and collapsed members and the Kirschberg evaporite member of the Edwards Group. The MODFLOW Horizontal-Flow Barrier Package was used to simulate faults in the local-scale model. The assumption was made that the degree to which a fault acts as a barrier to groundwater flow is proportional to the fault displacement. The final calibrated hydraulic-conductance values ranged from 0.01 to 0.2 per day for fault displacements ranging from 0 to more than 100 percent of the total aquifer thickness. The calibrated steady-state simulation generally reproduces the spatial distribution of measured water-level altitudes. Simulated water-level altitudes were within 9.0 meters of measured water-level altitudes at 74 of the 84 wells used as targets for the local-scale model for the calibrated steady-state simulation. The overall mean absolute difference between simulated and measured water-level altitudes is 4.2 meters, and the mean algebraic difference is 1.9 meters. The simulated springflow for San Antonio Springs was 7.7 percent greater and for San Pedro Springs was 4.2 percent less than the median measured springflow. Simulated tritium concentrations were within 0.14 tritium units of measured tritium concentrations for 11 of the 13 local-scale study tritium observations from the 10 local-scale study wells used to calibrate the steady-state local-scale model, with a mean absolute difference between simulated and measured tritium concentrations of 0.11 tritium units and a mean algebraic difference of -0.04 tritium units. Simulated tritium concentrations in the selected public-supply well during November 2007 were within 0.09 tritium units of the measured concentrations, with the exception of the shallowest observation from the well. The steady-state simulation water budget indicates that recharge occurring in the local-scale study area accounts for 31.8 percent of the sources of water to the Edwards aquifer in the local-scale model area and that inflow through the model boundaries contributes 68.2 percent. Most of the flow into the local-scale model area through the model boundaries occurs through the western and southern boundaries, 58.2 and 39.6 percent, respectively. The largest discharges from the Edwards aquifer in the local-scale model area are boundary outflow (71.4 percent) and withdrawals by wells (24.9 percent). Most of the flow out of the local-scale model area through the model boundaries occurs through the southern and eastern boundaries, 54.2 and 39.6 percent, respectively. The simulated zones of contribution for the selected public-supply well, Timberhill well nest, and Zarzamora well nest extend to the north, northeast, and northwest from each site in the confined zone of the aquifer into the recharge zone, where all recharge to the aquifer occurs. The area contributing recharge for the selected public-supply well has the greatest extent. The area contributing recharge for the Timberhill well nest encompasses approximately the western one-half of the area contributing recharge for the selected public-supply well, and that for the Zarzamora well nest encompasses approximately the eastern two-thirds of the area contributing recharge for the selected public-supply well. Simulated particle ages ranged from less than 1 day to more than 1,900 years in the 10 local-scale study wells (13 local-scale study tritium observations) used to calibrate the local-scale model. The simulated mean particle ages for the tritium observations representing selected well depths (shallow, intermediate, and deep) ranged from 2.5 to 15 years. The minimum (youngest) mean particle ages for the selected public-supply well and the Timberhill monitoring wells were for the intermediate well depth, while the youngest mean particle age for the Zarzamora monitoring wells was for the intermediate and deep well depth. The maximum (oldest) mean particle ages for the selected public-supply well and the Zarzamora monitoring wells were for the shallow well depth. The mean of simulated particle ages for tritium observations representing well depths open to the simulated conduit hydrogeologic units was 3.8 years, whereas the mean of simulated particle ages for tritium observations representing well depths not open to the simulated conduit hydrogeologic units was 9.6 years. The effect of short-circuit pathways, for example karst conduits, in the flow system on the movement of young water to the selected public-supply well could greatly alter contaminant arrival times compared to what might be expected from advection in a system without short circuiting. In a forecasting exercise, the simulated concentrations showed rapid initial response at the beginning and end of chemical input, followed by more gradual response as older water moved through the system. The nature of karst groundwater flow, where flow predominantly occurs via conduit flow paths, could lead to relatively rapid water quality responses to land-use changes. Results from the forecasting exercise indicate that timescales for change in the quality of water from the selected public-supply well could be on the order of a few years to decades for land-use changes that occur over days to decades, which has implications for source-water protection strategies that rely on land-use change to achieve water-quality objectives.

Texas↗

USGS2011

Selected approaches to estimate water-budget components of the High Plains, 1940 through 1949 and 2000 through 2009

The High Plains aquifer, underlying almost 112 million acres in the central United States, is one of the largest aquifers in the Nation. It is the primary water supply for drinking water, irrigation, animal production, and industry in the region. Expansion of irrigated agriculture throughout the past 60 years has helped make the High Plains one of the most productive agricultural regions in the Nation. Extensive withdrawals of groundwater for irrigation have caused water-level declines in many parts of the aquifer and increased concerns about the long-term sustainability of the aquifer. Quantification of water-budget components is a prerequisite for effective water-resources management. Components analyzed as part of this study were precipitation, evapotranspiration, recharge, surface runoff, groundwater discharge to streams, groundwater fluxes to and from adjacent geologic units, irrigation, and groundwater in storage. These components were assessed for 1940 through 1949 (representing conditions prior to substantial groundwater development and referred to as "pregroundwater development" throughout this report) and 2000 through 2009. Because no single method can perfectly quantify the magnitude of any part of a water budget at a regional scale, results from several methods and previously published work were compiled and compared for this study when feasible. Results varied among the several methods applied, as indicated by the range of average annual volumes given for each component listed in the following paragraphs. Precipitation was derived from three sources: the Parameter-Elevation Regressions on Independent Slopes Model, data developed using Next Generation Weather Radar and measured precipitation from weather stations by the Office of Hydrologic Development at the National Weather Service for the Sacramento-Soil Moisture Accounting model, and precipitation measured at weather stations and spatially distributed using an inverse-distance-weighted interpolation method. Precipitation estimates using these sources, as a 10-year average annual total volume for the High Plains, ranged from 192 to 199 million acre-feet (acre-ft) for 1940 through 1949 and from 185 to 199 million acre-ft for 2000 through 2009. Evapotranspiration was obtained from three sources: the National Weather Service Sacramento-Soil Moisture Accounting model, the Simplified-Surface-Energy-Balance model using remotely sensed data, and the Soil-Water-Balance model. Average annual total evapotranspiration estimated using these sources was 148 million acre-ft for 1940 through 1949 and ranged from 154 to 193 million acre-ft for 2000 through 2009. The maximum amount of shallow groundwater lost to evapotranspiration was approximated for areas where the water table was within 5 feet of land surface. The average annual total volume of evapotranspiration from shallow groundwater was 9.0 million acre-ft for 1940 through 1949 and ranged from 9.6 to 12.6 million acre-ft for 2000 through 2009. Recharge was estimated using two soil-water-balance models as well as previously published studies for various locations across the High Plains region. Average annual total recharge ranged from 8.3 to 13.2 million acre-ft for 1940 through 1949 and from 15.9 to 35.0 million acre-ft for 2000 through 2009. Surface runoff and groundwater discharge to streams were determined using discharge records from streamflow-gaging stations near the edges of the High Plains and the Base-Flow Index program. For 1940 through 1949, the average annual net surface runoff leaving the High Plains was 1.9 million acre-ft, and the net loss from the High Plains aquifer by groundwater discharge to streams was 3.1 million acre-ft. For 2000 through 2009, the average annual net surface runoff leaving the High Plains region was 1.3 million acre-ft and the net loss by groundwater discharge to streams was 3.9 million acre-ft. For 2000 through 2009, the average annual total estimated groundwater pumpage volume from two soil-water-balance models ranged from 8.7 to 16.2 million acre-ft. Average annual irrigation application rates for the High Plains ranged from 8.4 to 16.2 inches per year. The USGS Water-Use Program published estimated total annual pumpage from the High Plains aquifer for 2000 and 2005. Those volumes were greater than those estimated from the two soil-water-balance models. Total groundwater in storage in the High Plains aquifer was estimated as 3,173 million acre-ft prior to groundwater development and 2,907 million acre-ft in 2007. The average annual decrease of groundwater in storage between 2000 and 2007 was 10 million acre-ft per year.

High Plains Aquifer↗

USGS2010

Geochemical evolution processes and water-quality observations based on results of the National Water-Quality Assessment Program in the San Antonio segment of the Edwards aquifer, Texas, 1996-2006

As part of the National Water-Quality Assessment Program, the U.S. Geological Survey collected and analyzed groundwater samples during 1996-2006 from the San Antonio segment of the Edwards aquifer of central Texas, a productive karst aquifer developed in Cretaceous-age carbonate rocks. These National Water-Quality Assessment Program studies provide an extensive dataset of groundwater geochemistry and water quality, consisting of 249 groundwater samples collected from 136 sites (wells and springs), including (1) wells completed in the shallow, unconfined, and urbanized part of the aquifer in the vicinity of San Antonio (shallow/urban unconfined category), (2) wells completed in the unconfined (outcrop area) part of the regional aquifer (unconfined category), and (3) wells completed in and springs discharging from the confined part of the regional aquifer (confined category). This report evaluates these data to assess geochemical evolution processes, including local- and regional-scale processes controlling groundwater geochemistry, and to make water-quality observations pertaining to sources and distribution of natural constituents and anthropogenic contaminants, the relation between geochemistry and hydrologic conditions, and groundwater age tracers and travel time. Implications for monitoring water-quality trends in karst are also discussed. Geochemical and isotopic data are useful tracers of recharge, groundwater flow, fluid mixing, and water-rock interaction processes that affect water quality. Sources of dissolved constituents to Edwards aquifer groundwater include dissolution of and geochemical interaction with overlying soils and calcite and dolomite minerals that compose the aquifer. Geochemical tracers such as magnesium to calcium and strontium to calcium ratios and strontium isotope compositions are used to evaluate and constrain progressive fluid-evolution processes. Molar ratios of magnesium to calcium and strontium to calcium in groundwater typically increase along flow paths; results for samples of Edwards aquifer groundwater show an increase from shallow/urban unconfined, to unconfined, to confined groundwater categories. These differences are consistent with longer residence times and greater extents of water-rock interaction controlling fluid compositions as groundwater evolves from shallow unconfined groundwater to deeper confined groundwater. Results for stable isotopes of hydrogen and oxygen indicate specific geochemical processes affect some groundwater samples, including mixing with downdip saline water, mixing with recent recharge associated with tropical cyclonic storms, or mixing with recharge water than has undergone evaporation. The composition of surface water recharging the aquifer, as well as mixing with downdip water from the Trinity aquifer or the saline zone, also might affect water quality. A time-series record (1938-2006) of discharge at Comal Springs, one of the major aquifer discharge points, indicates an upward trend for nitrate and chloride concentrations, which likely reflects anthropogenic activities. A small number of organic contaminants were routinely or frequently detected in Edwards aquifer groundwater samples. These were the pesticides atrazine, its degradate deethylatrazine, and simazine; the drinking-water disinfection byproduct chloroform; and the solvent tetrachloroethene. Detection of these contaminants was most frequent in samples of the shallow/urban unconfined groundwater category and least frequent in samples of the unconfined groundwater category. Results indicate that the shallow/urban unconfined part of the aquifer is most affected by anthropogenic contaminants and the unconfined part of the aquifer is the least affected. The high frequency of detection for these anthropogenic contaminants aquifer-wide and in samples of deep, confined groundwater indicates that the entire aquifer is susceptible to water-quality changes as a result of anthropogenic activities. L

Texas↗

USGS2009

Water-level altitudes 2009 and water-level changes in the Chicot, Evangeline, and Jasper Aquifers and compaction 1973-2008 in the Chicot and Evangeline Aquifers, Houston-Galveston Region, Texas

This report, done in cooperation with the Harris-Galveston Subsidence District, the City of Houston, the Fort Bend Subsidence District, and the Lone Star Groundwater Conservation District, is one in an annual series of reports that depicts water-level altitudes and water-level changes in the Chicot, Evangeline, and Jasper aquifers, and compaction in the Chicot and Evangeline aquifers in the Houston-Galveston region, Texas. The report (excluding appendixes) contains 16 sheets and 15 tables: 3 sheets are maps showing current-year (2009) water-level altitudes for each aquifer, respectively; 3 sheets are maps showing 1-year (2008-09) water-level changes for each aquifer, respectively; 3 sheets are maps showing 5-year (2004-09) water-level changes for each aquifer, respectively; 4 sheets are maps showing long-term (1990-2009 and 1977-2009) water-level changes for the Chicot and Evangeline aquifers, respectively; 1 sheet is a map showing long-term (2000-2009) water-level change for the Jasper aquifer; 1 sheet is a map showing site locations of borehole extensometers; and 1 sheet comprises graphs showing measured compaction of subsurface material at the sites from 1973 or later through 2008, respectively. Tables listing the data used to construct the aquifer-data maps and the compaction graphs are included.

Texas↗