Earth and Planetary Science Letters 506 (2019) 388–396

Contents lists available at ScienceDirect

Earth and Planetary Science Letters

www.elsevier.com/locate/epsl

Sedimentary loading–unloading cycles and faulting in intermontane

basins: Insights from numerical modeling and field observations
in the NW Argentine Andes
Paolo Ballato a,c,∗ , Sascha Brune b,c , Manfred R. Strecker c
a
Department of Science, University of Roma Tre, Largo San Leonardo Murialdo, 1, 00146, Rome, Italy
b
Helmholtz-Centre Potsdam, GFZ German Research Centre for Geosciences, Telegrafenberg, D-14473, Potsdam, Germany
c
Institute of Earth and Environmental Sciences, University of Potsdam, Karl-Liebknecht-Str. 24/25, D-14476 Potsdam, Germany

a r t i c l e i n f o a b s t r a c t

Article history: The removal, redistribution, and transient storage of sediments in tectonically active mountain belts
Received 19 September 2018 is thought to exert a first-order control on shallow crustal stresses, fault activity, and hence on the
Received in revised form 25 October 2018 spatiotemporal pattern of regional deformation processes. Accordingly, sediment loading and unloading
Accepted 28 October 2018
cycles in intermontane sedimentary basins may inhibit or promote intrabasinal faulting, respectively, but
Available online xxxx
unambiguous evidence for this potential link has been elusive so far. Here we combine 2D numerical
Editor: A. Yin
experiments that simulate contractional deformation in a broken-foreland setting (i.e., a foreland where
Keywords: shortening is diachronously absorbed by spatially disparate, reverse faults uplifting basement blocks)
sedimentary loading and unloading cycles with field data from intermontane basins in the NW Argentine Andes. Our modeling results suggest
intermontane basins that thicker sedimentary fills (>0.7–1.0 km) may suppress basinal faulting processes, while thinner fills
intrabasinal faulting (<0.7 km) tend to delay faulting. Conversely, the removal of sedimentary loads via fluvial incision and
Argentinean broken foreland basin excavation promotes renewed intrabasinal faulting. These results help to better understand the
2D numerical experiments
tectono-sedimentary history of intermontane basins that straddle the eastern border of the Andean
Andes
Plateau in northwestern Argentina. For example, the Santa María and the Humahuaca basins record
intrabasinal deformation during or after sediment unloading, while the Quebrada del Toro Basin
reflects the suppression of intrabasinal faulting due to loading by coarse conglomerates. We conclude
that sedimentary loading and unloading cycles may exert a fundamental control on spatiotemporal
deformation patterns in intermontane basins of tectonically active broken forelands.
© 2018 Elsevier B.V. All rights reserved.

1. Introduction lationships between shortening, sedimentary basin fills, and an

increase in loading stresses suppressing faulting, may be repre-
In tectonically active mountain belts the removal, redistribu- sented by the evolution of thick sedimentary fills in intermon-
tion or transient storage of sedimentary and ice loads may have tane basins and adjacent intraorogenic plateau basins (Sobel et al.,
a strong impact on the shallow crustal stress field and thus the 2003; Garcia-Castellanos, 2007; Strecker et al., 2009). This process
spatiotemporal pattern of earthquakes and deformation (Hetzel may be exacerbated by aridity. For example, the development of
and Hampel, 2005; Mey et al., 2016. This may apply to individ- orographic rainfall barriers along the flanks of an orogenic plateau
ual faults (e.g., Hampel and Hetzel, 2006; Hampel et al., 2007; of such intermontane basins may ultimately promote the aridifi-
Maniatis et al., 2009; Steer et al., 2014; Turpeinen et al., 2008) cation of the interior of the orogen, the defeat of fluvial systems,
and more complex, linked fault networks that constitute orogenic loss of fluvial connectivity with the foreland, and the establish-
wedges or extensional provinces (e.g., Davis et al., 1983; Friedrich ment of endorheic basins. Over time, erosion, shortening and pro-
et al., 2003; Hampel et al., 2009; Simpson, 2006; Strecker et al., tracted sedimentation in these basins will accumulate thick fills,
2007; Whipple, 2009; Calais et al., 2010; Vernant et al., 2013; which will increase lithostatic pressure affecting the basin cen-
Ballato et al., 2015). A prime example reflecting the dynamic re- ters and the fault-bounded basin margins; in turn this may result
in reduced, and finally aborted fault activity (Sobel et al., 2003;
Garcia-Castellanos, 2007; Strecker et al., 2009; Heidarzadeh et al.,
2017). Conversely, headward erosion toward the orogen interior,
*Corresponding author.
E-mail addresses: [email protected], [email protected] basin capture, and the re-establishment of fluvial connectivity with
(P. Ballato). the foreland will result in incision, basin exhumation, and a reduc-

https://doi.org/10.1016/j.epsl.2018.10.043
0012-821X/© 2018 Elsevier B.V. All rights reserved.
 P. Ballato et al. / Earth and Planetary Science Letters 506 (2019) 388–396 389

Fig. 1. A) Digital Elevation Model of NW Argentina showing the location of the studied Santa María, Quebrada del Toro and Humahuaca basins and associated drainage
networks. Inset shows location of Fig. A and major morphostructural provinces of the southern Central Andes: SAB – Subandean Belt; EC – Eastern Cordillera; SBS – Santa
Barbara System; SP – Sierras Pampeanas; AP – Altiplano-Puna Plateau; FPC – Frontal and Precordillera; and PC – Principal Cordillera. B) Simplified tectono-stratigraphic
evolution of studied basins. Note that all basins show characteristic patterns of sedimentary filling, incision, and deformation. C) Regional-scale cross-section showing the
distribution of fault-bounded mountain ranges and intermontane basins across the eastern Andean flank at ca. 24.5◦ S (see Fig. A for location, modified from Pearson et al.,
2013). Note that the Calchaquí Valley is the northern continuation of the Santa María Valley.

tion in loading stresses, which should promote fault reactivation. verse faults and constitute an integral part of the evolution of the
Although a similar concept has been proposed and numerically in- Eastern Cordillera and the broken foreland of the Andes. Collec-
vestigated to explain Quaternary slip-rate variations during glacial– tively, our study suggests that surface processes associated with
interglacial cycles in an extensional setting (i.e., Hetzel and Ham- different degrees of transient sediment loading may modulate the
pel, 2005; Hampel et al., 2009), or to quantify stress transfer in spatiotemporal pattern of faulting in the interior of tectonically ac-
eroding landscapes during the seismic cycle (i.e., Calais et al., 2010; tive mountain belts.
Steer et al., 2014), the importance of oscillating sediment loading
and unloading processes for controlling fault activity in intermon- 2. Modeling setup
tane, compressional basins has neither been fully recognized yet
in field studies nor has it been explored by means of numerical Our model approach aims at employing a low-complexity setup
studies. that nevertheless captures the key tectonic processes of intermon-
Here, we employ 2D numerical modeling experiments to inves- tane basins. The model simulates an area that is 200 km wide, ini-
tigate faulting within the upper continental crust during sedimen- tially 15 km thick, and consists of a single brittle crustal layer over-
tary loading–unloading cycles in intermontane basins. We discuss lying a 1-km-thick décollement horizon at the base of the model
our results in light of field observations from three tectonically ac- (Fig. 2A). The model is horizontally shortened by adding new ma-
tive sedimentary basins in the northwestern Argentine Andes along terial at a velocity of 4 mm/yr from the right-hand side. The setup
the eastern margin of the Andean Plateau (Puna), where several features a free surface at the top, and free-slip boundaries allowing
cycles of basin filling and excavation have occurred during the last for tangential motion at the bottom and lateral boundaries.
few Ma (e.g., Strecker et al., 2009 and references therein; Pingel We employ a 2D version of the elasto-visco-plastic finite-
et al., 2013; Fig. 1). These basins are bounded by high-angle re- element code SLIM3D (Popov and Sobolev, 2008) that solves the
390 P. Ballato et al. / Earth and Planetary Science Letters 506 (2019) 388–396

Fig. 2. 2D modeling results showing strain-rate distributions A) after 2.0 and 4.5 Myr relevant for all models. Active faults are indicated by yellow to red colors. B) and C)
Model results after 6.0, 7.5, and 9.0 Myr from model M1 without sedimentary load and model M2 with a 1.25-km-thick sedimentary load applied for a duration of 4 m.y.
For M1 the sampling point (black circle) starts to be uplifted after ca. 6 Myr, while in M2 uplift occurs during/after unloading only. (For interpretation of the colors in the
figure(s), the reader is referred to the web version of this article.)

conservation equations of momentum, mass, and energy by em- unloading. Instead of implementing mass conservative surface pro-
ploying an Arbitrary Lagrangian–Eulerian formulation with the cesses (e.g., Burov and Poliakov, 2001; Armitage et al., 2014), we
particle-in-cell technique (Gerya, 2010). This code has been exten- therefore approximate the impact of sediments simply as vertical
sively benchmarked and previously applied to a variety of conver- loads at the top boundary (Clift et al., 2015) that fill a basin to
gent (Quinteros and Sobolev, 2012; Fennell et al., 2018), divergent a specified horizontal base level (Figs. 2C and 3). This base level
(Brune et al., 2012, 2016, 2017a, 2017b), and transform (Popov et is defined in a time-dependent way and linearly increases from 0
al., 2012; Brune, 2014) plate-boundary settings. Due to the simplic- to a maximum value within 1 m.y.; the model stays at this maxi-
ity of the setup with a single brittle layer, the only relevant defor- mum sediment load for a given time, and finally resets to zero load
mation mechanism is parameterized through the Mohr–Coulomb (Fig. 4A). To convert sediment load to sediment thickness, we use
failure criterion. a compacted sediment density of 2300 kg/m3 . This setup enables
Our model focuses on the loading and unloading cycles in inter- us to investigate the response of brittle compressional tectonics to
montane basins and does not include inherited structures formed sediment loading–unloading cycles in a self-consistent modeling
during previous phases of compression or extension. However, we framework.
account for tectonic inheritance in a simplified way by employing
an initial friction coefficient of 0.2 that is considerably less than 3. Model results
for undisturbed, intra-continental crust (Zoback et al., 1993). For
the cohesion parameter we use a value of 5 MPa. Long-lived faults The reference model M1 without sedimentation shows initial
are generated self-consistently via the standard technique of a homogeneous layer-parallel shortening until ∼2 Myr (model time)
strain-dependent friction coefficient (e.g. Huismans and Beaumont, when two inward-dipping faults with a dip of ∼45◦ nucleate and
2003; Buiter et al., 2006), where the friction coefficient at each el- accommodate shortening (pop-up structure; Fig. 2; Supplementary
ement decreases linearly from 0.2 to 0.02 between local strains of Material SM1). The dip of these faults increases during shorten-
0 and 0.5. We seed a small component of random strain as an ad- ing to ∼60◦ and becomes shallower within the first few km of
ditional means to account for simplified inherited structures and depth when topography starts developing (Thieulot et al., 2014).
in order to break the initial symmetry of the homogeneous layer Subsequently, a new system of inward-dipping faults nucleates
(Brune et al., 2014). and accommodates shortening. Each pop-up structure accommo-
The main objective of this modeling study is to investigate the dates shortening for about 2 million years and then becomes in-
impact of sedimentary loading and unloading cycles on deeper active. Note that the associated uplift rates (Fig. 4) are slightly
crustal tectonics. This setup applies to the basin configurations of overestimated in our study since isostatic compensation has to
NW Argentina where mass redistribution is highly non-cylindrical. be expected in this area where the elastic lithospheric thickness
Sediment is hence transported out-of-plane of the 2D model, is <20 km (Tassara et al., 2007; Pérez-Gussinyé et al., 2008). The
which is why mass is not conserved during sediment loading and general pattern of uplift and subsidence, however, is independent
 P. Ballato et al. / Earth and Planetary Science Letters 506 (2019) 388–396 391

Fig. 4. Time-dependent sediment load and tectonic response for all models.
A) Loading–unloading functions prescribed to the models. B) Uplift rates at the cen-
ter of the modeled basins (sampling point in Figs. 2 and 3). Dashed lines show
discrete time steps shown in Figs. 2 and 3.

Models M2 and M3 are characterized by a cycle of sedimen-

tary loading and unloading starting from 5 Myr (Figs. 2C and 3A;
Supplementary Animations SM2 and SM3). Loading and unloading
phases are separated by a period of constant load, equivalent to
a 1.25-km-thick sedimentary cover. These cycles last 4 and 2 m.y.
for model M2 and M3, respectively. Sediments are accumulated in
a topographic depression 40 to 50 km long, and delimited by two
adjacent, fault-bounded mountain ranges, very similar to reverse-
fault bounded intermontane basins in NW Argentina (Fig. 1). When
sediment unloading starts, a complex system of inward-dipping
faults nucleates within the basin until a major pop-up structure
generates basin uplift. This starts at about 0.2 and 0.4 m.y. after the
onset of sediment unloading, which corresponds to the removal of
0.5 and 1 km of sediments for model M2 and M3, respectively
(Fig. 4).
Models M4 and M5 are characterized by 0.75 and 0.5 km of
sediment loading, respectively (Figs. 3B and 3C; Supplementary
Animations SM4 and SM5). Although the loading–unloading func-
tion is like that of model M2, albeit at lower amplitude, basin up-
lift for both models starts at 7 Myr, and thus 1.4 m.y. earlier than
Fig. 3. 2D modeling results showing strain-rate distributions for models M3 to M5
the onset of basin unloading. This generates basin uplift and the
after 6.0, 7.5, and 9.0 Myr. Active faults are indicated by yellow to red colors.
A) Model M3 with a 1.25-km-thick sedimentary load for 2 m.y. B) Model M4 with a progressive reduction of sediment load until the basin is entirely
0.75-km-thick sedimentary load for 4 m.y. and C) Model M5 with a 0.5-km-thick evacuated at 8.8 and 7.5 Myr, for models M4 and M5, respectively.
sedimentary load for 4 m.y. While model M3 (and M2, see Fig. 2) exhibits de- Once the sediment load has been removed, basin uplift occurs via
formation and uplift during/after sedimentary unloading, M4 and M5 reveal fault a new set of inward-dipping faults, and deformation shifts toward
activation before the prescribed unloading event.
less deformed areas (Fig. 4).
Our models document that sediment loading and unloading cy-
of isostasy, as long as individual mountain ranges are capable of cles can have a strong impact on fault activity in intermontane
trapping large amounts of sediments within intermontane basins. basins, and consequently on the spatiotemporal pattern of defor-
At ∼9 Myr, the entire modeling space is occupied by mountain mation in an orogen. The comparison between the reference model
ranges with a wavelength of 30 to 40 km and elevations of 2–3 km M1 and models M2 and M3, shows that 1.25 km of sediment load
(Fig. 2B). At this stage, all faults are maximally weakened so that suppress the nucleation of new faults, promoting deformation in
the pop-up faults become simultaneously active and accommodate undeformed areas. Conversely, the onset of sediment unloading is
shortening until model termination at 12 Myr. This setup repro- directly linked to subsequent fault nucleation and basin uplift in
duces the fundamental structures of the transition between the both models (Fig. 4). The thickness of the sediment load, however,
Eastern Cordillera and the broken foreland of NW Argentina, and seems to play a key role in hampering or promoting fault activity.
specifically of the Santa María, Humahuaca, and Quebrada del Toro In fact, for a sediment thickness of 0.75 and 0.5 km (models M4
basins, which are delimited by a system of basement cored, mostly and M5, respectively), basin uplift and exhumation occurs before
double-verging, high-angle reverse faults (Fig. 1). the onset of sediment unloading. This implies that such a load is
392 P. Ballato et al. / Earth and Planetary Science Letters 506 (2019) 388–396

Fig. 5. A), B) Satellite images from the Santa María Valley (Fig. 1) with major intrabasinal and bounding structures. Note large Quaternary pediment surfaces extending from
basin-bounding ranges on the east side. C) Southeast-facing photograph near Entre Ríos (south of Fig. B). Black arrows mark the highest regional pediment surface (ca. 1 km
above the local base level) formed in deformed Late Pliocene Punaschotter conglomerates. This surface is in thrust contact with the basement rocks of the Sierra Aconquija.
D) Geological profiles showing four generations of locally deformed pediment surfaces. E) Panoramic view to the south (location in Fig. B) highlighting a major reverse fault
juxtaposing basement rocks with 1.2-Myr-old pediment surface II. Note that the same pediment surface, farther to the west, is back tilted (see also Fig. D).

not sufficient for delaying fault nucleation as observed for models 4.1. Santa María Basin
M2 and M3 (Fig. 4).
In the Santa María Basin, intermontane sedimentation in associ-
ation with localized intrabasinal deformation lasted until ∼2.9 Ma
4. Observations from intermontane basins in Northwestern with the deposition of the syntectonic Corral Quemado Forma-
Argentina tion (Fig. 1B; Santa María Group; Bossi et al., 2001; Kleinert and
Strecker, 2001). This was followed by basin deformation and ero-
The arid Andean Plateau comprises reverse-fault bounded sion during intrabasinal west-verging thrusting (Figs. 1B and 5;
mountain ranges and sedimentary basins that record protracted Strecker et al., 1989). Subsequently, an episode of basin filling
endorheic conditions in association with ∼4 km of late Cenozoic, led to the accumulation of an at least ∼800-m-thick sequence of
basement-derived conglomerates sourced from the uplifting ranges
large-scale surface uplift (e.g., Alonso et al., 1991; Garzione et al.,
between 2.9 and 2.5 Ma (Figs. 1B and 5C; Punaschotter conglom-
2006; Rohrmann et al., 2016; Fig. 1). The southeastern margin
erates; Strecker et al., 1989). These gently deformed conglomerates
of the plateau, the Argentine Puna, includes the Cenozoic East-
are preserved only along the southernmost sector of the basin,
ern Cordillera, which transitions into the broken foreland of the
where they have been cut by a major west-verging fault uplift-
Sierra Pampeanas to the south and the Santa Barbara System of in-
ing the 5-km-high Sierra Aconquija. Their occurrence at ∼800 m
verted Cretaceous extensional basins to the east, respectively (Mon
above the present-day base level suggests that the entire basin
and Salfity, 1995). This region is characterized by a series of in-
and most likely its southern and northern continuation (El Ca-
termontane basins that repeatedly experienced transient episodes
jon and Calchaqui basins, respectively; Mortimer et al., 2007;
of basin filling, incision, and a re-establishment of fluvial connec- Bywater-Reyes et al., 2010; Pingel et al., 2016) were rapidly filled
tivity with the foreland (Fig. 1; Tchilinguirian and Pereyra, 2001; by a significant volume of sediments (Fig. 1). Similar events oc-
Sobel et al., 2003; Strecker et al., 2009; Hain et al., 2011). The curred diachronously between 6 and 3 Ma at least in five inter-
Santa María, Humahuaca, and Quebrada del Toro basins are lo- montane basins of this sector of the Eastern Cordillera and the
cated in this region of compressional basins and ranges. These Sierras Pampeanas (Penck, 1920); these conglomerates were inter-
basins were an integral part of a contiguous foreland until ∼5.2, preted to reflect enhanced tectonic deformation and increased sea-
4.2 and 6.0 Ma, respectively, when spatially and temporally dis- sonal precipitation associated with the growth of orographic bar-
parate rock uplift led to foreland fragmentation and intermontane riers (Schoenbohm et al., 2015 and references therein). Although
basin development (Bossi et al., 2001; Kleinert and Strecker, 2001; the causes for this episode of aggradation are not the focus of our
Marrett and Strecker, 2000; Strecker et al., 2009; Pingel et al., study, basin-wide deposition of thick conglomerates in an inter-
2013; Pearson et al., 2013). Below, we discuss the basin fill and montane basin suggests reduced fluvial connectivity with the fore-
faulting history of these intermontane basins in the context of our land, associated either with damming of the basin outlet triggered
modeling results. by landslides (but it generally lasts only a few thousand years; e.g.,
 P. Ballato et al. / Earth and Planetary Science Letters 506 (2019) 388–396 393

Fig. 6. A) Perspective satellite view of the Tilcara area in the Humahuaca Basin (Fig. 1). Note the occurrence of a complex basin-fill history where different generations of
Quaternary river terraces and conglomerates cover unconformably faulted Cretaceous to Paleogene deposits of the Salta Group and Late Miocene to Pliocene intermontane
deposits. B) Field photo (see Fig. A) of a thrust fault juxtaposing the Maimará Formation with Late Quaternary terrace gravels; this setting suggests horizontal and vertical
displacement rates of approximately 0.6 and 0.3 mm/yr over the last 63 ka, respectively (Robinson et al., 2005; Sancho et al., 2008; Pingel et al., 2013).

Hermanns et al., 2004), tectonic uplift and severed fluvial networks and erosion through the development of a system of east-verging
(e.g., Hilley and Strecker, 2005), or reduced sediment-transport intrabasinal thrust faults (Fig. 1B; Pingel et al., 2013). Within
capacity through to a combination of increased sediment supply the last 1 Ma, however, the basin was characterized by a major
and decreased river discharge (e.g., Ballato and Strecker, 2014; episode of filling with the deposition of conglomerates at least un-
Schildgen et al., 2016; Tofelde et al., 2017). In any case, basin fill- til ∼800 ka, (Pingel et al., 2013; Streit et al., 2015; Fig. 6A) and
ing was followed by regional incision through the establishment few short episodes of filling and incision (at least 3 events over
of a longitudinal fluvial network (Rio Santa María) that led to a the last 110 ka; Figs. 1B and 6; Tchilinguirian and Pereyra, 2001;
significant removal of the Punaschotter conglomerates and the de- Robinson et al., 2005; Schildgen et al., 2016). The widespread oc-
velopment of five generations of pediments at successively lower currence of Pleistocene deposits at ∼350 m above the current
elevations adjusted to the incising Río Santa María (Figs. 1B and 5; base level indicates reduced fluvial connectivity with temporary
Strecker et al., 1989). These multiple surfaces document alternat- basin-wide deposition of thick conglomerate fills that were sub-
ing sediment transport processes associated with an alluvial-fan sequently removed by fluvial incision. Importantly, these cycles
setting during transitory stable base levels of the Río Santa María were associated with renewed intrabasinal deformation as doc-
and subsequent incision that caused the termination of these areas ument by regional unconformities, folded and faulted terraces,
as transport surfaces. As the alluvial-fan gravels of these surfaces and cross-cutting relationships (Sancho et al., 2008; Pingel et al.,
are generally thin, laterally extensive, and cover the underlying 2013). Near the town of Tilcara, an up to 200-m-thick sequence
Teriary sedimentary units unconformably, these extensive surfaces of 89 to 63-ka-old terrace conglomerates was overthrusted by late
represent pediments. Pediment surfaces in the piedmont region of Miocene to Pliocene deposits of the Maimará Formation at hor-
the Aconquija range are in thrust contact with the adjacent base- izontal and vertical displacement rates of ∼0.6 and 0.3 mm/yr
ment, and the original westward-sloping surfaces were affected over the last 63 ka, respectively (Fig. 6B; Sancho et al., 2008;
by folding and thrusting, in places leading to a reversal of ped- Pingel et al., 2013). These observations suggest that renewed fault-
iment inclination (pediment surfaces I, II and III; Figs. 5D and ing occurred during or after unloading of the late Pleistocene basin
5E) between ∼2.5 and 0.6 Ma (Fig. 1B). In some cases, existing fill.
folds were re-folded during the Quaternary incision phase. Overall,
these observations document that the vestiges of the pediments 4.3. Quebrada del Toro Basin
are valuable strain markers that document renewed tectonic defor-
mation of piedmont areas during the downwearing and removal of In the Quebrada del Toro Basin, intermontane sedimentation
the conglomeratic basin fills. was active between ∼4.17 and 0.98 Ma (syntectonic Sola, Alfarcito
and Carachi conglomerates; Fig. 1B; Marrett and Strecker, 2000;
4.2. Humahuaca Basin Hilley and Strecker, 2005). Subsequently, the Plio-Pleistocene con-
glomerates were intensely deformed, eroded, and again covered
In the Humahuaca Basin, intermontane sedimentation lasted by at least 700-m-thick, coarse basement-derived conglomerates
until ∼2.5 Ma (Tilcara Formation) and 1.5 Ma for the northern and intercalated lake sediments (Figs. 1B and 7) most likely dur-
sectors of the basin (Uquia Formation; Pingel et al., 2013 and ing aggradation over a fault-bounded uplifted block that led to
references therein; Fig. 1B). This was followed by basin uplift a progressive increase in the local base level and steepening of
394 P. Ballato et al. / Earth and Planetary Science Letters 506 (2019) 388–396

Fig. 7. Perspective satellite view of the San Bernardo area in the Quebrada del Toro Basin (see Fig. 1) showing deformed Miocene to Pleistocene clastic and volcaniclastic
deposits covered in angular unconformity (see dashed line) by 1- to 0.5-Myr-old terrace forming conglomerates. B) Field photo taken few km S of Fig. 7C showing fluvial
terraces T6 to T2. C) Field photo (see Fig. 7A for location) showing the geometrical relationship between the Lower Alfarcito conglomerates and the ∼1- to 0.5-Myr-old
conglomerates in the footwall of the San Bernardo Fault.

the downstream channel (Hilley and Strecker, 2005). In turn, these al., 2003; Garcia-Castellanos, 2007). However, if such a basin is
units were later incised episodically, forming a staircase morphol- overfilled, it will be incised and exhumed (e.g., Heidarzadeh et
ogy of successively lower terraces (Marrett and Strecker, 2000; al., 2017; Garcia-Castellanos and O’Connor, 2018). Alternatively, if
Hilley and Strecker, 2005; Tofelde et al., 2017), with at least six tectonic rates necessary to sustain basin isolation are lower than
cut-and-fill episodes during the last 500 ka (Figs. 1B and 7; Tofelde incision rates associated with headward erosion and knick-point
et al., 2017). Except for the northwestern sectors of the basin, migration that might be enhanced during times of climate change
where <10◦ tilted conglomerates exist in the footwall of the basin- and changes in runoff, basin incision and exhumation will follow
bounding San Bernardo Fault, the conglomeratic basin fill appears (Hilley and Strecker, 2005). In either case, the removal of sedi-
to be rather undeformed (Fig. 7C). In contrast to the other two in- ments may cause a re-activation or the new formation of contrac-
termontane basins, these observations suggest that loading stresses tional faults in the intermontane basin. This is well recorded in the
associated with the deposition of the second-generation of con- Santa Maria and Humahuaca basins where ∼800- and 350-m-thick
glomeratic basin fill may have caused intrabasinal suppression of sedimentary basin fills where largely removed by fluvial erosion.
faulting in the Quebrada del Toro. The ongoing incision and lim- Conversely, the limited removal of basin-fill sediments in the up-
ited removal of the valley fill (up to 220 m of net incision, but only per Quebrada del Toro Basin (up to 220 m, but only along the
along the Rio Toro and the San Bernardo Fault which is located Toro river which flows at the western margin of the basin), has
at the western edge of the basin; Figs. 1B and 7; Tofelde et al., not resulted yet in a decrease in lithostatic pressure followed by
2017), however, has not resulted yet in sufficient basin-wide sedi- renewed intrabasinal faulting.
ment unloading to have caused pronounced intrabasinal faulting.
6. Conclusions
5. Discussion
We illustrate that sediment loading–unloading cycles may have
Our results highlight the potential impact that combined a strong impact on the tectonic activity of intermontane basins.
climate-driven surface processes and tectonics may have on the Accordingly, sediment loading may suppress intrabasinal faulting,
level of fault activity in a tectonically active mountain belt. A con- while sediment unloading promotes faulting. This is more evident
ceptual sketch elucidating the relationships between sedimen- for thicker basin fills, while for thinner fills (<0.5 km over a length
tary loading–unloading and intrabasinal deformation based on our scale of 30–40 km), faulting appears to be delayed rather than
combined modeling and field study is shown in Fig. 8. As the inhibited. Our modeling results agree with previous studies that
rock-uplift rate increases relative to the transport efficiency and focused on the effect of glacial loading/unloading cycles on fault
incision capability of the principal stream draining a basin, in- activity. In contrast to these studies where pre-defined faults are
ternal drainage conditions, basin-wide aggradation, and locking of submitted to a duration of glacial loads on the order of thousands
the basin-bounding and intrabasinal faults may ensue, followed of years, we focus on sedimentary loading spanning several Myr
by a migration of tectonic activity to another area (e.g., Sobel et and investigate self-consistently generated crustal fault networks.
 P. Ballato et al. / Earth and Planetary Science Letters 506 (2019) 388–396 395

– Surface processes, tectonics and georesources: the Andean fore-

land basin of Argentina; STR 373/38-1). S.B. was supported by the
Helmholtz Young Investigators Group CRYSTALS (VH-NG-1132). We
thank John Armitage and Victor Ramos for their constructive and
motivating reviews. We are grateful to Heiko Pingel for helping
with figure drafting and to Carolina Montero Lopez, Fernando Hon-
gon, Leonardo Escalante and Alejandro Aramayo for introducing
the geology of Argentina. Finally, PB is indebted to Carolina Mon-
tero Lopez and her family for logistical support.

Appendix A. Supplementary material

Supplementary material related to this article can be found on-

line at https://doi.org/10.1016/j.epsl.2018.10.043.

References

Alonso, R., Jordan, T.E., Tabbutt, K.T., Vandervoort, D.S., 1991. Giant evaporite belts
of the Neogene central Andes. Geology 19, 401–404.
Armitage, J.J., Duller, R.A., Schmalholz, S.M., 2014. The influence of long-wavelength
tilting and climatic change on sediment accumulation. Lithosphere 6, 303–318.
https://doi.org/10.1130/L343.1.
Ballato, P., Landgraf, A., Fox, M., Stockli, D., Schildgen, T.F., Ghassemi, M.R., Kirby, E.,
Strecker, M.R., 2015. The growth of a mountain belt forced by base-level fall:
tectonics and surface processes during the evolution of the Alborz Mountains,
N Iran. Earth Planet. Sci. Lett. 425, 204–218.
Ballato, P., Strecker, M.R., 2014. Assessing tectonic and climatic causal mecha-
nisms in foreland basin stratal architecture: insights from the Alborz Mountains,
northern Iran. Earth Surf. Process. Landf. 39, 110–125.
Bossi, G.E., Georgieff, S.M., Gavriloff, I.J.C., Ibañez, L.M., Muruaga, C.M., 2001. Ceno-
zoic evolution of the intramontane Santa María basin, Pampean Ranges, north-
western Argentina. J. South Am. Earth Sci. 14, 725–734.
Brune, S., 2014. Evolution of stress and fault patterns in oblique rift systems: 3-D
numerical lithospheric-scale experiments from rift to breakup. Geochem. Geo-
phys. Geosyst. 15, 3392–3415.
Brune, S., Corti, G., Ranalli, G., 2017b. Controls of inherited lithospheric hetero-
geneity on rift linkage: numerical and analog models of interaction between
the Kenyan and Ethiopian rifts across the Turkana depression. Tectonics 36,
2017TC004739.
Brune, S., Heine, C., Clift, P.D., Pérez-Gussinyé, M., 2017a. Rifted margin architecture
and crustal rheology: reviewing Iberia-Newfoundland, Central South Atlantic,
and South China Sea. Mar. Pet. Geol. 79, 257–281.
Brune, S., Heine, C., Perez-Gussinye, M., Sobolev, S.V., 2014. Rift migration explains
continental margin asymmetry and crustal hyper-extension. Nat. Commun. 5.
Brune, S., Popov, A.A., Sobolev, S.V., 2012. Modeling suggests that oblique extension
facilitates rifting and continental break-up. J. Geophys. Res. 117, B08402.
Fig. 8. Conceptual model depicting the relationships between sedimentary loading Brune, S., Williams, S.E., Butterworth, N.P., Müller, R.D., 2016. Abrupt plate accelera-
and unloading cycles on deformation processes in a fault-bounded intermontane tions shape rifted continental margins. Nature 536, 201–204.
basin (see text for details). Buiter, S.J.H., Babeyko, A.Y., Ellis, S., Gerya, T.V., Kaus, B.J.P., Kellner, A., Schreurs, G.,
Yamada, Y., 2006. The numerical sandbox: comparison of model results for a
shortening and an extension experiment. In: Buiter, S.J.H., Schreurs, G. (Eds.),
The ensemble of previous studies and our present work leads to Analogue and Numerical Modelling of Crustal-Scale Processes. In: Geol. Soc.
the robust conclusion that a thickness of several 100 m of sedi- (Lond.) Spec. Publ., vol. 253, pp. 29–64.
ment or ice are sufficient to delay the brittle tectonic regime while Burov, E., Poliakov, A., 2001. Erosion and rheology controls on synrift and postrift
evolution: verifying old and new ideas using a fully coupled numerical
layers of ∼1 km thickness may be capable to fully inhibit faulting.
model. J. Geophys. Res., Solid Earth 106, 16461–16481. https://doi.org/10.1029/
Combined, our modeling and field observations highlight the 2001JB000433.
rich interplay between tectonics and erosional surface processes, Bywater-Reyes, S., Carrapa, B., Clementz, M., Schoenbohm, L.M., 2010. Effect of
which dictate the removal, redistribution and transient storage of late Cenozoic aridification on sedimentation in the Eastern Cordillera of north-
sedimentary loads, and thus modulate loading stresses and ac- west Argentina (Angastaco basin). Geology 38, 235–238. https://doi.org/10.1130/
G30532.1.
tive deformation patterns in the intermontane basis of the south-
Calais, E., Freed, A.M., Van Arsdale, R., Stein, S., 2010. Triggering of New Madrid
ern Central Andes. We suggest that this process may occur in seismicity by late-Pleistocene erosion. Nature 466, 608–611.
other tectonically active intermontane basins as well. Basins with Clift, P.D., Brune, S., Quinteros, J., 2015. Climate changes control offshore crustal
similar characteristics where these processes may have impacted structure at South China Sea continental margin. Earth Planet. Sci. Lett. 420,
66–72.
basin evolution include the Cenozoic Tien Shan and Qilian Shan
Davis, D., Suppe, J., Dahlen, F.A., 1983. Mechanics of fold-and-thrust belts and accre-
mountains of Central Asia; the North American Laramide broken- tionary wedges. J. Geophys. Res. 88, 1153–1172.
foreland province may represent a related paleotectonic setting Fennell, L.M., Quinteros, J., Iannelli, S.B., Litvak, V.D., Folguera, A., 2018. The role
during the Cretaceous and early Eocene. of the slab pull force in the late Oligocene to early Miocene extension in the
Southern Central Andes (27◦ –46◦ S): insights from numerical modeling. J. South
Am. Earth Sci. 87, 174–187.
Acknowledgements
Friedrich, A.M., Wernicke, B.P., Niemi, N.A., Bennett, R.A., Davis, J.L., 2003. Compar-
ison of geodetic and geologic data from the Wasatch region, Utah, and impli-
P.B. and M.S. were financed by the German Science Foundation’s cations for the spectral character of Earth deformation at periods of 10 to 10
(DFG) International Research Training Group IGK2018 (StRATEGy million years. J. Geophys. Res. 108, 2199.
396 P. Ballato et al. / Earth and Planetary Science Letters 506 (2019) 388–396

Garcia-Castellanos, D., 2007. The role of climate in high plateau formation. Insights Popov, A.A., Sobolev, S.V., 2008. SLIM3D: a tool for threedimensional thermo me-
from numerical experiments. Earth Planet. Sci. Lett. 257, 372–390. chanical modeling of lithospheric deformation with elasto-visco-plastic rheol-
Garcia-Castellanos, D., O’Connor, J.E., 2018. Outburst floods provide erodibility esti- ogy. Phys. Earth Planet. Inter. 171, 55–75.
mates consistent with long-term landscape evolution. Sci. Rep. 8, 10573. Popov, A.A., Sobolev, S.V., Zoback, M.D., 2012. Modeling evolution of the San Andreas
Garzione, C.N., Molnar, P., Libarkin, J.C., MacFadden, B.J., 2006. Rapid late Miocene Fault system in northern and central California. Geochem. Geophys. Geosyst. 13,
rise of the Bolivian Altiplano: evidence for removal of mantle lithosphere. Earth Q08016.
Planet. Sci. Lett. 241, 543–556. Quinteros, J., Sobolev, S.V., 2012. Why has the Nazca plate slowed since the Neo-
Gerya, T.V., 2010. Introduction to Numerical Geodynamic Modelling. Cambridge Uni- gene? Geology 41, 31–34.
versity Press. 358 pp. Robinson, R.A.I., Spencer, J.Q.G., Strecker, M.R., Richter, A., Alonso, R.N., 2005. Lumi-
Hain, M.P., Strecker, M.R., Bookhagen, B., Alonso, R.N., Pingel, H., Schmitt, A.K., 2011. nescence dating of alluvial fans in intramontane basins of NW Argentina. In:
Neogene to Quaternary broken foreland formation and sedimentation dynam- Harvey, A., Mather, A., Stokes, M. (Eds.), Alluvial Fans: Geomorphology, Sedi-
ics in the Andes of NW Argentina (25◦ S). Tectonics 30. https://doi.org/10.1029/ mentology, Dynamics. In: Geol. Soc. (Lond.) Spec. Publ., vol. 251, pp. 153–168.
2010TC002703.
Rohrmann, A., Sachse, D., Mulch, A., Pingel, H., Tofelde, S., Alonso, R.N., Strecker,
Hampel, A., Hetzel, R., 2006. Response of normal faults to glacial–interglacial fluctu-
M.R., 2016. Miocene orographic uplift forces rapid hydrological change in the
ations of ice and water masses on Earth’s surface. J. Geophys. Res. 111, B06406.
southern Central Andes. Sci. Rep. 6, 35678.
https://doi.org/10.1029/2005JB004124.
Sancho, C., Pena, J.L., Rivelli, F., Rhodes, E., Muñoz, A., 2008. Geomorphological evo-
Hampel, A., Hetzel, R., Densmore, A.L., 2007. Postglacial slip rate increase on the
lution of the Tilcara alluvial fan (Jujuy Province, NW Argentina): tectonic im-
Teton normal fault, northern Basin and Range Province, caused by melting
plications and paleoenvironmental considerations. J. South Am. Earth Sci. 26,
of the Yellowstone ice cap and deglaciation of the Teton Range? Geology 35,
68–77.
1107–1110.
Hampel, A., Hetzel, R., Maniatis, G., Karow, T., 2009. Three-dimensional numerical Schildgen, T., Robinson, R.A.J., Savi, S., Phillips, W.M., Spencer, J.Q.G., Bookhagen,
modeling of slip rate variations on normal and thrust fault arrays during ice B., Scherler, D., Tofelde, S., Alonso, R.N., Kubik, P.W., Binnie, S.A., Strecker,
cap growth and melting. J. Geophys. Res. 114, 8406–8420. M.R., 2016. Landscape response to late Pleistocene climate change in NW Ar-
Heidarzadeh, G., Ballato, P., Hassanzadeh, J., Ghassemi, M.R., Strecker, M.R., 2017. gentina: sediment flux modulated by basin geometry and connectivity. J. Geo-
Tectonic and climatic controls on Pliocene sedimentary basin evolution, over- phys. Res. 121, 392–414.
spilling and incision of the northwestern Iranian Plateau, Mianeh Basin. Earth Schoenbohm, L.M., Carrapa, B., McPherson, H.M., Pratt, J.R., Bywater-Reyes, S., Mor-
Planet. Sci. Lett. 469, 135–147. timer, E., 2015. Climate and tectonics along the southern margin of the Puna
Hermanns, R.L., Niedermann, S., Ivy-Ochs, S., Kubik, P.W., 2004. Rock avalanching Plateau, NW Argentina: origin of the late Cenozoic Punaschotter conglomerates.
into a landslide-dammed lake causing multiple dam failure in Las Conchas val- Geol. Soc. Am. Mem. 212, 251–260.
ley (NW Argentina) – evidence from surface exposure dating and stratigraphic Simpson, G.D.H., 2006. Modelling interactions between fold-thrust belt deformation,
analyses. Landslides 1, 113–122. foreland flexure and surface mass transport. Basin Res. 18, 125–143.
Hetzel, R., Hampel, A., 2005. Slip rate variations on normal faults during glacial– Sobel, E.R., Hilley, G.E., Strecker, M.R., 2003. Formation of internally drained contrac-
interglacial changes in surface loads. Nature 435, 81–84. tional basins by aridity-limited bedrock incision. J. Geophys. Res. 108, 2344.
Hilley, G.E., Strecker, M.R., 2005. Processes of oscillatory basin filling and excavation Steer, P., Simoes, M., Cattin, R., Shyu, J.B.H., 2014. Erosion influences the seismicity
in a tectonically active orogen: Quebrada del Toro Basin, NW Argentina. Geol. of active thrust faults. Nat. Commun. 5, 5564.
Soc. Am. Bull. 117, 887–901. Strecker, M.R., Alonso, R.N., Bookhagen, B., Carrapa, B., Coutand, I., Hain, M.P., Hilley,
Huismans, R.S., Beaumont, C., 2003. Symmetric and asymmetric lithospheric exten- G.E., Mortimer, E., Schoenbohm, L., Sobel, E.R., 2009. Does the topographic dis-
sion: relative effects of frictional-plastic and viscous strain softening. J. Geophys. tribution of the central Andean Puna Plateau result from climatic or geodynamic
Res., Solid Earth 108, 2496. processes? Geology 37, 643–646.
Kleinert, K., Strecker, M.R., 2001. Climate change in response to orographic barrier
Strecker, M.R., Alonso, R., Bookhagen, B., Carrapa, B., Hilley, G.E., Sobel, E.R., Trauth,
uplift: paleosol and stable isotope evidence from the late Neogene Santa María
M.H., 2007. Tectonics and climate of the Southern Central Andes. Annu. Rev.
basin, northwestern Argentina. Geol. Soc. Am. Bull. 113, 728–742.
Earth Planet. Sci. 35, 747–787.
Maniatis, G., Kurfeß, D., Hampel, A., Heidbach, O., 2009. Slip acceleration on normal
Strecker, M.R., Cerveny, P., Bloom, A.L., Malizia, D., 1989. Late Cenozoic tectonism
faults due to erosion and sedimentation — results from a new three-dimensional
and landscape development in the foreland of the Andes: Northern Sierras Pam-
numerical model coupling tectonics and landscape evolution. Earth Planet. Sci.
peanas, Argentina. Tectonics 8, 517–534.
Lett. 284, 570–582.
Marrett, R., Strecker, M.R., 2000. Response of intracontinental deformation in the Streit, R.L., Burbank, D.W., Strecker, M.R., Alonso, R.N., Cottle, J.M., Kylander-Clark,
central Andes to late Cenozoic reorganization of South American plate motions. A.R.C., 2015. Controls on intermontane basin filling, isolation and incision on the
Tectonics 19, 452–467. margin of the Puna Plateau, NW Argentina (∼23◦ S). Basin Res. 29 (S1), 131–155.
Mey, J., Scherler, D., Wickert, A.D., Egholm, D.L., Tesauro, M., Schildren, T.F., Strecker, Tassara, A., Swain, C., Hackney, R., Kirby, J., 2007. Elastic thickness structure of
M.R., 2016. Glacial isostatic uplift of the European Alps. Nat. Commun. 7, 13382. South America estimated using wavelets and satellite-derived gravity data. Earth
https://doi.org/10.1038/ncomms13382. Planet. Sci. Lett. 253, 17–36.
Mon, R., Salfity, J.A., 1995. Tectonic evolution of the Andes of northern Argentina. In: Tchilinguirian, P., Pereyra, F.X., 2001. Geomorfologıa del sector Salinas Grandes-
Tankard, A.J., Suárez Soruco, R., Welsink, H.J. (Eds.), Petroleum Basins of South Quebrada de Humahuaca, Provincia de Jujuy. Rev. Asoc. Geol. Argent. 56, 3–15.
America. Am. Assoc. Petrol. Geol., Tulsa, Oklahoma, pp. 269–283. Thieulot, C., Steer, P., Huismans, R.S., 2014. Three-dimensional numerical simula-
Mortimer, E., Carrapa, B., Coutand, I., Schoenbohm, L., Sobel, E., Sosa Gomez, J., tions of crustal systems undergoing orogeny and subjected to surface processes.
Strecker, M.R., 2007. Fragmentation of a foreland basin in response to out- Geochem. Geophys. Geosyst. 15, 4936–4957.
of-sequence basement uplifts and structural reactivation: El Cajón Campo del Tofelde, S., Schildgen, T.F., Savi, S., Pingel, H., Wickert, A.D., Bookhagen, B., Wittmann,
Arenal Basin, NW Argentina. Geol. Soc. Am. Bull. 119, 637–653. https://doi.org/ H., Alonso, R.N., Cottle, J., Strecker, M.R., 2017. 100 kyr fluvial cut-and-fill ter-
10.1130/B25884.1. race cycles since the Middle Pleistocene in the southern Central Andes, NW
Pearson, D.M., Kapp, P., DeCelles, P.G., Reiners, P.W., Gehrels, G.E., Ducea, M.N., Argentina. Earth Planet. Sci. Lett. 473, 141–153.
Pullen, A., 2013. Influence of pre-Andean crustal structure on Cenozoic thrust Turpeinen, H., Hampel, A., Karow, T., Maniatis, G., 2008. Effect of ice sheet growth
belt kinematics and shortening magnitude: Northwestern Argentina. Geo- and melting on the slip evolution of thrust faults. Earth Planet. Sci. Lett. 269,
sphere 9, 1766–1782. 230–241.
Penck, W., 1920. Der Südrand der Puna de Atacama (NW-Argentinien). Abh.
Vernant, P., Hivert, F., Chéry, J., Steer, P., Cattin, R., Rigo, A., 2013. Erosion-induced
Math.-Phys. Kl. Sächs. Akad. Wiss. 37. 420 pp.
isostatic rebound triggers extension in low convergent mountain ranges. Geol-
Pérez-Gussinyé, M., Lowry, A.R., Morgan, J.P., Tassara, A., 2008. Effective elastic thick-
ogy 41, 467–470.
ness variations along the Andean margin and their relationship to subduction
Whipple, K.X., 2009. The influence on climate on the tectonic evolution of mountain
geometry. Geochem. Geophys. Geosyst. 9, 1–21.
belts. Nat. Geosci. 2, 97–104.
Pingel, H., Mulch, A., Alonso, R.N., Cottle, J., Hynek, S.A., Poletti, J., et al., 2016. Sur-
face uplift and convective rainfall along the southern Central Andes (Angastaco Zoback, M.D., Apel, R., Baumgartner, J., Brudy, M., Emmermann, R., Engeser, B.,
Basin, NW Argentina). Earth Planet. Sci. Lett. 440, 33–42. Fuchs, K., Kessels, W., Rischmuller, H., Rummel, F., Vernik, L., 1993. Upper-crustal
Pingel, H., Strecker, M.R., Alonso, R.N., Schmitt, A.K., 2013. Neotectonic basin and strength inferred from stress measurements to 6 km depth in the KTB borehole.
landscape evolution in the Eastern Cordillera of NW Argentina, Humahuaca Nature 365, 633–635.
Basin (∼24◦ S). Basin Res. 25, 554–573.