Exhumation history and landscape evolution of the Sierra de San Luis ( Sierras Pampeanas , Argentina )-new insights from low-temperature thermochronological data

This paper presents low-temperature thermochronological data and K-Ar fault gouge ages from the Sierra de San Luis in the Eastern Sierras Pampeanas in order to constrain its low-temperature thermal evolution and exhumation history. Thermal modelling based on (U-Th)/He dating of apatite and zircon and apatite fission track dating point to the Middle Permian and the Triassic/Early Jurassic as main cooling/exhumation phases, equivalent to ca. 40-50% of the total exhumation recorded by the applied methods. Cooling rates are generally low to moderate, varying between 2-10 °C/Ma during the Permian and Triassic periods and 0.5-1.5 °C/Ma in post-Triassic times. Slow cooling and, thus, persistent residence of samples in partial retention/partial annealing temperature conditions strongly influenced obtained ages. Thermochronological data indicate no significant exhumation after Cretaceous times, suggesting that sampled rocks were already at or near surface by the Cretaceous or even before. As consequence, Cenozoic cooling rates are low, generally between 0.2-0.5 °C/Ma which is, depending on geothermal gradient used for calculation, equivalent to a total Cenozoic exhumation of 0.6-1.8 km. K-Ar fault gouge data reveal long-term brittle fault activity. Fault gouge ages constrain the end of ductile and onset of brittle deformation in the Sierra de San Luis to the Late Carboniferous/Early Permian. Youngest K-Ar illite ages of 222-172 Ma are interpreted to represent the last illite formation event, although fault activity is recorded up to the Holocene.


Introduction
The Sierra de San Luis is one of the southernmost ranges of the Sierras Pampeanas region (Pampean ranges, Fig. 1).These ranges correspond to basement block uplifts surrounded by basins of flat topography, which widely crop out in centralwestern Argentina between 27° to 33°30ʼ S (Ramos et al., 2002 and references therein).This region is also regarded as an example of active thick-skinned crustal deformation related to flattening of the Nazca plate subduction since ~15 Ma (Jordan et al., 1983;Jordan and Allmendinger, 1986;Kay and Abbruzzi, 1996;Ramos et al., 2002).Although the Sierras Pampeanas show a different and longer evolutionary paths than the Andean orogen, they share a similar morphogenetic history after the flattening of the Nazca plate and are regarded as a morphotectonic component of the Andean building (Jordan and Allmendinger, 1986).2c and location of regional map depicted in figure 3; C. Cross-section of studied transect, sample location and location of inferred major fault zones (location of fault zones based on geological maps from San Luis and San Francisco del Monte de Oro (Costa et al., 2000(Costa et al., , 2001a)); note that major fault zones illustrated here does not have a clear, single surface expression).CS: crystalline basement; SC: sedimentary cover; dotted lines: paleosurfaces.
Basement blocks consist of late Precambrian to early Paleozoic metamorphic and igneous rocks showing a topographic asymmetry represented by a steep western and a gentle eastern slope (Fig. 2), being the latter usually characterized by erosional paleosurfaces remnants.This topographic asymmetry is interpreted to be linked to Neogene uplift along reverse faults with listric geometry (González Bonorino, 1950;Jordan and Allmendinger, 1986; among others), which are usually located at the steeper western sides of the ranges (Fig. 1).The Neogene uplift resulted from the incorporation of the Juan Férnandez Ridge into the subduction of the Nazca plate, causing southeastwardpropagating flat-slab subduction beneath the Pampean region (Yáñez et al., 2001) and low heat flow into the crust of the overriding plate as well (Dávila and Carter, 2013;Collo et al., 2015).
This investigation provides thermochronological data of the western part of the Sierra de San Luis, which are complemented with K-Ar illite fine-fraction dating on fault gouges.Based on this, our work aims to constrain the cooling and exhumation history of this part of the range since the Late Paleozoic.In addition, this work aims to check the possible morphotectonic evolution scenarios, especially in the context of the onset of brittle deformation and the formation of paleolandsurfaces.

Geological and Morphotectonic Setting
Geological studies of the Sierra de San Luis have traditionally been focused on the tectonometamorphic evolution of the crystalline basement, which is considered to be completed by Early Carboniferous times (Criado Roque et al., 1981;Ortiz Suárez et al., 1992;Von Gosen, 1998;Steenken et al., 2004).Deformation of the basement was dominated by brittle deformation ever since, leading to characteristic N-S trending fault blocks bounded by major reverse faults during Cenozoic crustal shortening (Fig. 1).
In the late Paleozoic, the San Luis range was part of the continental Paganzo basin, whose stratigraphic record is widespread along the Sierras Pampeanas (Salfity and Gorustovic, 1983;Mpodozis and Ramos, 1989;Ramos et al., 2002;Limarino and Spalletti, 2006;Ramos, 2009).In the study region, the record of Paganzo strata is scarce, although minor outcrops are present in the Bajo de Véliz depression in the northeastern Sierra de San Luis (Hünicken et al., 1981) (Fig. 2a).
Subsequent Mesozoic rifting along reactivated Paleozoic structures led to the development of continental basins and deposition of Mesozoic clastic sediments to the west and around the Sierra de San Luis (Criado Roque, 1972;Schmidt et al., 1995;Costa et al., 2000).Rifting was accompanied by alkaline intraplate volcanism, which is recorded around the Sierra de San Luis (Llambías and Brogioni, 1981).Additionally, López and Solá (1981) described isolated outcrops of rift-related volcanic rocks in the Sierra de San Luis with K-Ar ages of ca.83±5.85Ma.
Main Andean Cenozoic crustal shortening and flat-slab subduction of the Nazca Plate is inferred to be related to plate reorganization and collision of the Juan Fernández Ridge during the Miocene (Jordan et al., 1983;Yáñez et al., 2001;Ramos et al., 2002), giving rise to uplift of Pampean basement blocks including the Sierra de San Luis.These uplifts were usually controlled by former main bounding faults of the Mesozoic rift system and other relevant anisotropies of the basement, which were tectonically inverted as a broken foreland (Jordan et al., 1983;Jordan and Allmendinger, 1986;Schmidt et al., 1995;Costa, 1992;Ramos et al., 2002;).
Erosional surface remnants were envisaged as key markers to constrain the geometry and characteristics of the post-Paleozoic structural relief and the associated uplift and exhumation path, as no pre-Quaternary sedimentary cover was developed or preserved atop the ranges (Costa, 1992;Costa et al., 1999).These paleolandsurfaces have been regarded in two ways . González Díaz (1981) and Criado Roque et al. (1981) considered them as a formerly continuous and essentially synchronous surface, which was uplifted and disrupted into several minor surfaces and juxtaposed by faults during the Andean Orogeny.Alternatively, Carignano (1999) and Rabassa et al. (2010Rabassa et al. ( , 2014) ) suggested that erosional paleosurfaces represent diachronous planation episodes, which are confined to different topographic levels and separated by topographic scarps.Based on field observations, the latter authors propose that paleosurface ages should range between Late Paleozoic and Paleogene.Based on few thermochronological data, a similar proposal on diachronous regional paleosurfaces had been previously raised by Jordan et al. (1989).
Neogene faulting activity documented along the western steep scarp of the Sierra de San Luis suggests that this hillslope constitutes the neotectonic uplift of the range, where at least 1000 m of structural relief was built-up since planation surfaces were formed (Costa, 1992;Costa et al., 2001b).On the other hand, recent thermochronological studies in other regions of the Sierras Pampeanas near the Sierra de San Luis indicate that exhumation and cooling associated with uplift since the Cenozoic was considerably small (Löbens et al., 2011;Bense et al., 2013;Enkelmann et al., 2014).

Thermochronology
"Exhumation" and "uplift" are widespread terms used in the literature to refer to the vertical transport of rock masses.Nevertheless, misinterpretations may arise as it is frequently not clear what these terms refer to.England and Molnar (1990) defined exhumation as the displacement of rocks with respect to the surface.If a crustal thermal profile is assumed, exhumation rates can thus be derived from thermochronological data (England and Molnar, 1990;Stüwe and Barr, 1998;Ring et al., 1999).Uplift, in turn, is related to vertical movements with respect to the geoid, although the mean sea level can be considered as the reference level as well (England and Molnar, 1990).England and Molnar (1990) further differentiated between "surface uplift" and "uplift of rocks", depending on whether displacements of the surface or displacements of rocks are considered.In this work, exhumation and uplift are considered following these definitions, being the latter related to surface uplift sensu England and Molnar (1990).
(U-Th)/He ages from apatite (AHe) and zircon (ZHe) as well as apatite fission track ages (AFT) can be interpreted in terms of an exhumation-induced cooling through low temperature conditions and provide an important tool for quantifying the cooling of rocks as they pass through relatively shallow crustal levels.According to Donelick et al. (1999) and Ketcham et al. (1999), among others, the thermal sensitivity of the apatite fission track method, namely the partialannealing zone (PAZ; Gleadow and Fitzgerald, 1987), ranges between 130 °C and 60 °C.For the (U-Th)/He system of apatite and zircon, this temperature interval is referred to as the partial-retention zone (PRZ) and ranges between 65 °C and 30 °C and 185 °C and 135 °C, respectively (e.g., Baldwin and Lister, 1998;Wolf et al., 1998;Reiners and Brandon, 2006).
The analytical procedure on (U-Th)/He and apatite fission track dating is described in the Appendix.Based on dating results, a two-stage approach of forward and inverse modelling of thermochronological data using HeFTy software (Ketcham, 2005) was followed to model numerically possible t-T paths for individual samples.Especially the combination of fission track data (age and track length distribution) with (U-Th)/He data can provide a diagenetic and sensitive tool for evaluating low-temperature thermal history.
Two boundary conditions were set to the thermal modeling: 1. the beginning of the time-temperature path was constrained by the zircon (U-Th)/He data and 2. the end of the time-temperature paths was set to 17 °C, according to annual mean temperatures in the study area (Müller, 1996).

Fault gouge dating and interpretation
Under brittle conditions, tectonic slip causes the crushing of rocks and grain-size reduction along faults.In these localized fault zones, the increased area/volume ratio of rock fragments together with fluid circulation favors high chemical reactivity, allowing retrograde processes to produce fault gouges composed of authigenic hydrosilicates such as illite.Thus, the formation time of the authigenic illite in a fault gouge, can be correlated with periods of motion along the fault and thus constrains the timing of faulting where favorable conditions for illite formation are present (e.g., Lyons and Snellenberg, 1971;Kralik et al., 1987;Wemmer, 1991;Solum et al., 2005;Haines et al., 2008;Zwingmann et al., 2010;Surace et al., 2011;Wolff et al., 2012;Bense et al., 2014).Bense et al. (2014) suggested a concept to evaluate the timing of brittle deformation based on K-Ar illite fine-fraction ages from fault gouges, which are developed in non-sedimentary host rocks during retrograde cooling.Following this idea, the interpretation of K-Ar illite ages is constrained by evaluation of several independent parameters, e.g., illite crystallinity, illite polytype quantification, illite grain-size, clay mineralogical observations, K-Ar muscovite and biotite host-rock cooling ages, as well as low-temperature thermochronological data derived from, e.g., fission track or the (U-Th)/He dating.This allows a better error evaluation of individual methods by highlighting concordant data in multi-methodological datasets, as well as an easier combination of all observations into a consistent regional evolution model.For a detailed discussion of methodological, interpretation and analytical issues, the reader is referred to Bense et al. (2014).Further details on data interpretation and analytical procedures are given in the Appendix.

Thermochronology
Eight samples, mainly from granitic rocks of the crystalline basement were collected along a FIG. 3. Simplified DEM-based map of the study area with locations of thermochronological samples (green triangles) and K-Ar illite fault-gouge samples (red triangles).For location of map extent see figure 2.
transect encompassing the western slope and upper parts at the South of the Sierra de San Luis (Figs. 2, 3, Tables 4 and 5).Mean zircon (U-Th)/He ages range from the Late Carboniferous to Middle Triassic (313 to 229 Ma).Ages show no clear correlation with elevation (Fig. 7, Table 5) and, due to the high scatter of individual ages, single ZHe ages are not interpreted to represent distinct events.Instead, ZHe ages are used as constraints for the modelling of the thermal history of individual samples (see below).
FIG. 4. Scheme for the interpretation of K-Ar illite fine fraction ages and polytype quantification data in combination with regional cooling data derived from (U-Th)/He, apatite fission track and K-Ar biotite method.The illite polytype transition zone marks the irreversible transition from one respective polytype to another under prograde conditions.Thus, under retrograde conditions there will always be a mixture of polytypes in the respective metamorphic zone.The temperatures for the transition between polytype-deformation-and metamorphic regimes have to be considered as continuously with not exact constrained temperature boundaries.The ages for the fault gouge fine-fraction (2-6 µm; 2 µm and <0.2 µm) are not indicated in this figure but would lie between the maximum 100% 2M 1 age and the minimum 100% 1M d age (note that this scheme involve simplifications. Additionally scheme is mainly concerted on interpretation of retrograde illite from fault gouges in non-sedimentary rocks; IC=illite crystallinity; T C =effective closure temperature; temperatures for respective system see Reiners and Brandon (2006), approximated illite polytype transformation temperatures according to Hunziker et al. (1986), Yoder and Eugster (1955), Velde (1965), Weaver (1989); temperature ranges of brittle and ductile deformation behavior of quartz taken from Passchier and Trouw (2005), IC=values for diagenetic zone, anchi-and epizone according to Kübler (1967); closure temperature of K-Ar biotite system according to Purdy and Jäger (1976)).
Apatite fission track ages range from Early Triassic to Early Jurassic (251 to 192 Ma; Table 4 and Fig. 8).Again, no correlation between age and elevation as expected for undisturbed elevation profiles (Fitzgerald et al., 2006) is recognizable.The youngest apparent ages (196 Ma and 192 Ma) can be found in the middle part of both slopes (Table 4).
Investigated samples are characterized by distinct shortened tracks with an unimodal track length distribution and a mean track lengths between 11.9 µm and 12.9 µm with a standard deviation of 0.9-1.5 µm (Table 4 and Fig. 8).The mean etch-pit diameters (Dpar values) of the seven samples are between 1.79 µm and 2.09 µm (Table 4).
Mean apatite (U-Th)/He ages along the San Luis transect range between the Middle Triassic and the Middle Cretaceous (286 to 105 Ma; Table 5 and Fig. 7).There is no correlation between age and elevation on either the western slope or on the eastern flank (Fig. 2).The youngest age (105 Ma) was found at the base of the western slope.Generally, apatite (U-Th)/He ages are younger or overlap within their 1σ-error with their corresponding AFT age.An exception is given by APM 48-08 and 32-08, as AHe ages are older than corresponding AFT ages, contradicting the normal age trend of AFT>AHe (Table 5 and Fig. 7).

K-Ar dating
Four small-scale fault zones along the Nogolí-Río Grande transect (Fig. 3) were sampled and dated by the K-Ar illite method in several grainsize fractions.Additionally, illite polytype quantification, illite crystallinity determination and mineralogical classification of gouges were done by X-ray diffraction.
In total, twelve K-Ar ages for the grain-size fractions of <0.2 µm, <2 µm and 2-6 µm from four samples were analysed.Samples were taken from associated small-scale faults in the middle part of the western hillslope (Fig. 3), where the bestexposed fault gouges were observed.No reliable kinematic information of fault slip data sets could have been surveyed, just slickenlines geometry.
Samples were exclusively taken from clay fault gouges developed in macroscopically muscovite free crystalline host rocks.Sampled locations FIG. 5. K-Ar ages of all analyzed grain-size fractions and samples as well as extrapolated ages for the 2M 1 and 1M d polytypes.All samples were taken from well-developed clay-gouges from small-scale faults in the middle part of the western hillslope (Fig. 3).Note the increase in illite age with increase in grain-size.The 100% 2M 1 and 100% 1M d ages are interpreted to represent the begin of deformation along the fault, respectively to represent the end of deformation and/or cooling below illite formation temperatures (more details are discussed in the text).
yield well developed clay gouges, in some cases with small grains of residual rock material.The monomineralic grains are only a few millimetres in diameter and consist almost exclusively of quartz.Polymineralic rock fragments were not observed within sampled material.Gouge thicknesses range from several millimetres up to several centimetres, but mostly vary between 0.5-2 cm.Sampled gouges have monochrome brownish to reddish colours.
K-Ar ages range from Early Pennsylvanian to Early Jurassic times (315-178 Ma).The age-analysis plot of all analysed samples is depicted in figure 5 and data are listed in table 1.All samples show a time interval span between fractions ranging from 18.7 Ma up to 84.7 Ma (Fig. 5).No overlapping ages could be observed within one sample.Radiogenic 40 Ar content ranges from 80.3% to 97.7%, potassium (K 2 O) content ranges from 1.26% to 6.17% (Table 1), indicating reliable analytical conditions for all analyses samples.XRD analyses of all samples confirm that illite, smectite and kaolinite are the major clay mineral components in the various fractions (Table 3 and Fig. 6).The presence of quartz is restricted to the fractions <2 µm and 2-6 µm, and none is found in the <0.2 µm fraction.Restricted to the 2-6 µm fractions some samples show minor traces of potassium feldspar.These observations were confirmed by TEM analysis.Additionally, TEM analyses also reveal traces of halloysite (Table 2).Glycolated XRD analyses were carried out to investigate the potential occurrence of expandable mixed-layers of illite and smectite.Major amounts of illite/smectite were found in all fractions of sample APG 50-09.Except for sample APG 50-09, illite/ smectite mixed-layers are nearly absent in the 2-6 µm fractions (Table 1).Considerable amounts of smectite could only be found in the <0.2 µm fraction.

Illite crystallinity (IC)
The IC, expressed as Kübler Index (KI), of all analysed samples varies from 0.155 Δ°2θ to 0.530 Δ°2θ (Table1).KI values from the air-dried <0.2 µm fractions indicate that two samples developed under diagenetic conditions.In contrast, the fractions of <2 µm and 2-6 µm yielded anchi-to epimetamorphic values.Variations in the Δ°2θ between the glycolated and the air-dried measurements correspond to the presence of illite/smectite mixed-layers (Table 1).No systematic variation with respect to the sample location was observed.The absence of very low IC together with fact that gouges were taken from macroscopically muscovite free host-rock indicate that grain-size fractions are not contaminated by muscovite phases.

Illite polytypism
All samples used for the illite polytype quantification show a random orientation expressed by a low Extrapolated ages for hypothetical samples consisting of 100% and 0% 2M 1 illite respectively.FIG. 7. Compilation of zircon and apatite (U-Th)/He single-grain ages and AFT ages.Generally ages scatter, but correlate with the respective temperature range of the used system (ZHe>AFT>AHe ages).A correlation between age and elevation, as expected for undisturbed elevation profiles (Fitzgerald et al., 2006), can only be observed for some neighboring samples (e.g., between samples 35, 34 and 33).This is interpreted to be the result of segmentation of the transect into several, individual fault blocks (see also figure 2).
(002)/(020) illite ratio as well as low (004) illite intensity.XRD tracings of random powders from all samples presented in this study contain a mixture of illite polytypes 2M 1 , 1M and 1M d .The analysed illite fractions are composed mainly of 1M d and 2M 1 polytypes.The 1M d polytype is dominant in the <0.2 µm fractions throughout all analysed samples.In the <2 µm fraction, the 1M d and 2M 1 polytypes are observed as the dominant phases.The 2-6 µm fractions are mostly made up of 2M 1 illite.The 1M illite content is subordinate (maximum of 7%) and more or less constant for all analyzed grain-size fractions and samples (Table 1).
In TEM (transmission electron microscopy) and SAED (selected area electron diffraction) analyses on sample APM 59-09, it was not possible to identify 1M illite, which may be attributed to the small TEM sample volume.However, other authors also reported difficulties to find the 1M polytype proposed by XRD studies during TEM analysis (e.g., Peacor et al., 2002;Solumn et al., 2005).
The relationship of increasing K-Ar ages with increasing grain size (see above and Table 1 and Fig. 5) is consistent with an increasing content of the 2M 1 illite polytype, which was formed in the earlier fault history under higher temperatures.The increase in 2M 1 is accompanied by a decrease in 1M d illite, as it was formed during the late fault history under lower temperatures.Additionally, polytype content correlate with obtained KI values, showing smaller KI-values (higher crystallinity) for samples with higher 2M 1 content (Table 2).
Based on the calculated polytype compositions of the samples, we extrapolated the "end-member" age of the 1M d polytype and the 2M 1 polytype (hypothetically samples which consist of 0% 2M 1 illite and 100% 2M 1 illite, respectively) by plotting the age of each individual grain-size fraction of a fault gouge sample against the 2M 1 illite content (Fig. 9 and Table 1).These plots show a coefficient of determination (R²) always larger than 0.9, confirming a clear linear relationship between age and 2M 1 polytype content.
Potential errors in polytype quantification may have been caused by smectite.XRD reflections from smectite exhibit overlap with illite polytype specific Presented ages are central ages ±1σ (Galbraith and Laslett, 1993); ages were calculated using zeta calibration method (Hurford and Green, 1983); glass dosimeter CN-5, and zeta value of S.L. is 323.16±10.1A cm-2 ; zeta error was calculated using ZETAMEAN software (Brandon, 1996); n: number of dated apatite crystals; ρ s /ρ i : spontaneous/induced track densities (×105 tracks cm -2 ); ρ d : number of tracks counted on dosimeter; N s /N i : number of counted spontaneous/induced tracks; N d : number of tracks counted on dosimeter; P(Χ 2 ): probability obtaining chi-squared value (Χ 2 ) for n degree of freedom (where n is the number of crystals -1); MTL: mean track length; SD: standard deviation of track length distribution; N: number of tracks measured; Dpar: etch pit diameter.
reflections (Fig. 6; Grathoff and Moore, 1996).In the case of superposition of smectite reflections on illite peaks, the illite content might be overestimated.As a consequence, polytype quantification as well as extrapolated 100% 2M 1 and 0% 2M 1 illite ages are subjected to error, possibly greater than the conservatively assumed general methodological error for polytype quantification of 2-5% (Grathoff and Moore, 1996).However, polytype quantification is in good accordance with other parameters, such as grain-size age, illite crystallinity and K-Ar age, indicating the consistency of the data set.Even most of the extrapolated polytype end-member ages (Table 3) are in good accordance with K-Ar mica cooling ages (see discussion).Non-deformational illite formation by fluid percolation cannot be excluded but is unlikely due to the consistency of the data set.

Thermal modelling
Based on the individual cooling paths derived from HeFTy modelling, a regional thermal history for the entire transect was compiled and is shown in figure 10.Cooling below the PRZ Z temperatures (≈175 °C) over the whole transect started in Late  The temperature regime for the apatite fission track partial-annealing zone (PAZ A ≈110-90 °C) was passed in Middle Permian to Early Triassic times.The lower temperature boundary recorded by our data (PRZ A ≈65 °C) was reached in late Permian to Jurassic times, or even in the Cretaceous.We attribute this great span of cooling ages to a prolongation of very low cooling rates, as evident in all models (Fig. 10).Slow cooling and associated long-lasting residence time in the PRZ Z , PAZ A and PRZ A also led to a broad scattering of individual helium and fission-track ages (Fig. 7) and the development of a broad, unimodal fission track length distribution and distinct shortened tracks (Fig. 8).Assuming that these very slow cooling rates gave rise to low rates of surface exhumation (i.e., limited erosion), suitable conditions were provided for the development of erosional surfaces during one or several regional planation events.
Based on model data, conservative calculations of cooling rates to temperatures around 175 °C (Figs. 10, 11) yields rates below 5 °C/Ma.For the temperature range of ca.175 °C (PRZ Z ) to ca. 65 °C (PRZ A ), rates vary from around 2 °C/Ma to 10 °C/Ma.An exception of slow to moderate cooling rates is given by samples APM 34-08 and 49-08, which yield rates of 0.5 °C/Ma to 1.5 °C/Ma.These very low cooling rates, together with the observation of higher ages for cooling up to ca. 175 °C (see above), strongly indicate a different thermal history of samples APM 34-08 and 49-08 in comparison to all other samples, at least for the cooling above the PRZ A .Results yielded by sample APM 34-08 are puzzling, whereas sample APM 49-08 is located at the boundary of a low relief area within the western Sierra de San Luis hillslope.Field studies have highlighted that the western range hillslope is a complex uplift scarp in the study area, composed by several morphotectonic domains bounded by fault zones (Costa, 1992;Costa et al., 1999).These geological evidences together with the mismatch of thermochronological data thus indicate that differential rates of exhumation and/or uplift resulted from segmentation and presence of several fault-bounded blocks (Fig. 12).
Very low cooling rates of <0.5 °C/Ma are calculated for all models in the temperature range of 65 °C to around 30 °C.Final cooling to the present day mean surface temperature of 17 °C (Müller, 1996) is less constrained by models but most likely in the range of 0.5 to 1.5 °C/Ma.

.2. Timing of faulting constrained by K-Ar ages
Following the method suggested by Bense et al. (2014), we interpret the 100% 1M d end-member FIG. 9. K-Ar ages versus % 2M 1 illite content for three analyzed samples.Polytype composition was determined using XRD-methods.
Plots show a coefficient of determination (R 2 ) better than 0.9, confirming a clear linear relationship between age and 2M 1 polytype content.Based on this relationship, ages for hypothetical fractions consisting of 0% 2M 1 were calculated by linear extrapolation (see text for a detailed discussion).
(mathematical age of an hypothetical sample consisting of 100% 1M illite, see above) age to best represent the age of the youngest deformational movement resulting in illite formation along sampled faults because the age-increasing influence of the "older" 2M 1 polytypes is eliminated.In turn, the complementary 2M 1 end-member age may represent a) the oldest generation of neoformed illite, under retrograde (cooling) conditions (see Bense et al., 2014) and/or b) the age of "detrital" muscovite, meaning crushed muscovite from the host rock (Fig. 13).Concerning the latter, the absence of If the 2M 1 illite age represents the onset of brittle deformation, it should always be younger than a K-Ar biotite cooling age from the same rock or nearby location (Fig. 4).Biotite shows a closure temperature around 300 °C (McDougall and Harrison, 1999) and can be interpreted to date the cooling of the basement close to brittle-ductile transition temperatures.Thus, biotite K-Ar ages represent cooling ages that predate the onset of brittle deformation and formation of illite in fault gouges (Fig. 4).
In the Nogolí region, K-Ar biotite ages from basement samples yield Carboniferous ages (345-328 Ma;Steenken et al., 2008) and are older than illite finefraction ages.Except one sample, they overlap with extrapolated 100% 2M 1 illite ages within error.
Further constraint is given by K-Ar and Ar/Ar muscovite ages (Fig. 14), yielding Devonian ages (380-350 Ma;Sims et al., 1998;Steenken et al., 2008) and Carbonifeous K-Ar fine fraction ages (345-299 Ma; Wemmer et al., 2011).The latter ages were taken from the San Luis Formation, represented by two belts of low-grade phyllites and quartz arenites, several kilometres to the north of the study area.The former mentioned muscovite ages are interpreted to represent the last mylonitization event caused by the Achalian Orogenic Cycle before brittle deformation started (Sims et al., 1998;Steenken et al., 2008), whereas K-Ar fine fraction ages from the San Luis Formation are interpreted to coincide with deformation under ductile/brittle transition conditions and the end of ductile deformation (Wemmer et al., 2011).All illite grain-size fraction ages presented here are younger or, in case of extrapolated 100% 2M 1 ages, overlap within error with K-Ar and Ar/Ar cooling ages, supporting our interpretation.Hence, the western Sierra de San Luis basement achieved depth/temperature levels below the brittle/ductile transition zone due to active deformation and exhumation during the latest Carboniferous/earliest Permian.
Illite-generating fault gouge activity along the sampled faults is interpreted to have ceased between 222 Ma and 173 Ma, as shown by the majority of the <0.2 µm grain-size fractions as well as the calculated 100% 1M d illite fractions (Table 1).Only one sample (APG 60-09) shows a younger age for the calculated 100% 1M d fraction of 119 Ma.This comparatively young age might be related to an extrapolation error resulting from its-compared to the other samples-low 1 M d illite content in the <0.2 µm fraction (Table 1).The youngest age documented by fault gouge dating must not be considered to represent the cessation of fault activity but to represent the last illite-forming event and, thus, the cooling of the fault block below illite-forming temperatures (approximately between 75-110 °C; e.g., Hamilton et al., 1992;Fig. 4).Cooling below the illite-forming temperatures is constrained by AFT and AHe ages (Tables 4 and 5).The youngest illite must overlap with the apatite fission track ages (representing cooling below 110 °C), whereas the AHe ages (representing cooling below 60 °C) must always be younger than the K-Ar illite ages (Fig. 4).This can indeed be observed for all analysed samples (Fig. 14; Tables 1-3, see also Bense et al., 2014).

Implications for the evolution of the Sierra de San Luis
The onset of brittle deformation in the western part of the present-day Sierra de San Luis during the Late Devonian to Early Carboniferous, suggest that this rock massif was already emplaced at shallow crustal levels by that time.Although there are no landscape remains of this evolutionary stage, it is suggested that surface uplift of the area related to the present day range configuration, driven by faulted blocks uplift, started during that time.This crustal mobility may be related with the development of shear zones under brittle-ductile conditions in several parts of the range during the Achalian Cycle, whose effects according to Sims et al. (1997Sims et al. ( , 1998) ) can be documented up to 355 Ma.Such significant episode of exhumation and uplift may also be linked to the collision of the Chilenia exotic terrane at the proto-pacific margin of Gondwana (Chanic orogeny) during the late Devonian (Ramos, 1988;Mpodozis and Ramos, 1989;Davis et al., 1999;Willner et al., 2011).
During the Carboniferous to the Early Permian (Fig. 12a), a mountainous landscape developed as suggested by Jordan et al. (1989), associated with the down-wearing of the Early Carboniferous topography.Aiming to reconstruct possible evolutionary paths, this primary landscape is represented as an imaginary primary topographic enveloping surface (PTES; Fig. 12).
Several intermountain basins of unknown extent and dimensions developed during this phase, which record fluvial-lacustrine sedimentation (e.g., fault-bounded depressions of Bajo de Véliz and La Estanzuela; Hünicken et al., 1981;Limarino and Spalletti, 2006;Chernicoff and Zappettini, 2007).These depocenters were part of the eastern Paganzo basin, which contains basement-derived clasts with provenance mainly from the east, indicating surface uplift of the easternmost Sierras Pampeanas to some extent during basin development (López-Gamundí et al., 1994;Chernicoff and Zappettini, 2007).A similar situation with a Carboniferous landscape carved into older rocks and covered or filled by fluvial sediments was documented in the Sierra de Chepes (La Rioja province) by Enkelmann et al., 2014).
Based on thermal modelling results, the middle Permian to early Jurassic can be identified as main exhumation phase, comprising around 40-50% of the total exhumation recorded by thermochronological data and can be ascribed to the effects of the San Rafael orogeny and the subsequent extensional orogenic FIG. 13.Illite polytypism in fault gouges in relation to temperature conditions during prograde and retrograde metamorphism.Generally when interpreting illite polytypism from fault gouges prograde metamorphic conditions are assumed.As a consequence 2M 1 illite is considered a detrital component derived from source rock, because temperature conditions in fault zone are regarded as insufficient for the development of 2M 1 illite (i.e., ca.200 °C).However, sufficient temperature for 2M 1 illite development can be ensured when the host rock experienced brittle deformation during regional cooling and passage of epizonal conditions during retrograde metamorphism (Bense et al., 2014).Thus, in case of non-sedimentary host-rocks, the 2M 1 illite must not be excluded from consideration as meaningful fault gouge age.
collapse (Ramos, 1988;Mpodozis and Ramos, 1989;Kleiman and Japas, 2009) in the study range.Upper Paleozoic to Jurassic exhumation was also reported for the Sierras de Córdoba, Sierra de Pie de Palo and Sierra de Aconquija (Löbens et al., 2013a, b;Richardson et al., 2013).In combination with the fact that Gondwanic sediments were preserved in the Bajo de Véliz region, data presented herein evidence a differential exhumation in the Sierra de San Luis.While the eastern range area (Bajo de Véliz) was exhumed to a surface or a near-surface position, allowing interaction with surficial processes, the western part (study area) was at least at several kilometres below the present surface by that time.This fact suggests that exhumation related to differential block movement driven by faulting was active at least during Permian times.A mountainous landscape probably prevailed by then in this region, whereas planation processes, which led to the present-day erosional surfaces and flat topography, might not start before the Late Permian at the western San Luis range (Fig. 12b).
Triassic exhumation and deformation is recorded by K-Ar fault gouge and thermochronological data, resulting from rifting in the Andean foreland during extension (Kay et al., 1991;Giambiagi et al., 2011;Sato et al., 2015).Coeval volcaniclastic rocks are recorded in the Sierra de Varela area, ca. 100 km to the south of the Sierra de San Luis with a half-graben array, though the timing of rifting is not well known (Costa et al., 1998b).
Although final exhumation to the surface level is not well constrained by the data, thermal modelling (Fig. 11, Tables 4 and 5) indicates that cooling below temperatures of about 60 °C for the study area did not occur up to the Upper Cretaceous (Fig. 12c).On the other hand, thermal modelling shows cooling below temperatures of 50-30 °C for all samples during the Cenozoic (Fig. 10).Considering a surface temperature of 17 °C (Müller, 1996), Cenozoic cooling rates between 0.2-0.5 °C/Ma can be estimated for the study area.This may be equivalent to a total Cenozoic exhumation of 0.6-1.8km, if a geothermal gradient of 18 to 23 °C/km is assumed based on recent estimations (Dávila and Carter, 2013).Even if a lower geothermal gradient is considered, as indicated by Collo et al. (2011;2015) for the Bermejo and Vinchina basins (~400 km to the northwest of the study area), maximum possible Cenozoic exhumation is constrained at 0.8-2 km.If the exhumation since the Cretaceous would have exceeded this difference, samples from the footslope area of the range must yield younger, reset ages.Additionally, this would be represented by a break in the slope in the age-elevation plot (Gleadow and Fitzgerald, 1987).
Thermochronological data indicate no significant/ measurable exhumation after Cretaceous times, suggesting that sampled rocks were already at or near surface by the Cretaceous or even before.Upper Cretaceous basaltic rocks emplaced near or at the surface further indicate that at least during the late Cretaceous the upper parts of the Sierra de San Luis were already exposed as a positive area prone to denudation processes (López and Solá, 1981).Following England and Molnar (1990) definitions, Cenozoic morphotectonic evolution might thus result from surface uplift of exhumed surfaces and/ or subsidence rather than from exhumation.There is no doubt that range uplift leading to the present day landscape occurred during the Neogene, due to an eastward migration of the locus of deformation (Ramos, 1988;Costa, 1992;Costa and Vita-Finzi, 1996;Costa et al., 2001b;Ramos et al., 2002).As a result, morphogenetic processes related to base level adjustment of fluvial systems seem to have dominated the present day morphology of the western slope (Costa, 1992).Evidences of a significant landscape rejuvenation starting in the Neogene are suggested by preservation of wide planation surfaces atop the range; landforms association along the western hillslope and piedmont indicating significant channel downcutting; thrusting of Neogene and Quaternary strata by the range bounding San Luis Fault System (Costa, 1992;Costa et al., 1999Costa et al., , 2000Costa et al., , 2001b)).
K-Ar fault gouge ages do not shed light on the crustal stress regime, but constrain the onset of brittle deformation in the western part of the range to the Late Carboniferous/Early Permian.This time span matches with regional tectonic processes such as the development of the Paganzo Basin (Limarino and Spalletti, 2006;Gulbranson et al., 2010) and the San Rafael Orogeny (Ramos, 1988;Kleiman and Japas, 2009;Geuna et al., 2010;Japas et al., 2013).Active deformation during basin development was accompanied by exhumation of the study area, as indicated by K-Ar fine fraction data from Wemmer et al. (2011) and thermochronological data presented here (Figs. 7,10).Ongoing Permian deformation and associated exhumation is recorded by K-Ar fault gouge and thermochronological data (Figs. 5,7,10), and might be related to basin inversion associated with the San Rafael Orogeny (Ramos, 1988;Kleiman and Japas, 2009;Geuna et al., 2010;Japas et al., 2013).
Carboniferous exhumation and uplift related to active deformation was also reported for the Sierra de Chepes (Fig. 1; Enkelmann et al., 2014), whereas coeval fault activity was recognized for the Sierras de Córdoba (Whitmeyer, 2008) and the Sierra de Comechingones (Löbens et al., 2011) based on Ar/Ar pseudotachyllite and K-Ar fault gouge data, respectively.Permian exhumation and deformation, in turn, is also recorded in other areas of the Sierras Pampeanas (Coughlin et al., 1998;Bense et al., 2013Bense et al., , 2014;;Sato et al., 2015).

Development of planation surfaces
It is assumed that regional tectonic unrest during the Triassic did not favour the development of planation surfaces of regional significance.In addition, a Triassic emergence of the sampled area is not supported by AHe ages, which yield predominantly Jurassic and Cretaceous ages (see above).
During Late Jurassic-Early Cretaceous times, rifting took place east and west of the Sierra de San Luis (Gordillo andLencinas, 1967, 1979;González and Toselli, 1973;Yrigoyen, 1975;Schmidt et al., 1995).As a consequence, regional planation processes took place most likely during the Lower-Middle Jurassic and/or during the Late Cretaceous, although data presented here cannot discriminate if major erosion surfaces are diachronous in time, as suggested by Carignano (1999) and Rabassa et al. (2010Rabassa et al. ( , 2014)).AHe ages indicate that both time intervals contributed to the exhumation of the sampled area.Samples from higher elevations in the eastern part of the study area show predominantly Jurassic ages (Table 5) and were probably also emerged during this time.Instead, samples representing the footslope area of the western range passed through PRZ A conditions in middle Cretaceous times (AHe mean age of 105 Ma, see Table 5).Whether or not emergence occurred during this time cannot be concluded from the data.
Further constraints comprise basalts emplaced at higher altitudes in the Sierra de San Luis to the northeast of the study area (Fig. 12c), which yield K-Ar ages of 83±6 Ma (López and Solá, 1981).The effusive character of the dated basalts is unclear because no robust field evidence for effusion above erosional surface has been reported and they might alternatively comprise a subvolcanic intrusion.If these basalts were emplaced atop the erosional surface, the emergence and formation of the erosional surface should be considered to be older than 83±5.85Ma.In this case, final exhumation to surface of the present day top of the Sierra de San Luis occurred between the Jurassic and Upper Cretaceous.If, in contrast, the basalts represent subvolcanic bodies with no relationship to a paleosurface, no further interpretations regarding the final exhumation can be made with available data.
Ages of basaltic rocks are in contrast to AHe single grain ages of the lowermost sample from the western footslope (APM 49-08), which show a minimum age for passage through the PRZ A of 72 Ma (Table 5), coinciding with a depth of about 2.4 km (considering a geothermal gradient of 25 °C/km).This favours the idea of diachronous development of the erosional surfaces in this part of the range due to differential block uplift, with an older surface on top of the range (>83 Ma) and a younger erosional surface (<72 Ma) barely developed in secondary blocks located at the main western range hillslope.

Conclusions
Exhumation of the section of the Sierra de San Luis studied here started during Carboniferous times.The middle Permian and Triassic to early Jurassic can be identified as the main exhumation phase, comprising around 40-50% of the total exhumation recorded by the applied thermochronological methods.Cooling rates varied between 2-10 °C/Ma during the Permian and Triassic periods.Post-Triassic cooling yields lower rates of 0.5-1.5 °C/km.Generally, exhumation took place differentially across the range, propagating from east to the west.
K-Ar fault gouge and thermochronological data reveal that crustal deformation and associated exhumation during the main exhumation phase were related to main tectonic processes.The Carboniferous evolution could be related to the development of the southeasternmost Paganzo Basin, whereas Permian exhumation is here associated with basin inversion during the San Rafael Orogeny.Crustal extension and rifting, in turn, gave rise to Triassic deformation and exhumation.
Because the range already cooled to near surface temperatures in Jurassic to Cretaceous times, thermochronological data give no explicit evidence for any Cenozoic exhumation, although Neogene range uplift is evident by geological and geomorphological data.Relative steady conditions persisted in the area nowadays, documented by planation surface remnants since the Cretaceous, and Apatite ages, as samples passed through the PRZ A predominantly in Jurassic to Cretaceous times.This episode probably comprises the period when most of the remaining erosional surfaces were developed.The final structural relief and present topography of the Sierra de San Luis has been achieved as a result of Neogene regional crustal shortening and consequent differential uplift of surfaces which were already exhumed during Late Paleozoic and predominantly Mesozoic times.
assisting in the field.Field work done by F. Bense was financially supported by the DAAD (project number D/08/48018).

K-Ar fault gouge: data acquisition
Each sampled fault gouge consists of approximately 250-1,000 g of fresh material.After careful selection, about 200 g of clay material was dispersed in distilled water and sieved <63 µm.The grain-size fractions <2 µm and 2-6 µm are gained from <63 µm fraction by differential settling in distilled water (Atterberg method following Stoke's law).Enrichment of the grain-size fraction <0.2 µm was accelerated by ultracentrifugation.For details the reader is referred to Wemmer (1991).Following the concentration of the different grain-size fractions, the samples were subjected to isotope measurements for dating and X-ray diffraction for determining the mineralogy, IC and PT.
The XRD analysis was done using a Phillips PW 1800 X-ray diffractometer.For the identification of the mineral content a step scan (0.020°2θ) was performed in the range from 4 to 70°2θ.This was important to verify the existence of illite and to rule out the occurrence of any other potassium-bearing minerals.The mineral composition of all samples is shown in table 2.
For the determination of the IC ''thin'' texture compounds according to Weber (1972) were prepared using 1.5 to 2.5 mg/cm 2 sample material.The metamorphic grade of the samples has been inferred from the peak width at half-height of the 10 Å peak (Kübler, 1967) using a software algorithm developed at the University of Göttingen by Friedrich (1991), rewritten to FORTRAN by K. Ullemeyer (Geomar, Kiel) in 2005.
Digital measurement was carried out by a step scan with 301 points in a range of 7-10°2θ, using a scan step of 0.01°2θ and an integration time of 4 s per step and a receiving slit of 0.1 mm as well as an automatic divergence slit.The presence of mixed layer clays that may obliterate the 10 Å peak has been tested by duplicate determination of the material under air-dry and glycolated conditions.Results of the IC determination are shown in table 3.
To determine the polytypism of illite, powder compounds were prepared and scanned in 561 steps in a range of 16-44°2θ, using a scan step of 0.05°2θ and an integration time of 30 s per step.The allocation of peaks to the corresponding polytypes was done as suggested by Grathoff and Moore (1996) and Grathoff et al. (1998).The randomness of orientation for the powder sample preparation was checked using the Dollase factor (Dollase, 1986).We tried to quantify the amount of different polytypes by several methods as described by Reynolds (1963); Maxwell and Hower (1967);Caillère et al. (1982); Grathoff and Moore (1996) and Grathoff et al. (1998), but due to bad peak shapes none of the used methods yield reasonable results.For this reason the general abundance of the polytypes was estimated as suggested by Grathoff and Moore (1996).The abundance of different polytypes in the analysed samples is shown in table 3.
Potassium and argon were determined following two different procedures.The argon isotopic composition was measured in a Pyrex glass extraction and purification vacuum line serviced with an on line 38 Ar spike pipette and coupled to a VG 1200 C noble gas mass spectrometer operating in static mode.Samples were pre-heated under vacuum at 120 °C for 24 h to reduce the amount of atmospheric argon adsorbed onto the mineral surfaces during sample preparation.Argon was extracted from the mineral fractions by fusing samples using a low blank resistance furnace within the Pyrex glass extraction and purification line.
The amount of radiogenic 40 Ar was determined by isotope dilution method using a highly enriched 38 Ar spike (Schumacher, 1975), which was calibrated against the biotite standard HD-B1 (Fuhrmann et al., 1987;Hess and Lippolt, 1994).The released gases were subjected to a two-stage purification procedure via Tigetters and SORB-ACs getters.Blanks for the extraction line and mass spectrometer were systematically determined and the mass discrimination factor was monitored by airshots.The overall error of the argon analysis is below 1.00%.Potassium was determined in duplicate by flame photometry using an Eppendorf Elex 63/61.The samples were dissolved in a mixture of HF and HNO 3 .CsCl and LiCl were added as an ionisation buffer and internal standard, respectively.The pooled error of duplicate Potassium determination on samples and standards is better than 1%.
The K-Ar age were calculated based on the 40 K abundance and decay constants recommended by the IUGS quoted in Steiger and Jäger (1977).The analytical error for the K-Ar age calculations is given at a 95% confidence level (2σ).Analytical results are presented in table 3. Details of argon and potassium analyses for the laboratory in Göttingen are given in Wemmer (1991).

FIG. 1 .
FIG. 1. Digital elevation model of the Sierras Pampeanas region where fault-bounded blocks surrounded by intermontane basins stand out in central western Argentina (see inset for location).Red dotted line shows boundaries of this geologic province with the main Andean building.Red rectangle highlights the Sierra de San Luis area depicted in figure 2 (NASA SRTM data http:// www2.jpl.nasa.gov/srtm/southAmerica.htm#PIA03388) (last visit 01/09/2016).

FIG. 8 .
FIG. 8. Track length distribution of the investigated samples.The amount of measured confined tracks is given by n.The apparent age of each sample is shown in italic letters.Elevation and FT parameters of samples are shown in table 4.

FIG. 10 .
FIG. 10.Possible t-T paths for individual samples based on Fission track single grain ages and the confined track length distribution as well as the (U-Th)/He ages of apatite and zircon; light grey paths: acceptable fit, dark grey: good fit, black line: best fitting path.With exception of APM 35-08 all models are based on (U-Th)/He data from zircon and apatite and Apatite Fission track data.Model 35-08 is only based on (U-Th)/He data.

FIG. 11 .
FIG. 11.Thermal history for the San Luis profile based on individual t-T paths presented in figure 10.Horizontal error bars decipher the range of good fits, the position of the marker represents position of best fit path (Fig. 10).For a better illustration the individual samples are exhibited at slightly different temperatures.The approximate temperature ranges of the Partial Retention Zone for zircon (PRZ Z ), the Partial Annealing Zone for apatite (PAZ A ) as well as the Partial Retention Zone for apatite (PRZ A ) are shown as grey bars.The high range in ages of sample APM 35-08 is due to lack of Apatite Fission Track data (Fig. 10).

FIG. 12 .
FIG. 12. Schematic exhumation model for the Sierra de San Luis and adjacent areas.See text for details (BV: Bajo de Véliz; PTES: Previous Topographic Eroded Surface; SSL: Sierra de San Luis; CD: Conlara Depression; CR: Comechingones range; SCR: Sierra Chica Range; white arrow indicating area affected by uplift or subsidence).For further details on the Sierra de Comechingones see Löbens et al. (2011).

FIG. 14 .
FIG. 14. Compilation of available geochronological data of the study area in comparison to K -Ar illite ages from fault gouges.For better view geochronometers are displayed in the x-axis.K-Ar muscovite and biotite data are taken from Steenken et al. (2008) and Wemmer et al. (2011).Black lines indicate individual ages, bars indicate overlapping error of individual ages.