Geology topics

Graham A. Sexstone

Publications and source records attributed to Graham A. Sexstone.

At least 19 recordsLinked to original sources

USGS2025

Fine-resolution satellite remote sensing improves spatially distributed snow modeling to near real time

Given the highly variable distribution of seasonal snowpacks in complex mountainous environments, the accurate snow modeling of basin-wide snow water equivalent (SWE) requires a spatially distributed approach at a sufficiently fine grid resolution (<500 m) to account for the important processes in the seasonal evolution of a snowpack (e.g., wind redistribution of snow to resolve patchy snow cover in an alpine zone). However, even well-validated snow evolution models, such as SnowModel, are prone to errors when key model inputs, such as the precipitation and wind speed and direction, are inaccurate or only available at coarse spatial resolutions. Incorporating fine-spatial-resolution remotely sensed snow-covered area (SCA) information into spatially distributed snow modeling has the potential to refine and improve fine-resolution snow water equivalent (SWE) estimates. This study developed 30 m resolution SnowModel simulations across the Big Thompson River, Fraser River, Three Lakes, and Willow Creek Basins, a total area of 4212 km 2 in Colorado, for the water years 2000–2023, and evaluated the incorporation of a Moderate Resolution Imaging Spectroradiometer (MODIS) and Landsat SCA datasets into the model’s development and calibration. The SnowModel was calibrated spatially to the Landsat mean annual snow persistence (SP) and temporally to the MODIS mean basin SCA using a multi-objective calibration procedure executed using Latin hypercube sampling and a stepwise calibration process. The Landsat mean annual SP was also used to further optimize the SnowModel simulations through the development of a spatially variable precipitation correction field. The evaluations of the SnowModel simulations using the Airborne Snow Observatories’ (ASO’s) light detection and ranging (lidar)-derived SWE estimates show that the versions of the SnowModel calibrated to the remotely sensed SCA had an improved performance (mean error ranging from −28 mm to −6 mm) compared with the baseline simulations (mean error ranging from 69 mm to 86 mm), and comparable spatial patterns to those of the ASO, especially at the highest elevations. Furthermore, this study’s results highlight how a regularly updated 30 m resolution SCA could be used to further improve the calibrated SnowModel simulations to near real time (latency of 5 days or less).

Colorado↗

USGS2025

Latitudinal gradients of snow contamination in the Rocky Mountains associated with anthropogenic sources

Seasonal snow is an important source of drinking water and recreation, and for agriculture in the Rocky Mountain region. Monitoring snow-water quality can inform on the effects to the albedo and energy balance of the snowpack, and the sources of natural and anthropogenic aerosol and gases. This study analyzed metals in the seasonal snowpack from water year (WY) 2018 for 49 sites. Calcium, lanthanum, and cerium concentrations support the importance of mineral dust to the southern Rocky Mountains. Mercury (Hg), zinc (Zn), and cadmium (Cd) concentrations showed a similar spatial pattern to mineral dust, whereas antimony (Sb) concentrations were highest in the northern Rocky Mountains. To assess the relative contributions from dust versus anthropogenic contaminant sources, enrichment factors (EF) were calculated, with values above 10 indicating anthropogenic contamination. For Cd, Hg, Sb, and Zn, EF values exceeded 10 at northern sites. These observations were compared to spatial trends of EF values of Hg from WY2009 to WY2018, regional monitoring networks, and back trajectory analyses. The agreement between these datasets revealed temporally consistent contaminant sources and/or transport processes to the northern Rocky Mountains snowpack. Sources include current and historical mining and smelting in the region. Strategies to limit the emissions of these metals to the Northern Rockies could benefit from focusing on remediation of contaminated sites, and continued monitoring and mitigation of active mining and smelting.

Colorado , Idaho, Montana, New Mexico, Utah, Wyomi↗

USGS2024

Streamflow timing and magnitude during snow drought depend on snow drought type and regional hydroclimate

Communities around the world rely on snowmelt to meet water demands, and periods of lower than normal snow accumulation, snow droughts, can decrease water supplies. Leveraging 172 minimally disturbed and seasonally snow-covered watersheds, we developed an approach to examine the effects of cool & dry, warm & dry, and warm & wet snow droughts on streamflow timing and magnitude by hydrologic region. Our results showed all types of snow droughts in all regions correlate with lower annual streamflow, lower maximum and minimum flows, and lower runoff ratios, with more numerous low flow days and earlier streamflow timing. However, departures from non-snow drought conditions differed substantially between drought types and regions. Consecutive snow droughts further reduced runoff ratios and increased low flow days, likely due to additional subsurface storage depletion. With warm snow drought occurrence expected to increase, we discuss impacts for water management systems whose design specifications may not reflect the changing hydroclimate.

Hydrological Sciences Journal↗

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

Modeling forest snow using relative canopy structure metrics

Snow and watershed models typically do not account for forest structure and shading; therefore, they display substantial uncertainty when attempting to account for forest change or when comparing hydrological response between forests with varying characteristics. This study collected snow water equivalent (SWE) measurements in a snow-dominated forest in Colorado, the United States, with variable canopy structure. The SWE measurements were integrated with 1 m Lidar derived canopy structure metrics and incoming solar radiation to create empirical SWE offset equations for four canopy structure groupings (forest gaps, south-facing forest edges, north-facing forest edges, and the interior forest) that varied in size compared to an open area. These simple equations indirectly integrate terrain shading and canopy shading and were able to estimate 40 to 70% of SWE variation in a heterogenous forested environment. The equations were then applied to a snow melt model with a 100 m grid size by applying the area-weighted average of SWE offsets from the four canopy structure groupings in each model cell. This tiled model configuration allowed for the model to better represent the subgrid heterogeneity of a forest environment that can be seen through an ensemble or range of potential outputs rather than a singular estimate.

Water↗

USGS2023

Aquatic carbon export and dynamics in mountain headwater streams of the western U.S.

Mountain headwater streams actively cycle carbon, receiving it from terrestrial landscapes and exporting it through downstream transport and gas exchange with the atmosphere. Although their importance is now widely recognized, aquatic carbon fluxes in headwater streams remain poorly characterized. In this study, aquatic carbon fluxes were measured in 15 mountain headwater streams and were used in a geostatistical analysis to determine how landscape characteristics influence aquatic carbon fluxes. In-stream sensors were used to measure fluorescent dissolved organic matter (fDOM) (a surrogate for dissolved organic carbon (DOC)) at a subset of sites to characterize dynamic responses to hydroclimatic events. Wetlands have a positive influence on aquatic carbon fluxes, whereas perennial snow/ice has the opposite effect, reflecting differences in soil organic matter content in these landscapes. Mean annual temperature (MAT) has a complex influence on DOC, with peak DOC exports in basins with MAT of 0–2°C. Precipitation has a strong positive influence on aquatic carbon fluxes, and declining snowpacks in the western United States may reduce future aquatic carbon exports. fDOM (and by implication DOC) and

Journal of Geophysical Research: Biogeosciences.↗

USGS2023

Snowpack relative permittivity and density derived from near-coincident lidar and ground-penetrating radar

Depth-based and radar-based remote sensing methods (e.g., lidar, synthetic aperture radar) are promising approaches for remotely measuring snow water equivalent (SWE) at high spatial resolution. These approaches require snow density estimates, obtained from in-situ measurements or density models, to calculate SWE. However, in-situ measurements are operationally limited, and few density models have seen extensive evaluation. Here, we combine near-coincident, lidar-measured snow depths with ground-penetrating radar (GPR) two-way travel times ( twt ) of snowpack thickness to derive >20 km of relative permittivity estimates from nine dry and two wet snow surveys at Grand Mesa, Cameron Pass, and Ranch Creek, Colorado. We tested three equations for converting dry snow relative permittivity to snow density and found the Kovacs et al. (1995) equation to yield the best comparison with in-situ measurements (RMSE = 54 kg m −3 ). Variogram analyses revealed a 19 m median correlation length for relative permittivity and snow density in dry snow, which increased to >30 m in wet conditions. We compared derived densities with estimated densities from several empirical models, the Snow Data Assimilation System (SNODAS), and the physically based iSnobal model. Estimated and derived densities were combined with snow depths and twt to evaluate density model performance within SWE remote sensing methods. The Jonas et al. (2009) empirical model yielded the most accurate SWE from lidar snow depths (RMSE = 51 mm), whereas SNODAS yielded the most accurate SWE from GPR twt (RMSE = 41 mm). Densities from both models generated SWE estimates within ±10% of derived SWE when SWE averaged >400 mm, however, model uncertainty increased to >20% when SWE averaged <300 mm. The development and refinement of density models, particularly in lower SWE conditions, is a high priority to fully realize the potential of SWE remote sensing methods.

Colorado↗

USGS2023

Snow surface roughness across spatio-temporal scales

The snow surface is at the interface between the atmosphere and Earth. The surface of the snowpack changes due to its interaction with precipitation, wind, humidity, short- and long-wave radiation, underlying terrain characteristics, and land cover. These connections create a dynamic snow surface that impacts the energy and mass balance of the snowpack, blowing snow potential, and other snowpack processes. Despite this, the snow surface is generally considered a constant parameter in many Earth system models. Data from the National Aeronautics and Space Administration (NASA) Cold Land Processes Experiment (CLPX) collected in 2002 and 2003 across northern Colorado were used to investigate the spatial and temporal variability of snow surface roughness. The random roughness (RR) and fractal dimension (D) metrics used in this investigation are well correlated. However, roughness is not correlated across scales, computed here from snow roughness boards at a millimeter resolution and airborne lidar at a meter resolution. Process scale differences were found based on land cover at each of the two measurement scales, as appraised through measurements in the forest and alpine.

Colorado↗

USGS2023

High resolution SnowModel simulations reveal future elevation-dependent snow loss and earlier, flashier surface water input for the Upper Colorado River Basin

Continued climate warming is reducing seasonal snowpacks in the western United States, where >50% of historical water supplies were snowmelt-derived. In the Upper Colorado River Basin, declining snow water equivalent (SWE) and altered surface water input (SWI, rainfall and snowmelt available to enter the soil) timing and magnitude affect streamflow generation and water availability. To adapt effectively to future conditions, we need to understand current spatiotemporal distributions of SWE and SWI and how they may change in future decades. We developed 100-m SnowModel simulations for water years 2001–2013 and two scenarios: control (CTL) and pseudo-global-warming (PGW). The PGW fraction of precipitation falling as snow was lower relative to CTL, except for November–April at high elevations. PGW peak SWE was lower for low (−45%) and mid elevations (−14%), while the date of peak SWE was uniformly earlier in the year for all elevations (17–23 days). Currently unmonitored high elevation snow represented a greater fraction of total PGW SWE. PGW peak daily SWI was higher for all elevations (30%–42%), while the dates of SWI peaks and centroids were earlier in the year for all elevations under PGW. PGW displayed elevated winter SWI, lower summer SWI, and changes in spring SWI timing were elevation-dependent. Although PGW peak SWI was elevated and earlier compared to CTL, SWI was more evenly distributed throughout the year for PGW. These simulated shifts in the timing and magnitude of SWE and SWI have broad implications for water management in dry, snow-dominated regions.

Colorado↗

USGS2022

Evaluating hydrologic region assignment techniques for ungaged basins in Alaska, USA

Building continental-scale hydrologic models in data-sparse regions requires an understanding of spatial variation in hydrologic processes. Extending these models to ungaged locations requires techniques to group ungaged locations with gaged ones to make process importance and model parameter transfer decisions to ungaged locations. This analysis (1) tested the utility of fundamental streamflow statistics (FDSS) in defining hydrologic regions across Alaska, USA; (2) evaluated if the hydrologic regions represented different hydrologic processes; and (3) tested the ability of random forest and direct assignment techniques, informed by statistically estimated FDSS (FDSSest) and basin characteristics (BCs), to correctly assign ungaged locations to hydrologic regions. Six hydrologic regions were identified across the domain using FDSS. Differences in mean flow, phase shift of the seasonal cycle, and skewness were the primary characteristics defining each region. Two regions represented arctic and continental climates, generally in the northern portion of the domain; four regions represented the southern, maritime portion of the domain. Random forest modeling with BCs (67% success rate) outperformed FDSSest (58% success rate) suggesting that no statistically estimated streamflow was needed to assign ungaged locations to a region. For regions with many sites, most region assignment techniques performed similarly. Random forest modeling performance declined when BCs and FDSSest were both used to predict region membership, suggesting FDSSest had little information in addition to BCs. This analysis demonstrated that FDSS-based hydrologic regions discern process differences across a data-sparse and hydrologically diverse landscape. Process importance rankings from random forest-derived BCs provided model-independent information for making modeling decisions.

Alaska↗

USGS2022

Black carbon dominated dust in recent radiative forcing on Rocky Mountain snowpacks

The vast majority of surface water resources in the semi-arid western United States start as winter snowpack. Solar radiation is a primary driver of snowmelt, making snowpack water resources especially sensitive to even small increases in concentrations of light absorbing particles such as mineral dust and combustion-related black carbon (BC). Here we show, using fresh snow measurements and snowpack modeling at 51 widely distributed sites in the Rocky Mountain region, that BC dominated impurity-driven radiative forcing in 2018. BC contributed three times more radiative forcing on average than dust, and up to 17 times more at individual locations. Evaluation of 2015 to 2018 archived samples from most of the same sites yielded similar results. These findings, together with long-term observations of atmospheric concentrations and atmospheric model studies, indicate that BC rather than dust has dominated radiative forcing by light absorbing impurities on snow for decades, indicating that mitigation strategies to reduce radiative forcing on headwater snow-water resources would need to focus on reducing winter and spring BC emissions.

Colorado, Idaho, Montana, New Mexico, Utah, Wyomin↗

USGS2022

Snow depth retrieval with an autonomous UAV-mounted software-defined radar

We present results from a field campaign to measure seasonal snow depth at Cameron Pass, Colorado, using a synthetic ultrawideband software-defined radar (SDRadar) implemented in commercially available Universal Software Radio Peripheral (USRP) software-defined radio hardware and flown on a small hexacopter unmanned aerial vehicle (UAV). We coherently synthesize an ultrawideband signal from stepped frequency 50-MHz subpulses across 600-2100-MHz frequency bands using a novel nonuniform nonlinear synthetic wideband waveform reconstruction technique that minimizes sweep time and completely eliminates problematic grating lobes and other processing artifacts traditionally seen in stepped waveform synthesis. We image seasonal snow across two transects: a 400-m open Meadow Transect and a 380-m forested transect. We present a surface detection algorithm that fuses data from LiDAR, global navigation satellite system (GNSS)/global positioning system (GPS), and features in the radargram itself to obtain high precision estimates of both snow and ground surface reflections, and thus total snow depth, represented as two-way travel time. The measurements are validated against independent ground-based ground-penetrating radar measurements with correlations coefficients as high as ρ = 0.9 demonstrated. Finally, we compare backscattered radar data collected by the UAV-SDRadar while hovering proximal to a known snow pit with in situ measured snow dielectric profiles and demonstrate imaging of snow stratigraphy.

Colorado↗

USGS2022

Upper Rio Grande Basin water-resource status and trends: Focus area study review and synthesis

The Upper Rio Grande Basin (URGB) is a critical international water resource under pressure from a myriad of climatic, ecological, infrastructural, water-use, and legal constraints. The objective of this study is to provide a comprehensive assessment of the spatial distribution and temporal trends of selected water-budget components (snow processes, evapotranspiration (ET), streamflow processes, and groundwater storage) using integrated analyses, such as watershed modeling and water availability and use data in the URGB over the past three decades. A spatially distributed snow evolution modeling system simulated snowpack processes over 34 years (1984–2017). It highlighted snow water equivalent declines from -35 to -77 mm/decade with widespread variability across elevation zones and land cover types. Gridded actual ET data from the SSEBop model were developed and tested for the URGB and demonstrated that all land-cover types had significant decreasing trends (1986-2015) ranging from -14 to -80 mm/decade. Conductivity-mass-balance (CMB) hydrograph separation results found that baseflow forms a large component of total streamflow, ranging from 29 to 69% (49% average) of total streamflow at 17 URGB sites upstream of Albuquerque, NM. Three of 4 graphical hydrograph separation methods in the U.S. Geological Survey Groundwater Toolbox were found to be inappropriate for estimating baseflow in the URGB; the most promising method, baseflow index (BFI) Standard, was optimized using CMB data and tested at three URGB sites, with resulting overestimation of 0 to 47%. Simulated changes in groundwater storage were extracted from historical and recent groundwater-flow models of select alluvial basins (San Luis, Española, Middle Rio Grande, and Tularosa-Hueco). In general, decreases in groundwater storage were observed from 1903 to 2013 except for the San Luis alluvial basin (Colorado), where periods of recovery are observed. The PRMS hydrologic model was successfully calibrated for 9 near-native subbasins (Nash-Sutcliffe efficiency 0.47 to 0.85) and parameters translated to the remaining subbasins; compared to simulated near-native flows (with minimal influence of reservoirs or diversions), observed Rio Grande streamgage flows demonstrated reductions of 40% or more for New Mexico and Texas areas of the basin. Significant decreasing trends (1980-2015) in precipitation, snowmelt rate, streamflow, and baseflow were observed at many of the 12 streamgage basins studied, which suggests that the decreasing trends for actual ET may be related to overall decreasing water availability in the basin, with negative implications for agricultural production and groundwater abstraction. Water security concerns arise from our findings of higher fraction precipitation as rain, slower snowmelt rates leading to decreasing streamflow production, and an increasing fraction of baseflow, all of which will affect the timing and magnitude of water available for human needs in the basin.

Upper Rio Grande Basin↗

USGS2020

Changes in climate and land cover affect seasonal streamflow forecasts in the Rio Grande headwaters

Seasonal streamflow forecast bias, changes in climate, snowpack, and land cover, and the effects of these changes on relations between basin‐wide snowpack, SNOw TELemetry (SNOTEL) station snowpack, and seasonal streamflow were evaluated in the headwaters of the Rio Grande, Colorado. Results indicate that shifts in the seasonality of precipitation and changing climatology are consistent with periods of overprediction and underprediction in streamflow forecasts. Multiple linear regression of SNOTEL data, postcedent precipitation, and land‐cover changes explained 2%–18% more variability in streamflow prediction than using SNOTEL station data alone. Simulated basin‐wide snowpack from a physically based model had significant negative trends in snow water equivalent (−4.33 mm/yr) and snow‐covered area (−0.05%/yr) during the melt period April–June. Simulated streamflow from a precipitation‐runoff model increased an average 5% when the effects of bark beetle‐induced tree mortality were compared to a baseline simulation with static vegetation. The effects of a 2013 wildfire increased simulated seasonal streamflow an average 35% for 1–4 years postfire. The combined effects of climate and land‐cover changes on snowpack‐streamflow relations highlight the difficulty in seasonal streamflow forecasting, which has important implications for water‐resource management.

Colorado↗

USGS2020

Runoff sensitivity to snow depletion curve representation within a continental scale hydrologic model

The spatial variability of snow water equivalent (SWE) can exert a strong influence on the timing and magnitude of snowmelt delivery to a watershed. Therefore, the representation of subgrid or subwatershed snow variability in hydrologic models is important for accurately simulating snowmelt dynamics and runoff response. The U.S. Geological Survey National Hydrologic Model infrastructure with the Precipitation-Runoff Modeling System (NHM-PRMS) represents the subgrid variability of SWE with snow depletion curves (SDCs), which relate snow-covered area to watershed-average SWE during the snowmelt period. The main objective of this research was to evaluate the sensitivity of simulated runoff to SDC representation within the NHM-PRMS across the continental United States (CONUS). SDCs for the model experiment were derived assuming a range of SWE coefficient of variation (CV) values and a lognormal probability distribution function. The NHM-PRMS was simulated at a daily time step for each SDC over a 14-year period. Results highlight that increasing the subgrid snow variability (by changing the SDC) resulted in a consistently slower snowmelt rate and longer snowmelt duration when averaged across the hydrologic response unit scale. Simulated runoff was also found to be sensitive to SDC representation, as increases in the subgrid SWE CV by 1.0 resulted in decreases in runoff ratio by as much as 12 percent in snow-dominated regions of the CONUS. Simulated decreases in runoff associated with slower snowmelt rates were approximately inversely proportional to increases in simulated evapotranspiration. High snow persistence and peak SWE:annual precipitation combined with a water limited dryness index were associated with the greatest runoff sensitivity to changing snowmelt. Results from this study highlight the importance of carefully parameterizing SDCs for hydrologic modeling. Furthermore, improving model representation of snowmelt input variability and its relation to runoff generation processes is shown to be an important consideration for future modeling applications.

Hydrological Processes↗

USGS2020

A new sampler for the collection and retrieval of dry dust deposition

Atmospheric dust can influence biogeochemical cycles, accelerate snowmelt, and affect air, water quality, and human health. Yet, the bulk of atmospherically transported material remains poorly quantified in terms of total mass fluxes and composition. This lack of information stems in part from the challenges associated with measuring dust deposition. Here we report on the design and efficacy of a new dry deposition sampler (Dry Deposition Sampling Unit (DSU)) and method that quantifies the gravitational flux of dust particles. The sampler can be used alone or within existing networks such as those employed by the National Atmospheric Deposition Program (NADP). Because the samplers are deployed sterile and the use of water to remove trapped dust is not required, this method allows for the recovery of unaltered dry material suitable for subsequent chemical and microbiological analyses. The samplers were tested in the laboratory and at 15 field sites in the western United States. With respect to material retention, sampler performance far exceeded commonly used methods. Retrieval efficiency was >97% in all trials and the sampler effectively preserved grain size distributions during wind exposure experiments. Field tests indicated favorable comparisons to dust-on-snow measurement across sites ( r 2 0.70, p < 0.05) and within sites to co-located aerosol data ( r 2 0.57–0.99, p < 0.05). The inclusion of dust deposition and composition monitoring into existing networks increases spatial and temporal understanding of the atmospheric transport on materials and substantively furthers knowledge of the effects of dust on terrestrial ecosystems and human exposure to dust and associated deleterious compounds.

Arizona, California, Colorado, Idaho, Nevada, Utah↗

USGS2020

Relating hydroclimatic change to streamflow, baseflow, and hydrologic partitioning in the Upper Rio Grande Basin, 1980 to 2015

Understanding how changing climatic conditions affect streamflow volume and timing is critical for effective water management. In the Rio Grande Basin of the southwest U.S., decreasing snowpack, increasing minimum temperatures, and decreasing streamflow have been observed in recent decades, but the effects of hydroclimatic changes on baseflow, or groundwater discharge to streams, have not been investigated. In this study, we determine how trends in precipitation, snowpack accumulation, and snowmelt rate relate to total streamflow, baseflow, and the hydrologic partitioning of baseflow and runoff at 12 sites in the Upper Rio Grande Basin (URGB) during 1980 to 2015. Total streamflow was partitioned into baseflow and runoff components at a daily time step using conductivity-mass-balance hydrograph separation. Trends in annual total streamflow, baseflow, runoff, baseflow index, precipitation, snowmelt rate, and peak snow water equivalent (SWE) were evaluated from 1980 to 2015 using the non-parametric Mann-Kendall trend test. Results indicate that baseflow forms a large component of total streamflow, contributing an average of 49% of total discharge upstream of Albuquerque, NM. During 1980 to 2015, decreasing trends in total streamflow occurred at 9 of 12 sites and were almost always associated with decreases in baseflow, suggesting that baseflow volumes can respond to changing climatic and anthropogenic conditions within decades. Decreasing snowmelt rates were more frequently associated with decreases in baseflow and total streamflow than were decreases in precipitation and peak SWE, highlighting the importance of snowmelt rate as a process controlling streamflow generation. If snow accumulation and snowmelt rates continue to decrease in the future, results indicate that total streamflow and baseflow volumes will decline, and that baseflow will become a larger fraction of total streamflow in the URGB.

Colorado, New Mexico↗

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↗