CaCO3, foraminifer fragmentation, benthic carbon isotopes, coarse lithic counts and radiogenic isotopes (Sr, Nd, Pb) of the Pliocene and earliest Pleistocene terrigenous component (3.3-2.4 Ma) of IODP Site 306-U1313

We present Plio-Pleistocene records of sediment color, %CaCO3, foraminifer fragmentation, benthic carbon isotopes (d13C) and radiogenic isotopes (Sr, Nd, Pb) of the terrigenous component from IODP Site U1313, a reoccupation of benchmark subtropical North Atlantic Ocean DSDP Site 607. We show that (inter)glacial cycles in sediment color and %CaCO3 pre-date major northern hemisphere glaciation and are unambiguously and consistently correlated to benthic oxygen isotopes back to 3.3 million years ago (Ma) and intermittently so probably back to the Miocene/Pliocene boundary. We show these lithological cycles to be driven by enhanced glacial fluxes of terrigenous material (eolian dust), not carbonate dissolution (the classic interpretation). Our radiogenic isotope data indicate a North American source for this dust (~3.3-2.4 Ma) in keeping with the interpreted source of terrestrial plant wax-derived biomarkers deposited at Site U1313. Yet our data indicate a mid latitude provenance regardless of (inter)glacial state, a finding that is inconsistent with the biomarker-inferred importance of glaciogenic mechanisms of dust production and transport. Moreover, we find that the relation between the biomarker and lithogenic components of dust accumulation is distinctly non-linear. Both records show a jump in glacial rates of accumulation from Marine Isotope Stage, MIS, G6 (2.72 Ma) onwards but the amplitude of this signal is about 3-8 times greater for biomarkers than for dust and particularly extreme during MIS 100 (2.52 Ma). We conclude that North America shifted abruptly to a distinctly more arid glacial regime from MIS G6, but major shifts in glacial North American vegetation biomes and regional wind fields (exacerbated by the growth of a large Laurentide Ice Sheet during MIS 100) likely explain amplification of this signal in the biomarker records. Our findings are consistent with wetter-than-modern reconstructions of North American continental climate under the warm high CO2 conditions of the Early Pliocene but contrast with most model predictions for the response of the hydrological cycle to anthropogenic warming over the coming 50 years (poleward expansion of the subtropical dry zones).

Sedimentological and benthic foraminiferal data pertaining to ODP Holes 208-1262A and 208-1263C

Eocene Thermal Maximum 2 (ETM2) occurred ~1.8 Myr after the Paleocene Eocene Thermal Maximum (PETM) and, like the PETM, was characterized by a negative carbon isotope excursion coupled with warming. We combined benthic foraminiferal and sedimentological records for Southeast Atlantic Sites 1263 (1500 m paleodepth) and 1262 (3600 m paleodepth) to show that benthic foraminiferal diversity and accumulation rates declined more precipitously and severely at the shallower site during peak ETM2. The sites are in close proximity, so differences in surface productivity cannot have caused this differential effect. Instead, on the basis of an analysis of climate modelling experiments, we infer that changes in ocean circulation pattern across ETM2 may have resulted in more pronounced warming at intermediate depths (Site 1263). The effects of more pronounced warming include increased metabolic rates, leading to a decrease in effective food supply and increased deoxygenation, thus potentially explaining the more severe benthic impacts at Site 1263. In response to more severe benthic disturbance, bioturbation may have decreased at Site 1263 as compared to Site 1262, hence differentially affecting the bulk carbonate record. We use a sediment-enabled Earth system model to test whether a reduction in bioturbation and/or the likely reduced carbonate saturation of more poorly ventilated waters can explain the more extreme excursion in bulk d13C and sharper transition in wt% CaCO3 at Site 1263. We find that both enhanced acidification and reduced bioturbation during peak ELMO conditions are needed to account for the observed features. Our combined ecological and modelling analysis illustrates the potential role of ocean circulation changes in amplifying local environmental changes and driving temporary, but drastic, loss of benthic biodiversity and abundance.

Heat flow measurements on METEOR cruise M86/5

Mud volcanoes (MV) are sources of mass and energy, transported from deeper levels of the sediment pile to the surface. Together with fluid and gas, thermal energy is emitted through these structures. Therefore heat flow determination is a sensible tool to detect and quantify the amount of convective flow. In the Gulf of Cadiz several mud volcanoes can be found along major tectonic lines (SWIM faults). We employ geothermal measurements to observe the activity of mud volcanoes and possible leakage at the faults apart from pronounced structures.

Tomographic measurements of O. umbonatus and N. truempyi during the PETM and ETM-2 at ODP Sites 208-1262 and 208-1263

Predicting the impact of ongoing anthropogenic CO2 emissions on calcifying marine organisms is complex, owing to the synergy between direct changes (acidification) and indirect changes through climate change (e.g., warming, changes in ocean circulation, and deoxygenation). Laboratory experiments, particularly on longer-lived organisms, tend to be too short to reveal the potential of organisms to acclimatize, adapt, or evolve and usually do not incorporate multiple stressors. We studied two examples of rapid carbon release in the geological record, Eocene Thermal Maximum 2 (~53.2 Ma) and the Paleocene Eocene Thermal Maximum (PETM, ~55.5 Ma), the best analogs over the last 65 Ma for future ocean acidification related to high atmospheric CO2 levels. We use benthic foraminifers, which suffered severe extinction during the PETM, as a model group. Using synchrotron radiation X-ray tomographic microscopy, we reconstruct the calcification response of survivor species and find, contrary to expectations, that calcification significantly increased during the PETM. In contrast, there was no significant response to the smaller Eocene Thermal Maximum 2, which was associated with a minor change in diversity only. These observations suggest that there is a response threshold for extinction and calcification response, while highlighting the utility of the geological record in helping constrain the sensitivity of biotic response to environmental change.

Geochemistry and isotopes at DSDP Site 68-503

Nodules occur in the siliceous calcareous ooze and siliceous marl at Site 503 in the eastern equatorial Pacific. They are present below a depth of about 11 meters throughout the green-colored reduced part of the section down to 228 meters, although they are most abundant between 30 and 85 meters. They are cylindrical or barrel-shaped, up to 70 mm long, and usually have an axial channel through them or are hollow. They appear to have formed around and/or within burrows. XRD studies and microprobe analyses show that they are homogeneous and consist of calcian rhododrosite and minor calcite; Mn is present to the extent of about 30%. Isotopic analyses of the carbonate give carbon values which range from -1.2 per mil to -3.8 per mil, and oxygen isotope compositions vary from +4.0 per mil to +6.0 per mil. These values are different from those for marine-derived carbonates as exemplified by the soft sediment filling of a burrow: d13C, -0.26 per mil; d18O, +1.05 per mil. The carbon isotope data indicate that carbonate derived (possibly indirectly) from seawater was mixed with some produced by organic diagenesis to form the nodules. The d18O values suggest that although they formed near the sediment surface, some modification or the introduction of additional diagenetic carbonate occurred during burial.

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.