Extensional tectonics and two-stage crustal accretion at oceanic transform faults.

Oceanic transform faults are seismically and tectonically active plate boundaries1 that leave scars-known as fracture zones-on oceanic plates that can cross entire ocean basins2. Current descriptions of plate tectonics assume transform faults to be conservative two-dimensional strike-slip boundaries1,3, at which lithosphere is neither created nor destroyed and along which the lithosphere cools and deepens as a function of the age of the plate4. However, a recent compilation of high-resolution multibeam bathymetric data from 41 oceanic transform faults and their associated fracture zones that covers all possible spreading rates shows that this assumption is incorrect. Here we show that the seafloor along transform faults is systemically deeper (by up to 1.6 kilometres) than their associated fracture zones, in contrast to expectations based on plate-cooling arguments. Accretion at intersections between oceanic ridges and transform faults seems to be strongly asymmetric: the outside corners of the intersections show shallower relief and more extensive magmatism, whereas the inside corners have deep nodal basins and seem to be magmatically starved. Three-dimensional viscoplastic numerical models show that plastic-shear failure within the deformation zone around the transform fault results in the plate boundary experiencing increasingly oblique shear at increasing depths below the seafloor. This results in extension around the inside corner, which thins the crust and lithosphere at the transform fault and is linked to deepening of the seafloor along the transform fault. Bathymetric data suggest that the thinned transform-fault crust is augmented by a second stage of magmatism as the transform fault intersects the opposing ridge axis. This makes accretion at transform-fault systems a two-stage process, fundamentally different from accretion elsewhere along mid-ocean ridges.

Broad fault zones enable deep fluid transport and limit earthquake magnitudes.

Constraining the controlling factors of fault rupture is fundamentally important. Fluids influence earthquake locations and magnitudes, although the exact pathways through the lithosphere are not well-known. Ocean transform faults are ideal for studying faults and fluid pathways given their relative simplicity. We analyse seismicity recorded by the Passive Imaging of the Lithosphere-Asthenosphere Boundary (PI-LAB) experiment, centred around the Chain Fracture Zone. We find earthquakes beneath morphological transpressional features occur deeper than the brittle-ductile transition predicted by simple thermal models, but elsewhere occur shallower. These features are characterised by multiple parallel fault segments and step overs, higher proportions of smaller events, gaps in large historical earthquakes, and seismic velocity structures consistent with hydrothermal alteration. Therefore, broader fault damage zones preferentially facilitate fluid transport. This cools the mantle and reduces the potential for large earthquakes at localized barriers that divide the transform into shorter asperity regions, limiting earthquake magnitudes on the transform.

Lower oceanic crust formed by in situ melt crystallisation revealed by seismic layering.

Oceanic crust forms at mid-ocean spreading centres through a combination of magmatic and tectonic processes, with the magmatic processes creating two distinct layers: the upper and the lower crust. While the upper crust is known to form from lava flows and basaltic dikes based on geophysical and drilling results, the formation of the gabbroic lower crust is still debated. Here we perform a full waveform inversion of wide-angle seismic data from relatively young (7-12-million-year-old) crust formed at the slow spreading Mid-Atlantic Ridge. The seismic velocity model reveals alternating, 400-500 m thick, high and low velocity layers with ±200 m/s velocity variations, below ~2 km from the oceanic basement. The uppermost low-velocity layer is consistent with hydrothermal alteration, defining the base of extensive hydrothermal circulation near the ridge axis. The underlying layering supports that the lower crust is formed through the intrusion of melt as sills at different depths, that cool and crystallise in situ. The layering extends up to 5-15 km distance along the seismic profile, covering 300,000-800,000 years, suggesting that this form of lower crustal accretion is a stable process.

Megathrust reflectivity reveals the updip limit of the 2014 Iquique earthquake rupture.

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.

Evolution of the Crustal and Upper Mantle Seismic Structure From 0-27 Ma in the Equatorial Atlantic Ocean at 2° 43'S.

We present seismic tomographic results from a unique seismic refraction and wide-angle survey along a 600 km long flow-line corridor of oceanic lithosphere ranging in age from 0 to 27 Ma in the equatorial Atlantic Ocean at 2° 43'S. The velocities in the crust near the ridge axis rapidly increase in the first 6 Myr and then change gradually with age. The upper crust (Layer 2) thickness varies between 2 and 2.4 km with an average thickness of 2.2 km and the crustal thickness varies from 5.6 to 6 km along the profile with an average crustal thickness of 5.8 km. At some locations, we observe negative velocity anomalies (∼-0.3 km/s) in the lower crust which could be either due to chemical heterogeneity in gabbroic rocks and/or the effects of fault related deformation zones leading to an increase in porosities up to 1.6% depending on the pore/crack geometry. The existence of a low velocity anomaly beneath the ridge axis suggests the presence of partial melt (∼1.3%) in the lower crust. Upper mantle velocities also remain low (∼7.8 km/s) from ridge axis up to 5 Ma, indicating a high temperature regime associated with mantle melting zone underneath. These results suggest that the evolution of the crust and uppermost mantle at this location occur in the first 10 Ma of its formation and then remains unchanged. Most of the structures in the older crust and upper mantle are fossilized structures and could provide information about past processes at ocean spreading centers.

Large slip, long duration, and moderate shaking of the Nicaragua 1992 tsunami earthquake caused by low near-trench rock rigidity.

Large earthquake ruptures propagating up to areas close to subduction trenches are infrequent, but when they occur, they heavily displace the ocean seafloor originating destructive tsunamis. The current paradigm is that the large seafloor deformation is caused by local factors reducing friction and increasing megathrust fault slip, or prompting the activation of ancillary faults or energy sources. As alternative to site-specific models, it has been proposed that large shallow slip could result from depth-dependent rock rigidity variations. To confront both hypotheses, here, we map elastic rock properties across the rupture zone of the M S7.0-M W7.7 1992 Nicaragua tsunami earthquake to estimate a property-compatible finite fault solution. The obtained self-consistent model accounts for trenchward increasing slip, constrains stress drop, and explains key tsunami earthquake characteristics such as long duration, high-frequency depletion, and magnitude discrepancy. The confirmation that these characteristics are all intrinsic attributes of shallow rupture opens new possibilities to improve tsunami hazard assessment.

Discovery of flat seismic reflections in the mantle beneath the young Juan de Fuca Plate.

Crustal properties of young oceanic lithosphere have been examined extensively, but the nature of the mantle lithosphere underneath remains elusive. Using a novel wide-angle seismic imaging technique, here we show the presence of two sub-horizontal reflections at ∼11 and ∼14.5 km below the seafloor over the 0.51-2.67 Ma old Juan de Fuca Plate. We find that the observed reflectors originate from 300-600-m-thick layers, with an ∼7-8% drop in P-wave velocity. They could be explained either by the presence of partially molten sills or frozen gabbroic sills. If partially molten, the shallower sill would define the base of a thin lithosphere with the constant thickness (11 km), requiring the presence of a mantle thermal anomaly extending up to 2.67 Ma. In contrast, if these reflections were frozen melt sills, they would imply the presence of thick young oceanic lithosphere (20-25 km), and extremely heterogeneous upper mantle.

Seismic Crustal Structure and Morphotectonic Features Associated With the Chain Fracture Zone and Their Role in the Evolution of the Equatorial Atlantic Region.

Oceanic transform faults and fracture zones (FZs) represent major bathymetric features that keep the records of past and present strike-slip motion along conservative plate boundaries. Although they play an important role in ridge segmentation and evolution of the lithosphere, their structural characteristics, and their variation in space and time, are poorly understood. To address some of the unknowns, we conducted interdisciplinary geophysical studies in the equatorial Atlantic Ocean, the region where some of the most prominent transform discontinuities have been developing. Here we present the results of the data analysis in the vicinity of the Chain FZ, on the South American Plate. The crustal structure across the Chain FZ, at the contact between ∼10 and 24 Ma oceanic lithosphere, is sampled along seismic reflection and refraction profiles. We observe that the crustal thickness within and across the Chain FZ ranges from ∼4.6-5.9 km, which compares with the observations reported for slow-slipping transform discontinuities globally. We attribute this presence of close to normal oceanic crustal thickness within FZs to the mechanism of lateral dike propagation, previously considered to be valid only in fast-slipping environments. Furthermore, the combination of our results with other data sets enabled us to extend the observations to morphotectonic characteristics on a regional scale. Our broader view suggests that the formation of the transverse ridge is closely associated with a global plate reorientation that was also responsible for the propagation and for shaping lower-order Mid-Atlantic Ridge segmentation around the equator.

Earthquake crisis unveils the growth of an incipient continental fault system.

Large continental faults extend for thousands of kilometres to form boundaries between rigid tectonic blocks. These faults are associated with prominent topographic features and can produce large earthquakes. Here we show the first evidence of a major tectonic structure in its initial-stage, the Al-Idrissi Fault System (AIFS), in the Alboran Sea. Combining bathymetric and seismic reflection data, together with seismological analyses of the 2016 Mw 6.4 earthquake offshore Morocco - the largest event ever recorded in the area - we unveil a 3D geometry for the AIFS. We report evidence of left-lateral strike-slip displacement, characterise the fault segmentation and demonstrate that AIFS is the source of the 2016 events. The occurrence of the Mw 6.4 earthquake together with historical and instrumental events supports that the AIFS is currently growing through propagation and linkage of its segments. Thus, the AIFS provides a unique model of the inception and growth of a young plate boundary fault system.

The Alboran volcanic-arc modulated the Messinian faunal exchange and salinity crisis.

What process triggered the Mediterranean Sea restriction remains debated since the discovery of the Messinian Salinity Crisis (MSC). Recent hypotheses infer that the MSC initiated after the closure of the Atlantic-Mediterranean Betic and Rifean corridors, being modulated through restriction at the Gibraltar Strait. These hypotheses however, do not integrate contemporaneous speciation patterns of the faunal exchange between Iberia and Africa and several geological features like the evaporite distribution. Exchange of terrestrial biota occurred before, during and after the MSC, and speciation models support an exchange path across the East Alborán basin (EAB) located a few hundreds of km east of the Gibraltar Strait. Yet, a structure explaining jointly geological and biological observations has remained undiscovered. We present new seismic data showing the velocity structure of a well-differentiated 14-17-km thick volcanic arc in the EAB. Isostatic considerations support that the arc-crust buoyancy created an archipelago leading to a filter bridge across the EAB. Sub-aerial erosional unconformities and onlap relationships support that the arc was active between ~10-6 Ma. Progressive arc build-up leading to an archipelago and its later subsidence can explain the extended exchange of terrestrial biota between Iberia and Africa (~7-3 Ma), and agrees with patterns of biota speciation and terrestrial fossil distribution before the MSC (10-6.2 Ma). In this scenario, the West Alboran Basin (WAB) could then be the long-postulated open-marine refuge for the Mediterranean taxa that repopulated the Mediterranean after the MSC, connected to the deep restricted Mediterranean basin through a sill at the Alboran volcanic arc archipelago.