Abstract

Recognizing detrital contributions from sediment source regions is fundamental to provenance studies of active and ancient orogenic settings. Detrital zircon U-Pb geochronology of unconsolidated sands from modern rivers that have source catchments with contrasting bedrock signatures provides insight into the fidelity of U-Pb age signatures in discriminating tectonic provenance and downstream propagation of environmental signals. We present 1705 new detrital zircon U-Pb ages for 15 samples of unconsolidated river sands from 12 modern rivers over a large spatial extent of Ecuador (∼1°N–5°S and ∼79°–77°W).

Results show distinctive U-Pb age distributions with characteristic zircon age populations for various tectonic provinces along the Andean convergent margin, including the forearc, magmatic arc, and internal (hinterland) and external (foreland) segments of the fold-thrust belt. (1) Forearc and magmatic arc (Western Cordillera) river sands are characterized by Neogene–Quaternary age populations from magmatic sources. (2) Rivers in the hinterland (Eastern Cordillera) segment of the Andean fold-thrust belt have substantial populations of Proterozoic and Paleozoic ages, representing upper Paleozoic–Mesozoic sedimentary and metasedimentary rocks of ultimate cratonic origin. (3) River sands in the frontal fold-thrust belt (Subandean Zone to Oriente Basin) show distinctive bimodal Jurassic age populations, a secondary Triassic population, and subordinate Early Cretaceous ages representative of Mesozoic plutonic and metamorphic bedrock.

Detrital zircon U-Pb results from a single regional watershed (Rio Pastaza) spanning the magmatic arc to foreland basin show drastic downstream variations, including the downstream loss of magmatic arc and hinterland signatures and abrupt introduction and dominance of selected sources within the fold-thrust belt. Disproportionate contributions from Mesozoic crystalline metamorphic rocks, which form high-elevation, high-relief areas subject to focused precipitation and active tectonic deformation, are likely the product of focused erosion and high volumes of local sediment input from the frontal fold-thrust belt, leading to dilution of upstream signatures from the hinterland and magmatic arc.

INTRODUCTION

Detrital zircon U-Pb geochronology has become a fundamental provenance tool providing critical records of mountain building, large-scale drainage evolution related to tectonic and climatic drivers, and source-to-sink evolution of depositional systems and sedimentary basins (e.g., Dickinson and Gehrels, 2003, 2008; Weislogel et al., 2006; Cawood et al., 2012; Laskowski et al., 2013; Blum and Pecha, 2014; Sharman et al., 2015; Horton, 2018; Ducea et al., 2018; Wang et al., 2018). Discrimination of competing sedimentary source regions and terranes is essential to detrital zircon provenance studies (e.g., DeCelles et al., 2000; Horton et al., 2010a, 2010b, 2015; Gehrels, 2014; Saylor and Sundell, 2016; Sickmann et al., 2018). However, despite the utility of detrital zircon U-Pb data sets, there are important limitations involved in interpretation of such data. At a fundamental level, do detrital zircon data accurately reflect the distribution of source units exposed in their river catchments? The commonly implicit assumption that detrital contributions are proportional to source unit exposure area within a catchment (Saylor et al., 2013) can be compromised in situations where there are disproportionate contributions from different bedrock units, due to: (1) heterogeneous zircon fertility (Moecher and Samson, 2006; Dickinson, 2008; Spencer et al., 2018); (2) uneven bedrock erodibility (Capaldi et al., 2017); (3) contrasting lithologic breakdown by weathering, erosion, and climate (Nesbitt et al., 1996; Amidon et al., 2005); and (4) varied mechanical durability during transport (e.g., Garzanti et al., 2013).

In this study, we present a large data set of 1705 new zircon U-Pb ages from unconsolidated sand samples collected from modern rivers across Ecuador in order to characterize detrital zircon U-Pb signatures of the different tectonic provinces and explore downstream variations in detrital signatures within a single large river system. Sampled rivers span the coastal forearc to retroarc foreland basin of Ecuador (∼80°–77°W) and have Andean mountainous catchment areas ranging from small (∼200 km2) to large (∼13,500 km2), covering a cumulative surface area of ∼50,000 km2 (Figs. 1 and 2). First, the results show that U-Pb age distributions of modern rivers in Ecuador clearly discriminate the major tectonic provinces (forearc, magmatic arc, and hinterland and foreland segments of the fold-thrust belt). Second, U-Pb age distributions along the large Rio Pastaza catchment (∼13,500 km2) show a downstream loss of important age populations, likely due to pronounced input from singular sources due to focused erosion driven by high precipitation and/or active tectonic processes.

Figure 1.

Regional map of Ecuador showing simplified surface geology and tectonic provinces (Baldock, 1982; Litherland et al., 1993; Pratt et al., 2005; Cochrane, 2013; Spikings et al., 2015). Tectonic province boundaries are indicated with dashed black lines. Locations of samples 1–15 (this study) and samples I–IV (Pepper et al., 2016) are represented by filled circles and correspond to outlined catchment (drainage basin) areas, which were created from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) digital elevation model V002 (https://asterweb.jpl.nasa.gov/gdem.asp) with 30 m resolution. Inset map in the lower right corner shows cratonic provinces discussed in the text (after Chew et al., 2007).

Figure 1.

Regional map of Ecuador showing simplified surface geology and tectonic provinces (Baldock, 1982; Litherland et al., 1993; Pratt et al., 2005; Cochrane, 2013; Spikings et al., 2015). Tectonic province boundaries are indicated with dashed black lines. Locations of samples 1–15 (this study) and samples I–IV (Pepper et al., 2016) are represented by filled circles and correspond to outlined catchment (drainage basin) areas, which were created from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) digital elevation model V002 (https://asterweb.jpl.nasa.gov/gdem.asp) with 30 m resolution. Inset map in the lower right corner shows cratonic provinces discussed in the text (after Chew et al., 2007).

Figure 2.

Regional map of Ecuador showing digital elevation model (DEM) topography. Tectonic province boundaries are indicated with dashed black lines. Solid white line delineates the continental divide. Locations of samples 1–15 (this study) and samples I–IV (Pepper et al., 2016) are represented by filled circles and correspond to outlined catchment (drainage basin) areas, which were created from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) DEM V002 (https://asterweb.jpl.nasa.gov/gdem.asp) with 30 m resolution.

Figure 2.

Regional map of Ecuador showing digital elevation model (DEM) topography. Tectonic province boundaries are indicated with dashed black lines. Solid white line delineates the continental divide. Locations of samples 1–15 (this study) and samples I–IV (Pepper et al., 2016) are represented by filled circles and correspond to outlined catchment (drainage basin) areas, which were created from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) DEM V002 (https://asterweb.jpl.nasa.gov/gdem.asp) with 30 m resolution.

GEOLOGIC SETTING

Ecuador is divided into various tectonic provinces with distinctive rock units of different age and tectonic history. These provinces are parallel to the regional tectonic strike of the northern Andes, with roughly north-south orientations that broadly coincide with major topographic boundaries (Figs. 1 and 2). Here we classify these tectonic provinces, according to their relative position, as different components of a Cordilleran-style orogen along an ocean-continent convergent plate boundary. We recognize that the precise boundaries between successive forearc, magmatic arc, retroarc fold-thrust belt, and foreland basin provinces are likely subject to internal complexities wherein these boundaries may deviate such a simplified framework at finer spatial scales. From west to east, the major tectonic provinces include the coastal forearc region, the magmatic arc (Western Cordillera), the internal hinterland segment of the fold-thrust belt (Eastern Cordillera), the external frontal segment of the fold-thrust belt (Subandean Zone), and the Oriente foreland basin.

Forearc and Magmatic Arc (Western Cordillera)

The coastal forearc and magmatic arc (Western Cordillera) region (Figs. 1 and 2) is composed of allochthonous mafic and ultramafic Upper Cretaceous rocks of oceanic plateau and/or intraoceanic island arc origin that accreted to the margin ca. 75 Ma, overlain by Upper Cretaceous–Paleogene sedimentary rocks and Paleogene–Quaternary magmatic rocks of both volcanic and plutonic origin (Kerr et al., 2002; Jaillard et al., 2004; Vallejo et al., 2009; Schütte et al., 2010; Spikings et al., 2015). Paleocene–Miocene (66–5.3 Ma) magmatic rocks include dominantly granodioritic intrusions and andesitic to rhyolitic lava flows, ignimbrites, and tuffs (Schütte et al., 2010, and references therein). Pliocene–Quaternary (5.3–0 Ma) magmatism is manifest as a volcanic chain that persists from Ecuador (∼2.5°S) northward into Colombia. In Ecuador alone, >80 recognized volcanoes (mostly stratovolcanoes) with Pliocene–Quaternary activity have produced massive amounts of volcanic material ranging from andesitic to rhyolitic lava flows, ignimbrites, and tuffs (Hall and Wood, 1985; Barberi et al., 1988; Hall et al., 2008; Bernard and Andrade, 2011; Jackson et al., 2019).

Hinterland Fold-Thrust Belt (Eastern Cordillera)

The Eastern Cordillera represents the internal hinterland segment of the fold-thrust belt (Figs. 1 and 2) composed of sublinear north-trending belts of variably metamorphosed and variably deformed upper Paleozoic–Lower Cretaceous sedimentary, igneous, and volcanic rocks. In general, these include: (1) upper Paleozoic–Lower Cretaceous low- to intermediate-grade metasedimentary rocks formed from pelitic protoliths ultimately derived from cratonic sources; (2) Triassic–Jurassic low- to high-grade metaigneous rocks and Triassic S-type granitoids; and (3) Jurassic–Lower Cretaceous volcanic and low-grade metavolcanic rocks (Aspden and Litherland, 1992; Aspden et al., 1992; Pratt et al., 2005; Ruiz et al., 2007; Spikings et al., 2015). Mesozoic metamorphism is linked to a complex, spatially variable history of sedimentary burial, crustal shortening, and terrane collision (Litherland et al., 1994; Kerr et al., 2002). Finally, several isolated Quaternary volcanic centers have developed on top of metamorphic basement in the Eastern Cordillera (Hall et al., 2008), but represent a limited cumulative area (<5%–10%) of this tectonic province.

Frontal Fold-Thrust Belt (Subandean Zone to Oriente Basin)

The Subandean Zone to Oriente Basin region represents the outer, frontal segment of the fold-thrust belt and associated flexural foreland basin (Figs. 1 and 2). The Oriente Basin is part of the larger hydrocarbon-bearing Subandean retroarc foreland basin system and is continuous southward with the Peruvian Marañón Basin and northward with the Colombian Putumayo Basin (Dashwood and Abbotts, 1990; Higley, 2001). The Cretaceous–Cenozoic sedimentary succession of the Oriente Basin is flat-lying to the east and locally deformed in the west by active fold-thrust structures of the frontal Subandean wedge-top depozone (Bès de Berc et al., 2005; Baby et al., 2013).

Mesozoic basin-fill strata have cratonic derivation and include sandstones, siltstones, organic-rich shales, and subordinate limestones (Tschopp, 1953; Dashwood and Abbotts, 1990; Barragán et al., 2004; Ruiz et al., 2007; Horton et al., 2010a; Vallejo et al., 2017). Cenozoic strata include coarser clastic material (sandstones, siltstones, and conglomerates) derived from Andean sources to the west (Christophoul et al., 2002; Jaillard et al., 1997; Martin-Gombojav and Winkler, 2008; Horton, 2018; Gutiérrez et al., 2019). The Cretaceous–Cenozoic succession unconformably overlies Jurassic magmatic arc rocks, which are exposed exclusively in the Subandean Zone as non-metamorphosed granites, monzogranites, and granodiorites (Litherland et al., 1994; Romeuf et al., 1995; Spikings et al., 2001; Pratt et al., 2005; Spikings et al., 2015).

METHODS

Detrital Zircon Geochronology

Uranium-lead (U-Pb) geochronological analyses of detrital zircon grains were performed on 15 unconsolidated sand samples collected from 12 rivers in Ecuador (Figs. 1 and 2; Supplemental Table S11). To reduce potential intersample bias that may result from hydrodynamic sorting and grain-size dependence (e.g., Malusà et al., 2016; Ibañez-Mejia et al., 2018), modern fluvial bar deposits of medium sand size were routinely targeted for sampling. Individual zircon crystals were density separated using a Wilfley water table and Frantz magnetic separator, followed by heavy liquid treatment using methylene iodide (density 3.33 g/cm3). Samples containing residual pyrite were washed in nitric acid. The final heavy fractions were then poured onto tape (rather than picked by hand), mounted in epoxy resin, polished, and imaged by backscatter electron imaging. Individual zircon grains were then chosen at random for analysis. U-Pb analyses were performed at the Arizona LaserChron Center (Tucson, Arizona, USA) by laser ablation inductively coupled plasma mass spectrometry on an Element2 HR ICPMS (Gehrels et al., 2006, 2008; Gehrels, 2014). Zircon standards used were Sri Lanka, FC-1, and R33. Analyses exhibiting >20% discordance, >5% reverse discordance, or >10% internal uncertainty were removed from further consideration. The reported ages represent 206Pb/238U ages for zircons <900 Ma and 206Pb/207Pb ages for zircons >900 Ma. For very young zircons (<2 Ma), a U-Th disequilibrium correction was applied to 206Pb/238U ages following the approach of Crowley et al. (2007), assuming a magma Th/U value of 2.3. Uncertainties are reported at the 1σ level. Probability distribution plots for data visualization were created using detritalPy software (Sharman et al., 2018). Grain dimensions were acquired for a subset of analyzed zircons, with measurements of long and short axes, and confirm that diagnostic age populations are not restricted to particular grain sizes.

Major Age Populations

Major age populations in the northern Andes have been identified from geochronological studies of crystalline rock units and detrital zircon U-Pb results from sedimentary units defining the Andean belt and adjacent sectors of northern South America (e.g., Teixeira et al., 1989; Litherland et al., 1994; Cordani and Teixeira, 2007; Cardona et al., 2010; Horton et al., 2010b, 2015; Saylor et al., 2013). Considerations of both Andean rock units and cratonic signatures are important due to the prevalence of upper Paleozoic–Mesozoic sedimentary rocks units composed of clastic material derived from the Amazonian craton (e.g., Dashwood and Abbotts, 1990; Cooper et al., 1995; Horton et al., 2010b; Vallejo et al., 2017; Gutiérrez et al., 2019).

In this study, major age populations (Fig. 3) are defined on the basis of significant age groups recognized in Cretaceous–Cenozoic basin fill in forearc and retroarc basins (Witt et al., 2017; Gutiérrez et al., 2019). Undisputed exposures of Precambrian basement rocks are restricted to cratonic provinces farther east and do not exist in the study area; thus Precambrian ages present in studied samples must have been originally derived from the craton and were subsequently incorporated into younger rock units. Basement age groups may be divided into Andean (<650 Ma) and cratonic (>900 Ma) affinities, recognizing that older Precambrian ages are common in upper Paleozoic–Mesozoic sedimentary and metasedimentary rocks of the northern Andes. Andean populations are dominated by Mesozoic–Cenozoic (250–0 Ma) ages common in magmatic units with subordinate Paleozoic (450–250 Ma), Famatinian orogeny (550–450 Ma), and late Neoproterozoic (650–550 Ma) ages of crystalline basement origin (Chew et al., 2007, 2008; Cardona et al., 2010; Gutiérrez et al., 2019). Ages >900 Ma are of cratonic affinity, with an original derivation from Precambrian cratonic provinces to the east, including: >2300 Ma (Central Amazonian); 2250–2150 Ma (Maroni-Itacaiúnas); 2150–1800 Ma (Ventuari-Tapajós); 1800–1550 Ma (Rio Negro–Jurena); 1550–1300 Ma (Rondonia–San Ignacio); and 1300–900 Ma (Sunsás) ages (Fig. 1, lower right inset) (Teixeira et al., 1989; Cordani et al., 2005; Cordani and Teixeira, 2007; Chew et al., 2007, Santos et al., 2008; Bahlburg et al., 2009; Bettencourt et al., 2010; Cardona et al., 2010).

Figure 3.

Composite plots of detrital zircon U-Pb age distributions for all detrital zircon U-Pb results from modern river sand samples in Ecuador. Plots show 0–300 Ma distribution (upper) and 0–2150 Ma distribution (lower), with each distribution defined by an age probability density function (thick black line), kernel density estimate (thin red line), and age histogram (open rectangles). Color bars highlight distinctive age populations in the composite plots, and match the pie diagram (lower right) showing color-coded age proportions for diagnostic age categories recognized in Ecuador (after Gutiérrez et al., 2019). Pl-Q—Pliocene–Quaternary; FMT—Famatinian orogeny; NPR—Neoproterozoic; SA—Sunsás; RSI—Rondonia–San Ignacio; RNJ—Rio Negro–Jurena; VT—Ventuari-Tapajós.

Figure 3.

Composite plots of detrital zircon U-Pb age distributions for all detrital zircon U-Pb results from modern river sand samples in Ecuador. Plots show 0–300 Ma distribution (upper) and 0–2150 Ma distribution (lower), with each distribution defined by an age probability density function (thick black line), kernel density estimate (thin red line), and age histogram (open rectangles). Color bars highlight distinctive age populations in the composite plots, and match the pie diagram (lower right) showing color-coded age proportions for diagnostic age categories recognized in Ecuador (after Gutiérrez et al., 2019). Pl-Q—Pliocene–Quaternary; FMT—Famatinian orogeny; NPR—Neoproterozoic; SA—Sunsás; RSI—Rondonia–San Ignacio; RNJ—Rio Negro–Jurena; VT—Ventuari-Tapajós.

RESULTS

Detrital zircon U-Pb geochronological results for 15 river sand samples from 12 modern Ecuadorian rivers at 1°N–5°S are presented with a series of diagrams defining the age distributions (probability density functions) in the context of the aforementioned diagnostic age populations of the northern Andes. A composite age distribution (Fig. 3) is defined by 2057 zircon U-Pb analyses, the sum of results from our 15 river samples (1705 new analyses) and 352 published analyses of four additional river samples (Pepper et al., 2016). The composite results reveal diverse age populations spanning the Precambrian through the Cenozoic, with Andean ages constituting 79% of all analyses (65% Cenozoic–Mesozoic and 14% Andean basement), and cratonic ages (>900 Ma) constituting the remaining 21% (Fig. 3; Table S1 [footnote 1]).

Among Cenozoic ages, 51% are Pliocene–Quaternary, 12% are Miocene, and 36% are Paleogene (dominantly Eocene–Oligocene). Among Mesozoic ages, over half (54%) are Jurassic, and the remainder are equal contributions of Cretaceous and Triassic (each 23%). Older Precambrian ages (>900 Ma) decrease in proportion with increasing age, consistent with decreasing proximity to the oldest cratonic core in the distal easternmost regions (Fig. 1). Among all cratonic ages in the cumulative data set (n = 423 analyses), 50% are 1300–900 Ma (Sunsás), 18% are 1550–1300 Ma (Rondonia–San Ignacio), 13% are 1800–1550 Ma (Rio Negro–Jurena), 13% are 2150–1800 Ma (Ventuari-Tapajós), and only 6% are >2150 Ma. The age distributions are presented for the different tectonic provinces according to sample location, river name, and catchment size (Table S2 [footnote 1]).

Forearc and Magmatic Arc (Western Cordillera)

Six new samples of rivers with individual catchment areas ranging from 197 to 3555 km2 in the forearc and magmatic arc and three published samples represent a cumulative catchment area of 18,222 km2 (Figs. 1 and 2). New geochronological results (n = 653 analyses) together with results from three published samples (n = 283 analyses; Pepper et al., 2016) show an overall dominance of Cenozoic, primarily Neogene (23–0 Ma) and late Paleogene (dominantly ca. 40–23 Ma) ages, with limited pre-Cenozoic ages (Fig. 4). Jurassic ages are virtually absent (<1% of analyses), and Paleozoic and older ages are also nearly absent in seven of the nine samples (Fig. 4).

Figure 4.

Plots of detrital zircon U-Pb age distributions for modern river sand samples from the forearc and Western Cordillera of Ecuador (locations shown in Figs. 1 and 2), including new samples 1–6 (this study) and published samples I–III (Pepper et al., 2016). Plots show 0–300 Ma distribution (left) and 0–2150 Ma distribution (right), with each distribution defined by an age probability density function (thick black line), kernel density estimate (thin red line), and age histogram (open rectangles). Color bars highlight distinctive age populations in the composite plots, and match the pie diagram (far right) showing color-coded age proportions for diagnostic age categories (see Fig. 3, including for abbreviation definitions).

Figure 4.

Plots of detrital zircon U-Pb age distributions for modern river sand samples from the forearc and Western Cordillera of Ecuador (locations shown in Figs. 1 and 2), including new samples 1–6 (this study) and published samples I–III (Pepper et al., 2016). Plots show 0–300 Ma distribution (left) and 0–2150 Ma distribution (right), with each distribution defined by an age probability density function (thick black line), kernel density estimate (thin red line), and age histogram (open rectangles). Color bars highlight distinctive age populations in the composite plots, and match the pie diagram (far right) showing color-coded age proportions for diagnostic age categories (see Fig. 3, including for abbreviation definitions).

Hinterland Fold-Thrust Belt (Eastern Cordillera)

Six new samples from rivers with individual catchment areas ranging from 341 to 8371 km2 from the Eastern Cordillera hinterland segment of the fold-thrust belt represent a cumulative catchment area of 21,671 km2 (Figs. 1 and 2). Geochronological results (n = 578 analyses) exhibit overall cosmopolitan age distributions dominated by Paleozoic–Precambrian ages (constituting 61% of analyses) (Fig. 5). Cenozoic ages are generally restricted to the Paleogene (dominantly ca. 50–30 Ma), as Neogene ages (23–0 Ma) are nearly absent, constituting only 2% of analyses. Mesozoic ages are principally Triassic (200–250 Ma), with minor Jurassic–Cretaceous (200–66 Ma) ages. Paleozoic (450–250 Ma), Famatinian (550–450 Ma), and Neoproterozoic (650–550 Ma) ages are present in five of six samples in the Eastern Cordillera province (all but sample 11: EC1535). Older Precambrian age populations are present in all samples of this region, with the most prevalent, the Sunsás group (1300–900 Ma), appearing in all samples, followed by Rondonia–San Ignació (1550–1300 Ma), Rio Negro–Jurena (1800–1550 Ma), and Ventuari-Tapajós (2150–1800 Ma).

Figure 5.

Plots of detrital zircon U-Pb age distributions for modern river sand samples from the Eastern Cordillera of Ecuador (locations shown in Figs. 1 and 2), including new samples 7–12 (this study). Plots show 0–300 Ma distribution (left) and 0–2150 Ma distribution (right), with each distribution defined by an age probability density function (thick black line), kernel density estimate (thin red line), and age histogram (open rectangles). Color bars highlight distinctive age populations in the composite plots, and match the pie diagram (far right) showing color-coded age proportions for diagnostic age categories (see Fig. 3, including for abbreviation definitions).

Figure 5.

Plots of detrital zircon U-Pb age distributions for modern river sand samples from the Eastern Cordillera of Ecuador (locations shown in Figs. 1 and 2), including new samples 7–12 (this study). Plots show 0–300 Ma distribution (left) and 0–2150 Ma distribution (right), with each distribution defined by an age probability density function (thick black line), kernel density estimate (thin red line), and age histogram (open rectangles). Color bars highlight distinctive age populations in the composite plots, and match the pie diagram (far right) showing color-coded age proportions for diagnostic age categories (see Fig. 3, including for abbreviation definitions).

Frontal Fold-Thrust Belt (Subandean Zone to Oriente Basin)

Three sand samples from rivers with individual catchment areas ranging from 4069 to 13,473 km2 in the Subandean Zone represent a cumulative catchment area of 19,351 km2 (Figs. 1 and 2). The new geochronological results (n = 474 analyses) together with those from one published sample (n = 69 analyses; Pepper et al., 2016) show age populations generally restricted to the Mesozoic, and Jurassic ages in particular, which constitute 59% of all analyses for this region and exhibit a clear bimodal distribution with age peaks corresponding to the Early Jurassic (29% of analyses; 200–174 Ma) and Late Jurassic (29% of analyses; 165–145 Ma) (Fig. 6). Secondary populations are Triassic (250–200 Ma) followed by limited Early Cretaceous (145–110 Ma) ages. Late Cretaceous and Cenozoic ages are virtually absent, as are Paleozoic and older populations, except in one sample (sample IV: RTENEC) where limited pre-Mesozoic ages are present.

Figure 6.

Plots of detrital zircon U-Pb age distributions for modern river sand samples from the Subandean Zone of Ecuador (locations shown in Figs. 1 and 2), including new samples 13–15 (this study) and published sample IV (Pepper et al., 2016). Plots show 0–300 Ma distribution (left) and 0–2150 Ma distribution (right), with each distribution defined by an age probability density function (thick black line), kernel density estimate (thin red line), and age histogram (open rectangles). Color bars highlight distinctive age populations in the composite plots, and match the pie diagram (far right) showing color-coded age proportions for diagnostic age categories (see Fig. 3, including for abbreviation definitions).

Figure 6.

Plots of detrital zircon U-Pb age distributions for modern river sand samples from the Subandean Zone of Ecuador (locations shown in Figs. 1 and 2), including new samples 13–15 (this study) and published sample IV (Pepper et al., 2016). Plots show 0–300 Ma distribution (left) and 0–2150 Ma distribution (right), with each distribution defined by an age probability density function (thick black line), kernel density estimate (thin red line), and age histogram (open rectangles). Color bars highlight distinctive age populations in the composite plots, and match the pie diagram (far right) showing color-coded age proportions for diagnostic age categories (see Fig. 3, including for abbreviation definitions).

DISCRIMINATION OF TECTONIC PROVINCES

Comparison of cumulative U-Pb age distributions by tectonic province reveals several diagnostic populations (Fig. 7). (1) The coastal forearc to magmatic arc (Western Cordillera) is characterized by Pliocene–Quaternary (5.3–0 Ma) and Paleogene (dominantly ca. 40–23 Ma) age populations. (2) The hinterland fold-thrust belt (Eastern Cordillera) is characterized by Paleogene (dominantly ca. 50–30 Ma), Triassic (250–200 Ma), minor Early Jurassic (200–145 Ma), Andean basement (650–250 Ma), and older Precambrian (>900 Ma) populations; Paleogene ages in the Eastern Cordillera hinterland are older than those in the forearc and magmatic arc. (3) The frontal fold-thrust belt (Subandean Zone to Oriente Basin) is characterized by two prominent Jurassic age peaks and a secondary Triassic age population that exhibits a broad distribution spanning most of the Triassic (200–250 Ma).

Figure 7.

Comparative plots of detrital zircon U-Pb age distributions (below) and composite cumulative distributions (above) for all 19 modern river sand samples from Ecuador, grouped according to tectonic province (N = 9 forearc–Western Cordillera, N = 6 Eastern Cordillera, and N = 4 Subandean Zone samples). Plots show 0–300 Ma distribution (left) and 0–2150 Ma distribution (right), with each distribution defined by an age probability density function (thick black line), kernel density estimate (thin red line), and age histogram (open rectangles). Color bars highlight distinctive age populations in the composite plots, and match the pie diagram (far right) showing color-coded age proportions for diagnostic age categories (see Fig. 3, including for abbreviation definitions).

Figure 7.

Comparative plots of detrital zircon U-Pb age distributions (below) and composite cumulative distributions (above) for all 19 modern river sand samples from Ecuador, grouped according to tectonic province (N = 9 forearc–Western Cordillera, N = 6 Eastern Cordillera, and N = 4 Subandean Zone samples). Plots show 0–300 Ma distribution (left) and 0–2150 Ma distribution (right), with each distribution defined by an age probability density function (thick black line), kernel density estimate (thin red line), and age histogram (open rectangles). Color bars highlight distinctive age populations in the composite plots, and match the pie diagram (far right) showing color-coded age proportions for diagnostic age categories (see Fig. 3, including for abbreviation definitions).

In this study, the sediment sources are explicitly known, in the sense that zircon grains in the modern rivers are derived from surface waters within individual stream catchments with known surface geology. Nevertheless, to independently verify diagnostic age peaks within the sediment source areas, we applied the non-negative matrix factorization (NMF) approach of Saylor et al. (2019). This inverse technique generates a series of age distributions representative of potential sources. The NMF model results for the 19 modern river samples (Fig. S1 [footnote 1]) demonstrate the diagnostic Mesozoic–Cenozoic age populations within the Subandean Zone (two Jurassic age peaks), Eastern Cordillera (Eocene, Early Cretaceous, and broad Triassic age peaks), and Western Cordillera (Pliocene–Quaternary, Miocene, Oligocene, and Late Cretaceous age peaks).

A statistical test for similarity among the 19 modern river samples using multidimensional scaling (MDS; Vermeesch, 2012) shows clear clustering of sample U-Pb age distributions into zones representing the three contrasting tectonic provinces—the forearc to magmatic arc (Western Cordillera), the hinterland segment of the fold-thrust belt (Eastern Cordillera), and the foreland segment of the fold-thrust belt (Subandean Zone to Oriente Basin) (Fig. 8). Only two samples have age distributions that appear outside of their respective fields. The first is a sample (sample 11: EC1535) from the Rio Paute in the Eastern Cordillera that appears to share both Eastern and Western Cordillera signatures in the MDS plot. This case is readily explained by the spatial overlap of the catchment into both the Eastern and Western Cordillera regions (Figs. 1 and 2), each accounting for roughly half of the total 5682 km2 catchment area for this sample.

Figure 8.

Multidimensional scaling plot showing quantitative similarity among detrital zircon U-Pb age distributions for all 19 modern river sand samples, where spatial proximity indicates statistical similarity between age distributions for different samples. Samples are color coded by region, and color ellipses highlight tectonic provinces distinguishable by their U-Pb age distributions (see text for discussion).

Figure 8.

Multidimensional scaling plot showing quantitative similarity among detrital zircon U-Pb age distributions for all 19 modern river sand samples, where spatial proximity indicates statistical similarity between age distributions for different samples. Samples are color coded by region, and color ellipses highlight tectonic provinces distinguishable by their U-Pb age distributions (see text for discussion).

The second apparent outlier in the MDS plot is from the Rio Pilaló in Western Cordillera (sample 3: EC1505), which appears to cluster very well with results diagnostic of the Eastern Cordillera (Fig. 8). The U-Pb age distribution for this sample (Fig. 4; n = 111 analyses) consists of a dominant Precambrian population (>900 Ma, 60%), primarily defined by Sunsás (1300–900 Ma) ages that constitute 49% of analyzed grains (54 of 111 analyses), with a subordinate range of ages >1300 Ma. The second-largest population is defined by Paleozoic (450–250 Ma, 13%) grains followed by 6% Famatinian (550–450 Ma), 5% late Neoproterozoic (650–550 Ma), and 6% early Neoproterozoic (900–650 Ma) ages. Although Neogene ages are hallmarks of the overall age distributions for the Western Cordillera region (Figs. 4 and 7), the Rio Pilaló sample is nearly devoid of Neogene zircons (Fig. 4, sample 3). This sample from the Rio Pilaló has the second-smallest catchment area of rivers in this study, representing only 251 km2, and a catchment composed of only Paleogene rocks (Table S2 [footnote 1]; Figs. 1 and 2). Detrital zircon geochronological results from Vallejo et al. (2019) for Paleogene stratigraphic units equivalent to those in this catchment show similar detrital U-Pb age distributions. These distributions record sedimentary input from source units in the Eastern Cordillera (Vallejo et al., 2019).

DOWNSTREAM TRENDS IN A LARGE CATCHMENT

River sand samples in the Subandean Zone represent large watersheds with western headwaters in the high-elevation Andean hinterland and magmatic arc (Fig. 2). Tributaries and main-stem rivers cross the Eastern Cordillera hinterland where U-Pb age distributions have characteristic Paleozoic and Precambrian age signatures (Figs. 5 and 6). The westernmost headwaters of the Rio Pastaza drain considerable portions of the magmatic arc in which the catchments are dominated by Pliocene–Quaternary volcanic and volcaniclastic rocks (Fig. 1; Barberi et al., 1988) (e.g., Rio Pastaza and Rio Napo watersheds, samples 15: EC1601 and 13: EC1602).

Given this scenario, it is surprising that U-Pb distributions of Subandean rivers lack the Pliocene–Quaternary and Precambrian–Paleozoic age populations diagnostic of the magmatic arc and hinterland regions, respectfully (Figs. 6 and 7). Lacking evidence for any significant sediment storage or depositional sink within upstream sectors of the catchments, we attribute the unexpected U-Pb distributions to dilution by disproportionately large detrital contributions from crystalline Mesozoic units. Five samples from the Rio Pastaza catchment demonstrate this abrupt downstream loss of important U-Pb age populations and overwhelming dilution by Mesozoic ages (Fig. 9). Measurements of the long and short axes of individual analyzed grains (Table S3 [footnote 1]) demonstrate that variations in age populations among these samples cannot be attributed to grain-size fractionation.

Figure 9.

Plots of detrital zircon U-Pb age distributions for modern river sand samples from the Rio Pastaza catchment area in Ecuador (locations shown in Figs. 1 and 2). Samples are arranged from upstream (top) to downstream (bottom) sample locations (Fig. 10). Plots show 0–300 Ma distribution (left) and 0–2150 Ma distribution (right), with each distribution defined by an age probability density function (thick black line), kernel density estimate (thin red line), and age histogram (open rectangles). Color bars highlight distinctive age populations in the composite plots, and match the pie diagram (far right) showing color-coded age proportions for diagnostic age categories (see Fig. 3, including for abbreviation definitions).

Figure 9.

Plots of detrital zircon U-Pb age distributions for modern river sand samples from the Rio Pastaza catchment area in Ecuador (locations shown in Figs. 1 and 2). Samples are arranged from upstream (top) to downstream (bottom) sample locations (Fig. 10). Plots show 0–300 Ma distribution (left) and 0–2150 Ma distribution (right), with each distribution defined by an age probability density function (thick black line), kernel density estimate (thin red line), and age histogram (open rectangles). Color bars highlight distinctive age populations in the composite plots, and match the pie diagram (far right) showing color-coded age proportions for diagnostic age categories (see Fig. 3, including for abbreviation definitions).

Upper Catchment: Magmatic Arc

The Rio Pastaza is formed by two longitudinal tributaries (the Rio Chambo and Rio Patate) within the magmatic arc. The Rio Chambo drains the northern half and the Rio Patate the southern half of the Pastaza’s upper catchment. Outcrops within the upper catchment are almost entirely Neogene–Quaternary volcanic and volcaniclastic rocks, and >20 volcanoes with documented Pliocene–Quaternary activity (Bernard and Andrade, 2011). A sample from the Rio Chambo (sample 4: EC1509) represents 2016 km2 of the Rio Pastaza’s upper catchment. The U-Pb signature of the sample is characterized by a unimodal Pliocene–Quaternary (5.3–0 Ma) age distribution (Fig. 9), consistent with sediment sources in the catchment.

Middle Catchment: Hinterland

Two river sand samples from the middle part of the Rio Pastaza’s catchment (Fig. 9, samples 9: EC1510 and 10: EC1513) within the hinterland fold-thrust belt province show a near absence of Pliocene–Quaternary zircons, and an abundance of Precambrian zircons >900 Ma, which constitute 47% and 42% of analyses for the two samples, respectively. Secondary age populations in the two mid-catchment samples are late Neoproterozoic–Famatinian (650–450 Ma; 12% and 10%, respectively) and Paleozoic–Triassic (450–200; 11% and 13%, respectively) (Fig. 9). While these U-Pb age distributions are consistent with those of other river sand samples from the Eastern Cordillera hinterland region (Fig. 5), the sample distributions bear no resemblance to that of the upper catchment sample, and are severely disproportionate to the relative exposure areas. This relationship is particularly striking because the large upstream catchment area is dominated by magmatic arc rocks, with >85% Cenozoic rocks (including >50% Pliocene–Quaternary volcanic rocks).

Lower Catchment: Proximal and Distal Subandean Zone

In the lower part of the Pastaza catchment, two samples represent both the proximal (western) part of the Subandean Zone and the distal (eastern) part of the Subandean Zone to Oriente Basin. A lower catchment sample from the proximal Subandean Zone (Fig. 9, sample 14: EC15191606) represents a 9252 km2 catchment area in which Cenozoic rocks make up 78% of exposed rock units (including 47% Quaternary volcanic rocks). Paleozoic–Triassic rocks constitute 3%, Jurassic rocks 4%, and Cretaceous rocks the remaining 15% of exposures in the catchment area. The U-Pb age distribution for the sample shows dominantly Mesozoic (250–66 Ma) ages (84%) that are principally Jurassic (145–200 Ma; 74%), scarce Precambrian (>900 Ma) ages (14%), and a virtual absence of Cenozoic (66–0 Ma) ages. Jurassic populations in these samples exhibit a bimodal distribution of Early (200–174 Ma; 36%) and Late Jurassic (165–145 Ma; 38%) ages.

The most downstream sample (sample 15: EC1601) is located near the transition from distal Subandean Zone to foreland basin and represents 13,473 km2, or 98%, of the total 13,700 km2 catchment of the Rio Pastaza when it debouches onto the Amazon Plain in the Pastaza megafan (Figs. 1 and 2; Bès de Berc et al., 2005; Bernal et al., 2011). The detrital zircon U-Pb age distribution of this sample is comparable with those of the previous two lower catchment samples, with overwhelming proportions of Mesozoic (250–66 Ma) ages (67%). In this sample, despite Jurassic crystalline rocks accounting for only 7% of its total upstream drainage area, Jurassic (145–200 Ma) ages constitute almost half (45%) of the results (96 of 212 analyses). Furthermore, while Cenozoic rocks make up 72% of exposed geologic units in the catchment, and more specifically 50% Pliocene–Quaternary volcanic rocks, Cenozoic zircons are nearly absent in this sample.

DISCUSSION

Provenance Signatures of Tectonic Provinces

The detrital zircon U-Pb results presented here for modern river sands enable discrimination of diagnostic sediment source regions that correspond to various tectonic provinces in the Ecuadorian segment of the northern Andes. Specifically, the U-Pb age distributions clearly identify the geochronological fingerprints of the major tectonic provinces recognized in cordilleran-type orogenic systems along ocean-continental subduction margins—namely the forearc, magmatic arc, and hinterland and foreland segments of the retroarc fold-thrust belt.

The U-Pb results from modern rivers in the west indicate that catchments spanning the Ecuadorian forearc and Western Cordillera are dominated by Neogene–Quaternary age populations, as expected for rocks related to the Andean magmatic arc. In contrast, the hinterland segment of the Andean fold-thrust belt (Eastern Cordillera) is defined by Proterozoic–Paleozoic age signatures diagnostic of the exposed Mesozoic sedimentary and metasedimentary rocks of ultimate cratonic origin that define this elevated region. Farther east, in frontal segments of the fold-thrust belt, rivers draining the Subandean Zone to Oriente Basin show distinctive bimodal Jurassic age populations, a secondary Triassic population, and subordinate Early Cretaceous ages representative of the major exposures of Mesozoic plutonic and metamorphic bedrock.

Paleozoic and Precambrian age signatures are generally restricted to and characteristic of the Eastern Cordillera, the elevated hinterland segment of the Andean fold-thrust belt in Ecuador (Fig. 2). As mentioned, most Precambrian ages, particularly those >900 Ma, were originally sourced from South American cratonic provinces to the east (Figs. 1 and 2). Therefore, the >900 Ma zircons are direct signatures of Eastern Cordillera detritus, indicating that individual grains have a long history involving initial derivation from the craton, Paleozoic–Mesozoic deposition in basin areas that included the present Eastern Cordillera, followed by Cenozoic erosion and transport within modern rivers.

Collectively, these results underscore the broad applicability of detrital zircon U-Pb methods to the northern Andes in general, and the Ecuadorian Andes in particular, for assessing the introduction and exhumation of sediment source regions and changes in regional paleodrainage patterns (e.g., Horton et al., 2010a, 2010b, 2015; Gutiérrez et al., 2019). These data provide a baseline for future studies of the Mesozoic–Cenozoic tectonic evolution of the northern Andes, as well as source-to-sink sediment routing and identification of critical provenance signatures within sedimentary basin systems within northern South America (Horton et al., 2015; Anderson et al., 2016; Horton, 2018).

Complex Downstream Variations in Source Signals

Variations in detrital zircon U-Pb age distributions of river sands from the large Rio Pastaza catchment reveal complexity in downstream signal propagation (Fig. 10). The characteristic Pliocene–Quaternary age population of the magmatic arc region in the upper catchment is lost at the expense of hinterland signatures in the middle catchment. In the lower catchment of the Subandean Zone, hinterland signatures are subdued at the expense of Subandean signatures, specifically Jurassic ages, even though Jurassic rocks account for a rather small proportion (7%) of rock units within the catchment area that is represented by the most-downstream sample.

Figure 10.

(A–D) Maps of the Rio Pastaza catchment area (thick black line) showing five river sand samples (color circles 4, 9, 10, 14, 15) relative to: (A) simplified surface geology and major faults (Litherland et al., 1993; Pratt et al., 2005; Bès de Berc et al., 2005); (B) digital elevation model (DEM) topography; (C) surface slope (in degrees); and (D) mean multiannual precipitation (in m/yr) for years 1965–1999 (INAMHI, 2014) (PRECIP—precipitation). (E) Longitudinal profile of the Rio Pastaza, from the upstream confluence of the Rio Chambo and Rio Patate (left) to the easternmost sample downstream (right), created from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) DEM V002 (https://asterweb.jpl.nasa.gov/gdem.asp) with 30 m resolution.

Figure 10.

(A–D) Maps of the Rio Pastaza catchment area (thick black line) showing five river sand samples (color circles 4, 9, 10, 14, 15) relative to: (A) simplified surface geology and major faults (Litherland et al., 1993; Pratt et al., 2005; Bès de Berc et al., 2005); (B) digital elevation model (DEM) topography; (C) surface slope (in degrees); and (D) mean multiannual precipitation (in m/yr) for years 1965–1999 (INAMHI, 2014) (PRECIP—precipitation). (E) Longitudinal profile of the Rio Pastaza, from the upstream confluence of the Rio Chambo and Rio Patate (left) to the easternmost sample downstream (right), created from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) DEM V002 (https://asterweb.jpl.nasa.gov/gdem.asp) with 30 m resolution.

The bimodal Jurassic U-Pb age distribution seen in downstream samples matches the ages of Jurassic igneous units exposed in the catchment, including the variably foliated Late Jurassic Azafrán plutonic complex within the Eastern Cordillera and the nonmetamorphosed Early Jurassic Abitagua Granite in the Subandean Zone (Litherland et al., 1994; Pratt et al., 2005; Cochrane 2013; Spikings et al., 2015). High surface slopes and high topographic relief across Paleozoic–Mesozoic igneous and metamorphic rock exposures in the Eastern Cordillera and proximal (western) Subandean Zone suggest that large, focused contributions of clastic material from these crystalline units result in their overrepresentation within U-Pb age distributions and the dilution of competing age signatures (Fig. 10C). In the middle catchment, the Rio Pastaza deeply and narrowly incises the Eastern Cordillera, with several large knickpoints with over 60–80 m of vertical relief (Fig. 10) and a steep stream gradient (22 m/km on average), suggesting that intense, focused erosion of rock units in this reach is sufficient to drown out the Pliocene–Quaternary zircon age signatures of the upper catchment.

Focused incision and erosion could be a product of spatially concentrated precipitation and/or active tectonic processes. Tectonic activity in the Subandean Zone since the Pleistocene has caused significant changes to drainage patterns, including several diversions of the Rio Pastaza by thrust-related fold growth and uplift of the Abitagua Granite (Fig. 10), as well as back-tilting due to fold growth above frontal structures farther east. 14C dating of fluvial terraces at various locations along the Rio Pastaza indicate rapid and potentially increasing rates of incision within hinterland and Subandean regions, with minimum average rates ranging from 5.0 to 6.7 mm/yr since ca. 18,000 yr B.P. (Bès de Berc et al., 2005). Although the highest rates may be influenced by dating of the youngest terraces (e.g., Gallen et al., 2015), on the basis of spatial relationships to active or recently active structures, Bès de Berc et al. (2005) argued that terrace development is mostly governed by active tectonics rather than climatic changes.

In considering downstream variations in U-Pb age distributions, it is possible that the Jurassic igneous units have high zircon concentrations or enhanced “zircon fertility” (e.g., Moecher and Samson, 2006; Dickinson, 2008) relative to Cenozoic sedimentary and Paleozoic–Cretaceous metamorphic units in the catchment (Fig. 10A), thereby contributing disproportionately large amounts of zircon. However, this is unlikely to be the singular cause for the observed dilution, on the basis of U-Pb age results from the Rio Zamora catchment to the south, where ∼50% of the catchment area is occupied by the Early Jurassic nonmetamorphosed Zamora Granite, a southern equivalent to the Abitagua Granite exposed in the Pastaza catchment (Spikings et al., 2015). A sample from the Rio Zamora (sample EC1572) does contain a large Jurassic U-Pb age population, but the age distribution is not as diluted by Jurassic zircons as would be expected based on the large surface area of the igneous body in the catchment if zircon fertility was the only contributing factor.

Another possible factor of disproportionate sediment contribution is the extreme spatial heterogeneity in mean annual precipitation between the eastern Amazonian lowlands and the high Andes (Fig. 10). Easterly trade winds carrying considerable moisture over the Amazon Basin are blocked by the topography of the Andes, forcing adiabatic expansion, condensation, and rainout, creating rainfall hot spots along the eastern flank of the Andes (e.g., Bookhagen and Strecker, 2008; Chavez and Takahashi, 2017). This wringing out of air masses at the mountain front leads to a sharp precipitation gradient in the Rio Pastaza catchment, with wet conditions (>3 m/yr) along the lower reaches of the catchment (<2 km elevation) and semiarid conditions in the upper Rio Pastaza valley (>3 km elevation), which receives <0.75 m/yr (Fig. 10; INAMHI, 2014). This sharp divide in mean annual precipitation occurs at 2–3 km elevation and coincides with the exposure levels of the Late Jurassic Azafrán pluton (Eastern Cordillera) and a precipitation maximum within the Early Jurassic Abitagua Granite (Subandean Zone). This spatially variable climatic framework may contribute to concentrated erosion within the fold-thrust belt and disproportionately large volumetric contributions from these Jurassic granites (Fig. 10).

A further consideration is the possible influence of temporally variable climatic phenomena such as El Niño–Southern Oscillation (ENSO; average reoccurrence between 2 and 7 yr) on modern river detrital zircon age distributions. Historically, La Niña events have been correlated with increased sediment and water discharge on the Rio Pastaza beyond background variability (Bernal et al., 2012); such flooding events could have led to further heterogeneous erosion along the stream channel, or rare mobilization of sediment sources in upstream reaches of the catchment.

Finally, are the number of analyses sufficient to capture all of the age populations present in downstream segments of the Rio Pastaza? Considering a commonly used threshold in detrital zircon provenance studies, whereby a sample size of 117 individual grain ages is considered to produce 95% confidence that populations >5% are identified (Vermeesch, 2004), two Rio Pastaza samples fall below this threshold—Eastern Cordillera sample 9: EC1510 with 105 grains, and sample 10: EC1513 with 91 grains (Fig. 9). However, the first of these (sample 9: EC1510) yields 95% confidence that populations >10% are identified (threshold of n = 103 grains), and the second, that populations >20% are identified at the same confidence level (threshold of n = 89) (Vermeesch, 2004). Considering that the upstream sample age distribution (Fig. 9, sample 4: EC1509) contains 92% Neogene ages (23–0 Ma)—which is much greater than a 20% fraction—we consider it highly unlikely that minor variations in the number of grains analyzed can explain the downstream loss of critical age populations.

Overall, drastic downstream changes in U-Pb age distributions along the Rio Pastaza provide a cautionary example whereby if it is not explicitly known that these samples are derived from the same large river system, they would likely be interpreted as having drastically different provenance.

CONCLUSIONS

Detrital zircon U-Pb geochronological results for 15 new samples (a total of 1705 analyses) of modern unconsolidated sand characterize the detrital fingerprints of 12 Andean rivers across the forearc, magmatic arc, and fold-thrust belt of Ecuador (1°N–5°S). These results show that river sediments from different tectonic provinces have diagnostic U-Pb age distributions defined by characteristic age populations. (1) In the forearc to magmatic arc (Western Cordillera), U-Pb age distributions are characterized by Cenozoic populations with an overall dominance of Pliocene–Quaternary ages. (2) In contrast, the hinterland segment of the fold-thrust belt (Eastern Cordillera) has a diverse (cosmopolitan) U-Pb age distribution characterized by large relative proportions of Proterozoic and Paleozoic ages representing upper Paleozoic–Mesozoic sedimentary and pelitic metamorphic rocks of ultimate cratonic derivation, and a near absence of Neogene–Quaternary ages. (3) For the frontal segment of the fold-thrust belt (Subandean Zone to Oriente Basin), U-Pb age distributions show limited diversity restricted to the Mesozoic, and are characterized by a large and bimodal Jurassic population with Early and Late Jurassic age peaks, a secondary cosmopolitan Triassic population, and subordinate Early Cretaceous ages.

Although samples from the Rio Pastaza headwaters yield age distributions consistent with all three tectonic provinces spanning the upper reaches of this catchment, these characteristic signatures are lost downstream as the river is overwhelmed by contributions from local sources. This downstream dilution of age distributions in the large Rio Pastaza watershed may be caused by focused erosion of fertile crystalline rocks due to an interplay between active tectonics and orographically induced precipitation focused at the range front.

ACKNOWLEDGMENTS

This research was supported by a U.S. National Science Foundation (NSF) Graduate Research Fellowship (grant 2016197896) awarded to LJJ, and NSF grant EAR-1338694 and National Geographic Society grant 9909-16 awarded to BKH. We thank E. Gabriela Gutiérrez, Sarah W.M. George, and Bernardo O. Beate for assistance in the field and insightful conversations. We thank reviewers R. Cecil, A. Laskowski, and Associate Editor B. Yanites for comments that helped improve the manuscript.

1Supplemental Material. Table S1: U-Pb geochronological results for modern river sand samples of Ecuador. Table S2: Sample information for modern river sand samples of Ecuador. Figure S1: NMF model results for modern river sand samples of Ecuador. Table S3: Grain-size measurements for modern river sand samples of the Rio Pastaza, Ecuador. Please visit https://doi.org/10.1130/GES02126.S1 or access the full-text article on www.gsapubs.org to view the Supplemental Material.
Science Editor: Andrea Hampel
Associate Editor: Brian J. Yanites

REFERENCES CITED

Amidon
,
W.H.
,
Burbank
,
D.W.
, and
Gehrels
,
G.E.
,
2005
,
Construction of detrital mineral populations: Insights from mixing of U-Pb zircon ages in Himalayan rivers
:
Basin Research
 , v.
17
, p.
463
485
, https://doi.org/10.1111/j.1365-2117.2005.00279.x.
Anderson
,
V.J.
,
Horton
,
B.K.
,
Saylor
,
J.E.
,
Mora
,
A.
,
Tesón
,
E.
,
Breecker
,
D.O.
, and
Ketcham
,
R.A.
,
2016
,
Andean topographic growth and basement uplift in southern Colombia: Implications for the evolution of the Magdalena, Orinoco, and Amazon river systems
:
Geosphere
 , v.
12
, p.
1235
1256
, https://doi.org/10.1130/GES01294.1.
Aspden
,
J.A.
, and
Litherland
,
M.
,
1992
,
The geology and Mesozoic collisional history of the Cordillera Real, Ecuador
:
Tectonophysics
 , v.
205
, p.
187
204
, https://doi.org/10.1016/0040-1951(92)90426-7.
Aspden
,
J.A.
,
Fortey
,
N.
,
Litherland
,
M.
,
Viteri
,
F.
, and
Harrison
,
S.M.
,
1992
,
Regional S-type granites in the Ecuadorian Andes: Possible remnants of the breakup of western Gondwana
:
Journal of South American Earth Sciences
 , v.
6
, p.
123
132
, https://doi.org/10.1016/0895-9811(92)90002-G.
Baby
,
P.
,
Rivadeneira
,
M.
,
Barragán
,
R.
, and
Christophoul
,
F.
,
2013
,
Thick-skinned tectonics in the Oriente foreland basin of Ecuador
, in
Norris
,
S.
,
Neeft
,
E.A.C.
, and
Geet
,
M.
, eds.,
Multiple Roles of Clays in Radioactive Waste Confinement
 :
Geological Society of London Special Publication
377
, p.
59
76
, https://doi.org/10.1144/SP377.1.
Bahlburg
,
H.
,
Vervoort
,
J.D.
,
Du Frane
,
S.A.
,
Bock
,
B.
,
Augustsson
,
C.
, and
Reimann
,
C.
,
2009
,
Timing of crust formation and recycling in accretionary orogens: Insights learned from the western margin of South America
:
Earth-Science Reviews
 , v.
97
, p.
215
241
, https://doi.org/10.1016/j.earscirev.2009.10.006.
Baldock
,
J.W.
,
1982
,
Geology of Ecuador (Explanatory Bulletin of the National Geological Map of the Republic of Ecuador)
: Quito, Ecuador,
Ministerio de Resursos Naturales y Energeticos, Direccion General de Geologia y Minas
,
59
p.
Barberi
,
F.
,
Coltelli
,
M.
,
Ferrara
,
G.
,
Innocenti
,
F.
,
Navarro
,
J.M.
, and
Santacroce
,
R.
,
1988
,
Plio-Quaternary volcanism in Ecuador
:
Geological Magazine
 , v.
125
, p.
1
14
, https://doi.org/10.1017/S0016756800009328.
Barragán
,
R.
,
Christophoul
,
F.
,
White
,
H.
,
Baby
,
P.
,
Rivadeneira
,
M.
,
Ramirez
,
F.
, and
Rodas
,
J.
,
2004
,
Estratigrafía secuencial del Cretácico de la Cuenca Oriente del Ecuador
, in
Baby
,
P.
,
Rivadeneira
,
M.
, and
Barragán
,
R.
, eds.,
La Cuenca Oriente: Geología y Petróleo: Travaux de l’Institut Français des Études Andine 144
 , p.
45
68
.
Bernal
,
C.
,
Christophoul
,
F.
,
Darrozes
,
J.
,
Soula
,
J.-C.
,
Baby
,
P.
, and
Burgos
,
J.
,
2011
,
Late Glacial and Holocene avulsions of the Rio Pastaza Megafan (Ecuador-Peru): Frequency and controlling factors
:
International Journal of Earth Sciences
 , v.
100
, p.
1759
1782
, https://doi.org/10.1007/s00531-010-0555-9.
Bernal
,
C.
,
Christophoul
,
F.
,
Soula
,
J.-C.
,
Darrozes
,
J.
,
Bourrel
,
L.
,
Laraque
,
A.
,
Burgos
,
J.
,
Bès de Berc
,
S.
, and
Baby
,
P.
,
2012
,
Gradual diversions of the Rio Pastaza in the Ecuadorian piedmont of the Andes from 1906 to 2008: Role of tectonics, alluvial fan aggradation, and ENSO events
:
International Journal of Earth Sciences
 , v.
101
, p.
1913
1928
, https://doi.org/10.1007/s00531-012-0752-9.
Bernard
,
B.
and
Andrade
,
D.
,
2011
,
Volcanes Cuaternarios del Ecuador Continental
: Quito, Ecuador,
Instituto Geofísico de la Escuela Politécnica Nacional
, informative poster.
Bès de Berc
,
S.B.
,
Soula
,
J.C.
,
Baby
,
P.
,
Souris
,
M.
,
Christophoul
,
F.
, and
Rosero
,
J.
,
2005
,
Geomorphic evidence of active deformation and uplift in a modern continental wedge-top–foredeep transition: Example of the eastern Ecuadorian Andes
:
Tectonophysics
 , v.
399
, p.
351
380
, https://doi.org/10.1016/j.tecto.2004.12.030.
Bettencourt
,
J.S.
,
Leite
,
W.B.
, Jr
,
Ruiz
,
A.S.
,
Matos
,
R.
,
Payolla
,
B.L.
, and
Tosdal
,
R.M.
,
2010
,
The Rondonian–San Ignacio Province in the SW Amazonian Craton: An overview
:
Journal of South American Earth Sciences
 , v.
29
, p.
28
46
, https://doi.org/10.1016/j.jsames.2009.08.006.
Blum
,
M.
, and
Pecha
,
M.
,
2014
,
Mid-Cretaceous to Paleocene North American drainage reorganization from detrital zircons
:
Geology
 , v.
42
, p.
607
610
, https://doi.org/10.1130/G35513.1.
Bookhagen
,
B.
, and
Strecker
,
M.R.
,
2008
,
Orographic barriers, high‐resolution TRMM rainfall, and relief variations along the eastern Andes
:
Geophysical Research Letters
 , v.
35
,
L06403
, https://doi.org/10.1029/2007GL032011.
Capaldi
,
T.N.
,
Horton
,
B.K.
,
McKenzie
,
N.R.
,
Stockli
,
D.F.
, and
Odlum
,
M.L.
,
2017
,
Sediment provenance in contractional orogens: The detrital zircon record from modern rivers in the Andean fold-thrust belt and foreland basin of western Argentina
:
Earth and Planetary Science Letters
 , v.
479
, p.
83
97
, https://doi.org/10.1016/j.epsl.2017.09.001.
Cardona
,
A.
,
Chew
,
D.
,
Valencia
,
V.A.
,
Bayona
,
G.
,
Mišković
,
A.
, and
Ibañez-Mejía
,
M.
,
2010
,
Grenvillian remnants in the northern Andes: Rodinian and Phanerozoic paleogeogaphic perspectives
:
Journal of South American Earth Sciences
 , v.
29
, p.
92
104
, https://doi.org/10.1016/j.jsames.2009.07.011.
Cawood
,
P.A.
,
Hawkesworth
,
C.J.
, and
Dhuime
,
B.
,
2012
,
Detrital zircon record and tectonic setting
:
Geology
 , v.
40
, p.
875
878
, https://doi.org/10.1130/G32945.1.
Chavez
,
S.P.
, and
Takahashi
,
K.
,
2017
,
Orographic rainfall hot spots in the Andes‐Amazon transition according to the TRMM precipitation radar and in situ data
:
Journal of Geophysical Research: Atmospheres
 , v.
122
, p.
5870
5882
, https://doi.org/10.1002/2016JD026282.
Chew
,
D.M.
,
Schaltegger
,
U.
,
Košler
,
J.
,
Whitehouse
,
M.J.
,
Gutjahr
,
M.
,
Spikings
,
R.A.
, and
Miškovic
,
A.
,
2007
,
U-Pb geochronologic evidence for the evolution of the Gondwanan margin of the north-central Andes
:
Geological Society of America Bulletin
 , v.
119
, p.
697
711
, https://doi.org/10.1130/B26080.1.
Chew
,
D.M.
,
Magna
,
T.
,
Kirkland
,
C.L.
,
Mišković
,
A.
,
Cardona
,
A.
,
Spikings
,
R.
, and
Schaltegger
,
U.
,
2008
,
Detrital zircon fingerprint of the Proto-Andes: Evidence for a Neoproterozoic active margin?
:
Precambrian Research
 , v.
167
, p.
186
200
, https://doi.org/10.1016/j.precamres.2008.08.002.
Christophoul
,
F.
,
Baby
,
P.
, and
Dávila
,
C.
,
2002
,
Stratigraphic responses to a major tectonic event in a foreland basin: The Ecuadorian Oriente Basin from Eocene to Oligocene times
:
Tectonophysics
 , v.
345
, p.
281
298
, https://doi.org/10.1016/S0040-1951(01)00217-7.
Cochrane
,
R.S.J.
,
2013
,
U-Pb thermochronology, geochronology and geochemistry of NW South America: Rift to drift transition, active margin dynamics and implications for the volume balance of continents [Ph.D. thesis]
:
Geneva, Switzerland
,
University of Geneva
,
191
p., https://doi.org/10.13097/archive-ouverte/unige:30029.
Cooper
,
M.A.
,
Addison
,
F.T.
,
Alvarez
,
R.
,
Coral
,
M.
,
Graham
,
R.H.
,
Hayward
,
A.B.
,
Howe
,
S.
,
Martinez
,
J.
,
Naar
,
J.
,
Peñas
,
R.
,
Pulham
,
A.J.
, and
Taborda
,
A.
,
1995
,
Basin development and tectonic history of the Llanos Basin, Eastern Cordillera, and middle Magdalena Valley, Colombia
:
American Association of Petroleum Geologists Bulletin
 , v.
79
, p.
1421
1442
.
Cordani
,
U.G.
, and
Teixeira
,
W.
,
2007
,
Proterozoic accretionary belts in the Amazonian Craton
, in
Hatcher
R.D.
, Jr
,
Carlson
,
M.P.
,
McBride
,
J.H.
, and
Martínez Catalán
,
J.R.
, eds.,
4–D Framework of Continental Crust
 :
Geological Society of America Memoir
200, p.
297
320
, https://doi.org/10.1130/2007.1200(14).
Cordani
,
U.G.
,
Cardona
,
A.
,
Jiménez
,
D.M.
,
Liu
,
D.
, and
Nutman
,
A.P.
,
2005
,
Geochronology of Proterozoic basement inliers in the Colombian Andes: Tectonic history of remnants of a fragmented Grenville belt
, in
Vaughan
,
A.P.M.
,
Leat
,
P.T.
, and
Pankhurst
,
R.J.
, eds.,
Terrane Processes at the Margins of Gondwana
 :
Geological Society of London Special Publication
246
, p.
329
346
, https://doi.org/10.1144/GSL.SP.2005.246.01.13.
Crowley
,
J.L.
,
Schoene
,
B.
, and
Bowring
,
S.A.
,
2007
,
U-Pb dating of zircon in the Bishop Tuff at the millennial scale
:
Geology
 , v.
35
, p.
1123
1126
, https://doi.org/10.1130/G24017A.1.
Dashwood
,
M.F.
, and
Abbotts
,
I.L.
,
1990
,
Aspects of the petroleum geology of the Oriente Basin, Ecuador
, in
Brooks
,
J.
, ed.,
Classic Petroleum Provinces
 :
Geological Society of London Special Publication
50
, p.
89
117
, https://doi.org/10.1144/GSL.SP.1990.050.01.06.
DeCelles
,
P.G.
,
Gehrels
,
G.E.
,
Quade
,
J.
,
LaReau
,
B.
, and
Spurlin
,
M.
,
2000
,
Tectonic implications of U-Pb zircon ages of the Himalayan orogenic belt in Nepal
:
Science
 , v.
288
, p.
497
499
, https://doi.org/10.1126/science.288.5465.497.
Dickinson
,
W.R.
,
2008
,
Impact of differential zircon fertility of granitoid basement rocks in North America on age populations of detrital zircons and implications for granite petrogenesis
:
Earth and Planetary Science Letters
 , v.
275
, p.
80
92
, https://doi.org/10.1016/j.epsl.2008.08.003.
Dickinson
,
W.R.
, and
Gehrels
,
G.E.
,
2003
,
U-Pb ages of detrital zircons from Permian and Jurassic eolian sandstones of the Colorado Plateau, USA: Paleogeographic implications
:
Sedimentary Geology
 , v.
163
, p.
29
66
, https://doi.org/10.1016/S0037-0738(03)00158-1.
Dickinson
,
W.R.
, and
Gehrels
,
G.E.
,
2008
,
U-Pb ages of detrital zircons in relation to paleogeography: Triassic paleodrainage networks and sediment dispersal across southwest Laurentia
:
Journal of Sedimentary Research
 , v.
78
, p.
745
764
, https://doi.org/10.2110/jsr.2008.088.
Ducea
,
M.N.
,
Giosan
,
L.
,
Carter
,
A.
,
Balica
,
C.
,
Stoica
,
A.M.
,
Roban
,
R.D.
,
Balintoni
,
I.
,
Filip
,
F.
, and
Petrescu
,
L.
,
2018
,
U-Pb detrital zircon geochronology of the lower Danube and its tributaries: Implications for the geology of the Carpathians
:
Geochemistry Geophysics Geosystems
 , v.
19
, p.
3208
3223
, https://doi.org/10.1029/2018GC007659.
Gallen
,
S.F.
,
Pazzaglia
,
F.J.
,
Wegmann
,
K.W.
,
Pederson
,
J.L.
, and
Gardner
,
T.W.
,
2015
,
The dynamic reference frame of rivers and apparent transience in incision rates
:
Geology
  v.
43
, p.
623
626
, https://doi.org/10.1130/G36692.1.
Garzanti
,
E.
,
Padoan
,
M.
,
Andò
,
S.
,
Resentini
,
A.
,
Vezzoli
,
G.
, and
Lustrino
,
M.
,
2013
,
Weathering and relative durability of detrital minerals in equatorial climate: Sand petrology and geochemistry in the East African Rift
:
The Journal of Geology
 , v.
121
, p.
547
580
, https://doi.org/10.1086/673259.
Gehrels
,
G.
,
2014
,
Detrital zircon U-Pb geochronology applied to tectonics
:
Annual Review of Earth and Planetary Sciences
 , v.
42
, p.
127
149
, https://doi.org/10.1146/annurev-earth-050212-124012.
Gehrels
,
G.
,
Valencia
,
V.
, and
Pullen
,
A.
,
2006
,
Detrital zircon geochronology by laser-ablation multicollector ICPMS at the Arizona LaserChron Center
:
The Paleontological Society Papers
 , v.
12
, p.
67
76
, https://doi.org/10.1017/S1089332600001352.
Gehrels
,
G.E.
,
Valencia
,
V.A.
, and
Ruiz
,
J.
,
2008
,
Enhanced precision, accuracy, efficiency, and spatial resolution of U‐Pb ages by laser ablation–multicollector–inductively coupled plasma–mass spectrometry
:
Geochemistry Geophysics Geosystems
 , v.
9
,
Q03017
, https://doi.org/10.1029/2007GC001805.
Gutiérrez
,
E.G.
,
Horton
,
B.K.
,
Vallejo
,
C.
,
Jackson
,
L.J.
, and
George
,
S.W.M.
,
2019
,
Provenance and geochronological insights into Late Cretaceous–Cenozoic foreland basin development in the Subandean Zone and Oriente Basin of Ecuador
, in
Horton
,
B.K.
, and
Folguera
,
A.
, eds.,
Andean Tectonics
 :
Amsterdam
,
Elsevier
, p.
237
268
, https://doi.org/10.1016/B978-0-12-816009-1.00011-3.
Hall
,
M.L.
, and
Wood
,
C.A.
,
1985
,
Volcano-tectonic segmentation of the northern Andes
:
Geology
 , v.
13
, p.
203
207
, https://doi.org/10.1130/0091-7613(1985)13<203:VSOTNA>2.0.CO;2.
Hall
,
M.L.
,
Samaniego
,
P.
,
Le Pennec
,
J.L.
, and
Johnson
,
J.B.
,
2008
,
Ecuadorian Andes volcanism: A review of Late Pliocene to present activity
:
Journal of Volcanology and Geothermal Research
 , v.
176
, p.
1
6
, https://doi.org/10.1016/j.jvolgeores.2008.06.012.
Higley
,
D.K.
,
2001
.
The Putumayo-Oriente-Maranon province of Columbia, Ecuador, and Peru—Mesozoic-Cenozoic and Paleozoic petroleum systems: U.S
. Geological Survey Digital Data Series 63, CD-ROM, https://doi.org/10.3133/ds63.
Horton
,
B.K.
,
2018
,
Sedimentary record of Andean mountain building
:
Earth-Science Reviews
 , v.
178
, p.
279
309
, https://doi.org/10.1016/j.earscirev.2017.11.025.
Horton
,
B.K.
,
Parra
,
M.
,
Saylor
,
J.E.
,
Nie
,
J.
,
Mora
,
A.
,
Torres
,
V.
,
Stockli
,
D.F.
, and
Strecker
,
M.R.
,
2010a
,
Resolving uplift of the northern Andes using detrital zircon age signatures
:
GSA Today
 , v.
20
, no.
7
, p.
4
10
, https://doi.org/10.1130/GSATG76A.1.
Horton
,
B.K.
,
Saylor
,
J.E.
,
Nie
,
J.
,
Mora
,
A.
,
Parra
,
M.
,
Reyes-Harker
,
A.
, and
Stockli
,
D.F.
,
2010b
,
Linking sedimentation in the northern Andes to basement configuration, Mesozoic extension, and Cenozoic shortening: Evidence from detrital zircon U-Pb ages, Eastern Cordillera, Colombia
:
Geological Society of America Bulletin
, v.
122
, p.
1423
1442
, https://doi.org/10.1130/B30118.1.
Horton
,
B.K.
,
Anderson
,
V.J.
,
Caballero
,
V.
,
Saylor
,
J.E.
,
Nie
,
J.
,
Parra
,
M.
, and
Mora
,
A.
,
2015
,
Application of detrital zircon U-Pb geochronology to surface and subsurface correlations of provenance, paleodrainage, and tectonics of the Middle Magdalena Valley Basin of Colombia
:
Geosphere
 , v.
11
, p.
1790
1811
, https://doi.org/10.1130/GES01251.1.
Ibañez-Mejia
,
M.
,
Pullen
,
A.
,
Pepper
,
M.
,
Urbani
,
F.
,
Ghoshal
,
G.
, and
Ibañez-Mejia
,
J.C.
,
2018
,
Use and abuse of detrital zircon U-Pb geochronology—A case from the Río Orinoco delta
,
eastern Venezuela
:
Geology
 , v.
46
, p.
1019
1022
, https://doi.org/10.1130/G45596.1.
INAMHI (Instituto Nacional de Meteorología e Hidrología)
,
2014
,
Mapa de precipitación media multianual—Serie 1965–1999
: Quito, Ecuador,
Instituto Nacional de Meteorología e Hidrología
, scale 1:1,000,000.
Jackson
,
L.J.
,
Horton
,
B.K.
,
Beate
,
B.O.
,
Bright
,
J.
, and
Breecker
,
D.O.
,
2019
,
Testing stable isotope paleoaltimetry with Quaternary volcanic glasses from the Ecuadorian Andes
:
Geology
 , v.
47
, p.
411
414
, https://doi.org/10.1130/G45861.1.
Jaillard
,
E.
,
Caron
,
M.
,
Dhondt
,
A.
,
Ordoñez
,
M.
,
Andrade
,
R.
,
Bengtson
,
P.
,
Bulot
,
L.
,
Cappetta
,
H.
,
Dávila
,
C.
,
Díaz
,
R.
and
Huacho
,
J.
,
1997
,
Síntesis estratigráfica y sedimentológica del Cretáceo y Paleógeno de la cuenca oriental del Ecuador
: Quito, Ecuador,
Orstom-Petroproduccion Publication
,
164
p.
Jaillard
,
E.
,
Ordonez
,
M.
,
Suarez
,
J.
,
Toro
,
J.
,
Iza
,
D.
, and
Lugo
,
W.
,
2004
,
Stratigraphy of the late Cretaceous–Paleogene deposits of the Cordillera Occidental of central Ecuador: Geodynamic implications
:
Journal of South American Earth Sciences
 , v.
17
, p.
49
58
, https://doi.org/10.1016/j.jsames.2004.05.003.
Kerr
,
A.C.
,
Aspden
,
J.A.
,
Tarney
,
J.
, and
Pilatasig
,
L.F.
,
2002
,
The nature and provenance of accreted oceanic terranes in western Ecuador: Geochemical and tectonic constraints
:
Journal of the Geological Society
 , v.
159
, p.
577
594
, https://doi.org/10.1144/0016-764901-151.
Laskowski
,
A.K.
,
DeCelles
,
P.G.
, and
Gehrels
,
G.E.
,
2013
,
Detrital zircon geochronology of Cordilleran retroarc foreland basin strata, western North America
:
Tectonics
 , v.
32
, p.
1027
1048
, https://doi.org/10.1002/tect.20065.
Litherland
,
M.
,
Zamora
,
A.
,
Eguez
,
A.
,
Aspden
,
J.
,
Baez
,
N.
,
Beate
,
B.
,
Benitez
,
S.
,
Jaillard
,
E.
, and
Van Thournout
,
F.
,
1993
,
Mapa geológico de la República del Ecuador
:
British Geological Survey, scale
 
1
:1,000,000.
Litherland
,
M.
,
Aspden
,
J.A.
, and
Jemielita
,
R.A.
,
1994
,
The Metamorphic Belts of Ecuador
:
British Geological Survey Overseas Memoir
11
,
147
p.
Malusà
,
M.G.
,
Resentini
,
A.
, and
Garzanti
,
E.
,
2016
,
Hydraulic sorting and mineral fertility bias in detrital geochronology
:
Gondwana Research
 , v.
31
, p.
1
19
, https://doi.org/10.1016/j.gr.2015.09.002.
Martin‐Gombojav
,
N.
, and
Winkler
,
W.
,
2008
,
Recycling of Proterozoic crust in the Andean Amazon foreland of Ecuador: Implications for orogenic development of the Northern Andes
:
Terra Nova
 , v.
20
, p.
22
31
, https://doi.org/10.1111/j.1365-3121.2007.00782.x.
Moecher
,
D.P.
, and
Samson
,
S.D.
,
2006
,
Differential zircon fertility of source terranes and natural bias in the detrital zircon record: Implications for sedimentary provenance analysis
:
Earth and Planetary Science Letters
 , v.
247
, p.
252
266
, https://doi.org/10.1016/j.epsl.2006.04.035.
Nesbitt
,
H.W.
,
Young
,
G.M.
,
McLennan
,
S.M.
, and
Keays
,
R.R.
,
1996
,
Effects of chemical weathering and sorting on the petrogenesis of siliciclastic sediments, with implications for provenance studies
:
The Journal of Geology
 , v.
104
, p.
525
542
, https://doi.org/10.1086/629850.
Pepper
,
M.
,
Gehrels
,
G.
,
Pullen
,
A.
,
Ibanez-Mejia
,
M.
,
Ward
,
K.M.
, and
Kapp
,
P.
,
2016
,
Magmatic history and crustal genesis of western South America: Constraints from U-Pb ages and Hf isotopes of detrital zircons in modern rivers
:
Geosphere
 , v.
12
, p.
1532
1555
, https://doi.org/10.1130/GES01315.1.
Pratt
,
W.T.
,
Duque
,
P.
, and
Ponce
,
M.
,
2005
,
An autochthonous geological model for the eastern Andes of Ecuador
:
Tectonophysics
 , v.
399
, p.
251
278
, https://doi.org/10.1016/j.tecto.2004.12.025.
Romeuf
,
N.
,
Aguirre
,
L.
,
Soler
,
P.
,
Féraud
,
G.
,
Jaillard
,
E.
, and
Ruffet
,
G.
,
1995
,
Middle Jurassic volcanism in the Northern and Central Andes
:
Revista Geológica de Chile
 , v.
22
, p.
245
259
.
Ruiz
,
G.M.H.
,
Seward
,
D.
, and
Winkler
,
W.
,
2007
,
Evolution of the Amazon Basin in Ecuador with special reference to hinterland tectonics: Data from zircon fission-track and heavy mineral analysis
:
Developments in Sedimentology
 , v.
58
, p.
907
934
, https://doi.org/10.1016/S0070-4571(07)58036-2.
Santos
,
J.O.S.
,
Rizzotto
,
G.J.
,
Potter
,
P.E.
,
McNaughton
,
N.J.
,
Matos
,
R.S.
,
Hartmann
,
L.A.
,
Chemale
,
F.
, Jr
, and
Quadros
,
M.E.S.
,
2008
,
Age and autochthonous evolution of the Sunsás Orogen in West Amazon Craton based on mapping and U-Pb geochronology
:
Precambrian Research
 , v.
165
, p.
120
152
, https://doi.org/10.1016/j.precamres.2008.06.009.
Saylor
,
J.E.
, and
Sundell
,
K.E.
,
2016
,
Quantifying comparison of large detrital geochronology data sets
:
Geosphere
 , v.
12
, p.
203
220
, https://doi.org/10.1130/GES01237.1.
Saylor
,
J.E.
,
Knowles
,
J.N.
,
Horton
,
B.K.
,
Nie
,
J.
, and
Mora
,
A.
,
2013
,
Mixing of source populations recorded in detrital zircon U-Pb age spectra of modern river sands
:
The Journal of Geology
 , v.
121
, p.
17
33
, https://doi.org/10.1086/668683.
Saylor
,
J.E.
,
Sundell
,
K.E.
and
Sharman
,
G.R.
,
2019
,
Characterizing sediment sources by non-negative matrix factorization of detrital geochronological data
:
Earth and Planetary Science Letters
 , v.
512
, p.
46
58
, https://doi.org/10.1016/j.epsl.2019.01.044.
Schütte
,
P.
,
Chiaradia
,
M.
, and
Beate
,
B.
,
2010
,
Geodynamic controls on Tertiary arc magmatism in Ecuador: Constraints from U-Pb zircon geochronology of Oligocene–Miocene intrusions and regional age distribution trends
:
Tectonophysics
 , v.
489
, p.
159
176
, https://doi.org/10.1016/j.tecto.2010.04.015.
Sharman
,
G.R.
,
Graham
,
S.A.
,
Grove
,
M.
,
Kimbrough
,
D.L.
, and
Wright
,
J.E.
,
2015
,
Detrital zircon provenance of the Late Cretaceous–Eocene California forearc: Influence of Laramide low-angle subduction on sediment dispersal and paleogeography
:
Geological Society of America Bulletin
 , v.
127
, p.
38
60
, https://doi.org/10.1130/B31065.1.
Sharman
,
G.R.
,
Sharman
,
J.P.
, and
Sylvester
,
Z.
,
2018
,
detritalPy: A Python‐based toolset for visualizing and analysing detrital geo‐thermochronologic data
:
The Depositional Record
 , v.
4
, p.
202
215
, https://doi.org/10.1002/dep2.45.
Sickmann
,
Z.T.
,
Schwartz
,
T.M.
, and
Graham
,
S.A.
,
2018
,
Refining stratigraphy and tectonic history using detrital zircon maximum depositional age: An example from the Cerro Fortaleza Formation, Austral Basin, southern Patagonia
:
Basin Research
 , v.
30
, p.
708
729
, https://doi.org/10.1111/bre.12272.
Spencer
,
C.J.
,
Kirkland
,
C.L.
, and
Roberts
,
N.M.
,
2018
,
Implications of erosion and bedrock composition on zircon fertility: Examples from South America and Western Australia
:
Terra Nova
 , v.
30
, p.
289
295
, https://doi.org/10.1111/ter.12338.
Spikings
,
R.A.
,
Winkler
,
W.
,
Seward
,
D.
, and
Handler
,
R.
,
2001
,
Along-strike variations in the thermal and tectonic response of the continental Ecuadorian Andes to the collision with heterogeneous oceanic crust
:
Earth and Planetary Science Letters
 , v.
186
, p.
57
73
, https://doi.org/10.1016/S0012-821X(01)00225-4.
Spikings
,
R.
,
Cochrane
,
R.
,
Villagomez
,
D.
,
Van der Lelij
,
R.
,
Vallejo
,
C.
,
Winkler
,
W.
, and
Beate
,
B.
,
2015
,
The geological history of northwestern South America: From Pangaea to the early collision of the Caribbean large igneous province (290–75 Ma)
:
Gondwana Research
 , v.
27
, p.
95
139
, https://doi.org/10.1016/j.gr.2014.06.004.
Teixeira
,
W.
,
Tassinari
,
C.C.G.
,
Cordani
,
U.G.
, and
Kawashita
,
K.
,
1989
,
A review of the geochronology of the Amazonian Craton: Tectonic implications
:
Precambrian Research
 , v.
42
, p.
213
227
, https://doi.org/10.1016/0301-9268(89)90012-0.
Tschopp
,
H.J.
,
1953
,
Oil explorations in the Oriente of Ecuador
:
American Association of Petroleum Geologists Bulletin
 , v.
37
, p.
2303
2347
.
Vallejo
,
C.
,
Winkler
,
W.
,
Spikings
,
R.A.
,
Luzieux
,
L.
,
Heller
,
F.
, and
Bussy
,
F.
,
2009
,
Mode and timing of terrane accretion in the forearc of the Andes in Ecuador
, in
Kay
,
S.M.
,
Ramos
,
V.A.
, and
Dickinson
,
W.R.
, eds.,
Backbone of the Americas: Shallow Subduction, Plateau Uplift, and Ridge and Terrane Collision
 :
Geological Society of America Memoir 204
, p.
197
216
, https://doi.org/10.1130/2009.1204(09).
Vallejo
,
C.
,
Tapia
,
D.
,
Gaibor
,
J.
,
Steel
,
R.
,
Cardenas
,
M.
,
Winkler
,
W.
,
Valdez
,
A.
,
Esteban
,
J.
,
Figuera
,
M.
,
Leal
,
J.
, and
Cuenca
,
D.
,
2017
,
Geology of the Campanian M1 sandstone oil reservoir of eastern Ecuador: A delta system sourced from the Amazon Craton
:
Marine and Petroleum Geology
 , v.
86
, p.
1207
1223
, https://doi.org/10.1016/j.marpetgeo.2017.07.022.
Vallejo
,
C.
,
Spikings
,
R.A.
,
Horton
,
B.K.
,
Luzieux
,
L.
,
Romero
,
C.
,
Winkler
,
W.
, and
Thomsen
,
T.B.
,
2019
,
Late Cretaceous to Miocene stratigraphy and provenance of the coastal forearc and Western Cordillera of Ecuador: Evidence for accretion of a single oceanic plateau fragment
, in
Horton
,
B.K.
, and
Folguera
,
A.
, eds.,
Andean Tectonics
 :
Amsterdam
,
Elsevier
, p.
209
236
, https://doi.org/10.1016/B978-0-12-816009-1.00010-1.
Vermeesch
,
P.
,
2004
,
How many grains are needed for a provenance study?
:
Earth and Planetary Science Letters
 , v.
224
, p.
441
451
, https://doi.org/10.1016/j.epsl.2004.05.037.
Vermeesch
,
P.
,
2012
,
On the visualisation of detrital age distributions
:
Chemical Geology
 , v.
312–313
, p.
190
194
, https://doi.org/10.1016/j.chemgeo.2012.04.021.
Wang
,
C.
,
Liang
,
X.
,
Foster
,
D.A.
,
Tong
,
C.
,
Liu
,
P.
,
Liang
,
X.
, and
Zhang
,
L.
,
2018
,
Linking source and sink: Detrital zircon provenance record of drainage systems in Vietnam and the Yinggehai–Song Hong Basin, South China Sea
:
Geological Society of America Bulletin
, v.
131
, p.
191
204
, https://doi.org/10.1130/B32007.1.
Weislogel
,
A.L.
,
Graham
,
S.A.
,
Chang
,
E.Z.
,
Wooden
,
J.L.
,
Gehrels
,
G.E.
, and
Yang
,
H.
,
2006
,
Detrital zircon provenance of the Late Triassic Songpan-Ganzi complex: Sedimentary record of collision of the North and South China blocks
:
Geology
 , v.
34
, p.
97
100
, https://doi.org/10.1130/G21929.1.
Witt
,
C.
,
Rivadeneira
,
M.
,
Poujol
,
M.
,
Barba
,
D.
,
Beida
,
D.
,
Beseme
,
G.
, and
Montenegro
,
G.
,
2017
,
Tracking ancient magmatism and Cenozoic topographic growth within the Northern Andes forearc: Constraints from detrital U-Pb zircon ages
:
Geological Society of America Bulletin
, v.
129
, p.
415
428
, https://doi.org/10.1130/B31530.1.
Gold Open Access: This paper is published under the terms of the CC-BY-NC license.