Formation of seep carbonates along the makran convergent margin, northern Arabian sea and a molecular and isotopic approach to constrain the carbon isotopic composition of parent methane

Authigenic carbonate deposits have been sampled with the remotely operated vehicle ‘MARUM-QUEST 4000 m’ from five methane seeps between 731 and 1823 m water depth along the convergent Makran continental margin, offshore Pakistan (northern Arabian Sea). Two seeps on the upper slope are located within the oxygen minimum zone (OMZ; ca. 100 to 1100 m water depth), the other sites are situated in oxygenated water below the OMZ (below 1100 m water depth). The carbonate deposits vary with regard to their spatial extent, sedimentary fabrics, and associated seep fauna: Within the OMZ, carbonates are spatially restricted and associated with microbial mats, whereas in the oxygenated zone below the OMZ extensive carbonate crusts are exposed on the seafloor with abundant metazoans (bathymodiolin mussels, tube worms, galatheid crabs). Aragonite and Mg-calcite are the dominant carbonate minerals, forming common early diagenetic microcrystalline cement and clotted to radial-fibrous cement. The δ18Ocarbonate values range from 1.3 to 4.2‰ V-PDB, indicating carbonate precipitation at ambient bottom-water temperature in shallow sediment depth. Extremely low δ13Ccarbonate values (as low − 54.6‰ V-PDB) point to anaerobic oxidation of methane (AOM) as trigger for carbonate precipitation, with biogenic methane as dominant carbon source. Prevalence of biogenic methane in the seepage gas is corroborated by δ13Cmethane values ranging from − 70.3 to − 66.7‰ V-PDB, and also by back-calculations considering δ13Cmethane values of carbonate and incorporated lipid biomarkers. These calculations (Δδ13Cmethane–carbonate, Δδ13CANME–methane, Δδ13CMOX–methane) prove to be useful to assess the carbon stable isotope composition of seeping methane if this has not been determined in the first place; such an approach represents a useful tool to reconstruct fluid composition of ancient seeps. AOM is also revealed by lipid biomarkers of anaerobic methane oxidizing archaea such as crocetane, pentamethylicosane (PMI), and sn2-hydroxyarchaeol strongly depleted in 13C (δ13C values as low as − 127‰ V-PDB). Biomarkers of sulphate-reducing bacteria are also abundant, showing slightly less negative δ13C values, but still significantly 13C-depleted (average values as low as − 101‰). Other bacterial biomarkers, such as bacteriohopanepolyols (BHPs), hopanols, and hopanoic acids are detected in most carbonates, but are particularly common in seep carbonates from the non-OMZ sites. The BHP patterns of these carbonates and their low δ13C values resemble patterns of aerobic methanotrophic bacteria. In the shallower OMZ sites, BHPs revealed much lower contents and varying compositions, most likely reflecting other sources than aerobic methanotrophic bacteria. 230Th/U carbonate ages indicate that AOM-induced carbonate precipitation at the deeper non-OMZ seeps occurred mainly during the late Pleistocene-Holocene transition, i.e. between 19 and 15 ka before present, when the global sea level was lower than today.

Mineralogical assemblage in sediment core GeoB12309-5 and in area samples from the Makran subduction zone

In order to study late Holocene changes in sediment supply into the northern Arabian Sea, a 5.3 m long gravity core was investigated by high-resolution geochemical and mineralogical techniques. The sediment core was recovered at a water depth of 956 m from the continental slope off Pakistan and covers a time span of 5 kyr. During the late Holocene source areas delivering material to the sampling site did, however, not change and were active throughout the year.

The hazard potential of the western segment of the Makran subduction zone, northern Arabian Sea

Evaluating the hazard potential of the Makran subduction zone requires understanding the previous records of the large earthquakes and tsunamis. We address this problem by searching for earthquake and tectonic proxies along the Makran Coast and linking those observations with the available constraints on historical seismicity and the tell-tale characteristics of sea floor morphology. The earthquake of Mw 8.1 of 1945 and the consequent tsunami that originated on the eastern part of the Makran are the only historically known hazardous events in this region. The seismic status of the western part of the subduction zone outside the rupture area of the 1945 earthquake remains an enigma. The near-shore shallow stratigraphy of the central part of Makran near Chabahar shows evidence of seismically induced liquefaction that we attribute to the distant effects of the 1945 earthquake. The coastal sites further westward around Jask are remarkable for the absence of liquefaction features, at least at the shallow level. Although a negative evidence, this possibly implies that the western part of Makran Coast region may not have been impacted by near-field large earthquakes in the recent past-a fact also supported by the analysis of historical data. On the other hand, the elevated marine terraces on the western Makran and their uplift rates are indicative of comparable degree of long-term tectonic activity, at least around Chabahar. The offshore data suggest occurrences of recently active submarine slumps on the eastern part of the Makran, reflective of shaking events, owing to the great 1945 earthquake. The ocean floor morphologic features on the western segment, on the contrary, are much subdued and the prograding delta lobes on the shelf edge also remain intact. The coast on the western Makran, in general, shows indications of progradation and uplift. The various lines of evidence thus suggest that although the western segment is potentially seismogenic, large earthquakes have not occurred there in the recent past, at least during the last 600 years. The recurrence period of earthquakes may range up to 1,000 years or more, an assessment based on the age of the youngest dated coastal ridge. The long elapsed time points to the fact that the western segment may have accumulated sufficient slip to produce a major earthquake.

Methane seep carbonates yield clumped isotope signatures out of equilibrium with formation temperatures

Methane cold seep systems typically exhibit extensive buildups of authigenic carbonate minerals, resulting from local increases in alkalinity driven by methane oxidation. Here, we demonstrate that modern seep authigenic carbonates exhibit anomalously low clumped isotope values (Delta(47)), as much as similar to 0.2 parts per thousand lower than expected values. In modern seeps, this range of disequilibrium translates into apparent temperatures that are always warmer than ambient temperatures, by up to 50 degrees C. We examine various mechanisms that may induce disequilibrium behaviour in modern seep carbonates, and suggest that the observed values result from several factors including kinetic isotopic effects during methane oxidation, mixing of inorganic carbon pools, pH effects and rapid precipitation. Ancient seep carbonates studied here also exhibit potential disequilibrium signals. Ultimately, these findings indicate the predominance of disequilibrium clumped isotope behaviour in modern cold seep carbonates that must be considered when characterizing environmental conditions in both modern and ancient cold seep settings.

Isotopes content, mineralogical composition, Uranium and Thorium contents and determined activity ratios from METEOR cruise M74/3 in 2007

Authigenic carbonate deposits have been sampled with the remotely operated vehicle 'MARUM-QUEST 4000 m' from five methane seeps between 731 and 1823 m water depth along the convergent Makran continental margin, offshore Pakistan (northern Arabian Sea). Two seeps on the upper slope are located within the oxygen minimum zone (OMZ; ca. 100 to 1100 m water depth), the other sites are situated in oxygenated water below the OMZ (below 1100 m water depth). The carbonate deposits vary with regard to their spatial extent, sedimentary fabrics, and associated seep fauna: Within the OMZ, carbonates are spatially restricted and associated with microbial mats, whereas in the oxygenated zone below the OMZ extensive carbonate crusts are exposed on the seafloor with abundant metazoans (bathymodiolin mussels, tube worms, galatheid crabs). Aragonite and Mg-calcite are the dominant carbonate minerals, forming common early diagenetic microcrystalline cement and clotted to radial-fibrous cement. The delta18O carbonate values range from 1.3 to 4.2 per mil V-PDB, indicating carbonate precipitation at ambient bottom-water temperature in shallow sediment depth. Extremely low delta13Ccarbonate values (as low - 54.6per mil V-PDB) point to anaerobic oxidation of methane (AOM) as trigger for carbonate precipitation, with biogenic methane as dominant carbon source. Prevalence of biogenic methane in the seepage gas is corroborated by delta13C methane values ranging from - 70.3 to - 66.7per mil V-PDB, and also by back-calculations considering delta 13C methane values of carbonate and incorporated lipid biomarkers.

Fluid-related inclusions in Alpine high-pressure peridotite reveal trace element recycling during subduction-zone dehydration of serpentinized mantle (Cima di Gagnone, Swiss Alps).

Serpentinites release at sub-arc depths volatiles and several fluid-mobile trace elements found in arc magmas. Constraining element uptake in these rocks and defining the trace element composition of fluids released upon serpentinite dehydration can improve our understanding of mass transfer across subduction zones and to volcanic arcs. The eclogite-facies garnet metaperidotite and chlorite harzburgite bodies embedded in paragneiss of the subduction melange from Cima di Gagnone derive from serpentinized peridotite protoliths and are unique examples of ultramafic rocks that experienced subduction metasomatism and devolatilization. In these rocks, metamorphic olivine and garnet trap polyphase inclusions representing the fluid released during high-pressure breakdown of antigorite and chlorite. Combining major element mapping and laser-ablation ICP-MS bulk inclusion analysis, we characterize the mineral content of polyphase inclusions and quantify the fluid composition. Silicates, Cl-bearing phases, sulphides, carbonates, and oxides document post-entrapment mineral growth in the inclusions starting immediately after fluid entrapment. Compositional data reveal the presence of two different fluid types. The first (type A) records a fluid prominently enriched in fluid-mobile elements, with Cl, Cs, Pb, As, Sb concentrations up to 10(3) PM (primitive mantle), similar to 10(2) PM Tit Ba, while Rb, B, Sr, Li, U concentrations are of the order of 10(1) PM, and alkalis are similar to 2 PM. The second fluid (type B) has considerably lower fluid-mobile element enrichments, but its enrichment patterns are comparable to type A fluid. Our data reveal multistage fluid uptake in these peridotite bodies, including selective element enrichment during seafloor alteration, followed by fluid-rock interaction along with subduction metamorphism in the plate interface melange. Here, infiltration of sediment-equilibrated fluid produced significant enrichment of the serpentinites in As, Sb, B, Pb, an enriched trace element pattern that was then transferred to the fluid released at greater depth upon serpentine dehydration (type A fluid). The type B fluid hosted by garnet may record the composition of the chlorite breakdown fluid released at even greater depth. The Gagnone study-case demonstrates that serpentinized peridotites acquire water and fluid-mobile elements during ocean floor hydration and through exchange with sediment-equilibrated fluids in the early subduction stages. Subsequent antigorite devolatilization at subarc depths delivers aqueous fluids to the mantle wedge that can be prominently enriched in sediment-derived components, potentially triggering arc magmatism without the need of concomitant dehydration/melting of metasediments or altered oceanic crust.

Biogeochemical measurements of cold seep sediments inhabit by vesicomyid clams in the Japan Deep Sea Trench

Vesicomyidae clams harbor sulfide-oxidizing endosymbionts and are typical members of cold seep communities associated with tectonic faults where active venting of fluids and gases takes place. We investigated the central biogeochemical processes that supported a vesicomyid clam colony as part of a locally restricted seep community in the Japan Trench at 5346 m water depth, one of the deepest seep settings studied to date. An integrated approach of biogeochemical and molecular ecological techniques was used combining in situ and ex situ measurements. During the cruise YK06-05 in 2006 with the RV Yokosuka to the Japan Trench, we investigated a clam colony inhabited by Abyssogena phaseoliformis (former known as Calyptogena phaseoliformis) and Isorropodon fossajaponicum (former known as Calyptogena fossajaponica). The targeted sampling and precise positioning of the in situ instruments were achieved with the manned research submersible Shinkai 6500 (JAMSTEC, Nankoku, Kochi, Japan). Sampling was first performed close to the rim of the JTC colony and then at the center. Immediately after sample recovery onboard, the sediment core was sub-sampled for ex situ rate measurements or preserved for later analyses. In sediment of the clam colony, low sulfate reduction (SR) rates (max. 128 nmol ml**-1 d**-1) were coupled to the anaerobic oxidation of methane (AOM). They were observed over a depth range of 15 cm, caused by active transport of sulfate due to bioturbation of the vesicomyid clams. A distinct separation between the seep and the surrounding seafloor was shown by steep horizontal geochemical gradients and pronounced microbial community shifts. The sediment below the clam colony was dominated by anaerobic methanotrophic archaea (ANME-2c) and sulfate-reducing Desulfobulbaceae (SEEP-SRB-3, SEEP-SRB-4). Aerobic methanotrophic bacteria were not detected in the sediment and the oxidation of sulfide seemed to be carried out chemolithoautotrophically by Sulfurovum species. Thus, major redox processes were mediated by distinct subgroups of seep-related microorganisms that might have been selected by this specific abyssal seep environment. Fluid flow and microbial activity was low but sufficient to support the clam community over decades and to build up high biomasses. Hence, the clams and their microbial communities adapted successfully to a low-energy regime and may represent widespread chemosynthetic communities in the Japan Trench.

Mineralogy of cold seep sediment inhabited by vesicomyid clams in the Japan Deep Sea Trench

Sediment sampling was performed at the center of the clam colony. Immediately after sample recovery onboard, the sediment core was sub-sampled and preserved for later analyses. Pyrite and carbonate content of the sediment was measured by X-ray refraction analysis as previously described (Ertefai et al., 2010).

On the recurrence of great subduction zone earthquakes

The last decade has witnessed two unusually large tsunamigenic earthquakes. The devastation from the 2004 Sumatra Andaman and the 2011 Tohoku-Oki earthquakes (both of moment magnitude >= 9.0) and their ensuing tsunamis comes as a harsh reminder on the need to assess and mitigate coastal hazards due to earthquakes and tsunamis worldwide. Along any given subduction zone, megathrust tsunamigenic earthquakes occur over intervals considerably longer than their documented histories and thus, 2004-type events may appear totally `out of the blue'. In order to understand and assess the risk from tsunamis, we need to know their long-term frequency and magnitude, going beyond documented history, to recent geological records. The ability to do this depends on our knowledge of the processes that govern subduction zones, their responses to interseismic and coseismic deformation, and on our expertise to identify and relate tsunami deposits to earthquake sources. In this article, we review the current state of understanding on the recurrence of great thrust earthquakes along global subduction zones.

Elastic thickness structure of the Andaman subduction zone: Implications for convergence of the Ninetyeast Ridge

We use the Bouguer coherence (Morlet isostatic response function) technique to compute the spatial variation of effective elastic thickness (T-e) of the Andaman subduction zone. The recovered T-e map resolves regional-scale features that correlate well with known surface structures of the subducting Indian plate and the overriding Burma plate. The major structure on the India plate, the Ninetyeast Ridge (NER), exhibits a weak mechanical strength, which is consistent with the expected signature of an oceanic ridge of hotspot origin. However, a markedly low strength (0< T-e <3 km) in that region, where the NER is close to the Andaman trench (north of 10 N), receives our main attention in this study. The subduction geometry derived from the Bouguer gravity forward modeling suggests that the NER has indented beneath the Andaman arc. We infer that the bending stresses of the viscous plate, which were reinforced within the subducting oceanic plate as a result of the partial subduction of the NER buoyant load, have reduced the lithospheric strength. The correlation, T-e < T-s (seismogenic thickness) reveals that the upper crust is actively deforming beneath the frontal arc Andaman region. The occurrence of normal-fault earthquakes in the frontal arc, low Te zone, is indicative of structural heterogeneities within the subducting plate. The fact that the NER along with its buoyant root is subducting under the Andaman region is inhibiting the subduction processes, as suggested by the changes in trench line, interrupted back-arc volcanism, variation in seismicity mechanism, slow subduction, etc. The low T-e and thinned crustal structure of the Andaman back-arc basin are attributed to a thermomechanically weakened lithosphere. The present study reveals that the ongoing back-arc spreading and strike-slip motion along the West Andaman Fault coupled with the ridge subduction exerts an important control on the frequency and magnitude of seismicity in the Andaman region. (C) 2013 Elsevier Ltd. All rights reserved.

Using space geodesy to constrain variations in seismogenic behavior on subduction megathrusts

The concept of seismogenic asperities and aseismic barriers has become a useful paradigm within which to understand the seismogenic behavior of major faults. Since asperities and barriers can be thought of as defining the potential rupture area of large megathrust earthquakes, it is thus important to identify their respective spatial extents, constrain their temporal longevity, and to develop a physical understanding for their behavior. Space geodesy is making critical contributions to the identification of slip asperities and barriers but progress in many geographical regions depends on improving the accuracy and precision of the basic measurements. This thesis begins with technical developments aimed at improving satellite radar interferometric measurements of ground deformation whereby we introduce an empirical correction algorithm for unwanted effects due to interferometric path delays that are due to spatially and temporally variable radar wave propagation speeds in the atmosphere. In chapter 2, I combine geodetic datasets with complementary spatio-temporal resolutions to improve our understanding of the spatial distribution of crustal deformation sources and their associated temporal evolution – here we use observations from Long Valley Caldera (California) as our test bed. In the third chapter I apply the tools developed in the first two chapters to analyze postseismic deformation associated with the 2010 Mw=8.8 Maule (Chile) earthquake. The result delimits patches where afterslip occurs, explores their relationship to coseismic rupture, quantifies frictional properties associated with inferred patches of afterslip, and discusses the relationship of asperities and barriers to long-term topography. The final chapter investigates interseismic deformation of the eastern Makran subduction zone by using satellite radar interferometry only, and demonstrates that with state-of-art techniques it is possible to quantify tectonic signals with small amplitude and long wavelength. Portions of the eastern Makran for which we estimate low fault coupling correspond to areas where bathymetric features on the downgoing plate are presently subducting, whereas the region of the 1945 M=8.1 earthquake appears to be more highly coupled.

I. A seismotectonic study of the Middle America Subduction Zone. II. Lithosphere and upper mantle structure of the Canadian Shield and Eastern North America

This thesis consists of two separate parts. Part I (Chapter 1) is concerned with seismotectonics of the Middle America subduction zone. In this chapter, stress distribution and Benioff zone geometry are investigated along almost 2000 km of this subduction zone, from the Rivera Fracture Zone in the north to Guatemala in the south. Particular emphasis is placed on the effects on stress distribution of two aseismic ridges, the Tehuantepec Ridge and the Orozco Fracture Zone, which subduct at seismic gaps. Stress distribution is determined by studying seismicity distribution, and by analysis of 190 focal mechanisms, both new and previously published, which are collected here. In addition, two recent large earthquakes that have occurred near the Tehuantepec Ridge and the Orozco Fracture Zone are discussed in more detail. A consistent stress release pattern is found along most of the Middle America subduction zone: thrust events at shallow depths, followed down-dip by an area of low seismic activity, followed by a zone of normal events at over 175 km from the trench and 60 km depth. The zone of low activity is interpreted as showing decoupling of the plates, and the zone of normal activity as showing the breakup of the descending plate. The portion of subducted lithosphere containing the Orozco Fracture Zone does not differ significantly, in Benioff zone geometry or in stress distribution, from adjoining segments. The Playa Azul earthquake of October 25, 1981, M_s=7.3, occurred in this area. Body and surface wave analysis of this event shows a simple source with a shallow thrust mechanism and gives M_o=1.3x10²⁷ dyne-cm. A stress drop of about 45 bars is calculated; this is slightly higher than that of other thrust events in this subduction zone. In the Tehuantepec Ridge area, only minor differences in stress distribution are seen relative to adjoining segments. For both ridges, the only major difference from adjoining areas is the infrequency or lack of occurrence of large interplate thrust events.

Part II involves upper mantle P wave structure studies, for the Canadian shield and eastern North America. In Chapter 2, the P wave structure of the Canadian shield is determined through forward waveform modeling of the phases P_nl, P, and PP. Effects of lateral heterogeneity are kept to a minimum by using earthquakes just outside the shield as sources, with propagation paths largely within the shield. Previous mantle structure studies have used recordings of P waves in the upper mantle triplication range of 15-30°; however, the lack of large earthquakes in the shield region makes compilation of a complete P wave dataset difficult. By using the phase PP, which undergoes triplications at 30-60°, much more information becomes available. The WKBJ technique is used to calculate synthetic seismograms for PP, and these records are modeled almost as well as the P. A new velocity model, designated S25, is proposed for the Canadian shield. This model contains a thick, high-Q, high-velocity lid to 165 km and a deep low-velocity zone. These features combine to produce seismograms that are markedly different from those generated by other shield structure models. The upper mantle discontinuities in S25 are placed at 405 and 660 km, with a simple linear gradient in velocity between them. Details of the shape of the discontinuities are not well constrained. Below 405 km, this model is not very different from many proposed P wave models for both shield and tectonic regions.

Chapter 3 looks in more detail at recordings of P_nl in eastern North America. First, seismograms from four eastern North American earthquakes are analyzed, and seismic moments for the events are calculated. These earthquakes are important in that they are among the largest to have occurred in eastern North America in the last thirty years, yet in some cases were not large enough to produce many good long-period teleseismic records. A simple layer-over-a-halfspace model is used for the initial modeling, and is found to provide an excellent fit for many features of the observed waveforms. The effects on P_nl of varying lid structure are then investigated. A thick lid with a positive gradient in velocity, such as that proposed for the Canadian shield in Chapter 2, will have a pronounced effect on the waveforms, beginning at distances of 800 or 900 km. P_nl records from the same eastern North American events are recalculated for several lid structure models, to survey what kinds of variations might be seen. For several records it is possible to see likely effects of lid structure in the data. However, the dataset is too sparse to make any general observations about variations in lid structure. This type of modeling is expected to be important in the future, as the analysis is extended to more recent eastern North American events, and as broadband instruments make more high-quality regional recordings available.

Fatty-acids and their C-13 signatures in seep carbonates from hydrocarbon seeps on the upper (GC 185) and lower (AC 645) continental slope of the Gulf of Mexico

Here we reported the fatty-acids and their δ 13C values in seep carbonates collected from Green Canyon lease block 185 (GC 185; Sample GC-F) at upper continental slope (water depth: ∼540 m), and Alaminos Canyon lease block 645 (GC 645; Sample AC-E) at lower continental slope (water depth: ∼2200 m) of the Gulf of Mexico. More than thirty kinds of fatty acids were detected in both samples. These fatty acids are maximized at C16. There is a clear even-over-odd carbon number predominance in carbon number range. The fatty acids are mainly composed of n-fatty acids, iso-/anteiso-fatty acids and terminally branched odd-numbered fatty acids (iso/anteiso). The low δ 13C values (−39.99‰ to.32.36‰) of n-C12:0, n-C13:0, i-C14:0and n-C14:0 suggest that they may relate to the chemosynthetic communities at seep sites. The unsaturated fatty acids n-C18:2 and C18:1Δ9 have the same δ 13C values, they may originate from theBeggiatoa/Thioploca. Unlike other fatty acids, the terminally branched fatty acids (iso/anteiso) show lowerδ 13C values (as low as −63.95‰) suggesting a possible relationship to sulfate reducing bacteria, which is common during anaerobic oxidation of methane at seep sites.