Methane concentration in sediment cores during METEOR cruise M92

We studied the concurrence of methanogenesis and sulfate reduction in surface sediments (0-25 cm below sea floor, cmbsf) at six stations (70, 145, 253, 407, 770 and 1024 m) along the Peruvian margin (12° S). This oceanographic region is characterized by high carbon export to the seafloor, creating an extensive oxygen minimum zone (OMZ) on the shelf, both factors that could favor surface methanogenesis. Sediments sampled along the depth transect traversed areas of anoxic and oxic conditions in the bottom-near water. Net methane production (batch incubations) and sulfate reduction (35S-sulfate radiotracer incubation) were determined in the upper 0-25 cmbsf of multicorer cores from all stations, while deep hydrogenotrophic methanogenesis (> 30 cmbsf, 14C-bicarbonate radiotracer incubation) was determined in two gravity cores at selected sites (78 and 407 m). Furthermore, stimulation (methanol addition) and inhibition (molybdate addition) experiments were carried out to investigate the relationship between sulfate reduction and methanogenesis. Highest rates of methanogenesis and sulfate reduction in the surface sediments, integrated over 0-25 cmbsf, were observed on the shelf (70-253 m, 0.06-0.1 and 0.5-4.7 mmol m-2 d-1, respectively), while lowest rates were discovered at the deepest site (1024 m, 0.03 and 0.2 mmol m-2 d-1, respectively). The addition of methanol resulted in significantly higher surface methanogenesis activity, suggesting that the process was mostly based on non-competitive substrates, i.e., substrates not used by sulfate reducers. In the deeper sediment horizons, where competition was probably relieved due to the decline of sulfate, the usage of competitive substrates was confirmed by the detection of hydrogenotrophic activity in the sulfate-depleted zone at the shallow shelf station (70 m). Surface methanogenesis appeared to be correlated to the availability of labile organic matter (C / N ratio) and organic carbon degradation (DIC production), both of which support the supply of methanogenic substrates. A negative correlation of methanogenesis rates with dissolved oxygen in the bottom-near water was not obvious, however, anoxic conditions within the OMZ might be advantageous for methanogenic organisms at the sediment-water interface. Our results revealed a high relevance of surface methanogenesis on the shelf, where the ratio between surface to deep (below sulfate penetration) methanogenic activity ranged between 0.13 and 105. In addition, methane concentration profiles indicate a partial release of surface methane into the water column as well as a partial consumption of methane by anaerobic methane oxidation (AOM) in the surface sediment. The present study suggests that surface methanogenesis might play a greater role in benthic methane budgeting than previously thought, especially for fueling AOM above the sulfate-methane transition zone.

Additional CaCO3 measurements from IODP Site 306-U1313

We present further %CaCO3 data from Site U1313 across the Pliocene-Pleistocene intensification of Northern Hemisphere glaciation. This data was measured on the U1313 secondary splice. We also present tie points between the primary and secondary splice for this interval based on graphical tuning of L* (sediment lightness).

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.

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.

Deep water formation in the North Pacific and deglacial CO2 rise

Deep water formation in the North Atlantic and Southern Ocean is widely thought to influence deglacial CO2 rise and climate change; here we suggest that deep water formation in the North Pacific may also play an important role. We present paired radiocarbon and boron isotope data from foraminifera from sediment core MD02-2489 at 3640 m in the North East Pacific. These show a pronounced excursion during Heinrich Stadial 1, with benthic-planktic radiocarbon offsets dropping to ~350 years, accompanied by a decrease in benthic d11B. We suggest this is driven by the onset of deep convection in the North Pacific, which mixes young shallow waters to depth, old deep waters to the surface, and low-pH water from intermediate depths into the deep ocean. This deep water formation event was likely driven by an increase in surface salinity, due to subdued atmospheric/monsoonal freshwater flux during Heinrich Stadial 1. The ability of North Pacific Deep Water (NPDW) formation to explain the excursions seen in our data is demonstrated in a series of experiments with an intermediate complexity Earth system model. These experiments also show that breakdown of stratification in the North Pacific leads to a rapid ~30 ppm increase in atmospheric CO2, along with decreases in atmospheric d13C and D14C, consistent with observations of the early deglaciation. Our inference of deep water formation is based mainly on results from a single sediment core, and our boron isotope data are unavoidably sparse in the key HS1 interval, so this hypothesis merits further testing. However we note that there is independent support for breakdown of stratification in shallower waters during this period, including a minimum in d15N, younging in intermediate water 14C, and regional warming. We also re-evaluate deglacial changes in North Pacific productivity and carbonate preservation in light of our new data, and suggest that the regional pulse of export production observed during the Bølling-Allerød is promoted by relatively stratified conditions, with increased light availability and a shallow, potent nutricline. Overall, our work highlights the potential of NPDW formation to play a significant and hitherto unrealized role in deglacial climate change and CO2 rise.

Geochemistry of Lesser Antilles lavas

We present new U-series disequilibrium and radiogenic isotope data for 7 mafic lavas from the Lesser Antilles arc. These are combined with published data in an internally consistent model that quantitatively estimates the amount of sediment and fluid added to the source of the Lesser Antilles arc system. Some lavas form an array consistent with bulk sediment addition (0.2-2%) whereas others appear to require addition of 0.4-2% sediment melt, particularly in the south of the arc. Evidence for both bulk sediment and sediment melt addition can be found within both the northern and central sections of the arc suggesting a thermal structure whereby the upper portions of the subducted sediment pile lie close to their solidus beneath much of the arc. Addition of up to 5% fluid derived from altered oceanic crust to these sediment enriched mantle wedge source regions can simulate the majority of the lavas on a plot of 207Pb/204Pb versus Ce/Pb. By taking into account the range in calculated wedge compositions and allowing for some mobility of Th in the fluid, the same model can also account for much of the observed range of U-Th-Ra disequilibria, especially if the eclogitic residue contains trace amounts of rutile. The implication of this more complex model is that the time scales for fluid addition and differentiation could be significantly shorter than those estimated in some previous studies.

Geochemistry of sediment cores from the Rainbow vent site

A geochemical investigation was carried out on two sediment cores collected at 2 and 5 km from the Rainbow hydrothermal vent site. Bulk sediment compositions indicate that these cores record clear enrichments in Fe, Cu, Mn, V, P and As from hydrothermal plume fallout (Cave et al., 2002, doi:10.1016/S0016-7037(02)00823-2). Sequential dissolution of the bulk sediments has been used to discriminate between a leach (biogenic and oxy-hydroxide) component and a residual phase (detrital and sulphide/sulphate fractions). Major element data (Al, Fe, Ti, Mn, Mg, Ca, Si and index%) reveal that the hydrothermal input, as recorded in the leach phase, is much stronger than apparent from bulk sediment analyses alone. REE patterns for the leach phase record contributions from both biogenic carbonate (mimicking seawater REE patterns) and hydrothermal oxy-hydroxides, with the latter exhibiting positive Eu anomalies (hydrothermal derived) and negative Ce anomalies (seawater derived). Based on major element and REE data, the residue contains contributions from aeolian dust input, local MORB material and a hydrothermal component. Ternary REE mixing calculations indicate that most of the REE within the residual fraction (~80%) is derived from hydrothermal material, while detrital contributions to the REE budget, as deep-sea clay and volcanic debris, are <20%. By combining bulk and REE data for the various end-member components of the residue, we calculate that the chemical composition of the residue hydrothermal end-member is high in Ca (6-15%) and with a Nd/Sr ratio of 0.004. These characteristics indicate the presence of low-solubility hydrothermal sulphate (rather than sulphide) material within the residue component of Rainbow hydrothermal sediments.