Laboratory shearing experiments at low sliding velocities of IODP Hole 343-C0019E

The 2011 Tohoku-Oki earthquake demonstrated that the shallowest reaches of plate boundary subduction megathrusts can host substantial coseismic slip that generates large and destructive tsunamis, contrary to the common assumption that the frictional properties of unconsolidated clay-rich sediments at depths less than View the MathML source should inhibit rupture. We report on laboratory shearing experiments at low sliding velocities (View the MathML source) using borehole samples recovered during IODP Expedition 343 (JFAST), spanning the plate-boundary décollement within the region of large coseismic slip during the Tohoku earthquake. We show that at sub-seismic slip rates the fault is weak (sliding friction µs=0.2-0.26), in contrast to the much stronger wall rocks (µs>~0.5). The fault is weak due to elevated smectite clay content and is frictionally similar to a pelagic clay layer of similar composition. The higher cohesion of intact wall rock samples coupled with their higher amorphous silica content suggests that the wall rock is stronger due to diagenetic cementation and low clay content. Our measurements also show that the strongly developed in-situ fabric in the fault zone does not contribute to its frictional weakness, but does lead to a near-cohesionless fault zone, which may facilitate rupture propagation by reducing shear strength and surface energy at the tip of the rupture front. We suggest that the shallow rupture and large coseismic slip during the 2011 Tohoku earthquake was facilitated by a weak and cohesionless fault combined with strong wall rocks that drive localized deformation within a narrow zone.

Shear experiments from IODP Expedition 322 and 333 samples

The mechanical behavior of the plate boundary fault zone is of paramount importance in subduction zones, because it controls megathrust earthquake nucleation and propagation as well as the structural style of the forearc. In the Nankai area along the NanTroSEIZE (Kumano) drilling transect offshore SW Japan, a heterogeneous sedimentary sequence overlying the oceanic crust enters the subduction zone. In order to predict how variations in lithology, and thus mechanical properties, affect the formation and evolution of the plate boundary fault, we conducted laboratory tests measuring the shear strengths of sediments approaching the trench covering each major lithological sedimentary unit. We observe that shear strength increases nonlinearly with depth, such that the (apparent) coefficient of friction decreases. In combination with a critical taper analysis, the results imply that the plate boundary position is located on the main frontal thrust. Further landward, the plate boundary is expected to step down into progressively lower stratigraphic units, assisted by moderately elevated pore pressures. As seismogenic depths are approached, the décollement may further step down to lower volcaniclastic or pelagic strata but this requires specific overpressure conditions. High-taper angle and elevated strengths in the toe region may be local features restricted to the Kumano transect.

(Table 1) Experimental results of Kd at various conditions from ODP Hole 131-808C

A comprehensive experimental study, utilizing a rocking autoclave hydrothermal apparatus with isotope tracers, was applied to evaluate the temperature of squeezing artifacts on B contents and isotopic compositions in pore waters. The partition coefficient (KD) was determined at temperatures from 25 ° to 350 °C, at 800 bars, and this information was applied to reconstruct pore water B and d11B in ODP drill sites, where pH, T, and porosity are known. The partition coefficient of B is a function of temperature, pH, and sediment mineralogy. The solution pH exerts a dominant control at low temperatures; however, KD decreases to a value of essentially zero (compared to that of KD = ~3.5 at 25 °C) at high temperatures indicating no adsorption. Two empirical equations were derived to represent most of the available experimental results. For pelagic clay rich sediments, a KD = -3.84-0.020T + 0.88pH (R = 0.84; 1sigma = 0.25) is established. For sediments that have experienced progressive metamorphism, a KD = -1.38-0.008T + 0.59pH (R = 0.81; 1sigma = 0.37) can be applied. Similarly the effect on pore water d11B can be corrected if the fractionation factors at different temperatures are assumed. The corrected B and d11B in ODP Sites 671, 672, and 808 indicate significant mobilization of bulk B in sediment (exchangeable + lattice bound) at depth, especially near the décollement zone or other potential flow conduits. Tectonically expelled fluids from mud diapirs of Barbados Ridge Complex, hot springs of Rumsey Hills, California, and mud pot waters of the Salton Sea geothermal field, are enriched in B (up to 20 mM) with lower d11B, supporting the argument of B mobilization as a result of fluid expulsion in accretionary prisms.

Permeability measurements on sheared and unsheared samples from different IODP Holes of Leg 316 and from IODP Hole 315-C0001E and ODP Hole 190-1174B

At subduction zones, the permeability of major fault zones influences pore pressure generation, controls fluid flow pathways and rates, and affects fault slip behavior and mechanical strength by mediating effective normal stress. Therefore, there is a need for detailed and systematic permeability measurements of natural materials from fault systems, particularly measurements that allow direct comparison between the permeability of sheared and unsheared samples from the same host rock or sediment. We conducted laboratory experiments to compare the permeability of sheared and uniaxially consolidated (unsheared) marine sediments sampled during IODP Expedition 316 and ODP Leg 190 to the Nankai Trough offshore Japan. These samples were retrieved from: (1) The décollement zone and incoming trench fill offshore Shikoku Island (the Muroto transect); (2) Slope sediments sampled offshore SW Honshu (the Kumano transect) ~ 25 km landward of the trench, including material overriden by a major out-of-sequence thrust fault, termed the "megasplay"; and (3) A region of diffuse thrust faulting near the toe of the accretionary prism along the Kumano transect. Our results show that shearing reduces fault-normal permeability by up to 1 order of magnitude, and this reduction is largest for shallow (< 500 mbsf) samples. Shearing-induced permeability reduction is smaller in samples from greater depth, where pre-existing fabric from compaction and lithification may be better developed. Our results indicate that localized shearing in fault zones should result in heterogeneous permeability in the uppermost few kilometers in accretionary prisms, which favors both the trapping of fluids beneath and within major faults, and the channeling of flow parallel to fault structure. These low permeabilities promote the development of elevated pore fluid pressures during accretion and underthrusting, and will also facilitate dynamic hydrologic processes within shear zones including dilatancy hardening and thermal pressurization.

Boron contents and boron, carbon, and oxygen isotope data for ODP Sites 160-970 and 160-971

Products of two mud volcanoes from the distal part of the Mediterranean Ridge accretionary complex have been investigated regarding their B, C, and O stable isotope signatures. The mud breccias have been divided into mud matrix, lithified clasts, biogenic deposits, and authigenic cements and crusts related to fluid flow and cementation. Isotope geochemistry is used to evaluate the depth of mobilization of each phase in the subduction zone. B contents and isotope ratios of the mud and mud clasts show a general trend of B enrichment and decreasing d11B values with increasing consolidation (i.e., depth). However, the majority of the clast and matrix samples relate to moderate depths of mobilization within the wedge (1-2 km below seafloor). The carbonate cements of most of these clasts as well as the authigenic crusts, however, provide evidence for a deep fluid influence, probably associated with the décollement at 5-6 km depth. This interpretation is supported by d13C ratios of the crust, which indicate precipitation of C from thermogenic methane, and by the d11B ratios of pore-water samples of mud-breccia drill cores. Clams (Vesicomya sp.) living adjacent to fluid vents have d11B and d18O values corresponding to brines known in the area, which acted as the parent solution for shell precipitation. Such brines are most likely Miocene pore waters trapped at deep levels within the backstop to the accretionary prism, probably prior to desiccation of the Mediterranean in the Messinian (6-5 Ma). Combining all results, deep fluid circulation and expulsion are identified as the main processes triggering mud liquefaction and extrusion, whereas brines contribute only locally. Given the high B contents, mud extrusion has to be considered a major backflux mechanism of B into the hydrosphere.

Common depth points and associated geographic coordinates

Fluids in subduction zones can influence seismogenic behaviour and prism morphology. The Eastern Makran subduction zone, offshore Pakistan, has a very thick incoming sediment section of up to 7.5 km, providing a large potential fluid source to the accretionary prism. A hydrate-related bottom simulating reflector (BSR), zones of high amplitude reflectivity, seafloor seep sites and reflective thrust faults are present across the accretionary prism, indicating the presence of fluids and suggesting active fluid migration. High amplitude free gas zones and seep sites are primarily associated with anticlinal hinge traps, and fluids here appear to be sourced from shallow biogenic sources and migrate to the seafloor along minor normal faults. There are no observed seep sites associated with the surface expression of the wedge thrust faults, potentially due to burial of the surface trace by failure of the steep thrust ridge slopes. Thrust fault reflectivity is restricted to the upper 3 km of sediment and the deeper décollement is non-reflective. We interpret that fluids and overpressure are not common in the deeper stratigraphic section. Thermal modelling of sediments at the deformation front suggests that the deeper sediment section is relatively dewatered and not currently contributing to fluid expulsion in the Makran accretionary prism.

(Table 1) Chemical compositions of beryllium in sediments and pore waters of ODP Hole 131-808

Studies of Be distributions in subduction zone sediments will help to understand questions regarding the enrichments of cosmogenic Be-10 in arc volcanic rocks. Analyses of Be-10 and Be-9 in sediments of Ocean Drilling Program Site 808, Nankai Trough and Be-9 in porewaters of Site 808 and Sites 671 and 672, Barbados ridge complex, show significant decreases in solid phase Be-10 and large increases of porewater Be-9 at the location of the décollement zone and below or at potential flow conduits. These data imply the potential mobilization of Be during pore fluid expulsion upon sediment burial. Experiments involving reaction between a décollement sediment and a synthetic NaCl-CaCl2 solution at elevated pressure and temperatures were conducted in an attempt to mimic early subduction zone processes. The results demonstrate that Be is mobilized under elevated pressure and temperature with a strong pH dependence. The Be mobilization provides an explanation of Be-10 enrichment in arc volcanic rocks and supports the argument of the importance of the fluid processes in subduction zones at convergent margins.

Whole rock and elements concentrations from IODP Hole 343-C0019E

The microstructures, mineralogy and chemistry of four representative samples collected from cores extracted from the Japan Trench during Integrated Ocean Drilling Project Expedition 343, the Japan Trench Fast Drilling Project (JFAST) have been studied using optical microscopy, TEM, SEM, XRF, XRD and microprobe analyses. The samples provide a transect from relatively undeformed marine sediments in the hanging wall, to the undeformed footwall material, crossing the thrust interface between the Pacific and North American plate, where the fault slipped during the March 2011 Tohoku-Oki earthquake. Our preliminary results suggest that the low strength of JFAST fault gouge material is caused by the high amount of clay minerals (~ 60% smectite, ~ 14 illite). Clay minerals in the décollement (gouge) sample are partly replaced by newly formed manganese oxide, which precipitated from hydrothermal fluids. Dauphine twins were found in quartz grains of the décollement sample suggesting local high stress possible during seismic loading. Other microstructures cannot be assigned unambiguously to co-seismic or a-seismic faulting processes. The observed scaly clay fabric is consistent with observations in many other plate-boundary fault zones. Significant grain size reduction was found in the fault (decollement) zone sample. But a change in lithology of the fault material cannot be ruled out. Microstructures typical for a-seismic deformation like dissolution-precipitation features (e.g. dissolved grain boundaries, mineral alteration) occur in all JFAST core samples, but more frequently in the décollement sample.

Lithium isotopic composition of pore fluids and sediments in the Costa Rica subduction zone

Pore fluid and sediment Li concentrations and isotopic ratios provide important insights on the hydrology, sediment contribution to the arc volcanoes and fluid-sediment reactions at the dominantly non-accretionary Costa Rica subduction zone. Ocean Drilling Program Site 1039 in the trench axis provides a reference section of 400 m of the incoming sediments, and Site 1040, situated arcward from the trench, consists of a deformed sedimentary wedge and apron sediments, the décollement, and the partially dewatered underthrust sediment section. At the reference site, pore fluids show important isotopic variations (delta6Li=-21.7 to -37.8 per mil), reflecting the interplay of in situ alteration of volcanic material and ion exchange with clay minerals. In the basal section, a reversal of Li concentration and delta6Li toward seawater values is observed, providing supporting evidence for a lateral seawater flow system in the upper oceanic basement underlying this sediment section. At Site 1040, pore fluid of the lower deformed wedge sediments and within the décollement is enriched in Li and the isotopic compositions are relatively light, suggesting infiltration of a deep-seated fluid. The delta6Li value of -22 per mil of this Li-enriched fluid (261 µM), when compared with the delta6Li value of the subducted sediment section (-11 per mil), suggests that the deep source fluid originates from mineral fluid dehydration and transformation reactions at temperatures of 100 to 150°C, consistent with the temperature range of the up-dip seismogenic zone and of transformation of smectite to illite. The distribution of Li and its isotopes in the underthrust section are similar to those at the reference site, indicating near complete subduction of the incoming sediments and that early dewatering of the underthrust sediments occurs predominantly by lateral flow into the ocean. The hemipelagic clay-rich sediment section of the subducting plate carries most of the Li into this subduction zone, and the pelagic diatomaceous and nannofossil calcareous oozes contain little Li. The Li isotopes of both the clay-rich hemipelagic sediments and of the pelagic oozes are, however, similar, with delta6Li values of -9 to -12 per mil. The observations that (1) the delta6Li values of the underthrust sediments are distinctly lower than that of the mantle, and (2) the lavas of the Costa Rican volcanoes are enriched in Li and 7Li, provide an approximation of the contribution of the subducted sediments to the arc volcanoes. A first order mass balance calculation suggests that approximately half of the Li flux delivered by subducted sediments and altered oceanic crust into the Middle American Trench is recycled to the Costa Rican arc and at most a quarter of sedimentary Li is returned into the ocean through thrust faults, primarily the décollement thrust.

(Table 1) I, Br, B, d11B and Mn concentrations in ODP Site 131-808

This study of the interstitial water concentration-depth distributions of iodide, bromide, boron, d11B, and dissolved organic carbon, as represented by absorbance at 325 nm (yellow substance: YS) and laser-induced fluorescence (LIF), is a follow-up of the extensive shipboard program of interstitial water analysis during ODP Leg 131. Most of the components studied are associated with processes involving the diagenesis of organic matter in these sediments. Three zones of the sediment column are discussed separately because of the different processes involved in causing concentration changes: 1. The upper few hundreds of meters: In this zone, characterized by very high sedimentation rates (>1200 m/m.y.), interstitial waters show very sharp increases in alkalinity, ammonia, iodide, bromide, YS, and LIF, mainly as a result of the diagenesis of organic carbon; 2. Whereas below 200 mbsf concentration gradients all show a decreasing trend, the zone at ~ 365 mbsf is characterized by concentration reversals, mainly due to the recent emplacement of deeper sediments above this depth as a result of thrust-faulting; 3. The décollement zone (945-964 mbsf) is characterized by concentration anomalies in various constituents (bromide, boron, d11B, manganese, LIF). These data are interpreted as resulting from an advective input of fluids along the zone of décollement as recent as ~ 200 ka. Possibly periodic inputs of anomalous fluids still seem to occur along this décollement zone.

Diatom biostratigraphy of sediments from ODP Leg 170 sites

The following data paper summarizes diatom biostratigraphic data from sediments drilled in the Costa Rica accretionary wedge during Ocean Drilling Program Leg 170. Quaternary through lower Miocene diatom zones characteristic of the equatorial Pacific region are recognized in the reference section, Site 1039, which was drilled on the downgoing Cocos plate. At Sites 1040-1043, where the recovered silty clay units are primarily wedge and apron sediments that overlie the underthrust sections, diatoms are generally low in abundance, and complete zonation of the cores was not possible above the décollement surface.

(Table 1) Sr-isotope values of representative sediment samples from ODP sites at the Costa Rica margin

Two sealed borehole hydrologic observatories (CORKs) were installed in two active hydrogeochemical systems at the Costa Rica subduction zone to investigate the relationship between tectonics, fluid flow, and fluid composition. The observatories were deployed during Ocean Drilling Program (ODP) Leg 205 at Site 1253, ~ 0.2 km seaward of the trench, in the upper igneous basement, and at Site 1255, ~ 0.5 km landward of the trench, in the décollement. Downhole instrumentation was designed to monitor formation fluid flow rates, composition, pressure, and temperature. The two-year records collected by this interdisciplinary effort constitute the first co-registered hydrological, chemical, and physical dataset from a subduction zone, providing critical information on the average and transient state of the subduction thrust and upper igneous basement. The continuous records at ODP Site 1253 show that the uppermost igneous basement is highly permeable hosting an average fluid flow rate of 0.3 m/yr, and indicate that the fluid sampled in the basement is a mixture between seawater (~ 50%) and a subduction zone fluid originating within the forearc (~ 50%). These results suggest that the uppermost basement serves as an efficient pathway for fluid expelled from the forearc that should be considered in models of subduction zone hydrogeology and deformation. Three transients in fluid flow rates were observed along the décollement at ODP Site 1255, two of which coincided with stepwise increases in formation pressure. These two transients are the result of aseismic slip dislocations that propagated up-dip from the seismogenic zone over the course of ~ 2 weeks terminating before reaching ODP Site 1255 and the trench. The nature and temporal behavior of strain and the associated hydrological response during these slow slip events may be an analog for the response of the seaward part of the subduction prism during or soon after large subduction zone earthquakes.

Phospholipid fatty acid and quinone composition of sediments from ODP Leg 205 sites

Assessing the habitability of deep-sea sediments undergoing compaction, compression, and subduction at convergent margins adds to our understanding of the limits of the terrestrial biosphere. In this work, we report exploratory biomarker data on sediments obtained at Ocean Drilling Program (ODP) Sites 1253, 1254, and 1255 during drilling at the Costa Rica subduction trench and forearc sedimentary wedge. The samples selected for postcruise biomarker analyses were located within intervals of potentially enhanced fluid flow within the décollement and sedimentary wedge fault zones (Sites 1254 and 1255) and within basal carbonates at the reference site (Site 1253). The passage of fluids that are geochemically distinct from ambient interstitial water provides a disequilibrium setting that may enhance habitability. Biomarker data show low levels of microbial biomass in subseafloor sediments sampled at the Costa Rica convergent margin as deep as ~370 meters below seafloor.