Correction to: Scientific Reports https://doi.org/10.1038/s41598-020-67154-8, published online 09 July 2020


The original version of this Article contained errors.


The data presented in de la Vega, E., Chalk, T.B., Wilson, P.A. et al. (2020) is based on a combination of δ11B-derived CO2 proxy data from Martínez-Botí, M. et al. (2015)1 and newly generated data, both from contiguous material gathered from ODP Site 999 in the Caribbean Sea.


The original version of this Article contains an error where the CO2 data reported from Martínez-Botí, M. et al. (2015) was plotted and calculated on its original age model as published. However, the newly generated high resolution δ11B/CO2 data was plotted and calculated on an updated age model, leading to a small number of points falling out of strict stratigraphic sequence.


The age model of the data from Martínez-Botí, M. et al. (2015) has now been updated to our new age model and has been corrected in all Figures, calculations, and Tables affected. This has impacted Figures 1 and 2, the Supplementary Figures 2, 3, 5, 6 and 7, and Supplementary Table 1, 2 and 3. The original Figures 1 and 2, and Supplementary Information file are provided below.

Figure 1
figure 1

Top panel: Current CO2 estimates from boron isotopes across the Plio-Pleistocene boundary. G. ruber data in red circles from site 999 (Martinez-Boti et al.13), T. sacculifer in red squares from site 999 (Seki et al.14), blue squares from site 926 (Sosdian et al.20) and pale red squares from site 999 (Bartoli et al.12). Bottom Panel: δ18O from benthic foraminifer Cibicidoides wuellerstorfi at ODP Site 999 (blue circles) with a 5 point running mean (this study and ref.53) compared to the benthic isotope stack of ref.19. M2 glacial and early Pleistocene strong glacials are highlighted in blue for context. Interglacial KM5 and KM5c are highlighted in yellow and orange, respectively. Note that there are no estimates for the M2 glacial and very few across the mid-Piacenzian warm period (mPWP), the low resolution of previous studies makes pin-pointing individual interglacials such as the KM5c future analogue difficult. These studies also differ in their estimates of Mid-Piacenzian CO212,13,14,20.

Figure 2
figure 2

Top panel: Red circles and lines show δ11B-derived CO2 data from Globigerinoides ruber at ODP Site 999 (this study and Martinez-Boti et al.13, Chalk et al.18), red squares are Trilobatus sacculifer at ODP 999 (this study and Seki et al.14), purple squares are T. sacculifer from ODP 668 (Honisch et al.23) and blue squares are T. sacculifer from ODP 926 (Sosdian et al.20). Black solid line shows ice core-derived CO2 from ref.58. Left; Late Pleistocene CO2 from boron isotopes14,18,23,52 and ice core data. Also shown are CO2 projections in line with RCP8.5 at current emission rates to the year 2040 (black broken line). Middle column; MPT CO2 estimates18,23 including disturbed ice estimates24,25 (Note: age adjusted for scale). Right; mPWP estimates of CO2 (this study combined with Martinez-Boti et al.13), new data from T. sacculifer is shown in red squares and shows no offset from G. ruber estimates. Second panel: Time periods as above, LR04 and ODP 999 δ18O from C. wuellerstorfi18,19,53. Third panel: Iron mass accumulation rate from the Southern Atlantic ODP Site 109028. Fourth panel: % Northern Component Water (NCW) estimated from δ13C in benthic foraminifera (grey) and ɛNd from fish debris (dark green) in the deep North Atlantic (core U131331). Note the lag of ocean circulation and CO2 relative to the M2 glaciation31.


This change does not impact the overall conclusion of the article, however, the average CO2 values during KM5c and the CO2 range during the specific intervals of the study slightly vary from the update.


The original Article and accompanying Supplementary Information files have been corrected.