Geology topics

Research about Paradox Valley

Source-linked reports with geographic coverage including Paradox Valley.

9 recordsLinked to original sources

USGS2024

Hydrogeologic conceptual model of groundwater occurrence and brine discharge to the Dolores River in the Paradox Valley, Montrose County, Colorado

Salinity, or total dissolved solids (TDS), of the Colorado River is a major concern in the southwestern United States where the river provides water to about 40 million people for municipal and industrial use and is used to irrigate about 5.5 million acres of land. Much of the salinity in the Colorado River Basin is derived from natural interactions of surface water and groundwater with various geologic materials (rocks, soils, and alluvial deposits). The Dolores River in southwest Colorado is a major tributary of the Colorado River that historically accounts for about 6 percent of the salinity load to the Upper Colorado River Basin with the Paradox Valley being the primary source of salinity to the Dolores River. The Paradox Valley, one of several salt-anticline valleys in the region, is a fault-bounded topographic basin aligned with and exposing an underlying salt-anticline core. Salt deposits in the Pennsylvanian Paradox Formation of the Hermosa Group form an elongated salt diapir oriented northwest to southeast that is up to 12,000 feet (ft) thick beneath the present valley floor. Surface erosion, groundwater circulation, and weathering during Tertiary and Quaternary valley formation contributed to development of a cap rock, collapse features, breccia, and brine at the top of the exposed salt diapir. Today (2023), brine occurring in the brecciated cap rock and underlying salt deposits is in hydraulic connection with an overlying freshwater alluvial aquifer, and depending on seasonal river stage and hydrologic conditions, the brine discharges to the Dolores River causing the observed increase in salinity as the river crosses the Paradox Valley. To reduce salinity concentrations in the Dolores River, the Bureau of Reclamation (Reclamation) operates the Paradox Valley Unit (PVU). The PVU project consists of nine shallow brine pumping wells near the Dolores River and one deep disposal well where the brine is injected for disposal. When operational, the PVU pumping wells extract brine from the base of the alluvial aquifer that is piped and injected into a deep disposal well about 3 miles southwest of the PVU. The PVU became fully operational July 1, 1996, and by 2015, operation of the PVU had reduced salinity concentrations in the Dolores River by as much as 70 percent compared to pre-PVU conditions. In response to a 4.5 magnitude earthquake, injection operations, and thus PVU pumping, were ceased from March 2019 to June 2022. A trial period of PVU operation began in June 2022 with a reduced injection rate, and thus PVU pumping rate, of about two-thirds capacity to gather additional information and guide future operational decisions. In cooperation with Reclamation, the U.S. Geological Survey (USGS) developed this report to present the current (2023) understanding of groundwater and brine occurrence and discharge to the Dolores River in the Paradox Valley. Results from the compilation of spatial datasets, groundwater sampling and age dating, and aquifer tests are presented to provide improved understanding of the Paradox Valley hydrogeology, to supply datasets for a numerical groundwater-flow and brine-transport model, and to support future operations of the PVU. The hydrogeologic data provided herein, along with the most recent loading analysis for the Dolores River in the Paradox Valley, and a previous conceptual model for brine discharge to the river are used to present a conceptual understanding of groundwater occurrence in the Paradox Valley.

Colorado↗

USGS2024

Simulation of groundwater flow and brine discharge to the Dolores River in the Paradox Valley, Montrose County, Colorado

Salinity, or total dissolved solids (TDS), of the Colorado River affects agricultural, municipal, and industrial water users and is an important concern in the Western United States. In the Paradox Valley of southwestern Colorado, natural discharge of sodium-chloride brine to the Dolores River from the underlying core of a salt-valley anticline accounts for about 6 percent of the salinity load to the Colorado River. Formation of the Paradox Valley began during the Miocene, and subsequent erosion exposed the Pennsylvania Paradox Formation in the core of the anticline where a cap rock, collapse features, breccia, and sodium-chloride saturated brine developed at the top of the exposed salt diapir. The discharge of brine to the Dolores River is affected by these dissolution features, along with seasonal hydrologic conditions and density-dependent flow between older dense brine and the younger fresh groundwater in the overlying alluvial aquifer. To reduce TDS concentrations in the Dolores River through the Paradox Valley, the Bureau of Reclamation has pumped brine from a series of shallow wells adjacent to the river since July 1996. The pumped brine is collected and piped to a deep disposal well where it is injected into the Mississippian Leadville Limestone at a depth of about 4,570-meters below land surface. The pumping and injection operation is collectively known as the Paradox Valley Unit (PVU), and by 2015, the PVU had substantially reduced TDS concentrations in the Dolores River by about 70 percent. Since 2019, injection-pressure limits and related seismic activity have constrained deep-well injection and thus brine pumping at the PVU. In cooperation with the Bureau of Reclamation, the U.S. Geological Survey developed a MODFLOW-6 three-dimensional, variable-density groundwater flow and TDS transport model of the Paradox Valley to evaluate the effects of PVU pumping operations on brine discharge to the Dolores River and to guide additional research. The finite-difference model grid consists of 76 rows and 48 columns oriented from northwest to southeast in alignment with valley topography and groundwater-flow directions in the near-surface freshwater alluvial aquifer. A 7-layer hydrogeologic framework was developed from existing datasets to represent the alluvial aquifer, cap rock, collapse breccia, and groundwater flow and TDS transport from the underlying Paradox Formation salt to the Dolores River. The model represents a 33-year transient calibration period from 1987 through 2020 that includes pre-PVU conditions from 1987 through June 1996 and post-PVU conditions from July 1996 through 2020. A 1,000-year simulation of groundwater flow and coupled TDS transport computed the initial conditions for the subsequent 33-year transient simulation. Observations of precipitation, streamflow, evaporation, agricultural land use, and PVU brine pumping rates were used to specify appropriate boundary conditions to the model representing time-varying recharge, tributary streamflow, groundwater underflow, evapotranspiration (ET), and PVU pumping. Values for average monthly streamflow and TDS concentration at the upstream streamgage, the Dolores River at Bedrock (USGS streamgage 09169500), were specified as model input where the Dolores River enters Paradox Valley. Observed pumping from the PVU, water levels and TDS concentrations in groundwater, and streamflow and estimated TDS concentrations at the downstream streamgage, the Dolores River near Bedrock (USGS streamgage 09171100), were calibration targets that constrained the manual calibration of model parameters representing aquifer hydraulic conductivity, storage, streambed conductance, recharge, and (ET). Two primary model-calibration targets were the match between observed and simulated TDS mass flux from PVU pumping wells and the match between estimated and simulated TDS mass flux to the Dolores River. The simulated TDS mass withdrawn by pumping wells is calculated by the model as the product of the assigned pumping rate and simulated groundwater TDS concentrations. Because actual pumping rates were assigned as simulated values, the total simulated PVU pumping for the 33-year calibration is within 0.5 percent of the observed values. However, simulated concentrations and thus mass flux of TDS withdrawn by the PVU pumping wells were consistently about 26 percent less than observed values for all the simulated time periods (33-year simulation, pre-PVU, and post-PVU). The representation of brine inflow was explored through additional modeling to evaluate the effect of the simulated brine source on groundwater TDS concentrations. Results indicated that a saturated-salt constant-flux brine source best replicated the magnitude and transient pattern observed for TDS mass flux from PVU pumping wells. The simulated TDS mass flux to the Dolores River is compared to estimates based on observed streamflow and specific conductance (SC) data for the downstream streamgage. The calibrated model provided a close fit of simulated to measured streamflow at the downstream streamgage, and the calibrated model fit to estimated TDS concentrations at the downstream streamgage was reasonable. The greatest differences between simulated and estimated values occurred during drought periods from June 2000 to March 2003, May 2012 to June 2013, and October 2013 to October 2014, when simulated TDS concentrations in the river were greater than estimated concentrations. In general, simulated TDS mass flux to the river for the pre-PVU period is in good agreement with estimated values (2-percent difference), but the model overestimated TDS mass flux to the river by about 41 percent during the post-PVU period. The model uncertainty with respect to TDS mass flux to the river indicates other processes or model parameters not well represented by the model are affecting the system, especially during drought. During model calibration, the most sensitive parameters were identified as vertical hydraulic conductivity of the alluvial aquifer, conductance of the Dolores River streambed, ET extinction depth and rate, and recharge rate. Five 5-year scenarios of conditions for 2021–25 were simulated to assist evaluation of alternative strategies to manage the discharge of brine into the Dolores River. The first scenario simulates no PVU pumping and serves as a base case for comparison to the other scenarios. Two scenarios simulate the effects of varying withdrawal timing at an annual rate about one-third less than during 2010 through 2018. During high-flow spring snowmelt runoff periods when brine discharge is naturally minimized, PVU pumping does not substantially affect salinity in the Dolores River, and comparison of these two scenarios indicates that scheduling brine withdrawals during times of low river stage is nearly as effective at reducing TDS mass flux to the river as pumping brine year-round. Cessation of pumping during periods of high river stage may be advantageous for system maintenance, brine injection, and seismic-risk reduction. The fourth scenario tested the effect of reducing irrigation-return flow on brine discharge and predicted a slight reduction of TDS mass flux to the Dolores River, but not as great a reduction as that of using the PVU to remove brine. The fifth scenario simulated 5 years of drought conditions without PVU pumping and indicates brine discharge during drought about 15 percent greater than during average hydrologic conditions. Results from scenario 5 are consistent with the calibrated model results and indicate that aquifer properties and ET processes and parameters may be affecting simulation results during drought. The Paradox Valley groundwater model provides a reasonable overall match to observed conditions in the Dolores River. The model is useful for evaluating relative differences between brine management scenarios to inform PVU operational decisions and to identify gaps in data and process understanding. Representation of the brine source, hydraulic-conductivity parameters, and recharge and ET processes were identified as potential areas for additional field and modeling research. Additional research in the Paradox Valley might include field-data collection that provides additional information on the hydrogeologic framework, groundwater levels, groundwater TDS concentrations, stream characteristics, and aquifer properties. Additional modeling efforts could benefit from applying advanced tools for model development, calibration, and visualization including parameter-estimation and sensitivity analysis. Statistical evaluation of known model uncertainties such as hydraulic conductivity, streambed conductance, representations of the brine source, recharge, and ET could improve the match between simulated and estimated TDS mass flux from PVU pumping wells and to the Dolores River further informing model predictions and system understanding for the Paradox Valley.

Colorado↗

USGS2020

Evolution of faulting induced by deep fluid injection, Paradox Valley, Colorado

High‐pressure fluid injection into a subhorizontal confined aquifer at 4.3–4.6 km depth induced > 7000 "> > 7000 >7000 earthquakes between 1991 and 2012 within once seismically quiescent Paradox Valley in Colorado, with magnitudes up to M w "> M w Mw 3.9. Earthquake hypocenters expanded laterally away from the well with time, defining the margins of the aquifer pressurized by injection at the well. Within 5 km of the well, alignment of earthquake hypocenters defines strikes of nine vertical fault zones. Previous studies show that these fault zones predate injection, producing left‐stepping offsets in the normal faults of the Wray‐Mesa fault system that cradles Paradox Valley. Hypocenters, rakes, and strikes of 2041 well‐constrained focal mechanisms show that most injection‐related earthquakes occur where these vertical faults intersect the pressurized aquifer. Well‐defined focal mechanisms show that this induced seismicity consists of Riedel shear faults at acute angles to the strikes of these fault zones. These small faults develop an anastomosing fault structure of focal planes along each planar fault zone, as fluid injection continues, even as their hypocenters define a single planar fault zone. Failure conditions at each hypocenter are found using a fully coupled poroelastic analysis of stress induced by fluid injection, and this analysis indicates a minimum Coulomb failure condition of 0.1 MPa. This failure condition is primarily a result of aquifer pore‐fluid pressurization, as almost all well‐located seismicity is within the pressurized aquifer. Reducing the rate of injection and frequent well shutdowns in the second decade nearly eliminated induced seismicity, except very near the well where gradients in pressurization are the largest. Despite these decreases in failure conditions and seismicity, some fault zones continued to produce earthquakes larger than M 3 as injection continued.

California↗

USGS2020

High-resolution mapping of the freshwater-brine interface using deterministic and Bayesian inversion of airborne electromagnetic data at Paradox Valley, USA

Salt loads in the Colorado River Basin are a primary water quality concern. Natural groundwater brine discharge to the Dolores River where it passes through the collapsed salt anticline of the Paradox Valley in western Colorado is a significant source of salt to the Colorado River. An airborne electromagnetic survey of Paradox Valley has provided insights into the 3D distribution of brine in the surficial aquifer. A combination of stochastic and deterministic resistivity inversions were used to interpret the top of the freshwater-brine interface and to qualitatively describe the vertical salinity gradients across the interface. Low-resistivity regions indicative of brine occur near the land surface where brine discharges to the Dolores River and increase in depth several kilometers up-gradient along the axis of the valley. The most conductive parts of the brine plume are found in the areas below and adjacent to the river, suggesting that the brine becomes shallower and more concentrated as it reaches its natural discharge location. A significant freshwater lens overlying the brine west of the Dolores River is spatially correlated to the intermittent West Paradox Creek and agricultural irrigation. Below this lens, the transition from freshwater to brine appears to occur abruptly over a few meters and correlates to available well information. However, away from these regions and particularly with distance from the river, the freshwater-brine interface appears to be more diffuse.

Colorado↗

USGS2019

Controls on spatial and temporal variations of brine discharge to the Dolores River in the Paradox Valley, Colorado, 2016–18

The Paradox Valley in southwestern Colorado is a collapsed anticline formed by movement of the salt-rich Paradox Formation at the core of the anticline. The salinity of the Dolores River, a tributary of the Colorado River, increases substantially as it crosses the valley because of discharge of brine-rich groundwater derived from the underlying salts. Although the brine is naturally occurring, it increases the salinity of the Colorado River, which is a major concern to downstream agricultural, municipal, and industrial water users. The U.S. Geological Survey in cooperation with the Bureau of Reclamation conducted a study to improve the characterization of processes controlling spatial and temporal variations in brine discharge to the Dolores River. For the study, three geophysical surveys were conducted in March, May, and September 2017, and water levels were monitored in selected ponds and groundwater wells from November 2016 to May 2018. The study also utilized streamflow and specific conductance data from two U.S. Geological Survey streamflow-gaging stations on the Dolores River to estimate salt load to the river. River-based continuous resistivity profiling and frequency domain electromagnetic induction surveys made during low-flow conditions indicated a zone of brine-rich groundwater close to the riverbed along an approximately 4-kilometer reach of the river. Under high-flow conditions, the brine was depressed as much as 2 meters below the riverbed, and brine discharge to the river was reduced to a minimum. Direct current electrical resistivity surveys show that the freshwater lens overlying the brine is much thicker (up to 10 meters) on the west bank than on the east bank (less than 5 meters). A large low-conductivity anomaly at river distance 6,800 meters was observed in all surveys and may represent a freshwater discharge zone or a losing reach of the river. Filling and draining of the wildlife ponds on the west side of the river had a negligible effect on salt loads in the river during the study period. Groundwater monitoring showed there was active exchange of water between the river and the adjacent alluvial aquifer. When river stage was low, groundwater flowed towards the river, and brine discharge to the river increased. When the river stage was high, the gradient was reversed, and fresh surface water recharged the alluvial aquifer minimizing brine discharge. Most of the salt load to the river occurred during the winter and appeared to be enhanced by diurnal stage fluctuations. A conceptual model of brine discharge to the river is presented at three scales. Groundwater at the regional scale drives dissolution of salt in the Paradox Formation and flow of brine into the base of the alluvial aquifer. Surface water–groundwater interactions at the scale of the alluvial aquifer control brine discharge to the river seasonally and interannually. At the finest scale, diurnal fluctuations in river stage drive exchange of freshwater with saltier pore water in the hyporheic zone, which appears to increase brine discharge to the river during winter.

Colorado↗

USGS2005

Experimental investigation of CO2-brine-rock interactions at elevated temperature and pressure: Implications for CO2 sequestration in deep-saline aquifers

Deep-saline aquifers are potential repositories for excess CO 2 , currently being emitted to the atmosphere from anthropogenic activities, but the reactivity of supercritical CO 2 with host aquifer fluids and formation minerals needs to be understood. Experiments reacting supercritical CO 2 with natural and synthetic brines in the presence and absence of limestone and plagioclase-rich arkosic sandstone showed that the reaction of CO 2 -saturated brine with limestone results in compositional, mineralogical, and porosity changes in the aquifer fluid and rock that are dependent on initial brine composition, especially dissolved calcium and sulfate. Experiments reacting CO 2 -saturated, low-sulfate brine with limestone dissolved 10% of the original calcite and increased rock porosity by 2.6%. Experiments reacting high-sulfate brine with limestone, both in the presence and absence of supercritical CO 2 , were characterized by the precipitation of anhydrite, dolomitization of the limestone, and a final decrease in porosity of 4.5%. However, based on favorable initial porosity changes of about 15% due to the dissolution of calcite, the combination of CO 2 co-injection with other mitigation strategies might help alleviate some of the well-bore scale and formation-plugging problems near the injection zone of a brine disposal well in Paradox Valley, Colorado, as well as provide a repository for CO 2 . Experiments showed that the solubility of CO 2 is enhanced in brine in the presence of limestone by 9% at 25 °C and 6% at 120 °C and 200 bar relative to the brine itself. The solubility of CO 2 is enhanced also in brine in the presence of arkosic sandstone by 5% at 120 °C and 300 bar. The storage of CO 2 in limestone aquifers is limited to only ionic and hydraulic trapping. However, brine reacted with supercritical CO 2 and arkose yielded fixation and sequestration of CO 2 in carbonate mineral phases. Brine desiccation was observed in all experiments containing a discrete CO 2 phase, promoting porosity-reducing precipitation reactions in aquifers near saturation with mineral phases.

Colorado↗

USGS2003

Effect of the Paradox Valley Unit on the dissolved-solids load of the Dolores River near Bedrock, Colorado, 1988-2001

Discharge of brine with an average dissolved-solids concentration of about 256,000 milligrams per liter from alluvium in Paradox Valley, a collapsed salt anticline, substantially increases the dissolved-solids load of the Dolores River. In 1996, the Bureau of Reclamation began operation of the Paradox Valley Unit, a series of brine-withdrawal wells completed in alluvium along the Dolores River and a deep-injection well for the brine, to decrease flow of brine into the river. This report presents the findings of a study to determine the effectiveness of the Paradox Valley Unit from 1988 through September 2001. Differences in dissolved-solids load of the Dolores River between two gaging stations, one upstream and one downstream from the Paradox Valley Unit, indicate that an average dissolved-solids load of about 313 tons per day (an annual average of about 115,000 tons) was contributed by brine inflow to the Dolores River before operation of the Paradox Valley Unit began in July 1996. By September 30, 2001, the dissolved-solids load contributed by brine had declined to an average of about 29 tons per day? a decrease of about 90 percent. This decrease might have been facilitated by a decrease in precipitation and streamflow into the Paradox Valley during the last few years of the assessed period.

Colorado↗

USGS1992

Geochemical effects of deep-well injection of the Paradox Valley brine into Paleozoic carbonate rocks, Colorado, U.S.A.

Brine seepage into the Dolores River from ground water in Paradox Valley, Colorado constitutes a major source of salt to the Colorado River. Plants are enderway to remove this source of salt by drawing down the Paradox Valley brine (PVB) and forcibly injecting it into a deep disposal well (4.8 km). Experiments were conducted to determine the effects of deep-well injection of PVB. The results show that PVB is near saturation with anhydrite at 25°C, and that heating results in anhydrite precipitation. The amount and the rate at which anhydrite forms is temperature, pressure, and substrate dependent. Paradox Valley brine heated in the presence of Precambrian rocks from the drill core produces the same amount of anhydrite as PVB heated alone, but at a greatly accelerated rate. A 30% dilution of PVB with Dolores River water completely eliminates anhydrite precipitation when the fluid is heated with the Precambrian rocks. Interaction of PVB and Leadville Limestone is characterized by dolomitization of calcite by brine Mg which releases Ca to solution. This added Ca reacts with SO 4 to form increased amounts of anhydrite. A 20% dilution of PVB by Dolores River water has no effect on dolomitization and reduces the amount of anhydrite only slightly. A 65% dilution of PVB by Dolores River water still does not prevent dolomitization but does suppress anhydrite formation. Computer modeling of PVB by programs utilizing the Pitzer ion-interaction parameters is in general agreement with the experimental results. Ion-activity products calculated by both SOLMINEQ and PHRQPITZ are close to equilibrium with both anhydrite and dolomite whenever these phases are present experimentally, although the calculations over-estimate by a factor of 2 the degree of saturation. Some discrepancies in the calculated results between the two programs are due largely to differences in mineral solubility data.

Colorado↗

USGS1978

Evaluation of hydrogeologic aspects of proposed salinity control in Paradox Valley, Colorado

The salt load in the Dolores River increases by about 200,000 tons per year where it crosses Paradox Valley, Colorado, because of the discharge of a sodium chloride brine from an underlying aquifer. A ground-water management program to nearly eliminate this major source of salt, which eventually enters the Colorado River, can be designed on the basis of an accurate description of the hydrogeologic framework of Paradox Valley.

Colorado↗