# Exceptionally high biosphere productivity at the beginning of Marine Isotopic Stage 11

Margaux Brandon, Amaelle Landais, Stéphanie Duchamp-Alphonse, Violaine

Favre, Léa Schmitz, Héloïse Abrial, Frédéric Prié, Thomas Extier, Thomas

Blunier

► To cite this version:

Margaux Brandon, Amaelle Landais, Stéphanie Duchamp-Alphonse, Violaine Favre, Léa Schmitz, et al.. Exceptionally high biosphere productivity at the beginning of Marine Isotopic Stage 11. Nature Communications, 2020, 11 (1), 10.1038/s41467-020-15739-2. http://dx.doi.org/10.1038/s41467-020-15739-2

HAL Id: hal-02765978

https://hal.science/hal-02765978

Submitted on 17 Sep 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
ARTICLE

https://doi.org/10.1038/s41467-020-15739-2

OPEN

<figure><img src="image_2.png" /></figure>

Check for updates

# Exceptionally high biosphere productivity at the beginning of Marine Isotopic Stage 11

Margaux Brandon<sup>1,2✉</sup>, Amaelle Landais<sup>1</sup>, Stéphanie Duchamp-Alphonse<sup>2</sup>, Violaine Favre<sup>1</sup>, Léa Schmitz<sup>1</sup>, Héloïse Abrial<sup>1</sup>, Frédéric Prié<sup>1</sup>, Thomas Extier<sup>1</sup> & Thomas Blunier<sup>3</sup>

Significant changes in atmospheric CO₂ over glacial-interglacial cycles have mainly been attributed to the Southern Ocean through physical and biological processes. However, little is known about the contribution of global biosphere productivity, associated with important CO₂ fluxes. Here we present the first high resolution record of Δ¹⁷O of O₂ in the Antarctic EPICA Dome C ice core over Termination V and Marine Isotopic Stage (MIS) 11 and reconstruct the global oxygen biosphere productivity over the last 445 ka. Our data show that compared to the younger terminations, biosphere productivity at the end of Termination V is 10 to 30 % higher. Comparisons with local palaeo observations suggest that strong terrestrial productivity in a context of low eccentricity might explain this pattern. We propose that higher biosphere productivity could have maintained low atmospheric CO₂ at the beginning of MIS 11, thus highlighting its control on the global climate during Termination V.

<footnote n="1">¹Université Paris-Saclay, CNRS, CEA, UVSQ, Laboratoire des sciences du climat et de l'environnement, 91191 Gif-sur-Yvette, France. ²Université Paris-Saclay, CNRS, GEOPS, 91405 Orsay, France. ³Centre for Ice and Climate, Niels Bohr Institute, University of Copenhagen, Copenhagen, Denmark.
✉email: margaux.brandon@universite-paris-saclay.fr</footnote>
The largest pre-anthropogenic changes in atmospheric CO₂ concentration of the last 800,000 years are observed during deglaciations, with increases of up to 100 ppm in a few thousand years¹. Oceanic carbon reservoir is broadly suspected to play a central role in these atmospheric CO₂ increases. Leading hypotheses invoke CO₂ degassing from the ocean induced by more vigorous convection² and concomitant decrease of the net organic matter export in the Southern Ocean³, enhanced exchanges between ocean surface and atmosphere due to sea-ice melting⁴, and increased sea-surface temperature⁵. Modelling studies simulating changes in oceanic processes are quite controversial⁶,⁷ but in all cases, they are not able to explain the full increase of atmospheric CO₂ rises during deglaciations. In parallel, terrestrial primary productivity and carbon stocks increase during deglaciations, thus acting as a significant CO₂ land sink⁸,⁹ so that an important additional source of CO₂ is required to explain the entire deglacial CO₂ pattern.

Quantifying changes in the carbon cycle over deglaciations relies on data compilation and modelling studies to estimate both carbon stocks and carbon fluxes⁸. Some information on carbon stocks can be obtained from δ¹³C of carbonates using a wealth of data obtained from marine sediments¹⁰ and δ¹³C of atmospheric CO₂¹¹,¹². In parallel, information on the evolution of past global productivity, a major component of carbon flux, is very sparse and often limited to the last deglaciation⁸,⁹.

In the ocean, past changes in biological carbon pump are best represented by changes in buried organic biomarkers combined with marine Total Organic Carbon/Particulate Inorganic Carbon ratio (TOC/CaCO₃), suggested to reflect the C-rain ratio¹³. However, despite their accuracy to provide biological export production, only one site combine such records in the Southern Ocean, for the last 800 ka, so far¹⁴.

On continents, pollen counting and sedimentary TOC are useful for biosphere productivity reconstruction¹⁵,¹⁶ but they are unfortunately indirect and rely on the use of biosphere models⁹. Moreover, similarly to oceanic records, these observations provide regional records that are not easy to use for documenting the past global carbon cycle. Ciais et al.⁸ proposed to use the isotopic composition of oxygen of atmosphere (δ¹⁸O_atm) as a tracer for terrestrial biosphere productivity. However, this proxy is a complex tracer being influenced by hydrological cycle at first order¹⁷,¹⁸ and its use as a quantitative tool for productivity reconstruction depends on the exact determination of associated fractionation factors in the water and biosphere cycles¹⁹.

A total estimate of the global biospheric fluxes and their temporal variations can be obtained more directly from measurements of Δ¹⁷O of O₂ (ln(δ¹⁷O + 1) – 0.516*ln(δ¹⁸O + 1)) in ice cores²⁰⁻²². This method provides O₂ fluxes and the conversion from O₂ to CO₂ fluxes can be done from the stoichiometry of the biological processes of photosynthesis and respiration²³. Δ¹⁷O of O₂ measures the variation of the triple isotopic composition of atmospheric O₂ with respect to modern oxygen so that by definition, Δ¹⁷O of O₂ is nil today. Previous experimental studies showed that Δ¹⁷O of O₂ increases in a closed biospheric system when the exchanges with the stratosphere are prevented: the biological productivity leads to Δ¹⁷O of O₂ increase while photochemical reactions occurring in the stratosphere have the effect of decreasing Δ¹⁷O of O₂²⁰. For paleoproductivity reconstructions, measurements of the evolution of Δ¹⁷O of O₂ in ancient air trapped in the Vostok and GISP2 ice cores already provided information on the evolution of global biosphere productivity, over the last 400 ka²¹,²². The results show a systematic larger productivity during interglacial than during glacial periods with interglacial levels remaining close to the current biosphere productivity.

Over the last nine deglaciations, Termination V (433–426 ka on ice core records on the latest AICC2012 chronology²⁴) is probably the most intriguing. This Termination is framed by the particularly long and strong glacial Marine Isotopic Stage 12 (MIS 12), followed by the long and warm interglacial MIS 11 (426–398 ka on AICC2012). This is the first Termination after the Mid-Brunhes event marking a fundamental change in the climate system from mild to warm periods with associated lower to higher CO₂ concentrations. Termination V is also occurring in a particular orbital context of low eccentricity, which is known to have an influence on the carbon cycle as observed in δ¹³C oceanic records²⁵. On the continents, pollen data¹⁶,²⁶ suggest a strong and long increase in terrestrial productivity during MIS 11. In the Ocean, MIS 11 displays an unusual increase in carbonate storage in low²⁷ and high-latitude environments²⁸. Yet, while it is clearly associated with a major phase in coral reef expansion²⁹ and a climax in calcareous phytoplankton productivity respectively³⁰, the impact of this large carbonate production on atmospheric pCO₂ is not understood. Therefore, the biosphere productivity fluxes during Termination V and MIS 11 need to be investigated.

Here we present the first measurements of the triple isotopic composition of atmospheric oxygen (Δ¹⁷O of O₂) in the Antarctic EPICA Dome C ice core over Termination V. Using these measurements and new correction factors compared to the previous record of Blunier et al.²², we reconstruct the oxygen fluxes associated with biosphere productivity over this particularly strong Termination. The biosphere productivity over Termination V and the beginning of MIS 11 is found to be 10–30% higher than productivity over the pre-industrial period, an exceptional value never encountered over the last four interglacial periods. The ice core δ¹⁸O_atm (or δ¹⁸O of O₂) record, and terrestrial and oceanic records related to biosphere productivity are used to discuss the relative contribution of changes in oceanic vs terrestrial biosphere fluxes to the atmospheric CO₂ rise during Termination V. We show that this productivity peak is most probably due to an increase of the terrestrial productivity during this period favoured by a particular context of low eccentricity. We propose that such strong productivity occurring concomitantly with an exceptional productivity carbonate peak in the marine realm, plays a role in maintaining the CO₂ level at a relatively low level at the beginning of MIS 11.

## Results

**Record of Δ¹⁷O of O₂ during Termination V.** The Antarctic EPICA Dome C (EDC) ice core (75°06'S; 123°21'E, 3233 m above sea level) covers the last 800 ka. We measured Δ¹⁷O of O₂ over the depth range 2735–2797 m, corresponding to an age between 405.7 and 444.1 ka with an average resolution of 780 years (50 samples) (Fig. 1). The raw data were corrected for several effects (Methods): values of the atmospheric air Δ¹⁷O of O₂ were obtained by corrections for gravitational fractionation, air bubble trapping fractionation and gas loss fractionation. Compared to the previous records²⁰⁻²², several corrections were addressed for the first time here (gas loss and air bubble trapping effects, see Supplementary Table 1). To check the coherency between the previous records that did not take into account these corrections, and the new one obtained on EDC, we measured Δ¹⁷O of O₂ over Termination II on the EDC ice core and compared our corrected Δ¹⁷O of O₂ curve with the previous record obtained by Blunier et al.²² on the Vostok ice core (Fig. 1). The two records are in agreement, displaying the same 51 ppm decrease over Termination II without any shift in the mean Δ¹⁷O of O₂ value because of compensating effects in the different corrections (see Methods). We thus conclude that our new record over Termination V can directly be used to complete the previous Δ¹⁷O of O₂ record from
<figure><img src="image_3.png" /><figcaption>**Fig. 1 Record of Δ¹⁷O of O₂ over the last 450 ka compared to evolution of CO₂.** a Δ¹⁷O of O₂ (light green: record presented in Blunier et al.²¹,²² covering the last 400 ka; dark green: new record with error bars showing the standard deviation of ± 6 per meg). b atmospheric CO₂ variations over the last 450 ka⁶⁹. The ice core records are presented on the AICC2012 timescale²⁴,⁷⁰. The grey shadow bars represent the period of rapid increase in atmospheric CO₂.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Fig. 2 Comparison of CO₂ and Δ¹⁷O of O₂ evolutions over the last five terminations.** CO₂⁶⁹ (black); Δ¹⁷O of O₂ combining new data (Termination II and Termination V, dotted green line) and data from Blunier et al.²¹,²² (solid green line). The red arrows indicate the strongest differences between the CO₂ and Δ¹⁷O of O₂ evolutions during interglacial periods. The grey shadow bars correspond to the main increase in atmospheric CO₂ over terminations.</figcaption></figure>

Blunier et al.²² and provide a full record of Δ¹⁷O of O₂ over the last five deglaciations.

Termination V displays a decrease in Δ¹⁷O of O₂ (56 ppm) of similar amplitude than the younger deglaciations (53, 35, 51 and 55 ppm for Terminations IV, III, II and I, respectively). The main difference with previous terminations is on the duration and timing of the Δ¹⁷O of O₂ decreases. For the four youngest terminations, the Δ¹⁷O of O₂ decreases are in parallel to the CO₂ increases²¹,²². Over Termination V, the main decrease of Δ¹⁷O of O₂ is more than twice as long as the CO₂ increase (from 434.8 to 410.2 ka and from 434.8 to 424.8 ka, i.e. 24.6 ka and 10 ka respectively, Fig. 2). This result contrasts with previous studies covering the last four deglaciations²¹,²² where decrease of Δ¹⁷O of O₂ and increase of atmospheric CO₂ concentration are synchronous (Supplementary Fig. 1).

Focusing on the last five deglaciations, only two cases can be identified with a significant deviation between the CO₂ and the Δ¹⁷O of O₂ evolutions during the deglaciations (red arrows on Fig. 2): the first case at the very end of Termination III, between 240 and 245 ka, when the short peak in atmospheric CO₂ concentration does not correspond to a sharp decrease in Δ¹⁷O of O₂ and the second case over Termination V, between 430 and 415 ka, when the maximum in CO₂ concentration is reached at 425 ka, more than 10 ka before the minimum in Δ¹⁷O of O₂ recorded at 415 ka. Contrary to the case of Termination III where the deviation of the Δ¹⁷O of O₂ signal from the CO₂ optimum is

only based on three points, the mismatch between CO₂ and Δ¹⁷O of O₂ evolutions over Termination V relies on 40 data points of Δ¹⁷O of O₂. When looking at the anticorrelation between Δ¹⁷O of O₂ and CO₂ evolutions over the last five deglaciations, Termination V stands out with a significantly much lower slope for the relationship between Δ¹⁷O of O₂ and CO₂ than the last four Terminations (Supplementary Fig. 1). In addition, to highlight the specificity of Termination V, our new data hence show that the systematic use of Δ¹⁷O of O₂ to infer CO₂ values when CO₂ measurements are not available³¹ is not always a valid approach (Supplementary Fig. 1).

The general anti-correlations between CO₂ and Δ¹⁷O of O₂ over terminations are explained at first order by stratospheric reactions. These reactions are leading to formation of ozone from the dioxygen reservoir and are associated with mass independent isotopic fractionations. These fractionations lead to an isotopic anomaly of ozone that can be transmitted to CO₂. There is eventually a correspondence between the Δ¹⁷O of CO₂ and the Δ¹⁷O of O₂, and it has been assumed that this anomaly transmission is rate-limited by the abundance of CO₂²⁰. As a consequence, an increase in CO₂ in the absence of changes in biosphere productivity should be associated with a decrease in Δ¹⁷O of O₂. The increase of biosphere productivity comes as an opposite effect: it produces oxygen with mass dependent isotopic fractionation which has the effect of increasing Δ¹⁷O of O₂. The similar CO₂ and Δ¹⁷O of O₂ anticorrelation observed over the
latest terminations is hence an indication that biosphere productivity variations are parallel to the variations of atmospheric concentration in CO₂ and that a constant relationship exists between the amplitudes of biosphere productivity variations and the amplitudes of atmospheric CO₂ variations. As a consequence, the deviation of the classical CO₂ vs Δ¹⁷O of O₂ anticorrelation toward unexpectedly high Δ¹⁷O of O₂ values should be interpreted as an anomalous increase in biosphere productivity compared to the situation for other terminations. We thus propose that the anomalously strong difference in the CO₂ and Δ¹⁷O of O₂ evolutions over the first part of MIS 11 at the end of Termination V is associated with a stronger biosphere productivity compared to the onset of the 4 last interglacial periods. We are aware that alternative possibilities can be suggested such as a decrease in the stratospheric ozone in agreement with relatively high N₂O atmospheric concentration (sink of ozone) during MIS 11³². However, similarly high N₂O atmospheric concentration is observed over the first part of MIS 9 without associated decoupling between atmospheric CO₂ increase and Δ¹⁷O of O₂ decrease. In the following, we thus favour the interpretation in terms of an anomalous change of paleoproduction during Termination V compared to younger terminations.

**Reconstruction of the global oxygen biosphere productivity.**
For a more quantitative reconstruction of paleoproduction based on the classical use of Δ¹⁷O of O₂, we followed the approach detailed in Landais et al.³³ which relies on estimates of the Δ¹⁷O of O₂ values produced by the terrestrial and oceanic biospheres. Such values are based on previously determined values of oxygen fractionation coefficients for the different biologic processes of oxygen uptake and production³⁴. Because of different proportions of C3 vs C4 plants or terrestrial vs oceanic biosphere productivity ratios, the global Δ¹⁷O of O₂ value produced by the biosphere can display strong variations between glacial and interglacial periods.

The reconstructed biosphere productivity using the Landais et al.³³ formulation and uncertainties is detailed in Supplementary Table 2 for different time periods of interest within this study, i.e. the pre-industrial, the LGM and the MIS 11. On average, we obtain a biosphere productivity at the beginning of MIS 11 which is 17 % higher than during MIS 1 and a biosphere productivity during the LGM which is 31 % lower than during MIS 1, a result in agreement with output of the IPSL coupled model equipped with vegetation and ocean productivity modules (Supplementary Fig. 4).

As in Landais et al.³³, it is shown in Supplementary Table 2 that the uncertainty in fractionation coefficients as well as on the ratio of oceanic to terrestrial productivity is leading to uncertainties in the reconstructed past global productivity. The fractionation coefficients are based on physical properties and are hence not expected to vary with time so that the possible bias on these coefficients should apply to the different periods. On opposite, the ratio of oceanic to terrestrial productivity is expected to vary with time and is a large source of uncertainty. We thus present in Supplementary Table 2 calculations performed with the largest range of possible ratios of oceanic to terrestrial productivity estimated in Landais et al.³³.

Since the publication of Landais et al.³³, new estimates of fractionation coefficients within the oxygen cycle are available¹⁹,³⁵,³⁶ and influence of the water cycle organisation has also been suggested²². As a consequence, we performed three additional types of sensitivity tests to better estimate the uncertainties related to fractionation processes within the photosynthesis³⁷, respiration³⁶ and water cycle²² (see SOM for details).

First, we address the uncertainty in photosynthesis fractionation. Some marine species show fractionation during marine

photosynthesis³⁷ while it was assumed in Landais et al.³³ that photosynthesis does not fractionate³⁸. We performed sensitivity tests with the largest fractionation effects during marine photosynthesis observed in Eisenstadt et al.³⁷ and obtained a MIS 11 productivity level decreased by 3% with such fractionation effect.

The second uncertainty we tested is the uncertainty in respiration fractionation. Recent studies have highlighted large variations in the relationship between δ¹⁷O and δ¹⁸O during respiration linked to temperature variations³⁶. When taking into account the maximum effect, i.e. a decrease of 0.005 for the slope of the relationship between ln(δ¹⁷O + 1) and ln(δ¹⁸O + 1) during respiration, we end up with a resulting biosphere productivity reconstruction 16% higher during MIS 11 than the average situation.

Then, we estimated the uncertainty in the fractionation within water cycle. The general assumption in Landais et al.³³ was that the relationship between δ¹⁷O and δ¹⁸O of water (meteoric water line) remains the same over glacial-interglacial cycles. However, measurements in the Vostok ice core have shown that the ¹⁷O-excess, defined as ¹⁷O-excess = ln(1 + δ¹⁷O) - 0.528*ln(1 + δ¹⁸O), is lower by up to 20 ppm during the last glacial maximum³⁹. While it has been shown that this is a local effect⁴⁰, we still performed a sensitivity test with decreasing ¹⁷O-excess of all continental meteoric waters by 20 ppm during glacial periods and increasing ¹⁷O-excess by 10 ppm during MIS 11 with respect to our current interglacial period. Taking into account such variations, the reconstructed biosphere productivity does not differ by more than 6% from the case without change in meteoric water ¹⁷O-excess (Supplementary Table 2).

Finally, we also took into account the uncertainty in the model used by Landais et al.³³ by combining the above reconstruction with an alternative reconstruction of global oxygen productivity using an alternative model, the one of Blunier et al.²², forced by our Δ¹⁷O of O₂ data (Supplementary Fig. 2).

Combining our different sensitivity studies for a global reconstruction of biosphere productivity (Fig. 3), we find that the global productivity is reduced by 10–40% during glacial periods compared to interglacials. At the beginning of most interglacial periods, the global productivity remains close to the pre-industrial level. The only exception to this general behaviour is the strong oxygen biosphere productivity at the end of Termination V reaching values 10–30% higher than during the pre-industrial period.

## Discussion

Complementary information on the origin and specificity of the Δ¹⁷O of O₂ signal over Termination V can be obtained from the ice core δ¹⁸O_atm record over the last 800 ka¹⁸ (Fig. 3). δ¹⁸O_atm is a complex parameter resulting from both biosphere productivity and low-latitude water cycle⁴¹. In particular, it has been shown that Weak Monsoon Intervals observed during Heinrich events lead to increases in δ¹⁸O_atm via changes in the low-latitude water cycle⁴². Another way to increase the δ¹⁸O_atm is to increase the ratio between terrestrial and oceanic biosphere productivity⁸,⁴¹.

As for Δ¹⁷O of O₂, δ¹⁸O_atm shows a particular behaviour over Termination V. First, while δ¹⁸O_atm variations during other terminations of the last 800 ka do not exhibit variations larger than 1.5 ‰¹⁸, the δ¹⁸O_atm decrease over Termination V is of 1.8 ‰. The highest level of δ¹⁸O_atm (1.43 ‰) is reached during Termination V at the beginning of the main increase of productivity inferred from Δ¹⁷O of O₂. In contrast with the other interglacial periods where the maximum in biosphere productivity corresponds to a minimum of δ¹⁸O_atm, there is no minimum of δ¹⁸O_atm observed at the beginning of MIS 11 (~420 ka) (Fig. 3).
<figure><img src="image_5.png" /><figcaption>**Fig. 3** Global and regional productivity records since 450 ka. **a** Eccentricity<sup>71</sup>. **b** CO<sub>2</sub> record (ppmv)<sup>69</sup>. **c** Reconstructed envelop for the ratio between global biosphere productivity and pre-industrial biosphere productivity as inferred from Δ<sup>17</sup>O of O<sub>2</sub> at age t (interpolation to 200 years and 101 binomial smoothing with Igor software). **d** δ<sup>18</sup>O<sub>atm</sub> record from the EDC ice core<sup>18,72-76</sup>. **e** LR04 Sea level stack (m compared to present) on AICC2012 timescale<sup>77</sup>. **f** Trees, shrubs and Picea pollen (%)<sup>16</sup>; **g** Si/Ti ratio<sup>16</sup>; a proxy of biogenic silica normalised to detrital, reflecting the changes in diatom productivity in the lake. **h** Pollen abundance (grains.cm<sup>-3</sup>) of ODP Site 646<sup>26</sup>, **i** Arboreal and Quercus robur pollen records (%) from Lake Ohrid, Balkan Peninsula<sup>46,47</sup>. **j** Alkenone mass accumulation rate (MAR) (µg·m<sup>-2</sup>·year<sup>-1</sup>) and TOC MAR (mg·m<sup>-2</sup>·year<sup>-1</sup>) records at Site PS2489-2/ODP1090<sup>3</sup>. **k** TOC/CaCO<sub>3</sub> ratio at Site PS2489-2/ODP1090, Atlantic Southern Ocean<sup>14</sup>. The horizontal dotted line separates the ice core records presented on the AICC2012 timescale<sup>24,70</sup> (above the line) and the terrestrial and oceanic records presented on the age model of each core (under the line). The grey shadow bars represent the period of rapid increase in atmospheric CO<sub>2</sub>.</figcaption></figure>

The minimum in δ¹⁸O_atm actually occurs toward the end of MIS 11, i.e. around 405 ka. The anomalously high δ¹⁸O_atm at the time of peak biosphere productivity inferred from Δ¹⁷O of O₂ cannot be explained by a weak monsoon period in opposite to an explanation given for high δ¹⁸O_atm levels observed over the last 800 ka⁴³. Indeed, data from other archives rather show a relatively strong East Asia monsoon activity during MIS 11 similarly to what is observed over the last 4 interglacial periods⁴⁴,⁴⁵. Instead, this high δ¹⁸O_atm values should be interpreted as an increase of the terrestrial vs oceanic productivity ratio expressed in flux of oxygen.

In most cases, terrestrial records are too short and not well dated enough to study the dynamic of vegetation and TOC over this termination. Moreover, hiatus of pollen data from oceanic records often prevent to study Termination V with high resolution. The best regional information on the terrestrial productivity

covering the last 445 ka can be recovered from palynological studies of sediments from El'gygytgyn<sup>16</sup>, Ohrid<sup>46,47</sup>, and Tenaghi Philippon<sup>15</sup> lakes, and from the oceanic ODP Site 646<sup>26</sup>. Pollen records can be used as indirect information on terrestrial biosphere productivity knowing the primary productivity associated with each ecosystem<sup>48</sup>. From this, the aforementioned records from the Northern Hemisphere highlight a major increase in terrestrial productivity during Termination V until mid-MIS 11 through the growth of plants associated with high productivity (Picea, Tree, Shrubs, Arboreal pollen, Quercus robur) in place of steppic species, tundra or frozen soils. This increase is the highest of the last 445 ka when focusing on high latitudes records (ODP site 646 at 58°N<sup>26</sup>, Lake El'gygytgyn at 67°N)<sup>16</sup>. At lakes Ohrid<sup>46,47</sup> and Tenaghi Philippon<sup>15</sup>, that are located around 40°N, while the increase in pollen abundance during Termination V is comparable to those observed during younger ones, the
pollen peak lasts longer (Supplementary Fig. 3). In parallel, the very high and long-lasting Si/Ti ratio and TOC percentage recorded in lake sediments¹⁶,⁴⁹ are the consequences of a strong and longer deglacial productivity¹⁵,¹⁶,⁴⁶,⁴⁷ (Fig. 3). These regional indications of terrestrial productivity are in agreement with the increase in global productivity inferred from $\Delta^{17}\text{O}$ of O$_2$ and increase in the terrestrial vs oceanic productivity ratio inferred from $\delta^{18}\text{O}_{\text{atm}}$.

Several factors could explain the higher terrestrial productivity recorded during Termination V and MIS 11 compared to other younger terminations. First, the relatively slow sea level rise over Termination V⁵⁰ leaves some emerged vegetated terrestrial platforms (e.g. Sunda plate). Then, terrestrial biosphere productivity could be favoured in the context of low eccentricity leading to less pronounced seasonality and a longer duration of summer growing season¹⁶. Finally, MIS 11 biosphere productivity can be enhanced at high latitude of the northern Hemisphere because of the strong warming and melting of ice sheet volumes in Northern Hemisphere at mid-MIS 11 leading to more available lands¹⁶,⁴⁹, favouring plant development/productivity. Some also suggest a more humid and warm climate, with higher precipitation compared to MIS 1¹⁶,⁴⁹.

In the Ocean, little is known about deglacial oceanic productivity over the last 800 ka. Increased contents of organic matter (Total Organic Carbon (TOC) and biomarkers)³,¹⁴ in sediments from the Subantarctic Zone of the Southern Ocean have been used to reflect increased downward flux of carbon associated with the photosynthetic biomass and thus, increased soft-tissue pump efficiency. However, they remain sparse and cannot directly be used to infer oceanic primary productivity. Besides, they do not consider the relative contribution of the downward flux of Particulate Inorganic Carbon (mainly CaCO$_3$) produced by calcifying plankton in the sunlit ocean (mainly coccolithophores and planktonic foraminifera), that creates a surface-to-deep alkalinity gradient, causing CO$_2$ to be released back to the atmosphere⁵¹ and that represents the carbonate counter pump. While the biological carbon pump, including both the soft-tissue pump and the carbonate counter pump, is expected to impact the atmospheric pCO$_2$, the soft-tissue pump is the only one to affect O$_2$ fluxes and thus, the $\Delta^{17}\text{O}$ of O$_2$.

To the best of our knowledge, the best regional information on Southern Ocean productivity spanning the last 445 ka is obtained from site PS2489-2/ODP1090³,¹⁴ located within the Atlantic sector of the Subantarctic Zone. This site provides TOC and alkenone mass accumulation rates (MARs) as well as marine TOC/CaCO$_3$ ratio that are helpful to discuss both soft-tissue pump and biological carbon pump efficiencies. Records show decreases in TOC and alkenone MARs, correlated with decreased TOC/CaCO$_3$ during Terminations V to I (Fig. 3). TOC and alkenone records during Termination V display no different trend compared to younger terminations, i.e. they show a less efficient soft-tissue pump probably associated with a decreased biological productivity in the Subantarctic Zone of the Southern Ocean, hence in agreement with the increase of terrestrial vs oceanic productivity over Termination V inferred from $\delta^{18}\text{O}_{\text{atm}}$.

The most striking pattern comes from the TOC/CaCO$_3$ ratio directly linked to biological pump efficiency and exhibiting a peak of very large amplitude over Termination V. This productivity signal goes along the exceptional production and accumulation of the coccolith *Gephyrocapsa caribbeanica* observed worldwide²⁸,⁵², and particularly within the Subantarctic Zone of the Southern Ocean³⁰,⁵³ at the beginning of, and during Termination V. Such exceptional coccolithophore acme could be enhanced by the minimum of eccentricity during this period, causing a longer duration of summer growing season⁵². In parallel, in the low-latitude ocean, this period is marked by a strong increase in reef

carbonate accumulation, probably associated with the sea level rise documented at that time²⁹,⁵⁴. The singular carbonate production by coccolithophore and coral ecosystem is more difficult to interpret straightforwardly in term of increase in biological productivity but supports the exceptional biological context over Termination V. Indeed this worldwide increase in carbonate production and preservation is expected to considerably decrease ocean $[\text{CO}_3^{2-}]^{28,55}$, and thus, significantly increase the atmospheric CO$_2$ concentration³⁰,⁵⁴ over Termination V. However, and opposite to the evolution of atmospheric CO$_2$ concentration during the other interglacial periods, the onset of MIS 11 is not associated with the highest CO$_2$ level of the whole interglacial period which means that a competing effect is at play. Some processes can play a role in the moderation of the pCO$_2$ such as physical pumping of CO$_2$ in the Southern Ocean through solubility effect or ventilation. A recent study provided indication on the change in the surface temperature of the Southern Ocean on the same timescale than atmospheric CO$_2$ variations over the last 720 ka with a link to ocean ventilation⁵⁶. The changes in Southern Ocean temperature and ventilation are not different during Termination V than over the four youngest terminations, therefore showing no evidence of an enhanced physical pumping over this time period. We rather postulate that the competing effect for a moderate CO$_2$ level over the first part of MIS 11 is the exceptional productivity during Termination V leading to a fixation of excess atmospheric CO$_2$ in the biosphere. Moreover, the maximum in atmospheric CO$_2$ concentration is actually observed more than 10 ka after the onset of MIS 11, i.e. at the exact time when global biosphere productivity, as inferred from $\Delta^{17}\text{O}$ of O$_2$, is decreasing. We thus propose that the strong peak of terrestrial productivity over Termination V, mainly driven by terrestrial biosphere, might slow down the increase of atmospheric CO$_2$ concentration induced by the carbonate production.

Our new record of $\Delta^{17}\text{O}$ of O$_2$ over Termination V reveals a very different picture from the evolution observed over the last 4 terminations: $\Delta^{17}\text{O}$ of O$_2$ is not strongly anti-correlated with CO$_2$ which should be interpreted as an anomalous evolution of the biosphere productivity during Termination V and beginning of MIS 11. The global oxygen biosphere productivity is larger by about 10 to 30% at the beginning of MIS 11 compared to the productivity observed at the onset of the younger interglacial periods. Comparison with the $\delta^{18}\text{O}_{\text{atm}}$ record over the last 800 ka and the few available data on continent and ocean shows that this exceptional productivity increase is mainly due to an increase in terrestrial productivity. In parallel, the exceptional change in the carbonate cycle in oceanic environment also suggests that this peak in biosphere productivity may have contributed to maintain the atmospheric CO$_2$ concentration at relatively low level at the beginning of MIS 11 while strong carbonate production and export, would have the opposite tendency of increasing it. The exceptional strong biosphere productivity during this period can be linked to low eccentricity around 400 ka, probably influencing productivity in both the terrestrial and oceanic ecosystems, due to longer summer insolation duration. The slow sea level rise during Termination V followed by an exceptional warming at high latitude of the Northern hemisphere over mid-MIS 11 may also have favoured increasing terrestrial biosphere productivity through vegetation in emerged low-latitude platform during Termination V and vegetation on high latitude terrestrial regions usually characterised by frozen soils during MIS 11.

Finally, this important result on the global productivity calls for future studies, especially at the regional scale through measurements of productivity signals in oceanic and terrestrial reservoirs over Termination V. This study should also be extended to refine recent terminations (e.g. Termination III where $\Delta^{17}\text{O}$ of O$_2$ resolution is very low) and bring attention to older ones. The
comparisons of high-resolution records of CO₂ and Δ¹⁷O of O₂ should permits to study the coupled dynamic of biosphere productivity and carbon cycle associated with large climate changes, under different orbital configurations.

## Method

**Material.** The ice samples used for this study come from bag samples allocated to gas consortium (A cut) of the EPICA Dome C 2 ice core. They were stored after drilling (2002-2003) in France (Le Fontanil) at -20 °C before being transported to LSCE, Gif sur Yvette in December 2017 for analysis in first semester of 2018. In total, 53 ice samples were analysed but 3 samples were lost due to an excess of pressure in the line.

**Extraction of air trapped in bubble from ice.** Extraction of air trapped in bubbles of ice was performed using a semi-automatic extraction and separation line. A 5-mm-thick section of ice was removed from each sample prior to measurement to prevent any contamination from ambient air or drilling fluid. The sample is then introduced in a frozen flask (-20 °C) attached to the extraction line together with 5cc flasks of outside air (standard for isotopic composition of oxygen). For each daily sequence, two flasks of exterior air and three duplicated ice samples (i.e. 6 ice samples) were analysed. After pumping the line 40 min with samples kept at -20 °C, the ice samples are melted to extract the air trapped in the bubble from the melted ice. The air then goes through H₂O and CO₂ traps immersed in ethanol cooled at -90 °C and in liquid nitrogen (-196 °C) respectively. The gas sample is transferred in a molecular sieve at -196 °C, then in a chromatographic column (1 m*2 mm) to separate O₂ and Ar from N₂ and finally goes through another molecular sieve to purify O₂ and Ar from remaining He. The extracted air is then trapped in stainless-steel tubes in an 8 port collection manifold immersed in liquid helium (4.13 K). The method follow description from Barkan and Luz²⁷ except that the chromatographic column is longer (1 m instead of 20 cm) and kept at higher temperature (0 °C instead of -90 °C).

**Mass spectrometer measurements.** Measurements of ¹⁸O/¹⁶O, ¹⁷O/¹⁶O and ³²O₂/⁴⁰Ar (δO₂/Ar) were performed using a multi-collector mass spectrometer Thermo Scientific MAT253 run in dual inlet mode. For each extracted air sample, 3 runs of 24 dual inlet measurements were performed against a laboratory standard gas obtained by mixing commercial O₂ and commercial Ar in atmospheric concentration. The mean of the three runs were calculated for each sample and the standard mean deviation is calculated as the pooled standard deviation over all duplicate samples hence integrating the variability associated with extraction, separation and mass spectrometry analysis.

The resulting pooled standard deviations before the corrections are 0.05‰, 0.02‰, 6.4‰ and 6 per meg for δ¹⁸O_atm, δ¹⁷O_atm, δO₂/Ar measurements and Δ¹⁷O of O₂ calculation, respectively.

**Atmospheric air calibration.** Every day, δ¹⁸O, δ¹⁷O and δO₂/Ar of atmospheric air are measured and are then used to calibrate our measurements following:

$$ \delta^{18}\text{O}_{\text{ext air corr}} = \left[ \frac{\left( \frac{\delta^{18}\text{O}_{\text{sample}}}{1000} \right) + 1}{\left( \delta^{18}\text{O}_{\text{ext air}}/1000 \right) + 1} - 1 \right] * 1000, \quad (1) $$

$$ \delta^{17}\text{O}_{\text{ext air corr}} = \left[ \frac{\left( \frac{\delta^{17}\text{O}_{\text{sample}}}{1000} \right) + 1}{\left( \delta^{17}\text{O}_{\text{ext air}}/1000 \right) + 1} - 1 \right] * 1000. \quad (2) $$

The δ¹⁸O_atext air and δ¹⁷O_atext air were constant during the two measurement periods so that we used the average values over the two corresponding periods to correct the raw data. The daily correction was the same every day during each period.

**Correction due to fractionation in the firn column.** Several corrections were applied on the measured data due to fractionation in the firn. We describe here corrections due to gravitational and thermal fractionation as well as pore close-off effect.

Gravitational fractionation operates in firn due to Earth gravity field. δ¹⁸O and δ¹⁷O were obtained by corrections for this diffusive process using δ¹⁵N in neighbouring samples⁵⁸. The correction applied depends on the difference of mass between the two isotopes considered so that:

$$ \delta^{18}\text{O}_{\text{gravitational corr}} = \delta^{18}\text{O}_{\text{measured}} - 2 * \delta^{15}\text{N}, \quad (3) $$

$$ \delta^{17}\text{O}_{\text{gravitational corr}} = \delta^{17}\text{O}_{\text{measured}} - 1 * \delta^{15}\text{N}. \quad (4) $$

Diffusive processes operating in the firn column due to changes in temperature or because of the Earth gravity must also be taken into account. These processes lead to isotopic fractionation of O₂ that was taken into account for δ¹⁸O_atm reconstruction in the NGRIP ice core over abrupt temperature changes of the last glacial period⁵⁹. To check the effect of thermal fractionation on Δ¹⁷O of O₂, we performed measurements of Δ¹⁷O of O₂ over the top 17 m of the EastGRIP firn

were a strong seasonal gradient is present. However, the resulting Δ¹⁷O of O₂ was not showing any significant deviation from the atmospheric value hence suggesting that thermal fractionation does not modify the Δ¹⁷O of O₂ of the atmosphere in the firn column (Supplementary Table 3). Moreover, in Antarctica, surface temperature variations are much lower than in Greenland during deglaciations or climatic variability of the last glacial period so that thermal fractionation is not expected to have a significant effect on the isotopic composition of trapped oxygen.

Pore close-off at the bottom of the firn has been shown to affect δO₂/N₂, δAr/N₂ with potential effects on δ¹⁵N and δ⁴⁰Ar in certain cases⁶⁰. We checked this possible effect on Δ¹⁷O of O₂ by comparing Δ¹⁷O of O₂ in bubbly ice at the top of the NEEM ice core. After correction of gravitational effect, we found a systematic enrichment of 13 per meg which could potentially bias the reconstruction of atmospheric Δ¹⁷O of O₂ from Δ¹⁷O of O₂ in trapped air.

**Gas loss correction.** During storage, O₂ in ice samples is subject to gas loss fractionation due to diffusion processes⁶¹ and O₂/N₂ ratio is always lower by several % in ice samples stored several years at -20 °C than ice samples stored at -50 °C²⁶. Such gas loss effect is also associated with isotopic fractionation of oxygen, δ¹⁸O_atm trapped in the ice being higher when δO₂/N₂ decreases with a slope for the relationship of -0.01 (δ¹⁸O_atm vs δO₂/N₂)¹⁷,¹⁸,⁶³,⁶⁴. We thus expect that δ¹⁷O can also be affected by this gas loss process and that it may create an anomaly of Δ¹⁷O of O₂. To check this effect, measurements of δ¹⁷O, δ¹⁸O_atm, Δ¹⁷O of O₂ and δO₂/Ar have been performed on three samples of GRIP ice core (clathrate ice stored during more than 20 years at -20 °C). Each of the three ice samples have been cut in order to analyse the interior and the exterior of the sample separately. δO₂/N₂ measurements could not be performed on exactly the same samples so that we used the δO₂/Ar measurements to estimate the amount of oxygen loss. Argon is also known to be affected by gas loss¹⁷,⁶¹ but in a smaller extend than oxygen⁶⁵ so that the decrease of δO₂/Ar is still a good indication of larger gas loss.

Data show a systematic lower Δ¹⁷O of O₂ value in the exterior sample compared to the interior sample, paralleling the lower δO₂/Ar value in the exterior sample compared to the interior sample as expected by gas loss (Supplementary Table 1). This systematic relationship and the Δ¹⁷O of O₂ and δO₂/Ar values can be used to propose a correction for the Δ¹⁷O of O₂ that takes into account the gas loss effect:

$$ \Delta^{17}\text{O}_{\text{gas loss corr}} = \Delta^{17}\text{O}_{\text{sample}} - 0.3945 * \left[ \left( \delta\text{O}_2/\text{Ar} \right)_{\text{sample}} - \left( \delta\text{O}_2/\text{Ar} \right)_{\text{std}} \right]. \quad (5) $$

**Comparison of EDC Δ¹⁷O of O₂ with previous Δ¹⁷O of O₂ record.** Previous Δ¹⁷O of O₂ measurements only took into account correction linked with gravitational fractionation. As a consequence, we checked the consistency of previous records performed on the GISP2 and Vostok ice core with our new EDC data by measuring again Δ¹⁷O of O₂ over Termination II using EDC samples and integrating the aforementioned correction (Fig. 1).

**Flux of oxygen associated with biosphere productivity.** We follow here the calculation of the flux of oxygen associated with biosphere productivity described in Landais et al.³³.

Since the atmospheric Δ¹⁷O of O₂ (Δ¹⁷O_atm) is influenced by the exchanges with biosphere and stratosphere, it is possible to write the following equation in a 3-box system at equilibrium (biosphere—bio-, atmosphere—atm-, stratosphere—strat-):

$$ F_{\text{bio}} * (\Delta^{17}\text{O}_{\text{bio}} - \Delta^{17}\text{O}_{\text{atm}}) = F_{\text{start}} * (\Delta^{17}\text{O}_{\text{start}} - \Delta^{17}\text{O}_{\text{atm}}), \quad (6) $$

where Fbio is the flux of oxygen exchanged by photosynthesis and respiration (both fluxes being considered at equilibrium) and Fstrat is the flux of oxygen exchanged between the lower atmosphere (or troposphere) and the stratosphere (fluxes in and out of the stratosphere are assumed equal).

To reconstruct past biospheric fluxes between terrestrial/oceanic biosphere and atmosphere based on Δ¹⁷O of O₂, it is necessary to know the evolution of the stratospheric flux as well as of the Δ¹⁷O_strat. Luz et al.²⁰ and Blunier et al.²² showed that a good assumption is to consider that the production rate of depleted O₂ in the stratosphere can be related to the atmospheric CO₂ concentration. It is then possible to express the evolution of biosphere oxygen flux in the past compared to pre-industrial from the following equation:

$$ \frac{F_{\text{bio,t}}}{F_{\text{bio,pre-industrial}}} = \frac{(\text{CO}_2)_t}{(\text{CO}_2)_{\text{pre-industrial}}} * \frac{\Delta^{17}\text{O}_{\text{bio,pre-industrial}}}{\Delta^{17}\text{O}_{\text{bio,t}} - \Delta^{17}\text{O}_{\text{atm,t}}}, \quad (7) $$

where Fbio,t is the biosphere oxygen flux at a given time, Fbio,pre-industrial is the present time biosphere oxygen flux, (CO₂,t) and (CO₂)pre-industrial correspond to the atmospheric CO₂ concentration at a given time and at the pre-industrial period respectively, Δ¹⁷O_bio,pre-industrial and Δ¹⁷O_bio,t are the values of Δ¹⁷O of O₂ in an atmosphere that would only be influenced by exchanges with the biosphere at the present time and at a given time respectively. Δ¹⁷O_atm,t is the Δ¹⁷O of O₂ of the atmosphere at a given time measured in the air trapped in ice core. We refer to pre-industrial period because of the long residence time of oxygen in the atmosphere
(>1000 years), Δ¹⁷O of O₂ is not exhibiting any significant variation over the last centuries.

Δ¹⁷O<sub>bio</sub> can be calculated from the fractionation coefficients associated with the different processes leading to oxygen fluxes in the biosphere (mainly photosynthesis, dark respiration and photosynthesis). Detailed calculations of Δ¹⁷O<sub>bio</sub> at the LGM and pre-industrial periods were obtained in Landais et al.³³ taking into account uncertainties in the determination of the fractionation coefficients³⁴ as well as on the relative fluxes of oxygen (Supplementary Table 2).

From the calculation of Δ¹⁷O<sub>bio</sub>, pre-industrial and Δ¹⁷O<sub>bio</sub>, LGM (Supplementary Table 2), Δ¹⁷O<sub>biot</sub> is calculated through a scaling on the variations of CO₂ concentration between pre-industrial period (280 ppmv) and the LGM (190 ppmv)⁶⁶ such as

$$ \Delta^{17}\text{O}_{\text{bio,t}} = \Delta^{17}\text{O}_{\text{bio,pre-industrial}} + (\Delta^{17}\text{O}_{\text{bio,LGM}} - \Delta^{17}\text{O}_{\text{bio,pre-industrial}}) \cdot \left( \frac{280 - (\text{CO}_2)_t}{90} \right). \quad (8) $$

**Details on uncertainty in photosynthesis fractionation.** We performed sensitivity tests to compare biosphere productivity reconstructions without any fractionation and reconstruction with the different fractionation effects during marine photosynthesis as observed in Eisenstadt et al.³⁷. The largest change on the reconstructed biosphere productivity is obtained using the observed photosynthesis fractionation effect associated with *Emiliania huxleyi* (δ¹⁸O = 5.81‰ and slope between ln(δ¹⁷O + 1) and ln(δ¹⁸O + 1) of 0.5253). The reconstructed global biosphere productivity is lower with such fractionation effect: the MIS 11 productivity level is decreased by 3% with such fractionation effect compared to the reconstruction with no fractionation at photosynthesis (i.e. the average situation, see Supplementary Table 2).

**Details on uncertainty in respiration fractionation.** When taking into account the maximum effect, i.e. a decrease of 0.005 for the slope of the relationship between ln(δ¹⁷O + 1) and ln(δ¹⁸O + 1) during respiration, we end up with a decrease of Δ¹⁷O<sub>bio</sub> by 65 ppm. The resulting biosphere productivity reconstruction is about 16% higher during MIS 11 than the average situation (Supplementary Table 2) leading to an anomalously high O₂ flux associated with gross primary productivity during Termination V.

**Validation of the biosphere productivity reconstruction.** The three sensitivity tests displayed (uncertainty on photosynthesis fractionation, on respiration fractionation and on fractionation within the water cycle) suggest that the biosphere productivity during MIS 11 is significantly higher than for our current interglacial, leading to the envelop displayed on Supplementary Fig. 2. Actually, explaining the Δ¹⁷O of O₂ anomaly without invoking a significant increase of the global productivity at the beginning of MIS 11 requires huge change of the Δ¹⁷O of O₂ produced by terrestrial biosphere (through changes in fractionation factor or changes in water cycle, hypothesis 1) or of the ratio of oceanic to terrestrial biosphere productivity (hypothesis 2). These two possibilities are rather unrealistic as detailed below.

In hypothesis 1, an increase of 30 ppm for Δ¹⁷O of O₂ produced by terrestrial biosphere is needed. This can be obtained either by a huge change in the fractionation factors between our current interglacial and MIS 11 (which is not realistic since these are based on physical processes which do not vary with time) or by a change of water cycle during MIS 11 with respect to present-day value. To reach this 30 ppm increase through modification of the water cycle, one option is to increase the ¹⁷O-excess of meteoric water by 30 ppm during MIS 11 compared to our current interglacial. We do not have ¹⁷O-excess values for MIS 11 yet, but ice core ¹⁷O-excess values obtained for the last interglacial are very similar to values obtained for the present interglacial³⁹. Moreover, increasing ¹⁷O-excess by 30 ppm would require unrealistic decrease of relative humidity at evaporation (by 30%, Barkan and Luz⁶⁷) during MIS 11 compared to our present interglacial. Finally, an increase of the mean δ¹⁸O and δ¹⁷O of meteoric water used by the plant along the meteoric water line (i.e. without global change in the global ¹⁷O-excess) would as well be a solution (see Fig. 4 of Landais et al.³³) but to reach the expected increase of 30 ppm in Δ¹⁷O of O₂ produced by the terrestrial biosphere, it would require an increase of global mean δ¹⁸O of meteoric water by ~3 ‰ with respect to present-day that will be transmitted to the global δ¹⁸O<sub>atm</sub>. The δ¹⁸O<sub>atm</sub> during MIS 11 is between 0 and 0.5 ‰ higher than during MIS 1 so that it does not support this hypothesis.

Hypothesis 2 requires an increase of the oceanic vs terrestrial productivity by a factor of 2 during MIS 11 compared to present-day value. This does not go along available observations suggesting an increase in the ratio of terrestrial to oceanic productivity during MIS 11 (see main text) and is still well above the uncertainty for the variation of the ratio of oceanic to terrestrial productivity during our interglacial period (14%, Supplementary Table 2).

An additional validation of our reconstruction of exceptional productivity during the beginning of MIS 11 comes from the comparison of our reconstruction with the oxygen biosphere productivity variations obtained by Blunier et al.²² over the last four terminations using a different vegetation model and different formulation of the link between atmospheric CO₂ and flux of Δ¹⁷O of O₂ anomaly

from the stratosphere. Indeed, in our approach, we directly dealt with Δ¹⁷O of O₂ anomaly in the box model while it may be more appropriate to deal with the δ¹⁷O and δ¹⁸O of O₂ values as done in Blunier et al.²² and suggested by Prokopenko et al.⁶⁸. We thus compared our reconstruction with the one performed with the model of Blunier et al.²² (Supplementary Fig. 2). The smaller difference in the oxygen flux associated with oxygen productivity between glacial and interglacial periods in Blunier et al.²² mainly comes from the different formulations of the dependency of the stratospheric anomaly Δ¹⁷O of O₂ to CO₂ atmospheric concentration. Another difference comes from the variations in the triple isotopic composition of oxygen in water between glacial and interglacial periods not taken into account in Landais et al.³³. We took this possible variation into account in our approach through the third sensitivity test explained above. We ran here the model of Blunier et al.²² over the new Termination V data. The comparison between the two methods of reconstruction (Landais et al.³³ vs Blunier et al.²²) is significant. The biosphere productivity increase at the end of Termination V is less important using the Blunier et al.²² approach than using the Landais et al.³³: 10% larger than for other interglacial periods following Blunier et al.²² instead of 20% on average following Landais et al.³³. Still, both approaches show that the biosphere productivity increase over Termination V is exceptional compared to the younger terminations.

## Data availability

The corrected data of Δ¹⁷O of O₂ are available in Pangaea database https://doi.pangaea.de/10.1594/PANGAEA.914609. The data are also available from the authors.

Received: 28 August 2019; Accepted: 20 March 2020; Published online: 30 April 2020

## References

1. Lüthi, D. et al. High-resolution carbon dioxide concentration record 650,000–800,000 years before present. *Nature* **453**, 379–382 (2008).

2. Menviel, L. et al. Southern Hemisphere westerlies as a driver of the early deglacial atmospheric CO₂ rise. *Nat. Commun.* **9**, 2503 (2018).

3. Martínez-Garcia, A. et al. Links between iron supply, marine productivity, sea surface temperature, and CO₂ over the last 1.1 Ma. *Paleoceanography* **24**, PA1207 (2009).

4. Stephens, B. B. & Keeling, R. F. The influence of Antarctic sea ice on glacial-interglacial CO₂ variations. *Nature* **404**, 171–175 (2000).

5. Volk, T. & Hoffert, M. I. in *Geophysical Monograph Series* (eds. Sundquist, E. T. & Broecker, W. S.) 99–110 (American Geophysical Union, 1985).

6. Tagliabue, A. et al. Quantifying the roles of ocean circulation and biogeochemistry in governing ocean carbon-13 and atmospheric carbon dioxide at the last glacial maximum. *Clim. Past* **5**, 695–706 (2009).

7. Toggweiler, J. R. Variation of atmospheric CO₂ by ventilation of the ocean's deepest water. *Paleoceanography* **14**, 571–588 (1999).

8. Ciais, P. et al. Large inert carbon pool in the terrestrial biosphere during the Last Glacial Maximum. *Nature Geoscience* **5**, 74–79 (2012).

9. Prentice, I. C., Harrison, S. P. & Bartlein, P. J. Global vegetation and terrestrial carbon cycle changes after the last ice age. *New Phytologist* **189**, 988–998 (2011).

10. Peterson, C. D., Lisiecki, L. E. & Stern, J. V. Deglacial whole-ocean δ¹³C change estimated from 480 benthic foraminiferal records. *Paleoceanography* **29**, 549–563 (2014).

11. Eggleston, S., Schmitt, J., Bereiter, B., Schneider, R. & Fischer, H. Evolution of the stable carbon isotope composition of atmospheric CO₂ over the last glacial cycle. *Paleoceanography* **31**, 434–452 (2016).

12. Schmitt, J. et al. Carbon Isotope Constraints on the Deglacial CO₂ Rise from Ice Cores. *Science* **336**, 711–714 (2012).

13. Duchamp-Alphonse, S. et al. Enhanced ocean-atmosphere carbon partitioning via the carbonate counter pump during the last deglacial. *Nat. Commun.* **9**, 2396 (2018).

14. Diekmann, B. & Kuhn, G. Sedimentary record of the mid-Pleistocene climate transition in the southeastern South Atlantic (ODP Site 1090). *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **182**, 241–258 (2002).

15. Tzedakis, P. C., Hooghiemstra, H. & Pálike, H. The last 1.35 million years at Tenaghi Philippon: revised chronostratigraphy and long-term vegetation trends. *Quaternary Sci. Rev.* **25**, 3416–3430 (2006).

16. Melles, M. et al. 2.8 Million Years of Arctic Climate Change from Lake El ggytgyn, NE Russia. *Science* **337**, 315–320 (2012).

17. Severinghaus, J. P., Beaudette, R., Headly, M. A., Taylor, K. & Brook, E. J. Oxygen-18 of O₂ records the impact of abrupt climate change on the terrestrial biosphere. *Science* **324**, 1431–1434 (2009).

18. Extier, T. et al. On the use of δ¹⁸Oatm for ice core dating. *Quaternary Sci. Rev.* **185**, 244–257 (2018).
19. Luz, B. & Barkan, E. The isotopic composition of atmospheric oxygen. *Global Biogeochem. Cycles* **25**, GB3001 (2011).

20. Luz, B., Barkan, E., Bender, M. L., Thiemens, M. H. & Boering, K. A. Triple-isotope composition of atmospheric oxygen as a tracer of biosphere productivity. *Nature* **400**, 547–550 (1999).

21. Blunier, T., Barnett, B., Bender, M. L. & Hendricks, M. B. Biological oxygen productivity during the last 60,000 years from triple oxygen isotope measurements. *Global Biogeochem. Cycles* **16**, https://doi.org/10.1029/2001GB001460 (2002).

22. Blunier, T., Bender, M. L., Barnett, B. & von Fischer, J. C. Planetary fertility during the past 400 ka based on the triple isotope composition of O₂ in trapped gases from the Vostok ice core. *Clim. Past* **8**, 1509–1526 (2012).

23. Hoffmann, G. et al. A model of the Earth's Dole effect. *Global Biogeochem. Cycles* **18**, GB1008 (2004).

24. Bazin, L. et al. An optimized multi-proxy, multi-site Antarctic ice and gas orbital chronology (AICC2012): 120-800 ka. *Clim. Past* **9**, 1715–1731 (2013).

25. Paillard, D. The Plio-Pleistocene climatic evolution as a consequence of orbital forcing on the carbon cycle. *Clim. Past* **13**, 1259–1267 (2017).

26. de Vernal, A. & Hillaire-Marcel, C. Natural variability of greenland climate, vegetation, and ice volume during the past million years. *Science* **320**, 1622–1625 (2008).

27. Cabioch, G. et al. The chronology and structure of the western New Caledonian barrier reef tracts. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **268**, 91–105 (2008).

28. Barker, S. et al. Globally increased pelagic carbonate production during the Mid-Brunhes dissolution interval and the CO₂ paradox of MIS 11. *Quaternary Sci. Rev.* **25**, 3278–3293 (2006).

29. Husson, L. et al. Reef carbonate productivity during quaternary sea level oscillations. *Geochem. Geophys. Geosyst.* **19**, 1148–1164 (2018).

30. Saavedra-Pellitero, M., Baumann, K.-H., Ullermann, J. & Lamy, F. Marine Isotope Stage 11 in the Pacific sector of the Southern Ocean; a coccolithophore perspective. *Quaternary Sci. Rev.* **158**, 1–14 (2017).

31. Yau, A. M., Bender, M. L., Marchant, D. R. & Mackay, S. L. Geochemical analyses of air from an ancient debris-covered glacier, Antarctica. *Quaternary Geochronol.* **28**, 29–39 (2015).

32. Schilt, A. et al. Glacial–interglacial and millennial-scale variations in the atmospheric nitrous oxide concentration during the last 800,000 years. *Quaternary Sci. Rev.* **29**, 182–192 (2010).

33. Landais, A., Lathiere, J., Barkan, E. & Luz, B. Reconsidering the change in global biosphere productivity between the Last Glacial Maximum and present day from the triple oxygen isotopic composition of air trapped in ice cores. *Global Biogeochem. Cycles* **21**, GB1025 (2007).

34. Helman, Y., Barkan, E., Eisenstadt, D., Luz, B. & Kaplan, A. Fractionation of the three stable oxygen isotopes by oxygen-producing and oxygen-consuming reactions in photosynthetic organisms. *Plant Physiol.* **138**, 2292 (2005).

35. Angert, A. et al. The contribution of respiration in tree stems to the Dole Effect. *Biogeosciences* **9**, 4037–4044 (2012).

36. Stolper, D. A., Fischer, W. W. & Bender, M. L. Effects of temperature and carbon source on the isotopic fractionations associated with O₂ respiration for ¹⁷O/¹⁶O and ¹⁸O/¹⁶O ratios in E. coli. *Geochim. Cosmochim. Acta* **240**, 152–172 (2018).

37. Eisenstadt, D., Barkan, E., Luz, B. & Kaplan, A. Enrichment of oxygen heavy isotopes during photosynthesis in phytoplankton. *Photosynthesis Res.* **103**, 97–103 (2010).

38. Guy, R. D., Fogel, M. L. & Berry, J. A. Photosynthetic fractionation of the stable isotopes of oxygen and carbon. *Plant Physiol.* **101**, 37–47 (1993).

39. Landais, A., Barkan, E. & Luz, B. Record of δ¹⁸O and ¹⁷O-excess in ice from Vostok Antarctica during the last 150,000 years. *Geophys. Res. Lett.* **35**, https://doi.org/10.1029/2007GL032096 (2008).

40. Winkler, R. et al. Deglaciation records of ¹⁷O-excess in East Antarctica: reliable reconstruction of oceanic normalized relative humidity from coastal sites. *Clim. Past* **8**, 1–16 (2012).

41. Bender, M., Sowers, T. & Labeyrie, L. The Dole Effect and its variations during the last 130,000 years as measured in the Vostok Ice Core. *Global Biogeochem. Cycles* **8**, 363–376 (1994).

42. Reutenauer, C. et al. Quantifying molecular oxygen isotope variations during a Heinrich stadial. *Clim. Past* **11**, 1527–1551 (2015).

43. Landais, A. et al. What drives the millennial and orbital variations of δ¹⁸Oatm? *Quaternary Sci. Rev.* **29**, 235–246 (2010).

44. Cheng, H. et al. The Asian monsoon over the past 640,000 years and ice age terminations. *Nature* **534**, 640–646 (2016).

45. Gallagher, S. J. et al. East Asian Monsoon History and Paleoeceanography of the Japan Sea Over the Last 460,000 Years. *Paleoceanogr. Paleoclimatol.* https://doi.org/10.1029/2018PA003331 (2018).

46. Sadori, L. et al. Pollen-based paleoenvironmental and paleoclimatic change at Lake Ohrid (south-eastern Europe) during the past 500 ka. *Biogeosciences* **13**, 1423–1437 (2016).

47. Kousis, I. et al. Centennial-scale vegetation dynamics and climate variability in SE Europe during Marine Isotope Stage 11 based on a pollen record from Lake Ohrid. *Quaternary Sci. Rev.* **190**, 20–38 (2018).

48. Whittaker, R. *Communities and Ecosystems* 2nd edn (Macmillan, New York, 1975).

49. Prokopenko, A. A., Bezrukova, E. V., Khursevich, G. K., Solotchina, E. P. & Kuzmin, M. I. Climate in continental interior Asia during the longest interglacial of the past 500 000 years: the new MIS 11 records from Lake Baikal, SE Siberia. *Clim. Past Discussions* **5**, 1951–1988 (2009).

50. Spratt, R. M. & Lisiecki, L. E. A Late Pleistocene sea level stack. *Clim. Past* **12**, 1079–1092 (2016).

51. Salter, I. et al. Carbonate counter pump stimulated by natural iron fertilization in the Polar Frontal Zone. *Nat. Geosci.* **7**, 885–889 (2014).

52. Rickaby, R. E. M. et al. Coccolith chemistry reveals secular variations in the global ocean carbon cycle? *Earth Planet. Sci. Lett.* **253**, 83–95 (2007).

53. Flores, J.-A., Filippelli, G. M., Sierro, F. J. & Latamir, J. The “White Ocean” Hypothesis: A Late Pleistocene Southern Ocean governed by Coccolithophores and driven by Phosphorus. *Front. Microbiol.* **3**, https://doi.org/10.3389/fmicb.2012.00233 (2012).

54. Droxler, A. W., Alley, R. B., Howard, W. R., Poore, R. Z. & Burckle, L. H. Unique and exceptionally long interglacial marine isotope stage 11: Window into Earth warm future climate. *Geophys. Monogr. Ser.* **137**, 1–14 (2003).

55. Qin, B., Li, T., Xiong, Z., Algeo, T. J. & Jia, Q. Deep-water carbonate ion concentrations in the Western Tropical Pacific since the Mid-Pleistocene: a major perturbation during the Mid-Brunhes. *J. Geophys. Res.: Oceans* **123**, 6876–6892 (2018).

56. Uemura, R. et al. Asynchrony between Antarctic temperature and CO₂ associated with obliquity over the past 720,000 years. *Nat. Commun.* **9**, 961 (2018).

57. Barkan, E. & Luz, B. High-precision measurements of ¹⁷O/¹⁶O and ¹⁸O/¹⁶O of O₂ and O₂/Ar ratio in air. *Rapid Commun. Mass Spectrom.* **17**, 2809–2814 (2003).

58. Dreyfus, G. B. et al. Firn processes and δ15N: potential for a gas-phase climate proxy. *Quaternary Sci. Rev.* **29**, 28–42 (2010).

59. Landais, A. et al. The glacial inception as recorded in the NorthGRIP Greenland ice core: timing, structure and associated abrupt temperature changes. *Clim. Dyn.* **26**, 273–284 (2005).

60. Severinghaus, J. & Battle, M. Fractionation of gases in polar ice during bubble close-off: New constraints from firn air Ne, Kr and Xe observations. *Earth Planet. Sci. Lett.* **244**, 474–500 (2006).

61. Bender, M., Sowers, T. & Lipenkov, V. On the concentrations of O₂, N₂, and Ar in trapped gases from ice cores. *J. Geophys. Res.* **100**, 18651 (1995).

62. Ikeda-Fukazawa, T. et al. Effects of molecular diffusion on trapped gas composition in polar ice cores. *Earth Planet. Sci. Lett.* **229**, 183–192 (2005).

63. Landais, A. A tentative reconstruction of the last interglacial and glacial inception in Greenland based on new gas measurements in the Greenland Ice Core Project (GRIP) ice core. *J. Geophys. Res.* **108**, 4563 (2003).

64. Capron, E. L'air piégé dans les glaces polaires: Contraintes chronologiques et caractérisation de la variabilité climatique rapide. Université de Versailles-St Quantin en Yvelines (2010).

65. Kobashi, T. et al. Post-bubble close-off fractionation of gases in polar firn and ice cores: effects of accumulation rate on permeation through overloading pressure. *Atmos. Chem. Phys.* **15**, 13895–13914 (2015).

66. Barnola, J.-M., Raynaud, D., Korotkevich, Y. S. & Lorius, C. Vostok ice core provides 160,000-year record of atmospheric CO₂. *Nature* **329**, 408 (1987).

67. Barkan, E. & Luz, B. Diffusivity fractionations of H₂¹⁶O/H₂¹⁷O and H₂¹⁶O/H₂¹⁸O in air and their implications for isotope hydrology. *Rapid Commun. Mass Spectrom.* **21**, 2999–3005 (2007).

68. Prokopenko, M. G., Pauluis, O. M., Granger, J. & Yeung, L. Y. Exact evaluation of gross photosynthetic production from the oxygen triple-isotope composition of O₂: implications for the net-to-gross primary production ratios. *Geophys. Res. Lett.* **38**, L14603 (2011).

69. Siegenthaler, U. Stable carbon cycle-climate relationship during the Late Pleistocene. *Science* **310**, 1313–1317 (2005).

70. Veres, D. et al. The Antarctic ice core chronology (AICC2012): an optimized multi-parameter and multi-site dating approach for the last 120 thousand years. *Clim. Past* **9**, 1733–1748 (2013).

71. Laskar, J. et al. A long-term numerical solution for the insolation quantities of the Earth. *Astron. Astrophys.* **428**, 261–285 (2004).

72. Petit, J. R. et al. Climate and atmospheric history of the past 420,000 years from the Vostok ice core. *Antarctica* **399**, 8 (1999).

73. Dreyfus, G. B. et al. Anomalous flow below 2700 m in the EPICA Dome C ice core detected using δ18O of atmospheric oxygen measurements. *Clim. Past* **3**, 341–353 (2007).

74. Dreyfus, G. B. et al. An ice core perspective on the age of the Matuyama–Brunhes boundary. *Earth Planet. Sci. Lett.* **274**, 151–156 (2008).

75. Landais, A. et al. Two-phase change in CO₂, Antarctic temperature and global climate during Termination II. *Nat. Geosci.* **6**, 1062–1065 (2013).
76. Bazin, L. et al. Phase relationships between orbital forcing and the composition of air trapped in Antarctic ice cores. *Clim. Past* 12, 729-748 (2016).

77. Lisiecki, L. E. & Raymo, M. E. A Pliocene-Pleistocene stack of 57 globally distributed benthic δ¹⁸O records. *Paleoceanography* 20, PA1003 (2005).

## Acknowledgements

The authors first acknowledge the three anonymous reviewers and the editor for their useful help in improving the manuscript. It is a pleasure to thank Jean Jouzel, Jochen Schmitt, Ji-Woong Yang and Michael Bender for useful discussions as well as Laurent Bopp, Masa Kageyama, Chiara Guarnieri and Jean-Yves Peterschmitt for their contribution on model reconstructions. We also thank Gregory Teste involved in cutting EDC samples. This work benefited from funding from the ANR program HUMI17 as well as the INSU-LEFE-IMAGO program BIOCOD. M.B. receives an IDEX-IDI PhD grant from Paris-Saclay. It is a contribution to the European Project for Ice Coring in Antarctica (EPICA). This work was supported by the French state aid managed by the ANR under the "Investissements d'avenir" programme with the reference ANR-11-IDEX-0004-17-EURE-0006.

## Author contributions

A.L., S.D.A. and M.B. designed this study. M.B, F.P., V.F, H.A. and L.S performed the Δ¹⁷O of O₂ measurements. M.B and A.L. corrected the data and calculated the flux of oxygen biosphere productivity with the help of T.B. M.B., A.L. and S.D.A. wrote the manuscript with the contribution of T.B and T.E.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** is available for this paper at https://doi.org/10.1038/s41467-020-15739-2.

**Correspondence** and requests for materials should be addressed to M.B.

**Peer review information** *Nature Communications* thanks the anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_6.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020