Code and glacial isostatic adjustment model outputs associated with "Glacial isostatic adjustment driven by asymmetric ice sheet melt during the last interglacial causes multiple local sea level peaks"

### Access
NetCDF files be accessed and downloaded from the directory via: [http://arcticdata.io/data/10.18739/A2S756N0M](http://arcticdata.io/data/10.18739/A2S756N0M).

### Overview
This dataset comprises the code used to produce the results in Creel, R.C., Austermann, J.A., Glacial isostatic adjustment driven by asymmetric ice sheet melt during the Last Interglacial causes multiple local sea-level peaks. Geology. https://doi.org/10.1130/G52483.1.

Research abstract:

Global mean sea-level (GMSL) change during the Last Interglacial (LIG, 129−116 kiloannum (ka)) gives perspective on how ice sheets respond to warming. Observations of multiple peaks in LIG relative sea level (RSL) records, combined with an assumption that the Laurentide Ice Sheet (LIS) collapsed prior to the LIG, have been used to infer Greenland and Antarctic ice sheet melt histories as well as oscillations in LIG GMSL. However, evidence for an LIS outburst flood at ca. 125 ka and extensive early-LIG Antarctic melt suggests that Laurentide remnants may have persisted longer into the LIG than typically thought even as Antarctic melt accelerated. Here, we explore the effect of concurrent early-Holocene Laurentide persistence and Antarctic collapse on glacial isostatic adjustment and sea level. In our models, we hold GMSL constant at present levels (i.e., GMSL = 0) from 128 ka to 117 ka by balancing excess Laurentide ice with early-LIG Antarctic melt. We find that due to glacial isostatic adjustment, this synchronous but asymmetric ice change causes multiple RSL peaks, separated by ∼4.2 ± 2.5 m of RSL fall near North America and ∼1.3 ± 0.7 m around the Indian Ocean. This spatial pattern resembles observations. These results show that multiple peaks in LIG RSL could have occurred with asymmetric ice changes between the Northern and Southern Hemisphere that sum to little, if any, change in GMSL. Our work highlights the need for LIG modeling studies to consider that dynamic cryospheric changes can occur even with near-constant GMSL.

This research was conducted at Lamont Doherty Earth Observatory in 2022 and 2023. It is entirely composed of modeling and compilation of existing data; no new data were produced for the study. The methodologies employed include glacial isostatic adjustment modeling using MATLAB code written by Jacky Austermann and post-processing of the resulting models using python scripts written by Roger Creel.