Geology topics

Research about Colorado

Source-linked reports with geographic coverage including Colorado.

At least 109 records · Page 6Linked to original sources

USGS2024

Evaluation and review of ecology-focused stream studies to support cooperative monitoring, Fountain Creek Basin, Colorado

The U.S. Geological Survey, in cooperation with Colorado Springs Utilities and Colorado Springs Stormwater Enterprise, synthesized previous studies and evaluated recent monitoring data to understand the distribution of fish and invertebrates in the Fountain Creek Basin and documented response to streamflow, water temperature, and water quality. The goal was to identify opportunities for aligning data collection to help maximize information gained from additional monitoring. Fifty-two publications were compiled from the literature that were completed within the study area between 1964 and 2022. Of these publications, 19 were fish and invertebrate focused. Overall, the distribution of fish and invertebrates in the Fountain Creek Basin has changed since the early 1900s. The occurrence of several fish species has increased or decreased since 2003, and a few species have not been collected in more than 100 years. Several mayfly, stonefly, and caddisfly taxa once common at several locations before 2000 are now rarely encountered, and those that now occur more frequently are associated with warmer-water streams. Decreasing invertebrate multimetric index values were noted at six locations, and the invasive Potamopyrgus antipodarum (New Zealand mud snail) is now established at two locations and occurs at several others, but in low numbers. Various streamflow characteristics were frequently noted to affect spatial and temporal patterns in fish and invertebrate communities, including early development and recruitment of Platygobio gracilis (flathead chub). Water quality and temperature contributed to patterns in aquatic communities, but less is known about the direct effects as these data were inconsistently available. Reach-scale habitat also contributed to patterns in aquatic communities, especially measures associated with the streambank, stream channel, and composition of streambed substrate. Moving forward, aligning consistent streamflow, water temperature, water quality, and geomorphic data collection at fish-, invertebrate-, and habitat-monitoring locations could maximize information gained from monitoring efforts and potentially inform evolving management activities and interests within the basin.

Colorado↗

USGS2024

Quantifying aspect-dependent snowpack response to high-elevation wildfire in the southern Rocky Mountains

Increasing wildfire frequency and severity in high-elevation seasonal snow zones presents a considerable water resource management challenge across the western United States (U.S.). Wildfires can affect snowpack accumulation and melt patterns, altering the quantity and timing of runoff. While prior research has shown that wildfire generally increases snow melt rates and advances snow disappearance dates, uncertainties remain regarding variations across complex terrain and the energy balance between burned and unburned areas. Utilizing paired in situ data sources within the 2020 Cameron Peak burn area on the Front Range of Colorado, U.S., during the 2021–2022 winter, we found no significant difference in peak snow water equivalent (SWE) magnitude between burned and unburned areas. However, the burned south aspect reached peak SWE 22 days earlier than burned north. During the ablation period, burned south melt rates were 71% faster than unburned south melt rates, whereas burned north melt rates were 94% faster than unburned north aspects. Snow disappeared 7–11 days earlier in burned areas than unburned areas. Net energy differences at the burned and unburned weather station sites were seasonally variable, the burned area snowpack lost more net energy during the winter, but gained more net energy during the spring. Increased incoming shortwave radiation at the burned site was 6 x more impactful in altering the net shortwave radiation balance than the decline in surface albedo. These findings emphasize the need for post-wildfire water resource planning that accounts for aspect-dependent differences in energy and mass balance to accurately predict snowpack storage and runoff timing.

Colorado↗

USGS2024

Airborne lidar accuracy analysis for dual photogrammetric and lidar sensor pilot project in Colorado, 2019

This report presents accuracy assessment results of the light detection and ranging (lidar) data collected in Colorado during a pilot project in fall 2019. The purpose of the pilot project was to assess the accuracy of lidar and imagery data collected simultaneously for the U.S. Department of Agriculture (USDA) National Agriculture Imagery Program and the U.S. Geological Survey National Geospatial Program 3D Elevation Program (3DEP). A multiagency group consisting of U.S. Department of the Interior agencies and USDA agencies participated in the effort. Department of the Interior agencies included Bureau of Land Management, National Park Service, and U.S. Geological Survey; USDA agencies included the Farm Services Agency, the Natural Resource Conservation Service, and U.S. Forest Service. This pilot project was designed to help determine if a lidar sensor system has the potential to meet future 3DEP topographic lidar collection requirements, ideally at the same altitudes and leaf-on times that National Agriculture Imagery Program is flown. The airborne sensor system from Leica Geosystems (part of Hexagon) (hereafter referred to as dual sensor system) was used in the pilot project and can collect imagery and three-dimensional point cloud data concurrently. This report examines the characteristics of lidar data from a geometric accuracy perspective. Field surveys were performed to evaluate the three-dimensional absolute and relative accuracy of the airborne lidar data and to determine if the data met 3DEP specifications.

Colorado↗

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↗

USGS2024

Wildland fire effects on sediment, salinity, and selenium yields in a basin underlain by Cretaceous marine shales near Rangely, Colorado

Understanding and quantifying soil erosion from rangelands is a high priority for land managers, especially in areas underlain by Cretaceous Mancos Shale, which is a natural source of sediment, salinity, and selenium to surface waters in many areas of western Colorado and eastern Utah. The purpose of this report is to present the results of a U.S. Geological Survey study that assessed sediment, salinity, and selenium yields after the Dead Dog wildfire (fire began June 11, 2017) in northwestern Colorado. Two methodologies were used to quantify erosion, with different data requirements and analytical complexity. The first approach was the use of a process-based erosion model, the Watershed Erosion Prediction Project, which uses inputs of climate, topography, vegetation, and soils data from existing datasets to predict erosion, making this approach easily extensible to other areas. The second approach required more complex data collection and was used to measure erosion and deposition by differencing digital elevation models created from uncrewed aerial vehicle imagery collected in 2016 (pre-fire) and 2021 (post-fire). Sediment, salinity, and selenium yields were calculated from the volumetric estimates of erosion from both methods, and a discussion of factors that may have contributed to overall findings, including vegetation, fire effects, and soil characteristics, is included. The two approaches yielded different outputs. Results from the Watershed Erosion Prediction Project model indicated that almost no erosion occurred after the Dead Dog fire. However, morphological changes in the study basin after the Dead Dog fire were visible in the pre- and post-fire imagery and measured in the digital elevation model differencing technique, with net erosion occurring in channel and landscape extents, though calculated erosion rates and salinity and selenium yields were relatively small. Visible and measured morphological changes consisted primarily of incision and deposition within stream channels and rill incision and expansion on steeper slopes. Widespread sheet erosion was not evident. Much of the new erosion originated within, and immediately below, previously vegetated areas that were then burned by the wildfire. Greater erosion rates and salinity and selenium yields were measured in the channel extent relative to the landscape extent. Calculated erosion rates ranged from 0.24 to 0.45 megagrams per hectare per year. These results indicate that the Dead Dog fire resulted in increased erosion in the study basin, yet these effects were relatively small based on the overall magnitude of modeled and measured erosion from the Watershed Erosion Prediction Project and the digital elevation model differencing technique. Minimal erosion in the basin is likely due to local site characteristics typical of soils derived from Mancos Shale, including the presence of robust physical crusts and biological soil crusts, and limitations of the methods based on data availability. Focusing uncrewed aerial vehicle flights on key areas (individual steep slopes, high-intensity burn areas, specific stream reaches) could likely increase understanding of erosional process with less effort and error than doing landscape-level flights.

Colorado↗

USGS2024

Evaluating distributed snow model resolution and meteorology parameterizations against streamflow observations: Finer Is not always better

Estimating snow conditions is often done using numerical snowpack evolution models at spatial resolutions of 500 m and greater; however, snow depth in complex terrain often varies on sub-meter scales. This study investigated how the spatial distribution of simulated snow conditions varied across seven model spatial resolutions from 30 to 1,000 m and over two meteorological data sets, coarser (≈12 km) and finer (4 km). Simulated snow covered area (SCA) was compared to remotely sensed SCA and simulated watershed mean peak snow water equivalent (SWE) was compared to four streamflow statistics representing different water management-relevant aspects of the hydrograph using non-parametric correlations. April 1 SWE tended to increase with model resolution, particularly below 4,000 masl. Finer meteorology simulations produced deeper April 1 SWE than coarser meteorology simulations. Finer resolution snow simulations tended to produce longer snowmelt durations and slower snowmelt rates than coarser resolution simulations. Finer resolution simulations had better agreement with SCA for both meteorology data sets, particularly at high and low elevations. However, finer resolution simulations did not generally outperform coarser simulations in snow versus streamflow statistic correlations. Snow versus streamflow correlations were most sensitive to meteorology, watershed properties, and then resolution. Watershed physiographic properties such as wetness index may increase snow versus streamflow metric correlations while elevation and slope may decrease correlations. At watershed scales, these results suggest that simulation resolution and choice of meteorology is less important than the physiographic properties of the watershed; however, if resolving snow distribution across the landscape is important, finer-resolution simulations are useful.

Colorado↗

USGS2024

On connecting hydro-social parameters to vegetation greenness differences in an evolving groundwater-dependent ecosystem

Understanding groundwater-dependent ecosystems (i.e., areas with a relatively shallow water table that plays a major role in supporting vegetation health) is key to sustaining water resources in the western United States. Groundwater-dependent ecosystems (GDEs) in Colorado have non-pristine temporal and spatial patterns, compared to agro-ecosystems, which make it difficult to quantify how these ecosystems are impacted by changes in water availability. The goal of this study is to examine how key hydrosocial parameters perturb GDE water use in time and in space. The temporal approach tests for the additive impacts of precipitation, surface water discharge, surface water mass balance as a surrogate for surface–groundwater exchange, and groundwater depth on the monthly Landsat normalized difference vegetation index (NDVI). The spatial approach tests for the additive impacts of river confluences, canal augmentation, development, perennial tributary confluences, and farmland modification on temporally integrated NDVI. Model results show a temporal trend (monthly, 1984–2019) is identifiable along segments of the Arkansas River at resolutions finer than 10 km. The temporal impacts of river discharge correlate with riparian water use sooner in time compared to precipitation, but this result is spatially variable and dependent on the covariates tested. Spatially, areal segments of the Arkansas River that have confluences with perennial streams have increased cumulative vegetation density. Quantifying temporal and spatial dependencies between the sources and effects of GDEs could aid in preventing the loss of a vulnerable ecosystem to increased water demand, changing climate, and evolving irrigation methodologies.

Colorado↗

USGS2024

A purple-colored 1M mica clay from Silverton, Colorado

A purple-colored clay of 1M mica polymorph in association with lead-zinc ore was collected from an abandoned mine dump near Silverton, Colorado. Electron micrographs show that the crystallites of the clay are less than 2µ in size and have poorly developed hexagonal outlines. Differential thermal and i.r. absorption analyses indicate similarity with those of muscovite. The structural formula of the 1M mica polymorph is (K 1.45 Na 0.02 ) (Al 3.75 Mg 0.19 ) (Si 6.90 Al 1.10 ) O 20 (OH) 4 . The purple color may be related to trace amounts of manganese in the clay.

Colorado↗

USGS2024

Evaluating post-wildfire debris-flow rainfall thresholds and volume models at the 2020 Grizzly Creek Fire in Glenwood Canyon, Colorado, USA

As wildfire increases in the western United States, so do postfire debris-flow hazards. The U.S. Geological Survey (USGS) has developed two separate models to estimate (1) rainfall intensity thresholds for postfire debris-flow initiation and (2) debris-flow volumes. However, the information necessary to test the accuracy of these models is seldom available. Here, we studied how well these models performed over a 2-year period in the 2020 Grizzly Creek Fire burn perimeter in Glenwood Canyon, Colorado, USA, through the development of a debris-flow response inventory. The study area had the advantage of a network of 11 rain gauges for rainfall intensity measurements and repeat lidar data for volume estimates. Our observations showed that 89 % of observed debris flows in the first year postfire were triggered by rainfall rates higher than the fire-wide rainfall threshold produced by the current USGS operational model (M1). No debris flows were observed in the second year postfire, despite eight rainstorms with intensities higher than the modeled rainfall threshold. We found that the operational model for debris-flow initiation rainfall thresholds works well in this region during the first year but may be too conservative in year 2 due to vegetation recovery and sediment depletion. However, rainfall thresholds in the second year can be improved by using updated remote sensing imagery to recalculate the debris-flow initiation likelihood with the M1 model. The current volume model overestimates for this region by a median value of 4.4 times. However, the offset between the model estimates and observations is approximately linear, and the volumes from the Grizzly Creek debris flows had a similar magnitude to historic postfire debris flows in the region. Consequently, the current volume model could be adjusted with a regional correction factor.

Colorado↗

USGS2024

Using the TSS-RESTREND methodology to diagnose post-reclamation vegetation trends on the western slope of Colorado

We evaluated the use of the time series segmented residual trends (TSS-RESTREND) methodology to analyze plant community trends after oil and gas reclamation. We focused on reclaimed well pads managed by the Bureau of Land Management in northwestern Colorado. We assessed whether TSS-RESTREND could detect postreclamation changes in a plant community and if such changes corresponded with management actions. We used precipitation data and the greenness-to-cover index to calculate the residuals of the vegetation–precipitation relationship (VPR residuals). The VPR residuals represent plant community trends caused by disturbance or management actions and not by precipitation. We then used breaks for additive season and trend and the Chow test on the VPR residuals of each well pad to identify abrupt changes in plant community composition from 2000 to 2020. Afterward, we applied a segmented residual trend (RESTREND) analysis to the VPR residuals before and after an identified breakpoint or a singular RESTREND when no significant breakpoint was found to determine if reclamation had an effect on vegetation response to precipitation. We found a slight positive increase in VPR residuals over time since reclamation, indicating a more positive response to precipitation over time. In addition, well pads with lower aridity index values had a small positive trend in VPR residuals over time, suggesting the negative impact of aridity on plant community composition diminishes with increasing time since reclamation. To further understand the connection between management actions and outcomes, we compared findings from TSS-RESTREND with aerial imagery and well pad documentation. With this information, we categorized the well pads into six groups based on reclamation outcomes. This approach provided insights into the effects of management actions on recovery. Overall, TSS-RESTREND methodology can help identify changes in plant community composition over time, enhancing our understanding of plant community dynamics in these severely degraded areas.

Colorado↗

USGS2024

Potential climate and human water-use effects on water-quality trends in a semiarid, western U.S. watershed: Fountain Creek, Colorado, USA

Nutrients, total dissolved solids (TDS), and trace elements affect the suitability of water for human and natural needs. Here, trends in such water-quality constituents are analyzed for 1999–2022 for eight nested monitoring sites in the 24,000 km 2 Fountain Creek watershed in Colorado, USA, by using the weighted regressions on time, discharge, and season (WRTDS) methodology. Fountain Creek shares characteristics with other western U.S. watersheds: (1) an expanding but more water-efficient population, (2) a heavy reliance on imported water, (3) a semiarid climate trending towards warmer and drier conditions, and (4) shifts of water from agricultural to municipal uses. The WRTDS analysis found both upward and downward trends in the concentrations of nutrients that reflected possible shifts in effluent management, instream uptake, and water conservation by a watershed population that grew by about 40%. Selenium, other trace elements, and TDS can pose water-quality challenges downstream and their concentrations were found to have a downwards trend. Those trends could be driven by either a warming and drying of the local climate or decreased agricultural irrigation, as both would reduce recharge and subsequent mobilization from natural geologic sources via groundwater discharge. The patterns illustrate how changes in climate and water use may have affected water quality in Fountain Creek and demonstrate the patterns to look for in other western watersheds.

Colorado↗

USGS2024

Assessment and characterization of ephemeral stream channel stability and mechanisms affecting erosion in Grand Valley, western Colorado, 2018–21

The Grand Valley in western Colorado is in the semiarid Southwest United States. The north side of the Grand Valley has many ungaged ephemeral streams, which are of particular interest because (1) the underlying bedrock geology, Late Cretaceous Mancos Shale, is a sedimentary rock deposit identified as a major salinity contributor to the Colorado River and (2) despite infrequent streamflows of short duration, monsoon-derived floods in these ephemeral streams can carry substantial amounts of sediment downstream, affecting upstream and downstream banks and channel cross sections. The study area is of interest, because salinity, or the total dissolved solids concentration, in the Colorado River causes an estimated $300 million to $400 million per year in economic damages in the United States, and it is estimated 62 percent of the Upper Colorado River Basin’s total dissolved solid loads originate from geologic sources. In an effort to minimize salt contributions to the Colorado River from public lands administered by the Bureau of Land Management, a comprehensive salinity control approach is typically used to reduce nonpoint sources of salinity through land management techniques and practices. In 2018, the U.S. Geological Survey, in cooperation with the Bureau of Land Management, began an assessment of ephemeral streams located on the north side of the Grand Valley, western Colorado, to characterize stream channel stability and identify mechanisms affecting erosion. The U.S. Geological Survey developed a method for automatically extracting channel cross-section geometry from existing remotely sensed terrain models. Based on estimated flood stage and surrogate streamflows, hydraulic characteristics were calculated. Furthermore, the channel geometries and hydraulic characteristics were used to estimate channel stability using a statistical model. Cross-section stabilities were determined from a stream channel stability assessment for a subset of 1,406 visited (field observed) locations out of 13,415 cross sections, which were delineated from remotely sensed terrain models. The application of Manning’s resistance equation in combination with multiple logistic regression models demonstrated channel stability can be estimated with a 0.845 goodness of fit for a validation dataset when using a combination of drainage area, width-to-depth ratio, sinuosity, and shear stress as the explanatory variables. Stream channel stability was extrapolated for 13,415 unvisited (not field observed) cross sections using the multiple logistic regression model and defined explanatory variables. Mapping of the ephemeral streams and their associated stabilities may be used by the Bureau of Land Management to prioritize areas for remediation or changes in management strategies to reduce sediment and salinity loading to the Colorado River. The study found channel stability within the ephemeral streams to be spatially variable, longitudinally discontinuous, and dictated by changes in channel bed slope. The stable ephemeral streams were relatively wide and shallow and often had smaller drainage areas with less potential for producing shear stresses capable of overcoming channel adhesion. A change in channel bed slope can provide the means necessary to generate shear stresses appropriate to initiate erosion and a subsequent stability transition to incising channels. Channel widening happens when either or both banks of an incising channel reach a critical height for mass wasting, or when channel curvature causes higher sidewall stress. Regardless, widening channels can promote increases in sinuosity and subsequently reduce steep channel bed slopes. Consequently, stable and widening channels can have comparable bed slopes, making channel bed slope a poor explanatory variable to predict channel stability overall, despite its function to initiate channel instability. The results were based on a surrogate 0.10 annual exceedance probability (AEP; return period equal to the 10-year flood) interval streamflow, although it was recognized fluctuations in streamflow would also affect channel stability. Past and current changes within the study area affect streamflow; therefore, mechanisms affecting erosion include land use disturbances, soil compaction, loss of vegetation cover, drought, less frequent and more extreme precipitation, and fires—which all intensify the potential runoff and erosion within the study area.

Colorado↗

USGS2024

Thick- and thin-skinned contractional styles and the tectonic evolution of the northern Sangre de Cristo Mountains, Colorado, USA

The Sangre de Cristo Mountains of southern Colorado and northern New Mexico, USA, contain an unusual combination of thick- and thin-skinned contractional structures involving both basement and cover rocks in the Laramide Rocky Mountain foreland. These structures are truncated by down-faulted extensional basins to the east and west. Together with synorogenic sediments, these structures preserve a record of the rise of the Ancestral Rocky Mountains, the Laramide orogeny, and Rio Grande rifting. Laramide structures within the mountains provide clues to processes that link the three events and to necessary conditions for thin-skinned and thick-skinned contractional structures to form together in continental interiors. To examine the full variety of structural styles, a portion of the northern Sangre de Cristo fold- and-thrust belt in Colorado was described and interpreted using geologic maps and structural cross-sections. Stratigraphic relations of the Ancestral Rocky Mountain highlands and basin fill were reconstructed from existing maps. These relations allow identification of faults inherited from the Ancestral Rocky Mountains, differentiation of thrust sheets, and in some cases, estimation of the magnitude of displacement. To examine relations between Laramide thrusts and Rio Grande rifting, kinematic data were collected from a thrust fault adjacent to rift faults. Three thrust fault styles were recognized: thin-skinned basement, thin-skinned cover rocks, and thick-skinned basement. Thin-skinned thrusts arising from a hinterland beneath the present San Luis Valley carried sheets of Proterozoic basement rocks northeast over a Laramide foreland. These basement thrusts are interpreted to be faults of the Ancestral Rocky Mountains that reactivated during the Laramide orogeny. The Laramide foreland consists of thin-skinned thrusts and folds in sedimentary cover rocks as young as 49 Ma. Both thin-skinned thrusts in basement and cover rocks are bounded by thick-skinned basement thrusts that moved intermittently throughout the Laramide orogeny. We infer that thin-skinned thrusts form in continental interiors where deformation is focused in weak strata of thick basin fill and in fluid-reaction weakened preexisting faults in basement rocks. Both conditions are met in the Sangre de Cristo Mountains. Basement thrusts adjacent to the San Luis Valley contain evidence of plastic contractional microstructures overprinted by extensional microstructures that may record the transition from Laramide contraction to Rio Grande extension of the crust.

Colorado↗

USGS2024

Deep resistivity geophysics of the San Juan–Silverton caldera complex, San Juan County, Colorado (USA)

Magnetotelluric (MT) and audiomagnetotelluric (AMT) data are used to better understand the subsurface geology and mineral resources in the San Juan–Silverton caldera complex located near Silverton, Colorado, western United States, as part of the extensive southern Rocky Mountains volcanic field that covers much of southwestern Colorado and northern New Mexico. Seven MT and AMT profiles of varying lengths image resistivity structure to depths of ~5 km. The AMT inversion models characterize geophysical responses of near-surface lithologies, structures, and mineralized systems and also help corroborate airborne electromagnetic data at shallow levels. The MT inversion models extend our depth of investigation from near the surface to great depths (~5 km) and help to form hypotheses about roots of the hydrothermal plumbing that fed shallower mineralized systems. Subsurface high resistivities occur beneath intermediate-composition lava flows and Proterozoic units. Subsurface moderate- to low-resistivity values may reflect hydrothermal plumbing that served as flow paths for mineralizing fluids and metallic ore formation. The model interpretations presented in this study could be utilized in remediation planning or mineral resource applications. The methods used could be applied to other watersheds with similar volcanic environments containing acid-generating historical mines or hydrothermally altered and mineralized source rocks.

Colorado↗

USGS2024

Groundwater hydrology, groundwater and surface-water interactions, aquifer testing, and groundwater-flow simulations for the Fountain Creek alluvial aquifer, near Colorado Springs, Colorado, 2018–20

From 2018 through 2020, the U.S. Geological Survey, in cooperation with the Air Force Civil Engineering Center, conducted an integrated study of the Fountain Creek alluvial aquifer located near Colorado Springs, Colorado. The objective of the study was to characterize hydrologic conditions for the alluvial aquifer pertinent to the potential for transport of solutes. Specific goals of this report were to characterize the groundwater hydrology of the area, to quantify groundwater and surface-water interactions, to estimate hydraulic properties of the aquifer using aquifer testing, and to complete numerical simulations of groundwater flow. Synoptic groundwater-level elevation measurements completed throughout this study, and as part of other U.S. Geological Survey programs between 1994 and 2020, indicate groundwater-level elevations fluctuate on annual and interannual timeframes. Groundwater-level fluctuations likely were caused by temporally variable groundwater recharge and discharge components in the area, with many wells showing maximum groundwater-level elevations during the winter months (November through March). From an interannual perspective, groundwater-level fluctuations appear to have reached maximum values during 2000 to 2003, decreased during 2003 to 2006, and remained relatively constant since that time, with the exception of several wells which have displayed rising groundwater-level elevations since 2018. Spatial evaluation of groundwater-level elevations indicates groundwater flow is generally from northeast to southwest within the vicinity of several alluvial paleochannels occurring along the northeastern margin of the aquifer. Within the center of the aquifer along Fountain Creek, groundwater flow is generally from north to south, approximately paralleling surface-water flow. To quantitatively understand the potential effect of groundwater recharge and groundwater pumping on fluctuations in groundwater-level elevation, a statistical transfer-function-noise model was applied. Results of the statistical model indicate throughout most of the aquifer, fluctuations were primarily the result of recharge seasonality. In the main stem of the aquifer where groundwater pumping wells were more concentrated, however, groundwater-level elevation fluctuations were also attributable to groundwater pumping through time. Three-dimensional evaluation of the aquifer geometry near Fountain Creek was combined with synoptic streamflow measurement and accounting of stream gains and losses to evaluate groundwater and surface-water interactions in the study area. Streamflow gain or loss calculations indicate Fountain Creek both gains from and loses flow to the alluvial aquifer, and gaining or losing reaches of the stream may be partially controlled by the depth to bedrock near the stream. Reaches with streamflow gains tend to coincide with areas where the estimated depth to bedrock is decreasing, meaning the alluvial aquifer is likely thinning in these areas and groundwater-flow paths may be converging and discharging groundwater to the stream. Losing reaches tended to coincide with locally greater depth to bedrock where the alluvial aquifer is likely thicker and has greater storage potential for surface water lost from Fountain Creek. Results of aquifer testing indicate hydraulic conductivity, estimated from slug tests and single-well pumping tests, ranged from 0.32 to 1,410 feet per day (ft/d) and 4.13 to 664 ft/d, respectively. These results are similar to the range of values from previous aquifer tests in the study area. Hydraulic conductivities from aquifer testing for this study were generally greater than the estimates of previous slug tests and had a mean value less than the estimates from previous pumping tests. Spatial evaluation of aquifer testing results indicates hydraulic conductivity tends to be greater in the main stem of the alluvial aquifer and lower in paleochannels upgradient from the main stem of the aquifer. The spatial variation in hydraulic conductivity may be attributed to the geomorphologic processes that formed the alluvial aquifer. Compacted sediment in the paleochannels has not been potentially transported sufficient distance to cause grain-size sorting, resulting in a poorly sorted deposit and lower hydraulic conductivities. In the central portion of the alluvial aquifer, near Fountain Creek, the sediments have been transported farther from their source areas and are likely better sorted, removing finer grained sediments that would cause lower hydraulic conductivity. A numerical groundwater-flow model was calibrated for the Fountain Creek alluvial aquifer for 2000–19 to simulate water-budget components, groundwater-flow directions, and groundwater-flow paths. The model simulated precipitation recharge, groundwater and surface-water interactions, evapotranspiration, high-volume groundwater pumping by pumping wells, and external inflows and outflows occurring along the boundaries of the alluvial aquifer. Model calibration was completed using manual and automated approaches, the latter of which assisted in quantifying model results sensitivity to input parameters. The calibrated model corresponds well with groundwater-level elevation observations, with a mean residual (observed minus simulated groundwater-level elevation) equal to −0.60 feet. Simulated groundwater base flow to streams was typically within 10 percent of base flow estimated by independent methods. Groundwater and surface-water interactions represented the largest water-budget components of the aquifer, with the second largest groundwater discharge component coming from pumping wells. Groundwater and surface-water interactions represent both the largest gain and loss terms in the water budget, because these interactions differ spatially, meaning in some areas of the model domain groundwater is being recharged by streams, whereas in other areas, groundwater is discharged to streams. Estimates of advective groundwater-flow paths indicate pumping wells may capture groundwater recharged from losing streams and groundwater that flows into the main stem of the alluvial aquifer from paleochannels.

Colorado↗

USGS2024

Climate-driven increases in stream metal concentrations in mineralized watersheds throughout the Colorado Rocky Mountains, USA

Increasing stream metal concentrations apparently caused by climate warming have been reported for a small number of mountain watersheds containing hydrothermally altered bedrock with abundant sulfide minerals (mineralized watersheds). Such increases are concerning and could negatively impact downstream ecosystem health, water resources, and mine-site remediation efforts. However, the pervasiveness and typical magnitude of these trends remain uncertain. We aggregated available streamwater chemistry data collected from late summer and fall over the past 40 years for 22 mineralized watersheds throughout the Colorado Rocky Mountains. Temporal trend analysis performed using the Regional Kendall Test indicates significant regional upward trends of ∼2% of the site median per year for sulfate, zinc, and copper concentrations in the 17 streams affected by acid rock drainage (ARD; median pH ≤ 5.5), equivalent to concentrations roughly doubling over the past 30 years. An examination of potential load trends utilizing streamflow data from eight “index gages” located near the sample sites provides strong support for regionally increasing sulfate and metal loads in ARD-affected streams, particularly at higher elevations. Declining streamflows are likely contributing to regionally increasing concentrations, but increasing loads appear to be on average an equal or greater contributor. Comparison of selected site characteristics with site concentration trend magnitudes shows the highest correlation for mean annual air temperature and mean elevation (R 2 of 0.42 and 0.35, respectively, with all others being ≤0.14). Future research on climate-driven controlling mechanisms should therefore focus on processes such as melting of frozen ground directly linked to site mean temperature and elevation.

Colorado↗

USGS2024

Urbanization and water management control stream water quality along a mountain to plains transition

Urbanization can have substantial effects on water quality due to altered hydrology and introduction of constituents to water bodies. In arid and semi-arid environments, streams are further stressed by dewatering as a result of diversions. We conducted a high-resolution synoptic survey of two streams in Colorado, USA that transition abruptly from granitic/metamorphic forested mountains to sedimentary urbanized plains and are both highly managed for water supply, yet differ in degree of urbanization. A synoptic mass balance approach developed for mine drainage applications was adapted to elucidate effects of urbanization, geology, and diversions on stream chemistry during baseflow conditions. Urbanization was a more important driver of stream concentrations than geology. The urban area was a strong source of bromide, calcium, chloride, and manganese, while lanthanum and dissolved organic carbon were primarily sourced from the mountains. A majority of streamflow was removed by diversions near the mountains/plains interface. Groundwater accounted for 31% of the subsequent flow increase to the urbanized stream, and delivered at least 33% of chloride loading. Constituents that were primarily urban-derived (bromide, calcium, chloride, and manganese) were 2–3 times higher in the urban region due to diversions; without diversions, stream water quality would have largely retained characteristics of forested streams through the urban reach. This study provides insights into processes that affect water quality in highly managed streams of the semi-arid western USA.

Colorado↗