Bulk carbonate oxygen isotopes: Oxygen isotopes were measured using a ThermoFisher GasBench II equipped with a CTC Combi-Pal autosampler and interfaced via continuous flow with a ThermoFisher MAT 253 isotope ratio mass spectrometer (IRMS). About 200-300 mg of bulk sample was used. Analytical precision is estimated to be ±0.08% for d18O of an internal standard (Carrara marble). Age model: The age model was developed by aligning the abrupt transitions of d18O record of G. bulloides in core MD01-2444 to warming and cooling events in the NGRIP d18Oice record, presented on: (i) the GICC05 timescale [Svensson, A. et al. A 60,000 year Greenland stratigraphic ice core chronology. Clim. Past 4, 47–57 (2008)]; and (ii) the WD2014: modified timescale,...