Auxiliary material and basic geochemical and benthic foraminiferal stable isotope data for the interval bracketing Eocene Thermal Maximum 2 (ETM2) at DSDP Sites 401 and 550 in the northeastern Atlantic Ocean

Ever since its discovery, Eocene Thermal Maximum 2 (ETM2; ~53.7 Ma) has been considered as one of the "little brothers" of the Paleocene-Eocene Thermal Maximum (PETM; ~56 Ma) as it displays similar characteristics including abrupt warming, ocean acidification, and biotic shifts. One of the remaining key questions is what effect these lesser climate perturbations had on ocean circulation and ventilation and, ultimately, biotic disruptions. Here we characterize ETM2 sections of the NE Atlantic (Deep Sea Drilling Project Sites 401 and 550) using multispecies benthic foraminiferal stable isotopes, grain size analysis, XRF core scanning, and carbonate content. The magnitude of the carbon isotope excursion (0.85-1.10 per mil) and bottom water warming (2-2.5°C) during ETM2 seems slightly smaller than in South Atlantic records. The comparison of the lateral d13C gradient between the North and South Atlantic reveals that a transient circulation switch took place during ETM2, a similar pattern as observed for the PETM. New grain size and published faunal data support this hypothesis by indicating a reduction in deepwater current velocity. Following ETM2, we record a distinct intensification of bottom water currents influencing Atlantic carbonate accumulation and biotic communities, while a dramatic and persistent clay reduction hints at a weakening of the regional hydrological cycle. Our findings highlight the similarities and differences between the PETM and ETM2. Moreover, the heterogeneity of hyperthermal expression emphasizes the need to specifically characterize each hyperthermal event and its background conditions to minimalize artifacts in global climate and carbonate burial models for the early Paleogene.

(Table 1) Trace element/Ca ratios of Orbulina universa shells from late Neogene Mediterranean samples

A Mediterranean composite sedimentary record was analyzed for Ba/Ca ratios on carbonate shells of Orbulina universa planktonic foraminifer (Ba/Ca)carb providing the opportunity to study and assess the extent of freshwater inputs on the basin and possible impacts on its dynamics during the Tortonian to Recent period. A number of scanning electron microscope analyses and auxiliary trace element measurements (Mn, Sr, and Mg), obtained from the same samples, exclude important diagenetic effects on the studied biogenic carbonates and corroborate the reliability of (Ba/Ca)carb ratios in foraminifera calcite as indicators of seawater source components during the studied interval. A long-term trend with (Ba/Ca)carb values shifting from ~7 to 3 µmol/mol from the base of the Tortonian to the top of the Messinian is observed. The interval of the late Messinian salinity crisis, where biogenic carbonates are missing or strongly diagenized, represents a crucial passage not monitored in our record. At the base of the Pliocene, up to about 4.7 Ma, the (Ba/Ca)carb record shows a decreasing trend from ~4 µmol/mol stabilizing itself to an about constant value of 0.9 ± 0.3 µmol/mol for the whole Plio-Pleistocene interval. These results suggest a dramatic change in the continental runoff values, up to ~3-16 times higher during part of the late Neogene (Tortonian-early Pliocene), with a possible profound modification in the physical dynamics of the Mediterranean basin. First-order mass balance equations used to estimate barium and salinity budgets in the Mediterranean Sea during the late Miocene-early Pliocene interval support the hypothesis of an active connection of the basin with the Paratethys region and of a definitive restriction at the base of the Pliocene after about 0.7 Ma from the well-known Messinian Lagomare phase. They also open intriguing scenarios on possible circulation shifts during the Neogene.

(Table 2) Stratigraphic position of midpoints of reddish and beige layers of lithologic cycles in the Loulja-A Section and their astronomical ages after tuning to the ETP and Summer Insolation Curves of the La2004 solution

An integrated high-resolution stratigraphy and orbital tuning is presented for the Loulja sections located in the Bou Regreg area on the Atlantic side of Morocco. The sections constitute the upward continuation of the upper Messinian Ain el Beida section and contain a well-exposed, continuous record of the interval straddling the Miocene-Pliocene (M-P) boundary. The older Loulja-A section, which covers the interval from ~5.59 to 5.12 Ma, reveals a dominantly precession-controlled color cyclicity that allows for a straightforward orbital tuning of the boundary interval and for detailed cyclostratigraphic correlations to the Mediterranean; the high-resolution and high-quality benthic isotope record allows us to trace the dominantly obliquity-controlled glacial history. Our results reveal that the M-P boundary coincides with a minor, partly precession-related shift to lighter "interglacial" values in d18O. This shift and hence the M-P boundary may not correlate with isotope stage TG5, as previously thought, but with an extra (weak) obliquity-controlled cycle between TG7 and TG5. Consequently, the M-P boundary and basal Pliocene flooding of the Mediterranean following the Messinian salinity crisis are not associated with a major deglaciation and glacio-eustatic sea level rise, indicating that other factors, such as tectonics, must have played a fundamental role. On the other hand, the onset of the Upper Evaporites in the Mediterranean marked by hyposaline conditions coincides with the major deglaciation step between marine isotope stage TG12 and TG11, suggesting that the associated sea level rise is at least partly responsible for the apparent onset of intermittently restricted marine conditions following the main desiccation phase. Finally, the Loulja-A section would represent an excellent auxiliary boundary stratotype for the M-P boundary as formally defined at the base of the Trubi marls in the Eraclea Minoa section on Sicily.

Secchi depth maps for the Baltic Sea, 2003-2012, with links to GeoTIFs

Secchi depth is a measure of water transparency. In the Baltic Sea region, Secchi depth maps are used to assess eutrophication and as input for habitat models. Due to their spatial and temporal coverage, satellite data would be the most suitable data source for such maps. But the Baltic Sea's optical properties are so different from the open ocean that globally calibrated standard models suffer from large errors. Regional predictive models that take the Baltic Sea's special optical properties into account are thus needed. This paper tests how accurately generalized linear models (GLMs) and generalized additive models (GAMs) with MODIS/Aqua and auxiliary data as inputs can predict Secchi depth at a regional scale. It uses cross-validation to test the prediction accuracy of hundreds of GAMs and GLMs with up to 5 input variables. A GAM with 3 input variables (chlorophyll a, remote sensing reflectance at 678 nm, and long-term mean salinity) made the most accurate predictions. Tested against field observations not used for model selection and calibration, the best model's mean absolute error (MAE) for daily predictions was 1.07 m (22%), more than 50% lower than for other publicly available Baltic Sea Secchi depth maps. The MAE for predicting monthly averages was 0.86 m (15%). Thus, the proposed model selection process was able to find a regional model with good prediction accuracy. It could be useful to find predictive models for environmental variables other than Secchi depth, using data from other satellite sensors, and for other regions where non-standard remote sensing models are needed for prediction and mapping. Annual and monthly mean Secchi depth maps for 2003-2012 come with this paper as Supplementary materials.

Cenozoic oxygen isotopic values for benthic foraminifera from DSDP and ODP low and high latitude marine sediment cores

The climate during the Cenozoic era changed in several steps from ice-free poles and warm conditions to ice-covered poles and cold conditions. Since the 1950s, a body of information on ice volume and temperature changes has been built up predominantly on the basis of measurements of the oxygen isotopic composition of shells of benthic foraminifera collected from marine sediment cores. The statistical methodology of time series analysis has also evolved, allowing more information to be extracted from these records. Here we provide a comprehensive view of Cenozoic climate evolution by means of a coherent and systematic application of time series analytical tools to each record from a compilation spanning the interval from 4 to 61 Myr ago. We quantitatively describe several prominent features of the oxygen isotope record, taking into account the various sources of uncertainty (including measurement, proxy noise, and dating errors). The estimated transition times and amplitudes allow us to assess causal climatological-tectonic influences on the following known features of the Cenozoic oxygen isotopic record: Paleocene-Eocene Thermal Maximum, Eocene-Oligocene Transition, Oligocene-Miocene Boundary, and the Middle Miocene Climate Optimum. We further describe and causally interpret the following features: Paleocene-Eocene warming trend, the two-step, long-term Eocene cooling, and the changes within the most recent interval (Miocene-Pliocene). We review the scope and methods of constructing Cenozoic stacks of benthic oxygen isotope records and present two new latitudinal stacks, which capture besides global ice volume also bottom water temperatures at low (less than 30°) and high latitudes. This review concludes with an identification of future directions for data collection, statistical method development, and climate modeling.

Plant Functional Type (PFT) map over Siberia derived from GlobCover 2005 dataset, link to dataset in NetCDF format

High-latitude ecosystems play an important role in the global carbon cycle and in regulating the climate system and are presently undergoing rapid environmental change. Accurate land cover data sets are required to both document these changes as well as to provide land-surface information for benchmarking and initializing Earth system models. Earth system models also require specific land cover classification systems based on plant functional types (PFTs), rather than species or ecosystems, and so post-processing of existing land cover data is often required. This study compares over Siberia, multiple land cover data sets against one another and with auxiliary data to identify key uncertainties that contribute to variability in PFT classifications that would introduce errors in Earth system modeling. Land cover classification systems from GLC 2000, GlobCover 2005 and 2009, and MODIS collections 5 and 5.1 are first aggregated to a common legend, and then compared to high-resolution land cover classification systems, vegetation continuous fields (MODIS VCFs) and satellite-derived tree heights (to discriminate against sparse, shrub, and forest vegetation). The GlobCover data set, with a lower threshold for tree cover and taller tree heights and a better spatial resolution, tends to have better distributions of tree cover compared to high-resolution data. It has therefore been chosen to build new PFT maps for the ORCHIDEE land surface model at 1 km scale. Compared to the original PFT data set, the new PFT maps based on GlobCover 2005 and an updated cross-walking approach mainly differ in the characterization of forests and degree of tree cover. The partition of grasslands and bare soils now appears more realistic compared with ground truth data. This new vegetation map provides a framework for further development of new PFTs in the ORCHIDEE model like shrubs, lichens and mosses, to represent the water and carbon cycles in northern latitudes better. Updated land cover data sets are critical for improving and maintaining the relevance of Earth system models for assessing climate and human impacts on biogeochemistry and biophysics.

Alkenone-derived UK'37 data and sea surface temperatures (SST) of sediment cores from the Gulf of Taranto (Italy)

Auxiliary data include one file with alkenone-derived UK'37 data and sea surface temperatures (SST). On these data Figs. 7 and 8 of the manuscript are based. The SST are derived from UK'37 by using the transfer function: SST = 29.876 UK'37 - 1.334 of Conte et al. (2006). The data are against the ages (in A.D.) of samples derived from cores GT91-1 (39[deg]59'23"N, 17[deg]45'25"E), GT89-3 and GT90-3 (both 39[deg]45'43"N, 17[deg]53'55"E ). Also included are composite records for UK'37 and SST. For creating the composite records, GT-89-3 was taken as reference core. In the overlapping period the GT89-3 data seem in general lower than the GT91-1 data. To accommodate for this in the composite record, the average difference (0.0343 UK'37 units; equivalent to 1.023 [deg]C) was subtracted from the GT91-1 record. Hereafter, for each depth in the overlapping interval the respective values (UK'37 or SST) of GT89-3 and GT91-1 were averaged. We have also averaged with 16 additional alkenone measurements, from 1793 to 1851, performed in the GT90-3 core.