# A solution for constraining past marine Polar Amplification

A. Morley, E. de la Vega, M. Raitzsch, J. Bijma, U. Ninnemann, G. Foster, Thomas B. Chalk, J. Meilland, R. Cave, J. Büscher, et al.

► To cite this version:

A. Morley, E. de la Vega, M. Raitzsch, J. Bijma, U. Ninnemann, et al. A solution for constraining past marine Polar Amplification. Nature Communications, 2024, 15 (1), pp.9002. 10.1038/s41467-024-53424-w. hal-04745907

HAL Id: hal-04745907

https://hal.science/hal-04745907v1

Submitted on 24 Oct 2024

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
<figure><img src="image_2.png" /></figure>

# A solution for constraining past marine Polar Amplification

Received: 4 July 2024

Accepted: 10 October 2024

Published online: 18 October 2024

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

Check for updates

A. Morley Ⓔ<fnref n="1" />,²✉, E. de la Vega<fnref n="1" />, M. Raitzsch³, J. Bijma Ⓔ⁴, U. Ninnemann⁵,
G. L. Foster⁶, T. B. Chalk⁷, J. Meilland Ⓔ⁸, R. R. Cave⁹, J. V. Büscher⁹,¹⁰ &
M. Kucera Ⓔ⁸

Most climate proxies of sea surface temperatures suffer from severe limitations when applied to cold temperatures that characterize Arctic environments. These limitations prevent us from constraining uncertainties for some of the most sensitive climate tipping points that can trigger rapid and dramatic global climate change such as Arctic/Polar Amplification, the disruption of the Atlantic Meridional Overturning Circulation, sea ice loss, and permafrost melting. Here, we present an approach to reconstructing sea surface temperatures globally using paired Mg/Ca - δ¹⁸O_c recorded in tests of the polar to subpolar planktonic foraminifera *Neogloboquadrina pachyderma*. We show that the fidelity of Mg/Ca-based paleoclimate reconstructions is compromised by variations in seawater carbonate chemistry which can be successfully quantified and isolated from paleotemperature reconstructions using a multiproxy approach. By applying the calibration to the last glacial maximum, we show that marine polar amplification has been underestimated by up to 3.0 ± 1.0 °C in model-based estimates.

One of the key drivers of Arctic climate and Arctic Amplification is the ocean<fnref n="1" />,² because of its role in transporting heat to high latitudes, which influences not only marine temperatures but also terrestrial climates and sea ice extent³,⁴. As a result, temperature changes in the high-latitude oceans have global implications, impacting planetary albedo, atmosphere-ocean carbon dioxide exchange, atmospheric and ocean circulation, deep-water formation, freshwater budgets, productivity, and biological diversity⁵⁻⁸. However, large uncertainties remain about the sensitivity of Arctic climate to natural⁹ and anthropogenic forcing¹⁰ owing to the limited observational and historical network in this area¹¹ and uncertainties linked to the strength of cloud feedback in model simulations¹². In the absence of a direct proxy for clouds, large-scale surface temperature patterns are useful analogues¹³ and critical for

constraining equilibrium climate sensitivity in climate models¹². Efforts to assess the magnitude of past polar surface temperature patterns using climate proxies are thus crucial to further our understanding of how the Arctic will respond to, and influence continued global warming⁴,¹⁴. Ice-core records provide accurate and high-resolution surface air temperatures over high-altitude central Greenland for the past 120 thousand years (ka), however, the reconstruction of sea surface temperatures (SSTs) from high-latitude oceans has been challenging. This is because planktonic foraminifera species diversity collapses to one dominant species (*N. pachyderma*) at cold SSTs (< 4°), leading to large uncertainties associated with reconstructed SSTs based on assemblage counts¹⁵,¹⁶. Other proxies are affected by their own set of limitations in high-latitude oceans, for example, dinoflagellate cyst

<footnote n="1">¹University of Galway, School of Geography, Archaeology, and Irish Studies, and Ryan Institute, Galway, Ireland. ²ICRAG – Irish Centre for Research in Applied Geosciences, Belfield, Dublin 4, Ireland. ³Dettmer Group GmbH & Co. KG., Bremen, Germany. ⁴Alfred-Wegener-Institut, Helmholtz-Zentrum für Polar- und Meeresforschung, Bremerhaven, Germany. ⁵University of Bergen, Department of Earth Science and Bjerknes Centre for Climate Research, Bergen, Norway. ⁶University of Southampton, School of Ocean and Earth Science, National Oceanography Centre Southampton, Southampton, UK. ⁷Centre Européen de Recherche et d'enseignement des géosciences de l'environnement (CEREGE), Aix-en-Provence, France. ⁸MARUM – Center for Marine Environmental Sciences, University of Bremen, Bremen, Germany. ⁹University of Galway, School of Natural Sciences, Galway, Ireland. ¹⁰Ulster University, School of Geography and Environmental Sciences, Coleraine, UK. ✉e-mail: audrey.morley@universityofgalway.ie</footnote>
assemblages at high latitudes are controlled by multiple variables, and it is often difficult to ascertain which is responsible for the observed assemblage change¹⁷. In addition, standard geochemical proxies such as oxygen isotopes (δ¹⁸O), Mg/Ca-thermometry and Alkenone-based SST reconstructions cannot accurately constrain cold temperatures < 4 °C because secondary controls besides temperature influence all these proxies at low temperatures¹⁸⁻²⁰.

The Mg/Ca ratio of the CaCO₃ tests of foraminifera is an established and powerful tool for the reconstruction of past SSTs²¹. It is based on the principle that the partition of magnesium in calcite is endothermic and, therefore, leads to higher Mg/Ca values at higher temperatures²². In fact, the first Mg/Ca-temperature calibration, which ultimately led to the ‘acceptance’ of this temperature proxy within the scientific community, was based on measurements of *N. pachyderma*²³. However, this seminal work also excluded samples collected from “cold sites” annually covered by sea ice from which the Mg content appeared disconnected from temperature. This makes it difficult to use the otherwise powerful Mg/Ca proxy in the polar realm. Several attempts have been made since to calibrate the Mg/Ca - thermometer for *N. pachyderma*²⁰,²³⁻²⁷, however, to date, calibration efforts have not been able to resolve SSTs below 4 °C. The difficulty in solving cold SSTs is exacerbated by the low sensitivity of the exponential Mg/Ca temperature equation and because it appears that secondary controls have a dominant influence on Mg/Ca at low temperatures. These non-thermal controls lead to unrealistically high temperatures at high latitude sites causing reconstructed glacial maxima SSTs at polar and subpolar latitudes to be too warm by 6-7 °C, especially in the Greenland and Iceland Seas²⁰ relative to modelled SSTs (Supplementary Fig. S5). Importantly, these warm glacial maximum SSTs are at odds with paleoclimate simulations, causing significant data-model mismatch²⁸. To resolve this problem, calibration equations for *N. pachyderma* have been modified to ever more insensitive calibration equations²⁷,²⁹,³⁰, reducing the actual temperature estimates and the corresponding amplitude of glacial-interglacial SST change. Subsequently, these low-SST and low-amplitude records have been interpreted as ‘winter’ or subsurface (below 100–150 m) temperature signals. The latter interpretation is based on the introduction of a cold bias in the geochemical signature (both Mg/Ca and oxygen isotopes) recorded in the gametogenic crust *N. pachyderma* deposits at the end of its lifecycle that is frequently assumed to be formed in colder and deeper water³¹ (see also Methods). However, recent advances in our understanding of *N. pachyderma* ecology challenge the assumption of systematic ontogenic vertical migration to deeper or colder waters associated with reproduction or “crusting”³²,³³. Moreover, culture-grown *Neogloboquadrina incompta* and *N. pachyderma* recovered from sediment traps have lower crust Mg/Ca ratios than ontogenetic calcite formed at the same temperature, suggesting that a lower temperature is not responsible for the low Mg/Ca ratio of neogloboquadrinid crusts²⁵,³⁴.

The departure from the predicted Mg/Ca-temperature relationship at low temperatures has led to the hypothesis that, in addition to temperature, seawater carbonate ion concentration ([CO₃²⁻]) may exert a strong secondary control on Mg/Ca especially at lower concentrations often associated with cold SSTs²⁰,²⁴. Culture experiments where planktonic foraminifera grew under controlled conditions¹⁵⁻³⁷ and core-top calibration studies, for which modern hydrographic parameters are available, support this hypothesis, especially when [CO₃²⁻] values are < 200 µmol.kg⁻¹, which is common at polar and Arctic latitudes. To account for and correct for a possible [CO₃²⁻] influence on Mg/Ca-thermometry, the species-specific sensitivity must be fully resolved¹⁸,³⁹. Once the sensitivity is determined, mathematical correction schemes can isolate temperature from [CO₃²⁻] within the Mg/Ca ratios³⁹⁻⁴¹. The resulting tool would allow for accurate SST reconstructions at high latitudes and thereby provide the means to quantify polar amplification on multiple timescales under various boundary conditions. Furthermore, accurate estimates of past SSTs

are essential to constrain the oceanic carbonate system via boron isotope (δ¹¹B) geochemistry⁴², currently leaving a gap in our understanding of high-latitude ocean-atmosphere carbon fluxes.

Here, we define the Mg/Ca - [CO₃²⁻] sensitivity for *N. pachyderma* by measuring Mg/Ca ratios on specimens collected live (as determined by the presence of active cytoplasm) from plankton tows from the Labrador and Nordic Seas that can be paired directly with hydrographic conditions of the calcification environment (Supplementary Fig. S1 and Supplementary Tables 1–3). Plankton tows via multinet sampling allow for an accurate analysis of the hydrographic conditions of the calcification environment when assuming that the plankton tow depth is representative of the mean calcification depth of living specimens (for details on hydrographic data used, see Supplementary Tables 1–3). This approach allows us to circumvent the fundamental limitations and assumptions stemming from the traditional core-top approach to proxy establishment e.g., our analysis does not rely on assumptions that are required to link tests found in core-top sediments to hydrographic conditions during their formation. This is particularly important for Polar and Arctic Oceans, considering the scale of recent warming and low sedimentation/accumulation rates that result in the mixing of pre-industrial with modern samples in the top 0.5 cm of sediment.

## Results and discussion

Like previous core top calibrations for *N. pachyderma*, our plankton tow-based results show a departure from the expected Mg/Ca-temperature relationship (Fig. 1a). To estimate the influence of seawater [CO₃²⁻] on Mg/Ca values, we isolate the contribution of [CO₃²⁻] on measured Mg/Ca values using water-column data for temperature and carbonate-ion concentration, following a similar approach as described by Evans et al.³⁵, and Morley et al.³⁹. First, we normalised measured Mg/Ca values against predicted Mg/Ca values (Mg/Ca[pred]) to reveal by how much the measured Mg/Ca is an overprediction of the expected. Mg/Ca[pred] represents the Mg/Ca sensitive to temperature only and is determined as Mg/Ca[pred]=0.4exp(0.09T). This equation was derived for *N. incompta* calcifying at [CO₃²⁻] values <200 µmol.kg⁻¹ in Morley et al.³⁹. Given that *N. incompta* and *N. pachyderma* are genetical sister species⁴³ and share ecological preferences⁴⁴, we assume that the temperature sensitivity is shared between both species. We then derive a Mg/Ca-[CO₃²⁻] sensitivity for all samples included in this study by comparing normalised Mg/Ca[Norm] against measured [CO₃²⁻] (Fig. 1b). Mg/Ca[Corr] (Fig. 1c) then reveals the temperature-only component and is calculated by dividing Mg/Ca[meas] over the Mg/Ca-[CO₃²⁻] sensitivity derived in Fig. 1b. We further complement the tow dataset with new core tops and paired hydrographic profiles including temperature, salinity, Dissolved Inorganic Carbon (DIC) and Alkalinity (ALK) from the Nordic Seas, and also include the published Kozdon et al.²⁷ core top dataset in our analysis using depth estimate from foraminifera δ¹⁸Oc projected on equilibrium δ¹⁸Oc derived from hydrographic profiles (see Methods). The comparisons of plankton tows with core tops allow us to test if the tow-derived Mg/Ca-[CO₃²⁻] sensitivity can be applied to climate archives.

The Mg/Ca-[CO₃²⁻] sensitivity for *N. pachyderma* from plankton tows (Fig. 1b) is consistent with calibrations for the closely related planktonic foraminifera *N. incompta*³⁹ but is more pronounced than in temperate planktonic foraminifera³⁵. This can be explained by a non-linear response of [CO₃²⁻] on Mg incorporation at different temperatures. Specifically, the steeper response for subpolar and polar foraminifera would imply that the [CO₃²⁻] effect is greater at low temperatures indicating that the high sensitivity of Mg to low [CO₃²⁻] may be linked to the decreasing efficiency of Mg exclusion under conditions less favourable for calcification. Along this line of argument, Evans et al.⁴⁵, hypothesised that for low-Mg foraminifera like *N. pachyderma*, the mechanism by which Mg is removed during calcification competes energetically with the need to modify the carbonate chemistry of seawater for calcification at low [CO₃²⁻]. They propose
that it may be more effective for foraminifera to prioritise the removal of H⁺ over Mg²⁺ at low [CO₃²⁻] (this is also supported by ref. 46). Also, low metabolic/respiration rates hypothesised at low temperatures may further reduce the energy available (required) to exclude Mg²⁺ and modify the carbonate chemistry for calcification²⁷, leading to the negative relationship between Mg/Ca and [CO₃²⁻] at low values.

Our Mg/Ca data from plankton samples provide an opportunity to directly quantify the carbonate ion concentration effect on Mg/Ca measured in Arctic *N. pachyderma*. This is because the carbonate ion concentration for the calcification environment can be determined reliably, and our sampled environments cover the entire range of low carbonate ion conditions from 200 to less than 50 µmol.kgSW⁻¹ (Fig. 1b). The results reveal that the departure of the measured Mg/Ca[Meas] from the theoretical Mg/Ca[Pred], calculated assuming they reflect only ambient temperature, is exponentially correlated with [CO₃²⁻]. When deriving [CO₃²⁻] values at the apparent calcification depths for the new and previously published core top values included in this study (see Methods), we find that the Mg/Ca sensitivity to

[CO₃²⁻] for core top samples is within the 95% confidence intervals of the Mg/Ca sensitivity to [CO₃²⁻] derived for plankton tow samples (Fig. 1b). This is important, as it allows us to apply the Mg/Ca sensitivity to [CO₃²⁻] derived from plankton tows to *N. pachyderma* preserved at the seafloor and in climate archives.

Before we compared datasets, we took into account the difference in Mg/Ca between the crust and ontogenetic calcite for *N. pachyderma*²⁵,²⁶,⁴⁷. Recent findings show that the Mg/Ca ratios of the outer test of *N. pachyderma* are between 8.6 to 15% higher in living (uncrusted) specimens when compared to the outer test of dead specimens recovered from the same plankton tows⁴⁸. Similarly, the analysis of intra-test banding of Mg/Ca ratios in crusted *N. pachyderma* tests from North Atlantic sediment traps reveals that the addition of crust lowers overall Mg/Ca ratios of a specimen by 15%³⁴ when assuming an average crust volume of 50%. To compare Mg/Ca values from tows and core tops we therefore adjust uncorrected Mg/Ca values from plankton tow samples down by 15% to account for the difference in Mg between ontogenetic and crust calcite.

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

**Fig. 1** | Calibration. In all panels, closed black squares refer to *N. pachyderma* collected from plankton tow samples, open dark grey squares represent *N. pachyderma* from core tops²⁷, and open turquoise squares represent *N. incompta* from core tops³⁹. Inverted red triangles refer to *N. pachyderma* core top data from the South Atlantic and the Antarctic continental shelf⁶⁹, while yellow triangles refer to *N. pachyderma* sediment trap data from the Antarctic Peninsular²⁴. a shows measured Mg/Ca values. b Shows normalised Mg/Ca values against water-column [CO₃²⁻] for tows and δ¹⁸O_c derived [CO₃²⁻] values for core tops. The grey shading connotes the 95% confidence intervals associated with the regression. c shows Mg/Ca values corrected for [CO₃²⁻] seawater carbonate-ion concentration. The blue shading connotes the 95% confident intervals of fully propagated uncertainties for

the slope of the calibration relationship, while the grey shading connotes the fully predicted uncertainties which would apply to the full uncertainty on a single data point. d shows δ¹⁸O_c values for *N. pachyderma* vs water column [CO₃²⁻]. *N. pachyderma* from plankton tows are corrected for a + 2‰ offset to account for the difference between ontogenetic and crust calcite. The blue shading connotes the 95% confidence intervals for the regression combining both tow and core top data, while the yellow and red shading connotes the 95% confidence intervals for the tows and core top regressions, respectively. In addition, function plots are shown for other unicellular planktonic calcifying organisms⁵³,⁵⁴. Source data for all panels are provided as a Source Data file.
With plankton tows, new, and previously published North Atlantic core top samples combined (e.g., ref. 20,27,39), Mg/Ca[Corr] describes the following equation: $Mg/Ca_{[Corr]} = 0.372 \pm 0.014\exp^{(0.0936 \pm 0.0045T)}$, $n = 235$, $r^2 = 0.90$, which shares the same sensitivity (within error) of the predicted temperature calibration equation, showing that $[CO_3^{2-}]$ is the sole or at least the primary non-thermal control on the predicted Mg/Ca - temperature relationship (Fig. 1c). Our results also demonstrate that the influence of $[CO_3^{2-}]$ on Mg/Ca measured on *N. pachyderma* can be mathematically quantified and isolated from temperature as for other planktonic foraminifera species (e.g., ref. 40), when $[CO_3^{2-}]$ is known. Thus, for applications of the Mg/Ca paleothermometer on fossil material, an independent $[CO_3^{2-}]$ proxy is required. Here, we investigate the influence of $[CO_3^{2-}]$ on the fractionation of oxygen isotopes measured in *N. pachyderma* tests and evaluate its potential to serve as a proxy for the Mg/Ca correction scheme. $\delta^{18}O_c$ is a traditional palaeoceanographic tracer combining the influence of $\delta^{18}O_{sw}$ of seawater and temperature$^{49-51}$. However, in addition to temperature, a dependence of oxygen isotope fractionation on the carbonate chemistry of seawater has previously been demonstrated in planktonic foraminifera$^{52,53}$ and other unicellular planktonic calcifying organisms$^{54}$, suggesting the existence of a universal carbonate ion effect on stable oxygen isotope ratios measured in unicellular planktonic calcifying organisms.

To evaluate whether past seawater $[CO_3^{2-}]$ values can be derived from $\delta^{18}O_c$ in the past, we first tested if the $\delta^{18}O_c/[CO_3^{2-}]$ sensitivity observed in material from plankton tows and core tops (Fig. 1d) is compromised by interactions with $\delta^{18}O_{sw}$ and/or temperature in a generalised linear model including $\delta^{18}O_c$, $\delta^{18}O_{sw}$, $[CO_3^{2-}]$, and temperature. We found that $\delta^{18}O_{sw}$ does not significantly contribute to the total variance ($n=52, p > 0.01$, Supplementary Table 5) and therefore excluded it from considerations. Furthermore, we found that temperature only weakly contributes to the total variance ($p = 0.017$) observed in the model. Importantly, the interaction term between $[CO_3^{2-}]$ and temperature is not significant ($n = 52, p = 0.14$), suggesting that while temperature contributes to the total variance in $\delta^{18}O_c$, it does not influence or interact with the $\delta^{18}O_c/[CO_3^{2-}]$ relationship. Nevertheless, we evaluate the contribution of temperature and the potential of a non-linear effect of temperature on $\delta^{18}O_c$ further, using a regression tree (Supplementary Fig. S2). As expected, the regression tree shows that $[CO_3^{2-}]$ explains most of the variance and that temperature contributes to the variance in the $\delta^{18}O_c/[CO_3^{2-}]$ linear model only for values above 5 °C. When analysing the relationships between temperature, $[CO_3^{2-}]$, and $\delta^{18}O_c$ separately, we find that temperature and $[CO_3^{2-}]$ only covary above 5 °C or 150 µmol.kg⁻¹ but not at colder temperatures or lower concentrations (Supplementary Fig. S3). Similarly, the covariance between temperature and $\delta^{18}O_c$ breaks down below 4-5 °C (also observed in ref. 55), while the relationship between $\delta^{18}O_c$ and $[CO_3^{2-}]$ remains unchanged regardless of concentrations (Supplementary Fig. S3). In summary, while temperature contributes weakly to the variance in $\delta^{18}O_c$ above 5 °C, it does not influence the interaction between $\delta^{18}O_c$ and $[CO_3^{2-}]$. Furthermore, these results confirm that the $\delta^{18}O_c/[CO_3^{2-}]$ sensitivity does not rely on covariance with temperature, which varies over space and time.

Following, we posit that the sensitivity or slope of $\delta^{18}O_c/[CO_3^{2-}]$ is regulated by the calcite saturation product ($K_{sp}$), which in turn is determined by the conversion between (isotopically heavier) bicarbonate and (isotopically lighter) carbonate ions in the calcification environment, that is a function of pH, salinity, and to a lesser extent temperature$^{34}$. We observe that the slopes of $\delta^{18}O_c/[CO_3^{2-}]$ for *N. pachyderma* from tows and core tops are statistically the same at -0.0279 but offset by ~2‰ (Fig. 1d). The slopes in our dataset are also steeper than reported for *Orbulina universa* (-0.002) or *Globigerina bulloides* (e.g., -0.004)$^{53}$, but not outside of the slopes observed previously for other unicellular planktonic calcifying organisms$^{54}$ (Fig. 1d). The steeper slope for *N. pachyderma* suggests that the relative

contribution of the isotopically light $[CO_3^{2-}]$ to the calcifying fluid is low and that comparatively less bicarbonate is converted into carbonate ions in the calcification environment when compared to *O. universa* or *G. bulloides*. This hypothesis agrees well with the biomineralization model controlling the exclusion/uptake of Mg described earlier in *N. pachyderma*, where low cell energy due to unfavourable calcification environments and low metabolic/respiration rates at cold temperatures leads to reduced $[CO_3^{2-}]$ available for calcification$^{45,46,56}$. In addition, the release of CO₂ during respiration and the absence of photosynthetic activity from symbionts, further decrease the relative contribution of isotopically light $[CO_3^{2-}]$ in favour of the isotopically heavier bicarbonate and potentially further enhances the sensitivity of oxygen isotope fractionation to varying $[CO_3^{2-}]$ in low $[CO_3^{2-}]$ settings. The observed 2‰ offset in the intercepts between tows and core tops agrees well with the average 2.1 ‰ offset observed between ontogenetic calcite and a fully crusted *N. pachyderma*$^{57}$. A caveat of this offset is that specimens of a similar degree of encrustation must be analysed together to avoid the influence of varying degrees of encrustation on the $\delta^{18}O_c/[CO_3^{2-}]$ relationship. Accordingly, the larger scatter observed for the core top versus plankton tow dataset could potentially be explained by different degrees of crust formation$^{57}$.

To explore the application potential of the $\delta^{18}O_c/[CO_3^{2-}]$ proxy to fossil material, we compared $[CO_3^{2-}]$ derived from $\delta^{18}O_c$ with CARINA$^{38}$ derived $[CO_3^{2-}]$ values for the Kozdon core top dataset (see Methods). Both methods provide $[CO_3^{2-}]$ values that are within a root mean square error (RMSE) of 8.44 µmol.kg⁻¹ of each other (adjusted $r^2 = 0.95$ $n = 52, p < 0.001$). Residuals are symmetrically distributed around zero (Supplementary Fig. S4). When $\delta^{18}O_c$ derived $[CO_3^{2-}]$ values are used to estimate the Mg/Ca-$[CO_3^{2-}]$ sensitivity, the scatter is much lower (when compared to CARINA), resulting in more robust correlation statistics (adjusted $r^2 = 0.81$ vs $r^2 = 0.58$) and lower uncertainties. The reduced scatter is most likely due to the uncertainties linked to assigning modern water column hydrography from CARINA to pre-industrial/late Holocene core tops. By using $\delta^{18}O_c$ to derive $[CO_3^{2-}]$, these uncertainties are removed. To further test the validity of this approach, we also derived $[CO_3^{2-}]$ from $\delta^{18}O_c$ values published in the core top dataset (paired Mg/Ca - $\delta^{18}O_c$) by Meland et al.$^{20}$. and find that the Mg/Ca sensitivity to $[CO_3^{2-}]$ is statistically identical with the tows and the Kozdon datasets (Fig. 1b).

Finally, we combine all *N. pachyderma* datasets (e.g., this study, 20,27) with a previously published dataset for *N. incompta*$^{39}$. When all datasets are combined, we derive the following Mg/Ca - $[CO_3^{2-}]$ sensitivity: $Mg/Ca[CO_3^{2-}] = 3227.95 * [CO_3^{2-}]^{-1.53}$ ($n = 235, r^2 = 0.87, p < 0.01$). All uncertainties (including calibration uncertainties accounting for the covariance of error in calibration coefficients) are propagated through Monte Carlo simulation [e.g., 40]. Briefly, to account for dependent and independent uncertainties, 1000 potential calibration slopes are created from selecting random environmental values within the 2σ uncertainty ranges for individual parameters ($\delta^{18}O_c$, $[CO_3^{2-}]$, Mg/Ca, and Temperature, see Methods). The median of these slopes is taken as our central calibration estimate with the 2σ described by the 95th percentile of potential slopes. We also include the prediction interval for the Monte Carlo Ensemble (Fig. 1c, see Methods for full details). Typical uncertainty in calibrated temperatures, combining calibration uncertainty, and $[CO_3^{2-}]$ uncertainty, is ±1.0 °C (±1.2 °C at 0 °C and ±0.9 °C at 15 °C, 2σ) if $\delta^{18}O_c$ are used to correct for $[CO_3^{2-}]$.

To evaluate the suitability of the proposed equations to other regions and genotypes of *N. pachyderma*, we applied the correction scheme to the comprehensive core top dataset published in Vazquez-Riveiros et al.$^{39}$, that includes 38 core top samples from the subpolar South Atlantic to the Antarctic continental shelf and the sediment trap series from the Antarctic Peninsula published in Hendry et al.$^{24}$. These datasets are dominated by genotype IV, which is the only *N*.
*pachyderma* genotype adapted to overwintering in sea ice⁶⁰. Nineteen of the 38 core top sites, as well as the sediment trap series are situated on the Antarctic margin and are therefore subject to the possible influence of sea ice on Mg/Ca ratios²⁴,⁵⁹. When Mg/Ca are corrected for δ¹⁸O_c-derived [CO₃²⁻] values, the Mg/Ca-temperature correlation statistics for the South Atlantic core top dataset are much more robust (e.g., r² = 0.87 vs. r² = 0.30, n = 36, see also data plotted in Fig. 1a, c) and within the uncertainties of the North Atlantic derived dataset. Similarly, the sediment traps from Hendry et al.²⁴ agree well with the North Atlantic dataset (Fig. 1a, c). These applications, therefore, suggest that the proposed correction scheme can be applied to both hemispheres and that genotypes do not influence the geochemical climate signature recorded in *N. pachyderma* specimens as previously hypothesised⁶⁰.

The combination of the Mg/Ca correction scheme and δ¹⁸O_c derived [CO₃²⁻] values presents a unique opportunity to apply the correction to previously published Mg/Ca · δ¹⁸O_c datasets and estimate SSTs in polar environments. In this way, it allows us to reevaluate the magnitude and spatial expression of marine polar amplification to glacial-interglacial climate forcing as recorded by Mg/Ca and δ¹⁸O_c values. Using previously published Mg/Ca · δ¹⁸O_c datasets²⁹,³⁰,⁶¹, we show that uncorrected Mg/Ca-based SSTs systematically underestimate maximum glacial cooling by 2 to 3° ± 1.0 °C for four previous glacial maxima (Fig. 2 and Supplementary Table 6). Furthermore, the magnitude of warming from glacial to interglacial SSTs in the Labrador Sea at IODP Site U1305 is significantly higher based on Mg/Ca[corr] values for the past five glacial-interglacial cycles increasing from 2.3–4.2 °C for uncorrected Mg/Ca-derived values to 7.9-11.0 ± 1.0 °C

<figure><img src="image_5.png" /><figcaption>**Fig. 2** | **Paleo application of Mg/Ca correction scheme to four glacial-interglacial periods including Marine Isotope Stage 1 (MIS1), MIS5, MIS9 and MIS11**²⁹,³⁰,⁶¹. Datasets originate from Eirik Drift, Southeast of Greenland, from IODP site U1305, and sister site MD99-2227. Purple lines show three-point running averages of corrected Mg/Ca-based sea surface temperatures (SSTs), while grey lines show three-point running averages of uncorrected Mg/Ca-based SSTs</figcaption></figure>

reconstructions. The blue shading shows the 95% confidence interval associated with Mg/Ca[corr] SST reconstructions. In orange are Modern Analogue Technique (MAT)-based SSTs reconstructions as shown in Irvali et al.³⁰ for MIS 5 (panel b.), 9 (panel c.), and 11 (panel d.). Note that MAT-based SSTs are not available for the last glacial-interglacial transition (panel a.). Source data for all panels are provided as a Source Data file.
once corrected. This increase in gradient reflects both colder SSTs during glacial periods and warmer SSTs during interglacials, high-lighting the importance of applying the correction to both glacial and interglacial periods. Assemblage-based (e.g., modern analogue tech-nique; MAT⁶²) glacial-interglacial SST gradients for the same site range between 7–9 °C. MAT-based glacial-interglacial SST gradients are, on average, 1.5 °C smaller than our corrected Mg/Ca-based estimates, which is likely the result of low species diversity below 4 °C leading to higher SST values during glacial maxima⁶³.

A paleoceanographic map of corrected *N. pachyderma* Mg/Ca records covering the LGM and Late Holocene demonstrates that the underestimated cooling (e.g., ΔSST) is most pronounced across mid-latitudes where our estimates suggest cooling between 6–9 ± 1.0 °C instead of milder 4-6° C cooling (Fig. 3a). Our re-assessment suggests that polar amplification or high-latitude LGM cooling was of a greater magnitude and characterized by a more extensive spatial expression than previously estimated, which raises questions about our understanding of heat transport under glacial boundary conditions. Furthermore, absolute Holocene and LGM SSTs values based on corrected *N. pachyderma* Mg/Ca records (Fig. 3b, c) are much closer to model-based reconstructions of Holocene and LGM SSTs than previous estimates based on this species (Supplementary Fig. S5). Considering these examples, we believe that a systematic application of the calibration approach, as presented in our study, provides, an opportunity to directly estimate past temperatures in Polar Oceans, where low diversity of planktonic foraminifera prevented the use of transfer functions and where dinocysts are rare.

Furthermore, these corrected datasets demonstrate that paired Mg/Ca · δ¹⁸O_c records, including reanalysis of published data, can now answer the call by the modelling community to reduce uncertainties linked to high latitude SST reconstructions and offer a more robust means of data-model comparison. In particular, Mg/Ca-based SST reconstructions in *N. pachyderma* can now evaluate the ability of climate models to simulate polar amplified warming and cooling, which is especially important as climate model simulations targeting warmer than present climates have historically not captured the full extent of polar amplified warming¹²,⁶⁴. Further, more accurate high-latitude paleoclimate data will inform the latest generation of climate models that are using paleoclimate-calibrated Community Earth System Models (CESM2)¹²,⁶⁵ to simulate multiple paleoclimate states—warmer and colder than the present. Finally, the corrected polar paleoclimate datasets will be a critical tool to refine estimates on the sensitivity of ice sheets to melting (e.g., sea level rise during warm climates⁶⁶), and to quantify how AMOC strength affects ocean heat transport on orbital to millennial timescales through glacial-interglacial cycles⁶⁷,⁶⁸.

# Methods

## Sample selection & preparation

*N. pachyderma* specimens were picked from multinet (HydroBios, type Midi) plankton tows collected during four separate oceanographic cruises, MSM09 (2008), MSM44 (2015), and MSM66 (2017) to the Labrador Sea/Baffin Bay and CE20009 (2020) to the Nordic Seas (see Supplementary Tables 1–3). Foraminifera were either immediately picked and placed into slides on board and subsequently frozen or directly frozen at –80 °C to prevent degradation of the cytoplasm before picking. This also allowed for clear separation of specimens that were alive when collected (cytoplasm bearing) from specimens that were dead (empty tests) for subsequent geochemical analysis. Specimen were divided into four size fractions where possible: 100–150 µm, 150–200 µm, 200–250 µm and 250–300 µm. Weights per individual were determined using a Sartorius XM1000P microbalance (nominal resolution of 1 µg) at the University of Bremen and a Sartorius Semi-microbalance SECURA (nominal resolution of 2 µg) at the University of Galway before δ¹⁸O_c and Mg/Ca analysis. Given the low average weights of tow specimens (e.g., 2.12 µg per individual), at least 20 specimens

per size fraction were required for δ¹⁸O_c and, on average, ~320 individuals for Trace Element (TE) analysis. To determine individual test weights, the total number of individuals measured was divided by the final weight per sample. Core tops analysed in this study were washed and dried to separate fine from coarse fractions. Between 250 and 350 specimens were picked within the size fraction 200–250 µm. While we selected specimens from this narrow size fraction, we did not perform a systematic analysis of the degree of encrustation of *N. pachyderma* in the core top samples.

## Hydrographic data & processing

Water column carbonate chemistry for each tow station was derived from hydrographic surveys carried out as close to the date and station as possible. For example, for stations from MSM 09, the CLIVAR and Davies Strait Time Series were collected sometimes within days (at most 1 month) and within 1 degree latitude and longitude of tow samples. When more than one station was available within a 1-degree radius, average values were calculated. Supplementary Table 1 provides a detailed list of datasets and stations used for carbonate chemistry, including coordinates and dates. For some stations, carbonate chemistry data from the same year and area was not available; in these cases, only profiles collected during the same season were considered and evaluated against available CTD data. For CE20009, paired hydrographic data and water samples, including ALK and DIC, were collected. ALK and DIC were analysed simultaneously from the same sample at the University of Galway using a VINDTA-3C⁶⁹ following the methods outlined in Dickson et al.⁷⁰. The VINDTA-3C was calibrated by duplicate measurement of Certified Reference Materials (A. Dickson, Scripps Institute of Oceanography, USA) with every batch of samples. The remaining carbonate system parameters were calculated using the Seacarb package in R.

For the core top datasets, calcification depth and hydrographic parameters were estimated by comparing the measured foraminiferal δ¹⁸O_c values to the expected equilibrium values of calcite (δ¹⁸O_c) precipitated from seawater. For the Kozdon dataset, we note that specimens for Mg/Ca and δ¹⁸O_c analysis were picked from the same sample but are not paired i.e., they constitute two separate sets of specimens. To determine calcification depth, we selected temperature, salinity, and carbonate chemistry using CARINA⁵⁸, while paired hydrographic profiles were used for CE20009 core tops. Since the dominant flux of *N. pachyderma* occurs in the summer, we only selected profiles collected between June and October from stations as close to the core sites as possible for the Kozdon dataset. Specific stations, locations and dates of profiles are detailed in the Supplementary Table 8. To obtain δ¹⁸O_sw depth profiles, we used modelled water-column data from LeGrande and Schmidt⁷¹. We retrieved data from both the closest station to each core site and within a 0.5° radius to provide locally averaged δ¹⁸O_sw values. Despite these constraints, it is possible that the hydrographic profiles used here do not represent the calcification environment of specimens recovered from core tops. For example, unrealistically deep calcification habitat estimates below 300 m may result from modern hydrographic profiles that are “too warm” for samples that predate modern global warming. For the CE20009 core tops, paired water column δ¹⁸O_sw values are available. We then calculated site-specific δ¹⁸O_sw - salinity linear relationships for each core. Using the δ¹⁸O - paleotemperature equation developed by O'Neil et al.⁷², we determined expected δ¹⁸O_c depth profiles (0–400 m) based on seasonal temperature, salinity, and δ¹⁸O_sw data. Calculations were completed by subtracting 0.27‰ to compare measured δ¹⁸O values of CO₂ produced by the reaction of calcite with H₃PO₄ and CO₂ equilibrated with water⁷³. The depth at which the expected equilibrium δ¹⁸O_c values matched the measured foraminiferal δ¹⁸O_c values (determined by linear interpolation) is the calcification depth, and the calcification temperature, salinity, and carbonate chemistry were determined based on respective values recorded at this calcification depth.
<figure><img src="image_6.png" /></figure>

**Fig. 3 | Data model inter-comparison.** Panel (a) shows the difference in sea surface temperature cooling (ΔSST) between the Holocene and the Last Glacial Maximum calculated using the isotope-enabled Community Earth System Model (iCESM) as shown in Tierney et al. 28. In the background together with ΔSST calculated using corrected Mg/Ca values measured on *N. pachyderma* (symbols). For a list of all datasets used, please see Supplementary Table 7. Panels (b) and (c) show Holocene and LGM SSTs, respectively, with the modelled (iCESM) SSTs in the background and

corrected Mg/Ca-based SSTs shown as symbols. Maps were produced with Ocean Data View⁹³. Please note that modelled annual SSTs are likely too cold for a direct comparison with *N. pachyderma*-based SSTs, given their seasonal reproductive cycle. Maps displaying uncorrected Mg/Ca values against modelled paleo data are shown in Supplementary Fig. S5. Source data for all panels are provided as a Source Data file.
We did not apply a systematic correction for a "vital effect" on δ¹⁸O_c values because the application of such a vital effect when estimating calcification depth for core top *N. pachyderma* is complicated by: (1) the δ¹⁸O_c offset from equilibrium measured in *N. pachyderma* collected from plankton tows is highly variable within a region, across regions, seasons, and on interannual timescales²⁶,⁷⁴⁻⁷⁶; (2) the addition of crust, typical for specimens recovered from core tops, increases the δ¹⁸O_c signal because crust is isotopically heavier than ontogenetic calcite⁵⁷,⁷⁷. Variations in the degree of crusting can, therefore, either "lower", completely "mask", or shift the vital effect towards positive values depending on the degree of crusting⁵⁷. As a result, there remains considerable uncertainty in the exact value to use when correcting for the competing signals of vital effects and crust in *N. pachyderma*, with studies reporting both calcification in equilibrium⁷⁷⁻⁸⁰ and out of equilibrium with seawater⁵⁵,⁸¹⁻⁸³.

Since it is not possible to isolate the vital effect from crust in core top specimens the application of a fixed offset leads to assumptions of calcification habitats that are not grounded in ecological observations. Specifically, it has long been assumed that the isotopically heavier crust calcite results from systematic ontogenic vertical migration to deeper and colder waters associated with reproduction³¹. However, from comprehensive, cross-regional plankton tow studies, we know that heavily crusted specimens can be found at any depth, including the top 0–50 m³² suggesting that there is no systematic ontogenic vertical migration associated with reproduction or “crusting” for *N. pachyderma*³²,³³. Furthermore, culture-grown crust calcite in *N. incompta* has a lower Mg/Ca ratio than ontogenetic calcite formed at the same temperature, suggesting that a lower temperature is not responsible for the low Mg/Ca ratio of neogloboquadrinid crusts²⁵. Finally, cross-sections of *N. pachyderma* tests from North Atlantic sediment traps also show that tests growing in a well-mixed water column with constant temperatures exhibit lower Mg/Ca ratios in their crust³⁴, suggesting that alongside temperature, other factors control the difference between crust and ontogenetic Mg/Ca.

By not applying a correction for vital effects or the contribution of crust calcite we avoid biasing calcification depth estimates towards deeper depth. This approach is also consistent with sediment trap studies from both the North and South Atlantic, where δ¹⁸O_c values measured on crusted *N. pachyderma* are in equilibrium with measured water column δ¹⁸O_c⁷⁷,⁷⁸. Except for two stations that yield anomalously deep estimates, all calcification depths fall between 7 and 67 ± 22 m 2σ water depth (Supplementary Table 8). Compared to the existing literature, these are at the shallow end of previously reported depth estimates⁸² but are in line with a seasonal calcification depth of 50 m⁸⁰ and observations that place the main seasonal bloom of *N. pachyderma* at the base of the chlorophyll max following the breakup of sea-ice⁸⁴.

For South Atlantic sediment traps included here, we note that crusted specimens were specifically excluded from the analysis to avoid the contribution of secondary or gametogenic crust calcite²⁴. We, therefore, applied the vital effect as derived in Hendry et al²⁴, when estimating [CO₂⁻] values and lowered Mg/Ca ratios by 15% to account for the difference in Mg between ontogenetic and crust calcite³⁴,⁴⁸. No vital effect or correction on Mg/Ca values was applied to the core top dataset from the South Atlantic⁵⁹ consistent with the processing of North Atlantic core top datasets analysed here.

## Analytical methods

To prepare foraminifera from plankton tows from MSM cruises for TE analysis, tests were gently crushed between two glass plates and transferred into acid-leached vials. The crushed foraminiferal tests were cleaned using a cleaning protocol focussing on the removal of organic matter (i.e., the oxidative step was repeated twice) followed by a final dissolution in dilute HNO₃. We did not pre-rinse samples with Milli-Q or include a reductive cleaning step since clay or metal oxide

removal is not required for tow samples. Samples from MSM cruises were analysed using a Nu AttoM high-resolution double-focusing inductively coupled plasma mass spectrometer at the Alfred Wegner Institute in Bremerhaven, Germany. Five replications were carried out for each sample, where the average repeatability for Mg/Ca among the five replicates was 0.95 ± 0.47% (relative standard deviation). Long-term reproducibility of Mg/Ca analyses is monitored by regular analysis of the carbonate standard JCp-1, yielding 3.94 ± 0.15 (2σ) mmol.mol⁻¹ over two years. Samples from cruise CE20009 were processed for elemental analysis at the University of Southampton. They were cleaned using established methods following Barker et al.⁸⁵, clay was removed via multiple Milli-Q rinses (for core top samples only) and the tests were oxidatively cleaned. Elemental analysis was performed on a Thermo Element Inductively Couple Plasma Mass Spectrometer (ICPMS) following Henehan et al.⁸⁶. Stable isotope analysis on MSM samples was performed on a Finnigan MAT 252 mass spectrometer equipped with an automated carbonate preparation device at MARUM, University of Bremen. Isotopic results were calibrated relative to the Vienna Pee Dee Belemnite (VPDB) using the NBS19 standard. The standard deviation of the laboratory standard was better than 0.07‰ for δ¹⁸O_c over the measuring period. Samples from CE20009 were analysed at FARLAB, using a Thermo Scientific MAT 253 IRMS coupled to a Kiel IV automated carbonate preparation device. Results are reported on the VPDB scale calibrated with NBS19, together with NBS18, and in-house synthetic carbonate standards CM12 and Riedal used in daily operations at FARLAB, University of Bergen. Long-term external precision (1 σ SD) is better than 0.08‰ for δ¹⁸O_c based on replicate analysis of the in-house standard CM12.

## PaleDATA correction

For all previously published datasets, δ¹⁸O_c data were corrected for glacial-interglacial sea-level changes using sea-level reconstructions by Spratt and Lisiecki⁸⁷. Here we use the commonly cited value of 0.11‰/ 10 m SL change for the correction of δ¹⁸O_c⁸⁸⁻⁹¹. We also applied a +10–15% offset to Mg/Ca[meas] values that underwent a reductive cleaning step (e.g., ref. 29,30,) or repeated acid leaches⁹² since these additional steps have been shown to lower Mg/Ca values in planktonic foraminifera⁸⁵. We then derived paleo [CO₂⁻] values from δ¹⁸O_c and corrected Mg/Ca[meas] using the correction scheme presented here. For the Late Holocene – LGM compilation, we included as many available records from literature as possible. Supplementary Table 7 details core sites, time intervals and references used for each estimate. For the calculation of mean glacial and interglacial values covering the last five glacial-interglacial cycles from IODP site U1305, please see Supplementary Table 6.

## Calibration uncertainty calculations

For the *N. pachyderma* data, first the relationship between δ¹⁸O and carbonate ion concentration is assessed (Fig. 1d). We performed a Monte Carlo simulation (n = 1000) on the slope of this relationship using the analytical uncertainty on the oxygen isotope measurements (typically 0.2 at 2σ) and carbonate ion concentration from hydrographic data (GLODAP) with a conservative prescribed uncertainty of ± 20 µmol/kg (2σ). This gives a range of potential slopes which can be used to determine a central estimate and uncertainties defined by the distribution of potential slopes (inside 95% of the potential slopes). The central coefficients are used to calculate the CO₂⁻:Mg/Ca[norm] relationship alongside the error on this slope to determine the uncertainties on the x-axis (Fig. 1b). These uncertainties are then propagated into the CO₂⁻:Mg/Ca[norm] relationship, where CO₂⁻ as given by δ¹⁸O is plotted against normalised Mg/Ca (Mg/Ca[meas]/Mg/Ca[pred]) to obtain a Mg/Ca[corr] for non-thermal effects alongside its uncertainties. Based on the previous relationships (Fig. 1a, b), Mg/Ca[corr] has an uncertainty of ± 0.2 (2σ), which we then calibrate against the in situ temperatures obtained from CTDs and/or the CARINA database.
These errors are related to the depth estimation of the foraminifera from the isotope data but typically range from 0 - 3 °C (minimum 0 °C, maximum 10 °C, average ± 2.8 °C with a standard deviation of 0.9 °C). Note, for the data from Meland et al.²⁰, we are unable to calculate uncertainty, so we use the average of the uncertainties from Kozdon et al.²⁷, which is 1.6 °C at 2σ). For the *N. incompta* data, all the steps to calculate the uncertainties are identical to the *N. pachyderma* dataset except, as there are no δ¹⁸O data, we assigned a temperature uncertainty of ± 0.8 °C (2σ)³⁹. The Mg/Ca$_{[corr]}$ uncertainty is then assumed to be the same as for *N. pachyderma*. Both species are combined to make the final calibration but using one or the other only gives consistent relationships within error, strengthening the argument against a species-specific vital effect.

The calibration Monte Carlo simulation finally is conducted on 1000 iterations of a slope using the propagated Mg/Ca[corr] and the individual temperature uncertainties (as above). We take the median slope of this curve to define the Mg/Ca[corr] to temperature relationship, and the 68 and 95 % confidence intervals (1 and 2σ respectively) to describe its uncertainties. These uncertainties represent the uncertainty in the defined relationship, assuming it to be consistent and are associated with a temperature reconstruction uncertainty of ±1.2 °C at temperatures of around 0 °C and ±0.9 °C at temperatures of around 15 °C. These confidence intervals give a conservative uncertainty on the relationship between Mg/Ca[corr] and temperature and are used for any relative change calculations presented in this study. We also plot a prediction uncertainty, which can be considered the uncertainty from a single future or additional data point but is overly conservative for the whole population of data.

**Data availability**

All new data and recalculated datasets shown here are available in the Supplementary Information/Source Data file. Source data are provided in this paper.

**Code availability**

The fully commented code for the error propagation calculation and input files have been deposited in the Zenodo database under DOI: 10.5281/zenodo.13833954.

**References**

1. Holland, M. M. & Bitz, C. M. Polar amplification of climate change in coupled models. *Clim. Dyn.* **21**, 221–232 (2003).

2. Manabe, S. & Stouffer, R. J. Sensitivity of a global climate model to an increase of CO2 concentration in the atmosphere. *J. Geophys. Res. Oceans* **85**, 5529–5554 (1980).

3. Årthun, M. et al. Skillful prediction of northern climate provided by the ocean. *Nat. Commun.* **8**, 1–11 (2017).

4. Årthun, M., Eldevik, T. & Smedsrud, L. H. The role of Atlantic heat transport in future Arctic winter sea ice loss. *J. Clim.* **32**, 3327–3341 (2019).

5. Parmentier, F.-J. W. et al. The impact of lower sea-ice extent on Arctic greenhouse-gas exchange. *Nat. Clim. Change* **3**, 195–202 (2013).

6. Jansen, E. et al. Past perspectives on the present era of abrupt Arctic climate change. *Nat. Clim. Change* **10**, 714–721 (2020).

7. Vavrus, S. J. The influence of Arctic amplification on mid-latitude weather and climate. *Curr. Clim. Change Rep.* **4**, 238–249 (2018).

8. Boetius, A. et al. Microbial ecology of the cryosphere: sea ice and glacial habitats. *Nat. Rev. Microbiol.* **13**, 677–690 (2015).

9. Huang, J. et al. Recently amplified arctic warming has contributed to a continual global warming trend. *Nat. Clim. Change* **7**, 875–879 (2017).

10. Meredith, M. et al. Polar regions. chapter 3, ipcc special report on the ocean and cryosphere in a changing climate. (2019).

11. Cowtan, K. & Way, R. G. Coverage bias in the HadCRUT4 temperature series and its impact on recent temperature trends. *Q. J. R. Meteorol. Soc.* **140**, 1935–1944 (2014).

12. Burls, N. & Sagoo, N. Increasingly sophisticated climate models need the out-of-sample tests paleoclimates provide. *J. Adv. Modelling Earth Syst.* **14**, e2022MS003389 (2022).

13. Fedorov, A. V. et al. Tightly linked zonal and meridional sea surface temperature gradients over the past five million years. *Nat. Geosci.* **8**, 975–980 (2015).

14. Polyakov, I. V. et al. Greater role for Atlantic inflows on sea-ice loss in the Eurasian Basin of the Arctic Ocean. *Science* **356**, 285–291 (2017).

15. Kucera, M. et al. Reconstruction of sea-surface temperatures from assemblages of planktonic foraminifera: multi-technique approach based on geographically constrained calibration data sets and its application to glacial Atlantic and Pacific Oceans. *Quat. Sci. Rev.* **24**, 951–998 (2005).

16. Weinelt, M. et al. Paleogeographic proxies in the northern North Atlantic. *Northern North Atlantic Changing Env.* 319–352 (2001).

17. Telford, R. J. Limitations of dinoflagellate cyst transfer functions. *Quat. Sci. Rev.* **25**, 1375–1382 (2006).

18. Bendle, J., Rosell-Melé, A. & Ziveri, P. Variability of unusual distributions of alkenones in the surface waters of the Nordic seas. *Paleoceanography* **20**, (2005).

19. 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).

20. Meland, M. Y. et al. Mg/Ca ratios in the planktonic foraminifer Neogloboquadrina pachyderma (sinistral) in the northern North Atlantic/Nordic Seas. *Geochem. Geophys. Geosyst.* **7**, https://doi.org/10.1029/2005GC001078 (2006).

21. Katz, M. E. et al. Traditional and emerging geochemical proxies in foraminifera. *J. Foraminifer. Res.* **40**, 165–192 (2010).

22. Katz, A. The interaction of magnesium with calcite during crystal growth at 25–90 C and one atmosphere. *Geochim. Cosmochim. Acta* **37**, 1563IN31579-15781586 (1973).

23. Nuernberg, D. Magnesium in tests of Neogloboquadrina pachyderma sinistral from high northern and southern latitudes. *J. Foraminifer. Res.* **25**, 350–368 (1995).

24. Hendry, K. R. et al. Controls on stable isotope and trace metal uptake in Neogloboquadrina pachyderma (sinistral) from an Antarctic sea-ice environment. *Earth Planet. Sci. Lett.* **278**, 67–77 (2009).

25. Davis, C. V. et al. Relationships between temperature, pH, and crusting on Mg/Ca ratios in laboratory-grown neogloboquadrina foraminifera. *Paleoceanography* **32**, 1137–1152 (2017).

26. Livsey, C. M. et al. High-resolution Mg/Ca and δ18O patterns in modern Neogloboquadrina pachyderma from the Fram Strait and Irminger Sea. *Paleoceanogr. Paleoclimatol.* **35**, e2020PA003969 (2020).

27. Kozdon, R. et al. Reassessing Mg/Ca temperature calibrations of Neogloboquadrina pachyderma (sinistral) using paired δ44/40Ca and Mg/Ca measurements. *Geochem. Geophys. Geosyst.* **10**, (2009).

28. Tierney, J. E. et al. Glacial cooling and climate sensitivity revisited. *Nature* **584**, 569–573 (2020).

29. Irvali, N. et al. Evidence for regional cooling, frontal advances, and East Greenland Ice Sheet changes during the demise of the last interglacial. *Quat. Sci. Rev.* **150**, 184–199 (2016).

30. Irvali, N. et al. A low climate threshold for south Greenland Ice Sheet demise during the Late Pleistocene. *Proc. Natl. Acad. Sci. USA* **117**, 190–195 (2020).
31. Emiliani, C. Depth habitats of some species of pelagic foraminifera as indicated by oxygen isotope ratios. *Am. J. Sci.* **252**, 149-158 (1954).

32. Tell, F. et al. Upper-ocean flux of biogenic calcite produced by the Arctic planktonic foraminifera Neogloboquadrina pachyderma. *Biogeosciences* **19**, 4903-4927 (2022).

33. Manno, C. & Pavlov, A. Living planktonic foraminifera in the Fram Strait (Arctic): absence of diel vertical migration during the mid- night sun. *Hydrobiologia* **721**, 285-295 (2014).

34. Jonkers, L. et al. Chamber formation leads to Mg/Ca banding in the planktonic foraminifer Neogloboquadrina pachyderma. *Earth Planet. Sci. Lett.* **451**, 177-184 (2016).

35. Evans, D. et al. Revisiting carbonate chemistry controls on planktic foraminifera Mg/Ca: implications for sea surface temperature and hydrology shifts over the Paleocene-Eocene Thermal Maximum and Eocene-Oligocene transition. *Climate* **12**, 819-835 (2016).

36. Kisakürek, B. et al. Controls on shell Mg/Ca and Sr/Ca in cultured planktonic foraminiferan, Globigerinoides ruber (white). *Earth Planet. Sci. Lett.* **273**, 260-269 (2008).

37. Russell, A. D. et al. Effects of seawater carbonate ion concentration and temperature on shell U, Mg, and Sr in cultured planktonic foraminifera. *Geochim. Cosmochim. Acta* **68**, 4347-4361 (2004).

38. Gray, W. R. et al. The effects of temperature, salinity, and the carbonate system on Mg/Ca in Globigerinoides ruber (white): A global sediment trap calibration. *Earth Planet. Sci. Lett.* **482**, 607-620 (2018).

39. Morley, A. et al. Environmental Controls on Mg/Ca in Neogloboquadrina incompta: A Core-Top Study From the Subpolar North Atlantic. *Geochem. Geophys. Geosyst.* **18**, 4276-4298 (2017).

40. Gray, W. R. & Evans, D. Nonthermal influences on Mg/Ca in planktonic foraminifera: a review of culture studies and application to the Last Glacial Maximum. *Paleoceanogr. Paleoclimatol.* **34**, 306-315 (2019).

41. Tierney, J. E. et al. Bayesian calibration of the Mg/Ca paleothermometer in planktic foraminifera. *Paleoceanogr. Paleoclimatol.* **34**, 2005-2030 (2019).

42. Foster, G. L. & Rae, J. W. Reconstructing ocean pH with boron isotopes in foraminifera. *Annu. Rev. Earth Planet. Sci.* **44**, 207-237 (2016).

43. Darling, K. F. et al. A resolution for the coiling direction paradox in Neogloboquadrina pachyderma. *Paleoceanography* **21**, https://doi.org/10.1029/2005PA001189 (2006).

44. Schiebel, R. & Hemleben, C. *Planktic Foraminifers in the Modern Ocean*. (2017).

45. Evans, D., Müller, W. & Erez, J. Assessing foraminifera biomimeralisation models through trace element data of cultures under variable seawater chemistry. *Geochim. Cosmochim. Acta* **236**, 198-217 (2018).

46. Zeebe, R. E. & Sanyal, A. Comparison of two potential strategies of planktonic foraminifera for house building: Mg2+ or H+ removal? *Geochim. Cosmochim. Acta* **66**, 1159-1169 (2002).

47. Jonkers, L. et al. Encrustation and trace element composition of Neogloboquadrina dutertrei assessed from single chamber analyses-implications for paleotemperature estimates. *Bio-geosciences* **9**, 4851-4860 (2012).

48. Hupp, B. N. & Fehrenbacher, J. S. Geochemical differences between alive, uncrusted and dead, crusted shells of Neogloboquadrina pachyderma: Implications for paleoreconstruction. *Paleoceanogr. Paleoclimatol.* **38**, e2023PA004638 (2023).

49. Ravelo, A. C. & Hillaire-Marcel, C. Chapter eighteen the use of oxygen and carbon isotopes of foraminifera in paleoceanography. *Dev. Mar. Geol.* **1**, 735-764 (2007).

50. Emiliani, C. Pleistocene temperatures. (1955).

51. Epstein, S. et al. Carbonate-water isotopic temperature scale. *Geol. Soc. Am. Bull.* **62**, 417-426 (1951).

52. Spero, H. J. et al. Effect of seawater carbonate concentration on foraminiferal carbon and oxygen isotopes. *Nature* **390**, 497-500 (1997).

53. Bijma, J., Spero, H. & Lea, D. in *Use of Proxies in Paleooceanography*. 489-512 (Springer,1999).

54. Ziveri, P. et al. A universal carbonate ion effect on stable oxygen isotope ratios in unicellular planktonic calcifying organisms. *Bio-geosciences* **9**, 1025-1032 (2012).

55. Bauch, D., Carstens, J. & Wefer, G. Oxygen isotope composition of living Neogloboquadrina pachyderma (sin.) in the Arctic Ocean. *Earth Planet. Sci. Lett.* **146**, 47-58 (1997).

56. Lombard, F. et al. Modelling the temperature dependent growth rates of planktic foraminifera. *Mar. Micropaleontol.* **70**, 1-7 (2009).

57. Kozdon, R. et al. Intratest oxygen isotope variability in the planktonic foraminifer N. pachyderma: Real vs. apparent vital effects by ion microprobe. *Chem. Geol.* **258**, 327-337 (2009).

58. Key, R. et al. The CARINA data synthesis project: introduction and overview. *Earth Syst. Sci. Data* **2**, 105-121 (2010).

59. Vázquez Riveiros, N. et al. Mg/Ca thermometry in planktic foraminifera: Improving paleotemperature estimations for G. bulloides and N. pachyderma left. *Geochem. Geophys. Geo-syst.* **17**, 1249-1264 (2016).

60. Darling, K. F. et al. Molecular evidence links cryptic diversification in polar planktonic protists to Quaternary climate dynamics. *Proc. Natl. Acad. Sci. USA* **101**, 7657-7662 (2004).

61. Thornalley, D. J., I. McCave, N. & Elderfield, H. Freshwater input and abrupt deglacial climate change in the North Atlantic. *Paleoceanography* **25**, https://doi.org/10.1029/2009PA001772 (2010).

62. Pflaumann, U. et al. Glacial North Atlantic: Sea-surface conditions reconstructed by GLAMAP 2000. *Paleoceanography* **18**, 1065 (2003).

63. De Vernal, A. et al. Comparing proxies for the reconstruction of LGM sea-surface conditions in the northern North Atlantic. *Quat. Sci. Rev.* **25**, 2820-2834 (2006).

64. Masson-Delmotte, V. et al. Past and future polar amplification of climate change: climate model intercomparisons and ice-core constraints. *Clim. Dyn.* **26**, 513-529 (2006).

65. Zhu, J. et al. LGM paleoclimate constraints inform cloud parameterizations and equilibrium climate sensitivity in CESM2. *J. Adv. Modeling Earth Syst.* **14**, e2021MS002776 (2022).

66. Dutton, A. et al. Sea-level rise due to polar ice-sheet mass loss during past warm periods. *Science* **349**, aaa4019 (2015).

67. Menviel, L. C. et al. An ice-climate oscillatory framework for Dansgaard-Oeschger cycles. *Nat. Rev. Earth Environ.* **1**, 677-693 (2020).

68. Kageyama, M. et al. Modeling the climate of the Last Glacial Maximum from PMIP1 to PMIP4. *Glob. Chang. Mag.* **29**, 80-81 (2021).

69. Mintrop, L. et al. Alkalinity determination by potentiometry: Inter-calibration using three different methods. **26**, https://doi.org/10.7773/cm.v26i1.573 (2000).

70. Dickson, A., Sabine, C. & Christian, J. Guide to best practices for ocean CO2 measurements. *PICES Spec. Publ.* **3**, IOCCP report (2007).

71. LeGrande, A. N. & Schmidt, G. A. Global gridded data set of the oxygen isotopic composition in seawater. *Geophys. Res. Lett.* **33**, L12604 (2006).

72. O'Neil, J. R., Clayton, R. N. & Mayeda, T. K. Oxygen isotope fractionation in divalent metal carbonates. *J. Chem. Phys.* **51**, 5547-5558 (1969).

73. Friedman, I. & O'Neil, J. R. *Data of geochemistry: Compilation of stable isotope fractionation factors of geochemical interest.* 440: US Government Printing Office. (1977).
74. Stangeew, E., Distribution and Isotopic Composition of Living Planktonic Foraminifera N. pachyderma (sinistral) and T. quinqueloba in the High Latitude North Atlantic. (2001).

75. Pados, T. et al. Oxygen and carbon isotope composition of modern planktic foraminifera and near-surface waters in the Fram Strait (Arctic Ocean)-a case study. *Biogeosciences* **12**, 1733-1752 (2015).

76. Volkmann, R. & Mensch, M. Stable isotope composition (δ18O, δ13C) of living planktic foraminifers in the outer Laptev Sea and the Fram Strait. *Mar. Micropaleontol.* **42**, 163-188 (2001).

77. Mikis, A. et al. Temporal variability in foraminiferal morphology and geochemistry at the West Antarctic Peninsula: a sediment trap study. *Biogeosciences* **16**, 3267-3282 (2019).

78. Jonkers, L. et al. Seasonal patterns of shell flux, δ18O and δ13C of small and large N. pachyderma (s) and G. bulloides in the subpolar North Atlantic. *Paleoceanography* **28**, 164-174 (2013).

79. Jonkers, L. et al. Seasonal stratification, shell flux, and oxygen isotope dynamics of left-coiling N. pachyderma and T. quinqueloba in the western subpolar North Atlantic. *Paleoceanography* **25**, https://doi.org/10.1029/2009PA001849 (2010).

80. Jonkers, L. et al. Variability in Neogloboquadrina pachyderma stable isotope ratios from isothermal conditions: implications for individual foraminifera analysis. *Climate* **18**, 89-101 (2022).

81. Kohfeld, K. E. et al. Neogloboquadrina pachyderma (sinistral coiling) as paleoceanographic tracers in polar oceans: Evidence from Northeast Water Polynya plankton tows, sediment traps, and surface sediments. *Paleoceanography* **11**, 679-699 (1996).

82. Simstich, J., Sarnthein, M. & Erlenkeuser, H. Paired δ18O signals of Neogloboquadrina pachyderma (s) and Turborotalita quinqueloba show thermal stratification structure in Nordic Seas. *Mar. Micropaleontol.* **48**, 107-125 (2003).

83. King, A. L. & Howard, W. R. δ18O seasonality of planktonic foraminifera from Southern Ocean sediment traps: Latitudinal gradients and implications for paleoclimate reconstructions. *Mar. Micropaleontol.* **56**, 1-24 (2005).

84. Greco, M. et al. Depth habitat of the planktonic foraminifera Neogloboquadrina pachyderma in the northern high latitudes explained by sea-ice and chlorophyll concentrations. *Biogeosciences* **16**, 3425-3437 (2019).

85. Barker, S., Greaves, M. & Elderfield, H. A study of cleaning procedures used for foraminiferal Mg/Ca paleothermometry. *Geochem. Geophys. Geosyst.* **4**, 8407 (2003).

86. Henehan, M. J. et al. Evaluating the utility of B/C ratios in planktic foraminifera as a proxy for the carbonate system: A case study of G lobigerinoides ruber. *Geochem., Geophys. Geosyst.* **16**, 1052-1069 (2015).

87. Spratt, R. M. & Lisiecki, L. E. A Late Pleistocene sea level stack. *Climate* **12**, 1079-1092 (2016).

88. Fairbanks, R. G. A 17,000 year glacio-eustatic sea level record: influence of glacial melting rate on the Younger Dryas event and deep-ocean circulation. *Nature* **342**, 637-642 (1989).

89. Schrag, D. P., Hampt, G. & Murray, D. W. Pore fluid constraints on the temperature and oxygen isotopic composition of the glacial ocean. *Science* **272**, 1930-1932 (1996).

90. Adkins, J. F., McIntyre, K. & Schrag, D. P. The salinity, temperature, and δ18O of the glacial deep ocean. *Science* **298**, 1769-1773 (2002).

91. Rohling, E. et al. Sea-level and deep-sea-temperature variability over the past 5.3 million years. *Nature* **508**, 477-482 (2014).

92. Ezat, M. M., Rasmussen, T. L. & Groeneveld, J. Reconstruction of hydrographic changes in the southern Norwegian Sea during the past 135 kyr and the impact of different foraminiferal Mg/Ca cleaning protocols. *Geochem., Geophys. Geosyst.* **17**, 3420-3436 (2016).

93. Schlitzer, R., Ocean data view. 2022.

## Acknowledgements

This research was funded by MSCA-IF Project ARCTICO [838529], the Marine Institute of Ireland Research Programme 2014-2020 (PDOC/19/05/02), and Science Foundation Ireland Frontiers for the Future Project 21/FFP-P/10261, awarded to A.M. In addition, A.M. acknowledges Grant in Aid funding from the Marine Institute for research expedition CE20009 on the RV Celtic Explorer with special thanks to the crew sailing under Master Anthony Hobin and Research Assistant Aedin McAlear for the analysis of Ocean ALK and DIC samples from the CE20009 survey. M.R. acknowledges DFG (German Research Foundation) for funding through research grant no. RA 2068/4-1. R.C. acknowledges funding for Project CE2COAST, funded by ANR (FR), BELSPO (BE), FCT (PT), IZM (LV), MI (IE), MIUR (IT), Rannis (IS) and RCN (NO) through the 2019 "Joint Transnational Call on Next Generation Climate Science in Europe for Oceans" initiated by JPI Climate and JPI Oceans.

## Author contributions

A.M. designed the project, carried out palaeoceanographic sample preparation of plankton tows from the Labrador Sea, compiled and analysed all datasets, wrote the manuscript, and coordinated the input of all co-authors. E.D.L.V. carried out palaeoceanographic sample preparation and geochemical analysis of plankton tows and core tops from CE20009 and contributed to the writing of the final manuscript. M.R. performed trace element and boron isotope analysis on *N. pachyderma* plankton tow samples from the Labrador Sea. J.B. supervised A.M. at the Alfred Wegner Institute and oversaw trace element and boron isotope analysis. U.N. oversaw all stable isotope analyses of water and foraminifera samples from CE20009 and contributed to the final version of this manuscript. G.F. oversaw trace element analysis on *N. pachyderma* plankton tow samples from CE20009 and contributed to the final version of this manuscript. T.C. performed uncertainty calculations for the Mg/Ca correction and contributed to the final version of this manuscript. J.M. collected all plankton tows during CE20009 and contributed to the final version of this manuscript. R.C. oversaw Alkalinity and DIC analysis on seawater samples collected during CE20009 and contributed to the final version of this manuscript. J.B. performed Alkalinity and DIC analysis on seawater samples collected during CE20009 and contributed to the final version of this manuscript. M.K. advised A.M. during project design, supervised her during her MSCA-IF at MARUM, University of Bremen, and contributed to the final version of this manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-024-53424-w.

**Correspondence** and requests for materials should be addressed to A. Morley.

**Peer review information** *Nature Communications* thanks Reinhard Kozdon and the other anonymous reviewers for their contribution to the peer review of this work. A peer review file is available.

**Reprints and permissions 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.
**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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024