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...

In south-central Texas the lower Guadalupe River has incised into the outcrop of the Carrizo-Wilcox aquifer. The river and the aquifer are hydraulically connected across the outcrop, although the connectivity is obscured at the surface by alluvium and surface-water and groundwater exchange dynamics are currently poorly understood. To investigate surface-water and groundwater exchange dynamics between the lower Guadalupe River and the Carrizo-Wilcox aquifer, a geophysical study was completed along a 14.86 km reach of the river by using water-borne gradient self-potential (SP) profiling and two-dimensional direct-current electric resistivity tomography. This paper explores the applicability of these water-borne geoelectric methods in delineating gaining and losing channel reaches, and demonstrates that geoelectric signals in the form of total electric field strength can be logged with an electric dipole and decomposed into component SP signals depicting regional and local groundwater flow patterns attributable to regional and localized hydraulic gradients. Localized SP anomalies of several tens of millivolts, indicative of hyporheic exchange flows, are observed and superimposed upon a 124 mV regional SP anomaly indicative of ambient groundwater exchange flows between the river and the aquifer. The observed SP signals are interpreted through two-dimensional finite-element modeling of streaming potentials attributable to ambient groundwater exchange and hyporheic exchange flow patterns. Variables of the channel environment such as temperature and concentration gradients, depth, and velocity are considered and subsequently eliminated as alternative sources of the SP signals that are presented.

One of the largest rechargeable groundwater systems by total available volume in the Rio Grande/Río Bravo Basin (hereinafter referred to as the “Rio Grande”) region of the United States and Mexico, the Mesilla Basin/Conejos-Médanos aquifer system, supplies water for irrigation as well as for cities of El Paso, Texas; Las Cruces, New Mexico; and Ciudad Juárez, Chihuahua, Mexico. The U.S. Geological Survey in cooperation with the Bureau of Reclamation assessed the groundwater resources in the Mesilla Basin and surrounding areas in Doña Ana County, N. Mex., and El Paso County, Tex., by using a combination of geophysical and geochemical methods. The study area consists of approximately 1,400 square miles in Doña Ana County, N. Mex., and 100 square miles in El Paso County, Tex. The Mesilla Basin composes most of the study area and can be divided into three parts: the Mesilla Valley, the West Mesa, and the East Bench. The Mesilla Valley is the part of the Mesilla Basin that was incised by the Rio Grande between Selden Canyon to the north and by a narrow valley (about 4 miles wide) to the southeast near El Paso, Tex., named the Paso del Norte, which is sometimes referred to in the literature as the “El Paso Narrows.”Previously published geophysical data for the study area were compiled and these data were augmented by collecting additional geophysical and geochemical data. Geophysical resistivity measurements from previously published helicopter frequency domain electromagnetic data, previously published direct-current resistivity soundings, and newly collected (2012) time-domain electromagnetic soundings were used in the study to detect spatial changes in the electrical properties of the subsurface, which reflect changes that occur within the hydrogeology. The geochemistry of the groundwater system was evaluated by analyzing groundwater samples collected in November 2010 for physicochemical properties, major ions, trace elements, nutrients, pesticides (reported but not used in the assessment), and environmental tracers. The data obtained from these samples (with the exception of the pesticide data) were used to gain insights into processes controlling the groundwater movement through the groundwater system in the study area. Results from the geophysical and geochemical assessments facilitated the interpretation of the geochemical characteristics of the groundwater sources and geochemical groups within the groundwater system.The groundwater-flow system in the study area consists primarily of the Mesilla Basin aquifer system, which can be divided into four hydrogeologic units by using an informal classification scheme based on basin-fill stratigraphy and sedimentology with an emphasis on aquifer characteristics. The four hydrogeologic units are (1) the Rio Grande alluvium, which is the shallow aquifer of the Mesilla Basin within the confines of the Mesilla Valley, and the three hydrogeologic units that compose the Santa Fe Group: (2) the lower part of the Santa Fe Group, which is the least productive zone, (3) the middle part of the Santa Fe Group, which is the primary water-bearing hydrogeologic unit in the basin and is generally saturated, and (4) the upper part of the Santa Fe Group, which is the most productive water-bearing unit within the Santa Fe Group but is only partially saturated in the north and largely unsaturated in the south and western parts of the Mesilla Basin.The helicopter frequency domain electromagnetic survey results indicated that approximately half of the resistivity values were less than 10 ohm-meters at depths of 50 and 100 feet with a transition where the resistivity values changed from relatively high values (greater than 20 ohm-meters) to relatively low resistivity values (less than 10 ohm-meters) near Vado, New Mexico. Slightly more than 25 percent of the gridded resistivity values from the three-dimensional grid of the combined inverse modeling results of the direct-current resistivity and time-domain electromagnetic soundings were equal to or less than 10 ohm-meters with large regions of low resistivity becoming apparent in the southernmost part of the study area near the Paso Del Norte where these low resistivity features are spatially the widest at or below the top of the bedrock. These low resistivity values might represent clayey deposits, sediments composed largely of sand and gravel saturated with saline water, or both. Historical dissolved-solids-concentration data within the surface geophysical subset area of the study area were compiled and compared to the inverse modeling results of the combined direct-current resistivity and time-domain soundings; this comparison was done to strengthen the interpretation made from the combined inverse modeling results that the low resistivity features were representative of sand and gravel deposits saturated with saline water and not clayey deposits.Water-level altitudes within the Rio Grande alluvium generally decreased from north to south, with a west to east decrease in water-level altitudes near Las Cruces, New Mexico, as a result of groundwater pumping. Groundwater flow within the Santa Fe Group is more complex than the groundwater flow within the Rio Grande alluvium because of the larger lateral and vertical extent of the Santa Fe Group compared to the Rio Grande alluvium. Groundwater from the Organ Mountains flows directly south towards the Paso del Norte. Groundwater from the Robledo Mountains, the Rough and Ready Hills, and the Sleeping Lady Hills generally flows to the southeast. Groundwater flowing near the north end of the midbasin uplift generally continues east towards the Rio Grande and then flows south on the east side of the midbasin uplift. Groundwater flowing near the west side of the midbasin uplift generally continues south parallel to the faults that make up the midbasin uplift and then flows east towards the Paso del Norte when it reaches the...

Synthesizing binational data to characterize shared water resources is critical to informing binational management. This work uses binational hydrogeology and water resource data in the Mesilla/Conejos-Médanos Basin (Basin) to describe the hydrologic conceptual model and identify potential research that could help inform sustainable management. The Basin aquifer is primarily composed of continuous basin-fill Santa Fe Group sediments, allowing for transboundary throughflow. Groundwater flow, however, may be partially or fully restricted by intrabasin uplifts and limited recharge. The shallow groundwater in the Rio Grande alluvium receives recharge from the Rio Grande and responds to changes in water supply and demand. About 11% of Rio Grande alluvial groundwater volume is recharged annually, an amount that is less than recent withdrawals. Potentially recoverable fresh to slightly brackish groundwater was estimated at 82,600 cubic hectometers in the U.S. portion of the Basin and 69,100 cubic hectometers in the Mexican portion. Alluvial groundwater geochemistry is governed by the evaporative concentration of the Rio Grande and agricultural diversions, whereas deeper groundwater geochemistry is governed by mixing and geochemical processes. Continued refinements to storage estimates, the water budget, and deep groundwater extent and geochemistry can improve estimates of sustainable use and inform alternative water sources.

This study presents a geochemical framework and geographic information system (GIS) method for assessing the intrinsic potential of surface water and groundwater to mobilize arsenic, molybdenum, selenium, uranium, and vanadium. The method was created using published groundwater and surface water geochemical data from the National Uranium Resource Evaluation database for 2302 groundwater and 915 surface water samples. The method was evaluated using published groundwater geochemical data from the Texas Water Development Board. Geochemical data were analyzed in GIS. Samples were categorized by environmental condition, which was determined by using reduction–oxidation—as indicated by pe—and pH ranges for each sample based on geochemical mobility frameworks developed by Smith (2007) and Perel’man (1986). Reduction–oxidation and pH influence the occurrence, persistence, and mobility of arsenic, molybdenum, selenium, uranium, and vanadium in groundwater and surface water. Reduction–oxidation categories were assigned to water samples using concentrations of redox-active constituents, including dissolved oxygen, iron, manganese, and sulfur. The presence of iron substrates and hydrogen sulfides were considered in relation to mobility mechanisms. Twelve-digit hydrologic unit code (HUC) boundaries were used in GIS as analysis areas to determine the most commonly occurring environmental condition in each HUC. The resulting maps identify the environmental conditions in different areas that can be used to identify where the elements are mobile. This methodology provides a systematic approach to identify areas where elements in groundwater and surface water may occur and persist and may be transferable to other locations.

The Rio Grande/Río Bravo del Norte (hereinafter referred to as the “Rio Grande”) is the primary source of recharge to the Mesilla Basin/Conejos-Médanos aquifer system in the Mesilla Valley of New Mexico and Texas. The Mesilla Basin aquifer system is the U.S. part of the Mesilla Basin/Conejos-Médanos aquifer system and is the primary source of water supply to several communities along the United States–Mexico border in and near the Mesilla Valley. Identifying the gaining and losing reaches of the Rio Grande in the Mesilla Valley is therefore critical for managing the quality and quantity of surface and groundwater resources available to stakeholders in the Mesilla Valley and downstream. A gradient self-potential (SP) logging survey was completed in the Rio Grande across the Mesilla Valley between 26 June and 2 July 2020, to identify reaches where surface-water gains and losses were occurring by interpreting an estimate of the streaming-potential component of the electrostatic field in the river, measured during bankfull flow. The survey, completed as part of the Transboundary Aquifer Assessment Program, began at Leasburg Dam in New Mexico near the northern terminus of the Mesilla Valley and ended ~72 kilometers (km) downstream at Canutillo, Texas. Electric potential data indicated a net losing condition for ~32 km between the Leasburg Dam and Mesilla Diversion Dam in New Mexico, with one ~200-m long reach showing an isolated saline-groundwater gaining condition. Downstream from the Mesilla Diversion Dam, electric-potential data indicated a neutral-to-mild gaining condition for 12 km that transitioned to a mild-to-moderate gaining condition between 12 and ~22 km downstream from the dam, before transitioning back to a losing condition along the remaining 18 km of the survey reach. The interpreted gaining and losing reaches are substantiated by potentiometric surface mapping completed in hydrostratigraphic units of the Mesilla Basin aquifer system between 2010 and 2011, and corroborated by surface-water temperature and conductivity logging and relative median streamflow gains and losses, quantified from streamflow measurements made annually at 16 seepage-measurement stations along the survey reach between 1988 and 1998 and between 2004 and 2013. The gaining and losing reaches of the Rio Grande in the Mesilla Valley, interpreted from electric potential data, compare well with relative median streamflow gains and losses along the 72-km long survey reach.