Search

Recordings of earthquakes by a temporary deployment of 84 short period seismometers in northern Chile were used to derive regional 3D seismic velocity models for the Taltal segment. We used the Regressive ESTimator (REST) package for event detection and automatic onset estimation of P‐ and S‐wave arrival times to create an earthquake catalog with 23,985 hypocenters. We followed standard acceptability criteria (i.e., azimuthal gap and residual cutoff) to create a high‐quality data set and inverted for 3D Vp, Vs and Vp/Vs models using local earthquake tomography. Plots of hypocenters from the catalog, comprising 16,349 earthquakes, reveal active structures in the upper crust, dip changes along the slab and fracturing within the oceanic crust. Moreover, the wavespeed models illuminate anomalies in both the Nazca and South American plates that correlate with the observed seismicity distribution, including variations from low (1.75) to high (>1.80) Vp/Vs near the Atacama fault system on the coastline and the Domeyko Fault System in the forearc. The seismic velocity models also provide evidence for fluid circulation caused by the subducting Taltal ridge on the coast and partial melting feeding a volcanic complex close to the Andes. Finally, the observed low Vp/Vs ratios (∼1.75) are associated with copper mining operations in the area, suggesting that this kind of imaging can be used to characterize the distribution of potential ore deposits in the area. Plain Language Summary: We recorded earthquakes in northern Chile with a network of 84 seismometers and used the arrival times of P and S waves to generate 3D wavespeed models of the region. These models reveal several structures in the area, including changes in the angle of the subducting Nazca plate and fractures in the oceanic crust. Among features observed in both the Nazca and South American plates are the Atacama and Domeyko fault systems. We also infer fluid circulation caused by the subducting Taltal ridge and partial melting that feeds a volcanic complex near the Andes. Low values of the Vp/Vs ratio are associated with copper mining operations in the area and could be used to identify new ore deposits. Key Points: Seismic catalog reveals forearc activity and slab dip variations. Vp anomalies in the oceanic plate are related to mid‐depth seismic eventsVelocity models uncover anomalies in Salar de Atacama and Taltal ridge that might influence seismicity distribution and hydration changesShallow low Vp/Vs (<1.75) correlate with ore deposits; deep high Vp/Vs (>1.80) suggest fluids and melting for the Lastarria volcanic complex [ABSTRACT FROM AUTHOR]

This study provides a comprehensive characterization of various hydrothermal systems in Southern Peru ranging from the faulted Precordillera's steep topography up to the volcanic High Cordillera (>4,000 m asl). The objective is to investigate thermal anomalies that may potentially serve as new geothermal resources. Our integrated approach combines: (a) geochemistry from 14 hot springs sampled throughout the Tacna region, and (b) 3D numerical modeling of coupled groundwater and heat transfer considering topography and faults embedded in homogeneous permeability. Water and gas analysis indicates that the springs located near volcanoes discharge Na‐K‐Cl waters with high temperatures (>87°C), high Total Dissolved Solid concentrations (TDS >3,452 mg/L), and free gases dominated by CO2 (>90 vol%). Springs located along the regional faults in the Precordillera discharge Ca‐SO4 and Na‐K‐Cl waters with moderate temperatures (27–53°C), intermediate TDS concentrations (464–2,458 mg/L), radiocarbon ages between 1.4 and 7.9 kyr, and free gases dominated by N2 (>95 vol%). The Aruma springs, which are located at the transition between the High and the Precordillera, display intermediate characteristics. Numerical models accurately replicate the locations and temperatures of the fault‐related springs only for permeable faults (>10−14 m2), revealing the creation of 100‐km long thermal plumes along faults, locally rising up the 150°C‐isotherm to about ∼1,000 m below the surface. This approach clearly distinguishes the spring origins, which are volcanic in High Cordillera and tectonic in Precordillera. Moreover, we highlight that steep topographic gradient and permeable reverse faults in the Andean forearc may generate considerable thermal anomalies, opening perspectives for the geothermal exploration. Plain Language Summary: Geothermal energy is an essential part of the transition to green energies. It consists of using the Earth' natural heat to generate electricity or provide direct heating. In the Peruvian Andes there is significant geothermal potential near volcanoes that constitute an efficient heat source. However, these are distant from many of the population hubs, mostly located along the Pacific Coast. In the Tacna region, Southern Peru, we observe hot springs aligned along faults near the city. We hypothesize that these faults act as pipes for the hot fluids making them suitable for geothermal energy extraction. The geochemical compositions of various hot springs in the region are compared with 3D numerical models simulating hydrothermal flows. We found that springs along volcanoes and those along faults have different geochemical signatures, evidencing different subsurface pathways. Numerical models accurately replicate the positions and temperatures of the observed and sampled hot springs and reveal significant thermal anomalies around compressional faults. Our findings open up promising perspectives for the geothermal exploration in the Andean forearc. Key Points: The geochemical properties of Southern Peru's hot springs are influenced by their spatial correlation with volcanic or tectonic featuresBy combining permeable faults and topography, 3D hydrothermal modeling replicates the temperatures and locations of observed springsDetection of regional topography‐driven thermal plumes along faults in the Andean forearc opens perspectives for geothermal exploration [ABSTRACT FROM AUTHOR]

The interaction of the northern Nazca and southwestern Caribbean oceanic plates with northwestern South America (NWSA) and the collision of the Panama‐Choco arc (PCA) have significant implications on the evolution of the northern Andes. Based on a quantitative kinematic reconstruction of the Caribbean and Farallon/Farallon‐derived plates, we reconstructed the subducting geometries beneath NWSA and the PCA accretion to the continent. The persistent northeastward migration of the Caribbean plate relative to NWSA in Cenozoic time caused the continuous northward advance of the Farallon‐Caribbean plate boundary, which in turn resulted in its progressive concave trench bending against NWSA. The increasing complexity during the Paleogene included the onset of Caribbean shallow subduction, the PCA approaching the continent, and the forced shallow Farallon subduction that ended in the fragmentation of the Farallon Plate into the Nazca and Cocos plates and the Coiba and Malpelo microplates by the late Oligocene. The convergence tectonics after late Oligocene comprised the accretional process of the PCA to NWSA, which evolved from subduction erosion of the forearc to collisional tectonics by the middle Miocene, as well as changes of convergence angle and slab dip of the Farallon‐derived plates, and the attachment of the Coiba and Malpelo microplates to the Nazca plate around 9 Ma, resulting in a change of convergence directions. During the Pliocene, the Nazca slab broke at 5.5°N, shaping the modern configuration. Overall, the proposed reconstruction is supported by geophysical data and is well correlated with the magmatic and deformation history of the northern Andes. Plain Language Summary: The tectonic reconstruction in convergent triple junctions is a particularly challenging task as the relative motion between plates could define highly changing boundaries. Indeed, the resulting interaction between these convergent plates may induce important changes in the disposition of the trenches, and in turn in the three‐dimensional geometry of the subducting plates. Therefore, these highly dynamic conditions throughout geological time may be accommodated by different phases of plate fragmentation and reorganization. These factors could explain the complex spatial‐temporal distribution of subduction‐related magmatism and the different episodes of deformation in the upper plates. This reasoning is validated in the northwestern corner of South America (SA), where the continent has been converging against the Caribbean and Farallon‐derived oceanic plates since Cretaceous time. Additionally, we study the effects of the collision and accretion of the Panama‐Choco arc with SA. To accomplish that, we review the kinematic history of the Farallon/Nazca and Caribbean oceanic plates relative to stable Guiana Craton (SA) and integrate these results with the magmatic and deformation evolution of the northern Andes, which allow us to propose a model of the geometrical evolution of the subducting slabs. The obtained model is additionally constrained by seismological data and published velocity anomalies. Key Points: The tectonics of convergent triple junctions is complicated by the relative plate motion and interaction of the involved platesWe propose a model for the kinematic and geometric evolution of the Farallon/Nazca and Caribbean plates throughout the CenozoicThe interaction between the Caribbean, Nazca and South American plates is closely related to the deformation history in the Northern Andes [ABSTRACT FROM AUTHOR]

Reconstructions of motions of the Nazca, South American, and Indian plates record short‐duration (≲10 Myr) variations in angular velocity, which enable a vector‐based test of the hypothesis that mountain uplift can cause changes in plate motion. Reductions in velocity of Nazca and South America between ∼12 and 6 Ma coincide with a phase of rapid surface uplift in the Central Andes. Decrease in the rate of India's convergence with Eurasia between ∼20 and 10 Ma corresponds to an increase in gravitational potential energy per unit area (GPE) within Tibet, marked by a transition from crustal thickening to thinning. The vectorial test shows that, in each case, the only change in driving force capable of balancing the change in basal drag is an increased resistance along the convergent boundary to the plate. Changes in GPE associated with mountain uplift provide a calibration for basal drags on plates. Basal tractions of ∼0.1–1 MPa provide resisting forces comparable in magnitude to driving forces from GPE variation in ocean lithosphere. The rapid change in motion of the Indian plate between about 48 and 41 Ma is explained by the juxtaposition of the Indian continent against the Andean‐type margin of the Transhimalaya and reduction in driving force due to loss of the slab. The net slab driving force lost was ∼2–4 TN m−1, in agreement with previous studies suggesting that forces resisting slabs' penetration into the mantle largely offset their negative buoyancy. Key Points: Calculation of forces required to cause rapid changes in plate motion test the hypothesis that some stem from increase in gravitational potential energy (GPE) of mountainsRapid decreases in velocity of Nazca, South American, and Indian plates require increases in resistance at their convergent boundariesGeological evidence shows that the Miocene decreases in velocity were contemporaneous with increases in GPE of the Andes and Tibet [ABSTRACT FROM AUTHOR]

Tephra layers preserved in marine sediments are strong tools to study the frequency, magnitude and source of past major explosive eruptions. Thirty‐seven volcanoes from the Ecuadorian and Colombian arc, in the northern Andes, experienced at least one eruption during the Holocene. The volcanic hazard is therefore particularly high for the populated areas of the Andes and in particular cases for the coastal region, and it is crucial to document such events to improve hazard assessment. The age and distribution of deposits from major Holocene eruptions have been studied in the Cordillera, but no descriptions of distal fallouts have been published. In this study, we focused on 28 Holocene tephra layers recorded in marine sediment cores collected along the northern Ecuador—Southern Colombia margin. New lithological, geochemical and isotope data together with 14C datings on foraminifers allow us to determine the age and volcanic source of marine tephra, and to propose a first land‐sea correlation of distal tephra fallouts. We show that at least seven explosive eruptions from Guagua Pichincha, Atacazo‐Ninahuilca, Cotopaxi, and Cerro Machín volcanoes left tephra deposits recorded in marine cores over 250 km away from their source. Volume estimates of emitted tephra range between 1.3 and 6.0 km3 for the tenth century Guagua Pichincha, ∼5 ka Atacazo‐Ninahuilca, ∼6.7 and ∼7.9 ka Cotopaxi events, suggesting that they were eruptions of Volcanic Explosivity Index of 5. The distribution of these deposits also brings new constraints for a better evaluation of the volcanic hazard in Ecuador. Plain Language Summary: During major explosive eruptions, large volumes of gases and tephra (lapilli and ash particles) are thrown into the atmosphere and can be spread by winds over 100 km and more. Tephra fallouts can impact the population, infrastructures and climate. It is therefore essential to document the age and magnitude of past major eruptions to better assess the volcanic hazards. In this study, we use the mineralogy, glass shard morphology, and the geochemical composition of tephra settled in marine sediments off Ecuador and Colombia to investigate their source. Thickness of tephra layers and radiocarbon ages performed on under‐ and over‐lying marine fauna allow us to determine the age of the eruptions, whereas the distribution of tephra yields constraints on the volume of fallout deposits. We show that the largest explosive eruptions from Ecuadorian and Colombian volcanoes reached the Pacific Ocean with a recurrence rate of about 1.5 events per millennium over the past 8 kyr. Key Points: We propose a first land‐sea correlation of distal Holocene tephra off Ecuador based on 14C age and geochemical dataProducts from at least seven explosive Holocene eruptions in Ecuador and south Colombia reached the Pacific OceanVolumes of tephra emitted by largest eruptions vary between 1.3 and 6.0 km3, suggesting they were VEI‐5 eruptions [ABSTRACT FROM AUTHOR]

New earthquake focal mechanism and centroid depth estimates show that the deformation style in the forelands of the Andes is spatially correlated with rift systems that stretched the South American lithosphere in the Mesozoic. Where the rifts trend sub‐parallel to the Andean range front, normal faults inherited from the rifts are being reactivated as reverse faults, causing the 30–45 km thick seismogenic layer to break up. Where the rift systems are absent from beneath the range front, the seismogenic layer is bending and being thrust beneath the Andes like a rigid plate. Force‐balance calculations show that the faults inerhited from former rift zones have an effective coefficient of static friction μ′ < 0.2. In order for these frictionally weak faults to remain seismogenic in the lower crust, their wall rocks are likely to be formed of dry granulite. Xenolith data support this view, and suggest that parts of the lower crust are now mostly metastable, having experienced temperatures at least 75–250°C hotter than present. The conditions in the lower crust make it unlikely that highly pressurized free water, or networks of intrinsically weak phyllosilicate minerals, are the cause of their low effective friction, as, at such high temperatures, both mechanisms would cause the faults to deform through viscous creep and not frictional slip. Therefore pre‐existing faults in the Andean forelands have remained weak and seismogenic after reactivation, and have influenced the style of mountain building in South America. However, the controls on their mechanical properties in the lower crust remain unclear. Plain Language Summary: This study is concerned with the controls on how mountain ranges grow. I show that the locations and types of earthquakes along the margins of the Andes mountains, which are generated when the mountain range grows, vary systematically with the positions of ancient fault zones. Where ancient faults lie along the margins of the Andes the entire crust is breaking up through slip on these faults. Where the same faults are not present along the margins of the Andes, the crust is being pushed beneath the mountains like a rigid plate. Therefore, the strength of faults along the margins of the Andes play a critical role in the growth of the mountain range. Notably, the earthquake‐generating faults along the margins of the Andes are much weaker than predicted by laboratory experiments, and the physical reason for their weakness remains unclear. Key Points: New earthquake data show the Andean forelands are breaking up in compression to 30–45 km depth in areas that experienced Mesozoic riftingForce‐balance calculations demonstrate that the effective coefficients of static friction on faults inherited from the rifts is <0.2These frictionally weak, seismogenic faults control the style of active mountain building in the forelands [ABSTRACT FROM AUTHOR]

We examined the relationship between the mechanical strength of the lithosphere and the distribution of seismicity within the overriding continental plate of the southern Central Andes (SCA, 29°–39°S), where the oceanic Nazca Plate changes its subduction angle between 33°S and 35°S, from subhorizontal in the north (<5°) to steep in the south (∼30°). We computed the long‐term lithospheric strength based on an existing 3D model describing variations in thickness, density, and temperature of the main geological units forming the lithosphere of the SCA and adjacent forearc and foreland regions. The comparison between our results and seismicity within the overriding plate (upper‐plate seismicity) shows that most of the events occur within the modeled brittle domain of the lithosphere. The depth where the deformation mode switches from brittle frictional to thermally activated ductile creep provides a conservative lower bound to the seismogenic zone in the overriding plate of the study area. We also found that the majority of upper‐plate earthquakes occurs within the realm of first‐order contrasts in integrated strength (12.7–13.3 log Pam in the Andean orogen vs. 13.5–13.9 log Pam in the forearc and the foreland). Specific conditions characterize the mechanically strong northern foreland of the Andes, where seismicity is likely explained by the effects of slab steepening. Plain Language Summary: The southern Central Andes (29°–39°S) are one of the most seismically active regions along the South‐American subduction zone. However, the causative dynamics behind localization of earthquakes in the South American plate are still not well understood. In this study, we computed the strength of the lithosphere (i.e., the resistance to deformation) with a specialized computer code and showcase that there exist a causative relationship with seismicity in the area. We found that earthquakes in the upper plate occur: (a) above the transition between domains with brittle to ductile deformation; and (b) where large and sharp horizontal strength contrasts exist. Specific conditions characterize the mechanically strong northern foreland of the Andes, where seismicity is likely explained by the additional effects of the oceanic subducted plate, where it transitions from a subhorizontal to a steep subduction angle. These results show the importance of quantifying the strength of the lithosphere and its spatial variations to better understand the distribution of seismicity records in the area. Key Points: The yield strength of the lithosphere was computed and compared with the upper‐plate seismicity distribution of the southern Central AndesThe brittle‐ductile transition effectively bounds the depth extent of the seismogenic zoneSeismicity localizes at the boundaries between weak and strong domains; slab steepening likely explains seismicity in the strong foreland [ABSTRACT FROM AUTHOR]