Geology Reports⌕ Search

Geology topics

Research about Colorado

Source-linked reports with geographic coverage including Colorado.

At least 127 records · Page 7Linked to original sources

Trajectories and tipping points of piñon–juniper woodlands after fire and thinning

Piñon–juniper (PJ) woodlands are a dominant community type across the Intermountain West, comprising over a million acres and experiencing critical effects from increasing wildfire. Large PJ mortality and regeneration failure after catastrophic wildfire have elevated concerns about the long-term viability of PJ woodlands. Thinning is increasingly used to safeguard forests from fire and in an attempt to increase climate resilience. We have only a limited understanding of how fire and thinning will affect the structure and function of PJ ecosystems. Here, we examined vegetation structure, microclimate conditions, and PJ regeneration dynamics following ~20 years post-fire and thinning treatments. We found that burned areas had undergone a state shift that did not show signs of returning to their previous state. This shift was characterized by (1) distinct plant community composition dominated by grasses; (2) a lack of PJ recruitment; (3) a decrease in the sizes of interspaces in between plants; (4) lower abundance of late successional biological soil crusts; (5) lower mean and minimum daily soil moisture values; (6) lower minimum daily vapor pressure deficit; and (7) higher photosynthetically active radiation. Thinning created distinct plant communities and served as an intermediate between intact and burned communities. More intensive thinning decreased PJ recruitment and late successional biocrust cover. Our results indicate that fire has the potential to create drier and more stressful microsite conditions, and that, in the absence of active management following fire, there may be shifts to persistent ecological states dominated by grasses. Additionally, more intensive thinning had a larger impact on community structure and recruitment than less intensive thinning, suggesting that careful consideration of goals could help avoid unintended consequences. While our results indicate the vulnerability of PJ ecosystems to fire, they also highlight management actions that could be adapted to create conditions that promote PJ re-establishment.

Colorado↗

Geoelectric evidence for a wide spatial footprint of active extension in central Colorado

Three-dimensional magnetotelluric (MT) imaging in central Colorado reveals a set of north-striking high-conductivity tracks at lower-crustal (50–20 km) depths, with conductive finger-like structures rising off these tracks into the middle crust (20–5 km depth). We interpret these features to represent saline aqueous fluids and partial melt that are products of active extensional tectonomagmatism. These conductors are distributed over a wider region than the narrow corridor along which Rio Grande rift structures are traditionally mapped at the surface, and they consequently demarcate regions of the lower crust where accommodation of bulk extensional strain has concentrated conductive phases. Our observations reveal limitations in existing models of Rio Grande rift activity and may reflect unrecognized spatiotemporal variations in rift system evolution globally.

Colorado↗

Resurvey of cross sections on the Yampa and Little Snake Rivers in Lily and Deerlodge Parks, Colorado

Resurveys of seven geomorphologic cross sections located in the Lily Park and Deerlodge Park, Colorado, reaches of the Yampa and Little Snake Rivers were conducted in October 2017. These cross sections extend from Lily Park, at the confluence of the two rivers, to Deerlodge Park within Dinosaur National Monument. Four cross sections were first surveyed in 1983 and then resurveyed in 1997. The remaining three cross sections were first surveyed in 1997. Analysis of historical aerial photographs (taken from 1961 to 2015) was conducted to contextualize the measured changes in the cross sections, confirm cross-section longitudinal positions along the rivers, and verify the timing of artificial realignment and straightening of the Little Snake River. Erosion occurred between 1983 and 1997 in all four cross sections first surveyed in 1983, largely through channel widening. Continued erosion occurred between 1997 and 2017 in six of the seven cross sections, also largely by channel widening with only minor changes in channel depth. Though erosion occurred over a longer time period, the net erosion observed at these cross sections over three decades is consistent with the net erosion documented by a sediment-transport-based monitoring program on the Yampa River and Little Snake Rivers from 2013 to 2020.

Colorado↗

Rethinking cost-share programs in consideration of economic equity: A case study of wildfire risk mitigation assistance for private landowners

Public agencies and organizations often deliver financial assistance through cost sharing, in which recipients contribute some portion toward total costs. However, cost sharing might raise equity concerns if it reduces participation among populations with lower incomes. Here, we revisit a past study using a richer dataset (n=1,689) to assess whether stated income levels affect survey respondents' willingness to participate in a cost share program for vegetation reduction to mitigate wildfire risk in western Colorado. Results show that residents with lower incomes are less likely to participate even though they can choose to contribute 0% toward a cost share. Residents reporting incomes less than $50,000 are 11 percentage points less likely to participate than those reporting incomes of $200,000 or more. They also are willing to pay a lower share (26 percentage points less) if they do participate. Results indicate potential economic equity concerns from the use of such programs.

Colorado↗

Use of Doppler velocity radars to monitor and predict debris and flood wave velocities and travel times in post-wildfire basins

The magnitude and timing of extreme events such as debris and floodflows (collectively referred to as floodflows) in post-wildfire basins are difficult to measure and are even more difficult to predict. To address this challenge, a sensor ensemble consisting of noncontact, ground-based (near-field), Doppler velocity (velocity) and pulsed (stage or gage height) radars, rain gages, and a redundant radio communication network was leveraged to monitor flood wave velocities, to validate travel times, and to compliment observations from NEXRAD weather radar. The sensor ensemble (DEbris and Floodflow Early warNing System, DEFENS) was deployed in Waldo Canyon, Pike National Forest, Colorado, USA, which was burned entirely (100 percent burned) by the Waldo Canyon fire during the summer of 2012 ( MTBS, 2020 ). Surface velocity, stage, and precipitation time series collected during the DEFENS deployment on 10 August 2015 were used to monitor and predict flood wave velocities and travel times as a function of stream discharge (discharge; streamflow). The 10 August 2015 event exhibited spatial and temporal variations in rainfall intensity and duration that resulted in a discharge equal to 5.01 cubic meters per second (m 3 /s). Discharge was estimated post-event using a slope-conveyance indirect discharge method and was verified using velocity radars and the probability concept algorithm. Mean flood wave velocities – represented by the kinematic celerity ( c k = 2.619 m e t e r s p e r s e c o n d , m / s ± 0.556 p e r c e n t "> 𝑐 𝑘 = 2.619 𝑚𝑒𝑡𝑒𝑟𝑠 𝑝𝑒𝑟 𝑠𝑒𝑐𝑜𝑛𝑑, m/s ± 0.556 𝑝𝑒𝑟𝑐𝑒𝑛𝑡) and dynamic celerity ( c d = 3.533 m / s ± 0.181 p e r c e n t a n d t h e i r u n c e r t a i n t i e s "> 𝑐 𝑑 = 3.533 m/s ± 0.181 𝑝𝑒𝑟𝑐𝑒𝑛𝑡) and their uncertainties were computed. L-moments were computed to establish probability density functions (PDFs) and associated statistics for each of the at-a-section hydraulic parameters to serve as a workflow for implementing alert networks in hydrologically similar basins that lack data. Measured flood wave velocities and travel times agreed well with predicted values. Absolute percent differences between predicted and measured flood wave velocities ranged from 1.6 percent to 49 percent and varied with water slope, hydraulic radius, and depth. The kinematic celerity was a better predictor for steep slopes and wide flood plains associated with the Upper Waldo and Middle Waldo radar streamgages; whereas, the dynamic celerity was a better surrogate for shallow slopes and incised channels such as the Lower Waldo radar streamgage. The method demonstrates the potential extensibility of a post-wildfire warning system by (1) leveraging multiple systems (i.e., weather radar, near-field velocity and stage radars, and rain gages) for accurate and timely warnings of debris and floodflows, (2) establishing an order of operations to site, install, and operate near-field radars and conventional rain gages to record floodflows, forecast travel times, and document geomorphic change in this basin and hydrologically similar basins that lack data, and (3) communicating data operationally with the Colorado Department of Transportation engineering staff, National Weather Service forecasters, and emergency managers.

Colorado↗

Estimates of Southern White-tailed Ptarmigan daily nest survival from multiple sites in the Southern Rocky Mountains of Colorado

Estimating vital rates of avian species is important to understand population dynamics and develop potential conservation strategies that target rates for management. Avian species have reduced potential for high annual fecundity in alpine ecosystems due to a short breeding window and harsh weather conditions. We located nests from Southern White-tailed Ptarmigan ( Lagopus leucura altipetens ) across six study sites in the Southern Rocky Mountains of Colorado to estimate daily nest survival from 2013–2017. We used a known-fate hierarchical nest survival model and fit several covariates, including environmental conditions representing daily weather events and shrub cover, to describe variation in daily survival and derive estimates of nest success. We located and monitored 198 nests from 128 radio-marked ptarmigan hens. The mean nest success estimated as a derived parameter from daily nest survival was 45.6% (95% credible interval [CI]: 31.2–59.6%) and ranged from 40.3% to 50.3% across sites. Variation in daily nest survival was poorly described by the covariates we fit (95% CI of most slope coefficients overlapped 0), although there was some support for a negative effect of relative elevation (nests at lower elevations within a site survived at higher rates) and a positive effect of nest age (older nests survived at higher rates). We examined how variation in nest success was likely to influence the finite rate of population growth using a simple simulation with an age-transition matrix parameterized with previously reported fecundity and survival estimates. We found that the finite growth rate was predicted to increase 18.7% when evaluated from the lower to upper 95% CI estimated values of nest success, conditional on the other vital rates used in our simulation. We discuss the broader implications of these findings in the context of managing for nest survival of Southern White-tailed Ptarmigan.

Colorado↗

Geophysical identification of potential groundwater – surface water interactions in the Bonita Peak Mining District, San Juan County, Colorado

This report includes the results and interpretation of geophysical investigations to aid in the location, identification, and conceptual site model (CSM) development of the near surface groundwater and surface water hydrology of the Upper Cement Creek and California Gulch of the Animas Rivers within the Bonita Peak Mining Districts (BPMD). Mine influenced water (MIW) from numerous mine adits within the BPMD has altered the geochemistry and ecology of these ground and surface waters. Many studies of varying scales have investigated the geology, hydrogeology, geochemistry, and aqueous chemistry of these adits and impacted areas. Field investigations using non-invasive geophysical methods were conducted in the summer of 2019 to aid in identifying potential groundwater inputs to these streams and further the CSM of the near surface within and along these streams. Two field deployments were completed in August and September of 2019. In August 2019, Phase 1 was completed using multifrequency electromagnetic induction (EMI) and thermal infrared imaging (TIR) to measure the bulk electrical conductivity, magnetic susceptibility, and the temperature response at known and suspected groundwater seeps into the Upper Cement Creek and California Gulch stretches of these alpine streams. Results and interpretation from Phase 1 were used to co-locate stretches of these streams for further detailed fiber optic distributed temperature system (FODTS) monitoring of the temperature dynamics for submeter identification of potential groundwater seeps during Phase 2. Phase 2 included additional EMI measurements, FODTS measurements, and repeated some Phase 1 measurements for improved data quality. All these data were used for the combined qualitative interpretation identifying seep locations. The data from both Phase 1 and Phase 2 are available in the USGS data release at: https://www.sciencebase.gov/catalog/item/5f18a30082cef313ed84879b. These data include detailed near-surface electrical conductivity and magnetic susceptibility mapping, discrete TIR imaging, and several kilometers of FODTS data. These data were georeferenced and used to identify potential groundwater seeps into these streams. The results indicate a complex near surface geophysical and hydrogeologic regime. Seeps are classified and spatially correlated to the EMI, TIR, and geochemical sampling results. This report is a thorough compilation of the geophysical field efforts in the summer of 2019, the data processing, and interpretations of these data. These data alone should not be used for final interpretations of the complex hydrogeology within these alpine systems. Rather, these data should be combined with other known data from various investigations to further refine the CSM and guide decisions for future investigations, watershed protective measures, and other site decisions. Geographical information system (GIS) layers of these results are available to project managers for such integration and combined interpretations.

Colorado↗

Geophysical characterization of mine influenced groundwater and surface water in the Mayflower section of the Animas River, Bonita Peak Mining District, Silverton Colorado

This report details findings from geophysical investigations to identify possible groundwater - surface water interactions near the Mayflower Section of the Animas River in Silverton, Colorado. The Mayflower Section is a mine influenced Superfund Site in Colorado. This investigation utilized electromagnetic induction (EM or EMI), magnetic, and fiber optic distributed temperature system (FODTS) geophysical methods to measure the bulk earth electrical conductivity, magnetic susceptibility, and temperature of specific surveyed volumes of the earth near the Mayflower Section of the Animas River. These physical parameters are used to understand the groundwater – surface water interactions, which can guide decision makers in their assessment of mine-impacted surface water. This report details the results from characterization and monitoring technologies to provide high data density and continuous monitoring of bulk earth electrical conductivity and temperature in the Mayflower section of the Animas River to identify zones of groundwater – surface water interactions and potential metal loading from mine influenced water. The investigation separated right and left bank characterization for each method and indicates more groundwater is entering from the right bank than the left bank and these predominantly right bank discharges potentially contain metal-rich water compared to the left bank. Results also indicate mineral veins facilitate preferential groundwater discharge to the river due to possible jointing, fractures, and permeability differences sometimes occurring along veins relative to host rock. For example, Boulder Gulch is likely groundwater dominated and may be receive water impacted by Mayflower Impoundments #1 and #2. Additionally, the beaver ponds near Blair Gulch may influence groundwater discharge and Mayflower Impoundment #4 is possibly impacting groundwater and surface water and may be connected to the wetlands to the west of the impoundment. These data could be further analyzed for smaller spatial scale analysis within areas of interest. The identification of these locations along sections of the river likely impacted by mine influenced groundwater potentially entering the Animas River and can be used by site investigators, decision makers, and stakeholders in mitigation decisions and strategies.

Colorado↗

Precursors to a continental-arc ignimbrite flare-up: Early central volcanoes of the San Juan Mountains, Colorado, USA

Our newly acquired and recently published map, geochronologic, and compositional data for early intermediate-composition central volcanoes in the northeastern San Juan Mountains provide insights about the broad magmatic precursors to the large continental-arc ignimbrite flare-up in the mid-Cenozoic Southern Rocky Mountain volcanic field (SRMVF). Initial volcanism migrated from central Colorado to northern New Mexico ca. 40–29 Ma, as part of a more regional trend of southward-progressing mid-Cenozoic magmatism in the U.S. segment of the North American Cordillera. Within the San Juan locus, which represents the largest preserved erosional remnant of the SRMVF and site of most intense eruptive activity, new 40 Ar/ 39 Ar and U-Pb zircon ages show that eruptions at many individual edifices began nearly concurrently, at ca. 35 Ma, with peak activity at 34–32 Ma. Broadly similar precursor effusive volcanism characterizes other major loci of continental-arc ignimbrite magmatism along the western American cordilleras, but none of these sites records early volcanism as voluminous, spatially widespread, well exposed, or compositionally diverse as the San Juan locus in Colorado.

Colorado↗

The 2022 Chaos Canyon landslide in Colorado: Insights revealed by seismic analysis, field investigations, and remote sensing

An unusual, high-alpine, rapid debris slide originating in ice-rich debris occurred on June 28, 2022, at 16:33:16 MDT at the head of Chaos Canyon, a formerly glacier-covered valley in Rocky Mountain National Park, CO, USA. In this study, we integrate eyewitness videos and seismic records of the event with meteorological data, field observations, pre- and post-event satellite imagery, and uncrewed aircraft vehicle imagery to characterize the event and future hazards it may pose. Deformation of the eventual slide mass preceded rapid failure by decades, starting in the early to mid-2000s, accelerating in 2018 (the warmest year on record), and reaching ~ 20 m/year in 2021. The main event, which was preceded by smaller sliding episodes earlier that day, had a volume of ~ 2.1 million m 3 , reached peak velocities of about 5 m/s, slid on a surface up to 80 m deep, and moved up to ~ 245 m downslope in < 2 min. We observed blocks of frozen debris (permafrost) in the landslide deposits. Within ~ 2 weeks, these blocks had melted and became dry, conical debris mounds (molards). We hypothesize that the rapid slide was induced by gradually increasing long-term air temperatures that thawed ice and increased pore pressures. The presence and suspected influence of permafrost on the occurrence of this landslide indicate other slopes in the park, and other moderate-to-low latitude alpine regions may experience similar slope stability issues as temperatures continue to warm.

Colorado↗

Effects of hunting on mating, relatedness, and genetic diversity in a puma population

Hunting mortality can affect population abundance, demography, patterns of dispersal and philopatry, breeding, and genetic diversity. We investigated the effects of hunting on the reproduction and genetic diversity in a puma population in western Colorado, USA. We genotyped over 11,000 single nucleotide polymorphisms (SNPs), using double-digest, restriction site-associated DNA sequencing (ddRADseq) in 291 tissue samples collected as part of a study on the effects of hunting on puma population abundance and demography in Colorado from 2004 to 2014. The study was designed with a reference period (years 1–5), during which hunting was suspended, followed by a treatment period (years 6–10), in which hunting was reinstated. Our objectives were to examine the effects of hunting on: (1) paternity and male reproductive success; (2) the relatedness between pumas within the population, and (3) genetic diversity. We found that hunting reduced the average age of male breeders. The number of unique fathers siring litters increased each year without hunting and decreased each year during the hunting period. Mated pairs were generally unrelated during both time periods, and females were more closely related than males. Hunting was also associated with increased relatedness among males and decreased relatedness among females in the population. Finally, genetic diversity increased during the period without hunting and decreased each year when hunting was present. This study demonstrates the utility of merging demographic data with large-scale genomic datasets in order to better understand the consequences of management actions. Specifically, we believe that this study highlights the need for long-term experimental research in which hunting mortality is manipulated, including at least one non-harvested control population, as part of a broader adaptive, zone management scheme.

Colorado↗

Proximity to roads does not modify inorganic nitrogen deposition in a topographically complex, high traffic, subalpine forest

Vehicles are an important source for N deposition that may negatively impact roadside ecosystems. While elevated roadside N deposition has been found in many locations, it is not yet known if vehicle emissions cause measurable increases of N deposition in complex, mountainous terrain adjacent to roads. To address this, this study investigated the effect of vehicle N emissions on throughfall (through trees) and open N deposition in a high traffic corridor in mountainous terrain of Rocky Mountain National Park, Colorado, USA. We measured bulk (wet + dry) atmospheric N deposition in throughfall and open samplers along two transects of 750 (throughfall) and 225 (open) m moving away from the road using ion exchange resin (IER) collectors. Contrary to most studies of roadside N deposition, we found no influence of road proximity on inorganic N deposition in throughfall or open sites, possibly due to terrain complexity. Interactions with vegetation modified regional N deposition; throughfall sites had 69% higher nitrate (NO 3 − ) deposition than open sites and larger trees were associated with higher ammonium (NH 4 + ) deposition as compared to smaller trees. When comparing to regional sites that are part of national monitoring networks, we confirmed that our estimates were unaffected by vehicle emissions as our throughfall IER collectors had similar total inorganic N deposition as wet + dry deposition from regional sites (8.64–13.56 vs 10.72–12.14 g N ha −1 day −1 , respectively). These findings do not negate vehicles as a local source of N emissions but suggest elevated N deposition adjacent to busy roads cannot be assumed for complex terrains. Instead, environmental variables may be more important drivers than proximity to roads in topographically complex ecosystems.

Colorado↗

Conservation genomics of an endangered montane amphibian reveals low population structure, low genomic diversity and selection pressure from disease

Wildlife diseases are a major global threat to biodiversity. Boreal toads ( Anaxyrus [Bufo] boreas ) are a state-endangered species in the southern Rocky Mountains of Colorado and New Mexico, and a species of concern in Wyoming, largely due to lethal skin infections caused by the amphibian chytrid fungus Batrachochytrium dendrobatidis ( Bd ). We performed conservation and landscape genomic analyses using single nucleotide polymorphisms from double-digest, restriction site-associated DNA sequencing in combination with the development of the first boreal toad (and first North American toad) reference genome to investigate population structure, genomic diversity, landscape connectivity and adaptive divergence. Genomic diversity ( π = 0.00034–0.00040) and effective population sizes ( N e = 8.9–38.4) were low, likely due to post-Pleistocene founder effects and Bd -related population crashes over the last three decades. Population structure was also low, likely due to formerly high connectivity among a higher density of geographically proximate populations. Boreal toad gene flow was facilitated by low precipitation, cold minimum temperatures, less tree canopy, low heat load and less urbanization. We found >8X more putatively adaptive loci related to Bd intensity than to all other environmental factors combined, and evidence for genes under selection related to immune response, heart development and regulation and skin function. These data suggest boreal toads in habitats with Bd have experienced stronger selection pressure from disease than from other, broad-scale environmental variations. These findings can be used by managers to conserve and recover the species through actions including reintroduction and supplementation of populations that have declined due to Bd .

Colorado↗

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↗

A general approach for evaluating of the coverage, resolution, and representation of streamflow monitoring networks

Streamflow monitoring networks provide information for a wide range of public interests in river and streams. A general approach to evaluate monitoring for different interests is developed to support network planning and design. The approach defines three theoretically distinct information metrics (coverage, resolution, and representation) based on the spatial distribution of a variable of interest. Coverage is the fraction of information that a network can provide about a variable when some areas are not monitored. Resolution is the information available from the network relative to the maximum information possible given the number of sites in the network. Representation is the information that a network provides about a benchmark distribution of a variable. Information is defined using Shannon entropy where the spatial discretization of a variable among spatial elements of a landscape or sites in a network indicates the uncertainty in the spatial distribution of the variable. This approach supports the design of networks for monitoring of variables with heterogeneous spatial distributions (“hot spots” and patches) that might otherwise be unmonitored because they occupy insignificant portions of the landscape. Areas where monitoring will maintain or improve the metrics serve as objective priorities for public interests in network design. The approach is demonstrated for the streamflow monitoring network operated by the United States Geological Survey during water year 2020 indicating gaps in the coverage of coastal rivers and the resolution of low flows.

Colorado↗

Annotated bibliography of scientific research on Gunnison sage-grouse published from January 2005 to September 2022

Integrating recent scientific knowledge into management decisions supports effective natural resource management and can lead to better resource outcomes. However, finding and accessing scientific knowledge can be time consuming and costly. To assist in this process, the U.S. Geological Survey is creating a series of annotated bibliographies on topics of management concern for western lands. Previously published reports introduced a methodology for preparing annotated bibliographies to facilitate the integration of recent, peer-reviewed science into resource management decisions. Therefore, relevant text from those efforts is reproduced here to frame the presentation. Centrocercus minimus (Gunnison sage-grouse; hereafter GUSG) has been a focus of scientific investigation since the early 2000s. The U.S. Fish and Wildlife Service listed GUSG as threatened under the Endangered Species Act in 2014 because of declining populations and increasing habitat loss. The U.S. Fish and Wildlife Service, Bureau of Land Management, and Colorado Parks and Wildlife have sought to increase the conservation of this species by adapting management and recovery plans to reduce threats and increase population resiliency. GUSG are studied less than the closely related Centrocercus urophasianus (greater sage-grouse); however, research efforts have recently increased to understand the life history, genetics, and habitat suitability of this sagebrush-obligate species. We compiled and summarized peer-reviewed journal articles, data products, and formal technical reports (such as U.S. Department of Agriculture Forest Service General Technical Reports and U.S. Geological Survey Open-File Reports) on GUSG, published between January 2005 and September 2022. We first systematically searched three reference databases and three government databases using the following search phrases: “Gunnison sage-grouse” or “lesser sage-grouse” or “Gunnison grouse” or “Gunnison sage grouse” or “lesser sage grouse” or “ Centrocercus minimus .” We refined the initial list of products by removing (1) duplicates, (2) publications that were not published as research, data products, or scientific review articles in peer-reviewed journals or as formal technical reports, and (3) products that did not have GUSG as a research focus or products that did not present new data or findings about GUSG. We summarized each product using a consistent structure (background, objectives, methods, location, findings, and implications) and identified the management topics (for example, population estimates or targets, habitat, and management efforts) addressed by each product. We also noted which publications included new geospatial data. The review process for this annotated bibliography included two initial internal colleague reviews of each summary, requesting input on each summary from an author of the original publication, and a formal peer review. Our initial searches resulted in 80 total products, of which 63 met our criteria for inclusion of which 53 were products that had not been summarized before. Across products summarized in the annotated bibliography, broad-scale habitat characteristics; population estimates or targets; behavior or demographics; and genetics were the most commonly addressed management topics. The bibliographies are available on the Science for Resource Managers tool ( https://apps.usgs.gov/science-for-resource-managers ) and are searchable by topic, location, and year, and the search tool includes links to each original publication. The studies compiled and summarized in this annotated bibliography may inform planning and management actions that seek to maintain and restore sagebrush landscapes and GUSG populations across the GUSG range.

Colorado↗

Local topography and streambed hydraulic conductivity influence riparian groundwater age and groundwater-surface water connection

The western U.S. is experiencing increasing rain to snow ratios due to climate change, and scientists are uncertain how changing recharge patterns will affect future groundwater-surface water connection. We examined how watershed topography and streambed hydraulic conductivity impact groundwater age and stream discharge at eight sites along a headwater stream within the Manitou Experimental Forest, CO USA. To do so, we measured: (a) continuous stream and groundwater discharge/level and specific conductivity from April to November 2021; (b) biweekly stream and groundwater chemistry; (c) groundwater chlorofluorocarbons and tritium in spring and fall; (d) streambed hydraulic conductivity; and (e) local slope. We used the chemistry data to calculate fluorite saturation states that were used to inform end-member mixing analysis of streamflow source. We then combined chlorofluorocarbon and tritium data to estimate the age composition of riparian groundwater. Our data suggest that future stream drying is more probable where local slope is steep and streambed hydraulic conductivity is high. In these areas, groundwater source shifted seasonally, as indicated by age increases, and we observed a high fraction of groundwater in streamflow, primarily interflow from adjacent hillslopes. In contrast, where local slope is flat and streambed hydraulic conductivity is low, streamflow is more likely to persist as groundwater age was seasonally constant and buffered by storage in alluvial sediments. Groundwater age and streamflow paired with characterization of watershed topography and subsurface characteristics enabled identification of likely controls on future stream drying patterns.

Colorado↗