Geology Reports⌕ Search

Geology topics

Research about Texas

Source-linked reports with geographic coverage including Texas.

At least 73 records · Page 4Linked to original sources

Watershed hydrology assessment for the Lower Colorado River Basin. Appendix A: Statistical hydrology

Statistical analysis of the observational record from U.S. Geological Survey (USGS) streamgages and period of historical flow observations prior to the gage installation provides an informative means of estimating flood flow frequency. The U.S. Geological Survey contributed to the InFRM team’s efforts by performing the statistical analysis of the gaged record and authored this Appendix to the Lower Watershed Hydrology Assessment. Flood flow frequency is defined by values or quantiles of streamflow for selected annual exceedance probabilities (AEPs) (England and others, 2019). The annual peak streamflow data collected as part of the systematic operation of a streamgage provides the foundation for a detailed analysis of peak streamflow, but additional historical information pertaining to peak streamflows that predates the installation of a streamgage also can be used. An annual peak streamflow is defined as the maximum instantaneous streamflow for a streamgage for a given water year, and annual peak streamflow data for USGS streamgages can be acquired through the USGS National Water Information System (NWIS) database (USGS, 2022). The statistical analyses are based on water-year increments. A water year is the 12-month period from October 1 of a given year through September 30 of the following year designated by the calendar year in which it ends. For the statistical hydrology portion of a multifaceted analysis, InFRM team members from the USGS analyzed annual peak streamflow records for the 45 USGS streamgages (gages) and 21 Lower Colorado River Authority (LCRA) streamgages (gages) in the lower Colorado River Basin listed in Table A.1 and Table A.8. The locations of USGS gages are also shown on Figure A.1, Figure A.2, and Figure A.3, and the locations of LCRA gages are shown in Section 1.4. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government.

Texas↗

Watershed hydrology assessment for the Lower Colorado River Basin. Appendix D: RiverWare analyses

RiverWare is a river system modeling tool developed by CADSWES (Center of Advanced Decision Support for Water and Environmental Systems) that allows the user to simulate complex reservoir operations and perform period-of-record analyses for different scenarios. For the InFRM hydrology studies, RiverWare is used to generate a homogeneous regulated POR by simulating the basin as if the reservoirs and their current rule sets had been present in the basin for the entire time period. Statistical analyses can then be performed on the extended records at the gages. This report summarizes the RiverWare portion of the hydrologic analysis being completed for the InFRM Hydrology study of the Colorado River Basin. The RiverWare model described in this chapter presents development of the Colorado River Basin hydrology, which mimics current operational conditions. The use of the RiverWare program allows for data extension to periods prior to dam construction. The utilization of longer gage record improves discharge frequency results and increases the confidence of the analysis being performed. The modeling evaluation criteria are: (1) evaluate output based on validating policies and functions, and (2) prioritize operation based on surcharge and flood control. A detailed explanation of the Colorado River Basin POR hydrology will be in a later section. Calibration results will also be shown that illustrate the overall model performance for the POR. The time window simulation run is for January 01, 1930 – September 30, 2019. This time window captures all big events occurred over the Colorado River basin. Each simulated water year was inspected individually to better validate the results. Historical pool elevations along with observed inflows and outflows were compared against the model simulated results.

Texas↗

Geologic framework and hydrostratigraphy of the Edwards and Trinity aquifers within northern Bexar and Comal Counties, Texas

During 2020–22, the U.S. Geological Survey, in cooperation with the Edwards Aquifer Authority, revised a previous publication that described the geologic framework and hydrostratigraphy of the Edwards and Trinity aquifers within northern Bexar and Comal Counties, Texas. This report presents the refined maps and descriptions of geologic framework and hydrostratigraphy of the Edwards and Trinity aquifers within northern Bexar and Comal Counties that resulted from additional field data. Two informal geologic units and their corresponding informal hydrostratigraphic unit (HSU) names are introduced in this report; these informal units were identified during geologic mapping work done in counties adjoining the study area. Hydrostratigraphically, the rocks exposed in the study area represent a section of the upper confining unit to the Edwards aquifer, the Edwards aquifer, the upper zone of the Trinity aquifer, the middle zone of the Trinity aquifer, and the lower confining unit to the middle zone of the Trinity aquifer. The Washita, Eagle Ford, Austin, and Taylor Groups are generally considered to be the upper confining unit to the Edwards aquifer. The Edwards aquifer was subdivided into nine informally named HSUs (from top to bottom) as follows: I, II, III, IV, V, VI, VII, Seco Pass, and VIII. The upper zone of the Trinity aquifer was subdivided into five informal HSUs and two subunits (from top to bottom) as follows: cavernous, Camp Bullis, upper evaporite, fossiliferous (subunits: upper and lower), and lower evaporite. The middle zone of the Trinity aquifer was subdivided into nine named HSUs (from top to bottom) as follows: Bulverde, Little Blanco, Twin Sisters, Doeppenschmidt, Herff Falls (where present), Rust, Honey Creek, Hensell, and Cow Creek. The middle zone of the Trinity aquifer is underlain by the confining Hammett HSU. Groundwater recharge and flow paths in the study area are influenced not only by the hydrostratigraphic characteristics of the individual HSUs but also by faults and fractures.

Texas↗

Bathymetric, hydrodynamic, biological, and water-quality characteristics of a nearshore area of the Laguna Madre near South Padre Island, Texas, 2021–22

A variety of data were collected by the U.S. Geological Survey, in cooperation with the City of South Padre Island, to better understand the physical and biological habitat in Tompkins Channel and adjacent seagrass beds in the lower Laguna Madre, Texas, where the construction of berms has been proposed in the City of South Padre Island’s Shoreline Master Plan. These berms would be used to create living shorelines, defined as shorelines that “connect the land and water to stabilize the shoreline, reduce erosion, and provide ecosystem services, like [maintaining] valuable habitat, that enhances coastal resilience.” Bathymetric surveys show variability in the lagoon depth within the study area and clearly delineate the relatively shallow bay from the deeper Tompkins Channel. There were seasonal patterns in daily mean tidal water-surface elevations, with higher mean tides in the fall and spring and lower mean tides in the winter and summer. Seasonal differences in mean tidal water-surface elevations were primarily caused by astronomical forces. Shorter-term (daily and weekly) variations in tidal water-surface elevations were produced by meteorological forcing resulting from high- and low-pressure systems, the passage of fronts, local storms, and tropical storms. The tide driven currents typically were toward the north-northwest during flood tide or toward the south-southeast during ebb tide, roughly parallel to the general South Padre Island shoreline orientation. Prevailing water current velocity (current speed and direction) indicates that there usually is a net inflow into the lower Laguna Madre from the Gulf of Mexico and a net outflow at locations north of the study area. The overall seagrass density was generally lowest in the northeast section of the transect. Areas of low seagrass density were also most commonly observed close to the channel edge of the transect. Overall, Syringodium filiforme was considerably more abundant than Thalassia testudinum and Halodule wrightii . Amphipods from the families Ampithoidae and Melitidae were the most abundant benthic invertebrates, followed by gastropods in the family Calyptraeidae. Higher abundances of benthic invertebrates were recorded at sites farther away from Tompkins Channel than at sites closer to it, perhaps because of the sensitivity of these organisms to disturbances caused by erosional forces in the channel. Changes in water-quality measured during a 2-day survey were consistent with patterns observed in prior studies, indicating biological processes likely play an important role in the physicochemical environment of the study area. Collectively, the data provide baseline information that can be used as reference conditions prior to any shoreline enhancement activities in the lower Laguna Madre.

Texas↗

Longer holding times decrease dispersal but increase mortality of translocated scaled quail

Scaled quail ( Callipepla squamata ) decline caused by habitat loss and fragmentation increased interest in translocation to reestablish populations. Yet factors determining translocation success are poorly understood. We tested hypotheses concerning the influence of source population and variation in delayed release strategy (1–9 weeks) on mortality and dispersal of wild-caught, translocated scaled quail. We trapped and translocated quail from 2016–2017 from source populations in the Edwards Plateau and Rolling Plains ecoregions to a large contiguous (>40,000 ha) release site in Knox County, Texas, USA. We evaluated mortality and dispersal of translocated females as a function of source population, holding time prior to release, age, release location, and year using a multi-state mark-recapture model with state uncertainty. Scaled quail translocated within the Rolling Plains were more likely to exhibit philopatry to the release site. Quail with longer holding times had higher mortality but lower dispersal rates. The Edwards Plateau is a suitable source site for translocation in the Rolling Plains. The reduced dispersal but higher mortality of translocated scaled quail associated with longer holding times creates a decision tradeoff for managers.

Texas↗

Karst groundwater vulnerability determined by modeled age and residence time tracers

Karst aquifers are a vital groundwater resource globally, but features such as rapid recharge and conduit flow make them highly vulnerable to land-surface contamination. We apply environmental age tracers to the south-central Texas Edwards aquifer, a karst resource in a rapidly urbanizing and drought-prone region, to assess vulnerability to land-surface contamination and risks unique to karst aquifers. We show that vulnerability of Edwards aquifer groundwater follows similar spatial and depth patterns common to porous-media type aquifers, despite complicated karst hydrogeologic features. Shallow and unconfined parts are more vulnerable to land-surface contamination than the deeper and confined parts, although even the oldest groundwater is mixed with some recent recharge. When modeled age-tracer results are coupled with other independent geochemical tracers of water-rock interaction specific to karst settings, they can yield insight into residence time and associated vulnerability.

Texas↗

Discrete streamflow measurements and waterborne self-potential logging of a 43-kilometer-long reach of the Elm Fork Trinity River upstream from Dallas, Texas

Continuous and discrete streamflow data were combined with waterborne self-potential (WaSP), surface-water temperature and surface-water conductivity surveys obtained along an approximately 43-kilometer (26.7 mile) surveyed reach of the Elm Fork Trinity River (hereinafter referred to as “Elm Fork”) upstream from Dallas, Texas, to investigate areas of gaining and losing streamflow under various streamflow and seasonal climatic conditions. Discrete streamflow measurements were made at 17 locations on October 12, 2021, and January 25, 2022, at 19 locations on May 17, 2022, and at 18 locations on August 9, 2022. WaSP data were measured from a kayak in January 2022 during a period of base flow along three individually surveyed reaches between the Lake Lewisville Dam and Frasier Dam on the Elm Fork. Together, these data indicated different parts of the Elm Fork functioned as either a gaining or losing stream depending on streamflow and seasonal climatic conditions. Overall, there were estimated net gains in streamflow during the first two discrete-measurement events of about 107 cubic feet per second (ft 3 /s) and 2 ft 3 /s in October 2021 and January 2022, respectively, and estimated net losses in streamflow in May 2022 and August 2022 of about 24 ft 3 /s and 18 ft 3 /s, respectively.

Texas↗

Occurrence of pharmaceutical compounds in the San Antonio segment of the Edwards (Balcones fault zone) aquifer, south-central Texas, June 2018–August 2020

The occurrence of pharmaceutical compounds (hereinafter referred to as “pharmaceuticals”) in surface water that recharges the San Antonio segment of the Edwards aquifer in south-central Texas is of concern with respect to potential effects on groundwater quality. This study, conducted during June 2018–August 2020 by the U.S. Geological Survey in cooperation with the San Antonio Water System, used a source-to-sink approach (that is, from aquifer recharge to aquifer discharge) to compare the occurrences and concentrations of pharmaceuticals in recharging (losing) streams to those in relatively shallow and relatively deep groundwater. Individual pharmaceutical concentrations in surface water were all less than 200 nanograms per liter, and most were less than 50 nanograms per liter. The two most common pharmaceuticals in surface water were metformin (50-percent detection frequency) and caffeine (25-percent detection frequency). In contrast to surface-water sites, few detections of pharmaceuticals above the detection limit were quantified at groundwater sites. On the basis of studied sample sites that are considered to be representative of the Edwards aquifer recharge zone, groundwater results collectively indicate that pharmaceuticals currently (2020) do not substantially impair water quality in the shallow unconfined part of the Edwards aquifer. Although groundwater pharmaceutical detections were not common, two detections of acetaminophen occurred in samples collected from a relatively deep confined part of the aquifer, indicating that this part of the aquifer is potentially vulnerable to pharmaceutical contamination. The near absence of wastewater treatment plant discharges and Texas Land Application Permit facilities within the drainage areas of streams in this study could explain reduced inputs of pharmaceuticals relative to pesticides on the Edwards aquifer recharge zone and their reduced occurrence in unconfined groundwater. In the western part of the study area, numerous pharmaceutical detections in the Frio River indicate that pharmaceutical sources in rural areas also exist with potential water-quality effects. Although reported pharmaceutical concentrations for the San Antonio segment of the Edwards aquifer are currently (2020) low, ongoing development in the region has the potential to increase contaminant loads, including for pharmaceuticals.

Texas↗

Seed banks of rare Physostegia correllii (Lamiaceae) in Lady Bird Lake, Austin, Texas, U.S.A.

Rare species threatened by climate and land-use change may harbor seeds in soil seed banks for periods of time even if adults have disappeared from the site. Soil samples were collected from sites with current Phyostegia correllii populations and from sites with former populations in Lady Bird Lake (a reservoir of the Colorado River, Austin, Texas. A seedling emergence study was conducted under greenhouse conditions, and the presence/absence of seedling emergence was recorded for two years. Seeds germinated from the seed banks of all current and former colonies tested. The presence of seed banks in a historical site (Blunn Creek) of Physostegia correllii suggests that management to encourage the germination of seeds might help to encourage the establishment of populations of this species. The re-establishment of disturbance fugitives might be facilitated by removing overhanging ground vegetation or imposing water management regimes that mimic natural floodplain dynamics.

Texas↗

Mapping abandoned uranium mine features using Worldview-3 imagery in portions of Karnes, Atascosa and Live Oak Counties, Texas

Worldview-3 (WV3) 16-band multispectral data were used to map exposed bedrock and mine waste piles associated with legacy open-pit mining of sandstone-hosted roll-front uranium deposits along the South Texas Coastal Plain. We used the “spectral hourglass” approach to extract spectral endmembers representative of these features from the image. This approach first requires calibrating the imagery to reflectance, then masking for vegetation, followed by spatial and spectral data reduction using a principal component analysis-based procedure that reduces noise and identifies homogeneous targets which are “pure” enough to be considered spectral endmembers. In this case, we used a single WV3 image which covered an ~11.5 km by ~19.5 km area of Karnes, Atascosa and Live Oak Counties, underlain by mined rocks from the Jackson Group and Catahoula Formation. Up to 58 spectral endmembers were identified using a further multi-dimensional class segregation method and were used as inputs for spectral angle mapper (SAM) classification. SAM classification resulted in the identification of at least 117 mine- and mine waste-related features, most of which were previously unknown. Class similarity was further evaluated, and the dominant minerals in each class were identified by comparison to spectral libraries and measured samples of actual Jackson Group uranium host rocks. Redundant classes were eliminated, and SAM was run a second time using a reduced set of 23 endmembers, which were found to map these same features as effectively as using the full 58 set of endmembers, but with significantly reduced noise and spectral outliers. Our classification results were validated by evaluating detailed scale mapping of three known mine sites (Esse-Spoonamore, Wright-McCrady and Garbysch-Thane) with published ground truth information about the vegetation cover, extent of erosion and exposure of waste pile materials and/or geologic information about host lithology and mineralization. Despite successful demonstration of the utility of WV3 data for inventorying mine features, additional landscape features such as bare agricultural fields and oil and gas drill pads were also identified. The elimination of such features will require combining the spectral classification maps presented in this study with high-quality topographic data. Also, the spectral endmembers identified during the course of this study could be useful for larger-scale mapping efforts using additional well-calibrated WV3 images beyond the coverage of our initial study area.

Texas↗

The use of historical data and global climate models to assess historical and future surface water and groundwater availability in the Trinity River Basin in Texas

This paper describes the results of a study that was done by the USGS to assess recent (2017) water availability, forecast long-term trends in water availability, assess changes in water availability, and forecast future water availability in the Trinity River Basin in Texas. The Trinity River Basin surface water model and Trinity River alluvium aquifer (TRAA) groundwater model were created to evaluate future conditions under different global climate models (GCM). The results of this study show minimal overall changes in water availability for both surface water and groundwater. Trend analyses using historical data (1900–2017) indicated an increase of annual precipitation on the watersheds that drain into the reservoirs in Regional Water Planning Group C. However, the Trinity River Basin surface water model GCM ensemble mean annual precipitation indicates a downward trend, resulting in a downward trend in surface runoff. Additionally, the GCM ensemble mean for the Trinity River Basin surface water model and the TRAA groundwater model both indicate a downward trend in recharge while the TRAA model GCM ensemble mean indicates an upward trend in the amount of groundwater leaving the aquifer to rivers and streams resulting in an upward trend of cumulative storage change.

Texas↗

Status of water-level altitudes and long-term and short-term water-level changes in the Chicot and Evangeline (undifferentiated) and Jasper aquifers, greater Houston area, Texas, 2022

Since the early 1900s, groundwater withdrawn from the primary aquifers that compose the Gulf Coast aquifer system—the Chicot, Evangeline, and Jasper aquifers—has been an important source of water in the greater Houston area, Texas. 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 the status of water-level altitudes and water-level changes in these aquifers in the greater Houston area. In this report, the Chicot and Evangeline aquifers are treated as a single hydrogeologic unit for the purposes of providing annual assessments of regional-scale water-level altitudes and changes over time. In 2022, shaded depictions of water-level altitudes for the Chicot and Evangeline aquifers (undifferentiated) ranged from about 270 feet (ft) below the North American Vertical Datum of 1988 (NAVD 88) to about 195 ft above NAVD 88. The largest decline in water-level altitudes indicated by the 1977–2022 long-term water-level-change map was southeast of The Woodlands. In comparison, the 1990–2022 long-term water-level-change map depicts declines in water-level altitudes in localized areas at or near certain wells in parts of northwestern Harris County and southern Montgomery County. The largest rise in water-level altitudes for 1977–2022 was observed in a relatively large area in southeastern Harris County, whereas the largest rise in water-level altitudes for 1990–2022 was observed in a relatively small area in central Harris County. The 5-year short-term water-level-change map depicts the largest declines at three wells in northern Fort Bend County and the largest rises at five wells in southwestern and central Harris County. The 1-year short-term water-level-change map depicts the largest decline at a well in north-central Harris County and the largest rises at a well in north-central Fort Bend County and a well in southwestern Harris County. In 2022, shaded depictions of water-level altitudes for the Jasper aquifer ranged from about 213 ft below NAVD 88 to about 285 ft above NAVD 88. The 2000–22 long-term water-level-change map depicts water-level declines throughout most of the study area where water-level-measurement data from the aquifer were collected, with the largest decline in north-central Harris County and south-central Montgomery County south of The Woodlands. The 5-year short-term water-level-change map depicts the largest decline at a well in Conroe and the largest rise at a well in south-central Montgomery County. The 1-year short-term water-level-change map depicts the largest decline at a well in south-central Montgomery County east of The Woodlands and the largest rise at a well in Conroe.

Texas↗

Flood warning toolset for the Sabinal River near Utopia, Texas

Introduction Floods are one of the most frequent and expensive natural disasters that occur across the United States. Rapid, high-water events that occur in local areas—flash floods—are especially difficult for emergency managers to predict and provide advance warning to the public, and insufficient data can hamper postflood recovery efforts. Central Texas is hilly, and it is known as a “flash flood alley” because of its high-intensity rains, shallow soils, and steep terrain, all of which combined can result in loss of life and property damage. For example, the flash flood event during July 2002 claimed 12 lives in central Texas, including 1 in the town of Utopia, which is on the east bank of the Sabinal River in a flash-flood-prone area along the Balcones Escarpment. During the flood event, the peak discharge recorded on July 5, 2002, at U.S. Geological Survey (USGS) streamgage 08198000 Sabinal River near Sabinal, Tex. (hereinafter referred to as the “Sabinal gage”), was 108,000 cubic feet per second (corresponding to a stream stage [also called gage height] of 33.74 feet). To put the 2002 flood into context, during a typical year the median daily discharge in the Sabinal River at the Sabinal gage is only about 23 cubic feet per second. In 2021, the USGS, in cooperation with the Bandera County River Authority and Groundwater District and the Texas Water Development Board, developed a flood warning toolset for the Sabinal River near Utopia. This study builds on earlier USGS flood work on the Medina River in Bandera County. The newly developed toolset consists of a newly installed USGS streamgage to collect continuous stream stage data (streamgage 08197970 Sabinal River at Utopia, Tex.; hereinafter referred to as the “Utopia gage”) 13 miles upstream from the Sabinal gage, a hydraulic model developed for the Sabinal River near Utopia, and an online library of digital flood-inundation maps referenced to the stream stage at the Utopia gage.

Texas↗

Flood-inundation maps created using a synthetic rating curve for a 10-mile reach of the Sabinal River and a 7-mile reach of the West Sabinal River near Utopia, Texas, 2021

In 2021, the U.S. Geological Survey (USGS), in cooperation with the Bandera County River Authority and Groundwater District and the Texas Water Development Board, studied floods to produce a library of flood-inundation maps for the Sabinal River near Utopia, Texas. Digital flood-inundation maps were created for a 10-mile reach of the Sabinal River from USGS streamgage 08197936 Sabinal River below Mill Creek near Vanderpool, Tex., at the upstream boundary of the study reach, to USGS streamgage 08197970 Sabinal River at Utopia, Tex. (hereinafter referred to as the “Utopia gage”), at the downstream boundary of the study reach, and for a 7-mile reach of the West Sabinal River. The flood-inundation maps depict estimates of the areal extent and depth of flooding corresponding to selected gage heights (the water-surface elevation at a streamgage, commonly referred to as “stage”) at the Utopia gage. Water-surface elevations were computed for the stream reach by means of a two-dimensional unsteady-state diffusion wave model with the U.S. Army Corps of Engineers Hydrologic Engineering Center River Analysis System program. A synthetic stage-discharge rating curve at the Utopia gage was developed using a regional regression equation to construct the model boundary condition inputs, and the upper bound of the stage-discharge relation was matched to a major flood event in July 2002. The hydraulic model was used to compute water-surface elevations for 35 stages at 0.5-foot (ft) increments referenced to the Utopia gage datum and ranging from 11 ft (near bankfull) to 28 ft (estimated peak stage during the July 2002 flood event). These flood-inundation maps, in conjunction with the real-time stage data from the Utopia gage, are intended to help guide the public in taking individual safety precautions and provide emergency management personnel with a tool to efficiently manage emergency flood operations and postflood recovery efforts.

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↗

Distribution of rare earth and other critical elements in lignites from the Eocene Jackson Group, Texas

Coal is increasingly evaluated as a source of rare earth elements (REEs) in the United States to address the overreliance on imported REEs. The objective of this study was to assess the distribution of REEs in lignites from selected mining areas in the Texas Gulf Coastal Plain region. Thirty-one archived lignite and rock samples previously collected by the U.S. Geological Survey were analyzed for their rare earth element and critical mineral content. These include samples from one core (5400 and 5500 lignite horizons) and two opencast lignite mines (Gibbons Creek 3500 and 4500 horizons, and San Miguel horizons A to D) in the Eocene Jackson Group of the Texas Gulf of Mexico Coastal Plain. Some lithologies in the Gibbons Creek 3500 and 4500 lignite-bearing sections have high total rare earth, yttrium (Y), and scandium (Sc) (REYSc) values, up to 7800 ppm (ash basis) REYSc. The lignite lithologies show an enrichment in rare earths, [samarium (Sm) through gadolinium (Gd)]. The basal Gibbons Creek 3500 lignite bench shows a heavy rare earth element enrichment pattern resembling that often seen in peats through high volatile A bituminous coals. The 5500 lignite sequence, overlying the latter lignite sections, shows a light rare earth enrichment. The San Miguel lignite benches have heavy rare earth enrichments with a negative europium (Eu) anomaly.

Texas↗

Can spatial food web subsidies associated with river hydrology and lateral connectivity be detected using stable isotopes?

During and following lateral connections, aquatic organisms residing in the river channel may assimilate material from sources imported from oxbows, and oxbow residents may consume and assimilate material imported from the channel. Hydrology, lateral connectivity, and stable isotope ratios of fishes and mussels were analyzed for evidence of spatial food web subsidies between the active channel and oxbow lakes in the floodplain of the Guadalupe River, Texas. During surveys conducted between March 2016 and April 2017, fish, mussel, periphyton, seston, and riparian plant samples were collected in and around two oxbows and adjacent channel sites for analysis of stable isotope ratios. Biplots of δ 13 C and δ 15 N were graphed for basal sources and specimens of six common fish species, four sunfish species ( Lepomis spp. combined), and two mussel species (Unionidae combined) captured from oxbows and the channel. Within each graph, polygons were drawn to indicate the space occupied by animals that could have assimilated feasible combinations of source materials originating from either oxbows or the river channel. Based on positions of animals within source polygons, riparian C4 grasses were not an important source of organic matter supporting biomass of fishes and mussels within the channel or oxbows. Overall, 84% of organisms had isotopic signatures consistent with assimilation of in situ sources, but also 76% of all organisms were inconclusive with regards to cross-habitat exchanges. Outliers that may have assimilated ex situ source material were observed for only 4% of 313 organisms from oxbows and 9% of 232 organisms from the channel, and some but not all of these cases followed high flow pulses that connected oxbows for extended periods. Several issues that compromise inferences from stable isotope analysis were identified, and estimation of spatial food web subsidies in fluvial systems could be enhanced by analyzing additional biomarkers, such as isotopic ratios of other elements and compound-specific stable isotopes, as well as additional sources, time-specific biotracers, and experimental approaches that directly track movement of sources and organisms in spatially structured food webs.

Texas↗

A residual oil zone (ROZ) assessment methodology with application to the central basin platform (Permian Basin, USA) for enhanced oil recovery (EOR) and long-term geologic CO2 storage

Residual oil zones (ROZ) form due to various geologic conditions and are located below the oil/water contact (OWC) of main pay zones (MPZ). Since ROZs usually contain immobile oil, they have not typically been considered commercially attractive for development by conventional primary recovery methods used in the initial phases of oil production. However, during the last decade some operators of these viable fields that also contain ROZs have extended carbon dioxide enhanced oil recovery (CO 2 -EOR) to below the MPZ to commercially recover oil from the associated ROZ. Increased interest in ROZs is also due to the application of anthropogenic CO 2 for oil recovery, leading to the subsurface sequestration of CO 2 , which can be part of the current net-zero carbon oil and climate change objectives. Several detailed studies of selected formations in the Permian Basin of the United States have shown that ROZs can be as common as traditional conventional oil reservoir traps, suggesting significant resources for potential additional hydrocarbon recovery and subsurface CO 2 sequestration via CO 2 -EOR. However, applications of CO 2 -EOR to ROZs have been limited despite the estimation of significant oil resources considered recoverable through CO 2 -EOR, and the benefit of concurrent geologic CO 2 storage that would help offset carbon emissions from the produced oil. The combination of insufficient economic incentives and technical reasons related to data scarcity, such as lack of penetration of wells and well logs, for locating ROZs has limited development of their resource potential when compared to known fields. This paper presents a probabilistic methodology for identifying and evaluating ROZ resources for CO 2 -EOR and CO 2 sequestration potential with the use of public and proprietary data sources. The methodology was developed during a pilot study that focused on the ROZ in the San Andres Formation of a nine-county area in the Permian Basin in West Texas. The pilot study estimated a mean oil in place of 25 × 10 9 barrels (bbl) of oil and a mean potential incremental oil recovery and CO 2 utilization of 2.6 × 10 9 bbl and 28.2 Tcf (1.46 × 10 9 tons), respectively, with 1 hydrocarbon pore volume (HCPV) of injection using the water alternating gas (WAG) method. The results of this pilot study are consistent with reported volumes in the literature for a similar area in the Permian Basin. The pilot study demonstrated that this methodology could be used to identify and assess the recoverable oil and coincident CO 2 storage volumes of ROZs in other formations and regions.

Texas↗