An official website of the United States government
Abstract. Geologic archives of the Laurentide Ice Sheet (LIS) provide abundant constraints regarding the size and extent of the ice sheet during the Last Glacial Maximum (LGM) and throughout the deglaciation. Direct observations of LGM LIS thickness are non-existent, however, due to ice surface elevations likely exceeding those of even the tallest summits in the northeastern United States (NE USA). Geomorphic and isotopic data from mountains across the NE USA can inform basal conditions, including the presence of warm- or cold-based regimes, while covered by ice. While warm-based ice and erosive conditions likely existed on the flanks of these summits and throughout neighboring valleys, cosmogenic nuclide inheritance and frost-riven blockfields on summits suggest ineffective glacial erosion and cold-based ice conditions. Geologic reconstructions indicate that a complex erosional and thermal regime likely existed across the NE USA sometime during and after the LGM, although this has not been confirmed by ice sheet models. Instead, current ice sheet models simulate warm-based ice conditions across this region, with disagreement likely arising from the use of low-resolution meshes (e.g., > 20 km) which are unable to resolve the high bedrock relief across the NE USA that strongly influenced overall ice flow and the complex LIS thermal state. Here, we use a newer-generation ice sheet model, the Ice-sheet and Sea-level System Model (ISSM), to simulate the LGM conditions of the LIS across the NE USA and in three localities with high bedrock relief (Adirondack Mountains, White Mountains, and Mount Katahdin), with results confirming the existence of a complex thermal regime as interpreted from the geologic data. The model uses a small ensemble of LGM climate boundary conditions and a high-resolution horizontal mesh that resolves bedrock features down to 30 m to reconstruct LGM ice flow, ice thickness, and thermal conditions. These results indicate that, across the NE USA, polythermal conditions existed during the LGM. While the majority of this domain is simulated to be warm-based, cold-based ice persists where ice velocities are slow (< 15 m yr−1), particularly across regional ice divides (e.g., Adirondack Mountains). Additionally, sharp thermal boundaries are simulated where cold-based ice across high-elevation summits (White Mountains and Mount Katahdin) flanks warm-based ice in adjacent valleys. We find that the elevation of this simulated thermal boundary ranges between 800–1500 m, largely supporting geologic interpretations that polythermal ice conditions existed across the NE USA during the LGM; however, this boundary varies geographically. In general, we show that a model using a finer spatial resolution compared to older models can simulate the polythermal conditions captured in the geologic data, with the model output being of potential utility for site selection in future geologic studies and for geomorphic interpretations of landscape evolution.
more »
« less
Ice sheet model simulations reveal polythermal ice conditions existed across the NE USA during the Last Glacial Maximum
Abstract. Geologic evidence of the Laurentide Ice Sheet (LIS) provides abundant constraints on the areal extent of the ice sheet during the Last Glacial Maximum (LGM). Direct observations of LGM LIS thickness are non-existent, however, with most geologic data across high elevation summits in the Northeastern United States (NE USA) often showing signs of inheritance, indicative of weakly erosive ice flow and the presence of cold-based ice. While warm-based ice and erosive conditions likely existed on the flanks of these summits and throughout neighboring valleys, summit inheritance issues have hampered efforts to constrain the timing of the emergence of ice-free conditions at high elevation summits. These geomorphic reconstructions indicate that a complex erosional and thermal regime likely existed across the southeasternmost extent of the LIS sometime during the LGM, although this has not been confirmed by ice sheet models. Instead, current ice sheet models simulate warm-based ice conditions across this region, with disagreement likely arising from the use of low resolution meshes (e.g., >20 km) which are unable to resolve the high bedrock relief across this region that strongly influenced overall ice flow and the complex LIS thermal state. Here we use a newer generation ice sheet model, the Ice-sheet and Sea-level System Model (ISSM), to simulate the LGM conditions of the LIS across the NE USA and at 3 localities with high bedrock relief (Adirondack Mountains, White Mountains, and Mount Katahdin), with results confirming the existence of a complex thermal regime as interpreted by the geologic data. The model uses higher-order physics, a small ensemble of LGM climate boundary conditions, and a high-resolution horizontal mesh that resolves bedrock features down to 30 meters to reconstruct LGM ice flow, ice thickness, and thermal conditions. These results indicate that across the NE USA, polythermal conditions existed during the LGM. While the majority of this domain is simulated to be warm-based, cold-based ice persists where ice velocities are slow (<15 m/yr) particularly across regional ice divides (e.g., Adirondacks). Additionally, sharp thermal boundaries are simulated where cold-based ice across high elevation summits (White Mountains and Mount Katahdin) flank warm-based ice in adjacent valleys. Because geologic data is geographically limited, these high-resolution simulations can help fill gaps in our understanding of the geographical distribution of the polythermal ice during the LGM. We find that the elevation of this simulated thermal boundary ranges between 800–1500 meters, largely supporting geologic interpretations that polythermal ice conditions existed across NE USA during the LGM, however this boundary varies geographically. In general, we show that a model with finer spatial resolution and higher order physics is able to simulate the polythermal conditions captured in the geologic data, with model output being of potential utility for site selection in future geologic studies and geomorphic interpretation of landscape evolution.
more »
« less
- Award ID(s):
- 2133699
- PAR ID:
- 10524718
- Publisher / Repository:
- The Cryosphere
- Date Published:
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Accurate reconstruction of Laurentide Ice Sheet volume changes following the Last Glacial Maximum is critical for understanding ice sheet contribution to sea-level rise, the resulting influence of meltwater on oceanic circulation, and the spatial and temporal patterns of deglaciation. Here, we provide empirical constraints on Laurentide Ice Sheet thinning during the last deglaciation by measuring in situ cosmogenic 10Be in 81 samples collected along vertical transects of nine mountains in the northeastern United States. In conjunction with 107 exposure age samples over five vertical transects from previous studies, we reconstruct ice sheet thinning history. At peripheral sites (within 200 km of the terminal moraine), we find evidence for ∼600 m of thinning between 19.5 ka and 17.5 ka, which is coincident with the slow initial margin retreat indicated by varve records. At locations >400 km north of the terminal moraine, exposure ages above and below 1200 m a.s.l. exhibit different patterns. Ages above this elevation are variable and older, while lower elevation ages are indistinguishable over 800−1000 m elevation ranges, a pattern that suggests a subglacial thermal boundary at ∼1200 m a.s.l. separating erosive, warm-based ice below and polythermal, minimally erosive ice above. Low-elevation ages from up-ice mountains are between 15 ka and 13 ka, which suggests rapid thinning of ∼1000 m coincident with Bølling-Allerød warming. These rates of rapid paleo-ice thinning are comparable to those of other vertical exposure age transects around the world and may have been faster than modern basin-wide thinning rates in Antarctica and Greenland, which suggests that the southeastern Laurentide Ice Sheet was highly sensitive to a warming climate.more » « less
-
We apply geologic evidence from ice-free areas in Antarctica to evaluate model simulations of ice sheet response to warm climates. This is important because such simulations are used to predict ice sheet behaviour in future warm climates, but geologic evidence of smaller-than-present past ice sheets is buried under the present ice sheet and therefore generally unavailable for model benchmarking. We leverage an alternative accessible geologic dataset for this purpose: cosmogenic-nuclide concentrations in bedrock surfaces of interior nunataks. These data produce a frequency distribution of ice thickness over multimillion-year periods, which is also simulated by ice sheet modelling. End-member transient models, parameterized with strong and weak marine ice sheet instability processes and ocean temperature forcings, simulate large and small sea-level impacts during warm periods and also predict contrasting and distinct frequency distributions of ice thickness. We identify regions of Antarctica where predicted frequency distributions reveal differences in end-member ice sheet behaviour. We then demonstrate that a single comprehensive dataset from one bedrock site in West Antarctica is sufficiently detailed to show that the data are consistent only with a weak marine ice sheet instability end-member, but other less extensive datasets are insufficient and/or ambiguous. Finally, we highlight locations where collecting additional data could constrain the amplitude of past and therefore future response to warm climates.more » « less
-
Abstract In the Ross Sea sector of Antarctica, periodic large-scale marine ice-sheet fluctuations since the mid-Miocene are recorded by drill core and seismic data, revealing a dynamic ice-sheet response to past increases in temperature and atmospheric CO2. In the adjacent, predominantly ice-free McMurdo Dry Valleys (MDVs), preserved terrestrial landscapes reflect persistent cold conditions and have been interpreted as indicators of a stable polar ice sheet, implying that the Antarctic Ice Sheet was largely insensitive during past warm periods. These disparate data-based perspectives highlight a long-standing debate around the past stability of the Antarctic Ice Sheet, with direct implications for the future ice-sheet response to ongoing climate warming. We reconcile marine records of dynamic ice-sheet behavior and episodic open-marine conditions with nearby ancient terrestrial landscapes recording consistent cold-polar conditions. Coupled ice-sheet and regional climate models nested at a high resolution are used to investigate surface temperatures in the MDVs during past warm periods. We find that high-elevation regions of the MDVs remain below freezing even when ice-free conditions prevail in the nearby Ross Sea. We compare observed landscapes with the spatial extent of modeled persistent cold conditions required for preservation of these ancient features, demonstrating that frozen MDVs landscapes could have coexisted with receded or collapsed ice sheets during past warm periods.more » « less
-
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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Posted Content:
- https://doi.org/10.5194/egusphere-2024-2091
