The transient simulation was performed with the Max Planck Institute for Meteorology Earth System Model version 1.2 in coarse resolution (MPI-ESM-CR) coupled to the ice-sheet model mPISM and the solid-earth model VILMA. MPI-ESM includes the spectral atmospheric model ECHAM6.3 at T31 horizontal resolution (approx. 3.75°) and 31 vertical levels, the land surface vegetation model JSBACH3.2, and the primitive equation ocean model MPIOM1.6 with a nominal resolution of 3°. Embedded into MPIOM1.6 is an Eulerian iceberg model (Erokhina and Mikolajewicz 2024). mPISM is based on PISM0.7.3. and was employed using a polar stereographic grid with a resolution of 10 km in the northern hemisphere and 15 km for Antarctica. VILMA was used in its 1D configuration and discretized on a Gaussian F128 grid with a nominal resolution of 0.7°. The applied setup was introduced in detail in Mikolajewicz et al. (2024).

Following an asynchronous spin-up from 46 ka to 26ka during which MPI-ESM was run with an acceleration factor of 10, the model was integrated without acceleration from a glacial state at 26 ka to the year 1950 with prescribed atmospheric greenhouse gas concentrations (Köhler et al., 2017) and insolation (Berger & Loutre, 1991). This simulation used annualy varying volcanic forcing (Schindlbeck-Belo et al., 2024). Forcing fields between MPI-ESM, mPISM, and VILMA as well as changes in the topography, glacier mask, river pathways, ocean bathymetry, and land-sea mask were updated every 10 years. The experiment was performed as part of a model ensemble that contains simulations differing in terms of their tuning parameters, parameterizations, bug fixes and volcanic boundary forcing. These differences are described according to the CMIP6 nomenclature, where r denotes the realization, i the initialization method, p differences in the physics and f in the forcing. The ensemble is described in detail in Mikolajewicz et al. (2024). This particular simulation corresponds to the simulation D2.1 in Mikolajewicz et al. (2024).

Note that the time runs with an offset of +25001 years, meaning that the year range 1-25000 of the data set equals the years 25000-1 BP (Before Present). All data are global datasets with the exception of output from the ice sheet model where regional northern hemispheric and southern hemispheric domains are used.