Description: Water and carbon are transferred from the ocean to the mantle in a process that alters mantle peridotite to create serpentinite and supports diverse ecosystems1. Serpentinized mantle rocks are found beneath the sea floor at slow- to ultraslow-spreading mid-ocean ridges1 and are thought to be present at about half the world’s rifted margins2, 3. Serpentinite is also inferred to exist in the downgoing plate at subduction zones4, where it may trigger arc magmatism or hydrate the deep Earth. Water is thought to reach the mantle via active faults3, 4. Here we show that serpentinization at the rifted continental margin offshore from western Spain was probably initiated when the whole crust cooled to become brittle and deformation was focused along large normal faults. We use seismic tomography to image the three-dimensional distribution of serpentinization in the mantle and find that the local volume of serpentinite beneath thinned, brittle crust is related to the amount of displacement along each fault. This implies that sea water reaches the mantle only when the faults are active. We estimate the fluid flux along the faults and find it is comparable to that inferred for mid-ocean ridge hydrothermal systems. We conclude that brittle processes in the crust may ultimately control the global flux of sea water into the Earth.

Description: We present the first application of Stochastic Heterogeneity Mapping based on the band-limited von Kármán function to a seismic reflection stack of a Mediterranean water eddy (meddy), a large salt lens of Mediterranean water. This process extracts two stochastic parameters directly from the reflectivity field of the seismic data: the Hurst number, which ranges from 0 to 1, and the correlation length (scale length). Lower Hurst numbers represent a richer range of high wavenumbers and correspond to a broader range of heterogeneity in reflection events. The Hurst number estimate for the top of the meddy (0.39) compares well with recent theoretical work, which required values between 0.25 and 0.5 to model internal wave surfaces in open ocean conditions based on simulating a Garrett-Munk spectrum (GM76) slope of −2. The scale lengths obtained do not fit as well to seismic reflection events as those used in other studies to model internal waves. We suggest two explanations for this discrepancy: (1) due to the fact that the stochastic parameters are derived from the reflectivity field rather than the impedance field the estimated scale lengths may be underestimated, as has been reported; and (2) because the meddy seismic image is a two-dimensional slice of a complex and dynamic three-dimensional object, the derived scale lengths are biased to the direction of flow. Nonetheless, varying stochastic parameters, which correspond to different spectral slopes in the Garrett-Munk spectrum (horizontal wavenumber spectrum), can provide an estimate of different internal wave scales from seismic data alone. We hence introduce Stochastic Heterogeneity Mapping as a novel tool in physical oceanography.

Description: The nature of the Ionian Sea crust has been the subject of scientific debate for more than 30 years, mainly because seismic imaging of the deep crust and upper mantle of the Ionian Abyssal Plain (IAP) has not been conclusive to date. The IAP is sandwiched between the Calabrian and Hellenic subduction zones in the central Mediterranean. To univocally confirm the proposed oceanic nature of the IAP crust as a remnant of the Tethys ocean and to confute its interpretation as a strongly thinned part of the African continental crust, a NE-SW oriented 131 km long seismic refraction and wide-angle reflection profile consisting of eight ocean bottom seismometers and hydrophones was acquired in 2014. A P-wave velocity model developed from travel time forward modelling is refined by gravimetric data and synthetic modelling of the seismic data. A roughly 6km thick crust with velocities ranging from 5.1km/s to 7.2km/s, top to bottom, can be traced throughout the IAP. In the vicinity of the Medina Seamounts at the southern IAP boundary, the crust thickens to about 9km and seismic velocities decrease to 6.8km/s at the crust-mantle boundary. The seismic velocity distribution and depth of the crust-mantle boundary in the IAP document its oceanic nature, and support the interpretation of the IAP as a remnant of the Tethys oceanic lithosphere formed during the Permian and Triassic period.

Description: Spreading processes associated with slow-spreading ridges are a complex interplay of volcanic accretion and tectonic dismemberment of the oceanic crust, resulting in an irregular seafloor morphology made up of blocks created by episodes of intense volcanic activity or tectonic deformation. These blocks undergo highly variable evolution, such as tilts or dissection by renewed tectonic extension, depending on their positions with respect to the spreading axis, core complexes, detachment or transform faults. Here, we use near-seafloor magnetic and bathymetric data and seismic profiles collected over the TAG Segment of the Mid-Atlantic Ridge to constrain the tectonic evolution of these blocks. Our study reveals that the presence and evolution of oceanic core complexes play a key role in triggering block movements. The deep subvertical detachment fault roots on the plate boundary, marked by a thermal anomaly and transient magma bodies. Thermal and magmatic variations control the structure and morphology of the seafloor above the subhorizontal detachment surface, occasionally leading to relocating the detachment.

Description: Offshore south central Chile (35° S–42° S), the morphology of the lowermost continental slope and trench floor witnesses a voluminous submarine mass-wasting event. The blocky slide body deposited in the Chile Trench at 73°46´ W 35°35´ S was targeted for study during RRS JAMES COOK Cruise JC23 and termed Reloca Slide. Its size of about 24 km3, its steep and high headscarp, the spatial distribution of slide deposits and the cohesive nature of major slide blocks make it interesting to address the issue of tsunami generation. We have obtained seismic reflection data that partly reveal the internal structure of the slide body. Gravity core samples were retrieved that will allow the slide to be dated and linked to the history of sedimentation and slope stability along this particular segment of the Chilean convergent margin. At present we assume a Holocene age for the sliding event.

Description: Accurate subsurface velocity models are crucial for geological interpretations based on seismic depth images. Seismic reflection tomography is an effective iterative method to update and refine a preliminary velocity model for depth imaging. Based on residual move-out analysis of reflectors in common image point gathers an update of the velocity is estimated by a ray-based tomography. To stabilize the tomography, several preconditioning strategies exist. Most critical is the estimation of the depth error to account for the residual move-out of the reflector in the common image point gathers. Because the depth errors for many closely spaced image gathers must be picked, manual picking is extremely time-consuming, human biased, and not reproducible. Data-driven picking algorithms based on coherence or semblance analysis are widely used for hyperbolic or linear events. However, for complex-shaped depth events, pure data-driven picking is difficult. To overcome this, the warping method named Non-Rigid Matching is used to estimate a depth error displacement field. Warping is used, e.g., to merge photographic images or to match two seismic images from time-lapse data. By calculating the displacements between an offset to its neighbouring offset in the common image point domain, a locally smooth-shaped displacement field is defined for each data sample. Depending on the complexity of the subsurface, sample tracking through the displacement field along predefined horizons or on a simple regular grid yields discrete depth error values for the tomography. The application to a multi-channel seismic line across the Sunda subduction zone offshore Lombok island, Indonesia, illustrates the approach and documents the advantages of the method to estimate a detailed velocity structure in a complex tectonic regime. By incorporating the warping scheme into the reflection tomography, we demonstrate an increase in the velocity resolution and precision by improving the data-driven accuracy of depth error picks with arbitrary shapes. This approach offers the possibility to use the full capacities of tomography and further leads to more accurate interpretations of complex geological structures.

Description: We resolve a previously unrecognized shallow subducting seamount from a re-processed multichannel seismic depth image crossing the 1994 M7.8 Java tsunami earthquake slip area. Seamount subduction is related to the uplift of the overriding plate by lateral shortening and vertical thickening, causing pronounced back-thrusting at the landward slope of the forearc high and the formation of splay faults branching off the landward flank of the subducting seamount. The location of the seamount in relation to the 1994 earthquake hypocentre and its co-seismic slip model suggests that the seamount acted as a seismic barrier to the up-dip co-seismic rupture propagation of this moderate size earthquake. The wrapping of the co-seismic slip contours around the seamount indicates that it diverted rupture propagation, documenting the control of forearc structures on seismic rupture.

Description: The updip limit of seismic rupture during a megathrust earthquake exerts a major control on the size of the resulting tsunami. Offshore Northern Chile, the 2014 Mw 8.1 Iquique earthquake ruptured the plate boundary between 19.5° and 21°S. Rupture terminated under the mid-continental slope and did not propagate updip to the trench. Here, we use state-of-the-art seismic reflection data to investigate the tectonic setting associated with the apparent updip arrest of rupture propagation at 15 km depth during the Iquique earthquake. We document a spatial correspondence between the rupture area and the seismic reflectivity of the plate boundary. North and updip of the rupture area, a coherent, highly reflective plate boundary indicates excess fluid pressure, which may prevent the accumulation of elastic strain. In contrast, the rupture area is characterized by the absence of plate boundary reflectivity, which suggests low fluid pressure that results in stress accumulation and thus controls the extent of earthquake rupture. Generalizing these results, seismic reflection data can provide insights into the physical state of the shallow plate boundary and help to assess the potential for future shallow rupture in the absence of direct measurements of interplate deformation from most outermost forearc slopes.