@article{ArayaVargasMeqbelRitteretal.2019,
  author    = {Araya Vargas, Jaime Andr{\´e}s and Meqbel, Naser M. and Ritter, Oliver and Brasse, H. and Weckmann, Ute and Yanez, Gonzalo and Godoy, B.},
  title     = {Fluid Distribution in the Central Andes Subduction Zone Imaged With Magnetotellurics},
  series = {Journal of geophysical research : Solid earth},
  volume    = {124},
  journal   = {Journal of geophysical research : Solid earth},
  number    = {4},
  publisher = {American Geophysical Union},
  address   = {Washington},
  issn      = {2169-9313},
  doi       = {10.1029/2018JB016933},
  pages     = {4017 -- 4034},
  year      = {2019},
  abstract  = {We present a model of the electrical resistivity structure of the lithosphere in the Central Andes between 20 degrees and 24 degrees S from 3-D inversion of 56 long-period magnetotelluric sites. Our model shows a complex resistivity structure with significant variability parallel and perpendicular to the trench direction. The continental forearc is characterized mainly by high electrical resistivity (>1,000m), suggesting overall low volumes of fluids. However, low resistivity zones (LRZs, <5m) were found in the continental forearc below areas where major trench-parallel faults systems intersect NW-SE transverse faults. Forearc LRZs indicate circulation and accumulation of fluids in highly permeable fault zones. The continental crust along the arc shows three distinctive resistivity domains, which coincide with segmentation in the distribution of volcanoes. The northern domain (20 degrees-20.5 degrees S) is characterized by resistivities >1,000m and the absence of active volcanism, suggesting the presence of a low-permeability block in the continental crust. The central domain (20.5 degrees-23 degrees S) exhibits a number of LRZs at varying depths, indicating different levels of a magmatic plumbing system. The southern domain (23 degrees-24 degrees S) is characterized by resistivities >1,000m, suggesting the absence of large magma reservoirs below the volcanic chain at crustal depths. Magma reservoirs located below the base of the crust or in the backarc may fed active volcanism in the southern domain. In the subcontinental mantle, the model exhibits LRZs in the forearc mantle wedge and above clusters of intermediate-depth seismicity, likely related to fluids produced by serpentinization of the mantle and eclogitization of the slab, respectively.},
  language  = {en}
}
@article{CastinoBookhagenStrecker2017,
  author    = {Castino, Fabiana and Bookhagen, Bodo and Strecker, Manfred},
  title     = {Rainfall variability and trends of the past six decades (1950-2014) in the subtropical NW Argentine Andes},
  series = {Climate dynamics : observational, theoretical and computational research on the climate system},
  volume    = {48},
  journal   = {Climate dynamics : observational, theoretical and computational research on the climate system},
  publisher = {Springer},
  address   = {New York},
  issn      = {0930-7575},
  doi       = {10.1007/s00382-016-3127-2},
  pages     = {1049 -- 1067},
  year      = {2017},
  abstract  = {The eastern flanks of the Central Andes are characterized by deep convection, exposing them to hydrometeorological extreme events, often resulting in floods and a variety of mass movements. We assessed the spatiotemporal pattern of rainfall trends and the changes in the magnitude and frequency of extreme events (ae<yen>95th percentile) along an E-W traverse across the southern Central Andes using rain-gauge and high-resolution gridded datasets (CPC-uni and TRMM 3B42 V7). We generated different climate indices and made three key observations: (1) an increase of the annual rainfall has occurred at the transition between low (< 0.5 km) and intermediate (0.5-3 km) elevations between 1950 and 2014. Also, rainfall increases during the wet season and, to a lesser degree, decreases during the dry season. Increasing trends in annual total amounts characterize the period 1979-2014 in the arid, high-elevation southern Andean Plateau, whereas trend reversals with decreasing annual total amounts were found at low elevations. (2) For all analyzed periods, we observed small or no changes in the median values of the rainfall-frequency distribution, but significant trends with intensification or attenuation in the 95th percentile. (3) In the southern Andean Plateau, extreme rainfall events exhibit trends towards increasing magnitude and, to a lesser degree, frequency during the wet season, at least since 1979. Our analysis revealed that low (< 0.5 km), intermediate (0.5-3 km), and high-elevation (> 3 km) areas respond differently to changing climate conditions, and the transition zone between low and intermediate elevations is characterized by the most significant changes.},
  language  = {en}
}
@article{CastinoBookhagenStrecker2017,
  author    = {Castino, Fabiana and Bookhagen, Bodo and Strecker, Manfred},
  title     = {Oscillations and trends of river discharge in the southern Central Andes and linkages with climate variability},
  series = {Journal of hydrology},
  volume    = {555},
  journal   = {Journal of hydrology},
  publisher = {Elsevier},
  address   = {Amsterdam},
  issn      = {0022-1694},
  doi       = {10.1016/j.jhydrol.2017.10.001},
  pages     = {108 -- 124},
  year      = {2017},
  abstract  = {This study analyzes the discharge variability of small to medium drainage basins (10(2)-10(4) km(2)) in the southern Central Andes of NW Argentina. The Hilbert-Huang Transform (HHT) was applied to evaluate non-stationary oscillatory modes of variability and trends, based on four time series of monthly normalized discharge anomaly between 1940 and 2015. Statistically significant trends reveal increasing discharge during the past decades and document an intensification of the hydrological cycle during this period. An Ensemble Empirical Mode Decomposition (EEMD) analysis revealed that discharge variability in this region can be best described by five quasi-periodic statistically significant oscillatory modes, with mean periods varying from 1 to 20 y. Moreover, we show that discharge variability is most likely linked to the phases of the Pacific Decadal Oscillation (PDO) at multi-decadal timescales (similar to 20 y) and, to a lesser degree, to the Tropical South Atlantic SST anomaly (TSA) variability at shorter timescales (similar to 2-5 y). Previous studies highlighted a rapid increase in discharge in the southern Central Andes during the 1970s, inferred to have been associated with the global 1976-77 climate shift. Our results suggest that the rapid discharge increase in the NW Argentine Andes coincides with the periodic enhancement of discharge, which is mainly linked to a negative to positive transition of the PDO phase and TSA variability associated with a long-term increasing trend. We therefore suggest that variations in discharge in this region are largely driven by both natural variability and the effects of global climate change. We furthermore posit that the links between atmospheric and hydrologic processes result from a combination of forcings that operate on different spatiotemporal scales. (C) 2017 Elsevier B.V. All rights reserved.},
  language  = {en}
}
@article{GuzmanPetrinovicBrodetal.2011,
  author    = {Guzman, S. and Petrinovic, I. A. and Brod, J. A. and Hongn, Fernando D. and Seggiaro, R. E. and Montero, C. and Carniel, Roberto and Dantas, E. L. and Sudo, Masafumi},
  title     = {Petrology of the Luingo caldera (SE margin of the Puna plateau) a middle Miocene window of the arc-back arc configuration},
  series = {Journal of volcanology and geothermal research},
  volume    = {200},
  journal   = {Journal of volcanology and geothermal research},
  number    = {3-4},
  publisher = {Elsevier},
  address   = {Amsterdam},
  issn      = {0377-0273},
  doi       = {10.1016/j.jvolgeores.2010.12.008},
  pages     = {171 -- 191},
  year      = {2011},
  abstract  = {We describe the petrographic characteristics, whole-rock geochemistry and mineral chemistry of rocks from the Pucarilla-Cerro Tipillas Volcanic Complex with emphasis on the rocks belonging to the middle Miocene Luingo caldera, located in the south-eastern portion of the Central Volcanic Zone (CVZ) of the Andes. We modelled the petrogenesis of the Luingo caldera rocks as a mixture of ca. 20\% crustal magmas and 80\% of mantle magmas by AFC with recharge processes. A comparison of Luingo geochemical data with the composition of Miocene-Pliocene volcanic rocks from the broad area, points to major thickening events during the middle Miocene for the western portion and during the upper Miocene for the eastern portion of the Southern CVZ. In the eastern sector (similar to 66 degrees W) the mantle source appears to change from a spinel-lherzolite type for the middle Miocene to a garnet-lherzolite type for the upper Miocene-Pliocene magmas. The areal distribution of the volcanic products led to the recognition of approximately equivalent areas covered by volcanic rocks both in the eastern and in the western Puna borders. This indicates a broad arc, which was structurally controlled at the proto-Puna/Puna margins, whose geochemical differences are related with variations in crustal thicknesses and heterogeneous mantle sources from west to east.},
  language  = {en}
}
@article{IbarraLiuMeessenetal.2019,
  author    = {Ibarra, Federico and Liu, Sibiao and Meeßen, Christian and Prezzi, Claudia Beatriz and Bott, Judith and Scheck-Wenderoth, Magdalena and Sobolev, Stephan Vladimir and Strecker, Manfred},
  title     = {3D data-derived lithospheric structure of the Central Andes and its implications for deformation: Insights from gravity and geodynamic modelling},
  series = {Tectonophysics : international journal of geotectonics and the geology and physics of the interior of the earth},
  volume    = {766},
  journal   = {Tectonophysics : international journal of geotectonics and the geology and physics of the interior of the earth},
  publisher = {Elsevier},
  address   = {Amsterdam},
  issn      = {0040-1951},
  doi       = {10.1016/j.tecto.2019.06.025},
  pages     = {453 -- 468},
  year      = {2019},
  abstract  = {We present a new three-dimensional density model of the Central Andes characterizing the structure and composition of the lithosphere together with a geodynamic simulation subjected to continental intraplate shortening. The principal aim of this study is to assess the link between heterogeneities in the lithosphere and different deformation patterns and styles along the orogen-foreland system of the Central Andes. First, we performed a 3D integration of new geological and geophysical data with previous models through forward modelling of Bouguer anomalies. Subsequently, a geodynamic model was set-up and parametrized from the previously obtained 3D structure and composition. We do not find a unambigous correlation between the resulting density configuration and terrane boundaries proposed by other authors. Our models reproduce the observed Bouguer anomaly and deformation patterns in the foreland. We find that thin-skinned deformation in the Subandean fold-and thrust belt is controlled by a thick sedimentary layer and coeval underthrusting of thin crust of the foreland beneath the thick crust of the Andean Plateau. In the adjacent thick-skinned deformation province of the inverted Cretaceous extensional Santa Barbara System sedimentary strata are much thinner and crustal thickness transitions from greater values in the Andean to a more reduced thickness in the foreland. Our results show that deformation processes occur where the highest gradients of lithospheric strength are present between the orogen and the foreland, thus suggesting a spatial correlation between deformation and lithospheric strength.},
  language  = {en}
}
@article{LiuSobolevBabeykoetal.2022,
  author    = {Liu, Sibiao and Sobolev, Stephan and Babeyko, Andrey and Pons, Micha{\"e}l},
  title     = {Controls of the foreland deformation pattern in the orogen-foreland shortening system},
  series = {Tectonics},
  volume    = {41},
  journal   = {Tectonics},
  number    = {2},
  publisher = {American Geophysical Union},
  address   = {Washington},
  issn      = {0278-7407},
  doi       = {10.1029/2021TC007121},
  pages     = {18},
  year      = {2022},
  abstract  = {Controls on the deformation pattern (shortening mode and tectonic style) of orogenic forelands during lithospheric shortening remain poorly understood. Here, we use high-resolution 2D thermomechanical models to demonstrate that orogenic crustal thickness and foreland lithospheric thickness significantly control the shortening mode in the foreland. Pure-shear shortening occurs when the orogenic crust is not thicker than the foreland crust or thick, but the foreland lithosphere is thin (<70-80 km, as in the Puna foreland case). Conversely, simple-shear shortening, characterized by foreland underthrusting beneath the orogen, arises when the orogenic crust is much thicker. This thickened crust results in high gravitational potential energy in the orogen, which triggers the migration of deformation to the foreland under further shortening. Our models present fully thick-skinned, fully thin-skinned, and intermediate tectonic styles in the foreland. The first tectonics forms in a pure-shear shortening mode whereas the others require a simple-shear mode and the presence of thick (>similar to 4 km) sediments that are mechanically weak (friction coefficient <similar to 0.05) or weakened rapidly during deformation. The formation of fully thin-skinned tectonics in thick and weak foreland sediments, as in the Subandean Ranges, requires the strength of the orogenic upper lithosphere to be less than one-third as strong as that of the foreland upper lithosphere. Our models successfully reproduce foreland deformation patterns in the Central and Southern Andes and the Laramide province.},
  language  = {en}
}
@article{NouryBernetSchildgenetal.2016,
  author    = {Noury, M. and Bernet, M. and Schildgen, Taylor F. and Simon-Labric, T. and Philippon, M. and Sempere, T.},
  title     = {Crustal-scale block tilting during Andean trench-parallel extension: Structural and geo-thermochronological insights},
  series = {Tectonics},
  volume    = {35},
  journal   = {Tectonics},
  publisher = {American Geophysical Union},
  address   = {Washington},
  issn      = {0278-7407},
  doi       = {10.1002/2016TC004231},
  pages     = {2052 -- 2069},
  year      = {2016},
  abstract  = {Despite a long history of plate convergence at the western margin of the South American plate that has been ongoing since at least the Early Paleozoic, the southern Peruvian fore-arc displays little to no evidence of shortening. In the light of this observation, we assess the deformation history of the southern Peruvian fore-arc and its geodynamic implications. To accomplish this, we present a new structural and geo-thermochronological data set (zircon U-Pb, mica Ar-40/Ar-39, apatite and zircon fission-track and zircon (U-Th)/He analyses) for samples collected along a 400km long transect parallel to the trench. Our results show that the Mesoproterozoic gneissic basement was mainly at temperatures 350 degrees C since the Neoproterozoic and was later intruded by Jurassic volcanic arc plutons. Along the coast, a peculiar apatite fission-track age pattern, coupled with field observations and a synthesis of available geological maps, allows us to identify crustal-scale tilted blocks that span the coastal Peruvian fore-arc. These blocks, bounded by normal faults that are orthogonal to the trench, suggest post-60Ma trench-parallel extension that potentially accommodated oroclinal bending in this region. Block tilting is consistent with the observed and previously described switch in the location of sedimentary sources in the fore-arc basin. Our data set allows us to estimate the cumulative slip on these faults to be less than 2km and questions the large amount of trench-parallel extension suggested to have accommodated this bending.},
  language  = {en}
}
@article{OrtizSaezAlvaradoetal.2022,
  author    = {Ortiz, Gustavo and Saez, Mauro and Alvarado, Patricia and Rivas, Carolina and Garc{\´i}a, V{\´i}ctor Hugo and Alonso, Ricardo and Zullo, Fernando Morales},
  title     = {Seismotectonic characterization of the 1948 (M-W 6.9) Anta earthquake Santa Barbara System, central Andes broken foreland of northwestern Argentina},
  series = {Journal of South American earth sciences},
  volume    = {116},
  journal   = {Journal of South American earth sciences},
  publisher = {Elsevier},
  address   = {Oxford},
  issn      = {0895-9811},
  doi       = {10.1016/j.jsames.2022.103822},
  pages     = {15},
  year      = {2022},
  abstract  = {The region of the Andean back-arc of northwestern Argentina has been struck by several magnitude >= 6 crustal earthquakes since the first historically recorded event in 1692. One of these events corresponds to the Anta earthquake on 25 August 1948, with epicenter in the Santa Barbara System causing three deaths and severe damage in Salta and Jujuy provinces with maximum Modified Mercalli seismic intensities (MMI) of IX. We collected and digitized analog seismograms of this earthquake from worldwide seismic observatories in order to perform first-motion analysis and modeling of long-period teleseismic P-waveforms. Our results indicate a simple seismic source of M0 = 2.85 x 1019 N m consistent with a moment magnitude Mw = 6.9. We have also tested for the focal depth determining a shallow source at 8 km with a reverse focal mechanism solution with a minor dextral strike-slip component (strike 20 degrees, dip 30 degrees, rake 120 degrees) from the best fit of waveforms. Using magnitude size empirical relationships, the comparison of the obtained Mw 6.9 magnitude value and the ca. 10,000 km2 area of MMI >= IX from our seismic intensity map, which was obtained from newspaper and many historical reports, indicates a rupture length of 42 +/- 8 km for the Anta earthquake. We show our results in a 3D geological model around the epicentral area, which integrates modern seismicity, geological data, and information of a previously studied east-west cross section located a few kilometers south of the 1948 epicenter. The integration of all available information provides evidence of the re-activation of the Pie de la Sierra del Gallo fault during the 1948 Mw 6.9 shallow earthquake; this thrust fault bounds the Santa Barbara System along its western foothill.},
  language  = {en}
}
@article{SchildgenvanderBeekD'Arcyetal.2022,
  author    = {Schildgen, Taylor F. and van der Beek, Peter A. and D'Arcy, Mitch and Roda-Boluda, Duna N. and Orr, Elizabeth N. and Wittmann, Hella},
  title     = {Quantifying drainage-divide migration from orographic rainfall over geologic timescales},
  series = {Earth \& planetary science letters},
  volume    = {579},
  journal   = {Earth \& planetary science letters},
  publisher = {Elsevier},
  address   = {Amsterdam},
  issn      = {0012-821X},
  doi       = {10.1016/j.epsl.2021.117345},
  pages     = {13},
  year      = {2022},
  abstract  = {Drainage-divide migration, controlled by rock-uplift and rainfall patterns, may play a major role in the geomorphic evolution of mountain ranges. However, divide-migration rates over geologic timescales have only been estimated by theoretical studies and remain empirically poorly constrained. Geomorphological evidence suggests that the Sierra de Aconquija, on the eastern side of the southern Central Andes, northwest Argentina, is undergoing active westward drainage-divide migration. The mountain range has been subjected to steep rock trajectories and pronounced orographic rainfall for the last several million years, presenting an ideal setting for using low-temperature thermochronometric data to explore its topographic evolution. We perform three-dimensional thermal-kinematic modeling of previously published thermochronometric data spanning the windward and leeward sides of the range to explore the most likely structural and topographic evolution of the range. We find that the data can be explained by scenarios involving drainage-divide migration alone, or by scenarios that also involve changes in the structures that have accommodated deformation through time. By combining new Be-10-derived catchment-average denudation rates with geomorphic constraints on probable fault activity, we conclude that the evolution of the range was likely dominated by west-vergent faulting on a high-angle reverse fault underlying the range, together with westward drainage-divide migration at a rate of several km per million years. Our findings place new constraints on the magnitudes and rates of drainage-divide migration in real landscapes, quantify the effects of orographic rainfall and erosion on the topographic evolution of a mountain range, and highlight the importance of considering drainage-divide migration when interpreting thermochronometer age patterns.},
  language  = {en}
}