Geology topics
Research about Colorado
Source-linked reports with geographic coverage including Colorado.
Macroinvertebrate community sample collection methods and data collected from Sand Creek and Medano Creek, Great Sand Dunes National Park and Preserve, Colorado, 2005–07
This report provides a table of site descriptions, sample information, and semiquantitative aquatic macroinvertebrate data from 105 samples collected between 2005 and 2007 from 7 stream sites within the Sand Creek and Medano Creek watersheds in Great Sand Dunes National Park and Preserve, Saguache County, Colorado. Additionally, a short description of sample collection methods and laboratory sample processing procedures is presented. These data were collected in anticipation of assessing the potential effects of fish toxicants on macroinvertebrates.
2014 annual summary of the lower Gunnison River Basin Selenium Management Program water-quality monitoring, Colorado
Dissolved-selenium loading analyses of data collected at 18 water-quality sites in the lower Gunnison River Basin in Colorado were completed through water year (WY) 2014. A WY is defined as October 1–September 30. Selenium is a trace element that bioaccumulates in aquatic food chains and can cause reproductive failure, deformities, and other harmful effects. This report presents information on the dissolved-selenium loads at 18 sites in the lower Gunnison River Basin for WYs 2011–2014. Annual dissolved-selenium loads were calculated at 5 sites with continuous U.S. Geological Survey (USGS) streamflow gages, whereas instantaneous dissolved-selenium loads were calculated for the remaining 13 sites using water-quality samples that had been collected periodically during WYs 2011–2014. Annual dissolved-selenium loads for WY 2014 ranged from 336 pounds (lb) at Uncompahgre River at Colona to 13,300 lb at Gunnison River near Grand Junction (Whitewater). Most sites in the basin had a median instantaneous dissolved-selenium load of less than 20.0 lb per day. In general, dissolved-selenium loads at Gunnison River main-stem sites showed an increase from upstream to downstream. The State of Colorado water-quality standard for dissolved selenium of 4.6 micrograms per liter (µg/L) was compared to the 85th percentiles for dissolved selenium at selected water-quality sites. Annual 85th percentiles for dissolved selenium were calculated for the five core USGS sites having streamflow gages using estimated dissolved-selenium concentrations from linear regression models. These annual 85th percentiles in WY 2014 ranged from 0.97 µg/L at Uncompahgre River at Colona to 16.7 µg/L at Uncompahgre River at Delta. Uncompahgre River at Delta and Whitewater were the only core sites where water samples exceeded the State of Colorado water-quality standard for dissolved selenium of 4.6 µg/L. Instantaneous 85th percentiles for dissolved selenium were calculated for sites with sufficient data using water-quality samples collected during WYs 2011–2014. The instantaneous 85th percentiles for samples for WY 2014 ranged from 1.1 µg/L at Uncompahgre River at Colona to 125 µg/L at Loutzenhizer Arroyo at North River Road. A trend analysis was completed for Whitewater to determine if dissolved-selenium loads are increasing or decreasing. The trend analysis indicates a decrease of 8,000 lb from WY 1986 to WY 2014, a 34.8 percent reduction during the time period, and an additional 6.2 percent reduction from a reported 28.6 percent reduction during WYs 1986–2008. The trend analysis for WY 1992 to WY 2014 indicates a decrease of 5,800 lb per year, or 27.9 percent.
Evaluation of National Atmospheric Deposition Program measurements for colocated sites CO89 and CO98 at Rocky Mountain National Park, water years 2010–14
Atmospheric wet-deposition monitoring in Rocky Mountain National Park included precipitation depth and aqueous chemical measurements at colocated National Atmospheric Deposition Program/National Trends Network (NADP/NTN) sites CO89 and CO98 (Loch Vale) during water years 2010–14 (study period). The colocated sites were separated by approximately 6.5 meters horizontally and 0.5 meter in elevation, in accordance with NADP siting criteria. Assessment of the 5-year record of colocated data is intended to inform man-agement decisions pertaining to the achievement of nitrogen deposition reduction goals of the Rocky Mountain National Park Nitrogen Deposition Reduction Plan. The data at site CO98 met NADP completeness criteria for the first time in 29 years of operation in 2011 and then again in 2012. During the study period, data at site CO89 met completeness criteria in 2012. Median weekly relative precipitation-depth differences between sites CO89 and CO98 ranged from 0 to 0.25 millimeter during the study period. Median weekly absolute percent differences in sample volume ranged from 5 to 10 percent. Median relative concentration differences for weekly ammonium (NH 4 + ) and nitrate (NO 3 - ) concentrations were near the NADP Central Analytical Laboratory’s method detection limits and thus were considered small. Absolute percent differences for water-year 2010–14 precipitation-weighted mean concentrations of NH 4 + , NO 3 - , and inorganic nitrogen (N inorg ) ranged from 0.0 to 25.7 percent. Absolute percent differences for water-year 2010–14 NH 4 + , NO 3 - , and N inorg deposition ranged from 2.1 to 18.9 percent, 3.3 to 24.5 percent, and 0.3 to 17.4 percent, respectively.
Effects of seasonal weather on breeding phenology and reproductive success of alpine ptarmigan in Colorado
Animal populations occurring at high elevations are often assumed to be in peril of extinctions or local extirpations due to elevational-dispersal limitations and thermoregulatory constraints as habitats change and warm. However, long-term monitoring of high-elevation populations is uncommon relative to those occurring at lower elevations, and evidence supporting this assumption is limited. We analyzed 45 years of reproductive data for two Colorado populations of white-tailed ptarmigan ( Lagopus leucura ), an alpine-endemic species with restricted distribution in western North America. Seasonal temperatures measured by the number of growing degree days warmed significantly at our study sites for pre-nesting, nesting, and brood-rearing seasonal periods (mean advance of 8 growing degree days per decade), and both populations advanced their reproductive phenology over the study period based on median hatch dates (median advance of 3.7 and 1.9 days per decade for the northern and southern sites, respectively). Reproductive performance measured by the number of chicks per hen declined significantly at one study site but not the other, and differences between sites may have been due to habitat degradation at one study area. Annual variability in chicks per hen was large at both sites but only weakly related to seasonal weather. An index of precipitation and temperature during the brood-rearing period was the best predictor for reproductive success with warm and dry conditions relating positively to number of chicks per hen. Our results provide evidence for two alpine ptarmigan populations that are remarkably invariant to fluctuations in seasonal weather with respect to reproductive success as measured by number of chicks per hen in the breeding population. These results are surprising given the general perception of alpine animal populations as being highly sensitive to warming temperatures.
Rapid response, monitoring, and mitigation of induced seismicity near Greeley, Colorado
On 1 June 2014 (03:35 UTC), an M w 3.2 earthquake occurred in Weld County, Colorado, a historically aseismic area of the Denver–Julesburg basin. Weld County is a prominent area of oil and gas development, including many high‐rate class II wastewater injection wells. In the days following the earthquake, the University of Colorado, with support from the U.S. Geological Survey and Incorporated Research Institutions for Seismology–Portable Array Seismic Studies of the Continental Lithosphere, rapidly deployed six seismic stations to characterize the seismicity associated with the 1 June earthquake (the Greeley sequence) and to investigate its possible connection to wastewater disposal. The spatial and temporal proximity of earthquakes to a high‐rate wastewater disposal well strongly suggests these earthquakes were induced. Scientific communication between the university, state agencies, and the energy industry led to rapid mitigation strategies to reduce the occurrence of further earthquakes. Mitigation efforts included implementing a temporary moratorium on injection at the well, cementing the bottom portion of the disposal well to minimize hydrologic connectivity between the disposal formation and the underlying crystalline basement, and subsequently allowing injection to resume at lower rates. Following the resumption of wastewater disposal, microseismicity was closely monitored for both increases in earthquake rate and magnitude. Following mitigation efforts, between 13 August 2014 and 29 December 2015, no earthquakes larger than M 1.5 occurred near the Greeley sequence. This study demonstrates that a detailed and rapid characterization of a seismic sequence in space and time relative to disposal, combined with collaboration and communication between scientists, regulators, and industry, can lead to objective and actionable mitigation efforts that potentially reduced the rate of earthquakes and the possible generation of larger earthquakes.
Climate regulates alpine lake ice cover phenology and aquatic ecosystem structure
High-elevation aquatic ecosystems are highly vulnerable to climate change, yet relatively few records are available to characterize shifts in ecosystem structure or their underlying mechanisms. Using a long-term dataset on seven alpine lakes (3126 to 3620 m) in Colorado, USA, we show that ice-off dates have shifted seven days earlier over the past 33 years and that spring weather conditions – especially snowfall – drive yearly variation in ice-off timing. In the most well-studied lake, earlier ice-off associated with increases in water residence times, thermal stratification, ion concentrations, dissolved nitrogen, pH, and chlorophyll-a. Mechanistically, low spring snowfall and warm temperatures reduce summer stream flow (increasing lake residence times) but enhance melting of glacial and permafrost ice (increasing lake solute inputs). The observed links among hydrological, chemical, and biological responses to climate factors highlight the potential for major shifts in the functioning of alpine lakes due to forecasted climate change.
Response of selenium concentrations in groundwater to seasonal canal leakage, lower Gunnison River Basin, Colorado, 2013
Selenium is a water-quality concern in the lower Gunnison River Basin because irrigation water interacting with seleniferous soils derived from the Mancos Shale Formation has mobilized selenium and increased its concentrations in surface water. Understanding the occurrence of elevated selenium concentrations in groundwater is necessary because groundwater discharge is an important source of selenium in surface water in the basin. In 2013, the U.S. Geological Survey, in cooperation with the Bureau of Reclamation and the Colorado Water Conservation Board, began a study to understand how changes in groundwater levels attributed to canal leakage affected the concentrations and speciation of dissolved selenium in groundwater. The purpose of this report is to characterize the groundwater adjacent to an unlined leaky canal. Two locations, near the East Canal (W-N1 and W-N2) and farther from the East Canal (W-M1 and W-M2), were selected for nested monitoring well installations. The pressure exerted by changes in canal stage was more readily transferred to the deep groundwater measured in the W-N1 near the canal than the shallow groundwater at the W-N2 well. No definitive relation could be made between canal water-level elevation and water-level elevations in monitoring wells farther from the canal (W-M1 and W-M2). Water flowing through the East Canal before the irrigation season had much higher selenium concentrations (140 micrograms per liter) than water in the canal during the irrigation season (3.02 micrograms per liter). Total selenium concentrations in the monitoring wells near the canal initially increased to 51.8 micrograms per liter in W-N1 and 1.66 micrograms per liter in W-N2. The initial increase in groundwater selenium concentrations presumably resulted from the dissolution of salts in the unsaturated zone by rising groundwater levels associated with canal leakage. The subsequent decrease in total selenium concentrations resulted from a combination of dilution by canal leakage and selenium reduction processes. Total selenium concentrations in monitoring wells located farther from the canal were not directly affected by canal leakage. Selenite/total selenium mass ratios in the East Canal samples ranged from about 0.02 to 0.13, indicating that about 2 to 13 percent of the total selenium in canal samples was composed of selenite. The increase in total selenium at W-N1 from before the irrigation season to the early irrigation season was accompanied by a decrease in the percentage of selenite from about 10 to 1 percent, indicating that selenate was added to the groundwater. A nitrate pulse occurred with the selenate pulse in W-N1 at the beginning of the irrigation season but apparently dissipated to a low enough concentration during the irrigation season to allow for selenate reduction to occur, as indicated by the relatively high percentages of selenite in W-N1 during the late irrigation season. W-N2 generally contained higher percentages of selenite than W-N1. Percentages of selenite in W-M1 did not change in response to filling the canal and generally composed less than 1 percent of the total selenium in that well. The predominance of selenate in W-M1, and apparent lack of selenate reduction, cannot be explained by a lack of anoxic conditions in the groundwater because all the available dissolved-oxygen data indicate that concentrations were less than 0.5 milligrams per liter. The most likely explanation for the lack of selenate reduction in W-M1 is that the exceptionally high concentrations of nitrate in the groundwater (about 340 to 390 milligrams per liter as nitrogen) inhibited selenate reduction. These high nitrate concentrations presumably come from the Mancos Shale and its weathering products because there was no evidence for a human source of nitrate at the lower Gunnison River Basin wetland. The high concentrations of selenate in W-M1 may persist and eventually discharge to surface water unless nitrate concentrations are reduced to low enough levels to permit substantial selenate reduction to occur. Well W-M2 contained relatively low concentrations of total selenium and high percentages of selenite before and at the onset of the irrigation season. An increase in total selenium concentration associated with a drying and wetting period later in the summer was accompanied by a decrease in the percentage of selenite to near 0 percent, indicating that selenate was added to the groundwater. This pattern is consistent with the examples of increasing concentrations of total selenium in the other wells and presumably resulted from the dissolution of selenate-bearing salts in the unsaturated zone by rising water levels in W-M2.
Three-dimensional surface deformation derived from airborne interferometric UAVSAR: Application to the Slumgullion Landslide
In order to provide surface geodetic measurements with “landslide-wide” spatial coverage, we develop and validate a method for the characterization of 3-D surface deformation using the unique capabilities of the Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) airborne repeat-pass radar interferometry system. We apply our method at the well-studied Slumgullion Landslide, which is 3.9 km long and moves persistently at rates up to ∼2 cm/day. A comparison with concurrent GPS measurements validates this method and shows that it provides reliable and accurate 3-D surface deformation measurements. The UAVSAR-derived vector velocity field measurements accurately capture the sharp boundaries defining previously identified kinematic units and geomorphic domains within the landslide. We acquired data across the landslide during spring and summer and identify that the landslide moves more slowly during summer except at its head, presumably in response to spatiotemporal variations in snowmelt infiltration. In order to constrain the mechanics controlling landslide motion from surface velocity measurements, we present an inversion framework for the extraction of slide thickness and basal geometry from dense 3-D surface velocity fields. We find that the average depth of the Slumgullion Landslide is 7.5 m, several meters less than previous depth estimates. We show that by considering a viscoplastic rheology, we can derive tighter theoretical bounds on the rheological parameter relating mean horizontal flow rate to surface velocity. Using inclinometer data for slow-moving, clay-rich landslides across the globe, we find a consistent value for the rheological parameter of 0.85 ± 0.08.
Episodic bedrock erosion by gully-head migration, Colorado High Plains, USA
This study explores the frequency of bedrock exposure in a soil-mantled low-relief (i.e. non-mountainous) landscape. In the High Plains of eastern Colorado, gully headcuts are among the few erosional features that will incise through the soil mantle to expose bedrock. We measured the last time of bedrock exposure using optically stimulated luminescence dating of alluvial sediment overlying bedrock in gully headcuts. Our dating suggests that headcuts in adjacent gullies expose bedrock asynchronously, and therefore, the headcuts are unlikely to have been triggered by a base-level drop in the trunk stream. This finding supports the hypothesis that headcuts can develop locally in gullies as a result of focused scour in locations where hydraulic stress during a flash flood is sufficiently high, and/or ground cover is sufficiently weak, to generate a scour hole that undermines vegetation. Alluvium dating also reveals that gullies have been a persistent part of this landscape since the early Holocene.
Debris flow occurrence and sediment persistence, Upper Colorado River Valley, CO
Debris flow magnitudes and frequencies are compared across the Upper Colorado River valley to assess influences on debris flow occurrence and to evaluate valley geometry effects on sediment persistence. Dendrochronology, field mapping, and aerial photographic analysis are used to evaluate whether a 19th century earthen, water-conveyance ditch has altered the regime of debris flow occurrence in the Colorado River headwaters. Identifying any shifts in disturbance processes or changes in magnitudes and frequencies of occurrence is fundamental to establishing the historical range of variability (HRV) at the site. We found no substantial difference in frequency of debris flows cataloged at eleven sites of deposition between the east (8) and west (11) sides of the Colorado River valley over the last century, but four of the five largest debris flows originated on the west side of the valley in association with the earthen ditch, while the fifth is on a steep hillslope of hydrothermally altered rock on the east side. These results suggest that the ditch has altered the regime of debris flow activity in the Colorado River headwaters as compared to HRV by increasing the frequency of debris flows large enough to reach the Colorado River valley. Valley confinement is a dominant control on response to debris flows, influencing volumes of aggradation and persistence of debris flow deposits. Large, frequent debris flows, exceeding HRV, create persistent effects due to valley geometry and geomorphic setting conducive to sediment storage that are easily delineated by valley confinement ratios which are useful to land managers.
Estimates of peak flood discharge for 21 sites in the Front Range in Colorado in response to extreme rainfall in September 2013
Extreme rainfall in September 2013 caused destructive floods in part of the Front Range in Boulder County, Colorado. Erosion from these floods cut roads and isolated mountain communities for several weeks, and large volumes of eroded sediment were deposited downstream, which caused further damage of property and infrastructures. Estimates of peak discharge for these floods and the associated rainfall characteristics will aid land and emergency managers in the future. Several methods (an ensemble) were used to estimate peak discharge at 21 measurement sites, and the ensemble average and standard deviation provided a final estimate of peak discharge and its uncertainty. Because of the substantial erosion and deposition of sediment, an additional estimate of peak discharge was made based on the flow resistance caused by sediment transport effects. Although the synoptic-scale rainfall was extreme (annual exceedance probability greater than 1,000 years, about 450 millimeters in 7 days) for these mountains, the resulting peak discharges were not. Ensemble average peak discharges per unit drainage area (unit peak discharge, [ Q u ]) for the floods were 1–2 orders of magnitude less than those for the maximum worldwide floods with similar drainage areas and had a wide range of values (0.21–16.2 cubic meters per second per square kilometer [m 3 s -1 km -2 ]). One possible explanation for these differences was that the band of high-accumulation, high-intensity rainfall was narrow (about 50 kilometers wide), oriented nearly perpendicular to the predominant drainage pattern of the mountains, and therefore entire drainage areas were not subjected to the same range of extreme rainfall. A linear relation (coefficient of determination [ R 2 ]=0.69) between Q u and the rainfall intensity ( I Tc , computed for a time interval equal to the time-of-concentration for the drainage area upstream from each site), had the form: Q u =0.26( I Tc -8.6), where the coefficient 0.26 can be considered to be an area-averaged peak runoff coefficient for the September 2013 rain storms in Boulder County, and the 8.6 millimeters per hour to be the rainfall intensity corresponding to a soil moisture threshold that controls the soil infiltration rate. Peak discharge estimates based on the sediment transport effects were generally less than the ensemble average and indicated that sediment transport may be a mechanism that limits velocities in these types of mountain streams such that the Froude number fluctuates about 1 suggesting that this type of floodflow can be approximated as critical flow.
Illuminating wildfire erosion and deposition patterns with repeat terrestrial lidar
Erosion following a wildfire is much greater than background erosion in forests because of wildfire-induced changes to soil erodibility and water infiltration. While many previous studies have documented post-wildfire erosion with point and small plot-scale measurements, the spatial distribution of post-fire erosion patterns at the watershed scale remains largely unexplored. In this study lidar surveys were collected periodically in a small, first-order drainage basin over a period of 2 years following a wildfire. The study site was relatively steep with slopes ranging from 17° to > 30°. During the study period, several different types of rain storms occurred on the site including low-intensity frontal storms (2.4 mm h −1 ) and high-intensity convective thunderstorms (79 mm h −1 ). These storms were the dominant drivers of erosion. Erosion resulting from dry ravel and debris flows was notably absent at the site. Successive lidar surveys were subtracted from one another to obtain digital maps of topographic change between surveys. The results show an evolution in geomorphic response, such that the erosional response after rain storms was strongly influenced by the previous erosional events and pre-fire site morphology. Hillslope and channel roughness increased over time, and the watershed armored as coarse cobbles and boulders were exposed. The erosional response was spatially nonuniform; shallow erosion from hillslopes (87% of the study area) contributed 3 times more sediment volume than erosion from convergent areas (13% of the study area). However, the total normalized erosion depth (volume/area) was highest in convergent areas. From a detailed understanding of the spatial locations of erosion, we made inferences regarding the processes driving erosion. It appears that hillslope erosion is controlled by rain splash (for detachment) and overland flow (for transport and quasi-channelized erosion), with the sites of highest erosion corresponding to locations with the lowest roughness. By contrast, in convergent areas we found erosion caused by overland flow. Soil erosion was locally interrupted by immobile objects such as boulders, bedrock, or tree trunks, resulting in a patchy erosion network with increasing roughness over time.
Groundwater quality, age, and susceptibility and vulnerability to nitrate contamination with linkages to land use and groundwater flow, Upper Black Squirrel Creek Basin, Colorado, 2013
The Upper Black Squirrel Creek Basin is located about 25 kilometers east of Colorado Springs, Colorado. The primary aquifer is a productive section of unconsolidated deposits that overlies bedrock units of the Denver Basin and is a critical resource for local water needs, including irrigation, domestic, and commercial use. The primary aquifer also serves an important regional role by the export of water to nearby communities in the Colorado Springs area. Changes in land use and development over the last decade, which includes substantial growth of subdivisions in the Upper Black Squirrel Creek Basin, have led to uncertainty regarding the potential effects to water quality throughout the basin. In response, the U.S. Geological Survey, in cooperation with Cherokee Metropolitan District, El Paso County, Meridian Service Metropolitan District, Mountain View Electric Association, Upper Black Squirrel Creek Groundwater Management District, Woodmen Hills Metropolitan District, Colorado State Land Board, and Colorado Water Conservation Board, and the stakeholders represented in the Groundwater Quality Study Committee of El Paso County conducted an assessment of groundwater quality and groundwater age with an emphasis on characterizing nitrate in the groundwater. Groundwater-quality samples were collected from 50 randomly selected wells between May and June 2013. The samples were analyzed for major ions, nutrients, dissolved gases, tritium ( 3 H), chlorofluorocarbons (CFC-11, CFC-12, and CFC-113), and fuel products (such as benzene, toluene, ethylbenzene, and xylenes). None of the groundwater samples exceeded the U.S. Environmental Protection Agency (EPA) National Primary Drinking Water Regulations for primary maximum contaminant levels (MCL) for major ions. Secondary maximum contaminant levels, which are not health concerns and affect mainly taste, color, or odor of the water, were observed in rare instances for pH (2 samples), chloride (1 sample), iron (3 samples), and manganese (8 samples). The secondary maximum contaminant level for total dissolved solids was also exceeded for two samples. Nitrate (nitrite plus nitrate as nitrogen in groundwater) was elevated above the estimated background concentration of natural recharge waters of 1 milligram per liter (mg/L) in 44 of the 50 wells sampled and showed a median concentration of 5.4 mg/L. Nitrate concentrations were above the MCL of 10 mg/L in 5 of the 50 wells sampled and above half of the EPA MCL (5 mg/L) in 27 of the 50 wells sampled, which included samples above the MCL. Dissolved-oxygen concentrations exceeded 0.5 mg/L in 95 percent of reported values (40 of 42 samples) and exceeded 2.0 mg/L in 90 percent of reported values (38 of 42 samples). The oxidized conditions observed in most areas indicate that nitrate from fertilizers and animal or human waste was geochemically stable and could persist in the groundwater for decades or perhaps longer. A historical analysis of median nitrate concentrations over nearly three decades showed an increase in nitrate of approximately 1 mg/L from 4.3 to 5.4 mg/L, although the increase was not determined to be significantly different using nonparametric statistical methods. Major-ion data indicate that groundwater representative of the primary aquifer was classified as calcium-sodium bicarbonate type water. Other water samples from wells located mainly along the periphery of the primary aquifer had cation-anion compositions consistent with distinct water sources, including groundwater contributions from the underlying bedrock aquifers. The areas with differentiable water sources were located mainly where alluvial deposits were thin and geologic contacts to the underlying bedrock aquifers were relatively shallow. Nitrate concentrations in the groundwater were evaluated for relations to land use. An agricultural region was defined using a sequence of land satellite imagery. Groundwater flow directions interpreted from median water-table elevations measured from 2000 to 2013 were used in conjunction with cropland locations to define the agricultural region boundaries by encompassing potential pathways of nitrate transport in the groundwater from nitrogen-based fertilizers. A statistically significant higher median nitrate concentration was observed for areas inside the agricultural region (6.7 mg/L) compared to areas outside the agricultural region (2.3 mg/L), although median concentrations in both areas were below the MCL (10 mg/L). Median nitrate concentration was also significantly greater in land parcels with septic use (4.9 mg/L) compared to nonseptic parcels (1.7 mg/L). In general, agriculture or septic use was identified as the primary source of nitrate, depending on location, while commercial, county, grazing, and residential land uses were generally secondary sources of nitrate. Apparent groundwater ages were estimated from chlorofluorocarbons (CFC-11, CFC-12, and CFC-113) and tritium ( 3 H) data using models that assumed piston flow and binary mixing (dilution of a young component with old, tracer-free water). The mean and median groundwater ages were about 30 years and the standard deviation was 6 years, indicating that most groundwater in the primary aquifer was “young” water that had recharged to the aquifer over the last few decades (post-1950s). The median fraction of young water was about 71 percent, and the standard deviation was 29 percent. The remaining water predated the 1950s, which may have originated from deeper geologic formations or may represent slow moving groundwater within the primary aquifer. Some of the oldest groundwater ages (older than 30 years) were observed in the upper reaches of the aquifer to the northwest where the primary aquifer is thin and intersects bedrock, supporting the hypothesis of geochemically distinct groundwater entering the primary aquifer from below. Groundwater that had reached the central part of the aquifer from upgradient areas of the basin was variable in age because of differences in flow paths and travel velocities. The groundwater age analysis showed that current (2013) land-use practices could affect water quality over decades to come, and that responses to remedial actions could be slow, especially for constituents, such as nitrate, that are stable under oxidized conditions. Fuel products (including acetone, benzene, diisopropyl ether, ethylbenzene, methyl acetate, methyl tertiary butyl ether (MTBE), methyl tert-pentyl ether, m- + p-xylene, o-xylene, tert-amyl alcohol, tert-butyl alcohol, tert-butyl ethyl ether, and toluene) were analyzed in groundwater from 49 of the 50 wells. Water from seven sites had detections for fuel compounds; all concentrations were below MCL. The results provided assurance of water quality and a valuable baseline to evaluate future trends of fuel constituents as the region is further developed. Probability maps were developed from logistic regression models to examine the likelihood that nitrate concentrations in groundwater exceeded specified levels. Susceptibility analysis examined relations between mid-level (5.0 mg/L) nitrate concentrations and climatic, hydrologic, and geologic variables; the significant variables were identified as depth to groundwater, soil organic matter, and soil water storage to 25-centimeter (cm) depth. The vulnerability assessments included natural factors driving susceptibility but also human factors related to land use and septic use. Vulnerability to low-level (2.5 mg/L) nitrate was related to depth to groundwater, septic zoning, and soil organic matter. The results highlighted that septic zoning affected low-level nitrate concentrations. Vulnerability to mid-level (5.0 mg/L) nitrate was examined using all 50 samples and also with two data outliers removed, which showed relatively high nitrate concentrations but also anomalous water chemistry or were located beyond the primary study area. Vulnerability to mid-level (5.0 mg/L) nitrate using all 50 samples was related to depth to groundwater, land use, septic use within a 500-meter (m) radius, soil water storage to a 25-cm depth, soil organic matter, and whether a location was within the agricultural region. The mid-level (5.0 mg/L) vulnerability model using 48 samples (two outliers removed) produced the best overall fit and was related to the same variables as when using all samples except septic use. The results for mid-level vulnerability provided additional support that septic use was associated with low levels of nitrate in the groundwater. Soil properties and land use were identified as the main drivers of moderate nitrate concentrations. Probabilities of exceeding low-level nitrate concentrations were high in most areas with the lowest probabilities usually to the northwest along thin geologic deposits in the upper part of the basin. The results of this investigation offer the foundational information needed for developing best management practices to mitigate nitrate contamination, basic concepts on water quality to aid public education, and information to guide regulatory measures if policy makers determine this is warranted. Science-based decision making will require continued monitoring and analysis of water quality in the future.
Rock-avalanche dynamics revealed by large-scale field mapping and seismic signals at a highly mobile avalanche in the West Salt Creek valley, western Colorado
On 25 May 2014, a rain-on-snow–induced rock avalanche occurred in the West Salt Creek valley on the northern flank of Grand Mesa in western Colorado (United States). The avalanche mobilized from a preexisting rock slide in the Green River Formation and traveled 4.6 km down the confined valley, killing three people. The avalanche was rare for the contiguous United States because of its large size (54.5 Mm3) and high mobility (height/length = 0.14). To understand the avalanche failure sequence, mechanisms, and mobility, we conducted a forensic analysis using large-scale (1:1000) structural mapping and seismic data. We used high-resolution, unmanned aircraft system imagery as a base for field mapping, and analyzed seismic data from 22 broadband stations (distances < 656 km from the rock-slide source area) and one short-period network. We inverted broadband data to derive a time series of forces that the avalanche exerted on the earth and tracked these forces using curves in the avalanche path. Our results revealed that the rock avalanche was a cascade of landslide events, rather than a single massive failure. The sequence began with an early morning landslide/debris flow that started ∼10 h before the main avalanche. The main avalanche lasted ∼3.5 min and traveled at average velocities ranging from 15 to 36 m/s. For at least two hours after the avalanche ceased movement, a central, hummock-rich core continued to move slowly. Since 25 May 2014, numerous shallow landslides, rock slides, and rock falls have created new structures and modified avalanche topography. Mobility of the main avalanche and central core was likely enhanced by valley floor material that liquefied from undrained loading by the overriding avalanche. Although the base was likely at least partially liquefied, our mapping indicates that the overriding avalanche internally deformed predominantly by sliding along discrete shear surfaces in material that was nearly dry and had substantial frictional strength. These results indicate that the West Salt Creek avalanche, and probably other long-traveled avalanches, could be modeled as two layers: a thin, liquefied basal layer, and a thicker and stronger overriding layer.
Mountains, glaciers, and mines—The geological story of the Blue River valley, Colorado, and its surrounding mountains
This report describes, in a nontechnical style, the geologic history and mining activity in the Blue River region of Colorado, which includes all of Summit County. The geologic story begins with the formation of ancient basement rocks, as old as about 1700 million years, and continues with the deposition of sedimentary rocks on a vast erosional surface beginning in the Cambrian Period (about 530 million years ago). This deposition was interrupted by uplift of the Ancestral Rocky Mountains during the late Paleozoic Era (about 300 million years ago). The present Rocky Mountains began to rise at the close of the Mesozoic Era (about 65 million years ago). A few tens of millions years ago, rifting began to form the Blue River valley; a major fault along the east side of the Gore Range dropped the east side down, forming the present valley. The valley once was filled by sediments and volcanic rocks that are now largely eroded. During the last few hundred-thousand years, at least two periods of glaciation sculpted the mountains bordering the valley and glaciers extended down the Blue River valley as far south as present Dillon Reservoir. Discovery of deposits of gold, silver, copper, and zinc in the late 1800s, particularly in the Breckenridge region, brought an influx of early settlers. The world-class molybdenum deposit at Climax, mined since the First World War, reopened in 2012 after a period of closure. The report includes a glossary to explain geologic terms used in the text, and numerous photos, maps, and diagrams illustrate the geologic principles discussed. References for further reading are also included.
Manual hierarchical clustering of regional geochemical data using a Bayesian finite mixture model
Interpretation of regional scale, multivariate geochemical data is aided by a statistical technique called “clustering.” We investigate a particular clustering procedure by applying it to geochemical data collected in the State of Colorado, United States of America. The clustering procedure partitions the field samples for the entire survey area into two clusters. The field samples in each cluster are partitioned again to create two subclusters, and so on. This manual procedure generates a hierarchy of clusters, and the different levels of the hierarchy show geochemical and geological processes occurring at different spatial scales. Although there are many different clustering methods, we use Bayesian finite mixture modeling with two probability distributions, which yields two clusters. The model parameters are estimated with Hamiltonian Monte Carlo sampling of the posterior probability density function, which usually has multiple modes. Each mode has its own set of model parameters; each set is checked to ensure that it is consistent both with the data and with independent geologic knowledge. The set of model parameters that is most consistent with the independent geologic knowledge is selected for detailed interpretation and partitioning of the field samples.
High resolution mapping of development in the wildland-urban interface using object based image extraction
The wildland-urban interface (WUI), the area where human development encroaches on undeveloped land, is expanding throughout the western United States resulting in increased wildfire risk to homes and communities. Although census based mapping efforts have provided insights into the pattern of development and expansion of the WUI at regional and national scales, these approaches do not provide sufficient detail for fine-scale fire and emergency management planning, which requires maps of individual building locations. Although fine-scale maps of the WUI have been developed, they are often limited in their spatial extent, have unknown accuracies and biases, and are costly to update over time. In this paper we assess a semi-automated Object Based Image Analysis (OBIA) approach that utilizes 4-band multispectral National Aerial Image Program (NAIP) imagery for the detection of individual buildings within the WUI. We evaluate this approach by comparing the accuracy and overall quality of extracted buildings to a building footprint control dataset. In addition, we assessed the effects of buffer distance, topographic conditions, and building characteristics on the accuracy and quality of building extraction. The overall accuracy and quality of our approach was positively related to buffer distance, with accuracies ranging from 50 to 95% for buffer distances from 0 to 100 m. Our results also indicate that building detection was sensitive to building size, with smaller outbuildings (footprints less than 75 m 2 ) having detection rates below 80% and larger residential buildings having detection rates above 90%. These findings demonstrate that this approach can successfully identify buildings in the WUI in diverse landscapes while achieving high accuracies at buffer distances appropriate for most fire management applications while overcoming cost and time constraints associated with traditional approaches. This study is unique in that it evaluates the ability of an OBIA approach to extract highly detailed data on building locations in a WUI setting.