Last Interglacial synthesis of high-latitude temperature: temperature anomalies and associated errors for 4 time slices

The Last Interglacial (LIG, 129-116 thousand of years BP, ka) represents a test bed for climate model feedbacks in warmer-than-present high latitude regions. However, mainly because aligning different palaeoclimatic archives and from different parts of the world is not trivial, a spatio-temporal picture of LIG temperature changes is difficult to obtain. Here, we have selected 47 polar ice core and sub-polar marine sediment records and developed a strategy to align them onto the recent AICC2012 ice core chronology. We provide the first compilation of high-latitude temperature changes across the LIG associated with a coherent temporal framework built between ice core and marine sediment records. Our new data synthesis highlights non-synchronous maximum temperature changes between the two hemispheres with the Southern Ocean and Antarctica records showing an early warming compared to North Atlantic records. We also observe warmer than present-day conditions that occur for a longer time period in southern high latitudes than in northern high latitudes. Finally, the amplitude of temperature changes at high northern latitudes is larger compared to high southern latitude temperature changes recorded at the onset and the demise of the LIG. We have also compiled four data-based time slices with temperature anomalies (compared to present-day conditions) at 115 ka, 120 ka, 125 ka and 130 ka and quantitatively estimated temperature uncertainties that include relative dating errors. This provides an improved benchmark for performing more robust model-data comparison. The surface temperature simulated by two General Circulation Models (CCSM3 and HadCM3) for 130 ka and 125 ka is compared to the corresponding time slice data synthesis. This comparison shows that the models predict warmer than present conditions earlier than documented in the North Atlantic, while neither model is able to produce the reconstructed early Southern Ocean and Antarctic warming. Our results highlight the importance of producing a sequence of time slices rather than one single time slice averaging the LIG climate conditions.

Water temperature reconstructions based on paired Uk'37 and Tex86 records, and paired surface and subsurface Mg/Ca records

This dataset contains the collection of available published paired Uk'37 and Tex86 records spanning multi-millennial to multi-million year time scales, as well as a collection of Mg/Ca-derived temperatures measured in parallel on surface and subsurface dwelling foraminifera, both used in the analyses of Ho and Laepple, Nature Geoscience 2016. As the signal-to-noise ratios of proxy-derived Holocene temperatures are relatively low, we selected records that contain at least the last deglaciation (oldest sample >18kyr BP).

Mg/Ca and Sr/Ca ratios of planktonic foraminifera from North Atlantic surface sediments

Reliable temperature estimates from both surface and subsurface ocean waters are needed to reconstruct past upper water column temperature gradients and past oceanic heat content. This work examines the relationships between trace element ratios in fossil shells and seawater temperature for surface-dwelling foraminifera species, Globigerinoides ruber (white) and Globigerina bulloides, and deep-dwelling species, Globorotalia inflata, Globorotalia truncatulinoides (dextral and sinistral) and Pulleniatina obliquiloculata. Mg/Ca and Sr/Ca ratios in shells picked in 29 modern core tops from the North Atlantic Ocean are calibrated using calculated isotopic temperatures. Mg/Ca ratios on G. ruber and G. bulloides agree with published data and relationships. For deep-dwelling species, Mg/Ca calibration follows the equation Mg/Ca = 0.78 (±0.04) * exp (0.051 (±0.003) * T) with a significant correlation coefficient of R**2 = 0.74. Moreover, there is no significant difference between the different deep-dwellers analyzed. For the Sr/Ca ratio, the surface dwellers and P. obliquiloculata do not record any temperature dependence. For the Globorotalia species, the thermo dependence of Sr/Ca ratio can be described by a single linear relationship: Sr/Ca = (0.0182 (±0.001) * T) + 1.097 (±0.018), R**2 = 0.85. Temperature estimates with a 1 sigma error of ±2.0°C and ±1.3°C can be derived from the Mg/Ca and Sr/Ca ratios, respectively, as long as the Sr geochemistry in the ocean has been constant through time.

(Table 1) Temperature and foraminiferal Mg/Ca ratios of samples off Japan

We studied Mg/Ca in high-Mg, shallow-water benthic foraminifera in culture and in samples from natural environments, in order to evaluate the expression of latitudinal and seasonal temperature variability in Mg/Ca in their tests. We cultured Planoglabratella opercularis (d'Orbigny) and Quinqueloculina yabei Asano under controlled temperature (10°-25°C) and salinity (30-38) conditions. Both species show a linear correlation between Mg/Ca and temperature, but they differ in temperature sensitivity. Salinity does not significantly influence Mg/Ca. In the samples collected in nature, Mg/Ca and seawater temperatures are positively correlated, but there are more complexities than in the records for cultured specimens due to such factors as seasonal fluctuations in temperature. We conclude that Mg/Ca ratios in monospecific benthic foraminiferal samples may be used as a reliable temperature proxy, if the lifetime of the species is taken into account.

(Table 1) Radiocarbon dating of a lake sediment core, southern West Greenland

The dominant processes determining biological structure in lakes at millennial timescales are complex. In this study, we used a multi-proxy approach to determine the relative importance of in-lake versus indirect processes on the Holocene development of an oligotrophic lake in SW Greenland (66.99°N, 50.97°W). A 14C and 210Pb-dated sediment core covering approximately 8500 years BP was analyzed for organic-inorganic carbon content, pigments, diatoms, chironomids, cladocerans, and stable isotopes (d13C, d18O). Relationships among the different proxies and a number of independent controlling variables (Holocene temperature, an isotope-inferred cooling period, and immigration of Betula nana into the catchment) were explored using redundancy analysis (RDA) independent of time. The main ecological trajectories in the lake biota were captured by ordination first axis sample scores (18-32% variance explained). The importance of the arrival of Betula (ca. 6500 years BP) into the catchment was indicated by a series of partial-constrained ordinations, uniquely explaining 12-17% of the variance in chironomids and up to 9% in pigments. Climate influences on lake biota were strongest during a short-lived cooling period (identified by altered stable isotopes) early in the development of the lake when all proxies changed rapidly, although only chironomids had a unique component (8% in a partial-RDA) explained by the cooling event. Holocene climate explained less variance than either catchment changes or biotic relationships. The sediment record at this site indicates the importance of catchment factors for lake development, the complexity of community trends even in relatively simple systems (invertebrates are the top predators in the lake) and the challenges of deriving palaeoclimate inferences from sediment records in low-Arctic freshwater lakes.

Mg/Ca and stable oxygen isotope ratios of Uvigerina spp. from the Southern Ocean

The sensitivity to temperature of Mg/Ca ratios in the shallow-infaunal benthic foraminifera Uvigerina spp. has been assessed. Core-top calibrations over ~1-20 °C show a range in sensitivity of 0.065-0.084 mmol/mol/°C but few data are available spanning the temperature range anticipated in deep-sea records over glacial-interglacial cycles. In contrast to epibenthic foraminiferal species, carbonate ion saturation appears not to affect Mg/Ca significantly. A method based on estimating the ratio of the temperature sensitivity of foraminiferal Mg/Ca to that of d18Ocalcite shows that sensitivity for Mg/Ca at the high end of the observed core-top range (~0.1 mmol/mol/°C) is required for consistency with LGM-Holocene differences in each property as constrained by independent proxy data. This is supported by a Mg/Ca record for Uvigerina spp. generated for the Southern Ocean over the past 440,000 years from Ocean Drilling Program Site 1123 (Chatham Rise, New Zealand). The record shows variability that correlates with climate oscillations. The LGM deep ocean temperature derived from the Mg/Ca record is -1.1 ± 0.3 °C. Transformation to temperature allows estimates to be made of changes in bottom water temperature and seawater d18O and comparison made with literature records. Analysis reveals a ~2.5-kyr lead in the record of temperature over calcite d18O and a longer lead over seawater d18O. This is a reflection of larger phase offsets at eccentricity periods; phase offsets at tilt and precession are within error zero.

Mg/Ca ratio of Globorotalia tumida and subsurface temperature reconstruction for ODP Site 130-806

The tropical Pacific thermocline strength, depth, and tilt are critical to tropical mean state and variability. During the early Pliocene (~3.5 to 4.5 Ma), the Eastern Equatorial Pacific (EEP) thermocline was deeper and the cold tongue was warmer than today, which resulted in an mean state with a reduced zonal sea surface temperature gradient, or El Padre. However, it is unclear whether the deep thermocline was a local feature of the EEP or a basin-wide condition with global implications. Our measurements of Mg/Ca of Globorotalia tumida in a western equatorial Pacific site indicate Pliocene subsurface temperatures warmer than today; thus, El Padre included a basin-wide thermocline that was relatively warm, deep, and weakly tilted. At ~4 Ma, thermocline steepening was coupled to cooling of the cold tongue. Since ~4 Ma, the basin-wide thermocline cooled/shoaled gradually, with implications for thermocline feedbacks in tropical dynamics and the interpretation of TEX86-derived temperatures.

Sr/Ca and Cd/Ca ratios of benthic foraminifera in surface deep-sea sediments

Measurements of Sr/Ca of benthic foraminifera show a linear decrease with water depth which is superimposed upon significant variability identified by analyses of individual foraminifera. New data for Cd/Ca support previous work in defining a contrast between waters shallower and deeper than ~2500 m. Measured element partition coefficients in foraminiferal calcium carbonate relative to sea water (D) have been described by means of a one-box model in which elements are extracted by Rayleigh distillation from a biomineralization reservoir that serves for calcification with a constant fractionation factor (alpha), such that D = (1 - f**alpha)/(l - f), where f is the proportion of Ca remaining after precipitation. A modification to the model recognises differences in element speciation. The model is consistent with differences between D[Sr], D[Ba], and D[Cd] in benthic but not planktonic foraminifera. Depth variations in D for Sr and Ba are consistent with the model, as are differences in depth variation of D[Cd] in calcitic and aragonitic benthic foraminifera. The shallower depth variations may reflect increasing calcification rates with increasing water depth to an optimum of about 2500 m. Observations of unusually lower DCd for some deep waters, not accompanied by similar [Sr], or D[Ba] may be because of dissolution or a calcification response to a lower carbonate saturation state.

Great Barrier Reef coral element/Ca and stable isotope data and U-Th ages from IODP Expedition 325

Tropical south-western Pacific temperatures are of vital importance to the Great Barrier Reef (GBR), but the role of sea surface temperatures (SSTs) in the growth of the GBR since the Last Glacial Maximum remains largely unknown. Here we present records of Sr/Ca and d18O for Last Glacial Maximum and deglacial corals that show a considerably steeper meridional SST gradient than the present day in the central GBR. We find a 1-2 °C larger temperature decrease between 17° and 20°S about 20,000 to 13,000 years ago. The result is best explained by the northward expansion of cooler subtropical waters due to a weakening of the South Pacific gyre and East Australian Current. Our findings indicate that the GBR experienced substantial meridional temperature change during the last deglaciation, and serve to explain anomalous deglacial drying of northeastern Australia. Overall, the GBR developed through significant SST change and may be more resilient than previously thought.

Potassium and argon concentrations and age determination of ODP Site 125-782 and 125-786

K-Ar whole-rock ages have been obtained for 30 samples from Sites 782 and 786, Ocean Drilling Program Leg 125 in the Izu-Bonin (Ogasawara) forearc region. They form a trimodal spread of ages between 9 Ma and 44 Ma and are, with a few exceptions, consistent with the inferred lithostratigraphy. The ages have been interpreted in terms of at least two distinct episodes of magmatic and/or hydrothermal activity. A group of ten samples, including the lava flows, gave an isochron age of 41.3 ± 0.5 Ma (middle-late Eocene). This is thought to represent the age of the principal magmatic development of the volcanic forearc basement, and is comparable to published ages on equivalent rocks from other parts of the forearc basement high (e.g., the Ogasawara Islands). It may be significant that this age is slightly younger than the timing of major plate reorganization in the Western Pacific at about 43 Ma. This was followed by a minor episode of intrusive magmatism at 34.6 ± 0.7 Ma (early Oligocene) which appears to have reset the ages of some of the earlier units. This event probably corresponds to the initiation of rifting of the "proto-arc" to form the Parece Vela Basin. Boninitic samples were erupted during both episodes of magmatism, the earlier being of low-Ca boninite type and the later being of medium- and high-Ca types. It is also possible that a third episode of intrusive magmatism affected the Izu-Bonin forearc region at both Sites 782 and 786 at about 17 Ma. This would be consistent with magmatic activity elsewhere in the region during the Miocene, associated with the end of active spreading in the Parece Vela Basin and the start of arc activity in the West Mariana Ridge.

Planktonic foraminiferal Cd/Ca ratios

Records of Cd/Ca in planktonic foraminiferal calcite of Globigerinoides bulloides in cores from the Subantarctic region of the Southern Ocean show large glacial-interglacial variations with lower Cd/Ca (by 0.06-0.10 µmol/mol) at glacial times. Interpretation of these records in terms of lower dissolved phosphate and inferred higher glacial nutrient utilization has significant implications for glacial atmospheric carbon dioxide (pCO2) draw-down. However, box core-top data for G. bulloides in the North Atlantic suggest that the incorporation of Cd into planktonic foraminifera relative to seawater (DCd) is temperature sensitive (DCd=0.637 exp 0.15T). When the Subantarctic planktonic Cd/Ca records are corrected for this temperature dependence, they show little or no glacial-interglacial diferences. If, as seems likely, this observation can be interpreted to indicate a minimal change (< 0.5 µmol/kg) in surface water phosphate concentrations, then the explanation for lowered glacial pCO2 must be looked for elsewhere.