Physical properties, pore water geochemistry, 210Pb-dating of sediment cores GeoB16426-1, GeoB16427-1, GeoB16429-1

We present differential bathymetry and sediment core data from the Japan Trench, sampled after the 2011 Tohoku-Oki (offshore Japan) earthquake to document that prominent bathymetric and structural changes along the trench axis relate to a large (~27.7 km**2) slump in the trench. Transient geochemical signals in the slump deposit and analysis of diffusive re-equilibration of disturbed SO4**2- profiles over time constrain the triggering of the slump to the 2011 earthquake. We propose a causal link between earthquake slip to the trench and rotational slumping above a subducting horst structure. We conclude that the earthquake-triggered slump is a leading agent for accretion of trench sediments into the forearc and hypothesize that forward growth of the prism and seaward advance of the deformation front by more than 2 km can occur, episodically, during a single-event, large mega-thrust earthquake.

Micro-facies and µXRF element scanning data of Dead Sea sediment cores DSEn and 5017-1

Laminated lake sediments from the Dead Sea basin provide high-resolution records of climatic variability in the eastern Mediterranean region, which is especially sensitive to changing climatic conditions. In this study, we aim on detailed reconstruction of climatic fluctuations and related changes in the frequency of flood and dust deposition events at ca. 3300 and especially at 2800 cal. yr BP from high-resolution sediment records of the Dead Sea basin. A ca. 4-m-thick, mostly varved sediment section from the western margin of the Dead Sea (DSEn - Ein Gedi profile) was analysed and correlated to the new International Continental Scientific Drilling Program (ICDP) Dead Sea Deep Drilling Project core 5017-1 from the deep basin. To detect even single event layers, we applied a multi-proxy approach of high-resolution microscopic thin section analyses, micro-X-ray fluorescence (µ-XRF) element scanning and magnetic susceptibility measurements, supported by grain size data and palynological analyses. Based on radiocarbon and varve dating, two pronounced dry periods were detected at ~3500-3300 and ~3000-2400 cal. yr BP which are differently expressed in the sediment records. In the shallow-water core (DSEn), the older dry period is characterised by a thick sand deposit, whereas the sedimentological change at 2800 cal. yr BP is less pronounced and characterised mainly by an enhanced frequency of coarse detrital layers interpreted as erosion events. In the 5017-1 deep-basin core, both dry periods are depicted by halite deposits. The onset of the younger dry period coincides with the Homeric Grand Solar Minimum at ca. 2800 cal. yr BP. Our results suggest that during this period, the Dead Sea region experienced an overall dry climate, superimposed by an increased occurrence of flash floods caused by a change in synoptic weather patterns.

Fig. S1. Loci of possible event locations based on differential ts - tp arrival times at the single functioning geophone

Geological storage of CO2 that has been captured at large, point source emitters represents a key potential method for reduction of anthropogenic greenhouse gas emissions. However, this technology will only be viable if it can be guaranteed that injected CO2 will remain trapped in the subsurface for thousands of years or more. A signi?cant issue for storage security is the geomechanical response of the reservoir. Concerns have been raised that geomechanical deformation induced by CO2 injection will create or reactivate fracture networks in the sealing caprocks, providing a pathway for CO2 leakage. In this paper, we examine three large-scale sites where CO2 is injected at rates of ab. 1 megatonne/y or more: Sleipner, Weyburn, and In Salah. We compare and contrast the observed geomechanical behavior of each site, with particular focus on the risks to storage security posed by geomechanical deformation. At Sleipner, the large, high-permeability storage aquifer has experienced little pore pressure increase over 15 y of injection, implying little possibility of geomechanical deformation. At Weyburn, 45 y of oil production has depleted pore pressures before increases associated with CO2 injection. The long history of the ?eld has led to complicated, sometimes nonintuitive geomechanical deformation. At In Salah, injection into the water leg of a gas reservoir has increased pore pressures, leading to uplift and substantial microseismic activity. The differences in the geomechanical responses of these sites emphasize the need for systematic geomechanical appraisal before injection in any potential storage site.

Global compilation of benthic data sets I

Approaches to quantify the organic carbon accumulation on a global scale generally do not consider the small-scale variability of sedimentary and oceanographic boundary conditions along continental margins. In this study, we present a new approach to regionalize the total organic carbon (TOC) content in surface sediments (<5 cm sediment depth). It is based on a compilation of more than 5500 single measurements from various sources. Global TOC distribution was determined by the application of a combined qualitative and quantitative-geostatistical method. Overall, 33 benthic TOC-based provinces were defined and used to process the global distribution pattern of the TOC content in surface sediments in a 1°x1° grid resolution. Regional dependencies of data points within each single province are expressed by modeled semi-variograms. Measured and estimated TOC values show good correlation, emphasizing the reasonable applicability of the method. The accumulation of organic carbon in marine surface sediments is a key parameter in the control of mineralization processes and the material exchange between the sediment and the ocean water. Our approach will help to improve global budgets of nutrient and carbon cycles.

Geochemistry and stratigraphy on sediment cores from Leg198 and Leg 208

The first complete cyclic sedimentary successions for the early Paleogene from drilling multiple holes have been retrieved during two ODP expeditions: Leg 198 (Shatsky Rise, NW Pacific Ocean) and Leg 208 (Walvis Ridge, SE Atlantic Ocean). These new records allow us to construct a comprehensive astronomically calibrated stratigraphic framework with an unprecedented accuracy for both the Atlantic and the Pacific Oceans covering the entire Paleocene epoch based on the identification of the stable long-eccentricity cycle (405-kyr). High resolution X-ray fluorescence (XRF) core scanner and non-destructive core logging data from Sites 1209 through1211 (Leg 198) and Sites 1262, 1267 (Leg 208) are the basis for such a robust chronostratigraphy. Former investigated marine (ODP Sites 1001 and 1051) and land-based (e.g., Zumaia) sections have been integrated as well. The high-fidelity chronology is the prerequisite for deciphering mechanisms in relation to prominent transient climatic events as well as completely new insights into Greenhouse climate variability in the early Paleogene. We demonstrate that the Paleocene epoch covers 24 long eccentricity cycles. We also show that no definite absolute age datums for the K/Pg boundary or the Paleocene - Eocene Thermal Maximum (PETM) can be provided by now, because of still existing uncertainties in orbital solutions and radiometric dating. However, we provide two options for tuning of the Paleocene which are only offset by 405-kyr. Our orbitally calibrated integrated Leg 208 magnetostratigraphy is used to revise the Geomagnetic Polarity Time Scale (GPTS) for Chron C29 to C25. We established a high-resolution calcareous nannofossil biostratigraphy for the South Atlantic which allows a much more detailed relative scaling of stages with biozones. The re-evaluation of the South Atlantic spreading rate model features higher frequent oscillations in spreading rates for magnetochron C28r, C27n, and C26n.

Mean dynamic topography and oceanographic parameters estimated from an inverse model and satellite geodesy, with link to model result in one single NetCDF file (392 MB), including the inverse data error covariance

Geostrophic surface velocities can be derived from the gradients of the mean dynamic topography-the difference between the mean sea surface and the geoid. Therefore, independently observed mean dynamic topography data are valuable input parameters and constraints for ocean circulation models. For a successful fit to observational dynamic topography data, not only the mean dynamic topography on the particular ocean model grid is required, but also information about its inverse covariance matrix. The calculation of the mean dynamic topography from satellite-based gravity field models and altimetric sea surface height measurements, however, is not straightforward. For this purpose, we previously developed an integrated approach to combining these two different observation groups in a consistent way without using the common filter approaches (Becker et al. in J Geodyn 59(60):99-110, 2012, doi:10.1016/j.jog.2011.07.0069; Becker in Konsistente Kombination von Schwerefeld, Altimetrie und hydrographischen Daten zur Modellierung der dynamischen Ozeantopographie, 2012, http://nbn-resolving.de/nbn:de:hbz:5n-29199). Within this combination method, the full spectral range of the observations is considered. Further, it allows the direct determination of the normal equations (i.e., the inverse of the error covariance matrix) of the mean dynamic topography on arbitrary grids, which is one of the requirements for ocean data assimilation. In this paper, we report progress through selection and improved processing of altimetric data sets. We focus on the preprocessing steps of along-track altimetry data from Jason-1 and Envisat to obtain a mean sea surface profile. During this procedure, a rigorous variance propagation is accomplished, so that, for the first time, the full covariance matrix of the mean sea surface is available. The combination of the mean profile and a combined GRACE/GOCE gravity field model yields a mean dynamic topography model for the North Atlantic Ocean that is characterized by a defined set of assumptions. We show that including the geodetically derived mean dynamic topography with the full error structure in a 3D stationary inverse ocean model improves modeled oceanographic features over previous estimates.

Analyses of shocked quartz from the K-P boundary

The precise cause and timing of the Cretaceous-Paleocene (K-P) mass extinction 65 Ma ago remains a matter of debate. Many advocate that the extinction was caused by a meteorite impact at Chicxulub, Mexico, and a number of potential kill-mechanisms have been proposed for this. Although we now have good constraints on the size of this impact and chemistry of the target rocks, estimates of its environmental consequences are hindered by a lack of knowledge about the obliquity of this impact. An oblique impact is likely to have been far more catastrophic than a sub-vertical one, because greater volumes of volatiles would have been released into the atmosphere. The principal purpose of this study was to characterize shocked quartz within distal K-P ejecta, to investigate whether the quartz distribution carried a signature of the direction and angle of impact. Our analyses show that the total number, maximum and average size of shocked quartz grains all decrease gradually with paleodistance from Chicxulub. We do not find particularly high abundances in Pacific sites relative to Atlantic and European sites, as has been previously reported, and the size-distribution around Chicxulub is relatively symmetric. Ejecta samples at any one site display features that are indicative of a wide range of shock pressures, but the mean degree of shock increases with paleodistance. These shock- and size-distributions are both consistent with the K-P layer having been formed by a single impact at Chicxulub. One site in the South Atlantic contains quartz indicating an anomalously high average shock degree, that may be indicative of an oblique impact with an uprange direction to the southeast +/- 45°. The apparent continuous coverage of proximal ejecta in this quadrant of the crater, however, suggests a relatively high impact angle of >45°. We conclude that some of the more extreme predictions of the environmental consequences of a low-angle impact at Chicxulub are probably not applicable.

Stable carbon and oxygen isotope compositions of benthic foraminifera N. truempyi samples and coarse fraction analyses from ODP Site 198-1209

We have compiled the first stratigraphically continuous high-resolution benthic foraminiferal stable isotope record for the Paleocene from a single site utilizing cores recovered at Pacific ODP Site 1209. The long-term trend in the benthic isotope record suggests a close coupling of volcanic CO2 input and deep sea warming. Over the short-term the record is characterized by slow excursions with a pronounced periodic beat related to the short (100-kyr) and long (405-kyr) eccentricity cycle. The phase relationship between the benthic isotope record and eccentricity is similar to patterns documented for the Oligocene and Miocene confirming the role of orbital forcing as the pace maker for paleoclimatic variability on Milankovitch time scales. In addition, the record documents an unusual transient warming of 2°C coeval with a 0.6 per mil carbon isotope excursion and a decrease in carbonate content at 61.75 Ma. This event, which bears some resemblance to Eocene hyperthermals, marks the onset of a long-term decline in d13C. The timing indicates it might be related to the initiation of volcanism along Greenland margin.

Physical oceaongraphy and hydrochemistry measured on water bottle samples during METEOR cruise M84/3 in 2001

Here we report on data from an oceanographic cruise on the German research vessel Meteor covering large parts of the Mediterranean Sea during spring of 2011. The main objective of this cruise was to conduct measurements of physical, chemical and biological variables on a section across the Mediterranean Sea with the goal of producing a synoptic picture of the distribution of relevant physical and biogeochemical properties, in order to compare those to historic data sets. During the cruise, a comprehensive data set of relevant variables following the guide lines for repeat hydrography outlined by the GO-SHIP group (http://www.go-ship.org/) was collected. The measurements include salinity and temperature (CTD), an over-determined carbonate system, inorganic nutrients, oxygen, transient tracers (CFC-12, SF6), helium isotopes and tritium, and carbon isotopes. The cruise sampled all major basins of the Mediterranean Sea following roughly an east-to-west section from the coast of Lebanon through to the Strait of Gibraltar, and to the coast of Portugal. Also a south-to-north section from the Ionian Sea to the Adriatic Sea was carried out. Additionally, sampling in the Aegean, Adriatic and Tyrrhenian Seas were carried out. The sections roughly followed lines and positions that have been sampled previously during other programs, thus providing the opportunity for comparative investigations of the temporal development of various parameters.

Planktic foraminiferal Mg/Ca and oxygen isotopes, plus estimated d18Oseawater and ice volume corrected d18Oseawater for the past 35000 years of sediment core GeoB10069-3

Evidence from geologic archives suggests that there were large changes in the tropical hydrologic cycle associated with the two prominent northern hemisphere deglacial cooling events, Heinrich Stadial 1 (HS1; ~19 to 15 kyr BP; kyr BP = 1000 yr before present) and the Younger Dryas (~12.9 to 11.7 kyr BP). These hydrologic shifts have been alternatively attributed to high and low latitude origin. Here, we present a new record of hydrologic variability based on planktic foraminifera-derived d18O of seawater (d18Osw) estimates from a sediment core from the tropical Eastern Indian Ocean, and using 12 additional d18Osw records, construct a single record of the dominant mode of tropical Eastern Equatorial Pacific and Indo-Pacific Warm Pool (IPWP) hydrologic variability. We show that deglacial hydrologic shifts parallel variations in the reconstructed interhemispheric temperature gradient, suggesting a strong response to variations in the Atlantic Meridional Overturning Circulation and the attendant heat redistribution. A transient model simulation of the last deglaciation suggests that hydrologic changes, including a southward shift in the Intertropical Convergence Zone (ITCZ) which likely occurred during these northern hemisphere cold events, coupled with oceanic advection and mixing, resulted in increased salinity in the Indonesian region of the IPWP and the eastern tropical Pacific, which is recorded by the d18Osw proxy. Based on our observations and modeling results we suggest the interhemispheric temperature gradient directly controls the tropical hydrologic cycle on these time scales, which in turn mediates poleward atmospheric heat transport.

Geochemistry at DSDP Leg 82 Holes

The nine holes (556-564) drilled during DSDP Leg 82 in a region west and southwest of the Azores Platform (Fig. 1) exhibit a wide variety of chemical compositions that indicate a complex petrogenetic history involving crystal fractionation, magma mixing, complex melting, and mantle heterogeneity. The major element chemistry of each hole except Hole 557 is typical of mid-ocean ridge basalts (MORBs), whereas the trace element and rare earth element (REE) abundances and ratios are more variable, and show that both depleted Type I and enriched Type II basalts have been erupted in the region. Hole 556 (30-34 Ma), located near a flow line through the Azores Triple Junction, contains typically depleted basalts, whereas Hole 557 (18 Ma), located near the same flow line but closer to the Azores Platform, is a highly enriched FeTi basalt, indicating that the Azores hot-spot anomaly has existed in its present configuration for at least 18 Ma, but less than 30-34 Ma. Hole 558 (34-37 Ma), located near a flow line through the FAMOUS and Leg 37 sites, includes both Type I and II basalts. Although the differences in Zr/Nb and light REE/heavy REE ratios imply different mantle sources, the (La/Ce)ch (>1) and Nd isotopic ratios are almost the same, suggesting that the complex melting and pervasive, small-scale mantle heterogeneity may account for the variations in trace element and REE ratios observed in Hole 558 (and FAMOUS sites). Farther south, Hole 559 (34-37 Ma), contains enriched Type II basalts, whereas Hole 561 (14-17 Ma), located further east near the same flow line, contains Type I and II basalts. In this case, the (La/Ce)ch and Nd isotopic ratios are different, indicating two distinct mantle sources. Again, the existence along the same flow line of two holes exhibiting such different chemistry suggests that mantle heterogeneity may exist on a more pervasive and transient smaller scale. (Hole 560 was not sampled for this study because the single basalt clast recovered was used for shipboard analysis.) All of the remaining three holes (562, 563, 564), located along a flow line about 100 km south of the Hayes Fracture Zone (33°N), contain only depleted Type I basalts. The contrast in chemical compositions suggests that the Hayes Fracture Zone may act as a "domain" boundary between an area of fairly homogeneous, depleted Type I basalts to the south (Holes 562-564) and a region of complex, highly variable basalts to the north near the Azores hot-spot anomaly (Holes 556-561).