Miocene surface uplift and orogenic evolution of the southern Andean Plateau (central Puna), northwestern Argentina

We present stable hydrogen-isotope analyses of volcanic glass ([Formula: see text]Dg) and radiometric ages (U-Pb zircon, U-Th calcite, AMS14C) from deformed sedimentary deposits in the vicinity of the intermontane Pocitos Basin in the central Puna of the Andean Plateau at about 24.5°S. Our results demonstrate 2-km surface uplift since the middle to late Miocene and protracted shortening that persists until the present day, while other sectors of the Puna show evidence for tectonically neutral and/or extensional settings. These findings are at odds with previous studies suggesting near-modern elevations (4 km) of the Puna Plateau since the late Eocene and formation of the intermontane Miocene Arizaro-Pocitos Basin associated with gravitational foundering of a dense lithosphere. Geophysical and geochemical data support the removal of continental lithosphere beneath the Puna, but the timing and mechanisms by which this removal occurs have remained controversial. We hypothesize that intermontane basin formation in the central Puna is the result of crustal shortening since about 20 Ma, followed by rapid surface uplift, likely related to lithospheric delamination.

Late Cretaceous through Cenozoic Paleoenvironmental History of the Bagua Basin, Peru: Paleoelevation Comparisons with the Central Andean Plateau

Located in northern Peru, at the lowest segment of the Central Andes, the Bagua Basin contains a Campanian to Pleistocene sedimentary record that archives the local paleoenvironmental and tectonic history. We present new δ18O and δ13C signatures of pedogenic carbonate nodules from paleosols in the Campanian–Maastrichtian Fundo El Triunfo Formation and in the upper Eocene–middle Miocene Sambimera Formation to reconstruct the isotopic composition of paleo-meteoric water and the floristic biome. We compare these results to modern isotopic values from a newly obtained modern water transect to interpret the environmental evolution of this area and its relationship with the neighboring Eastern Cordillera. A ~2‰ δ18O depletion between the latest Cretaceous and the latest Eocene reflects a shift from a coastal to inland environment. A negative δ18O shift of ~3‰ from the middle Miocene to the present day reveals the establishment of the Eastern Cordillera as an orographic barrier for the moisture traveling westward, sometime after deposition of the top of the Sambimera Formation at ~13 Ma. A shift in the δ13C signature from ~−25‰ in the Campanian–Miocene deposits to ~−23‰ in modern–Holocene times suggests a change in biome from dominant C3 plants to a mixture of C3 and C4 plants. This environmental shift reflects both the late Miocene global C4 expansion and the transition to more arid conditions in the basin. The Campanian–middle Miocene environmental reconstruction of the Bagua Basin indicates a steady paleoelevation setting in the northernmost Central Andes during most of the Cenozoic and constrains the uplift of the Eastern Cordillera to the late Miocene–Pleistocene. This paleoelevation history contrasts with that of the Central Andean Plateau, which is characterized by two major episodes of surface uplift: early–middle Miocene and late Miocene–Pliocene. The contrasting modern topographic configuration of the Central Andean Plateau and the northernmost Central Andes gives rise to the question of what factors created such a dramatic difference in topographic evolution of the two regions that shared an overall common tectonic history. We discuss the possible factors responsible for this contrasting topographic configuration and suggest that the diachronous flat slab episodes are likely a major factor, resulting in greater shortening and crustal thickness and, ultimately, in earlier surface uplift episodes occurring in the Central Andean Plateau.

Chemical Mohometry: Assessing Crustal Thickness of Ancient Orogens Using Geochemical and Isotopic Data

Convergent plate boundaries are key sites for continental crustal formation and recycling. Quantifying the evolution of crustal thickness and paleoelevation along ancient convergent margins represents a major goal in orogenic system analyses. Chemical and in some cases isotopic compositions of igneous rocks formed in modern supra-subduction arcs and collisional belts are sensitive to Moho depths at the location of magmatism, implying that igneous suites from fossil orogens carry information about crustal thickness from the time they formed. Several whole-rock chemical parameters correlate with crustal thickness, some of which were calibrated to serve as "mohometers," that is, quantitative proxies of paleo-Moho depths. Based on mineral-melt partition coefficients, this concept has been extended to detrital zircons, such that combined chemical and geochronological information extracted from these minerals allows us to reconstruct the crustal thickness evolution using the detrital archive. We discuss here the mohometric potential of a variety of chemical and isotopic parameters and show that their combined usage improves paleocrustal thickness estimates. Using a MATLAB® app developed for the underlying computations, we present examples from the modern and the deeper time geologic record to illustrate the promises and pitfalls of the technique. Since arcs are in isostatic equilibrium, mohometers are useful in reconstructing orogenic paleoelevation as well. Our analysis suggests that many global-scale correlations between magma composition and crustal thickness used in mohometry originate in the sub-arc mantle; additional effects resulting from intracrustal igneous differentiation depend on the compatible or incompatible behavior of the involved parameters.

Growth of Neogene Andes linked to changes in plate convergence using high-resolution kinematic models

The Andean cordillera was constructed during compressive tectonic events, whose causes and controls remain unclear. Exploring a possible link to plate convergence has been impeded by the coarse temporal resolution of existing plate kinematic models. Here we show that the Neogene evolution of the Andean margin is primarily related to changes in convergence as observed in new high-resolution plate reconstructions. Building on a compilation of plate finite rotations spanning the last 30 million years and using noise-mitigation techniques, we predict several short-term convergence changes that were unresolved in previous models. These changes are related to main tectono-magmatic events and require forces that are compatible with a range of geodynamic processes. These results allow to revise models of ongoing subduction orogeny at its type locality, emphasizing the role of upper plate deformation in the balance between kinematic energy associated with plate motion and gravitational potential energy stored in orogenic crustal roots.

A high-resolution model of the motion between Nazca and South American plates is presented. The work shows rapid changes that help explaining tectono-magmatic events via a balance between kinematic energy and gravitational potential energy stored in the roots of the Andes.

What controls the remobilization and deformation of surficial sediment by seismic shaking? Linking lacustrine slope stratigraphy to great earthquakes in South-Central Chile

Remobilization and deformation of surficial subaqueous slope sediments create turbidites and soft sediment deformation structures, which are common features in many depositional records. Palaeoseismic studies have used seismically-induced turbidites and soft sediment deformation structures preserved in sedimentary sequences to reconstruct recurrence patterns and - in some cases - allow quantifying rupture location and magnitude of past earthquakes. However, current understanding of earthquake-triggered remobilization and deformation lacks studies targeting where these processes take place, the subaqueous slope and involving direct comparison of sedimentary fingerprint with well-documented historical earthquakes. This study investigates the sedimentary imprint of six megathrust earthquakes with varying rupture characteristics in 17 slope sediment cores from two Chilean lakes, Riñihue and Calafquén, and evaluates how it links to seismic intensity, peak ground acceleration, bracketed duration and slope angle. Centimetre-scale stratigraphic gaps ranging from ca 1 to 20 cm - caused by remobilization of surficial slope sediment - were identified using high-resolution multi-proxy core correlation of slope to basin cores, and six types of soft sediment deformation structures ranging from ca 1 to 25 cm thickness using high-resolution three-dimensional X-ray computed tomography data. Stratigraphic gaps occur on slope angles of ≥2.3°, whereas deformation already occurs from slope angle 0.2°. The thickness of both stratigraphic gaps and soft sediment deformation structures increases with slope angle, suggesting that increased gravitational shear stress promotes both surficial remobilization and deformation. Seismic shaking is the dominant trigger for surficial remobilization and deformation at the studied lakes. Total remobilization depth correlates best with bracketed duration and is highest in both lakes for the strongest earthquakes (M _w ca 9.5). In lake Riñihue, soft sediment deformation structure thickness and type correlate best with peak ground acceleration providing the first field-based evidence of progressive soft sediment deformation structure development with increasing peak ground acceleration for soft sediment deformation structures caused by Kelvin-Helmholtz instability. The authors propose that long duration and low frequency content of seismic shaking favours surficial remobilization, whereas ground motion amplitude controls Kelvin-Helmholtz instability-related soft sediment deformation structure development.

Out-of-sequence skeletal growth causing oscillatory zoning in arc olivines

Primitive olivines from the monogenetic cones Los Hornitos, Central-South Andes, preserve dendritic, skeletal, and polyhedral growth textures. Consecutive stages of textural maturation occur along compositional gradients where high Fo–Ni cores of polyhedral olivines (Fo_92.5, Ni ~3500 ppm) contrast with the composition of dendritic olivines (Fo < 91.5, Ni < 3000 ppm), indicating sequential nucleation. Here we present a new growth model for oscillatory Fo–Ni olivine zoning that contrasts with the standard interpretation of continuous, sequential core-to-rim growth. Olivine grows rapidly via concentric addition of open-structured crystal frames, leaving behind compositional boundary layers that subsequently fill-in with Fo–Ni-depleted olivine, causing reversals. Elemental diffusion modeling reveals growth of individual crystal frames and eruption at the surface occurred over 3.5–40 days. Those timescales constrain magma ascent rates of 40–500 m/h (0.011 to 0.14 m/s) from the deep crust. Compared to ocean island basalts, where dendritic and skeletal olivines have been often described, magmas erupted at arc settings, experiencing storage and degassing, may lack such textures due to fundamentally different ascent histories.

Arc olivines are commonly explained through a paradigm of core-to-rim sequential growth and oscillatory zoning is interpreted to represent magma mixing. Here the authors show Fo–Ni–P oscillatory zoned olivines can grow as out-of-sequence crystal frames and complex zoning can occur in closed systems.

The role of ridges in the formation and longevity of flat slabs

Flat-slab subduction occurs when the descending plate becomes horizontal at some depth before resuming its descent into the mantle. It is often proposed as a mechanism for the uplifting of deep crustal rocks ('thick-skinned' deformation) far from plate boundaries, and for causing unusual patterns of volcanism, as far back as the Proterozoic eon. For example, the formation of the expansive Rocky Mountains and the subsequent voluminous volcanism across much of the western USA has been attributed to a broad region of flat-slab subduction beneath North America that occurred during the Laramide orogeny (80-55 million years ago). Here we study the largest modern flat slab, located in Peru, to better understand the processes controlling the formation and extent of flat slabs. We present new data that indicate that the subducting Nazca Ridge is necessary for the development and continued support of the horizontal plate at a depth of about 90 kilometres. By combining constraints from Rayleigh wave phase velocities with improved earthquake locations, we find that the flat slab is shallowest along the ridge, while to the northwest of the ridge, the slab is sagging, tearing, and re-initiating normal subduction. On the basis of our observations, we propose a conceptual model for the temporal evolution of the Peruvian flat slab in which the flat slab forms because of the combined effects of trench retreat along the Peruvian plate boundary, suction, and ridge subduction. We find that while the ridge is necessary but not sufficient for the formation of the flat slab, its removal is sufficient for the flat slab to fail. This provides new constraints on our understanding of the processes controlling the beginning and end of the Laramide orogeny and other putative episodes of flat-slab subduction.

2010 Maule earthquake slip correlates with pre-seismic locking of Andean subduction zone

The magnitude-8.8 Maule (Chile) earthquake of 27 February 2010 ruptured a segment of the Andean subduction zone megathrust that has been suspected to be of high seismic potential. It is the largest earthquake to rupture a mature seismic gap in a subduction zone that has been monitored with a dense space-geodetic network before the event. This provides an image of the pre-seismically locked state of the plate interface of unprecedentedly high resolution, allowing for an assessment of the spatial correlation of interseismic locking with coseismic slip. Pre-seismic locking might be used to anticipate future ruptures in many seismic gaps, given the fundamental assumption that locking and slip are similar. This hypothesis, however, could not be tested without the occurrence of the first gap-filling earthquake. Here we show evidence that the 2010 Maule earthquake slip distribution correlates closely with the patchwork of interseismic locking distribution as derived by inversion of global positioning system (GPS) observations during the previous decade. The earthquake nucleated in a region of high locking gradient and released most of the stresses accumulated in the area since the last major event in 1835. Two regions of high seismic slip (asperities) appeared to be nearly fully locked before the earthquake. Between these asperities, the rupture bridged a zone that was creeping interseismically with consistently low coseismic slip. The rupture stopped in areas that were highly locked before the earthquake but where pre-stress had been significantly reduced by overlapping twentieth-century earthquakes. Our work suggests that coseismic slip heterogeneity at the scale of single asperities should indicate the seismic potential of future great earthquakes, which thus might be anticipated by geodetic observations.