Abstract

Western South America provides an outstanding laboratory for studies of magmatism and crustal evolution because it contains Archean–Paleoproterozoic cratons that amalgamated during Neoproterozoic supercontinent assembly, as well as a long history of Andean magmatism that records crustal growth and reworking in an accretionary orogen. We have attempted to reconstruct the growth and evolution of western South America through U-Pb geochronologic and Hf isotopic analyses of detrital zircons from 59 samples of sand mainly from modern rivers. Results from 5524 new U-Pb ages and 1199 new Hf isotope determinations are reported. Our data are integrated with previously published geochronologic and Hf isotopic information, yielding a combined record that includes >42,000 ages and >1900 Hf isotope analyses. These large data sets yield five main conclusions: (1) South America has an age distribution that is similar to most other continents, presumably reflecting processes of crustal generation and/or preservation related to the supercontinent cycle, with age maxima at 2.2–1.8 Ga, 1.6–0.9 Ga, 700–400 Ma, and 360–200 Ma; (2) <200 Ma magmatism in western South America has age maxima at ca. 183, 166, 149, 125, 110, 88, 65, 35, 21, and 4 Ma (with significant north-south and east-west variations), yielding an average cyclicity of ∼33 m.y.; (3) for the past 200 m.y., no correlation exists between magmatism and the velocity of convergence between central South America and Pacific oceanic plates, the age of the downgoing plate, or the absolute motion of South America; (4) Hf isotopes record reworking of older crustal materials during most time periods, with incorporation of juvenile crust at ca. 1.6–1.0 Ga, 500–300 Ma, and ca. 175–35 Ma; and (5) the Hf isotopic signature of <200 Ma magmatism is apparently controlled by the generation of evolved crust during crustal thickening and eastward arc migration, versus juvenile magmas during extensional tectonism and westward and/or outboard migration of arc magmatism.

INTRODUCTION

Understanding the generation and reworking of continental crust has been a long-standing goal in Earth science (Armstrong, 1968). Western South America provides an excellent laboratory for studying crustal evolution because it contains older cratons that record Precambrian crustal genesis and assembly and also early Paleozoic through modern magmatic assemblages that formed along the Andean convergent and/or accretionary plate boundary (Cordani et al., 2000; Aceñolaza et al., 2002; Cawood, 2005; Franz et al., 2006; Cawood and Buchan, 2007; Ramos, 2009). This history has been reconstructed through a large number of detailed studies of specific areas (e.g., Bahlburg and Hervé, 1997; Loewy et al., 2004; Hervé et al., 2007; Augustsson and Bahlburg, 2008; Chew et al., 2008; Bahlburg et al., 2009, 2011; Miskovic et al., 2009), as well as several syntheses that have focused on assembly of South America (e.g., Goodwin, 1996) and evolution of the Andean margin (e.g., Aceñolaza et al., 2002; Haschke et al., 2002; Cediel et al., 2003; Cawood, 2005; Kay et al., 2005; Franz et al., 2006; Cawood and Buchan, 2007).

This study attempts to reconstruct the magmatic history and crustal evolution of western South America through U-Pb geochronologic and Hf isotopic analysis of detrital-zircon grains from 59 samples of modern sediment (Fig. 1). These new data are combined with previously published age and Hf isotope information, resulting in a compilation with >42,000 ages and >1900 Hf isotope determinations. In an effort to analyze these large data sets quantitatively, we (1) use maximum likelihood estimation to determine the main age groups present in age distributions, (2) evaluate possible cyclicity of magmatic records using a Lomb-Scargle least-squares spectral analysis and a continuous wavelet transform (CWT) analysis, (3) investigate possible connections between plate motions and magmatic records using cross wavelet transform and wavelet coherence analyses, and (4) use least-squares regression analysis to determine Hf isotope evolution trajectories.

Quantitative analysis of these U-Pb and Hf isotopic data sets provides new insights into the magmatic history and crustal evolution of western South America. All of the records reveal cyclic patterns that relate to the supercontinent cycle for longer time scales and, for the past 200 m.y., record changes in the tectonic setting of Andean arc magmatism. Interestingly, our data do not record cycles of Andean crustal evolution on a 30–50 m.y. time scale, which have been fundamental to recent discussions of orogenic cyclicity (e.g., Haschke et al., 2002, 2006; DeCelles et al., 2009, 2015; Ramos, 2009).

Combined U-Pb geochronology and Hf isotope analysis of zircon is a powerful method for evaluating regional patterns of magmatic history and crustal evolution due to the power of sedimentary averaging (Link et al., 2005), the durability of zircon under most crustal pressure-temperature (P-T) conditions (Speer, 1982), the ability of U-Pb geochronology to record ages of felsic to intermediate-composition magmatism (e.g, Condie et al., 2011), and the petrogenetic information provided by Hf isotopes (e.g, Augustsson et al., 2006; Bahlburg et al., 2009, 2011; Hawkesworth et al., 2010; Collins et al., 2011; Dhuime et al., 2011). It is important to note, however, that the age and Hf isotopic signature of detrital zircons provides a somewhat biased record of magmatic history and crustal evolution. A more complete analysis would require consideration of (1) variations in zircon abundance related to rock composition (e.g., SiO2 and/or Zr content) (Moecher and Samson, 2006; Dickinson, 2008), (2) three-dimensional analysis of the volumes of igneous rock present in the Andes through time (e.g., Paterson and Ducea, 2015), (3) igneous rocks that lack zircon, and (4) the importance of crustal additions and removals arising from tectonic processes (e.g., terrane accretion, subduction erosion, and delamination and/or relamination) and erosion. Such analyses are beyond the scope of this paper but will be facilitated by the new and compiled information presented herein.

Our study complements the work of Rino et al. (2004, 2008), Condie et al. (2005, 2009), Mapes et al. (2005), and Iizuka et al. (2010), who presented ages of detrital zircons from modern sand samples along the Amazon and Paraná Rivers. To complement the information from large rivers, we include samples from numerous smaller river systems (e.g., Milliman and Syvitski, 1992) that are present along the eastern and western flanks of the modern Andes.

GEOLOGICAL FRAMEWORK OF SOUTH AMERICA

The South American continent consists, to first order, of several Archean–Paleoproterozoic cratons that were effectively welded together during Mesoproterozoic through early Paleozoic orogenesis (Goodwin, 1996). This continental cratonic nucleus served as the backstop for accretionary processes that formed the Andean orogen along the paleo-Pacific margin (Coney, 1992; Cawood, 2005; Cawood and Buchan, 2007). The following sections outline the geologic evolution of cratonal South America and the northern, central, and southern Andes, with an emphasis on patterns of magmatism and crustal generation and/or reworking (Fig. 1).

Cratonal Nucleus

The nucleus of South America consists mainly of the Amazon, Sao Francisco, Parana, and Rio de la Plata cratons (Fig. 1). The Central Amazon craton (CAm) has Archean ages of 3.3–2.5 Ga. The Maroni-Itacaiunas province (MI) occurs east of the CAm with radiometric ages between 2.2 and 1.9 Ga. Westward, the Ventuari-Tapajos (VT) mobile belt (2.0–1.8 Ga) and the Rio Negro Juruena mobile belt (RNJ; 1.8–1.6 Ga) formed during the accretion of intraoceanic arcs. The Rondonian San Ignacio (RSI) province (1.6–1.3 Ga) consists of accreted arc terranes and microcontinents that experienced high-grade metamorphism during Mesoproterozoic time (Teixeira et al., 1989; Tassinari and Macambira, 1999; Santos et al., 2000, 2008; Cordani and Teixeira, 2007; Cordani et al., 2009). Along the western margin of the Amazon craton is the Sunsas orogen, which yields ages between 1.2 and 0.9 Ga (Litherland and Power, 1989; Teixeira et al., 2010; Ibanez-Mejia et al., 2011). Separating these cratons are the Brasiliano mobile belts, with dominant magmatism during Neoproterozoic–earliest Paleozoic (0.7–0.5 Ga) time (Goodwin, 1996; Rapela et al., 1998, 2007; Kröner and Cordani, 2003; Valeriano et al., 2008; Cordani, 2009; Cordani et al., 2009).

Northern Andes

The northern Andes (NA) are divided into three separate mountain ranges that consist mainly of Precambrian basement (Cordani et al., 2005; Ordóñez-Carmona et al., 2006; Ibanez-Mejia et al., 2011). The western Cordillera (“w” on Fig. 1) comprises oceanic magmatic and deep-marine sedimentary rocks of mainly Cretaceous age (McCourt et al., 1984; Cortés et al., 2005; Kennan and Pindell, 2009; Borrero et al., 2012). The main volcanic arc along the central Cordillera (“c” on Fig. 1) is composed mainly of Jurassic–Cretaceous igneous rocks intruded into and resting on both oceanic and continental crystalline basement (Aspden and McCourt, 1986; Aspden et al., 1987; Cediel et al., 2003). The eastern Cordillera (“e” on Fig. 1) consists of Archean–Mesoproterozoic inliers surrounded by 1.1–1.0 Ga rocks (Restrepo-Pace et al., 1997; Ruiz et al., 1999; Cordani et al., 2005; Jiménez-Mejía et al., 2006; Ordóñez-Carmona et al., 2006; Chew et al., 2007a, 2007b, 2008; Cardona et al., 2009; Miskovic et al., 2009; Ibanez-Mejia et al., 2011).

Central Andes

The central Andes (CA) consist of arc-type and basinal assemblages that formed in a convergent-margin setting from early Paleozoic time to the recent. Some of these assemblages may be exotic to South America (e.g., Cuyania terrane; Astini and Thomas, 1999; Ramos, 2004; Thomas et al., 2004; Finney, 2007), whereas most formed during long-lived accretionary tectonism along the proto-Andean margin (Cawood, 2005). Some terranes contain crystalline basement of Paleo- or Meso-Proterozoic age (Shackleton et al., 1979; Loewy et al., 2004; Chew et al., 2007a) or of Neoproterozoic age (Jordan and Allmendinger, 1986; Strecker et al., 1989; Allmendinger et al., 1997; Escayola et al., 2007).

Southern Andes

The southern Andes (SA) consist primarily of Proterozoic and Paleozoic rocks exposed in basement uplifts (Pankhurst et al., 2000, 2006; Augustsson et al., 2006), rimmed to the west by post-Triassic magmatic rocks of the Andean orogen (Kay et al., 2005; Hervé et al., 2007). The southern portion of SA belongs to the Patagonia microcontinent, which consists of widespread Paleozoic magmatic arc and basinal assemblages that formed in the broad accretionary Terra Australis orogen (Coney, 1992; Stern and De Wit, 2003; Cawood and Buchan, 2007; Ramos, 2008; Barbeau et al., 2009).

PRESENT STUDY

Sampling Methodology

Samples were collected from active channels of modern rivers and from beaches along shorelines. Most of the rivers sampled do not have dams or major construction projects that could have influenced zircon age distributions. Each sample consisted of a mixture of sediment from multiples sites (within a ∼10 m2 area) with a wide range of grain sizes (silt to coarse sand). Material coarser than 2 mm was removed during sample collection. This strategy was followed in an effort to maximize the number of age groups recognized and reduce biases from hydrodynamic sorting (Garzanti et al., 2009; Lawrence et al., 2011). Heavy minerals from each composite sample were initially concentrated with a gold pan, and then zircons were extracted using a Frantz isodynamic separator and heavy liquids (techniques described by Gehrels, 2000).

Sample Descriptions

Fifty-nine samples were collected from rivers draining the northern, central, and southern Andes and from beach deposits near river deltas (Fig. 1 and Table 1).

Northern Andes

Thirteen samples were collected from the northern Andes (NA) as follows. Sample NA01 was collected from the Caribbean coast. Seven samples were collected from the Magdalena River system in the Magdalena Valley (separating the eastern and central Cordilleras of Fig. 1), with NA02 from lowlands 140 km upstream, NA03 and NA04 from the Magdalena delta, NA05 along the coast 70 km west of the Magdalena delta, and NA06–NA09 from the flanks of the eastern and western Cordillera (EC and WC) in the Magdalena highlands. Four samples were collected farther south, with NA10 from the Andean highlands, NA11 from a tributary of the Amazon River system, and NA12 and NA13 along the western flank of the western Cordillera (“w” on Fig. 1).

Central Andes

Twenty-five samples were collected from both the western and eastern flanks of the central Andes (CA), with a spacing of ∼300–400 km. Along the eastern flank of the Andes, CA01–CA04 were collected from tributaries of the Rio de la Plata river system, CA05 and CA06 are from tributaries of the Amazon river system, CA07 and CA08 are from the Andean highlands, and CA09 and CA10 are from southern Amazon tributaries. Along the western flank of the Andes, samples CA11–CA22 were collected from small perennial rivers near the Atacama Desert. CA16 and CA23 are from beach sands to the west.

Southern Andes

Five southern Andes (SA) samples (SA01–SA05) were collected from the western flank of the orogen from rivers with a spacing of ∼300 km. Samples SA01 and SA02 are from perennial rivers within 10 km of the coast, whereas sample SA03 was collected from the coast and SA04 was collected farther upstream. Samples SA06–SA22 are from the eastern flank of the Andes, with a spacing of ∼200 km. These samples are from lowland portions of rivers, ∼30–50 km from the Atlantic coast.

U-Pb Analytical Methods

Uranium-lead geochronology of zircons was conducted by laser ablation–multicollector–inductively coupled plasma mass spectrometry (LA-MC-ICPMS) at the Arizona LaserChron Center (www.laserchron.org) using analytical methods described by Gehrels (2000), Gehrels et al. (2006, 2008), and Gehrels and Pecha (2014). Cathodoluminescence (CL) images, acquired with a Gatan Chroma CL2 system (www.geoarizonasem.org), were used to ensure that analyses were conducted on homogeneous (generally interior) portions of zircon grains.

The U-Pb analytical data are reported in Supplemental Data 1–31,2,3. These tables present the isotopic ratios, concentrations, calculated ages, and uncertainties for each analysis. Uncertainties are reported at 1-sigma, and include only internal uncertainty components (measurement of 206Pb*/238U, 206Pb*/207Pb*, 206Pb*/204Pb) and a factor to incorporate the overdispersion of primary standard analyses. External (systematic) uncertainties are reported for each sample (at 2-sigma) in Supplemental Data 1–3 (see footnotes 1–3) and include contributions from the scatter of primary standard analyses, composition of common Pb, age of the primary standard, and U decay constants. Also reported for each analysis is the “best age,” which is based on 206Pb*/238U for <900 Ma grains and 206Pb*/207Pb* for >900 Ma grains.

Following standard practice (e.g., Gehrels, 2012, 2014), the U-Pb data tables do not include analyses that are highly (>20%) discordant or highly (>10%) reverse discordant and analyses with >10% uncertainty. These filters are not applied to <600 Ma 206Pb*/207Pb* ages and <100 Ma 206Pb*/238U ages. Also included in Supplemental Data 1–3 (see footnotes 1–3) are Pb*/U concordia diagrams (from Ludwig, 2008) and age-distribution diagrams that sum all ages (and associated uncertainties) from a sample (or set of samples) into a single curve (from www.laserchron.org).

For interpretive purposes, U/Th ratios for each analysis are reported in Supplemental Data 1–3 (see footnotes 1–3). Most (>90%) analyses have U/Th in the range of typical igneous values (e.g., 1–10; Rubatto, 2002), which are interpreted to record episodes of igneous activity. Values >10 are interpreted as possibly of metamorphic origin (e.g., Gehrels et al., 2009).

Hf Analytical Methods

Hafnium analyses were also conducted by LA-MC-ICPMS at the Arizona LaserChron Center, using methods described by Cecil et al. (2011) and Gehrels and Pecha (2014). Laser ablation analyses were conducted with a laser beam diameter of 40 microns, with ablation pits located either on top of the U-Pb analysis pits or adjacent to the U-Pb pits but in the same domain (using CL images as guides). Each acquisition consisted of one 40-second integration on backgrounds (on peaks with no laser firing) followed by 60 one-second integrations with the laser firing. The isotope ratios and uncertainties reported are based on all measurements, with only a 2-sigma filter to remove outliers. Critical isotope ratios used in data reduction are as follows: 179Hf/177Hf = 0.73250 (Patchett and Tatsumoto, 1980); 173Yb/171Yb = 1.132338 (Vervoort et al., 2004); 176Yb/171Yb = 0.901691 (Vervoort et al., 2004; Amelin and Davis, 2005); 176Lu/175Lu = 0.02653 (Patchett, 1983); and the decay constant of 176Lu = 1.867e-11 (Scherer et al., 2001; Söderlund et al., 2004). Hafnium-depleted mantle model ages are not calculated because of uncertainties in the 176Hf/177Hf and 176Lu/177Hf of source material(s) from which the zircons crystallized.

Analyses were calibrated using solution analyses of JMC475, solutions of JMC475 mixed with varying concentrations of 176Lu and 176Yb, and laser ablation analyses of standard zircons Mud Tank, 91500, Temora2, R33, FC1, and Plesovice. Solutions of JMC475 yielded a weighted mean value of 0.282156 ± 0.000007 (2σ, n = 17). Solutions mixed with 176Yb and 176Lu [with 176(Yb + Lu)/176Hf up to 1.6] yielded a similar weighted mean value of 0.282150 ± 0.000008 (2σ, n = 16) with 173Yb/171Yb adjusted to 1.132423 (from value of 1.132338 suggested by Vervoort et al., 2004) to alleviate a slight correlation between measured 176Hf/177Hf and 176(Yb + Lu)/176Hf.

Zircon standards yielded measured 176Hf/177Hf values that were slightly lower than the reported values. With measured values increased by 0.8 epsilon units, the following results were obtained:

  • Mud Tank: weighted mean value of 0.282518 ± 0.000007 (2σ, n = 138) compared with reported value of 0.282507 (Woodhead and Hergt, 2005). This yields an offset of +0.4 epsilon units.

  • 91500: weighted mean value of 0.282292 ± 0.000021 (2σ, n = 33) compared with reported value of 0.282313 (Fisher et al., 2014). This yields an offset of –0.7 epsilon units.

  • Temora2: weighted mean value of 0.282677 ± 0.000009 (2σ, n = 33) compared with reported value of 0.282686 (Fisher et al., 2014). This yields an offset of –0.3 epsilon units.

  • Plesovice: weighted mean value of 0.282480 ± 0.000007 (2σ, n = 125) compared with reported value of 0.282484 (Slama et al., 2008). This yields an offset of –0.1 epsilon units.

  • FC1: weighted mean value of 0.282170 ± 0.000009 (2σ, n = 132) compared with reported value of 0.282183 (Fisher et al., 2014). This yields an offset of –0.4 epsilon units.

  • R33: weighted mean value of 0.282735 ± 0.000009 (2σ, n = 148) compared with reported value of 0.282764 (Fisher et al., 2014). This yields an offset of –1.0 epsilon units.

These offsets suggest that the accuracy (external reproducibility) of a set of analyses is 1–2 epsilon units. The precision (internal reproducibility) of individual analyses is ∼3 epsilon units, as indicated by an average measurement uncertainty of 2.8 epsilon units (average standard deviation of all standard measurements, at 2σ), and an average standard deviation for measurements of each standard of 2.8 epsilon units (at 2σ).

Data from all unknowns have been adjusted using a correction factor of +0.8 epsilon units and the revised 173Yb/171Yb value reported above.

Statistical Analyses

Three different statistical methods have been applied to analyze patterns and potential correlations among the geochronologic, isotopic, and plate motion data presented herein.

In an effort to characterize patterns of magmatism, we use a maximum likelihood analysis (from Unmix routine of Isoplot; Ludwig, 2008) to characterize the age and variation of overlapping ages. This method fits a Gaussian distribution to a set of ages and returns the mean and standard deviation of the constituent ages. The standard deviation is expressed at the 2σ level to characterize the range of included ages.

To investigate the cyclicity of magmatic records and plate convergence rates, we apply a Lomb-Scargle least-squares spectral analysis (Lomb, 1976; Scargle, 1982) to the geographically organized age distributions and the relative plate convergence rates of Maloney et al. (2013). Although the Lomb-Scargle least-squares spectral analysis is typically used to estimate periodicity of irregularly sampled data, the Lomb-Scargle method has some advantages (e.g., data vector can be irregularly sampled, data vector can have any length, and periodograms have known statistical properties) over the more common Fourier power spectral density (PSD) analysis. We note that although we prefer the Lomb-Scargle analysis, our data sets have a regular sampling interval (1 m.y.), and a Fourier PSD analysis produces a similar estimate of the PSD. Since we are comparing time series that measure different physical properties, we are not interested in the absolute power of each time series but in the relative peaks in periodicity between time series. To highlight any potential cyclical correlations between these time series, we normalize each PSD estimate setting the maximum power observed for each time series equal to one.

We also explore potential correlations between plate convergence rates and the Andean magmatic history by first employing a sophisticated continuous wavelet transform (CWT) analysis to each of the individual records. Traditional Fourier analysis implicitly assumes any physical relationships are stationary in time. The wavelet method has several advantages over traditional Fourier analysis and is particularly well suited for analyzing geological and geophysical time series that may exhibit periodic behavior that varies in time with localized intermittent periodicities (Torrence and Compo, 1998). For each individual time series, we calculate the continuous (local) wavelet power spectrum. Because the wavelet is not everywhere localized in time, we define a cone of influence (COI), outside of which edge effects are significant. We also define a ≥95% confidence interval using a first-order autoregressive process (Allen and Smith, 1996) that models red noise (random walk noise) to identify statistically significant areas of signal-generated wavelet power. The continuous (local) wavelet power spectrum results for each individual time series show signal power at similar periods as the Lomb-Scargle periodograms and also contain additional information about localized intermittent periodicities.

The CWT analyses described above provide statistically significant information about localized intermittent periodicities. To explore potential mechanistic explanations of the magmatic records from South America, we employ two approaches of analyzing pairs of CWTs generated from individual time series—the cross wavelet transform (XWT) and the wavelet coherence (WTC). Here we briefly describe these two approaches and refer the interested reader to Grinsted et al. (2004) for a more detailed formal mathematical description. The XWT provides a statistically significant measure (≥95% confidence interval against red noise) of the common power between two CWTs as a function of period and time. Both the XWT and WTC have the additional property of returning the phase difference between signals with common power and/or coherence, also as a function of period and time. This information is helpful in establishing a mechanistic connection as any cause and effect relationship would be phase locked, meaning the phase (time lag) would be constant for any period with statistically significant common power and/or coherence. The WTC extends the analysis of XWT to time series with significant coherence even if the common power is low. The WTC is analogous to traditional correlation coefficients with the added property of quantifying correlation within localized intermittent periodicities.

RESULTS

U-Pb Geochronology

A total of 5524 new U-Pb analyses are reported in Supplemental Data 1–3 (see footnotes 1–3). Figure 1 shows the location of each sample on a geologic map of South America, and Supplemental Item 14 provides a description of the zircon grains present in each sample. Highlights of the results from samples from northern, central, and southern regions are described in following sections and shown on Figures 2–4. Figures 5 and 6 summarize the U-Pb age distributions from each region, with all age distributions smoothed using a sliding window that averages each set of five closest ages.

Northern Andes

Samples from the northern Andes (NA) contain detrital zircons that are quite variable in age (Fig. 2). Most samples contain a few ages between ca. 2.1 and ca. 1.7 Ga. The only age group present in most samples ranges from ca. 1.6 Ga to ca. 950 Ma, with dominant ages near ca. 1.0 Ga. Several samples contain grains of 700–400 Ma, and most samples contain abundant <400 Ma grains. The dominant ages for <200 Ma grains are shown on Figures 5 and 6.

Central Andes

For >200 Ma ages, the age patterns in the central Andes (CA; Fig. 3) generally resemble age patterns from regions to the north (Fig. 2), although proportions of the main age groups are somewhat different (Figs. 5 and 6). Namely, most samples are missing the ca. 1.6 Ga to ca. 950 Ma and 400–200 Ma age groups. Similarly, for <200 Ma ages, some age groups (e.g., those with peak ages of ca. 183 and ca. 126 Ma) are also missing from CA, and the proportions of other age groups are somewhat different (Figs. 5 and 6). The 600–400 Ma age group is much more prominent in the CA than the NA.

Southern Andes

The southern Andes (SA) yield age distributions that are quite different from NA and CA, with predominance of 400–250 Ma and <25 Ma grains in most samples and abundant 200–25 Ma ages in only three samples (Figs. 4–6). Southern Andes samples contain very few ages between 900 and 700 Ma and >1.2 Ga.

East-West Comparison of U-Pb Ages

As a complement to the above north-south comparison of ages, Figures 5 and 6 compare age distributions from rivers that are located along the eastern versus western flanks of the modern Andes. Omitted from this compilation are samples from rivers in the northernmost Andes (which drain northward) and in southern Patagonia (which have no sampled westerly counterparts). This comparison highlights the dominance of older basement sources in rivers of the eastern Andes versus the prevalence of Phanerozoic igneous rocks in the western Andes. This compilation also records a strong cyclicity to all records of Andean magmatism, which is explored in detail below.

Hf Isotope Geochemistry

The results of our Hf isotope analyses are described in the following sections for samples from the northern, central, and southern Andes. A total of 1199 new analyses are reported in Supplemental Data 4–65,6,7. These tables report all details of individual analyses and present plots with results from each sample. Results are summarized on Figures 2–4.

Northern Andes

Hafnium isotope data from the northern Andes are quite variable and record interesting patterns through time. As shown on Figure 2, εHft values progressively decrease (become more evolved) from ca. 1.6 Ga through ca. 400 Ma, increase (become more juvenile) from ca. 300 to ca. 100 Ma, and then in most cases become more evolved again after ca. 100 Ma (Fig. 2).

Central Andes

Analyses from the central Andes reveal a generally similar pattern of changing εHft values (Fig. 3). Additional patterns include the presence of more abundant >1.6 Ga grains with variable εHft values, more juvenile values for 1.2–1.0 Ga grains, a distinctive pull-down of 650–500 Ma values, and considerably more negative εHft values for the youngest grains.

Southern Andes

The southern Andes (Fig. 4) yield few >1.2 Ga grains, mostly juvenile 1.2–1.0 Ga grains, abundant 600–300 Ma grains with intermediate εHft values, values that increase from moderately negative to positive from 300 to 100 Ma, and then mostly positive values for <100 Ma grains.

COMPARISON WITH EXISTING U-Pb AND Hf DATA FROM SOUTH AMERICA

U-Pb

In an attempt to compare our U-Pb geochronologic data set with previously published information, we have compiled ages reported from igneous and detrital samples from all of western South America (Figs. 5 and 6). All available U-Pb ages are included in the compilation, whereas only young (<200 Ma) volcanic K-Ar and Ar-Ar ages are included because of the likelihood that pre–200 Ma ages record cooling rather than igneous activity. Ages of individual analyses are included where reported; mean ages are reported where results from individual analyses are not available. Supplemental Data 78 includes all compiled ages and also a map that shows the distribution of sampled areas. The sources for the compiled data are listed in Supplemental Item 29. A total of 26,851 detrital ages and 9929 igneous ages are included in this compilation, in addition to our 5524 new ages.

Comparisons between our detrital data, previously published detrital data from both sedimentary rocks and modern sediment, and igneous rocks are provided on Figure 6. This comparison shows that our data and previously reported detrital data are quite similar in the dominant ages recognized and in the approximate proportions of different age groups. These two data sets are somewhat different from the age distribution defined by igneous samples, presumably because of temporal and/or spatial biases in the selection of studied igneous suites.

Included in this compilation are analyses of detrital-zircon ages from several South American rivers. Condie et al. (2005, 2009) report ages from rivers in Brazil and Peru, whereas Iizuka et al. (2010) and Rino et al. (2004, 2008) report ages from the Amazon and Parana rivers. Mapes et al. (2009) have also conducted analyses from Amazon River sediment, with analyses that are forthcoming. Campbell and Allen (2008) have included analyses from South American rivers in their database, but these ages have not been included in our compilation because sample locations were not reported.

Hf

Hafnium data were also compiled from previous studies that reported essential information about the analyses (e.g., instrumentation, methodology, and data reduction). The acceptable data are tabulated in Supplemental Data 810, with associated references reported in Supplemental Item 311. No acceptable analyses are available for the northern Andes; 688 were found from the central Andes and 85 from the southern Andes.

INTERPRETATION OF RESULTS

Our new and compiled data provide significant new insights into the magmatic history and crustal evolution of western South America. Each of these aspects is described below.

Magmatic History

The U-Pb data presented herein, combined with existing geochronologic information (Supplemental Data 7 [see footnote 8]; Supplemental Item 2 [see footnote 9]) and interpretations (e.g., Bahlburg and Hervé, 1997; Baldwin et al., 1997; Haschke et al., 2002, 2006; Cediel et al., 2003; Franz et al., 2006; Pankhurst et al., 2006; DeCelles et al., 2009, 2015; Miskovic et al., 2009; Ramos, 2009; Ducea et al., 2015; Paterson and Ducea, 2015) provide important new information about the history of magmatism in western South America. This record is most relevant for igneous activity of felsic to intermediate composition, given that most of the information is derived from ages of zircon, and zircon is generally of low abundance in mafic magmas. As noted above, the connection between age distribution and magmatic addition is also complicated by potential variations in Zr content (and hence zircon fertility) through space and time, the long and complex history of uplift and erosion in the Andes, and potential additions to and removals from the system by tectonic and erosional processes.

Given the objective of reconstructing the magmatic history of western South America, with emphasis on both ages and proportions of magmatism, it is important to recognize that all three sets of information (our data, published detrital data, and published igneous data) have biases. For igneous ages, it is clear that the age distribution is controlled in large part by the geographic and temporal focus of the constituent studies. This is apparent in the igneous age distribution (Figs. 5 and 6) by the large proportion of young igneous ages, which reflects the numerous studies of young Andean magmatism. In contrast, previously published detrital ages from sandstones are biased away from these young ages, given that detrital zircons in sandstones always predate deposition. The detrital ages reported herein are biased toward contributions from the modern Andes, given that most large rivers in South America have headwaters in the Andes and that most of our samples were collected from within or adjacent to the Andes. It is also important to note that detrital age distributions may overrepresent older sources given the tendency for zircon grains to be recycled or underrepresent older sources if older grains are too small to be analyzed or are compromised by Pb loss (Gehrels, 2014). For these reasons, the following discussions focus on general patterns of ages, with the assumption that short-term variations in the proportion of ages may reflect true variations in magmatic addition, whereas variations over longer time periods may be controlled mainly by geologic and/or sampling biases.

Figure 5 shows that there are notable similarities in all three records. Most impressive is the similarity in detrital data sets from this study and previous studies in terms of both presence or absence of ages and proportions of ages. Igneous ages are also similar for most areas and time periods but, as noted above, emphasize young Andean magmatism.

One of the first-order patterns apparent in all data sets is the relatively low abundance of pre–1.2 Ga ages, in spite of the fact that much of South America is underlain by Archean and Paleoproterozoic cratonal rocks (Fig. 1). This low abundance is due in part to biases in sampling (as noted above) and also documents the scarcity of Archean–Paleoproterozoic crystalline rocks in the modern Andes.

The oldest dominant age group includes 1.2–1.0 Ga grains, at least for the central and northern Andes (Fig. 5). These 1.2–1.0 Ga grains presumably record magmatism within the Sunsas belt (Fig. 1), which is exposed along the eastern flank of the northern and central Andes. Grains of these ages may also have been distributed widely across South America during Mesoproterozoic tectonism, as has been documented for sediment derived from the Grenville orogen in North America (e.g., Rainbird et al., 2012).

The next main phase of magmatism recorded in the three data sets begins at ca. 700 Ma and extends into early Paleozoic time (Fig. 5). This magmatism occurred within the Brasiliano belt (0.7–0.5 Ga) during assembly of the South American craton, which coincides with the Pan-African orogeny and the assembly of Gondwana (Windley, 1995; Brito Neves et al., 1999). This phase of continental assembly overlapped with the late Neoproterozoic onset of magmatism along the western margin of South America in a broad and long-lived accretionary orogen (Coney, 1992; Bahlburg and Hervé, 1997; Cawood, 2005; Cawood and Buchan, 2007). Igneous geochronology records magmatism along this convergent margin from ca. 610 Ma until ca. 470 Ma (Cawood and Buchan, 2007), which is also reflected in the record of detrital-zircon ages (Fig. 5).

Our three geochronologic records all show a decrease in magmatism during mid-Paleozoic time (ca. 400 to ca. 360 Ma), followed by significant magmatic addition between 360 and 200 Ma (Fig. 5). This pattern is quite similar to the magmatic history described by Bahlburg and Hervé (1997) and extended through the proto-Andes and to other portions of the Terra Australis orogen by Cawood and Buchan (2007).

The data provide a rich record of variations in magmatism during the past ∼200 m.y. As summarized on Figure 6, the record is somewhat different from east to west, from north to south, and in the three different data sets. As a starting point for discussion, we interpret the maxima in ages generated during this study (“This study (all)” on Fig. 6) to indicate times of higher magmatic addition, and we use maximum likelihood analysis to determine midpoints, ranges, and relative proportions (in percent) of magmatism at 183 ± 11 (9% of <200 Ma grains), 166 ± 10 (5%), 149 ± 9 (6%), 125 ± 18 (5%), 110 ± 9 (8%), 88 ± 7 (14%), 65 ± 6 (13%), 35 ± 5 (7%), 21 ± 6 (17%), and 4 ± 4 (16%) Ma (2σ). The shaded gray bars of Figure 6 correspond to these age ranges and midpoints and allow comparison of this interpreted record with other published age distributions and with specific records for NA, CA, SA, and western versus eastern portions of the Andes.

Temporal variations in age abundance in each of the different records shown in Figure 6 suggest that Andean magmatism is strongly cyclic. The period of this cyclicity can be interpreted from the times of higher magmatic addition noted above, which yields an average of ∼20 m.y. (10 maxima during 200 m.y.). The periodicity of this cyclicity is quantified with a Lomb-Scargle least-squares spectral analysis. As shown in Figure 7, the record of all previously published and new analyses has the highest signal power (dominant frequency) at ∼33 m.y., whereas each of the separate age records has the highest signal power between ∼20 and ∼34 m.y. Plate convergence velocities along the Andean margin have the highest power at ∼20 m.y. Possibilities for the origin of the apparent cyclicity in magmatism are explored below.

Crustal Evolution

The Hf data presented herein, combined with previously published Hf isotope analyses, provide new constraints on the crustal evolution of western South America. All of the available data are shown for NA, CA, and SA on Figures 8–10, and together for all areas on Figure 11. Accompanying these Hf isotope plots are age-distribution diagrams for our new data as well as for all available (igneous and detrital) data.

Yellow arrows on these diagrams show interpreted trajectories of Hf evolution (based on least-squares regression of relevant data). Arrows with slopes that are subparallel to average crustal evolution (black arrows) are interpreted to represent successive remelting of a relatively uniform suite of crustal materials. Arrays that project toward more positive (or less negative) values than average crust (black arrows) are interpreted to record a progressive increase in the proportion of juvenile magmas or slightly older juvenile crust. Arrays that project toward less positive (more negative) values than average crust (black arrows) may result from increased incorporation of old crustal material, involvement of increasingly negative crustal material, and/or reduced involvement of juvenile magmas or crust. Vertical arrays are interpreted to result from melting of crustal materials with highly variable epsilon Hf, along with possible incorporation of juvenile magmas.

Using this interpretive framework as a guide, the available U-Pb/Hf data support the following history of crustal genesis and reworking for western South America:

  1. Grains that are >1.6 Ga in age are scarce (especially in NA), and εHft values are highly variable, but there is a general trend toward more negative εHft values with time. This is interpreted to record Paleoproterozoic reworking of mostly Archean crustal materials. There are also numerous grains (especially in CA) with intermediate εHft values that include Paleoproterozoic juvenile material. This history pertains mainly to the cratonal rocks of South America (Fig. 1).

  2. Between ca. 1.6 Ga and ca. 1.0 Ga, εHft values in CA and possibly SA tend to become more positive, reflecting progressive incorporation of juvenile magmas. In contrast, εHft values in NA remain intermediate in value, reflecting lesser addition of juvenile crustal material. Given the proximity of the Sunsas belt with NA, the latter history may be relevant for the Mesoproterozoic crustal evolution of northwestern South America. In contrast, the more juvenile Mesoproterozoic zircons in CA and SA may have been shed from other regions of Rodinia—for example, North America, which yields positive εHft values for grains of this age (Gehrels and Pecha, 2014). Distribution of these zircons may have occurred as part of a global sediment dispersal system during assembly of Rodinia (e.g., Rainbird et al., 2012) or less likely, from early Paleozoic arrival of North American tectonic fragments with 1.2–1.0 Ga basement (e.g., Cuyania) along the Andean margin (e.g., Astini and Thomas, 1999; Ramos, 2004; Thomas et al., 2004).

  3. Between ca. 1.0 Ga and ca. 500 Ma, εHft values have trajectories that are generally similar to average crustal evolution (Figs. 8–11). This is interpreted to reflect mainly reworking of >1.0 Ga crustal materials, with little addition of juvenile magma. Franz et al. (2006), Willner et al. (2008), and Bahlburg et al. (2009, 2011) presented a similar conclusion on the basis of Sr, Nd, and Hf isotope signatures of igneous and sedimentary rocks from the central Andes. In addition to this overall trend, there are apparent pull-downs (negative excursions) in εHft at 650–550 and ca. 500 Ma (Fig. 11). The older of these excursions likely reflects reworking of older crustal components in Brasiliano orogens, whereas the younger excursion may record involvement of older crustal components within early Paleozoic arcs along the proto-Andean margin (blue dashed lines on Fig. 1). It is interesting that other portions of the broad Terra Australis orogen are also characterized by significant reworking of older crustal materials during this time period (e.g., Cawood et al., 2009, 2012; Kemp et al., 2009; Hawkesworth et al., 2010).

  4. Early Paleozoic (500–400 Ma) magmatism yields εHft values that remain intermediate or become slightly more positive in the CA and SA. This positive deviation from average crustal evolution suggests that the widespread early Paleozoic magmatism along the Andean margin (e.g., Rapela et al., 1992; Bahlburg and Hervé, 1997; Pankhurst et al., 1998; Coira et al., 1999; Ducea et al., 2010; Einhorn et al., 2015) may have occurred in an overall extensional tectonic regime with thinned continental crust and significant additions of mantle-derived melt (e.g., Ducea et al., 2015).

  5. Late Paleozoic (ca. 360 to ca. 250 Ma) magmatism is treated separately from early Paleozoic igneous activity because of the 400–360 m.y. magmatic lull and a significant change in Hf isotope evolution. As shown on Figure 11, there is an overall trend toward more negative εHft values during this time, reflecting increased involvement of older crustal materials in late Paleozoic magmatic arcs established along the Andean margin (green dashed lines on Fig. 1). This trend is apparent in NA (Fig. 8) and SA (Fig. 10) but less obvious in the data from CA (Fig. 9). Involvement of continental material may have increased as the late Paleozoic arcs approached and were firmly accreted to the proto-Andean margin (e.g., Willner et al., 2008; Bahlburg et al., 2009; Ramos, 2009).

  6. Between ca. 250 and ca. 100 Ma, εHft values become more positive in all areas (Figs. 8–11) and mirror the presence of juvenile Sr-Nd-Pb isotopes in rocks of the central Andes (Franz et al., 2006) and epsilon Nd values in Patagonia (Hervé et al., 2007). This trend is interpreted to result in large part from Triassic through Early Cretaceous extensional tectonism in the subduction system that is well documented all along the Andean margin (Fig. 12; Franz et al., 2006; Ramos, 2009; Maloney et al., 2013; Rossel et al., 2013). Such extension is also manifest as widespread rift-related strata and normal faults of Triassic and Jurassic age (Ramos, 2009), formation of the Late Jurassic Rocas Verdes backarc basin in the southern Andes (Dalziel et al., 1974), formation of the Salta rift system during Late Jurassic–Early Cretaceous time (Salfity and Marquillas, 1994), and opening of basalt-floored basins during Late Jurassic–Early Cretaceous time in the northern Andes (Maloney et al., 2013). The observed trend toward more positive εHft values may have also resulted from migration of the magmatic arc westward. As reported by Haschke et al. (2002, 2006) and Franz et al. (2006), the Late Jurassic–Early Cretaceous magmatic arc was situated along the far-western margin of the orogen, which is underlain by Paleozoic arc terranes containing little evolved continental material (Fig. 1).

  7. Between ca. 100 and ca. 35 Ma, most zircons lie along trends toward more negative εHft values that are steeper than average crustal evolution (Figs. 8–11). This is interpreted to reflect a progressive decrease in the proportion of juvenile magma generated within the arc systems and/or reworking of progressively older crustal materials. The former may be related to a change in subduction dynamics caused by an increase in the age of the slab being subducted and/or an increase in the convergence rate between South America and outboard ocean basins (e.g., Maloney et al., 2013). It is also likely that the average age of reworked crustal materials increased with time as the magmatic arc migrated back toward cratonal South America (e.g., Haschke et al., 2002, 2006; Franz et al., 2006). Although most zircons lie along this apparent trend, at least in NA and CA, it is important to note that many <100 Ma zircons (especially in SA) yield highly positive εHft values. Variations in the proportion of juvenile versus reworked evolved materials are not surprising given the length and complexity of the Andean magmatic arc system.

  8. Zircons that are younger than ca. 35 Ma yield highly variable εHft values, ranging from –15 to +15 (Fig. 11), which is similar to the broad range of epsilon Nd values in Patagonia (Pankhurst et al., 1999; Hervé et al., 2007). This variation is interpreted to reflect the broad range of tectonic settings along the Andean magmatic arc, from areas of crustal thickening where magmas are forming due to crustal melting, to areas of extension where mantle-derived melts are present (e.g., Kay et al., 2005; DeCelles et al., 2009, 2015; Ramos, 2009).

DISCUSSION

Magmatic Processes

There are relatively consistent patterns of magmatism and crustal evolution in all of the geochronologic and Hf isotopic data sets available for South America. In terms of ages of felsic to intermediate-composition (e.g., zircon-bearing) magmas, the history is apparently episodic, with long-term patterns that parallel the first-order global record (Fig. 12; Voice et al., 2011; Condie, 2014). Significant mismatches are the abundance of <200 Ma ages and the scarcity of >1.4 Ga ages in our South American record, which is not surprising given that most of our samples were collected in proximity to the Andean orogen or from rivers with headwaters in the Andes. This apparent episodicity, with a period of ca. 0.6 Ga, is widely ascribed to aspects of crustal generation and/or preservation associated with the supercontinent cycle (Hawkesworth et al., 2009, 2010; Belousova et al., 2010; Condie et al., 2011; Voice et al., 2011; Condie, 2014), with periods of low zircon production and/or preservation related to breakup of supercontinents (Fig. 12).

As noted above, there is also a strong episodicity in the ages of magmatism during the past 200 m.y. (Figs. 6 and 7). Generally similar times of high and low magmatism are apparent in all data sets, with age maxima at ca. 183, 166, 149, 125, 110, 88, 65, 35, 21, and 4 Ma (shaded gray bars on Fig. 6). There are significant variations in these records, however, with peak ages that are younger, older, or missing entirely along some transects, and the first-order periodicity ranges from ∼20 to ∼34 m.y. (average of ∼33 m.y.) depending on how the data are grouped (Fig. 7).

Possible causes of cyclicity in convergent margin systems have been explored by Davies and Stevenson (1992), Elliott et al. (1997), Haschke et al. (2002, 2006), DeCelles et al. (2009, 2015), Lee and King (2011), Gibert et al. (2012), Ramos et al. (2014), Ducea et al. (2015), and Paterson and Ducea (2015). Explanations generally focus on aspects of the downgoing plate (e.g., convergence velocity, age of downgoing slab, angle of downgoing slab, and dehydration providing fluids for the overlying mantle wedge) versus upper-plate controls such as tectonic thickening and/or thinning and processes of delamination. These recent syntheses have been constructed in large part on the basis of studies by Haschke et al. (2002, 2006) of magmatism in the central Andes, as well as magmatism in western North America.

Haschke et al. (2002, 2006) reported that the central Andes have experienced four phases of magmatism at 195–129, 123–82, 76–37, and 22–0 Ma, yielding an average periodicity of ∼50 m.y. (purple fields on Fig. 13). Magmatic gaps are reported to exist at 129–123, 82–76, and 37–22 Ma. Haschke et al. (2002) suggest that a correlation between magma addition and convergence rate may be due to a cycle of decreasing slab dip (periods of high magma addition) followed by slab breakoff or rollback (periods of low magma addition). Ramos (2009) and Ramos et al. (2014) propose an Andean orogenic and/or magmatic cycle that invokes changes in slab dip and delamination of lower crustal roots as a primary control on types and volumes of magmatism. DeCelles et al. (2009, 2015) provide an alternative explanation of cyclicity that involves periods of compressional tectonism, heating, hydration, eclogitization, and delamination of the upper plate, with less reliance on processes directly associated with the subducting plate (e.g., slab angle and subduction velocity).

Possible tectonic drivers for the cyclicity apparent from our data set are explored on Figure 14, which compares our U-Pb ages and Hf isotope values with convergence velocities (blue lines), rates of absolute westward motion of South America (green lines), and ages of the downgoing slab (purple lines) (from Maloney et al., 2013) for each of our three transects. Of the three parameters, there is little visual evidence of a correlation between magmatic history and either the age of the downgoing plate (purple lines) or the absolute westward motion of South America (green lines). Possible connections between magmatic history and convergence rate between South America and subducting oceanic plates (blue lines) are more difficult to evaluate visually because of the complexity of the various records.

Figure 15 shows the results of our attempts to statistically evaluate possible linkages between magmatism and convergence velocity using continuous wavelet transform (CWT) analysis for the central Andean section (similar results for SA and NA are shown in Supplemental Item 412). Although there appears to be some statistically significant common signal power with a periodicity between 16 and 24 m.y. (Figs. 15C–15E), the signals are inconsistently out of phase (not phase locked) across these periods (e.g., black arrows change direction for same period in Fig. 15E). This indicates that the time lag between the two signals with common power is changing with time. If there were a simple cause and effect relationship between magmatism and plate convergence rates (e.g., a plate convergence rate increases X m.y. before an age distribution increases), we would expect the phases to be locked (e.g., all black arrows pointing in the same direction for the same period in Fig. 15E). This interpretation is further supported by the lack of coherence between the two time series, where only two small areas centered at 20 and 40 Ma (period of ∼5 m.y.) show any statistically significant correlation (Fig. 15F). A similar lack of statistically significant correlation is observed for the southern and northern Andean sections (Supplemental Item 4 [see footnote 12]). Additionally, spectral analysis of the convergence velocities (from Maloney et al., 2013) yields highest signal power at ∼20 m.y. for all three regions (Fig. 7), which is somewhat shorter than the observed magmatic periodicities. We therefore conclude that the observed magmatic periodicities are not attributable to a simple genetic connection with plate convergence rates.

Crustal Generation and Reworking

As with magmatic processes, our data set provides new insights into processes of crustal generation, reworking, and preservation that operate at two different time scales. On a broader time scale, our data provide an interesting comparison with the Hf isotopic results presented by Belousova et al. (2010), Collins et al. (2011), Condie et al. (2011), and Cawood et al. (2012) (Fig. 12). Our data resemble the global patterns prior to ca. 800 Ma but are somewhat different for younger time periods. The first significant difference is the abundance of juvenile values from 800 to 650 Ma in the Belousova et al. (2010) compilation, which is not seen in South America. Next is the negative excursion from 650 to 500 Ma in the global data sets (due largely to Pan-African tectonism), which is not seen in our data set. A final difference is the abundance of juvenile values from 400 to 200 Ma in the Belousova et al. (2010) compilation (mainly due to juvenile terranes of Asia), which is not found in South America. In comparison with the internal versus external orogens recognized by Collins et al. (2011), as expected, our data from South America clearly resemble other external orogens (Fig. 12).

On a more recent time scale, our data record a dramatic change in crustal evolution during the past 200 m.y. that is not seen in older South American orogenic events or in the global record (Figs. 12–14). Magmatism in early Mesozoic magmatic arcs initiated with relatively evolved Hf signatures that reflect significant incorporation of older continental crust. Beginning at ca. 180 Ma, εHft values become more juvenile until ca. 100 Ma, and then become slightly more negative until ca. 35 Ma. The pre-100 Ma trend toward more juvenile values may result largely from extension within the magmatic arc, as has been proposed for the Tasmanides (Kemp et al., 2009; Hawkesworth et al., 2010). This excursion to more juvenile values may also have resulted from the Late Jurassic–Early Cretaceous magmatic arc having been located farther to the west, where there was little influence from older continental or continental margin assemblages.

As has been suggested in most previous syntheses (e.g., Haschke et al., 2002, 2006; Kay et al., 2005; DeCelles et al., 2009, 2015; Ramos, 2009; Ramos et al., 2014), the general trend toward more evolved values during the past ∼100 m.y. most likely resulted from the onset of shortening and crustal thickening within the Andes, causing more crustal melting, and also migration of magmatism eastward into older crust. The broad range of εHft values for <35 Ma zircons, from +15 to –15 (Fig. 11), results from profound variations in tectonic setting along the length of the Andes, with magmatism covering the full range of pure crustal recycling to mainly juvenile input from the mantle. Such a broad range of petrogenetic processes is rare during any time period in Earth history (Fig. 12).

Of particular interest for this study are the La/Yb, Sr isotope, and Nd isotope patterns reported by Haschke et al. (2002, 2006), which were interpreted to record changes in Andean subduction dynamics on a ∼50 m.y. time scale. As shown on Figure 13, the patterns appear to record increasing incorporation of crustal material toward the end of each cycle, as well as an overall trend toward more evolved values during the past 200 m.y. The data generated and compiled in this study from the same region (CA) show the general trend toward more evolved values, especially since ca. 100 Ma. However, the separate ∼50 m.y. pull-downs apparent from the Nd, Sr, and La/Yb data are not apparent in our Hf-in-zircon data set, and, as noted above and shown on Figure 13, our age data do not record the magmatic gaps interpreted by Haschke et al. (2002) to separate the magmatic cycles. A possible explanation of these discrepancies is that our record is dominated by felsic to intermediate composition magmatism, whereas the samples reported by Haschke et al. (2002) also include mafic components.

CONCLUSIONS

This study of U-Pb geochronology and Hf isotope geochemistry yields five major contributions to understanding the magmatic history and crustal evolution of western South America:

  1. Our U-Pb ages (n = 5524), combined with previously published ages (n = 36,764), record phases of high magmatic addition at 1.6–0.9 Ga, 700–400 Ma, 360–200 Ma, and 190 Ma to the present day (Fig. 12). Because this record has been assembled largely from studies of detrital zircons in rivers with headwaters in the modern Andes, the age distributions are biased toward felsic to intermediate composition magmatism, and toward younger events that are preserved in the Andean orogen.

  2. Magmatism during the past 200 m.y. has been variable north-to-south along the length of the Andes and east-west across the Andes, and igneous rocks versus detrital minerals yield slightly different patterns (Fig. 6). Most records are consistent, however, with age maxima at 183 (±11), 166 (±10), 149 (±9), 125 (±18), 110 (±9), 88 (±7), 65 (±6), 35 (±5), 21 (±6), and 4 (±4) Ma (2σ) (shaded gray bars on Fig. 6). Magmatic records for NA, CA, and SA have periodicities that range from ∼20 to ∼34 m.y. (Fig. 7). Although convergence rates between the South America and Pacific and Nazca oceanic plates have a somewhat similar periodicity of ∼20 m.y. (Maloney et al., 2013; Fig. 14), continuous wavelet transform analysis indicates that there is no temporal connection between magmatism and convergence rate (Fig. 15).

  3. The distribution of all ages from western South America has significant overlap with global age distributions, which are widely interpreted to record processes of crustal generation and/or preservation related to the supercontinent cycle (Hawkesworth et al., 2009, 2010; Belousova et al., 2010; Condie et al., 2011; Voice et al., 2011; Cawood et al., 2012; Condie, 2014). The main differences are the lesser abundance of >1.2 Ga ages and greater abundance of <200 Ma ages (Fig. 12) in our western South American record.

  4. Our Hf isotope determinations (n = 1199), combined with published Hf isotope data (n = 813), record eight phases in the crustal evolution of western South America (Fig. 11):

    • Prior to ca. 1.6 Ga, most grains record Paleoproterozoic reworking of older Archean crustal materials.

    • Between ca. 1.6 Ga and ca. 1.0 Ga, an increasing proportion of juvenile magma was mixed with reworked Paleoproterozoic–Archean crust.

    • Between ca. 1.0 Ga and ca. 500 Ma, most magmatic processes involved reworking of older crustal materials, with little generation of new crust.

    • Between ca. 500 and 400 Ma, increasing juvenile magmatism in CA resulted from extension and crustal thinning throughout the Famatina and related magmatic arcs (e.g., Ducea et al., 2015).

    • Following a magmatic lull between ca. 400 and ca. 360 Ma, there was significant juvenile magmatism in proto-Andean arc terranes, followed by a trend toward greater interaction with older crust until ca. 250 Ma. This increase in crustal reworking is presumably related to accretion and/or consolidation of oceanic arc terranes along the Andean margin.

    • Following a magmatic lull between ca. 240 and ca. 180 Ma, and continuing until ca. 100 Ma, magmatism increased, and significant volumes of juvenile magma were generated. This phase coincides with extensional tectonism (Ramos, 2009), a reduction in both convergence velocity and westward motion of South America (Maloney et al., 2013), and migration of the magmatic arc westward into more juvenile arc-type crust (e.g., Franz et al., 2006).

    • Between ca. 100 Ma and ca. 35 Ma, abundant magmatism continued, but with decreasing juvenile magma and increasing involvement of older continental crust. This phase presumably relates to contractional tectonism and initiation of uplift of the modern Andes (Kay et al., 2005; Franz et al., 2006; Ramos, 2009; DeCelles et al., 2015), an increase in both convergence velocity and westward motion of South America (Maloney et al., 2013), and migration of the magmatic arc eastward into more evolved crust (Haschke et al., 2002, 2006; Franz et al., 2006).

    • During the past ∼35 m.y., high-addition-rate magmatism continued with highly variable Hf isotope signatures. No evolutionary trends are apparent; although there is a spatial variation with more juvenile values in the northern and southern Andes and more evolved values in the central Andes (Figs. 11 and 14). These variations presumably reflect differing tectonic settings (contractional versus extensional) along and across the orogen in response to changes in both the subducting and overriding plates.

  5. The South American record of crustal evolution is quite similar to the global record (e.g., Belousova et al., 2010; Cawood et al., 2012; Condie, 2014), except that western South America lacks significant volumes of 800–650 Ma and 400–200 Ma juvenile crust and highly evolved 650–500 Ma crust (Fig. 12). The overall trend toward more juvenile magmatism is similar, however, to other external orogens of the world (Collins et al., 2011). The main pattern of increasing juvenile magmatism from 200 to 100 Ma and decreasing juvenile magmatism after 100 Ma (Fig. 12) may result from Jura-Cretaceous extension within and behind the Andean magmatic arc system, as well as migration of the magmatic arc first outboard into more juvenile arc terranes and then inboard into crust with greater proportions of older continental crust.

We are grateful for support from the Geological Society of America (student grant 9157-06) and the Chevron-Texaco summer support grant awarded to M. Pepper. Analyses were conducted in the Arizona LaserChron Center, which is supported by National Science Foundation grant EAR-1338583. We also thank LaserChron staff members Chen Li, Gayland Simpson, Chelsi White, Dominique Geisler, Mark Pecha, Clayton Loehn, and Glynis Jehle for assistance with sample preparation and analysis. Ned Brown provided valuable discussions of the work presented in this publication, and David Meko and Susan Beck provided assistance with our statistical analyses. We would also like to thank the ExxonMobil Convergent Orogenic Systems Analysis (COSA) project, which provided many discussions and covered costs of a central Andean field trip. We also thank Agustin Cardona, Franco Urbani, David Mendi, Camilo Bustamante, Giovani Jimenez, and Herman Bayona for assistance with sample collection during this project. A very thorough and insightful review was provided by H. Bahlburg.

1Supplemental Data 1. U-Pb data tables of northern Andes region. Please visit http://dx.doi.org/10.1130/GES01315.S1 or the full-text article on www.gsapubs.org to view the Supplemental Data.
2Supplemental Data 2. U-Pb data tables of central Andes region. Please visit http://dx.doi.org/10.1130/GES01315.S2 or the full-text article on www.gsapubs.org to view the Supplemental Data.
3Supplemental Data 3. U-Pb data tables of southern Andes region. Please visit http://dx.doi.org/10.1130/GES01315.S3 or the full-text article on www.gsapubs.org to view the Supplemental Data.
4Supplemental Item 1. Zircon properties. Please visit http://dx.doi.org/10.1130/GES01315.S4 or the full-text article on www.gsapubs.org to view the Supplemental Item.
5Supplemental Data 4. Hf data tables of northern Andes region. Please visit http://dx.doi.org/10.1130/GES01315.S5 or the full-text article on www.gsapubs.org to view the Supplemental Data.
6Supplemental Data 5. Hf data tables of central Andes region. Please visit http://dx.doi.org/10.1130/GES01315.S6 or the full-text article on www.gsapubs.org to view the Supplemental Data.
7Supplemental Data 6. Hf data tables of southern Andes region. Please visit http://dx.doi.org/10.1130/GES01315.S7 or the full-text article on www.gsapubs.org to view the Supplemental Data.
8Supplemental Data 7. Compilation of published U-Pb data. Please visit http://dx.doi.org/10.1130/GES01315.S8 or the full-text article on www.gsapubs.org to view the Supplemental Data.
9Supplemental Item 2. U-Pb compilation references. Please visit http://dx.doi.org/10.1130/GES01315.S9 or the full-text article on www.gsapubs.org to view the Supplemental Item.
10Supplemental Data 8. Compilation of published Hf data. Please visit http://dx.doi.org/10.1130/GES01315.S10 or the full-text article on www.gsapubs.org to view the Supplemental Data.
11Supplemental Item 3. Hf compilation references. Please visit http://dx.doi.org/10.1130/GES01315.S11 or the full-text article on www.gsapubs.org to view the Supplemental Item.
12Supplemental Item 4. Continuous wavelet transform analysis. Please visit http://dx.doi.org/10.1130/GES01315.S12 or the full-text article on www.gsapubs.org to view the Supplemental Item.