scholarly journals Subduction zone fluids and arc magmas conducted by lithospheric deformed regions beneath the central Andes

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
E. Contreras-Reyes ◽  
D. Díaz ◽  
J. P. Bello-González ◽  
K. Slezak ◽  
B. Potin ◽  
...  

AbstractDehydration of the oceanic subducting slab promotes the formation of magmatic arcs, intra-slab intermediate-depth seismicity, and hydration of the overlying mantle wedge. However, the complex permeability structure of the overriding plate controls the magma and fluid migration and their accumulation at shallower depths. In this regard, mapping the inner structure of the overriding crust and mantle is crucial to understand the magmatic and hydrological processes in subduction zones. We integrate 3-D P-wave, $$V_p/V_s$$ V p / V s , and electrical resistivity tomographic models of the northern Chilean subduction zone to map the magmatic and fluids derived from the subducting oceanic Nazca plate. Results show a continental crust relatively thick (50–65 km) characterized by a lower zone of high $$V_p$$ V p values (7.2–7.6 km/s), which is interpreted as the presence of plutonic rocks. The mantle lithospheric wedge is weakly hydrated ($$V_p/V_s$$ V p / V s = 1.75–1.8) while the forearc continental crust is traversed by regions of reduced electrical resistivity values ($$< 10^2$$ < 10 2 $$\Omega m$$ Ω m ) interpreted as zones of relatively high permeability/fracturing and fluid content. These regions spatially correlate with upper plate trans-lithospheric deformation zones. Ascending melts accumulate preferentially in the back-arc, whereas hydrothermal systems form trenchward of the volcanic arc. The results highlight the complex permeability structure of the upper South American plate.

2012 ◽  
Vol 4 (1) ◽  
pp. 745-781 ◽  
Author(s):  
C. J. Warren

Abstract. The exhumation of high and ultra-high pressure rocks is ubiquitous in Phanerozoic orogens created during continental collisions, and is common in many ocean-ocean and ocean-continent subduction zone environments. Three different tectonic environments have previously been reported, which exhume deeply buried material by different mechanisms and at different rates. However it is becoming increasingly clear that no single mechanism dominates in any particular tectonic environment, and the mechanism may change in time and space within the same subduction zone. In order for buoyant continental crust to subduct, it must remain attached to a stronger and denser substrate, but in order to exhume, it must detach (and therefore at least locally weaken) and be initially buoyant. Denser oceanic crust subducts more readily than more buoyant continental crust but exhumation must be assisted by entrainment within more buoyant and weak material such as serpentinite or driven by the exhumation of structurally lower continental crustal material. Weakening mechanisms responsible for the detachment of crust at depth include strain, hydration, melting, grain size reduction and the development of foliation. These may act locally or may act on the bulk of the subducted material. Metamorphic reactions, metastability and the composition of the subducted crust all affect buoyancy and overall strength. Subduction zones change in style both in time and space, and exhumation mechanisms change to reflect the tectonic style and overall force regime within the subduction zone. Exhumation events may be transient and occur only once in a particular subduction zone or orogen, or may be more continuous or occur multiple times.


2015 ◽  
Vol 112 (5) ◽  
pp. 1287-1291 ◽  
Author(s):  
David Bercovici ◽  
Gerald Schubert ◽  
Yanick Ricard

A simple model for necking and detachment of subducting slabs is developed to include the coupling between grain-sensitive rheology and grain-size evolution with damage. Necking is triggered by thickened buoyant crust entrained into a subduction zone, in which case grain damage accelerates necking and allows for relatively rapid slab detachment, i.e., within 1 My, depending on the size of the crustal plug. Thick continental crustal plugs can cause rapid necking while smaller plugs characteristic of ocean plateaux cause slower necking; oceanic lithosphere with normal or slightly thickened crust subducts without necking. The model potentially explains how large plateaux or continental crust drawn into subduction zones can cause slab loss and rapid changes in plate motion and/or induce abrupt continental rebound.


2020 ◽  
Vol 224 (3) ◽  
pp. 1553-1571
Author(s):  
Emily E Rodríguez ◽  
Daniel Evan Portner ◽  
Susan L Beck ◽  
Marcelo P Rocha ◽  
Marcelo B Bianchi ◽  
...  

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 &gt;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 &gt;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.


2021 ◽  
Author(s):  
Jungjin Lee ◽  
Mainak Mookherjee ◽  
Taehwan Kim ◽  
Haemyeong Jung ◽  
Reiner Klemd

&lt;p&gt;Subduction zones are often characterized by the presence of strong trench-parallel seismic anisotropy and large delay times. Hydrous minerals, owing to their large elastic anisotropy and strong lattice preferred orientations (LPOs) are often invoked to explain these observations. However, the elasticity and LPO of chloritoid, which is one such hydrous phases relevant in subduction zone settings, is poorly understood. In this study, we measured the LPO of polycrystalline chloritoid in natural rock samples and obtained the LPO-induced seismic anisotropy and evaluated the thermodynamic stability field of chloritoid in subduction zones. The LPO of chloritoid aggregates displayed a strong alignment of the [001] axes subnormal to the rock foliation, with a girdle distribution of the [100] axes and the (010) poles subparallel to the foliation. New elasticity data of single-crystal chloritoid showed a strong elastic anisotropy of chloritoid with 47% for S-waves (V&lt;sub&gt;S&lt;/sub&gt;) and 22% for P-waves (V&lt;sub&gt;P&lt;/sub&gt;), respectively. The combination of the LPO and the elastic anisotropy of the chloritoid aggregates produced a strong S-wave anisotropy of AV&lt;sub&gt;S&lt;/sub&gt; = 18% and a P-wave anisotropy of AV&lt;sub&gt;P&lt;/sub&gt; = 10%. Our results indicate that the strong LPO of chloritoid along the hydrated slab-mantle interface and in subducting slabs can influence trench-parallel seismic anisotropy in subduction zones with &amp;#8220;cold&amp;#8221; geotherms.&lt;/p&gt;


2021 ◽  
Vol 228 (2) ◽  
pp. 729-743
Author(s):  
Jiaqi Li ◽  
Min Chen ◽  
Jieyuan Ning ◽  
Tiezhao Bao ◽  
Ross Maguire ◽  
...  

SUMMARY The detailed structure near the 410-km discontinuity provides key constraints of the dynamic interactions between the upper mantle and the lower mantle through the mantle transition zone (MTZ) via mass and heat exchange. Meanwhile, the temperature of the subducting slab, which can be derived from its fast wave speed perturbation, is critical for understanding the mantle dynamics in subduction zones where the slab enters the MTZ. Multipathing, i.e. triplicated, body waves that bottom near the MTZ carry rich information of the 410-km discontinuity structure and can be used to constrain the discontinuity depth and radial variations of wave speeds across it. In this study, we systematically analysed the trade-off between model parameters in triplication studies using synthetic examples. Specifically, we illustrated the necessity of using array-normalized amplitude. Two 1-D depth profiles of the wave speed below the Tatar Strait of Russia in the Kuril subduction zone are obtained. We have observed triplications due to both the 410-km discontinuity and the slab upper surface. And, seismic structures for these two interfaces are simultaneously inverted. Our derived 410-km discontinuity depths for the northern and southern regions are at 420$\pm $15 and 425$\pm $15 km, respectively, with no observable uplift. The slab upper surface is inverted to be located about 50–70 km below the 410-km discontinuity. This location is between the depths of the 1 and 2 per cent P-wave speed perturbation contours of a regional 3-D full-waveform inversion (FWI) model, but we found twice the wave speed perturbation amplitude. A wave speed increase of 3.9–4.6 per cent within the slab, compared to 2.0–2.4 per cent from the 3-D FWI model, is necessary to fit the waveforms with the shortest period of 2 s, indicating that high-frequency waves are required to accurately resolve the detailed structures near the MTZ.


2021 ◽  
Vol 73 (1) ◽  
Author(s):  
Maria Jose Segovia ◽  
Daniel Diaz ◽  
Katarzyna Slezak ◽  
Felipe Zuñiga

AbstractTo analyze the process of subduction of the Nazca and South American plates in the area of the Southern Andes, and its relationship with the tectonic and volcanic regime of the place, magnetotelluric measurements were made through a transversal profile of the Chilean continental margin. The data-processing stage included the analysis of dimensional parameters, which as first results showed a three-dimensional environment for periods less than 1 s and two-dimensional for periods greater than 10 s. In addition, through the geomagnetic transfer function (tipper), the presence of structural electrical anisotropy was identified in the data. After the dimensional analysis, a deep electrical resistivity image was obtained by inverting a 2D and a 3D model. Surface conductive anomalies were obtained beneath the central depression related to the early dehydration of the slab and the serpentinization process of the mantle that coincides in location with a discontinuity in the electrical resistivity of a regional body that we identified as the Nazca plate. A shallow conductive body was located around the Calbuco volcano and was correlated with a magmatic chamber or reservoir which in turn appears to be connected to the Liquiñe Ofqui fault system and the Andean Transverse Fault system. In addition to the serpentinization process, when the oceanic crust reaches a depth of 80–100 km, the ascending fluids produced by the dehydration and phase changes of the minerals present in the oceanic plate produce basaltic melts in the wedge of the subcontinental mantle that give rise to an eclogitization process and this explains a large conductivity anomaly present beneath the main mountain range.


2019 ◽  
Vol 9 (1) ◽  
pp. 13-21
Author(s):  
Rezene Mahatsente

Abstract The Central Andes experienced major earthquake (Mw =8.2) in April 2014 in a region where the giant 1877 earthquake (Mw=8.8) occurred. The 2014 Iquique earthquake did not break the entire seismic gap zones as previously predicted. Geodetic and seismological observations indicate a highly coupled plate interface. To assess the locking mechanism of plate interfaces beneath Central Andes, a 2.5-D gravity model of the crust and upper mantle structure of the central segment of the subduction zone was developed based on terrestrial and satellite gravity data from the LAGEOS, GRACE and GOCE satellite missions. The densities and major structures of the gravity model are constrained by velocity models from receiver function and seismic tomography. The gravity model defined details of crustal and slab structure necessary to understand the cause of megathrust asperity generation. The densities of the upper and lower crust in the fore-arc (2970 – 3000 kg m−3) are much higher than the average density of continental crust. The high density bodies are interpreted as plutonic or ophiolitic structures emplaced onto continental crust. The plutonic or ophiolitic structures may be exerting pressure on the Nazca slab and lock the plate interfaces beneath the Central Andes subduction zone. Thus, normal pressure exerted by high density fore-arc structures and buoyancy force may control plate coupling in the Central Andes. However, this interpretation does not exclude other possible factors controlling plate coupling in the Central Andes. Seafloor roughness and variations in pore-fluid pressure in sediments along subduction channel can affect plate coupling and asperity generation.


2021 ◽  
Author(s):  
Thomas Alcock ◽  
Sergio Vinciguerra ◽  
Phillip Benson ◽  
Federico Vagnon

&lt;p&gt;Stromboli volcano has experienced four sector collapses over the past 13 thousand years, resulting in the formation of the Sciara del Fuoco (SDF) horseshoe-shaped depression and an inferred NE / SW striking rift zone across the SDF and the western sector of the island. These events have resulted in the formation of steep depressions on the slopes on the volcano where episodes of instability are continuously being observed and recorded. This study aims to quantify the fracture density inside and outside the rift zone to identify potential damaged zones that could reduce the edifice strength and promote fracturing. In order to do so we have carried out a multiscale analysis, by integrating satellite observations, field work and seismic and electrical resistivity analyses on cm scales blocks belonging to 11 lava units from the main volcanic cycles that have built the volcano edifice, ie. Paleostromboli, Nestromboli and Vancori. 0.5 m resolution Pleiades satellite data has been first used to highlight 23635 distinct linear features across the island. Fracture density has been calculated using Fracpaq based on the Mauldon et al (2001) method to determine the average fracture density of a given area on the basis of the average length of drawn segments within a predetermined circular area. 41.8 % of total fracture density is found around intrusions and fissures, with the summit area and the slopes of SDF having the highest average fracture density of 5.279&amp;#160; . Density, porosity, P- wave velocity in dry and wet conditions and electrical resistivity (in wet conditions) were measured&amp;#160; via an ultrasonic pulse generator and acquisition system (Pundit) and an on purpose built measuring quadrupole on cm scale blocks of lavas collected from both within and outside the proposed rift zone to assess the physical state and the crack damage of the different lava units.&amp;#160; Preliminary results show that P-wave velocity between ~ 2.25 km/s &lt; Vp &lt; 5km/s decreases with porosity while there is high variability electrical resistivity with 21.7 &lt; &amp;#961; &lt; 590 Ohm * m. This is presumably due to the lavas texture and the variable content of bubble/vesicles porosity and crack damage, that is reflected by an effective overall porosity between 0 and 9 %. Higher porosity is generally mirrored by lower p-wave velocity values. Neostromboli blocks show the most variability in both P-wave velocity and electrical resistivity. Further work will assess crack density throughout optical analyses and systematically investigate the UCS and elastic moduli. This integrated approach is expected to provide a multiscale fracture density and allow to develop further laboratory testing on how slip surfaces can evolve to a flank collapse at Stromboli.&lt;/p&gt;


There are well established differences in the chemical and isotopic characteristics of the calc-alkaline basalt—andesite-dacite-rhyolite association of the northern (n.v.z.), central (c.v.z.) and southern volcanic zones (s.v.z.) of the South American Andes. Volcanic rocks of the alkaline basalt-trachyte association occur within and to the east of these active volcanic zones. The chemical and isotopic characteristics of the n.v.z. basaltic andesites and andesites and the s.v.z. basalts, basaltic andesites and andesites are consistent with derivation by fractional crystallization of basaltic parent magmas formed by partial melting of the asthenospheric mantle wedge containing components from subducted oceanic lithosphere. Conversely, the alkaline lavas are derived from basaltic parent magmas formed from mantle of ‘within-plate’ character. Recent basaltic andesites from the Cerro Galan volcanic centre to the SE of the c.v.z. are derived from mantle containing both subduction zone and within-plate components, and have experienced assimilation and fractional crystallization (a.f.c.) during uprise through the continental crust. The c.v.z. basaltic andesites are derived from mantle containing subduction-zone components, probably accompanied by a.f.c. within the continental crust. Some c.v.z. lavas and pyroclastic rocks show petrological and geochemical evidence for magma mixing. The petrogenesis of the c.v.z. lavas is therefore a complex process in which magmas derived from heterogeneous mantle experience assimilation, fractional crystallization, and magma mixing during uprise through the continental crust.


2019 ◽  
Vol 9 (2) ◽  
pp. 152
Author(s):  
Rahmat Setyo Yuliatmoko ◽  
Telly Kurniawan

The amount of stress released by an earthquake can be calculated with a stress drop, the stress ratio before and after an earthquake where the stress accumulated in a fault or a subduction zone is immediately released during an earthquake. The purpose of this research is to calculate the amount of stress drop in faults and subduction in Maluku and Halmahera and their variations and relate them to the geological conditions in the area so that the tectonic characteristics in the area can be identified. This research employed mathematical analysis and the Nelder Mead Simplex nonlinear inversion methods. The results show that Maluku and Halmahera are the area with complex tectonic conditions and large earthquake impacts. The Maluku sea earthquake generated a stress drop of 0.81 MPa with a reverse fault mechanism in the zone of subduction, while for the Halmahera earthquake the stress drop value was 52.72 MPa, a typical strike-slip mechanism in the fault zone. It can be concluded that there is a difference in the stress drop between the subduction and fault zones; the stress drop in the fault was greater than that in the subduction zone due to different rock structure and faulting mechanisms as well as differences in the move slip rate that plays a role in the process of holding out the stress on a rock. This information is very important to know the amount of pressure released from the earthquake which has a very large impact as part of disaster mitigation measures.


Sign in / Sign up

Export Citation Format

Share Document