Search USGS⌕ Search

Geology topics

Research about Colorado, Wyoming

Source-linked reports with geographic coverage including Colorado, Wyoming.

At least 19 recordsLinked to original sources

Modeling observational uncertainty improves estimation of breeding-season parameters that inform lek-based monitoring for Centrocercus urophasianus (Greater Sage-Grouse)

Counts of males on leks are widely used as abundance indices to monitor populations of lek-breeding species, but they can be influenced by variation in male attendance, inter-lek movement, and mortality. When there is imperfect observability of males on leks, estimates of lek attendance and inter-lek movement are likely to be biased low. We compared traditional multi-state mark–recapture (MSMR) and MSMR models with state uncertainty (MSMR-SU) to determine how accounting for uncertainty in GPS-derived lek attendance improves inference of rates underlying lek-based abundance indices for male Centrocercus urophasianus (Greater Sage-Grouse). We found that modeling observational uncertainty with MSMR-SU models consistently reduced bias in lek attendance and inter-lek movement estimates for adult and yearling males compared to MSMR models. Estimates of daily maximum male lek attendance (adults: 0.79; yearlings: 0.68) and inter-lek movement rates (adults: 0.10; yearlings: 0.07) from MSMR-SU models were higher than daily maximum lek attendance (adults: 0.74; yearlings: 0.55) and inter-lek movement rates (adults: 0.01; yearlings: 0.04) from MSMR models. Consistent with previous studies, adults exhibited higher lek attendance and mortality rates than yearlings. However, unlike earlier studies, inter-lek movement rates were similar between age classes, with adults tending to move earlier in the lekking season than yearlings. Wind speed and precipitation had a weak negative effect on lek attendance. Greater distance to the nearest known lek strongly reduced inter-lek movement. Accounting for lekking-season mortality shifted the estimated date of maximum lek attendance earlier in the season by 8–9 days for adults and 1–3 days for yearlings. Accounting for state uncertainty due to error in state assignment should reduce bias in demographic estimates for use in wildlife conservation and management decisions.

Colorado, Wyoming↗

Ecology of reintroduced Rocky Mountain bighorn sheep in Dinosaur National Monument

Translocations have been widely used to restore and conserve bighorn sheep ( Ovis canadensis ) populations in North America. Some translocations have been successful, but many populations remain small and genetically isolated. Population structure can influence the viability and long-term success of reintroductions. Social ungulates often function as interconnected subpopulations (metapopulations); however, few studies evaluate subpopulation sizes, connectivity, and genetic diversity within metapopulations. To address this gap, we conducted a comprehensive study of a reintroduced Rocky Mountain bighorn sheep ( Ovis canadensis canadensis ) population in Dinosaur National Monument in Colorado and Utah, USA, between 2006–2020. We analyzed global positioning system (GPS) radio-collar data, genetic samples, and results of health testing to evaluate abundance, distribution, genetic structure and diversity, habitat use, movement and connectivity, and presence of or exposure to respiratory pathogens. We integrated these analyses to evaluate the outcomes of a reintroduction effort that began in 1952, over 70 years ago, and to inform management decisions in Dinosaur National Monument. We also provide a framework for evaluating metapopulation processes, including a non-invasive approach that links genetic structure with Bayesian spatial capture-recapture analyses to estimate subpopulation sizes. Despite models indicating continuous suitable habitat, we found a spatially structured population with at least 4 subpopulations with constrained connectivity. Evidence from step selection and density analyses suggested that movement among subpopulations may be limited by semi-permeable barriers including rivers and human disturbance, which could contribute to maintenance of spatial structure over time. In 2006, antibody to Mycoplasma ovipneumoniae was detected in all geographically and genetically distinct subpopulations. Widespread clinical signs of disease and a confirmed exposure to M. ovipneumoniae in 2019 indicate a long-term disease challenge. Proximity to domestic sheep creates repeated opportunities for introduction of new M. ovipneumoniae strains. We estimated abundance in 2019 at 109 (95% CrI = 87–133), composed of subpopulations ranging from 18–39 animals (95% CrIs from 11–50). Genetic diversity was relatively high compared to other reintroduced and native Rocky Mountain bighorn sheep populations, which is likely a consequence of multiple translocations from different sources. Three of 4 subpopulation centers generally aligned with the locations of original translocation release sites. Persistence in the presence of pathogens may be facilitated by metapopulation structure and moderately high genetic diversity. Conversely, metapopulation structure can also facilitate pathogen persistence. Our approach offers a path to advance understanding of the population ecology of reintroduced bighorn sheep and can inform effective conservation and management of their populations.

Colorado, Wyoming↗

Integrated stratigraphic and geochemical analysis of organic-rich intervals of the Lewis Shale in the eastern Washakie Basin, Wyoming

Geological studies in the Cretaceous Western Interior Seaway (KWIS) in North America reveal highly variable sedimentological conditions on its western shore caused by rapidly changing sea level and detrital input during the seaway’s closure. Here we examine a 601-foot (183 meters) continuous core through the lower part of the Maastrichtian Lewis Shale in the eastern Washakie Basin, Wyoming, through integrating stratigraphic and geochemical analyses to better understand organic matter deposition and preservation during the final marine transgression within the seaway. The core penetrates eight organic-rich flooding intervals (F1-F8) and a regional condensed section, informally known as the Asquith marker. The lower portions of the core record sediment input from a southern source, likely the Sierra Madre/Park Ranges, while the upper part records sediment input from a northern source, likely the Granite Mountains. This provenance transition is supported by regional stratigraphic analysis and changes in bulk mineral and trace metal composition. The Asquith marker and early flooding surfaces are enriched in oil-prone, marine organic matter deposited under dysoxic to euxinic conditions, whereas younger flooding surfaces show increased terrigenous input and poorer preservation conditions. The Asquith marker is identified as a prime oil-prone source rock. Although younger flooding surfaces also exhibit favorable source-rock properties, their generative potential is reduced due to increased clastic and terrigenous organic matter deposition caused by regression. These results emphasize the importance of sediment source variability, organic matter preservation, and changing redox conditions to provide insight into the sediment provenance and petroleum potential of the Lewis Shale in the context of the final transgression within the KWIS.

Colorado, Wyoming↗

Quantifying intraspecific variation in host resistance ad tolerance to a lethal pathogen

Testing for intraspecific variation for host tolerance or resistance in wild populations is important for informing conservation decisions about captive breeding, translocation, and disease treatment. Here, we test the importance of tolerance and resistance in multiple populations of boreal toads ( Anaxyrus boreas boreas ) against Batrachochytrium dendrobatidis (Bd), the amphibian fungal pathogen responsible for the greatest host biodiversity loss due to disease. Boreal toads have severely declined in Colorado (CO) due to Bd, but toad populations challenged with Bd in western Wyoming (WY) appear to be less affected. We used a common garden infection experiment to expose post-metamorphic toads sourced from four populations (2 in CO and 2 in WY) to Bd and monitored changes in mass, pathogen burden and survival for 8 weeks. We used a multi-state modelling approach to estimate weekly survival and transition probabilities between infected and cleared states, reflecting a dynamic infection process that traditional approaches fail to capture. We found that WY boreal toads are more tolerant to Bd infection with higher survival probabilities than those in CO when infected with identical pathogen burdens. WY toads also appeared more resistant to Bd with a higher probability of infection clearance and an average of 5 days longer to reach peak infection burdens. Our results demonstrate strong intraspecific differences in tolerance and resistance that likely contribute to why population declines vary regionally across this species. Our multi-state framework allowed us to gain inference on typically hidden disease processes when testing for host tolerance or resistance. Our findings demonstrate that describing an entire host species as ‘tolerant’ or ‘resistant’ (or lack thereof) is unwise without testing for intraspecific variation.

Colorado, Wyoming↗

Stratigraphic cross sections of the Lewis Shale in the eastern part of the southwestern Wyoming Province, Wyoming and Colorado

Three stratigraphic cross sections A–A', B–B', and C–C' were created for the Lewis Shale and associated strata in the eastern part of the Southwestern Wyoming Province of Wyoming and Colorado. The cross sections highlight 15 clinothems within the Lewis Shale, Fox Hills Sandstone, and Lance Formation progradational system (also referred to as the Lewis Shale system). Additionally, the cross sections indicate that multiple source areas were active at the same time during deposition of the Lewis Shale. Specifically, the north-to-southeast cross section A–A' demonstrates that the northern, sand-rich source was being deposited and onlapping onto an older, southern, mud-rich source.

Colorado, Wyoming↗

Testing the potential of streamflow data to predict spring migration of an ungulate herds

In mountainous and high latitude regions, migratory animals exploit green waves of emerging vegetation coinciding with rising daily mean temperatures initiating snowmelt across the landscape. Snowmelt also causes rivers and streams draining these regions to swell, a process referred to as to as the ‘spring pulse.’ Networks of streamgages measuring streamflow in these regions often have long-term and continuous periods of record available in real-time and at the daily time step, and thus produce data with potential to predict temporal migration patterns for species exploiting green waves. We tested the potential of models informed by streamflow data to predict timing of spring migration of mule deer ( Odocoileus hemionus ) herds in a headwater basin of the Colorado River. Models using streamflow data were compared with those informed by traditional temperature-derived measures of the onset of spring. Non-parametric linear-regression techniques were used to test for temporal stationarity in each variable, and logistic-regression models were used to produce probabilities of migration initiation. Our analysis indicates that models using daily streamflow data can perform as well as those using temperature-derived data to predict past-migration patterns, and nearly as well in potential to forecast future migrations. The best performing model was used to generate probabilities of onset of migration for mule deer herds over the 69-year period-of-record from a streamgage. That model indicated spring migration has been trending toward earlier initiations, with modeled median initiations shifting from a Julian day of 123 in the mid 20 th century to Julian day 115 over the most recent two decades. The period of 1960 to 1979 had the latest modeled median initiations with Julian day of 128. The analyses demonstrate promise for merging existing hydrologic and biological data collection platforms in these regions to explore timing of past migration patterns and predict migration onsets in real-time.

Colorado, Wyoming↗

Fingerprinting historical tributary contributions to floodplain sediment using bulk geochemistry

Sediment deposition on floodplains is essential for the development and maintenance of riparian ecosystems. Upstream erosion is known to influence downstream floodplain construction, but linking these disparate processes is challenging, especially over large spatial and temporal scales. Sediment fingerprinting is thus a robust tool to establish process linkages between downstream floodplain development and sediment production in distal headwater basins. Here we use sediment geochemistry to connect historical erosion in several tributaries of the Yampa River in Colorado and Wyoming, USA, to the construction of downstream floodplains on which extensive cottonwood forests established. Using a combination of conventional techniques and the relatively novel machine-learning random forest algorithm, we build multiple fingerprints of diagnostic geochemical tracers that are then input into a Bayesian mixing model to apportion provenance of floodplain sediment. Sediment samples for provenance analysis were collected from an excavated floodplain in Deerlodge Park on the Yampa River at the rooting surface of the surrounding cottonwood forest and dominantly comprised of very fine (4Φ) sand. Fingerprinting analysis of the 4Φ fraction of collected floodplain sink (n = 38) and tributary source (n = 218) samples revealed floodplain sediment to be dominantly sourced from the tributaries of Muddy Creek (45 ± 4%) and Sand Wash (42 ± 6%). Dendrochronology results moreover indicate the Deerlodge floodplain sediment was deposited in ∼1912, which falls squarely within the time (1880–1940) these tributaries were actively eroding. Taken together, study results indicate a demonstrable link between historical tributary erosion and downstream floodplain construction and concomitant forest establishment. Our findings suggest processes operating in tributary watersheds play an important role in the dynamics of large rivers and emphasize both the need for holistic, collaborative management of sediment as an essential resource and the potential to utilize sediment fingerprinting to inform and direct river ecosystem management.

Colorado, Wyoming↗

Sex-specific behaviors of hunted mule deer during rifle season

Animal populations face increased threats to mobility and access to critical habitat from a variety of human disturbances including roads, residential development, agriculture, and energy development. Disturbance from human hunting is known to alter habitat use in ungulates, but recent work suggests that hunting may also trigger the onset of migration. Whether this holds true across ungulate species and hunting systems warrants further empirical testing. We used global positioning system location data from mule deer ( Odocoileus hemionus ) in south-central Wyoming, USA, to evaluate the sex-specific effects of hunting on habitat selection and migratory behavior from 2016 to 2018. We modeled habitat selection before and during hunting season using a step selection function, and we used time-to-event models to evaluate if hunting triggered migration. We found habitat selection and migration timing to be sex specific. Males responded to hunting season by selecting security habitat away from motorized routes, whereas females used habitat through hunting season that retained higher forage quality. Weather, as indexed by temperature and precipitation (i.e., snowfall), influenced migration timing for males and females. Migration timing in males was influenced by migration distance, where individuals traveling >50 km tended to migrate earlier than individuals moving <50 km. For deer that survived to rifle season, hunting was less influential on migration timing than environmental factors. Rifle season increased the likelihood of migration by 2% in females and <0.01% in males compared to outside rifle season. Our findings suggest that roadless areas on mule deer summer ranges and within migration corridors reduce the effects of hunting disturbance. Consequently, managers may consider limiting the use of motorized vehicles as a method for reducing effects on migration from hunting disturbance.

Colorado, Wyoming↗

Spatial memory shapes migration and its benefits: Evidence from a large herbivore

From fine-scale foraging to broad-scale migration, animal movement is shaped by the distribution of resources. There is mounting evidence, however, that learning and memory also guide movement. Although migratory mammals commonly track resource waves, how resource tracking and memory guide long-distance migration has not been reconciled. We examined these hypotheses using movement data from four populations of migratory mule deer ( n = 91). Spatial memory had an extraordinary influence on migration, affecting movement 2–28 times more strongly than tracking spring green-up or autumn snow depth. Importantly, with only an ability to track resources, simulated deer were unable to recreate empirical migratory routes. In contrast, simulated deer with memory of empirical routes used those routes and obtained higher foraging benefits. For migratory terrestrial mammals, spatial memory provides knowledge of where seasonal ranges and migratory routes exist, whereas resource tracking determines when to beneficially move within those areas.

Colorado, Wyoming↗

Using social-context matching to improve spatial function-transfer performance for cultural ecosystem service models

Recreational and aesthetic enjoyment of public lands is increasing across a wide range of activities, highlighting the need to assess and adapt management to accommodate these uses. Despite a growing number of studies on mapping cultural ecosystem services, most are local-scale assessments that rely on costly and time-consuming primary data collection. As a result, the availability of spatial information on non-market values associated with cultural ecosystem services (social values) remains limited. Spatial function transfer, if it could be justified for social-value models, would expedite the development of social-value information and promote its more regular inclusion in ecosystem service assessments. We used survey data from six national forests in Colorado and Wyoming to explore the potential for transferring cultural ecosystem service models between forests and specifically to test the hypothesis that transfer performance increases with social-context similarity between transferring and receiving areas. Results confirm this relationship but fall just short of being able to predict with certainty when transferred models will meet the minimum performance criterion needed for defensible use by managers. Social values are highly variable and can be difficult to predict, but our results suggest that with the right combination of indicators spatial function transfer can become a defensible means of generating social-value information when primary data collection is not feasible.

Colorado, Wyoming↗

Life-history variation of two inland salmonids revealed through otolith microchemistry analysis

Increasingly, otolith microchemistry analysis is used as a tool to trace fish migrations, especially migrations of diadromous fishes. Yet, few studies have used otolith microchemistry to trace migrations in small inland watersheds, leaving major knowledge gaps in our understanding of inland fish spatial ecology. Here, we evaluate the use of tributary habitat for spawning and describe and compare fluvial brown trout ( Salmo trutta ) and rainbow trout ( Oncorhynchus mykiss ) natal origin distribution, time spent in natal streams, and spawning site fidelity. 63% of rainbow trout and 57% of brown trout migrated after hatching. Brown trout showed greater variation in time spent in natal tributaries, suggesting that individuals are temporally distributing risk among offspring. By contrast, rainbow trout showed greater variation in natal origin, suggesting that individuals are spatially distributing risk among offspring. Our results indicate there is high inter- and intraspecific migration variation in inland salmonid populations, which may be linked to access to a mosaic of spawning and rearing habitat types.

Colorado, Wyoming↗

Quantifying 87Sr/86Sr temporal stability and spatial heterogeneity for use in tracking fish movement

The specificity and accuracy of inferred fish origin and movement relies on describing spatial heterogeneity and temporal stability of environmental signatures. But the cost and logistics of sample collection often precludes the complete quantification of environmental signature temporal stability and spatial heterogeneity. We used repeated sampling and a novel approach (Bayesian ridge regression, BRR) to quantify the temporal stability and spatial heterogeneity of 87 Sr/ 86 Sr, respectively. We explained 86% of observed variation in 87 Sr/ 86 Sr using a BRR model and estimated 87 Sr/ 86 Sr throughout the Upper North Platte River Basin with high accuracy (±0.00106). Year to year variation in 87 Sr/ 86 Sr signatures ranged from 0.00007 to 0.00073 (SD), while seasonal variation ranged from 0.00091 to 0.00134 (SD). We then assessed the specificity and discussed the accuracy of inferring movement using three scenarios of described spatial heterogeneity. Our results indicate reliable inference of fish movement requires comprehensive quantification of spatial heterogeneity and temporal variation in environmental signatures.

Colorado, Wyoming↗

Geologic map of the Fort Collins 30'×60' quadrangle, Larimer and Jackson Counties, Colorado, and Albany and Laramie Counties, Wyoming

The rocks and landforms of the Fort Collins 30 ′ × 60 ′ 1:100,000-scale U.S. Geological Survey quadrangle reveals a particularly complete record of geologic history in the northern Front Range of Colorado. The Proterozoic basement rocks exposed in the core of the range preserve evidence of Paleoproterozoic marine sedimentation, volcanism, and regional soft-sediment deformation, followed by regional folding and gradational metamorphism. Mesoproterozoic time was marked by intrusion of the Berthoud Plutonic Suite into crust that was structurally neutral or moderately extending in an east-northeast direction. Evidence of the late Paleozoic Anasazi uplift (Ancestral Rocky Mountains uplift) within the quadrangle is recorded by removal of Permian and older sediments and deposition of proximal Pennsylvanian and Permian strata unconformably onto the exhumed Proterozoic basement rocks. The Phanerozoic sediments indicate a steady progression of fluvial, eolian, and lacustrine environments throughout most of the Mesozoic Era which was a time of relatively slow sediment accumulation. Early Cretaceous time was marked by incursion of the Cretaceous Western Interior Seaway, a shallow-water marine embayment that persisted throughout the latter part of the Mesozoic Era. Sedimentation rates increased significantly in the latter part of this period during down-warping related to distant crustal loading by thrusting along the western continental margin. With onset of the Laramide orogeny in latest Cretaceous time, mountain building resumed in this region. This deformation placed Proterozoic rock over Cretaceous and Paleocene strata along the western margin of the Front Range and Medicine Bow Mountains. Post-Laramide time was marked by a prolonged period of weathering, erosion, and planation of the basement-rock surface, extending perhaps into late Oligocene or early Miocene time. Erosion on the eastern slope of the Front Range in late Paleogene to early Neogene time produced a broad, rolling surface surrounding residual highlands and east-trending fluvial channels filled with coarse, boulder gravel. Significant global cooling during the Pliocene led to glaciation during the Quaternary. In the Rocky Mountain region, renewed uplift allowed erosion to accentuate the topographic relief across the high mountains of the map area and established the elevations necessary to trigger accumulation of persistent snow and ice. Mountain glaciers advanced and retreated during at least three glacial-interglacial cycles during the middle and late Pleistocene in this area. Erosion continues to this day on the High Plains east of the mountain front, and progressive incision of the drainage is recorded by at least five major gravel-clad terrace and pediment surfaces along the major fluvial channels that connect to the South Platte River system.

Colorado, Wyoming↗

Estimating soil respiration in a subalpine landscape using point, terrain, climate and greenness data

Landscape carbon (C) flux estimates are necessary for assessing the ability of terrestrial ecosystems to buffer further increases in anthropogenic carbon dioxide (CO2) emissions. Advances in remote sensing have allowed for coarse-scale estimates of gross primary productivity (GPP) (e.g., MODIS 17), yet efforts to assess spatial patterns in respiration lag behind those of GPP. Here, we demonstrate a method to predict growing season soil respiration at a regional scale in a forested ecosystem. We related field measurements (n=144) of growing season soil respiration across subalpine forests in the Southern Rocky Mountains ecoregion to a suite of biophysical predictors with a Random Forest model (30 m pixel size). We found that Landsat Enhanced Vegetation Index (EVI), growing season AI, temperature, precipitation, elevation, and slope aspect explained spatiotemporal variability in soil respiration. Our model had a psuedo-r2 of 0.45 and root mean squared error (RMSE) of roughly one-quarter of the mean value of respiration. Predicted growing season soil respiration across the region was remarkably consistent across 2004, 2005 and 2006 (150-d averages of 542.8, 544.3, and 536.5 g C m-2, respectively). Yet, we observed substantial variability in spatial patterns of soil respiration predictions that varied between years, suggesting that our method is sensitive to changes in respiration drivers. We compared our estimates to MODIS GPP and nocturnal net ecosystem exchange (NEE) derived from eddy covariance towers as a proxy for ecosystem respiration. Averaged across the predictive region, mean predicted growing season soil respiration was 73% of MODIS GPP, while predicted soil respiration was generally within 20% of nocturnal NEE from eddy covariance towers. This study demonstrated that geospatial and remotely-sensed datasets can be used in a statistical modeling framework to estimate soil respiration at landscape scales.

Colorado, Wyoming↗

Integrating physiological stress into the movement ecology of migratory ungulates: A spatial analysis with mule deer

Rapid climate and human land-use change may limit the ability of long-distance migratory herbivores to optimally track or “surf” high-quality forage during spring green-up. Understanding how anthropogenic and environmental stressors influence migratory movements is of critical importance because of their potential to cause a mismatch between the timing of animal movements and the emergence of high-quality forage. We measured stress hormones (fecal glucocorticoid metabolites; FGMs) to test hypotheses about the effects of high-quality forage tracking, human land use, and use of stopover sites on the physiological state of individuals along a migratory route. We collected and analyzed FGM concentrations from 399 mule deer (Odocoileus hemionus) samples obtained along a 241 km migratory route in western Wyoming, USA, during spring 2015 and 2016. In support of a fitness benefit hypothesis, individuals occupying areas closer to peak forage quality had decreased FGM levels. Specifically, for every 10-day interval closer to peak forage quality, we observed a 7% decrease in FGMs. Additionally, we observed support for both an additive anthropogenic stress hypothesis and a hypothesis that stopovers act as physiological refugia, wherein individuals sampled far from stopover sites exhibited 341% higher FGM levels if in areas of low landscape integrity compared to areas of high landscape integrity. Overall, our findings indicate that the physiological state of mule deer during migration is influenced by both anthropogenic disturbances and their ability to track high-quality forage. The availability of stopovers, however, modulates physiological responses to those stressors. Thus, our results support a recent call for the prioritization of stopover locations and connectivity between those locations in conservation planning for migratory large herbivores.

Colorado, Wyoming↗

Productivity and CO 2 exchange of Great Plains ecoregions. I. Shortgrass steppe: Flux tower estimates

The shortgrass steppe (SGS) occupies the southwestern part of the Great Plains. Half of the land is cultivated, but significant areas remain under natural vegetation. Despite previous studies of the SGS carbon cycle, not all aspects have been completely addressed, including gross productivity, ecosystem respiration, and ecophysiological parameters. Our analysis of 1998 − 2007 flux tower measurements at five Bowen ratio–energy balance (BREB) and three eddy covariance (EC) sites characterized seasonal and interannual variability of gross photosynthesis and ecosystem respiration. Identification of the nonrectangular hyperbolic equation for the diurnal CO 2 exchange, with vapor pressure deficit (VPD) limitation and exponential temperature response, quantified quantum yield α , photosynthetic capacity A max , and respiration rate r d with variation ranges (19 \< α \ < 51 mmol mol − 1 , 0.48 \< A max \ < 2.1 mg CO 2 m − 2 s − 1 , 0.15 \< r d \ < 0.49 mg CO 2 m − 2 s − 1 ). Gross photosynthesis varied from 1 100 to 2 700 g CO 2 m − 2 yr − 1 , respiration from 900 to 3,000 g CO 2 m − 2 yr − 1 , and net ecosystem production from − 900 to + 700 g CO 2 m − 2 yr − 1 , indicating that SGS may switch from a sink to a source depending on weather. Comparison of the 2004 − 2006 measurements at two BREB and two parallel EC flux towers located at comparable SGS sites showed moderately higher photosynthesis, lower respiration, and higher net production at the BREB than EC sites. However, the difference was not related only to methodologies, as the normalized difference vegetation index at the BREB sites was higher than at the EC sites. Overall magnitudes and seasonal patterns at the BREB and the EC sites during the 3-yr period were similar, with trajectories within the ± 1.5 standard deviation around the mean of the four sites and mostly reflecting the effects of meteorology.

Colorado, Wyoming↗