Nazca Plate Research Articles

H2 Production in Subduction Zone

Convergence zones where ophiolites outcrops, such as Oman and New Caledonia, allow the generation of natural H2 by redox phenomena. Nevertheless, the role of hydration of the mantle wedge between the overriding plate and the slab is often evoked to explain the H2 emanations. In the case of active subduction, such as that of the Nazca plate beneath the Andes, the presence of H2 emanation has been noted. In this study, we modelled the processes of (de)hydration and melting within the subduction zone. The quantification of water involved in mantle wedge hydration served as a proxy for H2 production. Multiple simulations have been conducted to investigate oceanic-continental subduction, wherein various parameters, especially the subduction angle, have been varied. It appears that in the case of flat subduction, the mantle hydration zone is large, extending up to 500 km from the trench. The extend is controlled by the velocity of the overriding plate. In the case of a steep subduction, the zone is narrower, and is located between the trench and the volcanic arc. Magma formation competes with H2 generation for the use of water expelled from the subducting plate. In the transition from a steep to a flat slab, the mantle hydration zone widens and the volcanic zone moves away from the trench. The oceanic crust may undergo melting, leading to a change in magma composition, before volcanism diminishes in intensity and then disappears.

Crustal seismicity and updated stress mapping in the Bolivian Orocline, Central Andes

The Central Andes region, with significant seismic activity, has distinct faulting mechanisms across its varied topography. Notably, the sub-Andean province predominantly features reverse faulting, whereas the high Andean plateau shows a predominance of normal and strike-slip faults. Despite the importance of this extensive orogenic plateau, Bolivia remains under-documented in seismic studies. We used data from recently installed seismic stations in Bolivia and regional stations in Brazil to determine 13 new focal mechanisms of shallow crustal earthquakes in Bolivia, some of them felt in major cities. Employing probabilistic full waveform moment tensor inversion and P-wave polarities, we mapped the stress distribution across the Central Andes. The main patterns of faulting mechanisms were: reverse faulting along the NW-SE trending sub-Andean belt north of the Orocline (north of Cochabamba) with NE-SW compression; reverse faulting along the N-S trending sub-Andean belt south of the Orocline (south of Santa Cruz) with ∼ E-W oriented compression; strike-slip faulting within the Eastern Cordillera with NE-SW P axes. The results indicate that the sub-Andean belt experiences compressional forces perpendicular to the front of the Andean plateau, arising from both the spreading stresses of the plateau and the broader plate-wide compression. On the other hand, the high plateau (Altiplano) is characterized by normal and strike-slip mechanisms, suggesting a dynamic equilibrium between local extensional gravitational stresses and regional compressional forces due to the Nazca plate convergence.

Deep‐Focus Earthquake Mechanisms at the Subducting Nazca Plate (Peru‐Brazil Border): Cold Slab Behavior in a Warm Plate

AbstractWe calculate focal mechanisms and centroid depths for deep‐focus earthquakes (DFEs) along the Peru‐Brazil border. We obtained a total of 28 focal solutions for events with magnitudes between 4.2 and 7.5 Mw and occurring between 2014 and 2022. Focal mechanisms indicate predominance of normal faulting, demonstrating a rather uniform down‐dip compression (DDC) regime within the plate. The orientations of the nodal planes suggest that earthquakes tend to occur along faults parallel to the local slab strike, although other fault types are documented. Stress orientations derived from the focal mechanisms agree with patterns expected if faulting were initiated by transformational faulting on a metastable olivine wedge (MOW) under DDC. Centroid depths range between 557 and 659 km, defining a narrow seismic zone within the lower portion of the subducting plate and an aseismic upper portion. We suggest that DFEs nucleate through transformational faulting within a narrow MOW preserved at a colder slab segment right above the lower mantle and juxtaposed to a shallower, warmer segment at around 500 km depth. This thermal complexity was possibly produced through flat subduction initiated by the subduction of the Nazca Ridge. We speculate that subduction of other aseismic ridges is possibly controlling the thermal state of the Nazca slab as a whole and, consequently, the depth distribution of DFEs along the South America subduction front.

Internal deformation of the North Andean Sliver in Ecuador-southern Colombia observed by InSAR

Summary In the Northern Andes, partitioning of oblique subduction of the Nazca plate beneath the South American continent induces a northeastward motion of the North Andean Sliver. The strain resulting from this motion is absorbed by crustal faults, which have produced magnitude 7 + earthquakes historically in the Andean Cordillera of Ecuador and southern Colombia. In order to quantify the strain in that area, we derive a high-resolution surface velocity map using InSAR time-series processing. We analyzed 6 to 8 years of Sentinel-1 data and combined different satellite line-of-sight directions to produce a reliable velocity map in the East direction. We use interpolated GNSS data to express the velocity map with respect to Stable South America and remove the long-wavelength pattern due to the post-seismic deformation following the 2016 Mw 7.8 Pedernales earthquake. The InSAR velocity map finds high E-W shortening strain rates along N-S trending structures within the Western Cordillera and the Interandean valley, with little deformation taking place east of them. This result strengthens the previous proposition of a ∼350 km long Quito-Latacunga tectonic block, forming a restraining bend in the overall right-lateral strike-slip fault system accommodating the northeastward escape motion of the North Andean Sliver. However, the high spatial resolution provided by InSAR indicates that previously proposed boundaries for this block need to be revised. In particular, InSAR results highlight high strain rate (>300 nstrain/yr) along undescribed active structures, south and west of the proposed limits for the Quito-Latacunga block, respectively in Peltetec and Ibarra regions. Interestingly, the two areas with the largest strain rates spatially correlate with the proposed areas of large historical earthquakes. Modeling of the InSAR and GNSS velocities in these areas suggests shallow coupling and high slip rates on structures which, previously, were not identified as active. We also demonstrate a slow-down of the shallow aseismic slip on the Quito fault after the Pedernales earthquake, suggesting that stress changes following large megathrust events might trigger transient slip behaviors on crustal faults. The high-resolution strain map provided by this work provides a new basis for future tectonic models in the Ecuadorian and southern Colombian Andes, and will contribute to the seismic hazard assessment in this highly populated area of the Andes.

Automatic relocation of intermediate-depth earthquakes using adaptive teleseismic arrays

SUMMARY Intermediate-depth earthquakes, accommodating intraslab deformation, typically occur within subduction zone settings at depths between 60–300 km. These events are in a unique position to inform us about the geodynamics of the subducting slab, specifically the geometry of the slab and the stress state of the host material. Improvements in the density and quality of recorded seismic data enhance our ability to determine precise locations of intermediate-depth earthquakes, in order to establish connections between event nucleation and the tectonic setting. Depth phases (near-source surface reflections, e.g. pP and sP) are crucial for the accurate determination of earthquake source depth using global seismic data. However, they suffer from poor signal-to-noise ratios in the P wave coda. This reduces the ability to systematically measure differential traveltimes to the direct P arrival, particularly for the frequent lower magnitude seismicity which highlights considerable seismogenic regions of the subducted slabs. To address this limitation, we have developed an automated approach to group globally distributed stations at teleseismic distances into ad-hoc arrays with apertures of 2.5$^\circ$, before optimizing and applying phase-weighted beamforming techniques to each array. Resultant vespagrams allow automated picking algorithms to determine differential arrival times between the depth phases and their corresponding direct P arrival. Using these differential times we can then determine the depths of earthquakes, which in turn can be used to create a catalogue of relocated events. This will allow new comparisons and insights into the governing controls on the distribution of earthquakes in subducted slabs. We demonstrate this method by relocating intermediate-depth events associated with northern Chile and the Peruvian flat slab regions of the subducting Nazca plate. The relocated Chilean catalogue contains comparable event depths to an established catalogue, calculated using a semi-automated global methodology, which serves to validate our fully automatic methodology. The new Peruvian catalogue we generate indicates three broad zones of seismicity approximately between latitudes 1–7$^\circ$S, 7–13$^\circ$S and 13–19$^\circ$S. These align with flat to steep slab dip transitions and the previously identified Pucallpa Nest. We also find a regionally deeper slab top than indicated by recent slab models, with intraslab events concentrated at points where the slab bends, suggesting a link between slab flexure and intermediate-depth earthquake nucleation.

Seismic observations of Nazca-plate crustal thicknesses providing constraints for a first-order asthenospheric-mantle potential-temperature anomalies assessment

The Nazca plate is populated with several magmatic tracks exhibiting complex morphologies responding to time-dependent positioning of spreading centers (East Pacific Rise, Cocos-Nazca). Mantle plumes beneath the Nazca plate are/were located on- or off-ridge (with ridge being spreading center), modulating the resulting hotspot-tracks’ crustal thickness accordingly. We use published seismic observational constraints of Moho depths to study asthenospheric-mantle potential-temperature anomalies (ΔTp) beneath the Nazca plate. We use a simple thermodynamics formulation for adiabatic cooling and decompression melting, and a depth-dependent, static analytical description presenting different options for hydrated-solidus-liquidus curves up to 8 GPa-lithostatic pressure as a function of bulk water content in the peridotite source. The calculated Tp-anomalies are relative to an average Nazca crust (∼6 km thick) formed at the East Pacific Rise with source hydration levels of about 0.01 wt%. As active-upwelling perturbing passive MOR environments, on-ridge plumes added mass (lower crustal intrusion, surface extrusion) to young crust creating OIB hotspot-tracks whose mean crustal thicknesses (≤18 km) are consistent with inferred ΔTp range-values of [-100, 75] °C (Iquique Ridge), [-50, 130] °C (Nazca Ridge), and [-20, 150] °C (Carnegie Ridge south), considering conservative bulk water contents of 0.005–0.08 wt%. The Juan Fernández OIB hotspot track was formed by an off-ridge active-upwelling plume impinging under a 27 Myr-old oceanic lithosphere with ∼7 km-thick pre-existing crust. Currently, the track presents 4-5 km-high isolated volcanic edifices and about 1 km crustal root, totaling about 12–13 km aggregate melt thickness, and suggesting ΔTp range-values of [-20, 160] °C.

Comparison between rupture parameters of intermediate and deep earthquakes at the Peru–Brazil–Bolivia border and northern Chile

We determined the main parameters of the source rupture process of intermediate- and deep-depth earthquakes occurring in the Peru–Brazil–Bolivia border region and northern Chile. The parameters of depth, fault-plane orientation, scalar seismic moment, slip distribution, and radiated seismic energy are obtained from seismograms. We selected 15 intermediate-depth earthquakes (100 km < h < 300 km) and 10 very deep earthquakes (h > 500 km) with magnitudes MW ≥ 6.0. For most events, the slip distribution over the rupture plane shows a single asperity, and the source time function presents a simple pulse. There are differences between intermediate-depth and deep earthquakes. The rupture areas, maximum slip and source time function (STF) duration are larger for intermediate-depth events than for deep events. Additionally, the STF’s show a sharper increase for deep earthquakes. The scaled radiated seismic energy shows larger values for deep depth events. The stress regime pattern derived from the obtained focal mechanism agrees with the geometry of the subduction of the Nazca plate. At intermediate depths, in the northern area up to 12°S, the stress pattern corresponds to a horizontal extension, while in the southern area, the tension axes dip at an angle of 30°. At deep depths, the stress regime corresponds to vertical compression in the north and dips of approximately 45° in the south.

The thermal structure of the Colombian lithosphere: A regional and basin-scale analysis

It is well-established that the thermal state of the lithosphere strongly influences various regional and local geological processes, including crustal deformation, hydrocarbon maturation, hydrogen generation, and geothermal phenomena. Moreover, the thermal structure exhibits high sensitivity to tectonic features, a property of particular significance in Colombia, where three main tectonic plates converge, and lithospheric tearing has been documented. In this contribution, we focus on elucidating the impact of plate architecture on the thermal field in central-eastern Colombia at both shallow and deep levels. To accomplish this, we constructed a series of two-dimensional profiles and derived a numerical solution of the heat equation using the conservative finite difference method. As constraints, we incorporate an inferred distribution of rocks in the deep crust and upper mantle based on global and local lithospheric thickness models. Material parameters for the various rocks, both exposed and inferred, were obtained from the literature. Additionally, we used superficial heat flow estimates and apparent geothermal gradients compiled by the Colombian Geological Survey.Our results suggest a significant influence of the lithospheric Caldas tear on the thermal state of Colombia, with the breaking off occurring in the Nazca plate under the Eastern Cordillera range around 5°N. The modeled asthenospheric heat flow remains approximately 25 mWm−2, except in the northern Eastern Cordillera range, where the background heat flow increases rapidly to 40 mWm−2. Consequently, our model predicts partial melting in the lower crust to the north and a thermally unstable lower crust to the south of the Caldas tear. The material parameters that best fit the surface data suggest the presence of a basement moderately enriched in radioactive elements in the Eastern Llanos basin. After accounting for compaction, we also confirm a strong influence of the tectonic setting on the thermal state of sedimentary basins.

Joint flexural-density modeling of the Taltal, Copiapó, and Iquique hotspot ridges and the surrounding oceanic plate, offshore Chile

Abstract Based on gravity and bathymetric data and using a novel two-dimensional joint flexural-density modeling approach, this work studies the physical properties of the oceanic Nazca plate around the Taltal, Copiapó, and Iquique hotspot ridges offshore northern Chile. The area is located westward of the Chilean Trench where the Taltal and Copiapó Ridges collide with the continental margin. The results show that the variability in density structure at different scales is a key factor in explaining the observed gravity signal, playing an important role in the lithospheric flexure and hence the elastic properties of the Nazca plate in this setting. The results can be interpreted as evidence of spatial and temporal heterogeneities in the plate-weakening process at the hotspots, magmatic underplating, and crustal and upper mantle fracturing and/or hydration. These processes might be relevant for the ascent of magma pathways of later (secondary) volcanism and influence the mechanical segmentation of the oceanic plate. The latter is critical in explaining the active seismogenic contact between the oceanic Nazca and overriding South America plates.

Plume‐Driven Subduction Termination in 3‐D Mantle Convection Models

AbstractThe effect of mantle plumes is secondary to that of subducting slabs for modern plate tectonics when considering plate driving forces. However, the impact of plumes on tectonics and planetary surface evolution may nonetheless have been significant. We use numerical mantle convection models in a 3‐D spherical chunk geometry with damage rheology to study some of the dynamics of plume‐slab interactions. Substantiating our earlier 2‐D results, we observe a range of interaction scenarios, and that the plume‐driven subduction terminations we had identified earlier persist in more realistic convective flow. We analyze the dynamics of plume affected subduction, including in terms of their geometry, frequency, and the overall effect of plumes on surface dynamics as a function of the fraction of internal to bottom heating. Some versions of such plume‐slab interplay may be relevant for geologic events, for example, for the inferred ∼183 Ma Karoo large igneous province formation and associated slab disruption. More recent examples may include the impingement of the Afar plume underneath Africa leading to disruption of the Hellenic slab, and the current complex structure imaged for the subduction of the Nazca plate under South America. Our results imply that plumes may play a significant role not just in kick‐starting plate tectonics, but also in major modifications of slab‐driven plate motions, including for the present‐day mantle.

The Volcanism of the Meseta del Lago Buenos Aires, Patagonia: the Transition from Subduction to Slab Window

Abstract The Meseta del Lago Buenos Aires (MLBA) in southern Patagonia, a volcanic plateau formed from ~12 Ma to present, provides an opportunity to investigate the temporal evolution in volcanism as this region transitions from the subduction of the Nazca plate to the formation of the slab window produced by the collision of the Chile Ridge and the Andean subduction zone. Here, we report new major, minor, and trace element contents, as well as Sr, Nd, Pb, and Hf isotopes of the MLBA lavas. Three distinct geochemical endmembers can be distinguished in the MLBA basalts: a subduction-influenced endmember, a transitional component similar to the South Atlantic enriched mid-ocean ridge basalts, and an enriched component akin to the EM1 mantle composition. Lavas older than ~1.5 Ma define a compositional continuum between the subduction-influenced and transitional endmembers; this trend is also present in many other southern Patagonian plateaus regardless of their distance to the trench, eruption age, and the composition of the continental blocks where they are located. In contrast, MLBA basalts younger than ~1.5 Ma uniquely define a transition into the EM1 mantle component at the time when this region was affected by the slab window. The estimated pressures and temperatures of mantle-melt equilibration for the MLBA basalts indicates an increase in both parameters after the formation of the slab window that roughly correlate with the changes in lava composition. The basalts’ composition from all southern Patagonia plateaus points to the presence of the South Atlantic mid-ocean ridge basalt mantle influenced by the Discovery, Shona, and Bouvet hotspots rather than the sub-slab mantle, as represented by the Chile Ridge basalts. This observation challenges the hypothesis that the sub-slab mantle within the slab window has had an important role in the composition of the erupted lavas. Instead, it suggests the presence of a South Atlantic mantle beneath southern Patagonia either within the mantle wedge, consistent with a long-lasting South Atlantic convection cell beneath South America, or in the subcontinental lithospheric mantle metasomatized before or just after the opening of the South Atlantic basin, as demonstrated by the composition of southern Patagonia mantle xenoliths. Although it is difficult to precisely distinguish the contributions of the asthenosphere from that of the metasomatized subcontinental lithospheric mantle beneath this region, our work suggests significant contributions from the latter in the composition of the MLBA lavas.

Stratigraphy and structure of Chachahuén volcanic complex, southern Mendoza province, Argentina

The Chachahuén Volcanic Complex is an eroded Upper Miocene arc to back-arc volcanic system that was emplaced in the north-eastern edge of the Neuquén Basin, southern Mendoza province, Argentina. It was formed as the result of shallowing followed by steepening of the Nazca plate during Miocene-Pliocene, and is characterized by a volcanism which ranges from trachydacites to basalts. Thanks to erosion, the stratigraphy and the shallow plumbing systems of the volcano are well-exposed. A new evolution for the volcano stratigraphy is proposed, based on new field observations, satellite imagery, and incorporating the radiometric and petrological data of previous works. The stratigraphy is constrained around one main explosive event, recorded by a thick dacitic pyroclastic density current (PDC) deposit (the Corrales Ignimbrite) (>150m). The main units forming the Chachahuén Volcanic Complex are (1) the pre-dacite PDC deposits (trachydacitic to rhyolitic) Vizcachas Formation, (2) the thick Corrales Ignimbrite, (3) a post-dacite trachyandesitic PDC deposits (dominated by block and ash (BAF) deposits) and thick lava flows, and (4) thin mafic basaltic lava flows. The collapse of a large elliptical caldera occurred during the Trachyandesite phase. The main magma transport channels and feeders are sub-vertical dykes, which exhibit a radial distribution centred on the main caldera depression. Moreover, cryptodomes mainly preserved outside the rim of the caldera accommodated the transport and emplacement of the more silicic magma during the Trachyandesite phase.

The relationship between GRACE gravity and the seismic b-value: a case study of the Northern Chile Triple Junction (25° S–40° S)

SUMMARY The northern Chile Triple Junction (CTJ) is characterized by the ongoing subduction of the Nazca plate beneath the South American plate. The geological structures within the subduction zone undergo complex changes, resulting in significant tectonic activities and intense seismicity along the western margin of South America. Based on the Gravity Recovery and Climate Experiment (GRACE) data and earthquake catalogues, this study selects the northern CTJ area (25° S–40° S, 75° W–65° W) as the research object, adopts the mathematical methods of independent component analysis (ICA) and principal component analysis (PCA) to separate the earthquake-related signals within the GRACE data, and fits the changes of seismic b-values through the frequency–magnitude relationship. The characteristics of gravity changes before and after seismic events, the seismic activity parameter b-values, and the relationship between the gravity signals and b-values are discussed. The results show that mathematical methods can effectively extract seismic-related gravity components from the GRACE data. ICA, compared to PCA, provides better results in capturing the temporal variations associated with b-value time-series, which exhibit good consistency in long-term trend changes. The average change of b-values in the study area is 0.66 ± 0.003, fluctuating over time. Generally, prior to larger seismic events, b-values tend to decrease. Along the western margin of South America, b-values are low; this aligns with the active tectonic activities between subducting plates. Additionally, a certain correlation between b-values and gravity changes is observed, but due to the influence of tectonic activities, the correspondence between b-values and gravity anomalies may not be consistent across different areas. The b-value is highly consistent with the strain rate model. Low b-values correspond to high strain rates along the western edge of South America, which is in line with the tectonic characteristics of frequent seismic activity in this area. A gradual concentration of gravity anomalies before major earthquakes is observed, accompanied by the gradual accumulation of smaller seismic events. Meanwhile, several months before the two major earthquakes, the spatial distribution of gravity appears to be similar to the coseismic signals, but the nature of its generation remains to be explored. These methods and results not only add to the applications of GRACE in seismic studies but also raise questions for further exploration.

Thermal modeling of subduction zones with prescribed and evolving 2D and 3D slab geometries

The determination of the temperature in and above the slab in subduction zones, using models where the top of the slab is precisely known, is important to test hypotheses regarding the causes of arc volcanism and intermediate-depth seismicity. While 2D and 3D models can predict the thermal structure with high precision for fixed slab geometries, a number of regions are characterized by relatively large geometrical changes over time. Examples include the flat slab segments in South America that evolved from more steeply dipping geometries to the present day flat slab geometry. We devise, implement, and test a numerical approach to model the thermal evolution of a subduction zone with prescribed changes in slab geometry over time. Our numerical model approximates the subduction zone geometry by employing time dependent deformation of a Bézier spline that is used as the slab interface in a finite element discretization of the Stokes and heat equations. We implement the numerical model using the FEniCS open source finite element suite and describe the means by which we compute approximations of the subduction zone velocity, temperature, and pressure fields. We compute and compare the 3D time evolving numerical model with its 2D analogy at cross-sections for slabs that evolve to the present-day structure of a flat segment of the subducting Nazca plate.

Implications of Flat‐Slab Subduction on Hydration, Slab Seismicity, and Arc Volcanism in the Pampean Region of Chile and Argentina

AbstractThe Pampean flat slab in central Chile and Argentina is characterized by the inland migration and subsequent cessation of arc volcanism since the mid‐Miocene. Slab flattening also affects the distribution and number of intermediate‐depth earthquakes and the evolution of the overlying continental thermal structure. In this study, we combine thermal‐mechanical models with petrological models to examine temporal changes in pressure, temperature, and composition during flat‐slab subduction and estimate water carrying capacity, predicted melt distributions and predicted changes in melt composition. Model results indicate that the present‐day flattened Nazca plate carries water to ∼700 km inland from the trench and could cause flux melting if the material above the slab remains fertile. Observed slab seismicity matches areas where hydrated materials have ∼&gt;3 wt% H2O in the oceanic crust and mantle lithosphere. Seismicity increases as slab water carrying capacity decreases (slab dehydration). As P‐T conditions and compositions of the rock trapped above the slab change during slab flattening, flux melting switches from a peridotite‐dominated early phase to a combined mid‐ocean ridge basalt/eclogite and peridotite melting at ∼8 Ma. The results provide broad consistency with known earthquake distributions, seismic velocities, and observed temporal and spatial changes in volcanic patterns above the Pampean flat slab and point toward the role of melt depletion in the decrease and ultimate cessation of arc volcanism in this region.

Comparative Analysis of Structural Reinforcement with Viscoelastic Energy Dissipators, Friction and Metal Creep in Tall Buildings

Seismic risk is a challenging problem in tall buildings due to the possibility of loss of human life and economic caused by seismic events. Peru is located at the interaction of the South American plate and the Nazca plate, which is why various seismic events of moderate to large magnitude occur. Today there are many ways to solve these problems and it is a very challenging case to reinforce tall buildings. In addition, technological advances in software facilitate and help through programmed models in tall buildings that analyze their structure characteristics such as drift, shear and others. This article proposes a comparative analysis of three types of dissipators: viscous fluid, friction, and metal creep through a Time-History analysis in a 15-story high-rise building located in Peru. The proposed methodology considers three stages: (i) definition of the characteristics and properties of the structure in accordance with Peruvian Standard E.030, in addition three accelerograms are used for the dynamic time-history analysis and maximum displacements and drifts are determined by ETABS software. (ii) calculate the design drift of the tall building and the properties of the viscous fluid, friction, and creep dissipator. In addition, calculations are made for the design parameters of each dissipator, and it is modeled as required for the case study. (iii) the new drifts and the damping values that the building presents for each dissipator are analyzed. According to the results obtained, the dissipator with the best results is of the flow type, since it has better performance in drifts and manages to produce an average damping of 96.87% for tall buildings. While the viscous dissipators obtain a 57.85% damping and the friction ones are estimated at 81.57%.

Tracing hotspot traces in the Andes

Two segments of subduction of the Nazca plate beneath the South American plate occur at low angles based on seismic hypocenter locations, approaching nearly horizontal below ~100 km in depth. In contrast with most of the rest of the subduction zone, the two segments, beneath central Chile, and central and northern Peru, lack active volcanoes along the crest of the Andes and have more subdued topography to the east of the Andean crest. Each low-angle subduction segment occurs to the east of the intersection of inferred mantle hotspot traces on the Nazca plate with the Peru-Chile Trench: the Nazca ridge (at the southern part of the Peruvian segment), and the Juan Fernández island-seamount chain (offshore the Chilean segment). A third inferred trace, the Galápagos-Carnegie ridge, may be correlated with a zone on incipient low-angle subduction beneath Colombia. &#x0D; The importance of such hotspot traces in contributing to low-angle subduction beneath the Andes is strengthened by updated South American-Nazca plate reconstructions, including three oceanic hotspot traces, in comparison with a new isotopic date compilation of igneous rocks from the mountain range. The Juan Fernández hotspot trace, reconstructed from Pacific-hotspot models to the Nazca-Farallon plate, encountered the subduction zone offshore southern Peru ~65 Ma, broadening arc volcanism to the east; the trace-trench intersection migrated gradually and then rapidly southward, widening the arc east to Bolivia and northern Argentina; it then stabilized about 13 Ma offshore central Chile, producing the contemporary low-angle Pampean segment. The Juan Fernández hotspot may also have been responsible for formation of the Manihiki Plateau on the Pacific plate much earlier, ~125 Ma. The Easter-Nazca hotspot trace intersected the subduction zone beneath Colombia before ~50 Ma and migrated southward beneath Ecuador beginning ~15 Ma, with progressive low-angle subduction implied by migrating volcanic cessation along the Andean crest to southern Peru. The Galápagos-Carnegie hotspot trace only recently encountered the subduction zone, apparently inducing a new low-angle segment and cessation of magmatism in Colombia. The reconstructions and magmatic history provided here strongly support a previously proposed genetic relationship of hotspot traces and low-angle subduction. Additionally, the reconstructions suggest remnants of older subducted traces in the asthenosphere may have sourced post-rift magmatism in eastern Brazil and Paraguay, which cannot be explained otherwise by simple hotspot mechanisms.

Plio-Pleistocene rear-arc volcanism in the Southern Volcanic Zone: Eruptive styles of the Varvarco Volcanic Field

The Varvarco Volcanic Field (VVF) is located in the southern part of the Las Loicas Trough, as part of the Late Pliocene-Early Pleistocene rear-arc volcanic belt in the Transitional Southern Volcanic Zone (34.5–37°S). Its volcanic products show an elliptical distribution, elongated parallel to NW-SE main structures that regionally controlled the Las Loicas Trough. A detailed field and petrographic study was carried out to identify main lithofacies and establish its eruptive styles.The VVF magmatic evolution is initially characterized by a voluminous explosive stage represented in the area by dense and dilute pyroclastic density currents (PDC) deposits (massive lapilli tuffs, cross-stratified lapilli tuffs and diffuse-stratified tuffs). Afterwards, it evolved into an effusive stage represented by basaltic lava flows (coherent basalts), associated with Hawaiian to Strombolian eruptive style, which constitutes most of the VVF volume. The final stage of the VVF history was linked to a stratovolcano-type activity where both effusive (coherent basalts and andesites, and rhyolitic coulees) and explosive lithofacies (such as massive lithic breccias and massive lapilli tuffs) are described. Within this stage, the uppermost effusive levels were intruded by dacitic and rhyolitic domes and basaltic dykes. Available ages allow to conclude that the VVF emplacement was developed during Plio-Pleistocene times, linked to the re-steepening of the Nazca plate, after the Late Miocene Payenia shallow subduction regime.

Are Volatiles From Subducted Ridges on the Pampean Flat Slab Fracking the Crust? Evidence From an Enhanced Seismicity Catalog

AbstractSeamounts and ridges are often invoked to explain subduction‐related phenomena such as flat slab generation, but the extent of their involvement remains controversial. An analysis of seismicity in the region of the Pampean flat slab through an application of an automated catalog generation algorithm resulted in 35,924 well constrained local earthquake hypocenters and a total of 12,172 focal mechanisms. Several new features related to the subduction of the Juan Fernandez Ridge (JFR) were discovered, including (a) a series of parallel lineaments of seismicity in the subducted Nazca plate separated by about 50 km and trending about 20°, and (b) a strong spatial correlation between these deeper (&gt;80 km depth) regions of intense seismicity and concentrations of activity in the crust almost directly above it. Focal mechanisms of the deeper events are almost exclusively (∼81%) normal, while those in the crust are predominantly (∼70%) reverse. The deeper lineaments mirror the orientation and spacing of several seamount chains seen on the Nazca plate, suggesting that these patterns are caused by the same types of features at depth. This would imply that relatively minor features persist as slab anomalies long after they are subducted. The correlation of the deeper seismicity that defines these features with seismicity in the mid to lower crust suggests a genetic relation between the two. We postulate that volatiles from the subducted ridges percolate into the South American crust and induce seismicity essentially by fracking it.

Cryptic trans-lithospheric fault systems at the western margin of South America: implications for the formation and localization of gold-rich deposit superclusters

We present a review of frontier research advances in the investigation of cryptic structures that transect the South American Andes at oblique strike directions. The intersections between these cryptic structures and the superimposed Andean belt correlate with the spatial distribution of gold-rich mineral deposit clusters. The deposit clusters can be described as superclusters, as they comprise various gold deposit types that formed at multiple times throughout the Phanerozoic, impinging repeatedly on the structural intersections. However, the cryptic inherited fault structures are difficult to detect, because their deeper-seated roots are often overlain by younger supracrustal successions, and/or their exposed surface manifestations are structurally obscured by subsequent tectonic-magmatic activity. Thus, it also remains a challenge to constrain the nature and timing of formation, and the respective subsequent evolutionary path, of these proposed pre-Andean structures. Based on various case studies, we demonstrate that the localization of identified Phanerozoic gold deposit superclusters along the western South American margin is fundamentally controlled by structural inheritance often dating back to at least the Mesoproterozoic. Integration of multi-approach observations and datasets allows insights into a larger-scale tectonic history that showcases the successive inheritance of major structures originating from the Amazonian Craton, over the Paleozoic Gondwana margin, into the Cenozoic magmatic belts of the Andes, and even into recent fractures within the subducting oceanic Nazca plate, recording &amp;gt;1.2-billion-years of progressive structural inheritance and growth at one of the longest-lived tectonic margins in Earth history. In contrast to previous models of the spatial distribution of gold deposits, based on statistical approaches and spatial periodicity in self-organized systems focusing on single subduction and/or accretion episodes and belts, we propose that the structural inheritance and intersections are key to the localization of gold deposits in the Andes. In combination with bulk-geochemical data from magmatic rocks, we suggest that inherited structures maintained a trans-lithospheric connectivity to pre-fertilized gold enriched upper mantle reservoirs, which were tapped during multiple tectono-magmatic reactivation episodes.