USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earth Interactions”

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 37 records · Page 2Linked to original sources
USGS2011

Simulations of historical and future trends in snowfall and groundwater recharge for basins draining to Long Island Sound

A regional watershed model was developed for watersheds contributing to Long Island Sound, including the Connecticut River basin. The study region covers approximately 40 900 km 2 , extending from a moderate coastal climate zone in the south to a mountainous northern New England climate zone dominated by snowmelt in the north. The input data indicate that precipitation and temperature have been increasing for the last 46 years (1961– 2006) across the region. Minimum temperature has increased more than maximum temperature over the same period (1961–2006). The model simulation indicates that there was an upward trend in groundwater recharge across most of the modeled region. However, trends in increasing precipitation and groundwater recharge are not significant at the 0.05 level if the drought of 1961–67 is removed from the time series. The trend in simulated snowfall is not significant across much of the region, although there is a significant downward trend in southeast Connecticut and in central Massachusetts. To simulate future trends, two input datasets, one assuming high carbon emissions and one assuming low carbon emissions, were developed from GCM forecasts. Under both of the carbon emission scenarios, simulations indicate that historical trends will continue, with increases in groundwater recharge over much of the region and substantial snowfall decreases across Massachusetts, Connecticut, southern Vermont, and southern New Hampshire. The increases in groundwater recharge and decreases in snowfall are most pronounced for the high emission scenario.

Connecticut, Massachusetts, New Hampshire, Rhode I↗
USGS2011

Impacts of climate change on the growing season in the United States

Understanding the effects of climate change on the vegetative growing season is key to quantifying future hydrologic water budget conditions. The U.S. Geological Survey modeled changes in future growing season length at 14 basins across 11 states. Simulations for each basin were generated using five general circulation models with three emission scenarios as inputs to the Precipitation-Runoff Modeling System (PRMS). PRMS is a deterministic, distributed-parameter, watershed model developed to simulate the effects of various combinations of precipitation, climate, and land use on watershed response. PRMS was modified to include a growing season calculation in this study. The growing season was examined for trends in the total length (annual), as well as changes in the timing of onset (spring) and the end (fall) of the growing season. The results showed an increase in the annual growing season length in all 14 basins, averaging 27–47 days for the three emission scenarios. The change in the spring and fall growing season onset and end varied across the 14 basins, with larger increases in the total length of the growing season occurring in the mountainous regions and smaller increases occurring in the Midwest, Northeast, and Southeast regions. The Clear Creek basin, 1 of the 14 basins in this study, was evaluated to examine the growing season length determined by emission scenario, as compared to a growing season length fixed baseline condition. The Clear Creek basin showed substantial variation in hydrologic responses, including streamflow, as a result of growing season length determined by emission scenario.

Earth Interactions↗
USGS2011

Effects of baseline conditions on the simulated hydrologic response to projected climate change

Changes in temperature and precipitation projected from five general circulation models, using one late-twentieth-century and three twenty-first-century emission scenarios, were downscaled to three different baseline conditions. Baseline conditions are periods of measured temperature and precipitation data selected to represent twentieth-century climate. The hydrologic effects of the climate projections are evaluated using the Precipitation-Runoff Modeling System (PRMS), which is a watershed hydrology simulation model. The Almanor Catchment in the North Fork of the Feather River basin, California, is used as a case study. Differences and similarities between PRMS simulations of hydrologic components (i.e., snowpack formation and melt, evapotranspiration, and streamflow) are examined, and results indicate that the selection of a specific time period used for baseline conditions has a substantial effect on some, but not all, hydrologic variables. This effect seems to be amplified in hydrologic variables, which accumulate over time, such as soil-moisture content. Results also indicate that uncertainty related to the selection of baseline conditions should be evaluated using a range of different baseline conditions. This is particularly important for studies in basins with highly variable climate, such as the Almanor Catchment.

California↗
USGS2012

Flood pulsing in the Sudd wetland: analysis of seasonal variations in 2 inundation and evapotranspiration in Southern Sudan

Located on the Bahr el Jebel in South Sudan, the Sudd is one of the largest floodplain wetlands in the world. Seasonal inundation drives the hydrologic, geomorphological, and ecological processes, and the annual flood pulse is essential to the functioning of the Sudd. Despite the importance of the flood pulse, various hydrological interventions are planned upstream of the Sudd to increase economic benefits and food security. These will not be without consequences, in particular for wetlands where the biological productivity, biodiversity, and human livelihoods are dependent on the flood pulse and both the costs and benefits need to be carefully evaluated. Many African countries still lack regional baseline information on the temporal extent, distribution, and characteristics of wetlands, making it hard to assess the consequences of development interventions. Because of political instability in Sudan and the inaccessible nature of the Sudd, recent measurements of flooding and seasonal dynamics are inadequate. Analyses of multitemporal and multisensor remote sensing datasets are presented in this paper, in order to investigate and characterize flood pulsing within the Sudd wetland over a 12-month period. Wetland area has been mapped along with dominant components of open water and flooded vegetation at five time periods over a single year. The total area of flooding (both rain and river fed) over the 12 months was 41 334 km2, with 9176 km2 of this constituting the permanent wetland. Mean annual total evaporation is shown to be higher and with narrower distribution of values from areas of open water (1718 mm) than from flooded vegetation (1641 mm). Although the exact figures require validation against ground-based measurements, the results highlight the relative differences in inundation patterns and evaporation across the Sudd.

Earth Interactions↗
USGS2013

The ancient blue oak woodlands of California: longevity and hydroclimatic history

Ancient blue oak trees are still widespread across the foothills of the Coast Ranges, Cascades, and Sierra Nevada in California. The most extensive tracts of intact old-growth blue oak woodland appear to survive on rugged and remote terrain in the south Coast Ranges and on the foothills west and southwest of Mt. Lassen. In our sampling of old-growth stands, most blue oak appear to have recruited to the canopy in the mid- to late-19 th century. The oldest living blue oak tree sampled was over 459-years old and several dead blue oak logs had over 500 annual rings. Precipitation sensitive tree-ring chronologies up to 700-years long have been developed from old blue oak trees and logs. Annual ring-width chronologies of blue oak are strongly correlated with cool season precipitation totals, streamflow in the major rivers of California, and the estuarine water quality of San Francisco Bay. A new network of 36 blue oak chronologies records spatial anomalies in growth that arise from latitudinal changes in the mean storm track and location of landfalling atmospheric rivers. These long, climate-sensitive blue oak chronologies have been used to reconstruct hydroclimatic history in California and will help to better understand and manage water resources. The environmental history embedded in blue oak growth chronologies may help justify efforts to conserve these authentic old-growth native woodlands.

California↗
USGS2013

Variability common to first leaf dates and snowpack in the western conterminous United States

Singular value decomposition is used to identify the common variability in first leaf dates (FLDs) and 1 April snow water equivalent (SWE) for the western United States during the period 1900–2012. Results indicate two modes of joint variability that explain 57% of the variability in FLD and 69% of the variability in SWE. The first mode of joint variability is related to widespread late winter–spring warming or cooling across the entire west. The second mode can be described as a north–south dipole in temperature for FLD, as well as in cool season temperature and precipitation for SWE, that is closely correlated to the El Niño–Southern Oscillation. Additionally, both modes of variability indicate a relation with the Pacific–North American atmospheric pattern. These results indicate that there is a substantial amount of common variance in FLD and SWE that is related to large-scale modes of climate variability.

Earth Interactions↗
USGS2014

Evaluation of statistically downscaled GCM output as input for hydrological and stream temperature simulation in the Apalachicola–Chattahoochee–Flint River Basin (1961–99)

The accuracy of statistically downscaled general circulation model (GCM) simulations of daily surface climate for historical conditions (1961–99) and the implications when they are used to drive hydrologic and stream temperature models were assessed for the Apalachicola–Chattahoochee–Flint River basin (ACFB). The ACFB is a 50 000 km 2 basin located in the southeastern United States. Three GCMs were statistically downscaled, using an asynchronous regional regression model (ARRM), to ⅛° grids of daily precipitation and minimum and maximum air temperature. These ARRM-based climate datasets were used as input to the Precipitation-Runoff Modeling System (PRMS), a deterministic, distributed-parameter, physical-process watershed model used to simulate and evaluate the effects of various combinations of climate and land use on watershed response. The ACFB was divided into 258 hydrologic response units (HRUs) in which the components of flow (groundwater, subsurface, and surface) are computed in response to climate, land surface, and subsurface characteristics of the basin. Daily simulations of flow components from PRMS were used with the climate to simulate in-stream water temperatures using the Stream Network Temperature (SNTemp) model, a mechanistic, one-dimensional heat transport model for branched stream networks. The climate, hydrology, and stream temperature for historical conditions were evaluated by comparing model outputs produced from historical climate forcings developed from gridded station data (GSD) versus those produced from the three statistically downscaled GCMs using the ARRM methodology. The PRMS and SNTemp models were forced with the GSD and the outputs produced were treated as “truth.” This allowed for a spatial comparison by HRU of the GSD-based output with ARRM-based output. Distributional similarities between GSD- and ARRM-based model outputs were compared using the two-sample Kolmogorov–Smirnov (KS) test in combination with descriptive metrics such as the mean and variance and an evaluation of rare and sustained events. In general, precipitation and streamflow quantities were negatively biased in the downscaled GCM outputs, and results indicate that the downscaled GCM simulations consistently underestimate the largest precipitation events relative to the GSD. The KS test results indicate that ARRM-based air temperatures are similar to GSD at the daily time step for the majority of the ACFB, with perhaps subweekly averaging for stream temperature. Depending on GCM and spatial location, ARRM-based precipitation and streamflow requires averaging of up to 30 days to become similar to the GSD-based output. Evaluation of the model skill for historical conditions suggests some guidelines for use of future projections; while it seems correct to place greater confidence in evaluation metrics which perform well historically, this does not necessarily mean those metrics will accurately reflect model outputs for future climatic conditions. Results from this study indicate no “best” overall model, but the breadth of analysis can be used to give the product users an indication of the applicability of the results to address their particular problem. Since results for historical conditions indicate that model outputs can have significant biases associated with them, the range in future projections examined in terms of change relative to historical conditions for each individual GCM may be more appropriate.

Florida, Georgia↗
USGS2015

Recent changes in annual area burned in interior Alaska: The impact of fire management

The Alaskan boreal forest is characterized by frequent extensive wildfires whose spatial extent has been mapped for the past 70 years. Simple predictions based on this record indicate that area burned will increase as a response to climate warming in Alaska. However, two additional factors have affected the area burned in this time record: the Pacific decadal oscillation (PDO) switched from cool and moist to warm and dry in the late 1970s and the Alaska Fire Service instituted a fire suppression policy in the late 1980s. In this paper a geographic information system (GIS) is used in combination with statistical analyses to reevaluate the changes in area burned through time in Alaska considering both the influence of the PDO and fire management. The authors found that the area burned has increased since the PDO switch and that fire management drastically decreased the area burned in highly suppressed zones. However, the temporal analysis of this study shows that the area burned is increasing more rapidly in suppressed zones than in the unsuppressed zone since the late 1980s. These results indicate that fire policies as well as regional climate patterns are important as large-scale controls on fires over time and across the Alaskan boreal forest.

Alaska↗
USGS2016

Urban effects on regional climate: a case study in the Phoenix and Tucson ‘sun’ corridor

Land use and land cover change (LULCC) due to urban expansion alter the surface albedo, heat capacity, and thermal conductivity of the surface. Consequently, the energy balance in urban regions is different from that of natural surfaces. To evaluate the changes in regional climate that could arise due to projected urbanization in the Phoenix-Tucson corridor, Arizona, we applied the coupled WRF-NOAH-UCM (which includes a detailed urban radiation scheme) to this region. Land cover changes were represented using land cover data for 2005 and projections to 2050, and historical North American Regional Reanalysis (NARR) data were used to specify the lateral boundary conditions. Results suggest that temperature changes will be well defined, reflecting the urban heat island (UHI) effect within areas experiencing LULCC. Changes in precipitation are less robust, but seem to indicate reductions in precipitation over the mountainous regions northeast of Phoenix and decreased evening precipitation over the newly-urbanized area.

Earth Interactions↗
USGS2017

Northern hemisphere jet stream positions indices as diagnostic tools for climate and ecosystem dynamics

The latitudinal position of the Northern Hemisphere jet stream (NHJ) modulates the occurrence and frequency of extreme weather events. Precipitation anomalies in particular are associated with NHJ variability; the resulting floods and droughts can have considerable societal and economic impacts. This study develops a new climatology of the 300-hPa NHJ using a bottom-up approach based on seasonally explicit latitudinal NHJ positions. Four seasons with coherent NHJ patterns were identified (January–February, April–May, July–August, and October–November), along with 32 longitudinal sectors where the seasonal NHJ shows strong spatial coherence. These 32 longitudinal sectors were then used as NHJ position indices to examine the influence of seasonal NHJ position on the geographical distribution of NH precipitation and temperature variability and their link to atmospheric circulation pattern. The analyses show that the NHJ indices are related to broad-scale patterns in temperature and precipitation variability, in terrestrial vegetation productivity and spring phenology, and can be used as diagnostic/prognostic tools to link ecosystem and socioeconomic dynamics to upper-level atmospheric patterns.

Earth Interactions↗
USGS2017

Evidence that recent warming is reducing upper Colorado River flows

The upper Colorado River basin (UCRB) is one of the primary sources of water for the western United States, and increasing temperatures likely will elevate the risk of reduced water supply in the basin. Although variability in water-year precipitation explains more of the variability in water-year UCRB streamflow than water-year UCRB temperature, since the late 1980s, increases in temperature in the UCRB have caused a substantial reduction in UCRB runoff efficiency (the ratio of streamflow to precipitation). These reductions in flow because of increasing temperatures are the largest documented temperature-related reductions since record keeping began. Increases in UCRB temperature over the past three decades have resulted in a mean UCRB water-year streamflow departure of −1306 million m 3 (or −7% of mean water-year streamflow). Additionally, warm-season (April through September) temperature has had a larger effect on variability in water-year UCRB streamflow than the cool-season (October through March) temperature. The greater contribution of warm-season temperature, relative to cool-season temperature, to variability of UCRB flow suggests that evaporation or snowmelt, rather than changes from snow to rain during the cool season, has driven recent reductions in UCRB flow. It is expected that as warming continues, the negative effects of temperature on water-year UCRB streamflow will become more evident and problematic.

Arizona, Colorado, New Mexico, Utah, Wyoming↗
USGS2018

Do downscaled general circulation models reliably simulate historical climatic conditions?

The accuracy of statistically downscaled (SD) general circulation model (GCM) simulations of monthly surface climate for historical conditions (1950–2005) was assessed for the conterminous United States (CONUS). The SD monthly precipitation (PPT) and temperature (TAVE) from 95 GCMs from phases 3 and 5 of the Coupled Model Intercomparison Project (CMIP3 and CMIP5) were used as inputs to a monthly water balance model (MWBM). Distributions of MWBM input (PPT and TAVE) and output [runoff (RUN)] variables derived from gridded station data (GSD) and historical SD climate were compared using the Kolmogorov–Smirnov (KS) test For all three variables considered, the KS test results showed that variables simulated using CMIP5 generally are more reliable than those derived from CMIP3, likely due to improvements in PPT simulations. At most locations across the CONUS, the largest differences between GSD and SD PPT and RUN occurred in the lowest part of the distributions (i.e., low-flow RUN and low-magnitude PPT). Results indicate that for the majority of the CONUS, there are downscaled GCMs that can reliably simulate historical climatic conditions. But, in some geographic locations, none of the SD GCMs replicated historical conditions for two of the three variables (PPT and RUN) based on the KS test, with a significance level of 0.05. In these locations, improved GCM simulations of PPT are needed to more reliably estimate components of the hydrologic cycle. Simple metrics and statistical tests, such as those described here, can provide an initial set of criteria to help simplify GCM selection.

Earth Interactions↗
USGS2020

Basinwide hydroclimatic drought in the Colorado River basin

The Colorado River basin (CRB) supplies water to approximately 40 million people and is essential to hydropower generation, agriculture, and industry. In this study, a monthly water balance model is used to compute hydroclimatic water balance components (i.e., potential evapotranspiration, actual evapotranspiration, and runoff) for the period 1901–2014 across the entire CRB. The time series of monthly runoff is aggregated to compute water-year runoff and then used to identify drought periods in the basin. For the 1901–2014 period, eight basinwide drought periods were identified. The driest drought period spanned years 1901–04, whereas the longest drought period occurred during 1943–56. The eight droughts were primarily driven by winter precipitation deficits rather than warm temperature anomalies. In addition, an analysis of prehistoric drought for the CRB—computed using tree-ring-based reconstructions of the Palmer drought severity index—indicates that during some past centuries drought frequency was higher than during the twentieth century and that some centuries experienced droughts that were much longer than those during the twentieth century. More frequent or longer droughts than those that occurred during the twentieth century, combined with continued warming associated with climate change, may lead to substantial future water deficits in the CRB.

Arizona, California, Colorado, Nevada, New Mexico,↗
USGS2021

Simulated atmospheric response to four projected land-use land-cover change scenarios for 2050 in the north-central United States

Land-use land-cover change (LULCC) has become an important topic of research for the central United States because of the extensive conversion of the natural prairie into agricultural land, especially in the northern Great Plains. As a result, shifts in the natural climate (minimum/maximum temperature, precipitation, etc.) across the north-central United States have been observed, as noted within the Fourth National Climate Assessment (NCA4) report. Thus, it is necessary to understand how further LULCC will affect the near-surface atmosphere, the lower troposphere, and the planetary boundary layer (PBL) atmosphere over this region. The goal of this work was to investigate the utility of a new future land-use land-cover (LULC) dataset within the Weather Research and Forecasting (WRF) modeling system. The present study utilizes a modeled future land-use dataset developed by the Forecasting Scenarios of Land-Use Change (FORE-SCE) model to investigate the influence of future (2050) land use on a simulated PBL development within the WRF Model. Three primary areas of LULCC were identified within the FORE-SCE future LULC dataset across Nebraska and South Dakota. Variations in LULC between the 2005 LULC control simulation and four FORE-SCE simulations affected near-surface temperature (0.5°–1°C) and specific humidity (0.3–0.5 g kg −1 ). The differences noted in the temperature and moisture fields affected the development of the simulated PBL, leading to variations in PBL height and convective available potential energy. Overall, utilizing the FORE-SCE dataset within WRF produced notable differences relative to the control simulation over areas of LULCC represented in the FORE-SCE dataset.

Nebraska, South Dakota↗
USGS2022

Temporal greenness trends in stable natural land cover and relationships with climatic variability across the conterminous United States

Assessment of temporal trends in vegetation greenness and related influences aids understanding of recent change in terrestrial ecosystems and feedbacks from weather, climate, and environment. We analyzed 1-km normalized difference vegetation index (NDVI) timeseries data (1989–2016) derived from the Advanced Very High Resolution Radiometer (AVHRR) and developed growing season time-integrated NDVI (GS-TIN) for estimating seasonal vegetation activity across stable natural land cover in the conterminous United States (CONUS). After removing areas from analysis that had experienced land cover conversion or modification, we conducted a monotonic trend analysis on the GS-TIN timeseries and found that significant positive temporal trends occurred over 35% of the area, while significant negative trends were observed over only 3.5%. Positive trends were prevalent in the forested lands of the eastern third of CONUS and far northwest, as well as in grasslands in the north central plains. We observed negative and nonsignificant trends mainly in the shrublands and grasslands across the northwest, southwest, and west central plains. To understand the relationship of climate variability with these temporal trends, we conducted partial and multiple correlation analyses on GS-TIN, growing season temperature, and water-year precipitation timeseries. The GS-TIN trends in northern forests were positively correlated with temperature. The GS-TIN trends in the central and western shrublands and grasslands were negatively correlated with temperature and positively correlated with precipitation. Our results revealed spatial patterns in vegetation greenness trends for different stable natural vegetation types across CONUS, enhancing understanding gained from prior studies based on coarser 8-km AVHRR data.

Earth Interactions↗
USGS2022

Extensive droughts in the conterminous United States during multiple centuries

Extensive and severe droughts have substantial effects on water supplies, agriculture, and aquatic ecosystems. To better understand these droughts, we used tree-ring-based reconstructions of the Palmer drought severity index (PDSI) for the period 1475–2017 to examine droughts that covered at least 33% of the conterminous United States (CONUS). We identified 37 spatially extensive drought events for the CONUS and examined their spatial and temporal patterns. The duration of the extensive drought events ranged from 3 to 12 yr and on average affected 43% of the CONUS. The recent (2000–08) drought in the southwestern CONUS, often referred to as the turn-of-the-century drought, is likely one of the longest droughts in the CONUS during the past 500 years. A principal components analysis of the PDSI data from 1475 through 2017 resulted in three principal components (PCs) that explain about 48% of the variability of PDSI and are helpful to understand the temporal and spatial variability of the 37 extensive droughts in the CONUS. Analyses of the relations between the three PCs and well-known climate indices, such as indices of El Niño–Southern Oscillation, indicate statistically significant correlations; however, the correlations do not appear to be large enough (all with an absolute value less than 0.45) to be useful for the development of drought prediction models.

Conterminous United States↗
USGS2024

Investigating the atmospheric conditions associated with impactful shallow landslides in California (USA)

Shallow landslides are often triggered during rainfall events, which can increase subsurface soil water pressure and destabilize hillslopes. The likelihood of regional shallow landslide initiation is often assessed through a comparison of rainfall intensity and duration to pre-established thresholds. While informative for landslide warning, this exclusive focus on rainfall exceeding thresholds does not consider the meteorological conditions producing the rainfall. Here, we ask the question, are there common meteorological characteristics that lead to landslide-triggering precipitation? We develop a catalog of 18 post-1995 widespread, impactful shallow landslide events occurring within 13 storms across California, USA, where initiation time could be constrained to a ≤6-h window. We examine storm characteristics during the landslide initiation window using atmospheric reanalysis products, radar observations, and quantitative precipitation estimates. We find that, while there are some common atmospheric characteristics across landslide events, they can occur under a range of atmospheric conditions. For example, all Northern California landslide events assessed are associated with moderate to strong atmospheric rivers (ARs), while Southern California landslides feature non-AR to strong AR conditions. The storm events evaluated herein share many characteristics of hydrologically important storms in California that did not necessarily result in landslides; thus, atmospheric characteristics alone may not be sufficient to determine whether landslides will occur. However, documenting the characteristics of landslide-triggering storms defines the conditions under which landslides tend to occur, provides analog events that can be useful in forecast applications, helps define future research directions relating to atmospheric conditions and landslides, and supports interdisciplinary research efforts.

California↗
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↗