Search USGS⌕ Search

Geology topics

M. Alisa Mast

Publications and source records attributed to M. Alisa Mast.

At least 19 recordsLinked to original sources
USGS2026

Peak-, mean-, and low-streamflow regional-regression equations for natural streamflow in central and western Colorado, 2019

The U.S. Geological Survey (USGS), in cooperation with the Colorado Department of Transportation, developed peak-, mean-, and low-streamflow regional-regression equations for estimating various statistics for natural streamflow in hydrologic regions of central and western Colorado. The peak-streamflow regression equations were developed using data from 418 streamgages, consisting of 15,202 years of record and a mean of approximately 36 years of record per streamgage. The mean- and low-streamflow regional-regression equations were developed using data from 323 streamgages where daily streamflow data were collected year-round. The annual exceedance-probability discharges for each streamgage were computed using the USGS software program PeakFQ. Mean monthly and 7-day minimum and maximum streamflows were computed using the USGS software program SWToolbox. Streamflow-duration values were computed using an R script. The regional-regression equations were determined using data for the period of record for a given streamgage through water year 2019. Geographic information systems datasets were used to develop 55 basin and 42 climatic characteristics, which were evaluated as candidate explanatory variables in the regression analysis. For the peak-streamflow regional-regression equations, the study area was divided into four hydrologic regions based on mean basin elevation, including the Plateau (less than 8,014 feet), Mid-Elevation (8,015 feet to 9,492 feet), Sub-Alpine (9,493 feet to 10,490 feet), and Alpine (greater than 10,490 feet) regions. For the peak-streamflow equations, the selection of basin and climatic characteristics was based on the 1-percent annual exceedance-probability discharge for each hydrologic region. For the mean streamflow, streamflow-duration values, and 7-day minimum and maximum streamflows, the study area was divided into four hydrologic regions based on river basin, including the (1) Colorado-East Slope Headwaters, (2) Green River, (3) Rio Grande, and (4) San Juan-Dolores. For mean streamflows, basin and climatic characteristics were evaluated separately for the annual period and each month for each hydrologic region. Regional regression equations published in this report are available for use in the USGS web-based program StreamStats.

Colorado↗
USGS2025

Environmental characterization of Blue Mesa Reservoir and potential causes of and management strategies for harmful algal blooms, 1970 through 2023, Curecanti National Recreation Area, Colorado

Blue Mesa Reservoir, in the Curecanti National Recreation Area, is the largest storage reservoir in Colorado and consists of three distinct basins: Iola (the shallowest), Cebolla, and Sapinero. After algal toxins were first documented in Iola basin in 2018, the U.S. Geological Survey began a study in cooperation with the National Park Service, Colorado River Water Conservation District, Upper Gunnison River Water Conservancy District, Gunnison County, Project 7 Water Authority, and Uncompahgre Valley Water Users Association to better understand occurrence of toxic cyanobacteria harmful algal blooms (cyanoHABs) and identify possible causal mechanisms to potentially inform management strategies. Toxic cyanoHABS occurred when the algal toxin microcystin exceeded a concentration of 8 micrograms per liter primarily in Iola basin in 2018 and 2020–22, years having some of the lowest reservoir water-level elevations (reservoir levels) since 1984. The toxic cyanoHABs started in mid-September and continued through the fall months. Algal abundance was greatest in Iola basin compared to Cebolla and Sapinero basins, with Aphanizomenon , a toxin-producing cyanobacterium, being the most abundant. During blooms, enhanced algal photosynthesis caused elevated pH and dissolved oxygen concentrations especially in Iola basin. Continuous monitor data in Iola basin indicated peaks in phycocyanin fluorescence, pH, and dissolved oxygen concentration that preceded the onset of toxic cyanoHABs by about 2 weeks potentially indicating a useful early warning monitoring strategy for future response to toxic cyanoHABs. Long-term trends showed increases in mean air and surface-water temperatures and chlorophyll- a concentrations in the reservoir but no change in nutrient inputs from major tributaries. In Iola basin, reservoir level was positively correlated with Secchi disk depth and inversely correlated with total phosphorus concentration. Because of its shallow depth, the effect of low reservoir levels may disproportionately affect Iola basin compared to other basins, resulting in algal blooms and toxin production especially at reservoir levels below about 7,470 feet above North American Vertical Datum of 1988. Elevated phosphorus at low reservoir level likely was primarily phosphorus contained in algal tissue. This report indicates that the main driver for recent toxic cyanoHABs in Iola basin is low reservoir level that likely causes favorable conditions (shallow and warm) for algal growth and increased recruitment of algae from bottom sediments such as during wind-driven turbulence. Control of external nutrients to the reservoir is unlikely to help control algal blooms because Aphanizomenon fixes nitrogen from the atmosphere, and there is an abundant geogenic source of phosphorus. Maintenance of reservoir levels greater than about 7,470 feet might help minimize the occurrence of toxic cyanoHABs. Additional data could help better understand how the timing and duration of reservoir levels below 7,470 feet contribute to toxic cyanoHABs.

Colorado↗
USGS2025

Multi-scale geophysical mapping of the brine and bedrock surfaces along the Dolores River, Paradox Valley, Colorado, December 2023

Total dissolved solids derived from salt dome–sourced brine in the underlying alluvial aquifer substantially increase with distance in the reach of the Dolores River that passes through Paradox Valley in southwestern Colorado. The area has been the site of salinity control operations since the 1990s to reduce salt loading to the downstream Colorado River. Previous airborne and ground/water-based electromagnetic (EM) geophysical data have successfully mapped the top of the brine surface, albeit with relatively coarse near-surface resolution and limited spatial coverage. This present December 2023 study used ground-based high-resolution EM and passive seismic (horizontal-to-vertical spectral ratio, HVSR) tools to map in detail the depth and thickness of the brine zone in the alluvial aquifer (top of the brine down to bedrock contact) in areas immediately surrounding the Dolores River where previous airborne EM (AEM) results indicated brine within 10 m of land surface. Results indicate the deepest bedrock is generally associated with the shallowest brine and local depressions in the collapse breccia (caprock to the Paradox Formation salt) may facilitate vertical migration of brine into the alluvial aquifer. Additionally, the ground-based EM mapping corroborated general patterns in depth to brine that were observed in previous AEM results while also revealing additional detail, including suspected focused brine discharge zones to the Dolores River. A river-based survey wherein EM data, channel depth and river water–specific conductance information were collected augmented these findings and indicated specific areas on both the western and eastern sides of the river where focused brine discharge may occur. This study comprises a large-scale, ground- and water-based geophysical mapping effort, including hundreds of HVSR soundings and 100s of kilometres of EM data, which were successfully translated into digital brine and bedrock surfaces that could be incorporated into groundwater modelling efforts, future well siting or other decision-making.

Colorado↗
USGS2025

Remote sensing of chlorophyll a and temperature to support algal bloom monitoring in Blue Mesa Reservoir, Colorado

We present methods to reconstruct historical chlorophyll a and surface water temperatures from satellite-based remote sensing products for Blue Mesa Reservoir, Colorado, to support algal bloom monitoring. A machine learning model was trained to construct chlorophyll a concentrations from Sentinel-2 satellite imagery and in situ measurements of chlorophyll a concentrations (out of bag RMSE = 1.9 μg/L, R 2 = 0.63) and reconstruct summertime chlorophyll a concentrations over the entire reservoir from 2016 through 2023. Concurrently, we developed an approach to retrieve remotely sensed water temperatures from the Landsat collection 2 provisional surface temperature product (MAE = 0.6°C) and reconstructed summertime surface water temperature records from 2000 through 2023. Finally, we demonstrate how the reconstructed chlorophyll a and temperature records can yield insight on reservoir dynamics. The chlorophyll a records indicate that algal blooms have a consistent spatial pattern across multiple years, initiating in the eastern end of the reservoir and spreading to the west over time. Water temperatures increased at a linearized rate of 0.3°C per decade from 2000 through 2023 and were inversely proportional to reservoir water surface elevation. Finally, mean summer remotely sensed chlorophyll a concentration had a moderately positive correlation with mean summer remotely sensed water temperature.

Colorado↗
USGS2024

Testing tree-ring cellulose δ18O with water isotopes for Holocene lake δ18O interpretations in the central Rocky Mountains USA

Stable isotopes of water preserved in geologic archives, primarily as oxygen (δ18O), have proven critical for documenting Earth’s climatic and hydrologic systems past and present. However, timescale differences of water isotope inputs to proxy systems and the signal embedded in long paleorecords often confound translation to observed hydroclimatic metrics. Here, a unique 20-year dataset of meteorology, hydrology, and the isotopic composition of weekly meteoric and surface water samples (δ18O, δ2 H) are combined with paleoclimate δ18O data from tree-ring cellulose and lake carbonate to better understand proxy signals of Upper Colorado river basin drought. Annual tree-ring cellulose δ18O from Picea engelmannii growing within a glacier-fed creek and a spring discharge area were used to derive annual source water δ18O using a cellulose source-water isotope model. Comparisons with the monitoring record indicates that tree-ring cellulose δ18O tracks variations in wet and dry hydroclimatic extremes. Source water isotopes are shown to reflect the hydroclimate of the current year and some number of previous years as an effective moisture-discharge proxy rather than a precipitation isotope proxy. Results contextualize Holocene lake carbonate δ18O data. The contemporary-to-paleo comparison identifies changes in seasonal precipitation extremes during recent millennia and several earlier arid and monsoon-dominated Holocene periods that exceed the arid maximum of the calibration period.

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↗
USGS2022

Estimating urban air pollution contribution to South Platte River nitrogen loads with National Atmospheric Deposition Program data and SPARROW model

Air pollution is commonly disregarded as a source of nutrient loading to impaired surface waters managed under the Clean Water Act per states’ 303(d) list programs. The contribution of air pollution to 2017–2018 South Platte River nitrogen (N) loads was estimated from the headwaters to the gage at Weldona, Colorado, USA (100 km downstream of Denver), using data from the National Atmospheric Deposition Program (NADP) and the SPAtially Referenced Regressions On Watershed attributes (SPARROW) model. The NADP offers wet-deposition raster created by spatial interpolation of data collected from regionally representative monitoring sites, excluding the influences from urban site data. For this study, NADP wet-deposition data obtained from sites within the Denver-Boulder, Colorado, urban corridor were included and excluded in new spatial interpolations of wet-deposition raster, which were used as input for SPARROW to model the influence of urban air pollution sources on South Platte River loads. Because urban air pollution is already incorporated into the NADP Total Deposition modeling methodology, dry N deposition was held constant for each SPARROW modeling scenario when dry deposition was included. By including the urban wet-deposition data in the model, estimated N loading to the South Platte River at Denver increased by 9–11 percent. Factoring in dry deposition at a 1:1.8 dry:wet ratio obtained from the results, urban air pollution was estimated to contribute as much as 20 percent of the nitrate Total Maximum Daily Load for Segment 14 of the South Platte River.

Colorado↗
USGS2021

Characterization of groundwater quality and discharge with emphasis on selenium in an irrigated agricultural drainage near Delta, Colorado, 2017–19

Selenium is a water-quality constituent of concern for aquatic ecosystems in the lower Gunnison River Basin. Selenium is derived from bedrock of the Mancos Shale and is mobilized and transported to groundwater and surface water by application of irrigation water. Although it is recognized that groundwater contributes an appreciable amount of selenium to surface water, few studies have addressed interactions between the two. The U.S. Geological Survey in cooperation with the Colorado Water Conservation Board conducted a study during 2017–19 to characterize the quality and quantity of groundwater discharging to an agricultural drainage near Delta, Colorado, locally known as Sunflower Drain. Water quality in the study area is characterized by high dissolved solids with elevated concentrations of selenium and nitrate resulting from dissolution of soluble salts in the Mancos Shale. Selenium concentrations have decreased by 50 percent since the early 2000s, possibly in response to irrigation system improvements. Stable water isotopes indicate streamflow is dominated by canal water during the irrigation season (April to October) and, during the nonirrigation season (November to March), is dominated by groundwater that has undergone some degree of evaporation. Pesticide and pharmaceutical compounds were infrequently detected, and results indicate they were derived from sources outside the study area such that they do not appear to be useful as tracers of groundwater sources. Stable isotopes of nitrate indicate that nitrate originates from the Mancos Shale, and the isotopic composition is enriched by denitrification in the groundwater system. Using a mass-balance approach, estimated groundwater discharge rates to Sunflower Drain ranged from 0.15 to 0.27 cubic feet per second per mile with one losing reach identified. Selenium, sulfate, and nitrate concentrations in groundwater estimated by mass-balance calculations were similar to concentrations measured in the Poly 17 observation well, located in a largely irrigated area in east tributary. One tributary reach had higher concentrations of selenium, sulfate, and nitrate likely reflecting localized inputs of more concentrated groundwater, similar to the concentrations in the Poly 7 observation well, which is downgradient from a residential area in the west tributary. Three pilot studies were conducted, including fiber optic distributed temperature sensing to detect groundwater discharge zones in the stream channel, a passive seismic technique to estimate depth to bedrock, and use of radon-222 as a geochemical tracer of groundwater discharge. All three techniques show promise as additional approaches for investigating groundwater discharge surface-water systems in irrigated drainage areas on Mancos Shale. The factors that affect groundwater movement mainly include when and where irrigation water is transported and applied, and the distribution of bedrock of the Mancos Shale and overlying alluvial deposits. The average groundwater recharge rate for the study area was estimated at 8.1 inches per year, based on mass balance calculations from synoptic survey data. Along the western tributary of Sunflower Drain, there was evidence that spills from the East Canal may recharge the groundwater aquifer adjacent to the stream channel. Groundwater movement to the stream channel may be controlled by the topography of the alluvial/bedrock interface or focused along human-made features, such as tile drains and ditches constructed around irrigated fields. On larger scales, the top of bedrock was also important, creating a topographic constriction that caused a zone of groundwater discharge. The groundwater system is complex, and further study could better define the system, possibly through application of a groundwater flow model and more extensive studies using some of the exploratory methods evaluated in this study.

Colorado↗
USGS2021

Water-quality change following remediation using structural bulkheads in abandoned draining mines, upper Arkansas River and upper Animas River, Colorado USA

Water-quality effects after remediating abandoned draining mine tunnels using structural bulkheads were examined in two study areas in Colorado, USA. A bulkhead was installed in the Dinero mine tunnel in 2009 to improve water quality in Lake Fork Creek, a tributary to the upper Arkansas River. Although bulkhead installation improved pH, and manganese and zinc concentrations and loads at the Dinero mine tunnel, water-quality degradation was observed at the nearby Nelson tunnel. Only manganese concentrations improved in Lake Fork Creek downstream from the tunnel. To improve water quality in Cement Creek, a tributary of the Animas River, multiple bulkheads were installed in mine tunnels during 1996–2003 and a water treatment plant operated from 1989 to 2003 to treat drainage from several draining tunnels. After bulkhead installation and cessation of active water treatment (about 2003), water quality (pH and dissolved copper, manganese, and zinc concentrations) degraded at the mouth of Cement Creek. The patterns and timing were similar to post-bulkhead increased discharge and trace-metal loads at non-bulkheaded tunnels indicating the bulkheads might have been the cause. Pre-1989 water-quality data for Cement Creek are scarce, although limited historical data indicate possible, slight improvement in only manganese concentrations after bulkhead installation. Increased zinc loads in Lake Fork Creek and decreased pH through time in Cement Creek may indicate increased groundwater discharge to the streams after bulkhead installation. In these two study areas, bulkheads did not substantially improve downstream water quality.

Colorado↗
USGS2020

Preferential elution of ionic solutes in melting snowpacks: Improving process understanding through field observations and modeling in the Rocky Mountains

The preferential elution of ions from melting snowpacks is a complex problem that has been linked to temporary acidification of water bodies. However, the understanding of these processes in snowpacks around the world, including the polar regions that are experiencing unprecedented warming and melting, remains limited despite being instrumental in supporting climate change adaptation. In this study, data collected from a snowmelt lysimeter and snowpits at meadow and forest-gap sites in a high elevation watershed in Colorado were combined with the PULSE multi-phase snowpack chemistry model to investigate the controls of meltwater chemistry and preferential elution. The snowdepth at the meadow site was 64% of that at the forest-gap site, and the snowmelt rate was greater there (meadow snowpit) due to higher solar irradiance. Cations such as Ca 2+ and NH 4 + "> 4+ were deposited mostly within the upper layers of both the meadow and forest-gap snowpacks, and acid anions such as NO 3 − "> 3− and SO 4 2 − "> 42− were more evenly distributed. The snow ion concentrations were generally greater at the forest-gap snowpit, except for NH 4 + "> 4+ , which indicates that wind erosion of wet and dry deposited ions from the meadow may have reduced concentrations of residual snow. Furthermore, at the forest-gap site, snow interception and scavenging processes such as sublimation, ventilation, and throughfall led to particular ion enrichment of Ca 2+ , Mg 2+ , K + , Cl − , SO 4 2 − "> 42− and NO 3 − "> 3− . Model simulations and observations highlight that preferential elution is enhanced by low snowmelt rates, with the model indicating that this is due to lower dilution rates and increased contact time and area between the percolating meltwater and the snow. Results suggest that low snowmelt rates can cause multiple early meltwater ionic pulses for ions subject to lower ion exclusion. Ion exclusion rates at the grain-size level have been estimated for the first time.

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↗
USGS2019

2017 Monitoring and tracking wet nitrogen deposition at Rocky Mountain National Park

The Colorado Department of Public Health and Environment (CDPHE), the National Park Service (NPS), and the U.S. Environmental Protection Agency (EPA) issued the Nitrogen Deposition Reduction Plan (NDRP) in 2007 to address the effects and trends of nitrogen deposition at Rocky Mountain National Park (RMNP). The agencies chose a glidepath approach to reduce wet nitrogen deposition to a level of 1.5 kilograms of nitrogen per hectare per year (kg N/ha/yr) by the year 2032 to protect sensitive ecosystems within RMNP from adverse effects. Another goal of the NDRP is to “reverse the trend of increasing nitrogen deposition at the park.” Trends in wet deposition data were analyzed at three sites in RMNP and three regional sites outside of the park. Wet nitrogen deposition (5-year rolling average) at Loch Vale in RMNP was 3.3 kg N/ha/yr, which is above the glidepath (2.4 kg N/ha/yr) in 2017. Wet nitrogen deposition has not decreased at RMNP or other sites in the region over the long-term. Ammonium concentrations showed a statistically significant increasing trend at all sites and nitrate concentrations showed a significant decreasing trend at four of the five sites over the period of record. In more recent years (2011-2017), wet nitrogen deposition showed no significant trend at monitoring sites in RMNP. Ammonium concentrations also showed no significant trend over the short-term, however nitrate concentrations did significantly decrease at two of the six sites.

Colorado↗
USGS2019

Science needs for continued development of total nitrogen deposition budgets in the United States

The objectives of this white paper are to describe the state of the science with respect to total Nr deposition budgets in North America and the research needed to improve these budgets from both a measurement and modeling standpoint. The document is intended to serve as a plan for TDep research activities but also, more broadly, to provide program managers, natural resource managers, policy makers and scientists with an understanding of the need for complete and accurate Nr deposition budgets to protect ecosystem health and human welfare, and the linkages between the underlying policy-relevant science questions and the specific knowledge and data gaps needed to improve Nr deposition budgets.

Report↗
USGS2018

Estimating metal concentrations with regression analysis and water-quality surrogates at nine sites on the Animas and San Juan Rivers, Colorado, New Mexico, and Utah

The purpose of this report is to evaluate the use of site-specific regression models to estimate metal concentrations at nine U.S. Geological Survey streamflow-gaging stations on the Animas and San Juan Rivers in Colorado, New Mexico, and Utah. Downstream users could use these regression models to determine if metal concentrations are elevated and pose a risk to water supplies, agriculture, and recreation. Multiple linear-regression models were developed by relating metal concentrations in discrete water-quality samples to continuously monitored streamflow and surrogate parameters (specific conductance, pH, turbidity, and water temperature) collected at the U.S. Geological Survey stations. Models were developed for dissolved and total concentrations of aluminum, arsenic, cadmium, copper, iron, lead, manganese, and zinc using water-quality samples collected from 2005 to 2017 by several Federal, State, Tribal, and local agencies using different collection methods and analytical laboratories. Model performance varied but, in general, models for dissolved metals did not perform as well as those for total metals. Dissolved metals generally were correlated to specific conductance or streamflow and total metals generally were better correlated with turbidity. Explanatory variables in the models reflected hydrologic and geochemical processes within the basin. A larger number of regression models were statistically significant for the most upstream sites, where metal concentrations were elevated by drainage from abandoned mines and mineralized bedrock. Models generally did not perform as well at downstream sites, especially for dissolved metals, which occurred at lower concentrations than at the upstream sites. In the lower reaches of the rivers, the input of more alkaline water from tributaries and groundwater reduced metal solubility and diluted metal concentrations. The number and distribution of samples in the calibration datasets also may have been a factor in model development. At some sites on the San Juan River, calibration datasets were more limited and did not represent the full range of observed hydrologic and water-quality conditions, especially during storm events in summer and fall. Recommendations for model use are given based on estimates of model precision, biases, and adequacy of the calibration datasets in terms of the number of samples and representativeness of the observed range of streamflow and water-quality conditions.

Colorado, New Mexico, Utah↗
USGS2018

Linking transit times to catchment sensitivity to atmospheric deposition of acidity and nitrogen in mountains of the western United States

Transit times are hypothesized to influence catchment sensitivity to atmospheric deposition of acidity and nitrogen (N) because they help determine the amount of time available for infiltrating precipitation to interact with catchment soil and biota. Transit time metrics, including fraction of young water ( F yw ) and mean transit time (MTT), were calculated for 11 headwater catchments in mountains of the western United States based on differences in the amplitude of the seasonal signal of δ 18 O in streamflow and precipitation. Results were statistically compared with catchment characteristics to elucidate controlling mechanisms. Transit times also were compared with stream solute concentrations to test the hypothesis that transit times are a primary influence on weathering rates and biological assimilation of atmospherically deposited N. Results indicate that transit times in the study catchments are strongly related to soil, vegetation, and topographic characteristics, with barren terrain (bare rock and talus) and steep slopes linked to high F yw and short MTT, whereas forest soil (hydrogroup B) was linked to low F yw and greater MTT. Concentrations of silicate weathering products (Na + and Si) were negatively related to F yw and barren terrain, and positively related to MTT and forest soil, supporting the concept that weathering fluxes and buffering capacity tend to be low in alpine areas due to short transit times. Nitrate concentrations were positively related to N deposition, catchment slope, and barren terrain, and negatively related to forest, indicating that hydrologic and/or biogeochemical processes associated with steep slopes limit uptake of atmospherically deposited N by biota. Interannual and seasonal variability in transit times and source water contributions in the study catchments was substantial, reflecting the influence of strong temporal variations in snowmelt inputs in high‐elevation catchments of the western United States. Results from this study confirm that short transit times in these areas are a key reason they are highly sensitive to atmospheric pollution and climate change.

Hydrological Processes↗
USGS2018

Spatial and temporal variation in sources of atmospheric nitrogen deposition in the Rocky Mountains using nitrogen isotopes

Variation in source areas and source types of atmospheric nitrogen (N) deposition to high-elevation ecosystems in the Rocky Mountains were evaluated using spatially and temporally distributed N isotope data from atmospheric deposition networks for 1995-2016. This unique dataset links N in wet deposition and snowpack to mobile and stationary emissions sources, and enhances understanding of the impacts of anthropogenic activities and environmental policies that mitigate effects of accelerated N cycling across the Rocky Mountain region. δ 15 N−NO 3 − at 50 U.S. Geological Survey Rocky Mountain Snowpack (Snowpack) sites ranged from −3.3‰ to +6.5‰, with a mean value of +1.4‰. At 15 National Atmospheric Deposition Program (NADP)/National Trends Network wet deposition (NADP Wetfall) sites, summer δ 15 N−NO 3 − is significantly lower ranging from −7.6‰ to −1.3‰ while winter δ 15 N−NO 3 − ranges from −2.6‰ to +5.5‰, with a mean value of +0.7‰ during the cool season. The strong seasonal difference in NADP Wetfall δ 15 N−NO 3 − is due in part to variation in the proportion of N originating from source regions at different times of the year due to seasonal changes in weather patterns. Snowpack NO 3 − and δ 15 N−NO 3 − are significantly related to NADP Wetfall (fall and winter) suggesting that bulk snowpack samples provide a reliable estimate at high elevations. Spatial trends show higher NO 3 − concentrations and δ 15 N−NO 3 − in the Southern Rocky Mountains located near larger anthropogenic N emission sources compared to the Northern Rocky Mountains. NADP Wetfall δ 15 N−NH 4 + ranged from −10‰ to 0‰, with no observed spatial pattern. However, the lowest δ 15 N−NH 4 + (−9‰), and the highest NH 4 + concentration (35 μeq/L) were observed at a Utah site dominated by local agricultural activities, whereas the higher δ 15 N−NH 4 + observed in Colorado and Wyoming are likely due to mixed sources, including fossil fuel combustion and agricultural sources. These findings show spatial and seasonal variation in N isotope data that reflect differences in sources of anthropogenic N deposition to high-elevation ecosystems and have important implications for environmental policy across the Rocky Mountain region.

Rocky Mountains↗