skip to main content


Title: Atmospheric River Contributions to Ice Sheet Hydroclimate at the Last Glacial Maximum
Abstract

Atmospheric rivers (ARs) are an important driver of surface mass balance over today's Greenland and Antarctic ice sheets. Using paleoclimate simulations with the Community Earth System Model, we find ARs also had a key influence on the extensive ice sheets of the Last Glacial Maximum (LGM). ARs provide up to 53% of total precipitation along the margins of the eastern Laurentide ice sheet and up to 22%–27% of precipitation along the margins of the Patagonian, western Cordilleran, and western Fennoscandian ice sheets. Despite overall cold conditions at the LGM, surface temperatures during AR events are often above freezing, resulting in more rain than snow along ice sheet margins and conditions that promote surface melt. The results suggest ARs may have had an important role in ice sheet growth and melt during previous glacial periods and may have accelerated ice sheet retreat following the LGM.

 
more » « less
Award ID(s):
1903600 1903528
PAR ID:
10395338
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
50
Issue:
1
ISSN:
0094-8276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. This dataset contains the atmospheric river catalogues and the associated precipitation and temperature data for the Preindustrial and Last Glacial Maximum CESM2 simulations presented in the GRL manuscript:  Atmospheric river contributions to ice sheet hydro climate at the Last Glacial Maximum. The atmospheric river catalogue files (zipped) are in netcdf format and organized by year. There are 100 years of data for both simulations.  The Preindustrial simulation catalogue begins in model year 41 and ends in model year 140.  The LGM simulation catalogue begins in model year 1 and ends in year 100. Each yearly file has a temporal resolution of 6 hours (1460 time steps each file) and a spatial resolution of 0.9° x 1.25° (the native resolution of the CESM simulation). A variable in the file called "ar_binary_tag" indicates whether an atmospheric river is present at each grid cell and each tilmestep: 1 indicates an atmospheric river is present; 0 indicates an atmospheric river is not present.  The precipitation and temperature files are 100-year annual or 100-year seasonal averages of atmospheric river precipitation/temperature. See the Methods section of the article for more details on the atmospheric river detection algorithm and precipitation/temperature calculations.

    Associated article abstract:

    Atmospheric rivers (ARs) are an important driver of surface mass balance over today’s Greenland and Antarctic ice sheets. Using paleoclimate simulations with the Community Earth System Model, we find ARs also had a key influence on the extensive ice sheets of the Last Glacial Maximum (LGM). ARs provide up to 53% of total precipitation along the margins of the eastern Laurentide ice sheet and up to 22-27% of precipitation along the margins of the Patagonian, western Cordilleran, and western Fennoscandian ice sheets. Despite overall cold conditions at the LGM, surface temperatures during AR events are often above freezing, resulting in more rain than snow along ice sheet margins and conditions that promote surface melt. The results suggest  ARs may have had an important role in ice sheet growth and melt during previous glacial periods and may have accelerated ice sheet retreat following the LGM.

     
    more » « less
  2. Abstract. Studying the retreat of the Patagonian Ice Sheet (PIS) during the last deglaciation represents an important opportunity to understand how ice sheets outside the polar regions have responded to deglacial changes in temperature and large-scale atmospheric circulation. At the northernmost extension of the PIS during the Last Glacial Maximum (LGM), the Chilean Lake District (CLD) was influenced by the southern westerly winds (SWW), which strongly modulated the hydrologic and heat budgets of the region. Despite progress in constraining the nature and timing of deglacial ice retreat across this area, considerable uncertainty in the glacial history still exists due to a lack of geologic constraints on past ice margin change. Where the glacial chronology is lacking, ice sheet models can provide important insight into our understanding of the characteristics and drivers of deglacial ice retreat. Here we apply the Ice Sheet and Sea-level System Model (ISSM) to simulate the LGM and last deglacial ice history of the PIS across the CLD at high spatial resolution (450 m). We present a transient simulation of ice margin change across the last deglaciation using climate inputs from the National Center for Atmospheric Research Community Climate System Model (CCSM3) Trace-21ka experiment. At the LGM, the simulated ice extent across the CLD agrees well with the most comprehensive reconstruction of PIS ice history (PATICE). Coincident with deglacial warming, ice retreat ensues after 19 ka, with large-scale ice retreat occurring across the CLD between 18 and 16.5 ka. By 17 ka, the northern portion of the CLD becomes ice free, and by 15 ka, ice only persists at high elevations as mountain glaciers and small ice caps. Our simulated ice history agrees well with PATICE for early deglacial ice retreat but diverges at and after 15 ka, where the geologic reconstruction suggests the persistence of an ice cap across the southern CLD until 10 ka. However, given the high uncertainty in the geologic reconstruction of the PIS across the CLD during the later deglaciation, this work emphasizes a need for improved geologic constraints on past ice margin change. While deglacial warming drove the ice retreat across this region, sensitivity tests reveal that modest variations in wintertime precipitation (∼10 %) can modulate the pacing of ice retreat by up to 2 ka, which has implications when comparing simulated outputs of ice margin change to geologic reconstructions. While we find that TraCE-21ka simulates large-scale changes in the SWW across the CLD that are consistent with regional paleoclimate reconstructions, the magnitude of the simulated precipitation changes is smaller than what is found in proxy records. From our sensitivity analysis, we can deduce that larger anomalies in precipitation, as found in paleoclimate proxies, may have had a large impact on modulating the magnitude and timing of deglacial ice retreat. This fact highlights an additional need for better constraints on the deglacial change in strength, position, and extent of the SWW as it relates to understanding the drivers of deglacial PIS behavior.

     
    more » « less
  3. Abstract

    Atmospheric rivers (ARs) are efficient mechanisms for transporting atmospheric moisture from low latitudes to the Antarctic Ice Sheet (AIS). While AR events occur infrequently, they can lead to extreme precipitation and surface melt events on the AIS. Here we estimate the contribution of ARs to total Antarctic precipitation, by combining precipitation from atmospheric reanalyses and a polar‐specific AR detection algorithm. We show that ARs contribute substantially to Antarctic precipitation, especially in East Antarctica at elevations below 3,000 m. ARs contribute substantially to year‐to‐year variability in Antarctic precipitation. Our results highlight that ARs are an important component for understanding present and future Antarctic mass balance trends and variability.

     
    more » « less
  4. Abstract

    Retreat or advance of an ice sheet perturbs the Earth's solid surface, rotational vector, and the gravitational field, which in turn feeds back onto the evolution of the ice sheet over a range of timescales. Throughout the last glacial cycle, ice sheets over the Northern Hemisphere have gone through multiple growth and retreat phases, but the dynamics during these phases are not well understood. In this study, we apply a coupled ice sheet‐glacial isostatic adjustment model to simulate the Northern Hemisphere Ice Sheets over the last glacial cycle. We focus on understanding the influence of solid Earth deformation and gravitational field perturbations associated with surface (ice and water) loading changes on the dynamics of terrestrial and marine‐based ice sheets during different phases of the glacial cycle. Our results show that solid Earth deformation enhances glaciation during growth phases and melting during retreat phases in terrestrial regions through ice‐elevation feedback, and gravitational field perturbations have a stabilizing influence on marine‐based ice sheets in regions such as Hudson Bay in North America and Barents and Kara Seas in Eurasia during retreat phases through sea‐level feedback. Our results also indicate that solid Earth deformation influences the relative sensitivity of the North American and Eurasian ice sheets to climate and thus the timing and magnitude of their fluctuations throughout the last glacial cycle.

     
    more » « less
  5. Abstract

    Quantifying glacial erosion contributes to our understanding of landscape evolution and topographic relief production in high altitude and high latitude areas. Combiningin situ10Be and26Al analysis of bedrock, boulder, and river sand samples, geomorphological mapping, and field investigations, we examine glacial erosion patterns of former ice caps in the Shaluli Shan of the southeastern Tibetan Plateau. The general landform pattern shows a zonal pattern of landscape modification produced by ice caps of up to 4000 km2during pre‐LGM (Last Glacial Maximum) glaciations, while the dating results and landforms on the plateau surface imply that the LGM ice cap further modified the scoured terrain into different zones. Modeled glacial erosion depth of 0–0.38 m per 100 ka bedrock sample located close to the western margin of the LGM ice cap, indicates limited erosion prior to LGM and Late Glacial moraine deposition. A strong erosion zone exists proximal to the LGM ice cap marginal zone, indicated by modeled glacial erosion depth >2.23 m per 100 ka from bedrock samples. Modeled glacial erosion depths of 0–1.77 m per 100 ka from samples collected along the edge of a central upland, confirm the presence of a zone of intermediate erosion in‐between the central upland and the strong erosion zone. Significant nuclide inheritance in river sand samples from basins on the scoured plateau surface also indicate restricted glacial erosion during the last glaciation. Our study, for the first time, shows clear evidence for preservation of glacial landforms formed during previous glaciations under non‐erosive ice on the Tibetan Plateau. As patterns of glacial erosion intensity are largely driven by the basal thermal regime, our results confirm earlier inferences from geomorphology for a concentric basal thermal pattern for the Haizishan ice cap during the LGM. © 2018 John Wiley & Sons, Ltd.

     
    more » « less