preindustrial Control experiment to be used in VolMIP analyses. The piControl experiment is the CMIP6-DECK piControl experiment described in Eyring et al. (2016). piControl provides initial climate states that are sampled to start most of VolMIP experiments (Zanchettin et al., 2016). The dataset contains monthly values of selected variables spatially averaged over four regions. These are the full globe (GL), the Northern Hemisphere extratropics (30°-90°N, NH), the tropics (30°S-30°N, TR), and the Southern Hemisphere (30°-90°S, hereafter SH). The considered variables have the following cmor names: hfls, hfss, pr, rlds, rldscs, rlus, rlut, rlutcs, rsds, rsdscs, rsdt, rsus, rsut, rsutcs, tas. Additionally, the climate indices NAO and Nino34 are part of the dataset. Considered models are CanESM5, IPSL-CM6A-LR, GISS-E2.1-G, MIROC-ES2L, MPI-ESM1.2-LR (named MPI-ESM-LR in the files of this dataset) and UKESM1. Considered experiments are piControl and volc-pinatubo-full, with initial date and final date as specified for each model in Zanchettin et al. (2021). Different realizations are considered for the participating models depending on availability.

Idealized volcanic-forcing coupled climate model experiment using the 1991 Pinatubo forcing as used in the CMIP6 historical simulations. It is a Tier 1 (mandatory) VolMIP experiment based on a large ensemble of short-term “Pinatubo” climate simulations aimed at accurately estimating simulated responses to volcanic forcing that may be comparable to the amplitude of internal interannual climate variability. Initialization is based on equally distributed predefined states of ENSO (cold/neutral/warm states) and of the North Atlantic Oscillation (NAO, negative/neutral/positive states). Sampling of an eastern phase of the Quasi-Biennial Oscillation (QBO), as observed after the 1991 Pinatubo eruption, is preferred for those models that spontaneously generate such mode of stratospheric variability. VIRF diagnostics must be calculated for this experiment for the whole integration and for all ensemble members, as these are required for the “volc-pinatubo-strat”/“surf” experiments. A minimum length of integration of 3 years is requested. Details about the experiment are provided by Zanchettin et al. (2016). The dataset contains monthly values of selected variables spatially averaged over four regions. These are the full globe (GL), the Northern Hemisphere extratropics (30°-90°N, NH), the tropics (30°S-30°N, TR), and the Southern Hemisphere (30°-90°S, hereafter SH). The considered variables have the following cmor names: hfls, hfss, pr, rlds, rldscs, rlus, rlut, rlutcs, rsds, rsdscs, rsdt, rsus, rsut, rsutcs, tas. Additionally, the climate indices NAO and Nino34 are part of the dataset. Considered models are CanESM5, IPSL-CM6A-LR, GISS-E2.1-G, MIROC-ES2L, MPI-ESM1.2-LR (named MPI-ESM-LR in the files of this dataset) and UKESM1. Considered experiments are piControl and volc-pinatubo-full, with initial date and final date as specified for each model in Zanchettin et al. (2021). Different realizations are considered for the participating models depending on availability.

Ensemble of MPI-ESM1-2-HR CMIP6 historical simulations with low-pass filtered solar and ozone variability (i.e., using a 33-years running-mean). The simulations are performed within the BMBF project "Solar contribution to climate change on decadal to centennial timescales" (SOLCHECK) of the "Role of the middle atmosphere in climate" (ROMIC II: https://romic2.iap-kborn.de/en/romic/strategy). The experimental setup is identical to the MPI-ESM1-2-HR historical CMIP6 simulations except for the solar and ozone variability.

For transient-deglaciation-prescribed-glac1d-methane (r1i1p1f1): The transient experiment was performed with the Max Planck Institute for Meteorology Earth System Model version 1.2 in coarse resolution (MPI-ESM-CR). The model 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°. The model has been extended with a full methane cycle, including a simplified atmospheric chemistry, as described in Kleinen et al. (2020), Kleinen et al. (2021) and Kleinen et al. (2022). For the experiment, the model was integrated from a glacial state at 23 ka to the year 1950 with prescribed atmospheric greenhouse gas concentrations (Köhler et al., 2017) and insolation (Berger & Loutre, 1991). Ice sheets and surface topographies were prescribed from GLAC-1D (Tarasov et al., 2012) reconstructions. All forcing fields are updated every 10 years of the simulations and initiate changes in the topography, glacier mask, river pathways, ocean bathymetry, and land-sea mask. For transient-deglaciation-prescribed-glac1d-methane (r1i1p1f2): This experiment is derived from the r1i1p1f1 experiment. It is branched off at 18 ka BP, and starting at 15.2 ka BP, meltwater from the Laurentide ice sheet was removed from the system and stored. The accumulated meltwater was released over a period of 1200 years starting at 12.8 ka BP. This induced a collapse of the AMOC, leading to climatic conditions very similar to the Younger Dryas cold period seen in climate reconstructions. After the end of the meltwater release, the circulation recovered quickly and climatic conditions converged with the r1i1p1f1 experiment. Important: Please be aware that due to CMORization constraints data set time values concerning to "Before Present(BP)" time got a year offset (+25001). Means for example year range 2001 to 25000 model time is equal to 23000 to 1 BP.

The transient simulation was performed with the Max Planck Institute for Meteorology Earth System Model version 1.2 in coarse resolution (MPI-ESM-CR). The model 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°. The applied setup was introduced in detail in Kapsch et al. (2021 and 2022). For the experiment, the model was integrated 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). Ice sheets and surface topographies were prescribed from GLAC-1D (Tarasov et al., 2012) reconstructions. All forcing fields are updated every 10 years of the simulations and initiate changes in the topography, glacier mask, river pathways, ocean bathymetry, and land-sea mask. Meltwater from ice sheets is calculated as the temporal derivative of ice thickness at grid points covered by grounded ice sheets. The derived meltwater is then distributed by the hydrological discharge model and finally released into the ocean as freshwater. The experiment was performed as part of a model ensemble that contains simulations differing in terms of their tuning parameters, parameterizations and bug fixes. The ensemble is described in detail in the Supplementary Material of Kapsch et al. (2022). The current simulation refers to GLAC1D_P1 in Kapsch et al. (2022) and is significantly warmer during the glacial than model versions P2 and P3, indicating a weaker sensitivity to greenhouse gas changes. The glacial cooling over the North Atlantic is more realistic in P2 and P3 than in P1. P2 and P3 are also more sensitive to small changes in the meltwater forcing than P1. The experiment was performed as part of a model ensemble that contains simulations differing in terms of their tuning parameters, parameterizations, bug fixes and 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 the Supplementary Material of Kapsch et al. (2022).

The transient simulation was performed with the Max Planck Institute for Meteorology Earth System Model version 1.2 in coarse resolution (MPI-ESM-CR). The model 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°. The applied setup was introduced in detail in Kapsch et al. (2021 and 2022). For the experiment, the model was integrated 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). Ice sheets and surface topographies were prescribed from ICE-6G (Peltier et al., 2015) reconstructions. All forcing fields are updated every 10 years of the simulations and initiate changes in the topography, glacier mask, river pathways, ocean bathymetry, and land-sea mask. Meltwater from ice sheets is calculated as the temporal derivative of ice thickness at grid points covered by grounded ice sheets. The derived meltwater is then distributed by the hydrological discharge model and finally released into the ocean as freshwater. The experiment was performed as part of a model ensemble that contains simulations differing in terms of their tuning parameters, parameterizations and bug fixes. The ensemble is described in detail in the Supplementary Material of Kapsch et al. (2022). The current simulation refers to ICE6G_P1 in Kapsch et al. (2022) and is significantly warmer during the glacial than model versions P2 and P3, indicating a weaker sensitivity to greenhouse gas changes. The glacial cooling over the North Atlantic is more realistic in P2 and P3 than in P1. P2 and P3 are also more sensitive to small changes in the meltwater forcing than P1. The experiment was performed as part of a model ensemble that contains simulations differing in terms of their tuning parameters, parameterizations, bug fixes and 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 the Supplementary Material of Kapsch et al. (2022).

The transient simulation was performed with the Max Planck Institute for Meteorology Earth System Model version 1.2 in coarse resolution (MPI-ESM-CR). The model 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°. The applied setup was introduced in detail in Kapsch et al. (2021 and 2022). For the experiment, the model was integrated 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). Ice sheets and surface topographies were prescribed from GLAC-1D (Tarasov et al., 2012) reconstructions. All forcing fields are updated every 10 years of the simulations and initiate changes in the topography, glacier mask, river pathways, ocean bathymetry, and land-sea mask. Meltwater from ice sheets is calculated as the temporal derivative of ice thickness at grid points covered by grounded ice sheets. The derived meltwater is then distributed by the hydrological discharge model and finally released into the ocean as freshwater. The experiment was performed as part of a model ensemble that contains simulations differing in terms of their tuning parameters, parameterizations and bug fixes. The ensemble is described in detail in the Supplementary Material of Kapsch et al. (2022). The current simulation refers to GLAC1D_P2 in Kapsch et al. (2022) and is significantly colder during the glacial than model version P1, indicating a stronger sensitivity to greenhouse gas changes. The glacial cooling over the North Atlantic is more realistic in P2 and P3 than in P1. P2 and P3 are also more sensitive to small changes in the meltwater forcing than P1. The experiment was performed as part of a model ensemble that contains simulations differing in terms of their tuning parameters, parameterizations, bug fixes and 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 the Supplementary Material of Kapsch et al. (2022).

The transient simulation was performed with the Max Planck Institute for Meteorology Earth System Model version 1.2 in coarse resolution (MPI-ESM-CR). The model 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°. The applied setup was introduced in detail in Kapsch et al. (2021 and 2022). For the experiment, the model was integrated 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). Ice sheets and surface topographies were prescribed from ICE-6G (Peltier et al., 2015) reconstructions. All forcing fields are updated every 10 years of the simulations and initiate changes in the topography, glacier mask, river pathways, ocean bathymetry, and land-sea mask. Meltwater from ice sheets is calculated as the temporal derivative of ice thickness at grid points covered by grounded ice sheets. The derived meltwater is then distributed by the hydrological discharge model and finally released into the ocean as freshwater. The experiment was performed as part of a model ensemble that contains simulations differing in terms of their tuning parameters, parameterizations and bug fixes. The ensemble is described in detail in the Supplementary Material of Kapsch et al. (2022). The current simulation refers to ICE6G_P2 in Kapsch et al. (2022) and is significantly colder during the glacial than model version P1, indicating a stronger sensitivity to greenhouse gas changes. The glacial cooling over the North Atlantic is more realistic in P2 and P3 than in P1. P2 and P3 are also more sensitive to small changes in the meltwater forcing than P1. The experiment was performed as part of a model ensemble that contains simulations differing in terms of their tuning parameters, parameterizations, bug fixes and 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 the Supplementary Material of Kapsch et al. (2022).

The transient simulation was performed with the Max Planck Institute for Meteorology Earth System Model version 1.2 in coarse resolution (MPI-ESM-CR). The model 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°. The applied setup was introduced in detail in Kapsch et al. (2021 and 2022). For the experiment, the model was integrated 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). Ice sheets and surface topographies were prescribed from ICE-6G (Peltier et al., 2015) reconstructions. All forcing fields are updated every 10 years of the simulations and initiate changes in the topography, glacier mask, river pathways, ocean bathymetry, and land-sea mask. Meltwater from ice sheets is calculated as the temporal derivative of ice thickness at grid points covered by grounded ice sheets. The derived meltwater is then distributed by the hydrological discharge model and finally released into the ocean as freshwater. The experiment was performed as part of a model ensemble that contains simulations differing in terms of their tuning parameters, parameterizations and bug fixes. The ensemble is described in detail in the Supplementary Material of Kapsch et al. (2022). The current simulation refers to ICE6G_P3 in Kapsch et al. (2022) and is significantly colder during the glacial than model version P1, indicating a stronger sensitivity to greenhouse gas changes. The glacial cooling over the North Atlantic is more realistic in P2 and P3 than in P1. P2 and P3 are also more sensitive to small changes in the meltwater forcing than P1. The experiment was performed as part of a model ensemble that contains simulations differing in terms of their tuning parameters, parameterizations, bug fixes and 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 the Supplementary Material of Kapsch et al. (2022).

The transient simulation was performed with the Max Planck Institute for Meteorology Earth System Model version 1.2 in coarse resolution (MPI-ESM-CR). The model 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°. The applied setup was introduced in detail in Kapsch et al. (2021 and 2022). For the experiment, the model was integrated 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). Ice sheets and surface topographies were prescribed from GLAC-1D (Tarasov et al., 2012) reconstructions. All forcing fields are updated every 10 years of the simulations and initiate changes in the topography, glacier mask, river pathways, ocean bathymetry, and land-sea mask. Meltwater from ice sheets is calculated as the temporal derivative of ice thickness at grid points covered by grounded ice sheets. The derived meltwater is then distributed by the hydrological discharge model and finally released into the ocean as freshwater. The experiment was performed as part of a model ensemble that contains simulations differing in terms of their tuning parameters, parameterizations and bug fixes. The ensemble is described in detail in the Supplementary Material of Kapsch et al. (2022). The current simulation refers to GLAC1D_P3 in Kapsch et al. (2022) and is significantly colder during the glacial than model version P1, indicating a stronger sensitivity to greenhouse gas changes. The glacial cooling over the North Atlantic is more realistic in P2 and P3 than in P1. P2 and P3 are also more sensitive to small changes in the meltwater forcing than P1. The experiment was performed as part of a model ensemble that contains simulations differing in terms of their tuning parameters, parameterizations, bug fixes and 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 the Supplementary Material of Kapsch et al. (2022).