Dynamic treeline and cryosphere response to pronounced mid-Holocene climatic variability in the US Rocky Mountains

Gregory T. Pederson; Daniel K. Stahle; David B McWethy; Matthew Toohey; Johann Jungclaus; Craig Lee; Justin T. Martin; Mio Alt; Nickolas E. Kichas; Nathan J. Chellman; Joseph R. McConnell; Cathy Whitlock

doi:10.1073/pnas.2412162121

USGS · 70261894

Dynamic treeline and cryosphere response to pronounced mid-Holocene climatic variability in the US Rocky Mountains

Gregory T. Pederson·Daniel K. Stahle·David B McWethy·Matthew Toohey·Johann Jungclaus·Craig Lee·Justin T. Martin·Mio Alt·Nickolas E. Kichas·Nathan J. Chellman·Joseph R. McConnell·Cathy Whitlock

View original source ↗Export metadata ↓

Abstract

Climate-driven changes in high-elevation forest distribution and reductions in snow and ice cover have major implications for ecosystems and global water security. In the Greater Yellowstone Ecosystem of the Rocky Mountains (United States), recent melting of a high-elevation (3,091 m asl) ice patch exposed a mature stand of whitebark pine ( Pinus albicaulis ) trees, located ~180 m in elevation above modern treeline, that date to the mid-Holocene (c. 5,950 to 5,440 cal y BP). Here, we used this subfossil wood record to develop tree-ring-based temperature estimates for the upper-elevation climate conditions that resulted in ancient forest establishment and growth and the subsequent regional ice-patch growth and downslope shift of treeline. Results suggest that mid-Holocene forest establishment and growth occurred under warm-season (May-Oct) mean temperatures of 6.2 °C (±0.2 °C), until a multicentury cooling anomaly suppressed temperatures below 5.8 °C, resulting in stand mortality by c. 5,440 y BP. Transient climate model simulations indicate that regional cooling was driven by changes in summer insolation and Northern Hemisphere volcanism. The initial cooling event was followed centuries later (c. 5,100 y BP) by sustained Icelandic volcanic eruptions that forced a centennial-scale 1.0 °C summer cooling anomaly and led to rapid ice-patch growth and preservation of the trees. With recent warming (c. 2000–2020 CE), warm-season temperatures now equal and will soon exceed those of the mid-Holocene period of high treeline. It is likely that perennial ice cover will again disappear from the region, and treeline may expand upslope so long as plant-available moisture and disturbance are not limiting.

Explore related subjects

Proceedings of the National Academy of Sciences Article

Source-reported bounding extent: 44.77582421595028° to 45.200185933759144° latitude; -110.40405049321484° to -109.75503610328678° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Gregory T. Pederson, Daniel K. Stahle, David B McWethy, Matthew Toohey, Johann Jungclaus, Craig Lee, Justin T. Martin, Mio Alt, Nickolas E. Kichas, Nathan J. Chellman, Joseph R. McConnell, Cathy Whitlock. 2024-12-30. Dynamic treeline and cryosphere response to pronounced mid-Holocene climatic variability in the US Rocky Mountains. https://doi.org/10.1073/pnas.2412162121

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

Deep critical zone controls on shallow landslides

The deep critical zone (CZ) has long been recognized for its importance in influencing shallow landslides but was not considered feasible to include in slope stability models at the watershed scale. In this study, we demonstrate that simple approximations of the CZ in a fully coupled hydrologic and soil slope stability model can effectively capture the location, timing, and likely size of shallow landslides. To achieve this, we use coupled, process-based models that incorporate the effects of 1) deep CZ structures, 2) three-dimensional transient hydrology, and 3) multidimensional slope stability, calibrated with data from an intensively monitored field site. Our results show that the hydrologically active deep CZ guides groundwater flow, influencing where it drains from or exfiltrates to the soil mantle, producing distinct patterns of soil saturation and seepage forces at the soil-bedrock boundary. Deep conductive weathered critical zone drains the soil mantle, reducing the likelihood of destabilizing pore pressures, while the downslope thinning of the CZ forces groundwater to the surface. This creates localized instability and a tendency for similar-sized landslides across the landscape. In contrast, the absence of conductive weathered bedrock results in more widespread destabilizing pore pressures, leading to larger landslides and the likelihood of landslides earlier in a storm than in landscapes underlain by a deep CZ. Our findings suggest that first-order variations of deep CZ can provide physical explanations for variations observed in the susceptibility, magnitude, and timing of shallow landslides, and that CZ structure may be inferred from patterns and timing of landsliding.

Proceedings of the National Academy of Sciences↗

USGS2026

Widespread terrestrial ecosystem disruption at the onset of the Paleocene–Eocene Thermal Maximum

The Paleocene–Eocene Thermal Maximum (PETM, ~56 Mya) interval was marked by massive 13 C-depleted carbon emissions into the ocean/atmosphere system, manifested as a negative carbon isotope excursion (CIE) in sedimentary components, and ~5 °C global average warming. Episodes of hydrological perturbations and soil-erosion have been widely documented for the PETM but their link with vegetation- and carbon cycle changes remain poorly constrained. Here, we present organic microfossil evidence showing a strong increase in fern-dominated pioneer vegetation that replaced coniferous forests on the margin of the Norwegian Sea during the first millennia of the CIE. With the present stratigraphic constraints, the “fern spike” occurred simultaneously in terrestrial settings along the North Sea, Arctic Ocean, the US east coast and in southern Australia, indicating that pioneer vegetation persisted for several millennia following a partial collapse of previously stable terrestrial ecosystems. Both the ferns and influx of microcharcoal imply recurrent physical disturbance, including soil destabilization and erosion, potentially linked to droughts, wildfires, and strong hydrological forcing resulting from extreme climate change. Together with evidence for reworked clay minerals and ancient organic matter (kerogen), these findings show that highly disturbed terrestrial ecosystems were widespread across mid- and high-latitude regions globally. Carbon cycle model simulations suggest that a substantial loss of standing and buried biomass, along with oxidation of soil organic matter, acted as important positive feedbacks during the onset of the CIE. Additionally, enhanced kerogen weathering likely contributed as another major positive feedback throughout both the onset and main phase of the CIE.

Proceedings of the National Academy of Sciences↗

USGS2025

Landscape changes elevate the risk of avian influenza virus diversification and emergence in the East Asian–Australasian Flyway

Highly pathogenic avian influenza viruses (HPAIV) persistently threaten wild waterfowl, domestic poultry, and public health. The East Asian–Australasian Flyway plays a crucial role in HPAIV dynamics due to its large populations of migratory waterfowl and poultry. Over recent decades, this flyway has undergone substantial landscape changes, including both losses and gains of waterfowl habitats. These changes can affect waterfowl distributions, increase contact with poultry, and consequently alter ecological conditions that favor avian influenza virus (AIV) evolution. However, limited research has assessed these likely impacts. Here, we integrated empirical data and an individual-based model to simulate AIV transmission in migratory waterfowl and domestic poultry, including wild-to-poultry spillover and reassortment dynamics in poultry, across landscapes representing the years 2000 and 2015. We used the reassortment incidence as a proxy for ecological and transmission conditions that support viral diversification and the emergence of novel subtypes. Our simulations show that landscape change reshaped the waterfowl distribution, facilitated bird aggregation at improved habitats, increased coinfection, and raised reassortment rate by 1,593%, indicating a substantially higher potential for viral diversification and emergence. Model-generated risk maps show expanded and increased reassortment risk in southeastern China, the Yellow River Basin, and northeastern China. These findings suggest the importance of landscape change as a driver of potential AIV diversification and subtype emergence. This underscores the need for interdisciplinary approaches that integrate landscape dynamics, host movement, and viral evolution to better assess and mitigate future risk.

East Asian–Australasian Flyway↗