Geochemistry, U-Pb SHRIMP zircon dating and Hf isotopes of the Gondwanan magmatism in NW Argentina: petrogenesis and geodynamic implications

We have carried out zircon U-Pb SHRIMP dating and Hf isotope determinations as well as geochemical analyses on three plutonic units of Gondwanan magmatism that crop out in NW Argentina. Two episodes of different age and genesis have been identified. The older one includes gabbros and diorites (Río Grande Unit) of 267±3 Ma and granitoids (belonging to the Llullaillaco Unit) of 263±1 Ma (late Permian, Guadalupian); the parent magmas were generated in an intraplate environment and derived from an enriched mantle but were subsequently contaminated by crustal components. The younger rocks are granodiorites with arc signature (Chuculaqui Unit) and an age of 247±2 Ma (middle Triassic-Anisian). Hf isotope signature of the units indicates mantle sources as well as crustal components. Hf model ages obtained are consistent with the presence of crustal Mesoproterozoic (mainly Ectasian to Calymnian (TDM(c)=1.24 to 1.44 Ga-negative Hf (T)) and juvenile Cryogenian sources (TDM=0.65 to 0.79 Ga-positive ɛHf (T)), supporting the idea of a continuous, mostly Mesoproterozoic, basement under the Central Andes, as an extension of the ArequipaAntofalla massif. The tectonic setting and age of the Gondwanan magmatism in NW Argentina allow to differentiate: a. Permian intra-plate magmatism developed under similar conditions to the upper section of the Choiyoi magmatism exposed in the Frontal Cordillera and San Rafael Block, Argentina; b. Triassic magmatism belonging to a poorly known subduction-related magmatic arc segment of mostly NS trend with evidence of porphyry type mineralization in Chile, allowing to extend this metallotect into Argentina.


Introduction
A long and discontinuous Permian/Triassic magmatic belt occurs along the western margin of Gondwana in South America, outcropping from central Perú as far as approximately 39ºS in Argentina. From west to east it extends approximately 600 km from the Gondwana margin into the foreland, reaching the longitude of the present Sierras Pampeanas and the western margin of the Río de la Plata craton (Fig. 1, inset). It comprises epizonal plutonic and volcanic rocks, the latter including pyroclastic facies. In Argentina and Chile the assemblage has been named Choiyoi magmatic province (Kay et al., 1989;Llambías et al., 1993;Llambías, 1999), even though originally the term Choiyoi was used to identify only the volcanic units (originally Groeber, 1946, 1951, andlater Choiyoi Group;e.g., Llambías et al., 1993).
In Argentina, the Choiyoi igneous rocks crop out in the basement of the Neuquén basin and the main Cordillera of southern Mendoza, in the Cordillera Frontal of Mendoza and San Juan, San Rafael block and its southern extension in the La Pampa province. To the NW the magmatic province extends into Chile and to the SE into northern Patagonia. The Cordillera Frontal outcrops are the most voluminous and have been the subject of study over the last 40 years (e.g., Rolleri and Criado Roque, 1970;Mpodozis and Kay, 1992;Llambías et al., 1993, Breitkreuz andZeil, 1994;Lucassen et al., 1999). This region is characterized by felsic rocks with subordinated basic and mesosilicic rocks. Volcanic sequences dominate in Argentina and plutonic rocks in Chile, although important plutons such as the Colangüil Batholith also occur in Argentina (Llambías andSato, 1990, 1995).
Isolated Permian-Triassic plutons have been recognized in the Puna region of Salta province, Argentina, suggesting that the magmatic event reached that latitude (Zappettini and Blasco, 1998;Page and Zappettini, 1999;Poma et al., 2009).
The aim of this paper is to contribute to the knowledge of the Gondwanan magmatism through the presentation and interpretation of chemical and isotopic data of previously poorly known units that represent the northernmost outcrops identified in Argentina. We present new zircon U-Pb SHRIMP and Hf isotope data and we explore the nature and characteristics of the magma sources to better constrain a petrogenetic model and characterize the crustal components of the Puna region.
in the western margin of Gondwana are almost continuous during the Phanerozoic, although there is evidence for a time of relative quiescence during the Devonian and Early Mississippian (390 to 340 Ma) from northern Perú to southern Chile due to the development of a passive margin at those times (Bahlburg et al., 2009). In central to southern Chile this scenario is related to the prior collision of Chilenia (Willner et al., 2009).
During the early Carboniferous magmatism was reinitiated along the Gondwana margin in the central and southern Andes (Kay et al., 1989;Brown, 1991;Breitkreuz et al., 1992;Breitkreuz and Van Schmus, 1996). In particular, intrusive activity began around 330 Ma (Lucassen et al., 1999 and references therein). I-type earliest subduction-related granitoids are recorded in the Eastern Cordillera of Perú with U-Pb zircon ages at ca. 350-325 Ma (Miskovic et al., 2009), between 21° and 22°S in Chile (Ujina, Rosario, El Colorado and Quebrada Blanca porphyries, Collahuasi Group, 300 Ma;cf. Munizaga et al., 2008) Alasino et al., 2012). A subduction-related setting has been also ascribed to Permian magmatism and Triassic units have been interpreted to be the result of partial rifting and transtension following the suture of the Arequipa-Antofalla terrane and the Amazonian craton, due to stresses originated during the Pangea break-up (Kontak et al., 1990;Atherton and Petford, 1990). Although, other authors have pointed out that an early stage the magmatism was related to subduction, followed by non orogenic magmatism related to active rifting (Kay et al., 1989;Mpodozis and Ramos, 1990;Llambías and Sato, 1990). Japas (2009) andRocha Campos et al. (2011) consider that the 31°S to 36°S segment could be regarded as a transitional zone between different subduction segments after 270 Ma. The northern segment would be normal in dip angle and the southern shallower of that. It should be noted that in northern Chile, there is evidence of arc-related magmatism of Triassic age with porphyry copperrelated mineralization (243.2±2.1 to 248.7±3.3 Ma U-Pb SHRIMP; Munizaga et al., 2008).
Associated volcanic rocks have been identified to the north of the Aracar volcano, named as the Laguna de Aracar Formation (Koukharsky, 1969 2 ) that extends southward along the NW border of the Arizaro Salar. It comprises acidic volcanic and pyroclastic rocks dated at 266±28 Ma (K-Ar whole rock; Zappettini and Blasco, 1998). Equivalent volcanic sequences in Chile have been dated at 259±8 Ma and 261±9 Ma (K-Ar whole rock; Ramírez et al., 1991).
Other plutonic outcrops of the region have been assigned to the Gondwanan magmatism as an outcome of the mapping and geochronological studies herein presented. They are herein identified as Río Grande and Chuculaqui Units and were originally assigned to the Famatianian magmatic arc (Zappettini and Blasco, 1998).
We have focused the work on the description of the Río Grande, Llullaillaco and Chuculaqui Units (Fig. 1).
The Río Grande Unit form small outcrops to the East of the Llullaillaco salar and also to the north of the Río Grande salar where it is intruded by the Chuculaqui granitoids. These rocks intrude Ordovician metasediments unit (Zappettini and Blasco, 1998).
The Llullaillaco Unit constitutes three main groups of outcrops located to the NE of the Aracar volcano, to the East of the Llullaillaco salar and to the North of the Río Grande salar. Outcrops are partially covered by Cenozoic lava f lows and intrude rocks of the Río Grande Unit (Zappettini and Blasco, 1998).
The Chuculaqui Unit is located to the West in the center south extreme of the studied region. It constitutes two main bodies and minor satellital outcrops partially cover by the Cenozoic volcanites.

Río Grande Unit
This unit includes gabbros, diabases and diorites with hypidiomorphic granular to porphyritic texture, sometimes exhibiting igneous lamination, as it is frequently displayed by plagioclase feldspar. The most abundant minerals are plagioclase, pyroxene and olivine. Plagioclase crystals are zoned, locally showing patchy extinction due to complex zoning patterns and variable Ca-Na compositions, these being characteristics compatible with mixing-mingling processes; labradorite-bytownite has been recognized in the basic varieties, and andesine-oligoclase in the mesosilicic types. Ortho-and clino-pyroxene are replaced by fibrous amphibole pseudomorphs with a few crystals showing preserved cores. Olivine is present in a few samples showing partial replacement by an orthopyroxene rim. Accessory minerals include apatite (0.5 to 5.4 mm in length), sphene, zircon and opaque minerals; some are euhedral but there is a predominance of rounded shapes like drops. Locally, plagioclase cumulate textures have been recognized.

Llullaillaco Unit
It consists of red granite, with subordinated microdioritic facies. The main facies is a leucocratic reddish alkali-feldspar granite. The textures are granular and microgranular allotriomorphic, with simultaneous growth of alkali feldspars and quartz producing granophyre intergrowth; miarolitic structures are also observed. These textures indicate rapid cooling conditions and shallow emplacement. Quartz crystals usually have rounded borders although some of them show original bipyramidal habit. Alkali feldspar is strongly perthitic and plagioclase is scarce or not visible. Biotite is subordinated (about 1%) and it is frequently replaced by chlorite and epidote. Accessory minerals are zircon, apatite needles, fluorite forming mosaic crystals grouped in cavities (upholster vugs without evidence of lined), and rounded grains of magnetite. The presence of one alkali perthitic feldspar is compatible with an hypersolvus granite (Bowen and Tuttle, 1950;Tuttle and Bowen, 1958).

Chuculaqui Unit
This unit comprises gray tonalites to granodiorites as the main facies, with subordinate granites and quartz-diorites. Poma et al. (2009) described associated mafic facies in the borders of this plutonic suite, including gabbroid lenses like mafic microgranular enclaves. These are microdioritic in composition and the material is hybridized in various degrees. Field relationships indicate that during emplacement the silicic melt incorporated mafic material; additional textural evidence indicates in some localities that the granitic intrusion was probably coeval with the mafic magmatism.
The rocks are hypidiomorphic, with inequigranular medium to coarse grains. In some granitic rocks monzonitic textures are common showing poikilitic quartz and K-feldspar enclose euhedral plagioclase and prismatic amphibole.
The most abundant mineral is zoned oligoclase (An 26-28 ) with light sericitic alteration; alkali feldspar crystals (orthoclase) are perthitic. In the granodiorite facies, plagioclase occasionally shows a myrmekitic intergrowth at grain boundaries with alkali feldspar. The mafic minerals (10% to 20%) are amphibole and scarce biotite, both of them partially replaced by epidote, chlorite and associated opaque minerals. Accessory minerals are idiomorphic sphene, apatite, zircon and opaque minerals.

Geochemistry
The studied units cover a wide compositional range (Fig 2a, Table 1) from gabbros (45% SiO 2 ) to high silica granites (up to 78% SiO 2 ). In the K 2 O versus SiO 2 diagram ( Fig. 2b) the Río Grande Unit rocks plot in the medium-to high-K field, the Llullaillaco Unit is constrained to the high-K field and Chuqulaqui Unit spans the medium-and high-K fields.
The Río Grande Unit gabbros and diorites are metaluminous subalkaline basic to mesosilic rocks with high values of FeO, CaO and MgO, and variable contents of Na 2 O and Ti 2 O. Zr shows a positive correlation with SiO 2 . In the MORB-normalized trace element spider diagrams (Pearce, 1983) the rocks (Fig. 3a) display a continuous variation in HFSE contents, in particular Nb, Ta, Zr, Hf and Ti depletion. The REE diagram (Fig. 4a) is characterized by a relatively low La/Yb slope and Eu anomalies (0.81 to 1.42) coincident with the presence of cumular plagioclase as was petrographically identified.
The red granites of the Llullaillaco Unit have characteristics of evolved high-silica magma with low FeO, MgO, CaO contents and are slightly peraluminous with A/CNK ≤1.1 (Table 1). A distinguishing characteristic is their high Rb, Th and U contents and low Sr. Zircon contents show a random distribution versus SiO 2 . The discrimination diagram of Pearce (1984) shows (Fig. 5a) the granitoids plot in volcanic arc field transitional to the within plate field. In a spider diagram normalized to MORB (Fig. 3b) rocks show a consistent pattern characterized by strong depletion in Ba and Ce, and high Th. The REE diagram (Fig. 4b) is characterized by having a sharp negative Eu anomaly related to plagioclase fractionation and low La/Yb ratio suggesting the presence of a low-pressure residual mineralogy in the source.
The Chuculaqui Unit rocks span between 65 and 70% SiO 2 and are metaluminous, plotting in the VAG field (Fig. 5a). These rocks have characteristically variable La/Ta ratios even higher than 30 and Ba/Ta ratios (Fig. 5b) greater than 450. Zr values decrease with increasing SiO 2 indicating a normal negative trend. In a MORB-normalized spider diagram (Fig. 3c) all samples exhibit enriched LILE (Ba) relative to light REE and both enriched relative to HFSE Nb, Ta, Zr, and Hf, corresponding to continental arc-related signatures. The REE diagrams (Fig. 4c) show small negative anomalies in Eu (0.66 to 0.93) indicative of plagioclase fractionation and La/Yb ratios <9.

U-Pb geochronology
Dating young zircons (Phanerozoic) faces the problems of low counts of 207 Pb and the difficulty to detect deviations from slightly older cores and from subtle amounts of Pb loss. To minimize the first problem the time of counting 207 Pb was increased from 10 seconds to 20 seconds, and grains and areas of grains poor in U (<100 ppm) were avoided. Additionally we have used the TuffZirc algorithm (Ludwig and Mundil, 2002), which is largely insensitive to both Pb loss and inheritance to plot the 206 Pb/ 238 U ages corrected using the 207 Pb counts. Most of the data group reasonably in the Concordia plots, with few analyses deviating from the Concordia line. All Concordia ages are within error with the TuffZirc ages. These ages are presented as insets in the Concordia plots.

Río Grande unit
Rocks of this unit are relatively rich in zircon. All grains are 100-300 µm prisms terminated in pyramids (aspect ratio 2:1 to 5:1). The zircon grains are simple and back-scattered electron images show no evidence of older inherited cores or of younger metamorphic rims or zones. They have characteristics of magmatic grains such as the zoning and the Th/U ratios averaging 1.6 ( Table 2). Examples of dated zircon are provided in figures 6a to d. Twelve analyses group at the 206 Pb/ 238 U age of 267±3 Ma (MSWD= 0.019; 2σ) (Fig. 6e). Grains b.2-4 (232±3 Ma), b.5-3 and b.2-1 may be affected by lead loss and are not included in the age calculation. Grain b.5-1 is also not used because its age is highly discordant (-343%). This age is considered the age of crystallization of the gabbro-diorite body.

Llullaillaco Unit
Zircons from sample SA10-03 are short prisms (aspect ratio 2:1 to 3:1), 80 µm to 200 µm long, clear and colorless. Zoning is subtle but present in all grains. Some grains such as b.4-1, g7-1, and g7-3 show recrystallization patches which are brighter (richer in U) in BSE images (Fig. 7a to d). These areas are sealing fractures and represent a later event in the rock. Only one of these patches was analyzed (b.4-1) and effectively yielded the highest U content of 1,807 ppm and the lowest Th/U ratio of the sample (0.05), much lower than the Th/U average of the other zircons (0.68) (cf. Table 2). This patch also has the youngest age of 258±3 Ma (albeit within error with the other ages from this sample). The recrystallization patches may have resulted from the activity of late-magmatic fluids, which may  have occurred shortly after the cooling of the granite (up to 1 or 2 Ma later). Most of grains (n=10) group at the age of 263±1 Ma (Upper Middle Permian, Capitanian) but a slightly older population of three grains is also present at 280±2 Ma (Lower Permian, Artinskian) (Fig. 7e). The age of 263±1 Ma is considered the age of crystallization of the granitic body that has incorporated components of a slightly older continental crust, probably inherited from intra-arc contamination.

Chuculaqui Unit
U-Pb ages have been determined for zircon and titanite from a granodiorite. The sample has magmatic titanite, which is dark reddish brown occurring as shards of 500-600 µm in diameter ( Fig. 8a and b). The mineral is U-rich (average is 631 ppm). The results (Table 2) are concordant. Common lead (i.e., nonradiogenic) is up to 3.71% and average 1.98%. The Concordia age of seven analyses is 247±2 Ma (MSWD=0.76) (Fig. 8e). This titanite is Th-poor and Th/U ratios are almost constant and relatively low (average is 0.85).
The zircon population consists of long prisms 100 to 300 µm long of magmatic origin showing zoning ( Fig. 8c and d) and relatively high Th/U ratios (between 0.61 and 1.33). The 13 dated zircons (Table 2) have the same Concordia age at 246±3 Ma (MSWD=1.7).
The age obtained for the titanite is within uncertainty of the age of the zircons, suggesting that the minerals were co-magmatic.

Lu-Hf isotopes
While the zircon U-Pb age for igneous rocks represents the timing of magma crystallization, Hf isotopes allow distinguishing juvenile, essentially mantle-derived crust of a given age, having positive ε Hf (t) , from contemporary crust derived from remelting of older crust, characterized by negative ε Hf (t) . Juvenile magmas are defined as those generated from the depleted mantle or by re-melting of material recently extracted from it (Belousova et al., 2010).
The single-stage Hf model age (T DM ) value can be used for zircons with positive ε Hf (t) as a proxy for the maximum age of magma extraction from the depleted mantle. Otherwise, the two-stage Hf model (T DM (c) ) provides a first approximation to the source age of host magma from which zircon with negative ε Hf (t) crystallized. Hf model (T DM ) based on a depleted mantle source, is calculated using ( 176 Hf/ 177 Hf) i=0.279718 at 4.56 Ga and 176 Lu/ 177 Hf=0.0384, and producing a present-day value of 176 Hf/ 177 Hf=0.28325 . The T DM (c) age in zircon is calculated from the initial Hf isotopic composition of the zircon, using an average crustal Lu/Hf ratio (0.015; Griffin et al., 2004). The initial Hf composition of zircon represents the 176 Hf/ 177 Hf value calculated at the time  the zircon crystallized, using the U-Pb age previously obtained in the same spot of the same crystal. Such model ages indicate the crustal residence time for the rocks that hosted the zircon.

Río Grande Unit
Eight zircons of the Río Grande unit already dated by U-Pb were selected for Hf analyses. The grains were selected according to the degree of concordance and lower common lead content. All 176 Hf/ 177 Hf ratios are similar and the data produced negative ε Hf (t) (from -0.76 to -3.66). The average Lu-Hf model age of the zircon assuming a crustal origin (T DM (c) in table 3) is about 1.4 Ga (Fig. 9a).

Llullaillaco Unit
Ten measurements for Hf isotopes were undertaken on zircons from red granite. The ε Hf(t) measured in this fraction varies between +1.55 and -3.64, with one outlier at -6.26. This range is mostly coincident with that obtained for Río Grande Formation zircons. The Hf T DM(c) ages obtained range between 1.16 and 1.64 Ga, mostly grouped at 1.24 and 1.36 Ga (Fig. 9b).

Chuculaqui Unit
Hf isotope determinations in the dated Chuculaqui granodiorite zircons all give positive ε Hf values between +1.92 and +5.66 (Fig. 9c). Because of this dominantly juvenile input the values of ε Hf (t) and Hf model ages were calculated using the one-stage depleted mantle model (Table 3). T DM ages obtained range between 0.65 and 0.79 Ga. These data indicate that the crustal source of the granodiorite has an important juvenile component of Cryogenian age.

Magma sources
Late Paleozoic-Early Triassic magmatism timing in northern Chile shows two main peaks at about 300 Ma and 244 Ma (Munizaga et al., 2008), with porphyry Cu-Mo type mineralization being related to the younger event. Results obtained in Argentina also show two main episodes of magmatism, with different geochemical and Hf isotope signatures. The younger age obtained (246±2 Ma for Chuculaqui igneous event) is coincident with the younger age reported in Chile, but the older ones (267±3 Ma for the Río Grande, 263±1 Ma for the Llullaillaco igneous events, as well as the inherited 280±2 Ma zircon population in the latter unit) are not matched by rocks of similar ages in Chile.
When comparing the Hf results for both ages' groups ( Fig. 10), the Río Grande and Llullaillaco events have no equivalent in northern Chile, but they share a common Hf evolution trend with the 300 Ma group with more evolved signatures (Cluster 2) suggesting that both magmatic events were sourced by crustal protoliths with similar isotopic signature and Hf model ages. The younger Chuculaqui event has ε Hf (t) and age values essentially coincident with those of Cluster 3 from Munizaga et al. (2008) -which includes the Characolla granite porphyry and El Colorado dacite-both sharing sources with common isotopic imprints.
Hf data obtained from the Río Grande unit suggest the presence of a Mesoproterozoic crust, or sedimentary material derived from such crust, beneath the Puna region. The negative ɛ Hf also would indicate that these mafic magmas were derived from mantle subsequently contaminated by crustal components.
The ε Hf (t) variations in the Llullaillaco unit suggest that the melt is derived from a mixed juvenile-crustal source. This range is mostly coincident with that obtained for Río Grande unit zircons, as well as Hf T DM (c) ages that point to a Mesoproterozoic (Sunsás) crust, accordingly suggesting a common origin for both units. The chemistry of high-silica Llullaillaco rocks indicates crustal participation in their parental melts and within-plate extensional affinities.
Otherwise, positive ε Hf (t) values and T DM ages of about 0.71 Ga in the Chuculaqui unit indicate that the crustal source of the granodiorite has an important juvenile component of Cryogenian age.
From a geochemical point of view the Río Grande unit magmas would appear to have been generated in a setting with intraplate or no-compressional affinities, subsequently contaminated by crustal components as indicted by the Hf isotope data. The lack of inherited xenocrysts should be noted and may result from high magmatic temperatures.
The Llullaillaco granites are representative of a prevailing crustal component without excluding some juvenile mantelic input at the time of magmatism. The dominant crustal source and shallow emplacement are displayed by their geochemical and textural features. The REE pattern is representative of high-silica evolved magmas with LREE depletion caused by minor accessory minerals and The 'crustal' model ages (T DM (C) ) assume that the zircon's parental magma was produced from a volume of average continental crust ( 176 Lu/ 177 Hf=0.015; Griffin et al., 2004), that was originally derived from the depleted mantle. 176 Lu decay is 1.876 x 10-11 yr-1 (Scherer et al., 2001). The granitoids of the Chuculaqui unit show typical Cordilleran-type magmatic arc characteristics considering their geochemistry and their evolutionary trend. The positive ε Hf (t) values for the zircons of this unit also point to relatively juvenile inputs with little presence of recycled crust sources.
The data obtained are consistent with the results presented by Munizaga et al. (2008) in Chile, confirming the observation that the magmas show contributions from inhomogeneous older crust material as well as variable magmatic sources, with mixture of crustal melts and mantle-derived magmas.

The age of the crust in the Puna region: Isotopic constraints
The two stage Depleted Mantle Mesoproterozoic model ages obtained for the Permian Río Grande and Llullaillaco units indicate a significant residence time in the crust for the magma sources and their emplacement on continental crust. This magmatism in northwest Argentina thus confirms the presence of Mesoproterozoic crustal components mainly Ectasian to Calymnian (1.24 to 1.44 Ga-negative ε Hf (t) ) supporting the idea of a continuous basement under the Central Andes, as an extension of the Arequipa-Antofalla massif (Fig. 11). Coincidentally, Zappettini and Santos (2011) have reported Ectasian and Calymnian T DM (c) Hf ages (1.36 and 1.48 Ga) with negative ε Hf (t) (between 0 and -22.3) from Cretaceous syenitic intrusions in Eastern Puna. Preliminary data obtained by the authors from Ordovician granitoids from Western Puna also point to Calymnian T DM (c) Hf ages (between 1.42 and 1.62 Ga) and negative ε Hf (t) (between 0 and -3.39).
The T DM (c) Hf ages of zircons at around 1.4 Ga obtained by Munizaga et al. (2008) for the Permian magmatism in Chile are also indicative of the presence of the aforementioned Mesoproterozoic (Ectasian to Calymnian) crustal source.
Additionally, the 650 to 790 Ma T DM zircon Hf ages for the Chuculaqui Formation as well as similar data reported by Munizaga et al. (2008), with positive ε Hf (t) and T DM Hf varying between 604 and 748 Ma (recalculated from reported data by Munizaga et FIG. 10. ε Hf (T)    and the data obtained suggest that such enriched mantle was active beneath the Puna region. It should be noted that the Chuculaqui Unit lacks relict or inherited zircon that would point to reworking of ancient crust. This excludes the obtained T DM Hf ages to be interpreted as a binary mix of such an ancient crust and mantle material at the time of the generation of the granodiorite, although some Triassic juvenile input is not precluded. The extension of the Arequipa-Antofalla massif beneath the Altiplano, as indicated by the reported Hf data, is consistent with previous Pb-isotopic signature of plutonic rocks from western Puna that confirms the presence of Sunsás-San Ignacio age crust in NW Argentina (cf. Ramos, 2008;Fig. 5 and references therein). Loewy et al. (2004) have defined two domains as constituting the Antofalla basement: the northern one that extends in Chile between Belén and Sierra Moreno includes juvenile magmatism at 1.5-1.4 Ga, evidence of metamorphism at 1.2-1.0 Ga and later FIG. 11. Main basement blocks with ages from basement outcrops and model ages from the Gondwanan igneous rocks (this paper and Munizaga et al., 2008). Other basement ages from Ramos (2008). magmatism at about 500-400 Ma; and the southern domain, exposed from Limón Verde in northern Chile to Antofalla in western Argentina Puna (22°-26°S), including juvenile material of 700-600 Ma and magmatism and metamorphism recorded at 500-400 Ma.
The studied area pertains to the southern domain, but none of these latter events have been identified, except for the 700 to 600 Ma juvenile magmatic episodes that could be correlated with the 710 Ma T DM zircon Hf ages identified in the Chuculaqui Unit. The 1.24 to 1.36 Ga and 1.0 Ga events would extend the tectonostratigraphic history of this region at least to the Mesoproterozoic and, considering the isolated Statherian T DM (c) Hf age (1.64 Ga) with negative ε Hf (t) , up to the late Paleoproterozoic. The Arequipa-Antofalla massif is interpreted to be allochthonous to Amazonia (cf. Loewy et al., 2004) and that its accretion to Amazonia would have taken place during the Grenville-Sunsás Orogeny (1.0-1.3 Ga) (Chew et al., 2007). The recorded evidence of Early Mesoproterozoic crust in its southern extension could provide additional evidence for its connection to the Maz and Río Apa blocks, constituting together the hypothetical MARA craton, following the model proposed by Casquet et al. (2009Casquet et al. ( , 2010. The 1.24 to 1.36 Ga crustal ages obtained are coincident with the Andean-type magmatic arc recorded from the Maz terrane at 1.26 to 1.33 Ga (Casquet et al., 2011 and references therein), as well as to the thermal episode at 1.3 Ga that affected the Río Apa block (Cordani et al., 2010), pointing to a common history for the three areas at least during the Mesoproterozoic.

Concluding remarks
Gondwanan magmatism developed between Early Carboniferous and Early Triassic times. In NW Argentina it comprises two episodes of different age and genesis: the oldest includes gabbros and diorites (Río Grande Unit) and granitoids (belonging to the Llullaillaco Unit) of late Permian age (Guadalupian) generated in an intraplate environment, hypothetically from an enriched mantle subsequently contaminated with crustal components; the youngest is represented by granodiorites (Chuculaqui Unit) of middle Triassic age (Anisian) with Cordilleran-type arc signature. The results obtained here are comparable to those presented by Munizaga et al. (2008) for northern Chile (see discussion above) and can be related to the Pre-Andean cycle as distinguished in Chile by Charrier et al. (2007).
The Choiyoi magmatic province, according to Llambías and Sato (1995) and Llambías (1999), developed in a tectonic setting evolving from a late Carboniferous to Permian subduction-related magmatic arc through a collisional regime and subsequent early Triassic post-orogenic granite magmatism, the latter developed in the Argentine side of the Frontal Cordillera. Japas (2009) andRocha Campos et al. (2011) subdivided the Choiyoi magmatism in a lower section (with ages between 280 and 265 Ma) and an upper section (with ages between 265 and 250 Ma). Triassic volcanic sequences are separately grouped in a synrift phase. Late Triassic to early Jurassic volcanic sequences are grouped further south, in the basement of the Neuquén Basin, into the Precuyano Cycle related also to a rift setting.
We consider that the upper section of the Choiyoi magmatism, with main outcrops in the Frontal Cordillera and San Rafael Block, reaches the NW Argentine Puna where it is represented by the Río Grande and Llullaillaco units. In fact, the tectonic conditions that originated both the Permian magmatism described in this paper and that of the upper section of the Choiyoi series, as well as their geochemical signatures, are analogous. The ages determined in the Puna units are slightly older than those known from further south in the Frontal Cordillera and San Rafael Block (Kleiman and Japas, 2009), suggesting that these common tectonic conditions were originated earlier in the north, with progressive migration to the south.
Conversely, the Triassic magmatism represented by the 246 Ma Chuculaqui Unit in the Puna region cannot be correlated with the upper section of the Choiyoi magmatism from Frontal Cordillera and San Rafael Block, the latter being older than 250 Ma (Rocha Campos et al., 2011). Furthermore, the Triassic magmatism from the Gondwanan Cycle of similar age (246 Ma and younger), identified in the San Rafael Block, has been grouped in the above mentioned synrift phase (Rocha Campos et al., 2011) with geochemical characteristics and tectonic setting emplacement that differ from those identified for the Chuculaqui unit. Rather, the latter could be ascribed to a poorly known continental magmatic arc segment of mostly NS trend. Although during the Pre-Andean tectonic cycle subduction along the continental margin was presumably interrupted or considerably diminished (Charrier et al., 2007) the presence of the Chuculaqui Unit and similar rocks of the same age in northern Chile (Munizaga et al., 2008) are indicative of, at least, some restricted subduction related magmatism during the Triassic.
Interestingly, it should be highlighted that in Chile there are porphyry Cu type deposits (La Profunda and Characolla) genetically related to these rocks (Munizaga et al., 2008). This implies the presence of an early Triassic metallogenic belt that continues into Argentina, in the areas where the Chuculaqui magmatic event has been defined.