Slow Seismic Velocities Research Articles

Inclined bending seismic reflection layer in the crust illuminated by distributed fibre-optic-sensing measurements in western Japan

Distributed acoustic sensing (DAS), which enables a single fibre-optic cable to function as multiple sensors, is a technique to measure the strain rate distributed along the cable. This technique is applied to record ground motions in western Japan via a 50 km-long fibre-optic cable beneath a national road. The measured values are strain changes along the cable every 5 m, corresponding to 9788 sensor deployments. This high-density measurement along the long cable successfully recorded the 2021 M2.8 and M3.2 earthquakes that occurred in the crust within the distance of the cable in southern Kyoto. The direct S waves were followed by seismic waves approximately 8–14 s later, which were reflected by lower crustal structures. These waveforms were previously reported by observing many earthquakes via multiple seismometers, but the DAS observations clearly illuminate reflected wavefields from single earthquake observations for the first time. The numerical simulation of the strain-rate wavefields of these earthquakes reveals the existence of a north-dipping thin layer with a slow seismic velocity in the lower crust, which becomes steeper in the shallower part. This layer might represent the path of slab-derived fluid to the shallow fault zone.

Seismic Imaging Beneath Cascadia Shows Shallow Mantle Flow Patterns Guide Lower Mantle Upwellings

AbstractThe mantle transition zone (MTZ) plays an important role in modulating material transport between the upper mantle and the lower mantle. Constraining this transport is essential for understanding geochemical reservoirs, hydration cycles, and the evolution of the Earth. Slabs and hotspots are assumed to be the dominant locations of transport. However, the degree of material transport in other areas is debated. We apply P‐to‐S receiver functions to an amphibious data set from Cascadia to image the MTZ discontinuities beneath mid‐ocean ridges, a hotspot, and a subduction zone. We find a MTZ thinned by 10 ± 6 km beneath the ridges and by 8 ± 4 km beneath the base of the slab, closely resembling the 660 discontinuity topography. Depressions on the 410 discontinuity are smaller, 5 ± 2 km on average, focused in the north and the south and accompanied by supra‐410 discontinuity melt phases. The depressions occur away from locations of uplifted 660 discontinuity, but near slow seismic velocity anomalies imaged in the upper mantle. This suggests lower mantle upwellings occur beneath ridges and beneath the base of slabs but stall in the transition zone, with upper mantle convection determining upward material transport from the transition zone. Therefore, upper mantle dynamics play a larger role in determining transfer than typically assumed.

Mapping structures on the core–mantle boundary using Sdiff postcursors: Part I. Method and Validation

SUMMARY Ultra-low velocity zones (ULVZs) are patches of extremely slow seismic velocities on the core–mantle boundary (CMB). Here, we target them using the postcursors to S core-diffracted phases (Sdiff) caused by ULVZs. We use traveltimes of these postcursors to make probabilistic maps using a reversible-jump Markov chain Monte Carlo inversion setup. For the forward model, we extend 2-D wave front tracking (2DWT) software, previously developed for surface wave multipathing studies, to the CMB. The 2DWT is able to model the full multipathing behaviour of Sdiff postcursors and compute arrival times for a given ULVZ input velocity structure on the order of a few CPU seconds, as opposed to 100s of CPU hours required for 3-D full waveform synthetics. We validate the method using synthetic data sets produced by the 2DWT, as well as 3-D full waveform synthetics, using a parametrisation formed from a collection of ellipses. We also test idealistic data coverage versus a case of more realistic coverage. We show ULVZ size and velocity reduction can typically be well recovered, and our maps show the inherent trade-off between these parameters around the edge of the ULVZ. Our method cannot directly constrain ULVZ height; tests show that it underestimates ULVZ velocity reductions and overestimates ellipticity for thinner ULVZs due to neglecting mantle effects.

Mineralogy, fabric and deformation domains in D″ across the southwestern border of the African LLSVP

SUMMARYRecent advances in seismic anisotropy studies that jointly use reflections and shear wave splitting have proven to place tight constraints on the plausible anisotropic and deformation scenarios in the D″ region. We apply this novel methodology to a large area of the D″ region beneath the South Atlantic, in proximity to and within the African large low seismic velocity province (LLSVP). This area of the mantle is characterized by a transition from fast to slow seismic velocity anomalies and it is thought to be the location of deep-seated plumes responsible for hotspot volcanism. Attempting to probe mantle composition and deformation along the LLSVP borders may provide key information on mantle dynamics. By analysing seismic phases sampling this region, we detect a D″ discontinuity over a large area beneath the South Atlantic, with inferred depth ranges ∼170 to ∼240 km above the core–mantle boundary. We find evidence for a D″ reflector within the area of the LLSVP. Shear wave splitting observations suggest that anisotropy is present in this region of the mantle, in agreement with previous studies that partially sampled this region. We model the observations considering lattice- and shape-preferred orientation of materials expected in the D″ region. A regional variation of mineralogy, phase transition boundaries, and deformation direction is required to explain the data. We infer two distinct domains of mineralogy and deformation: aligned post-perovskite outside the LLSVP and aligned bridgmanite within the LLSVP. The scenario depicted by this study agrees well with the current hypotheses for the composition of the LLSVP and with the prevalence of vertical deformation directions expected to occur along the LLSVPs borders.

Dense mantle flows periodically spaced below ocean basins

Understanding mantle flow is key to elucidate how deep Earth dynamics relate to tectonics at the global scale. The convective mass transport is reflected in lateral variations of the gravity field, seismic velocities, as well as deformations of the Earth's surface. Yet, upper to mid-mantle dynamics have been difficult to constrain at the medium scales of thousands of km. Here, we analyze the second-order horizontal derivatives of seafloor topography and of the gravity potential over the Pacific and Northern Indian ocean basins, and provide evidence for periodic undulations of 1600-2000 km wavelength in both signals, elongated along the direction of absolute plate motion. We investigate potential crustal and lithospheric sources and show that at least part of this signal must originate below the lithosphere, with alignments of sub-lithospheric upper mantle mass excess below seafloor lows. Furthermore, we find that these alignments coincide geographically over wide areas with similarly periodic slow seismic velocity fingers located at upper mantle depths. These two fields may thus record an intermediate scale of mantle convection below ocean basins, which cannot be explained by purely thermal convection and requires instead lateral variations in composition in the upper mantle. Elucidating the nature of the detected mass excess sources coincident with the slow seismic velocities calls for a joint dynamical modeling of all observations in a thermo-chemical context.

Dynamics of the North American Plate: Large‐Scale Driving Mechanism From Far‐Field Slabs and the Interpretation of Shallow Negative Seismic Anomalies

AbstractWith a small fraction of marginal subduction zones, the driving mechanism for the North American plate motion is in debate. We construct global mantle flow models simultaneously constrained by geoid and plate motions to investigate the driving forces for the North American plate motion. By comparing the model with only near‐field subducting slabs and that with global subducting slabs, we find that the contribution to the motion of the North American plate from the near‐field Aleutian, central American, and Caribbean slabs is small. In contrast, other far‐field slabs, primarily the major segments around western Pacific subduction margins, provide the dominant large‐scale driving forces for the North American plate motion. The coupling between far‐field slabs and the North American plate suggests a new form of active plate interactions within the global self‐organizing plate tectonic system. We further evaluate the extremely slow seismic velocity anomalies associated with the shallow partial melt around the southwestern North America. Interpreting these negative seismic shear‐velocity anomalies as purely thermal origin generates considerably excessive resistance to the North American plate motion. A significantly reduced velocity‐to‐density scaling for these negative seismic shear‐velocity anomalies must be incorporated into the construction of the buoyancy field to predict the North American plate motion. We also examine the importance of lower mantle buoyancy including the ancient descending Kula‐Farallon plates and the active upwelling below the Pacific margin of the North American plate. Lower mantle buoyancy primarily affects the amplitudes, as opposed to the patterns of both North American and global plate motions.

Spectral Element Modeling of Acoustic to Seismic Coupling Over Topography

AbstractAcoustic waves in the atmosphere are commonly recorded on seismometers as they couple into the ground. These signals, here called ground coupled airwaves, are not commonly considered in numerical modeling of infrasound propagation, which often assumes a rigid unmeshed boundary. Starting from an analytically‐tractable spherical wave model, we analyze how the coupling of an acoustic wave into a planar elastic halfspace changes over a wide range of scenarios. We use energy admittance to quantify acoustic to seismic coupling over both a planar elastic halfspace and meshed topography. Our spectral element and analytic calculations have different maxima as a function of incidence angle, with very high admittance values for near‐vertical incidence (maximum ≈78%). Energy admittance calculations at shallow incidence angles are much smaller (less than 1%). In simulations over the complex topography of Sakurajima Volcano, we attribute the variable spatial pattern of energy admittance to changes in earth parameters between each model. The observed pressure difference over the simulated 15 km region appears to be . While this value is relatively small, the cumulative addition over 100s of km and multiple acoustic bounce points may be significant. Acoustic to seismic coupling along the propagation path may bias long distance yield estimates, particularly when infrasound propagates over regions with steep topography or particularly slow seismic velocities, such as alluvial planes.

The seismic structure of the Antarctic upper mantle

Abstract The deployment of seismic stations and the development of ambient noise tomography as well as new analysis methods provide an opportunity for higher-resolution imaging of Antarctica. Here we review recent seismic structure models and describe their implications for the dynamics and history of the Antarctic upper mantle. Results show that most of East Antarctica is underlain by continental lithosphere to depths of approximately 200 km. The thickest lithosphere is found in a band 500–1000 km inboard from the Transantarctic Mountains, representing the continuation of cratonic lithosphere with Australian affinity beneath the ice. Dronning Maud Land and the Lambert Graben show much thinner lithosphere, consistent with Phanerozoic lithospheric disruption. The Transantarctic Mountains mark a sharp boundary between cratonic lithosphere and the warmer upper mantle of West Antarctica. In the southern Transantarctic Mountains, cratonic lithosphere has been replaced by warm asthenosphere, giving rise to Cenozoic volcanism and an elevated mountainous region. The Marie Byrd Land volcanic dome is underlain by slow seismic velocities extending through the transition zone, consistent with a mantle plume. Slow-velocity anomalies beneath the coast from the Amundsen Sea Embayment to the Antarctic Peninsula are likely to result from upwelling of warm asthenosphere during subduction of the Antarctic–Phoenix spreading centre.

Crustal and Upper‐Mantle Structure Below Central and Southern Mexico

AbstractThe Mexican subduction zone is a complex convergent margin characterized by a slab with an unusual morphology, an abnormal location of the volcanic arc, and the existence of slow earthquakes (slow slip events and tectonic tremors). The number of seismic imaging studies of the Mexican subduction zone has increased over the past decade. These studies have revealed the seismic wave velocity structure beneath Central and Southern Mexico. However, the existing tomographic models are either confined to specific areas or have coarse resolution. As a consequence, they do not capture the lateral and finer‐scale variations in the mantle and the crust structure. Here, we fit frequency‐dependent traveltime differences between observed and synthetic seismograms in a three‐dimensional model of Central and Southern Mexico to constrain crustal and upper mantle seismic velocity structures jointly. We use ∼3,300 seismic records, filtered between 5 and 50 s, from 74 regional earthquakes. We perform point spread function tests to assess the resolution and trade‐offs between model parameters. Our tomographic model reveals several seismic wave velocity features. These features include the geometry of the Cocos slab, the partial melting zone beneath the Trans‐Mexican Volcanic Belt, and the Yucatan slab diving below Los Tuxtlas Volcanic Field that are consistent with previous studies. We also identify horizontal structures in the lower crust and the shallow upper mantle that were not previously identified. These structures include slow seismic wave velocity anomalies that may be correlated with ultraslow seismic velocity zones and high conductivity regions, where slow earthquakes have been identified.

Continental Drift with Deep Cratonic Roots

The influence of the continental lithosphere and its root (or keel) on the continental drift of Earth is a key element in the history of plate tectonics. Previous geodynamic studies of mantle flow suggested that the cratonic root is moderately mechanically coupled with the underlying mantle, and stable continental drift on Earth's timescales occurs when the effective viscosity contrast between the continental lithosphere and the underlying mantle is approximately 103. Both geodynamics and seismological studies indicate that mechanically weak mobile belts (i.e., orogenic or suture zones) that surround cratons may play a role in the longevity of the cratonic lithosphere over geologically long timescales (i.e., over 1,000 million years) because they act as a buffer region against the high-viscosity cratons. Low-viscosity asthenosphere, characterized by slow seismic velocities, reduces the basal drag force acting on the cratonic root, which may also contribute to the longevity of the cratonic lithosphere. ▪ The role of the continental lithosphere and its root on the continental drift is reviewed from recent geodynamic and seismological studies. ▪ The cratonic root is moderately mechanically coupled with the underlying mantle and deformed by mantle flow over geological timescales. ▪ Orogenic belts or suture zones that surround cratons act as a buffer to protect cratons and are essential for their longevity. ▪ Low-viscosity asthenosphere may reduce the basal drag acting on the cratonic root and also contribute to its stability and longevity.

Mantle dynamics of the Andean Subduction Zone from continent-scale teleseismic S-wave tomography

SUMMARY The Andean Subduction Zone is one of the longest continuous subduction zones on Earth. The relative simplicity of the two-plate system has makes it an ideal natural laboratory to study the dynamics in subduction zones. We measure teleseismic S and SKS traveltime residuals at >1000 seismic stations that have been deployed across South America over the last 30 yr to produce a finite-frequency teleseismic S-wave tomography model of the mantle beneath the Andean Subduction Zone related to the Nazca Plate, spanning from ∼5°N to 45°S and from depths of ∼130 to 1200 km. Within our model, the subducted Nazca slab is imaged as a fast velocity seismic anomaly. The geometry and amplitude of the Nazca slab anomaly varies along the margin while the slab anomaly continues into the lower mantle along the entirety of the subduction margin. Beneath northern Brazil, the Nazca slab appears to stagnate at ∼1000 km depth and extend eastward subhorizontally for >2000 km. South of 25°S the slab anomaly in the lower mantle extends offshore of eastern Argentina, hence we do not image if a similar stagnation occurs. We image several distinct features surrounding the slab including two vertically oriented slow seismic velocity anomalies: one beneath the Peruvian flat slab and the other beneath the Paraná Basin of Brazil. The presence of the latter anomaly directly adjacent to the stagnant Nazca slab suggests that the plume, known as the Paraná Plume, may be a focused upwelling formed in response to slab stagnation in the lower mantle. Additionally, we image a high amplitude fast seismic velocity anomaly beneath the Chile trench at the latitude of the Sierras Pampeanas which extends from ∼400 to ∼1000 km depth. This anomaly may be the remnants of an older, detached slab, however its relationship with the Nazca–South America subduction zone remains enigmatic.

Multi-parameter imaging in the 2016 Mw 7.0 Kumamoto earthquake, Kyushu, Japan and its implications for fluids in seismogenesis

Fluids associated with structural heterogeneities in a seismogenic layer play a crucial role in triggering earthquakes and controlling their rupture processes, but the mechanism of fluid intrusion and the relationship with earthquake initiation are unclear. This study presents new seismic evidence of an earthquake mostly triggered by fluid intrusion associated with magmatism in Kyushu, SW Japan. Seismic imaging revealed two zones with slow seismic velocities beneath the Unzen and Aso Volcanoes with high Poisson's ratio, extending into the rupture belt of the 2016 Mw 7.0 Kumamoto earthquake. The following features, particularly the high Poisson's ratio anomalies, suggest the presence of fluids in the source area ascending from nearby magmatism. Anomalous crack density and saturation rate imply that the fractures/pores are not fully saturated with fluids in the crust of the Beppu-Shimabara Graben, possibly providing a path for fluids to penetrate into the source area, either through highly fracture-damaged faults or via incompletely saturated pores/cracks. Fluid penetrating into the hypocenter areas from deep magma increases the pore pressure and decreases the Coulomb failure function. We therefore speculate that a system such as this, with high pore pressure and fluids in the host faults might lower the threshold for fault failure, thus increasing the likelihood of failure of the mainshock fault. The results of this study suggest that the presence of fractures/pores within extensional rift zones around the world might be critical for fluid migration, and contribute to the initiation, development and migration of seismic activity.

Short‐ and Long‐Term Variations in the Reykjanes Geothermal Reservoir From Seismic Noise Interferometry

AbstractThe Reykjanes Geothermal System (RGS) is a high‐temperature geothermal system located on the Reykjanes peninsula, a transtensional plate‐boundary zone located on the southwestern tip of Iceland. The area is characterized by high seismicity, recent volcanism, and high‐temperature geothermal fields. We use seismic noise records from April 2014 to August 2015 to study stress changes and potential deformation of the subsurface caused by injection and production operations at RGS through seismic interferometry. We retrieve continuous time series of waveform similarity values and seismic velocity changes during this period. The S‐transform of the similarity values allows us to clearly identify three variations in the mechanical properties of the Reykjanes peninsula related to rapid changes of RGS production. In addition, we observe a slow seismic velocity decrease of 0.36%/year in the reservoir due to the water deficit and seasonal variations associated with the energy production demand.

Normal faulting and viscous buckling in the Tibetan Plateau induced by a weak lower crust

Flow of weak lower crust has been invoked to reconcile observed topographic gradients, uniform elevations, slow seismic velocity, and high conductivity measured in the Tibetan Plateau, with viscosity estimates of 1016–1021 Pa·s. Here we investigate the dynamic response resulting from a range of lower crust viscosities in a 3-D lithospheric-scale geodynamic simulation of the India–Eurasia collision zone to determine bounds of physically viable lower crustal strengths. We show that thickening of the plateau is accommodated through viscous buckling of the upper crust in response to lower crustal flow for a lower crustal viscosity on the order of 1020 Pa·s. This generates two east–west trending bands of surface subsidence and dilatation consistent with observed normal faulting and estimates of vertical velocity. These results suggest viscous buckling of the upper crust, induced by lower crustal flow from gravitational pressure gradients due to high topography, is responsible for the observed extension in Tibet.

Dynamics of microseismicity and its relationship with the active structures in the western Corinth Rift (Greece)

We analyse the complete earthquake archive of the western Corinth Rift using both cross-correlations between pairs of event waveforms and accurate differential traveltimes observed at common stations, in order to identify small-scale fault structures at depth. The waveform database was generated by the dense Corinth Rift Laboratory network and includes about 205 000 events between 2000 and 2015. Half of them are accurately relocated using double-difference techniques. The novelty of this relocated catalogue is the integration of the recent westernmost earthquakes due to the extension of the network in 2010 to the western extremity of the Corinth Rift and the consideration of the whole database over more than 15 yr. The total relocated seismicity exhibits well-defined clusters at the root of the main normal faults mainly between 5 and 10 km depth in the middle of the gulf and illuminates thin active structure planes dipping north about 20° under the northern coast. Some seismicity is observed in the footwall of the main active faults, along the West and East Helike faults. We also built a multiplet database based on waveform similarity taking into account cross-correlation coefficients weighted by signal-to-noise ratios. Short-term multiplets are concentrated in the middle of the gulf along the Kamarai fault system, in a 1–2 km thick layer at 6–8 km depth, interpreted as a highly fractured geological layer. They are often associated to slow seismic migration velocities occurring in this zone during strong swarm episodes and are thus likely to be triggered by pore pressure variations. On the other hand, most long-term and regular multiplets are located deeper (7–10 km), under the northern coast, within a layer less than 0.3 km thick. They occur at the border of nearly planar structures with low seismicity rate, which we identify as fault planes, and they may be explained by aseismic slip on the fault surface around them. This supports the existence of an immature structure growing downdip towards the north at the base of the active geological layer, which possibly connects to the ductile middle crust around 15 km depth, as suggested by the occurrence of deeper events in the continuity of the 1995–fault plane. The different migration velocities (from 0.05 km d−1 to several km d−1) highlighted during the western 2014–swarms indicate that both pore pressure and creep diffusion are operating in the fault zone. The fast migrations observed in the Psathopyrgos fault zone, where a slow slip event was detected by dilatometers in 2002, compare with that for creeping faults. To the west, from spatial distribution of events, we show that the Rion–Patras fault connecting the western extremity of the Corinth Rift fault system to the Patras Rift, is dipping around 60° north–west with a rake angle of −115°. Finally, we identified two new areas within the central active zone which may correspond to large scale, locked asperities on active fault surfaces, similar in size to the main asperity broken during the 1995, MW 6.3, Aigion earthquake.

Fat Plumes May Reflect the Complex Rheology of the Lower Mantle.

Recent tomographic imaging of the mantle below major hot spots shows slow seismic velocities extending down to the core‐mantle boundary, confirming the existence of mantle plumes. However, these plumes are much thicker than previously thought. Using new laboratory experiments and scaling laws, we show that thermal plumes developing in a visco‐plastic fluid present much larger diameters than plumes developing in a Newtonian fluid. Such a rheology requiring a yield stress is consistent with a lower mantle predominantly deforming by pure dislocation climb. Yield stress values between 1 and 10 MPa, implying dislocation densities between 108 and 1010 m−2, would be sufficient to reproduce the plumes morphology observed in tomographic images.

Southwestern Africa on the burner: Pleistocene carbonatite volcanism linked to deep mantle upwelling in Angola

The origin of intraplate carbonatitic to alkaline volcanism in Africa is controversial. A tectonic control, i.e., decompression melting associated with far-field stress, is suggested by correlation with lithospheric sutures, repeated magmatic cycles in the same areas over several million years, synchronicity across the plate, and lack of clear age progression patterns. Conversely, a dominant role for mantle convection is supported by the coincidence of Cenozoic volcanism with regions of lithospheric uplift, positive free-air gravity anomalies, and slow seismic velocities. To improve constraints on the genesis of African volcanism, here we report the first radiometric and isotopic results for the Catanda complex, which hosts the only extrusive carbonatites in Angola. Apatite (U-Th-Sm)/He and phlogopite 40Ar/39Ar ages of Catanda aillikite lavas indicate eruption at ca. 500–800 ka, more than 100 m.y. after emplacement of abundant kimberlites and carbonatites in this region. The lavas share similar high-μ (HIMU)–like Sr-Nd-Pb-Hf isotope compositions with other young mantle-derived volcanics from Africa (e.g., Northern Kenya Rift; Cameroon Line). The position of the Catanda complex in the Lucapa corridor, a long-lived extensional structure, suggests a possible tectonic control for the volcanism. The complex is also located on the Bie Dome, a broad region of fast Pleistocene uplift attributed to mantle upwelling. Seismic tomography models indicate convection of deep hot material beneath regions of active volcanism in Africa, including a large area encompassing Angola and northern Namibia. This is strong evidence that intraplate late Cenozoic volcanism, including the Catanda complex, resulted from the interplay between mantle convection and preexisting lithospheric heterogeneities.

The nature of subslab slow velocity anomalies beneath South America

AbstractSlow seismic velocity anomalies are commonly imaged beneath subducting slabs in tomographic studies, yet a unifying explanation for their distribution has not been agreed upon. In South America two such anomalies have been imaged associated with subduction of the Nazca Ridge in Peru and the Juan Fernández Ridge in Chile. Here we present new seismic images of the subslab slow velocity anomaly beneath Chile, which give a unique view of the nature of such anomalies. Slow seismic velocities within a large hole in the subducted Nazca slab connect with a subslab slow anomaly that appears correlated with the extent of the subducted Juan Fernández Ridge. The hole in the slab may allow the subslab material to rise into the mantle wedge, revealing the positive buoyancy of the slow material. We propose a new model for subslab slow velocity anomalies beneath the Nazca slab related to the entrainment of hot spot material.

A new approach to obtaining a 3D shear wave velocity model of the crust and upper mantle: An application to eastern Turkey

We present a new approach to the joint inversion of surface wave dispersion data and receiver functions by utilizing Common Conversion Point (CCP) stacking to reconcile the different sampling domains of the two datasets. Utilizing CCP stacking allows us to suppress noise in the data by waveform stacking, and correct for backazimuthal variations and complex crustal structure by mapping receiver functions back to their theoretical location. When applied to eastern Turkey, this approach leads to a higher resolution image of the subsurface and clearly delineates different tectonic features in eastern Turkey that were not apparent using other approaches. We observe that the slow seismic velocities near the Karliova Triple Junction correlate to moderate strain rates and high heat flow, which leads to a rheologically weak crust that has allowed for the upward propagation of Miocene and younger volcanics near the triple junction. We find seismically fast, presumably rigid blocks located in the southeastern Anatolian Plate and Arabian Plate are separated by a band of low shear wave velocities that correspond to the East Anatolian Fault Zone, which is consistent with the presence of fluids in the fault zone. We observe that the Arabian Plate has underthrust the Eurasian Plate as far as the northern boundary of the Bitlis Massif, which can explain the high exhumation rates in the Bitlis Massif as a result of slab break-off of the Arabian oceanic lithosphere. We also find a shallow (~33km) anomaly beneath eastern Turkey that we interpret as a localized wedge of mantle that was underthrust by a crustal fragment during the collision of Arabia and Eurasia. These observations are possible because of the high-resolution images obtained by combining common conversion point receiver function stacks with ambient noise dispersion data to create a data-driven three-dimensional shear wave velocity model.

Flat versus normal subduction zones: a comparison based on 3-D regional traveltime tomography and petrological modelling of central Chile and western Argentina (29°–35°S)

Our study compares the seismic properties between the flat and normal subduction regions in central Chile, to better understand the links between the slab geometry, surface deformation and the deeper structures. In comparison with previous studies, we show the most complete 3-D regional seismic tomography images for this region, in which we use (1) a larger seismic data set compiled from several short-term seismic catalogues, (2) a denser seismic array allowing a better resolution of the subduction zone from the trench to the backarc and into the upper ∼30 km of the slab and (3) a starting 1-D background velocity model specifically calculated for this region and refined over the years. We assess and discuss our tomography results using regional seismic attenuation models and estimating rock types on the basis of pressure and temperature conditions computed from thermomechanical models. Our results show significant seismic differences between the flat and normal subduction zones. As expected, the faster seismic velocities and increased seismicity within the flat slab and overriding lithosphere are generally consistent with a cooler thermal state. Our results are also consistent with dehydration of the mantle above the subducted Juan Fernandez Ridge at the eastern tip of the flat slab segment, indicating that the latter retains some fluids during subduction. However, fluids in the upper portion of the flat slab segment are not seismically detected, since we report instead fast slab seismic velocities which contradict the argument of its buoyancy being the cause of horizontal subduction. The forearc region, above the flat slab, exhibits high Vs and very low Vp/Vs ratios, uncorrelated with typical rock compositions, increased density or reduced temperature; this feature is possibly linked with the aftershock effects of the Mw7.1 1997 Punitaqui earthquake, the flat slab geometry and/or seismic anisotropy. At the surface, the seismic variations correlate with the geological terranes. The Andean crust is strongly reduced in seismic velocities along the La Ramada–Aconcagua deformation belt, suggesting structural damage. Slow seismic velocities along the Andean Moho match non-eclogitized hydrated rocks, consistent with a previous delamination event or a felsic composition, which in turn supports the extent of the Chilenia terrane at these depths. We confirm previous studies that suggest that the Cuyania terrane in the backarc region is mafic and contains an eclogitized lower crust below 50-km depth. We also hypothesize major Andean basement detachment faults (or shear zones) to extend towards the plate interface and canalize slab-derived fluids into the continental crust.