Search USGS⌕ Search

SEARCH · Search USGS

Results for “The Mountaineer”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31Linked to original sources

Biological and chemical recovery of acidified Catskill Mountain streams in response to the Clean Air Act Amendments of 1990

Decades of acidic deposition have adversely affected aquatic and terrestrial ecosystems in acid-sensitive watersheds in parts of the eastern United States. The national Acid Rain Program (Title IV of the 1990 Clean Air Act Amendments - CAAA) helped reduce emissions of sulfur dioxide (SO 2 ) and nitrogen oxides (NO x ) and resulted in sharp decreases in the acidity of atmospheric deposition. The decrease in acidic deposition produced a steady decline in the acidity of streams in many poorly buffered waters across the western Adirondacks and parts of the Catskill Mountains of New York. Until recently, however, there has been little evidence of biological recovery in most acid-sensitive streams in both regions. Long-term deposition and stream-chemistry records and fish-community data from quantitative surveys done during 1991–93 and again during 2012–19 at 13 sites in the upper Neversink River and its tributaries were evaluated to determine if chemical and biological recovery were evident in this Catskill Mountain watershed and if they could be linked to regional declines in acidic deposition. Between 1991 and 2019, large decreases in sulfate and nitrate deposition in the basin mirrored declines in total nationwide SO 2 and NO x emissions. There were corresponding decreases in sulfate and nitrate concentrations in deposition at a National Trends Network station at Frost Valley (NY68) and coincident declines in sulfate concentrations at four long-term monitoring sites in the Neversink River watershed. Mean acid neutralizing capacity and pH increased and inorganic aluminum (Al i ) concentrations from routine summertime samples decreased significantly at most moderately to severely acidified sites between the two study periods. Richness, density, and biomass of fish communities increased at most sites, while the density and biomass of brook trout Salvelinus fontinalis populations increased at fewer sites that were undergoing chemical recovery. Although recovery is far from complete, trends in deposition chemistry, water quality, and fish assemblages in streams of the upper Neversink watershed indicate that the 1990 CAAA is having positive impacts on aquatic ecosystems in the Catskill Mountain region, New York.

New York↗

Crustal migration of CO2-rich magmatic fluids recorded by tree-ring radiocarbon and seismicity at Mammoth Mountain, CA, USA

Unrest at Mammoth Mountain over the past several decades, manifest by seismicity, ground deformation, diffuse CO 2 emissions, and elevated 3 He/ 4 He ratios in fumarolic gases has been driven by the release of CO 2 -rich fluids from basaltic intrusions in the middle to lower crust. Recent unrest included the occurrence of three lower-crustal (32–19 km depth) seismic swarms beneath Mammoth Mountain in 2006, 2008 and 2009 that were consistently followed by peaks in the occurrence rate of shallow (≤10 km depth) earthquakes. We measured 14 C in the growth rings (1998–2012) of a tree growing in the largest (∼0.3 km 2 ) area of diffuse CO 2 emissions on Mammoth Mountain (the Horseshoe Lake tree kill; HLTK) and applied atmospheric CO 2 concentration source area modeling to confirm that the tree was a reliable integrator of magmatic CO 2 emissions over most of this area. The tree-ring 14 C record implied that magmatic CO 2 emissions from the HLTK were relatively stable from 1998 to 2009, nearly doubled from 2009 to 2011, and then declined by the 2012 growing season. The initial increase in CO 2 emissions was detected during the growing season that immediately followed the largest (February 2010) peak in the occurrence rate of shallow earthquakes. Migration of CO 2 -rich magmatic fluids may have driven observed patterns of elevated deep, then shallow seismicity, while the relationship between pore fluid pressures within a shallow (upper 3 km of crust) fluid reservoir and permeability structure of the reservoir cap rock may have controlled the temporal pattern of surface CO 2 emissions.

California↗

Unraveling the dynamics of magmatic CO2 degassing at Mammoth Mountain, California

The accumulation of magmatic CO 2 beneath low-permeability barriers may lead to the formation of CO 2 -rich gas reservoirs within volcanic systems. Such accumulation is often evidenced by high surface CO 2 emissions that fluctuate over time. The temporal variability in surface degassing is believed in part to reflect a complex interplay between deep magmatic degassing and the permeability of degassing pathways. A better understanding of the dynamics of CO 2 degassing is required to improve monitoring and hazards mitigation in these systems. Owing to the availability of long-term records of CO 2 emissions rates and seismicity, Mammoth Mountain in California constitutes an ideal site towards such predictive understanding. Mammoth Mountain is characterized by intense soil CO 2 degassing (up to ∼1000 t d −1 ) and tree kill areas that resulted from leakage of CO 2 from a CO 2 -rich gas reservoir located in the upper ∼4 km. The release of CO 2 -rich fluids from deeper basaltic intrusions towards the reservoir induces seismicity and potentially reactivates faults connecting the reservoir to the surface. While this conceptual model is well-accepted, there is still a debate whether temporally variable surface CO 2 fluxes directly reflect degassing of intrusions or variations in fault permeability. Here, we report the first large-scale numerical model of fluid and heat transport for Mammoth Mountain. We discuss processes (i) leading to the initial formation of the CO 2 -rich gas reservoir prior to the occurrence of high surface CO 2 degassing rates and (ii) controlling current CO 2 degassing at the surface. Although the modeling settings are site-specific, the key mechanisms discussed in this study are likely at play at other volcanic systems hosting CO 2 -rich gas reservoirs. In particular, our model results illustrate the role of convection in stripping a CO 2 -rich gas phase from a rising hydrothermal fluid and leading to an accumulation of a large mass of CO 2 (∼10 7 –10 8 t) in a shallow gas reservoir. Moreover, we show that both, short-lived (months to years) and long-lived (hundreds of years) events of magmatic fluid injection can lead to critical pressures within the reservoir and potentially trigger fault reactivation. Our sensitivity analysis suggests that observed temporal fluctuations in surface degassing are only indirectly controlled by variations in magmatic degassing and are mainly the result of temporally variable fault permeability. Finally, we suggest that long-term CO 2 emission monitoring, seismic tomography and coupled thermal–hydraulic–mechanical modeling are important for CO 2 -related hazard mitigation.

California↗

Forest disturbance interactions and successional pathways in the Southern Rocky Mountains

The pine forests in the southern portion of the Rocky Mountains are a heterogeneous mosaic of disturbance and recovery. The most extensive and intensive stress and mortality are received from human activity, fire, and mountain pine beetles (MPB; Dendroctonus ponderosae ). Understanding disturbance interactions and disturbance-succession pathways are crucial for adapting management strategies to mitigate their impacts and anticipate future ecosystem change. Driven by this goal, we assessed the forest disturbance and recovery history in the Southern Rocky Mountains Ecoregion using a 13-year time series of Landsat image stacks. An automated classification workflow that integrates temporal segmentation techniques and a random forest classifier was used to examine disturbance patterns. To enhance efficiency in selecting representative samples at the ecoregion scale, a new sampling strategy that takes advantage of the scene-overlap among adjacent Landsat images was designed. The segment-based assessment revealed that the overall accuracy for all 14 scenes varied from 73.6% to 92.5%, with a mean of 83.1%. A design-based inference indicated the average producer’s and user’s accuracies for MPB mortality were 85.4% and 82.5% respectively. We found that burn severity was largely unrelated to the severity of pre-fire beetle outbreaks in this region, where the severity of post-fire beetle outbreaks generally decreased in relation to burn severity. Approximately half the clear-cut and burned areas were in various stages of recovery, but the regeneration rate was much slower for MPB-disturbed sites. Pre-fire beetle outbreaks and subsequent fire produced positive compound effects on seedling reestablishment in this ecoregion. Taken together, these results emphasize that although multiple disturbances do play a role in the resilience mechanism of the serotinous lodgepole pine, the overall recovery could be slow due to the vast area of beetle mortality.

Colorado, New Mexico, Wyoming↗

Historical dominance of low-severity fire in dry and wet mixed-conifer forest habitats of the endangered terrestrial Jemez Mountains salamander ( Plethodon neomexicanus )

Anthropogenic alteration of ecosystem processes confounds forest management and conservation of rare, declining species. Restoration of forest structure and fire hazard reduction are central goals of forest management policy in the western United States, but restoration priorities and treatments have become increasingly contentious. Numerous studies have documented changes in fire regimes, forest stand structure and species composition following a century of fire exclusion in dry, frequent-fire forests of the western U.S. (e.g., ponderosa pine and dry mixed-conifer). In contrast, wet mixed-conifer forests are thought to have historically burned infrequently with mixed- or high-severity fire—resulting in reduced impacts from fire exclusion and low restoration need—but data are limited. In this study we quantified the current forest habitat of the federally endangered, terrestrial Jemez Mountains salamander (Plethodon neomexicanus) and compared it to dendroecological reconstructions of historical habitat (e.g., stand structure and composition), and fire regime parameters along a gradient from upper ponderosa pine to wet mixed-conifer forests. We found that current fire-free intervals in Jemez Mountains salamander habitat (116–165 years) are significantly longer than historical intervals, even in wet mixed-conifer forests. Historical mean fire intervals ranged from 10 to 42 years along the forest gradient. Low-severity fires were historically dominant across all forest types (92 of 102 fires). Although some mixed- or highseverity fire historically occurred at 67% of the plots over the last four centuries, complete mortality within 1.0 ha plots was rare, and asynchronous within and among sites. Climate was an important driver of temporal variability in fire severity, such that mixed- and high-severity fires were associated with more extreme drought than low-severity fires. Tree density in dry conifer forests historically ranged from open (90 trees/ha) to moderately dense (400 trees/ha), but has doubled on average since fire exclusion. Infill of fire-sensitive tree species has contributed to the conversion of historically dry mixedconifer to wet mixed-conifer forest. We conclude that low-severity fire, which has been absent for over a century, was a critical ecosystem process across the forest gradient in Jemez Mountains salamander habitat, and thus is an important element of ecosystem restoration, resilience, and rare species recovery.

Forest Ecology and Management↗

Using stream-side groundwater discharge for geochemical exploration in mountainous terrain

Groundwater chemistry has been predominantly used in geochemical exploration studies to identify mineralized targets concealed under transported cover in areas with gentle topography. Another potentially valuable ap-plication that has received little attention is using groundwater chemistry to identify deposits concealed within mountain ridges. A number of geochemical exploration studies have employed surface water chemistry in hilly and mountainous terrain and have demonstrated mixed success. However, groundwater chemistry is potentially a more refined and powerful tool because groundwater sample locations can be selected to target a specific hillslope zone, whereas surface water integrates the geochemical signature of the entire watershed area above the sample site. In this study, we propose an approach that utilizes shallow groundwater samples collected from near-stream, hand-installed sampling points. We tested this approach in Handcart Gulch, an unmined catchment in the Montezuma Mining District, Colorado, USA, where previous studies suggested the possible presence of a concealed mineral occurrence. Seventeen near-stream mini-boreholes were drilled in a stream-parallel 3.3-km-long transect to depths of ~2 m in ferricrete (iron-oxide cemented colluvium and alluvium) using a handheld rock-coring drill. Groundwater in the mini-boreholes was sampled for major and trace element chemistry, as well as environmental tracers including stable water isotopes, dissolved noble gases, helium isotopes, and tritium, to identify possible variations in groundwater age and recharge conditions/sources. Chemistry data delineate a well-defined 900-m-wide zone of elevated SO4, Cu, and other metal concentrations on the east side of the stream. In this zone, Cu concentrations up to 2670 μg/L are 1–3 orders of magnitude greater than concentrations in upstream and downstream areas. Environmental tracer data show no corresponding anomalies in groundwater age or recharge conditions/sources within the zone, suggesting that the source of the elevated metals is a change in bedrock mineralogy as opposed to other hydrologic factors. The groundwater samples therefore successfully identify and refine the location of a potential concealed high-Cu mineral occurrence in the watershed, de-monstrating that near-stream groundwater discharge may be a valuable medium for geochemical exploration in mountainous areas.

Colorado↗

Parsing complex terrain controls on mountain glacier response to climate forcing

Glaciers are a key indicator of changing climate in the high mountain landscape. Glacier variations across a mountain range are ultimately driven by regional climate forcing. However, changes also reflect local, topographically driven processes such as snow avalanching, snow wind-drifting, and radiation shading as well as the initial glacier conditions such as hypsometry and ice thickness. Here we assess the role of these various terrain influences on change to Little Ice Age (LIA) glaciers in Glacier National Park, U.S.A . With available data for LIA and modern glacier areas, we estimate glacier volumes using simple ice flow assumptions, and topographically driven processes using terrain proxies. At the LIA glacial maxima there were 82 glaciers larger than 0.1 km 2 ranging from 0.11 to 4.97 km 2 . Over the course of the 20 th century, every single LIA glacier decreased in area and 60% (49 glaciers) diminished to below the 0.1 km 2 threshold. Glaciers with large initial area (>1.5 km 2 ) at the end of LIA persisted. Within the intermediate size class (0.5 km 2 < area < 1.5 km 2 ), LIA glacier persistence is poorly explained by initial glacier volume, ice thickness, or elevation. Instead, wind exposure is an important explanatory factor. Our analysis demonstrates the complex response of cirque glaciers to post-LIA climate change in this region: individual glaciers have not necessarily undergone equivalent and synchronous change. Nevertheless, that all glaciers in this mountain range experienced retreat demonstrates that local processes mediated adjustments of some glaciers, but completely decoupled none from the regional climate forcing.

Montana↗

Magmatic unrest beneath Mammoth Mountain, California

Mammoth Mountain, which stands on the southwest rim of Long Valley caldera in eastern California, last erupted ∼57,000 years BP. Episodic volcanic unrest detected beneath the mountain since late 1979, however, emphasizes that the underlying volcanic system is still active and capable of producing future volcanic eruptions. The unrest symptoms include swarms of small ( M ≤ 3) earthquakes, spasmodic bursts (rapid-fire sequences of brittle-failure earthquakes with overlapping coda), long-period (LP) and very-long-period (VLP) volcanic earthquakes, ground deformation, diffuse emission of magmatic CO 2 , and fumarole gases with elevated 3 He/ 4 He ratios. Spatial-temporal relations defined by the multi-parameter monitoring data together with earthquake source mechanisms suggest that this Mammoth Mountain unrest is driven by the episodic release of a volume of CO 2 -rich hydrous magmatic fluid derived from the upper reaches of a plexus of basaltic dikes and sills at mid-crustal depths (10–20 km). As the mobilized fluid ascends through the brittle–plastic transition zone and into overlying brittle crust, it triggers earthquake swarm activity and, in the case of the prolonged, 11-month-long earthquake swarm of 1989, crustal deformation and the onset of diffuse CO 2 emissions. Future volcanic activity from this system would most likely involve steam explosions or small-volume, basaltic, strombolian or Hawaiian style eruptions. The impact of such an event would depend critically on vent location and season.

California↗

Long-term year-round observations of magmatic CO2 emissions on Mammoth Mountain, California, USA

Diffuse emission of magmatic CO 2 is one of the main indicators of volcanic unrest at Mammoth Mountain, but the presence of deep seasonal snowpack at the site has hindered year-round CO 2 flux observations. A permanent eddy covariance station was established at the largest area of diffuse CO 2 degassing on Mammoth Mountain (Horseshoe Lake tree kill) that measured CO 2 fluxes ( F c ) and meteorological parameters on a half-hourly basis. From July 22, 2014 to May 24, 2020, F c ranged from −35 to 10,546 g m −2 d −1 . F c decreased on average by 53% over the study period, tracking the long-term decline in CO 2 emissions following the last major increase that occurred at the Horseshoe Lake tree kill area from 2009 to 2011. Statistical and spectral analyses were applied to the F c and ancillary meteorological parameter time series to understand (1) relationships between these parameters, (2) their dominant periodicities, and (3) changes in F c that may be unexplained by meteorological forcing. Variations in detrended F c ( F cdt ) were most strongly correlated with wind direction and atmospheric temperature, followed by atmospheric pressure on diurnal to annual time scales, but wind direction likely exerted the most direct control on F cdt . Comparison of the smoothed (180-d span) F cdt time series to the time series of average-daily snow water equivalent measured ~1 km away suggested that snowpack may have suppressed CO 2 emissions. No evidence of a change in CO 2 emissions related to the last major seismic swarm beneath Mammoth Mountain on February 2–18, 2014 was observed.

California↗

Paleomagnetism of Proterozoic mafic dikes from the Tobacco Root Mountains, southwest Montana

Paleomagnetic data from Proterozoic mafic dikes in southwestern Montana provides evidence for two distinct episodes of subparallel dike emplacement at ca. 1450 and 780 Ma. Published geochemical data from dikes in the southern Tobacco Root Mountains has identified three distinct compositional groups, termed groups A, B, and C. Geochronological data from the group A dikes yielded a Sm-Nd age of 1448 ?? 49 Ma. Emplacement of these dikes is thought to reflect mafic magmatism associated with extension accompanying development of the adjacent Mesoproterozoic Belt Basin. Paleomagnetic results from these dikes and a group C dike yield antipodal magnetizations with a group-mean direction of D = 225.0??, I = 61.8?? (k = 27.9, ??95 = 7.7??, N = 14 independent means/24 sites). The average paleomagnetic pole (8.7??N, 216.1??E, A95 = 10.3??) is considered to be primary on the basis of positive baked contact tests and similarity to poles of ca. 1.45-1.4 Ga from intrusions elsewhere in North America, but is discordant with respect to poles from age equivalent sedimentary rocks of the Meosoproterozoic Belt Supergroup. 40Ar/39Ar dates from geochemical group B dikes are consistent with published U-Pb dates that demonstrate dike emplacement at 780 Ma as part of the regional Gunbarrel magmatic event. Hornblende concentrates from the group B dikes yield 40Ar/39Ar apparent ages of 778-772 Ma, whereas biotite from a baked contact zone yielded a plateau date of 788 Ma. Paleomagnetic results from the group B dikes yield a mean direction of D = 301.5??, I = -17.1?? (k = 65.7, ??95 = 4.0??, N = 12 independent means/23 sites) with a paleomagnetic pole at 14.6??N, 127.0??E (A95 = 3.2??). The combination of geochronologic data, results of a baked contact test, and spatial agreement of the paleomagnetic poles with poles of similar age elsewhere in North America indicates that this is also a primary magnetization associated with dike emplacement. Paleomagnetic data from some of the Tobacco Root Mountains dikes provide evidence that they were partially to completely remagnetized during latest Cretaceous to early Tertiary time, perhaps due to thermal affects associated with emplacement of the Late Cretaceous Tobacco Root Batholith. The overall agreement of paleomagnetic poles from the Proterozoic dikes with those of age equivalent rocks elsewhere in North America and agreement of the secondary magnetization with expected directions for the latest Cretaceous/early Tertiary indicate that the rocks of the Tobacco Root Mountains have not experienced significant tilting or vertical axis rotation since the Mesoproterozoic. The new paleomagnetic poles from this study thus provide key data for refining Meso- and Neoproterozoic parts of the North American APW path. ?? 2008 Elsevier B.V. All rights reserved.

Precambrian Research↗

The geochemistry and petrogenesis of the Paleoproterozoic Green Mountain arc: A composite(?), bimodal, oceanic, fringing arc

The inferred subduction affinity of the ∼1780-Ma Green Mountain arc, a dominantly bimodal igneous terrane (together with immature marine and volcaniclastic sedimentary rocks) accreted to the southern margin of the Wyoming province, is integral to arc-accretion models of the Paleoproterozoic growth of southern Laurentia. Conversely, the dominantly bimodal nature of many putative arc-related igneous suites throughout southern Laurentia, including the Green Mountain arc, has also been used to support models of growth by extension of pre-existing crust. We report new geochemical and isotopic data from ∼1780-Ma gabbroic and granodioritic to tonalitic rocks of the Big Creek Gneiss, interpreted as consanguineous with previously studied metavolcanic rocks of the Green Mountain Formation. The ∼1780-Ma Big Creek Gneiss mafic rocks show clear geochemical signatures of a subduction origin and provide no supporting evidence for extensional tectonism. The ∼1780-Ma Big Creek Gneiss felsic rocks are attributed to partial melting of mafic and/or mixed lower-crustal material. The bimodal nature of the suite results from the combination of arc basalts and felsic crustal melts. The lack of andesite is consistent with the observed tholeiitic differentiation trend of the mafic magmas. The lower ɛ Nd (1780 Ma) values for the felsic rocks vs. the mafic rocks suggest that the unexposed lower crust of the arc may be older than the arc and that Trans-Hudson- or Penokean-aged rocks possibly form the substratum of the arc. Our results reinforce previous interpretations that arc-related magmatism played a key role in the Paleoproterozoic crustal growth of southern Laurentia, but also support the possibility of unexposed older crust as basement to the arcs.

Wyoming, Colorado↗

Millennial-scale climatic and cultural impacts on vegetation and fire at the southern edge of the Rocky Mountains, USA

Mountain forests and woodlands in semiarid regions of the world are threatened by climatic change and other human impacts. In the southwestern USA, climate and culturally driven changes to the structure and fire regimes of dry coniferous forests over recent centuries are well documented by tree-ring archives. However, the roles of climate and people as drivers of millennial-scale changes are less understood. We present a new record of vegetation dynamics and regional fire activity inferred from pollen, microscopic charcoal, and sediment geochemistry from Santa Fe Lake, NM (3532 m a.s.l.), the southernmost natural lake in the Rocky Mountains. Declining elemental intensities (e.g., Ti, Fe, Si) record local deglaciation after 14,000 cal yr BP followed by upslope expansion of trees indicated by increasing Pinus and Picea pollen percentages as the climate warmed during the Late Glacial. Cool, dry growing seasons probably limited the establishment of dense forests until after 12,200 cal yr BP when we hypothesize that Pinus ponderosa (ponderosa pine) expanded regionally (i.e., within 10–100 km) due to continued warming and a possible increase in monsoonal precipitation. Subalpine forests established near Santa Fe Lake by 10,300 cal yr BP as indicated by abundant conifer stomates and increasing Picea and Pinus aristata pollen percentages, which are highest after 5100 cal yr BP. Increasing Juniperus -type, Quercus , and Pinus edulis -type pollen record changes in vegetation belts at lower elevations after 10,300 cal yr BP. Maize pollen, an unambiguous indicator of Puebloan agriculture, first appears at 800 cal yr BP (1150 cal yr CE) and is present through 1600 cal yr CE when Spanish colonization brought cultural upheaval and population loss to local indigenous communities. Regional fire activity, inferred from microscopic charcoal influx, remained relatively constant for most of our record. However, average microscopic charcoal influx more than doubled relative to the rest of the Holocene during the last 2000 years, with the highest influx after 1550 cal yr BP (400 cal yr CE). The increased influx of microscopic charcoal coincides with archaeological evidence for dramatic regional population growth in the northern Rio Grande region and probably reflects extensive cultural burning associated with the development of an agricultural economy. Although climate was an important driver of vegetation change and fire regimes in the mountains of the southwestern USA during the Holocene, our results suggest that human impacts on regional burning during the last 2000 years far exceeded earlier climate-driven variation.

Rocky Mountains↗

Mapping variations in weight percent silica measured from multispectral thermal infrared imagery - Examples from the Hiller Mountains, Nevada, USA and Tres Virgenes-La Reforma, Baja California Sur, Mexico

Remotely sensed multispectral thermal infrared (8-13 ??m) images are increasingly being used to map variations in surface silicate mineralogy. These studies utilize the shift to longer wavelengths in the main spectral feature in minerals in this wavelength region (reststrahlen band) as the mineralogy changes from felsic to mafic. An approach is described for determining the amount of this shift and then using the shift with a reference curve, derived from laboratory data, to remotely determine the weight percent SiO2 of the surface. The approach has broad applicability to many study areas and can also be fine-tuned to give greater accuracy in a particular study area if field samples are available. The approach was assessed using airborne multispectral thermal infrared images from the Hiller Mountains, Nevada, USA and the Tres Virgenes-La Reforma, Baja California Sur, Mexico. Results indicate the general approach slightly overestimates the weight percent SiO2 of low silica rocks (e.g. basalt) and underestimates the weight percent SiO2 of high silica rocks (e.g. granite). Fine tuning the general approach with measurements from field samples provided good results for both areas with errors in the recovered weight percent SiO2 of a few percent. The map units identified by these techniques and traditional mapping at the Hiller Mountains demonstrate the continuity of the crystalline rocks from the Hiller Mountains southward to the White Hills supporting the idea that these ranges represent an essentially continuous footwall block below a regional detachment. Results from the Baja California data verify the most recent volcanism to be basaltic-andesite. ?? 2005 Elsevier Inc. All rights reserved.

Remote Sensing of Environment↗

Detecting subtle change from dense landsat time series: Case studies of mountain pine beetle and spruce beetle disturbance

In contrast to abrupt changes caused by land cover conversion, subtle changes driven by a shift in the condition, structure, or other biological attributes of land often lead to minimal and slower alterations of the terrestrial surface. Accurate mapping and monitoring of subtle change are crucial for an early warning of long-term gradual change that may eventually result in land cover conversion. Freely accessible moderate-resolution datasets such as the Landsat archive have great potential to characterize subtle change by capturing low-magnitude spectral changes in long-term observations. However, past studies have reported limited success in accurately extracting subtle changes from satellite-based time series analysis. In this study, we introduce a supervised framework named ‘PIDS’ to detect subtle forest disturbance from a comprehensive Landsat data archive by leveraging disturbance-based calibration sites. PIDS consists of four components: (1) Parameter optimization; (2) Index selection; (3) Dynamic stratified monitoring; and (4) Spatial consideration. PIDS was applied to map the early stage of bark beetle infestations (i.e., a lower per-pixel fraction of trees cover that show visual signs of infestation), which are a typical example of subtle change in conifer forests. Landsat Analysis Ready Data were used as the time series inputs for mapping mountain pine beetle and spruce beetle disturbance between 2001 and 2019 in Colorado, USA. PIDS-detection map assessment showed that the overall performance of PIDS (namely ‘F1 score’) was 0.86 for mountain pine beetle and 0.73 for spruce beetle, making a substantial improvement (> 0.3) compared to other approaches/products including COntinuous monitoring of Land Disturbance, LandTrendr, and the National Land Cover Database forest disturbance product. A sub-pixel analysis of tree canopy mortality percentage was performed by linking classified high-resolution (0.3- and 1-m) aerial imagery and 30-m PIDS-detection maps. Results show that PIDS typically detects mountain pine beetle infestation when ≥56% of a Landsat pixel is occupied by red-stage canopy mortality (one year after initial infestation), and spruce beetle infestation when ≥55% is occupied by gray-stage mortality (two years after initial infestation). This study addresses an important methodological goal pertinent to the utility of event-based reference samples for detecting subtle forest change, which could be potentially applied to other types of subtle land change.

Colorado↗

Uncertainty of ICESat-2 ATL06- and ATL08-derived snow depths for glacierized and vegetated mountain regions

Seasonal snow melt dominates the hydrologic budget across a large portion of the globe. Snow accumulation and melt vary over a broad range of spatial scales, preventing accurate extrapolation of sparse in situ observations to watershed scales. The lidar onboard the Ice, Cloud, and land Elevation, Satellite (ICESat-2) was designed for precise mapping of ice sheets and sea ice, and here we assess the feasibility of snow depth-mapping using ICESat-2 data in more complex and rugged mountain landscapes. We explore the utility of ATL08 Land and Vegetation Height and ATL06 Land Ice Height differencing from reference elevation datasets in two end member study sites. We analyze ∼3 years of data for Reynolds Creek Experimental Watershed in Idaho's Owyhee Mountains and Wolverine Glacier in southcentral Alaska's Kenai Mountains. Our analysis reveals decimeter-scale uncertainties in derived snow depth and glacier mass balance at the watershed scale. Both accuracy and precision decrease as slope increases: the magnitudes of the median and median of the absolute deviation of elevation errors (MAD) vary from ∼0.2 m for slopes <5° to >1 m for slopes >20°. For glacierized regions, failure to account for intra- and inter-annual evolution of glacier surface elevations can strongly bias ATL06 elevations, resulting in under-estimation of the mass balance gradient with elevation. Based on these results, we conclude that ATL08 and ATL06 observations are best suited for characterization of watershed-scale snow depth and mass balance gradients over relatively shallow slopes with thick snowpacks . In these regions, ICESat-2 elevation residual-derived snow depth and mass balance transects can provide valuable watershed scale constraints on terrain parameter- and model-derived estimates of snow accumulation and melt.

Alaska, Idaho↗

206Pb-230Th-234U-238U and 207Pb-235U geochronology of Quaternary opal, Yucca Mountain, Nevada

U–Th–Pb isotopic systems have been studied in submillimeter-thick outermost layers of Quaternary opal occurring in calcite–silica fracture and cavity coatings within Tertiary tuffs at Yucca Mountain, Nevada, USA. These coatings preserve a record of paleohydrologic conditions at this site, which is being evaluated as a potential high-level nuclear waste repository. The opal precipitated from groundwater is variably enriched in 234 U (measured 234 U/ 238 U activity ratio 1.124–6.179) and has high U (30–313 ppm), low Th (0.008–3.7 ppm), and low common Pb concentrations (measured 206 Pb/ 204 Pb up to 11,370). It has been demonstrated that the laboratory acid treatment used in this study to clean sample surfaces and to remove adherent calcite, did not disturb U–Th–Pb isotopic systems in opal. The opal ages calculated from 206 Pb∗/ 238 U and 207 Pb∗/ 235 U ratios display strong reverse discordance because of excess radiogenic 206 Pb∗ derived from the elevated initial 234 U. The data are best interpreted using projections of a new four-dimensional concordia diagram defined by 206 Pb∗/ 238 U, 207 Pb∗/ 235 U, 234 U/ 238 U activity , and 230 Th/ 238 U activity . Ages and initial 234 U/ 238 U activity ratios have been calculated using different projections of this diagram and tested for concordance. The data are discordant, that is observed 207 Pb∗/ 235 U ages of 170 ± 32 (2σ) to 1772 ± 40 ka are systematically older than 230 Th/U ages of 34.1 ± 0.6 to 452 ± 32 ka. The age discordance is not a result of migration of uranium and its decay products under the open system conditions, but a consequence of noninstantaneous growth of opal. Combined U–Pb and 230 Th/U ages support the model of slow mineral deposition at the rates of millimeters per million years resulting in layering on a scale too fine for mechanical sampling. In this case, U–Pb ages provide more accurate estimates of the average age for mixed multiage samples than 230 Th/U ages, because ages based on shorter-lived isotopes are nonlinearly biased by younger mineral additions. Use of the combined U–Th–Pb technique to date Yucca Mountain Quaternary opals significantly extends the age range beyond that of the 230 Th/U dating method and shows that selected fracture pathways in the unsaturated zone felsic tuffs of Yucca Mountain have been active throughout the Quaternary.

Geochimica et Cosmochimica Acta↗

The saturated zone at Yucca Mountain: An overview of the characterization and assessment of the saturated zone as a barrier to potential radionuclide migration

The US Department of Energy is pursuing Yucca Mountain, Nevada, for the development of a geologic repository for the disposal of spent nuclear fuel and high-level radioactive waste, if the repository is able to meet applicable radiation protection standards established by the US Nuclear Regulatory Commission and the US Environmental Protection Agency (EPA). Effective performance of such a repository would rely on a number of natural and engineered barriers to isolate radioactive waste from the accessible environment. Groundwater beneath Yucca Mountain is the primary medium through which most radionuclides might move away from the potential repository. The saturated zone (SZ) system is expected to act as a natural barrier to this possible movement of radionuclides both by delaying their transport and by reducing their concentration before they reach the accessible environment. Information obtained from Yucca Mountain Site Characterization Project activities is used to estimate groundwater flow rates through the site-scale SZ flow and transport model area and to constrain general conceptual models of groundwater flow in the site-scale area. The site-scale conceptual model is a synthesis of what is known about flow and transport processes at the scale required for total system performance assessment of the site. This knowledge builds on and is consistent with knowledge that has accumulated at the regional scale but is more detailed because more data are available at the site-scale level. The mathematical basis of the site-scale model and the associated numerical approaches are designed to assist in quantifying the uncertainty in the permeability of rocks in the geologic framework model and to represent accurately the flow and transport processes included in the site-scale conceptual model. Confidence in the results of the mathematical model was obtained by comparing calculated to observed hydraulic heads, estimated to measured permeabilities, and lateral flow rates calculated by the site-scale model to those calculated by the regional-scale flow model. In addition, it was confirmed that the flow paths leaving the region of the potential repository are consistent with those inferred from gradients of measured head and those independently inferred from water-chemistry data. The general approach of the site-scale SZ flow and transport model analysis is to calculate unit breakthrough curves for radionuclides at the interface between the SZ and the biosphere using the three-dimensional site-scale SZ flow and transport model. Uncertainties are explicitly incorporated into the site-scale SZ flow and transport abstractions through key parameters and conceptual models. ?? 2002 Elsevier Science B.V. All rights reserved.

Journal of Contaminant Hydrology↗

The site-scale saturated zone flow model for Yucca Mountain: Calibration of different conceptual models and their impact on flow paths

This paper presents several different conceptual models of the Large Hydraulic Gradient (LHG) region north of Yucca Mountain and describes the impact of those models on groundwater flow near the potential high-level repository site. The results are based on a numerical model of site-scale saturated zone beneath Yucca Mountain. This model is used for performance assessment predictions of radionuclide transport and to guide future data collection and modeling activities. The numerical model is calibrated by matching available water level measurements using parameter estimation techniques, along with more informal comparisons of the model to hydrologic and geochemical information. The model software (hydrologic simulation code FEHM and parameter estimation software PEST) and model setup allows for efficient calibration of multiple conceptual models. Until now, the Large Hydraulic Gradient has been simulated using a low-permeability, east–west oriented feature, even though direct evidence for this feature is lacking. In addition to this model, we investigate and calibrate three additional conceptual models of the Large Hydraulic Gradient, all of which are based on a presumed zone of hydrothermal chemical alteration north of Yucca Mountain. After examining the heads and permeabilities obtained from the calibrated models, we present particle pathways from the potential repository that record differences in the predicted groundwater flow regime. The results show that Large Hydraulic Gradient can be represented with the alternate conceptual models that include the hydrothermally altered zone. The predicted pathways are mildly sensitive to the choice of the conceptual model and more sensitive to the quality of calibration in the vicinity on the repository. These differences are most likely due to different degrees of fit of model to data, and do not represent important differences in hydrologic conditions for the different conceptual models.

Journal of Contaminant Hydrology↗