Geology topics

Research about Montana

Source-linked reports with geographic coverage including Montana.

At least 343 records · Page 19Linked to original sources

USGS2011

Audiomagnetotelluric data to characterize the Revett-type copper-silver deposits at Rock Creek in the Cabinet Mountains Wilderness, Montana

The Revett-type deposits at Rock Creek are part of the concealed stratabound copper-silver deposits located in the Cabinet Mountains Wilderness of Montana. The U.S. Geological Survey is conducting a series of multidisciplinary studies as part of the Assessment Techniques for Concealed Mineral Resources project. Geologic, geochemical, geophysical, and mineral resources data are being evaluated with existing and new mineral deposit models to predict the possibility and probability of undiscovered deposits in covered terranes. To help characterize the size, resistivity, and depth of the mineral deposit concealed beneath thick overburden, a regional southwest-northeast audiomagnetotelluric sounding profile was acquired. Further studies will attempt to determine if induced polarization parameters can be extracted from the magnetotelluric data to determine the size of the mineralized area. The purpose of this report is to release the audiomagnetotelluric sounding data collected along that southwest-northeast profile. No interpretation of the data is included.

Montana↗

USGS2011

Investigation of pier scour in coarse-bed streams in Montana, 2001 through 2007

A primary goal of ongoing field research of bridge scour is improvement of scour-prediction equations so that pier-scour depth is predicted accurately-an important element of hydraulic analysis and design of highway bridges that cross streams, rivers, and other waterways. Scour depth for piers in streambeds with a mixture of sand, gravel, cobbles, and boulders (coarse-bed streams, which are common in Montana) generally is less than the scour depth in finer-grained (sandy) streambeds under similar conditions. That difference is attributed to an armor layer of coarser material. Pier-scour data from the U.S. Geological Survey were used in this study to develop a bed-material correction factor, which was incorporated into the Federal Highway Administration's recommended equation for computing pier scour. This report describes results of a study of pier scour in coarse-bed streams at 59 bridge sites during 2001-2007 in the mountain and foothill regions of western Montana. Respective drainage areas ranged from about 3 square miles (mi 2 ) to almost 20,000 mi 2 . Data collected and analyzed for this study included 103 pier-scour measurements; the report further describes data collection, shows expansion of the national coarse pier-scour database, discusses use of the new data in evaluation of relative accuracy of various predictive equations, and demonstrates how differences in size and gradation between surface bed material and shallow-subsurface bed material might relate to pier scour. Nearly all measurements were made under clear-water conditions with no incoming sediment supply to the bridge opening. Half of the measurements showed approach velocities that equaled or surpassed the critical velocity for incipient motion of bed material, possibly indicating that measurements were made very near the threshold between clear-water and live-bed scour, where maximum scour was shown in laboratory studies. Data collected in this study were compared to selected pier-scour data from the nationwide Bridge Scour Data Management System (BSDMS), to show the effect of bed-material size and gradation on scour depth. Unsteady field flow conditions and armoring by coarser material reduced scour relative to the clear-water/sandy-bed laboratory results at steady flow. The new correction factor and the standard scour equation produced the most accurate estimates of scour depth in armored, coarse-bed conditions. Maximum relative scour occurred at similar velocity across variations in bed material and gradation. Pier scour decreased with increased variation in particle size and gradation.

Montana↗

USGS2011

A digital model for planning water management at Benton Lake National Wildlife Refuge, west-central Montana

Benton Lake National Wildlife Refuge is an important area for waterfowl production and migratory stopover in west-central Montana. Eight wetland units covering about 5,600 acres are the essential features of the refuge. Water availability for the wetland units can be uncertain owing to the large natural variations in precipitation and runoff and the high cost of pumping supplemental water. The U.S. Geological Survey, in cooperation with the U.S. Fish and Wildlife Service, has developed a digital model for planning water management. The model can simulate strategies for water transfers among the eight wetland units and account for variability in runoff and pumped water. This report describes this digital model, which uses a water-accounting spreadsheet to track inputs and outputs to each of the wetland units of Benton Lake National Wildlife Refuge. Inputs to the model include (1) monthly values for precipitation, pumped water, runoff, and evaporation; (2) water-level/capacity data for each wetland unit; and (3) the pan-evaporation coefficient. Outputs include monthly water volume and flooded surface area for each unit for as many as 5 consecutive years. The digital model was calibrated by comparing simulated and historical measured water volumes for specific test years.

Montana↗

USGS2011

Diet overlap of top-level predators in recent sympatry: bull trout and nonnative lake trout

The establishment of nonnative lake trout Salvelinus namaycush in lakes containing lacustrine–adfluvial bull trout Salvelinus confluentus often results in a precipitous decline in bull trout abundance. The exact mechanism for the decline is unknown, but one hypothesis is related to competitive exclusion for prey resources. We had the rare opportunity to study the diets of bull trout and nonnative lake trout in Swan Lake, Montana during a concomitant study. The presence of nonnative lake trout in Swan Lake is relatively recent and the population is experiencing rapid population growth. The objective of this study was to evaluate the diets of bull trout and lake trout during the early expansion of this nonnative predator. Diets were sampled from 142 bull trout and 327 lake trout during the autumn in 2007 and 2008. Bull trout and lake trout had similar diets, both consumed Mysis diluviana as the primary invertebrate, especially at juvenile stages, and kokanee Oncorhynchus nerka as the primary vertebrate prey, as adults. A diet shift from primarily M. diluviana to fish occurred at similar lengths for both species, 506 mm (476–545 mm, 95% CI) for bull trout and 495 mm (470–518 mm CI) for lake trout. These data indicate high diet overlap between these two morphologically similar top-level predators. Competitive exclusion may be a possible mechanism if the observed overlap remains similar at varying prey densities and availability.

Montana↗

USGS2011

Shared bacterial and viral respiratory agents in bighorn sheep ( Ovis canadensis ), domestic sheep ( Ovis aries ), and goats ( Capra hircus ) in Montana

Transmission of infectious agents from livestock reservoirs has been hypothesized to cause respiratory disease outbreaks in bighorn sheep ( Ovis canadensis ), and land management policies intended to limit this transmission have proven controversial. This cross-sectional study compares the infectious agents present in multiple populations of bighorn sheep near to and distant from their interface with domestic sheep ( O. aries ) and domestic goat ( Capra hircus ) and provides critical baseline information needed for interpretations of cross-species transmission risks. Bighorn sheep and livestock shared exposure to Pasteurellaceae, viral, and endoparasite agents. In contrast, although the impact is uncertain, Mycoplasma sp. was isolated from livestock but not bighorn sheep. These results may be the result of historic cross-species transmission of agents that has resulted in a mosaic of endemic and exotic agents. Future work using longitudinal and multiple population comparisons is needed to rigorously establish the risk of outbreaks from cross-species transmission of infectious agents.

Montana↗

USGS2011

The distribution and abundance ofa nuisance native alga, Didymosphenia geminata,in streams of Glacier National Park: Climate drivers and management implications

Didymosphenia geminata (didymo) is a freshwater alga native to North America, including Glacier National Park, Montana. It has long been considered a cold-water species, but has recently spread to lower latitudes and warmer waters, and increasingly forms large blooms that cover streambeds. We used a comprehensive monitoring data set from the National Park Service (NPS) and USGS models of stream temperatures to explore the drivers of didymo abundance in Glacier National Park. We estimate that approximately 64% of the stream length in the park contains didymo, with around 5% in a bloom state. Results suggest that didymo abundance likely increased over the study period (2007–2009), with blooms becoming more common. Our models suggest that didymo abundance is positively related to summer stream temperatures and negatively related to total nitrogen and the distance downstream from lakes. Regional climate model simulations indicate that stream temperatures in the park will likely continue to increase over the coming decades, which may increase the extent and severity of didymo blooms. As a result, didymo may be a useful indicator of thermal and hydrological modification associated with climate warming, especially in a relatively pristine system like Glacier where proximate human-related disturbances are absent or reduced. Glacier National Park plays an important role as a sentinel for climate change and associated education across the Rocky Mountain region.

Montana↗

USGS2010

Low genotyping error rates and noninvasive sampling in Bighorn Sheep

Noninvasive DNA sampling allows studies of natural populations without disturbing the target animals. Unfortunately, high genotyping error rates often make noninvasive studies difficult. We report low error rates (0.0–7.5%/locus) when genotyping 18 microsatellite loci in only 4 multiplex polymerase chain reaction amplifications using fecal DNA from bighorn sheep ( Ovis canadensis ). The average locus-specific error rates varied significantly between the 2 populations (0.13% vs. 1.6%; P < 0.001), as did multi-locus genotype error rates (2.3% vs. 14.1%; P < 0.007). This illustrates the importance of quantifying error rates in each study population (and for each season and sample preservation method) before initiating a noninvasive study. Our error rates are among the lowest reported for fecal samples collected noninvasively in the field. This and other recent studies suggest that noninvasive fecal samples can be used in species with pellet-form feces for nearly any study (e.g., of population structure, gene flow, dispersal, parentage, and even genome-wide studies to detect local adaptation) that previously required high-quality blood or tissue samples.

Montana↗

USGS2010

Grizzly bear density in Glacier National Park, Montana

We present the first rigorous estimate of grizzly bear ( Ursus arctos ) population density and distribution in and around Glacier National Park (GNP), Montana, USA. We used genetic analysis to identify individual bears from hair samples collected via 2 concurrent sampling methods: 1) systematically distributed, baited, barbed-wire hair traps and 2) unbaited bear rub trees found along trails. We used Huggins closed mixture models in Program MARK to estimate total population size and developed a method to account for heterogeneity caused by unequal access to rub trees. We corrected our estimate for lack of geographic closure using a new method that utilizes information from radiocollared bears and the distribution of bears captured with DNA sampling. Adjusted for closure, the average number of grizzly bears in our study area was 240.7 (95% CI = 202–303) in 1998 and 240.6 (95% CI = 205–304) in 2000. Average grizzly bear density was 30 bears/1,000 km 2 , with 2.4 times more bears detected per hair trap inside than outside GNP. We provide baseline information important for managing one of the few remaining populations of grizzlies in the contiguous United States.

Montana↗

USGS2010

The dispersion of fibrous amphiboles by glacial processes in the area surrounding Libby, Montana, USA

Mining operations began at a world-class vermiculite deposit at Vermiculite Mountain near Libby, Montana, circa 1920 and ended in 1990. Fibrous and asbestiform amphiboles intergrown with vermiculite ore are suspected to be a causative factor in an abnormally high number of cases of respiratory diseases in former mine and mill workers, and in residents of Libby. The question addressed in this report is whether some of the amphibole from Vermiculite Mountain could have been dispersed by Pleistocene glacial processes rather than by human activity after vermiculite mining began. The history of Pinedale glaciation in the Libby area provides a framework for estimating the presence and distribution of asbestiform amphiboles derived from Vermiculite Mountain and found in naturally occurring sediments of Glacial Lake Kootenai that underlie the Libby Valley area. There were two situations where sediments derived from Vermiculite Mountain were deposited into Glacial Lake Kootenai: (1) as lake-bottom sediments derived from meltwater flowing down Rainy Creek when the valley south of Vermiculite Mountain was free of ice but active ice still covered Vermiculite Mountain; and (2) as lake-bottom sediments eroded from the Rainy Creek outwash and re-deposited during a re-advance of the Purcell Trench Glacier lobe near Moyie Springs, Idaho.

Montana↗

USGS2010

Potential economic benefits of adapting agricultural production systems to future climate change

Potential economic impacts of future climate change on crop enterprise net returns and annual net farm income (NFI) are evaluated for small and large representative farms in Flathead Valley in Northwest Montana. Crop enterprise net returns and NFI in an historical climate period (1960–2005) and future climate period (2006–2050) are compared when agricultural production systems (APSs) are adapted to future climate change. Climate conditions in the future climate period are based on the A1B, B1, and A2 CO 2 emission scenarios from the Intergovernmental Panel on Climate Change Fourth Assessment Report. Steps in the evaluation include: (1) specifying crop enterprises and APSs (i.e., combinations of crop enterprises) in consultation with locals producers; (2) simulating crop yields for two soils, crop prices, crop enterprises costs, and NFIs for APSs; (3) determining the dominant APS in the historical and future climate periods in terms of NFI; and (4) determining whether NFI for the dominant APS in the historical climate period is superior to NFI for the dominant APS in the future climate period. Crop yields are simulated using the Environmental/Policy Integrated Climate (EPIC) model and dominance comparisons for NFI are based on the stochastic efficiency with respect to a function (SERF) criterion. Probability distributions that best fit the EPIC-simulated crop yields are used to simulate 100 values for crop yields for the two soils in the historical and future climate periods. Best-fitting probability distributions for historical inflation-adjusted crop prices and specified triangular probability distributions for crop enterprise costs are used to simulate 100 values for crop prices and crop enterprise costs. Averaged over all crop enterprises, farm sizes, and soil types, simulated net return per ha averaged over all crop enterprises decreased 24% and simulated mean NFI for APSs decreased 57% between the historical and future climate periods. Although adapting APSs to future climate change is advantageous (i.e., NFI with adaptation is superior to NFI without adaptation based on SERF), in six of the nine cases in which adaptation is advantageous, NFI with adaptation in the future climate period is inferior to NFI in the historical climate period. Therefore, adaptation of APSs to future climate change in Flathead Valley is insufficient to offset the adverse impacts on NFI of such change.

Montana↗

USGS2010

Avalanche ecology and large magnitude avalanche events: Glacier National Park, Montana, USA

Large magnitude snow avalanches play an important role ecologically in terms of wildlife habitat, vegetation diversity, and sediment transport within a watershed. Ecological effects from these infrequent avalanches can last for decades. Understanding the frequency of such large magnitude avalanches is also critical to avalanche forecasting for the Going-to-the-Sun Road (GTSR). In January 2009, a large magnitude avalanche cycle occurred in and around Glacier National Park, Montana. The study site is the Little Granite avalanche path located along the GTSR. The study is designed to quantify change in vegetative cover immediately after a large magnitude event and document ecological response over a multi-year period. GPS field mapping was completed to determine the redefined perimeter of the avalanche path. Vegetation was inventoried using modified U.S. Forest Service Forest Inventory and Analysis plots, cross sections were taken from over 100 dead trees throughout the avalanche path, and an avalanche chronology was developed. Initial results indicate that the perimeter of this path was expanded by 30%. The avalanche travelled approximately 1200 vertical meters and 3 linear kilometers. Stands of large conifers as old as 150 years were decimated by the avalanche, causing a shift in dominant vegetation types in many parts of the avalanche path. Woody debris is a major ground cover up to 3 m in depth on lower portions of the avalanche path and will likely affect tree regrowth. Monitoring and measuring the post-avalanche vegetation recovery of this particular avalanche path provides a unique dataset for determining the ecological role of avalanches in mountain landscapes.

Montana↗

USGS2010

Characterizing wet slab and glide slab avalanche occurrence along the Going-to-the-Sun Road, Glacier National Park, Montana, USA

Wet slab and glide slab snow avalanches are dangerous and yet can be particularly difficult to predict. Both wet slab and glide slab avalanches are thought to depend upon free water moving through the snowpack but are driven by different processes. In Glacier National Park, Montana, both types of avalanches can occur in the same year and affect the Going-to-the-Sun Road (GTSR). Both wet slab and glide slab avalanches along the GTSR from 2003-2010 are investigated. Meteorological data from two high-elevation weather stations and one SNOTEL site are used in conjunction with an avalanche database and snowpit profiles. These data were used to characterize years when only glide slab avalanches occurred and those years when both glide slab and wet slab avalanches occurred. Results of 168 glide slab and 57 wet slab avalanches along the GTSR suggest both types of avalanche occurrence depend on sustained warming periods with intense solar radiation (or rain on snow) to produce free water in the snowpack. Differences in temperature and net radiation metrics between wet slab and glide slab avalanches emerge as one moves from one day to seven days prior to avalanche occurrence. On average, a more rapid warming precedes wet slab avalanche occurrence. Glide slab and wet slab avalanches require a similar amount of net radiation. Wet slab avalanches do not occur every year, while glide slab avalanches occur annually. These results aim to enhance understanding of the required meteorological conditions for wet slab and glide slab avalanches and aid in improved wet snow avalanche forecasting.

Montana↗

USGS2009

Simulations of cataclysmic outburst floods from Pleistocene Glacial Lake Missoula

Using a flow domain that we constructed from 30 m digital-elevation model data of western United States and Canada and a two-dimensional numerical model for shallow-water flow over rugged terrain, we simulated outburst floods from Pleistocene Glacial Lake Missoula. We modeled a large, but not the largest, flood, using initial lake elevation at 1250 m instead of 1285 m. Rupture of the ice dam, centered on modern Lake Pend Oreille, catastrophically floods eastern Washington and rapidly fills the broad Pasco, Yakima, and Umatilla Basins. Maximum flood stage is reached in Pasco and Yakima Basins 38 h after the dam break, whereas maximum flood stage in Umatilla Basin occurs 17 h later. Drainage of these basins through narrow Columbia gorge takes an additional 445 h. For this modeled flood, peak discharges in eastern Washington range from 10 to 20 × 10 6 m 3 /s. However, constrictions in Columbia gorge limit peak discharges to <6 × 10 6 m 3 /s and greatly extend the duration of flooding. We compare these model results with field observations of scabland distribution and high-water indicators. Our model predictions of the locations of maximum scour (product of bed shear stress and average flow velocity) match the distribution of existing scablands. We compare model peak stages to high-water indicators from the Rathdrum-Spokane valley, Walulla Gap, and along Columbia gorge. Though peak stages from this less-than-maximal flood model attain or exceed peak-stage indicators along Rathdrum-Spokane valley and along Columbia gorge, simulated peak stages near Walulla Gap are 10–40 m below observed peak-stage indicators. Despite this discrepancy, our match to field observations in most of the region indicates that additional sources of water other than Glacial Lake Missoula are not required to explain the Missoula floods.

Montana↗

USGS2009

A shallow subsurface controlled release facility in Bozeman, Montana, USA, for testing near surface CO2 detection techniques and transport models

A controlled field pilot has been developed in Bozeman, Montana, USA, to study near surface CO2 transport and detection technologies. A slotted horizontal well divided into six zones was installed in the shallow subsurface. The scale and CO2 release rates were chosen to be relevant to developing monitoring strategies for geological carbon storage. The field site was characterized before injection, and CO2 transport and concentrations in saturated soil and the vadose zone were modeled. Controlled releases of CO2 from the horizontal well were performed in the summers of 2007 and 2008, and collaborators from six national labs, three universities, and the U.S. Geological Survey investigated movement of CO2 through the soil, water, plants, and air with a wide range of near surface detection techniques. An overview of these results will be presented.

Montana↗

USGS2009

Changes in the chemistry of shallow groundwater related to the 2008 injection of CO2 at the ZERT field site, Bozeman, Montana

Approximately 300 kg/day of food-grade CO 2 was injected through a perforated pipe placed horizontally 2–2.3 m deep during July 9–August 7, 2008 at the MSU-ZERT field test to evaluate atmospheric and near-surface monitoring and detection techniques applicable to the subsurface storage and potential leakage of CO 2 . As part of this multidisciplinary research project, 80 samples of water were collected from 10 shallow monitoring wells (1.5 or 3.0 m deep) installed 1–6 m from the injection pipe, at the southwestern end of the slotted section (zone VI), and from two distant monitoring wells. The samples were collected before, during, and following CO 2 injection. The main objective of study was to investigate changes in the concentrations of major, minor, and trace inorganic and organic compounds during and following CO 2 injection. The ultimate goals were (1) to better understand the potential of groundwater quality impacts related to CO 2 leakage from deep storage operations, (2) to develop geochemical tools that could provide early detection of CO 2 intrusion into underground sources of drinking water (USDW), and (3) to test the predictive capabilities of geochemical codes against field data. Field determinations showed rapid and systematic changes in pH (7.0–5.6), alkalinity (400–1,330 mg/l as HCO 3 ), and electrical conductance (600–1,800 μS/cm) following CO 2 injection in samples collected from the 1.5 m-deep wells. Laboratory results show major increases in the concentrations of Ca (90–240 mg/l), Mg (25–70 mg/l), Fe (5–1,200 ppb), and Mn (5–1,400 ppb) following CO 2 injection. These chemical changes could provide early detection of CO 2 leakage into shallow groundwater from deep storage operations. Dissolution of observed carbonate minerals and desorption-ion exchange resulting from lowered pH values following CO 2 injection are the likely geochemical processes responsible for the observed increases in the concentrations of solutes; concentrations generally decreased temporarily following four significant precipitation events. The DOC values obtained are 5 ± 2 mg/l, and the variations do not correlate with CO 2 injection. CO 2 injection, however, is responsible for detection of BTEX (e.g. benzene, 0–0.8 ppb), mobilization of metals, the lowered pH values, and increases in the concentrations of other solutes in groundwater. The trace metal and BTEX concentrations are all significantly below the maximum contaminant levels (MCLs). Sequential leaching of core samples is being carried out to investigate the source of metals and other solutes.

Montana↗

USGS2009

Channel water balance and exchange with subsurface flow along a mountain headwater stream in Montana, United States

Channel water balances of contiguous reaches along streams represent a poorly understood scale of stream-subsurface interaction. We measured reach water balances along a headwater stream in Montana, United States, during summer base flow recessions. Reach water balances were estimated from series of tracer tests in 13 consecutive reaches delineated evenly along a 2.6 km valley segment. For each reach, we estimated net change in discharge, gross hydrologic loss, and gross hydrologic gain from tracer dilution and mass recovery. Four series of tracer tests were performed during relatively high, intermediate, and low base flow conditions. The relative distribution of channel water along the stream was strongly related to a transition in valley structure, with a general increase in gross losses through the recession. During tracer tests at intermediate and low flows, there were frequent substantial losses of tracer mass (>10%) that could not be explained by net loss in flow over the reach, indicating that many of the study reaches were concurrently losing and gaining water. For example, one reach with little net change in discharge exchanged nearly 20% of upstream flow with gains and losses along the reach. These substantial bidirectional exchanges suggest that some channel interactions with subsurface flow paths were not measurable by net change in flow or transient storage of recovered tracer. Understanding bidirectional channel water balances in stream reaches along valleys is critical to an accurate assessment of stream solute fate and transport and to a full assessment of exchanges between the stream channel and surrounding subsurface.

Montana↗

USGS2009

Structural development of high-temperature mylonites in the Archean Wyoming province, northwestern Madison Range, Montana

The Crooked Creek mylonite, in the northwestern Madison Range, southwestern Montana, is defined by several curved lenses of high non-coaxial strain exposed over a 7-km-wide, northeast-trending strip. The country rocks, part of the Archean Wyoming province, are dominantly trondhjemitic to granitic orthogneiss with subordinate amphibolite, quartzite, aluminous gneiss, and sills of metabasite (mafic granulite). Data presented here support an interpretation that the mylonite formed during a period of rapid, heterogeneous strain at near-peak metamorphic conditions during an early deformational event (D 1 ) caused by northwest–southeast-directed transpression. The mylonite has a well-developed L-S tectonite fabric and a fine-grained, recrystallized (granoblastic) texture. The strong linear fabric, interpreted as the stretching direction, is defined by elongate compositional “fish,” fold axes, aligned elongate minerals, and mullion axes. The margins of the mylonitic zones are concordant with and grade into regions of unmylonitized gneiss. A second deformational event (D 2 ) has folded the mylonite surface to produce meter- to kilometer-scale, tight-to-isoclinal, gently plunging folds in both the mylonite and country rock, and represents a northwest–southeast shortening event. Planar or linear fabrics associated with D 2 are remarkably absent. A third regional deformational event (D 3 ) produced open, kilometer-scale folds generally with gently north-plunging fold axes. Thermobarometric measurements presented here indicate that metamorphic conditions during D 1 were the same in both the mylonite and the country gneiss, reaching upper amphibolite- to lower granulite-facies conditions: 700 ± 50° C and 8.5 ± 0.5 kb. Previous geochronological studies of mylonitic and cross-cutting rocks in the Jerome Rock Lake area, east of the Crooked Creek mylonite, bracket the timing of this high-grade metamorphism and mylonitization between 2.78 and 2.56 Ga, nearly a billion years before the 1.78-Ga Big Sky orogeny, which overprinted the basement rocks exposed in adjacent ranges of the Wyoming province.

Montana↗

USGS2009

Long-duration drought variability and impacts on ecosystem services: A case study from Glacier National Park, Montana

Instrumental climate records suggest that summer precipitation and winter snowpack in Glacier National Park (Glacier NP), Montana, vary significantly over decadal to multidecadal time scales. Because instrumental records for the region are limited to the twentieth century, knowledge of the range of variability associated with these moisture anomalies and their impacts on ecosystems and physical processes are limited. The authors developed a reconstruction of summer (June–August) moisture variability spanning a.d. 1540–2000 from a multispecies network of tree-ring chronologies in Glacier NP. Decadal-scale drought and pluvial regimes were defined as any event lasting 10 yr or greater, and the significance of each potential regime was assessed using intervention analysis. Intervention analysis prevents single intervening years of average or opposing moisture conditions from ending what was otherwise a sustained moisture regime. The reconstruction shows numerous decadal-scale shifts between persistent drought and wet events prior to the instrumental period (before a.d. 1900). Notable wet events include a series of three long-duration, high-magnitude pluvial regimes spanning the end of the Little Ice Age ( a.d. 1770–1840). Though the late-nineteenth century was marked by a series of >10 yr droughts, the single most severe dry event occurred in the early-twentieth century ( a.d. 1917–41). These decadal-scale dry and wet events, in conjunction with periods of high and low snowpack, have served as a driver of ecosystem processes such as forest fires and glacial dynamics in the Glacier NP region. Using a suite of paleoproxy reconstructions and information from previous studies examining the relationship between climate variability and natural processes, the authors explore how such persistent moisture anomalies affect the delivery of vital goods and services provided by Glacier NP and surrounding areas. These analyses show that regional water resources and tourism are particularly vulnerable to persistent moisture anomalies in the Glacier NP area. Many of these same decadal-scale wet and dry events were also seen among a wider network of hydroclimatic reconstructions along a north–south transect of the Rocky Mountains. Such natural climate variability can, in turn, have enormous impacts on the sustainable provision of natural resources over wide areas. Overall, these results highlight the susceptibility of goods and services provided by protected areas like Glacier NP to natural climate variability, and show that this susceptibility will likely be compounded by the effects of future human-induced climate change.

Montana↗