SciELO - Scientific Electronic Library Online

 
vol.37 número2Lahares secundarios en el volcán Popocatépetl: El lahar Nexpayantla del 4 de febrero, 2010Caracterización de Necrosoles en el bosque seco tropical en Colombia: un acercamiento desde la geología forense índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Revista mexicana de ciencias geológicas

versión On-line ISSN 2007-2902versión impresa ISSN 1026-8774

Rev. mex. cienc. geol vol.37 no.2 Ciudad de México ago. 2020  Epub 17-Sep-2025

https://doi.org/10.22201/cgeo.20072902e.2020.2.1557 

Regular Articles

U-Pb geochronology of rutile: deciphering the cooling history of the Oaxacan Complex granulites, southern Mexico

Geocronología U-Pb en rutilo: una contribución para descifrar la historia de enfriamiento de granulitas del Complejo Oaxaqueño, sur de México

Miguel Gerardo Adame-Martínez1  # 

Luigi A. Solari2  * 

Carlos Ortega-Obregón2 

Fanis Abdullin2 

1Facultad de Ingeniería, Universidad Autónoma de San Luis Potosí, Dr. Manuel Nava # 8, Zona Universitaria Poniente, 78290, San Luis Potosí, S.L.P., Mexico.

2 Centro de Geociencias, Universidad Nacional Autónoma de México, Campus Juriquilla, 76230, Querétaro, Qro., Mexico.

# Current Address: Industrial Minera México, S.A. de C.V., Centro de Investigaciones Metalúrgicas, Centenario y Turquesa s/n, Col. La Prieta, C.P. 33800, Parral, Chihuahua, Mexico.


ABSTRACT

Rutile (TiO2) is a heavy mineral, commonly found as accessory in many lithologies, such as basic igneous rocks, high-grade metamorphic units, as well as a detritus in sedimentary clastic rocks. Its chemical composition is sensitive to the crystallization environment, allowing a characterization of either metabasic or metasedimentary protoliths in metamorphic rocks. Thanks to the capability to accept U in its crystalline network, at least in metasedimentary, high-grade protoliths, rutile can be dated by U-Pb geochronology. Furthermore, its closure temperature of ca. 600 °C for the U-Pb system makes rutile a suitable chronometer, complementary to zircon, to unravel provenance and exhumation paths in both sedimentary siliciclastic cover and basement units. Besides, the Zr-in thermometer allows for a very precise calculation of the rutile crystallization temperature. In the example case presented here, focused on granulite facies units of the Grenvillian Oaxacan Complex (OC), rutile crystallisation took place in the range 808-873 °C. Data for different localities indicate that cooling and exhumation after the Zapotecan granulite facies event (ca. 990 Ma) was heterogeneous among the different tectonic slices that constitute the OC. Cooling occurred in the central sector (Nochixtlán-Oaxaca) right after the granulite peak, with fast cooling rates of ca. 40 °C/Ma. To the north and south, the cooling to ca. 600 °C was much slower, with calculated cooling rates of ca. 3 °C/Ma for the northern OC outcrops in Coatepec (Puebla) to ca. 6 °C/Ma south of Ejutla (Oaxaca). This can be related to a combination of factors, such as an early collapse of some sectors of the orogen, a change of conditions in the subducing plate, or more in general, to a sudden change in the geodynamic conditions during the Zapotecan orogeny and Amazonia-Baltica amalgamation.

Key words: rutile; granulite metamorphism; U-Pb geochronology; chemical composition; Oaxacan Complex; Mexico

RESUMEN

El rutilo (TiO2) es un mineral pesado, accesorio común en muchas litologías tanto ígneas básicas como metamórficas de alto grado, o bien presente como detrito en rocas sedimentarias siliciclásticas Tiene unas importantes características petrogenéticas, como son la composición química sensible al ambiente de formación y una marcada y diferente afinidad en protolitos metabásicos y metasedimentarios. Gracias a su capacidad para aceptar U en su red cristalina, al menos en los protolitos metasedimentarios de alto grado metamórfico, el rutilo puede datarse por medio de geocronología por U-Pb. Además, siendo la temperatura de cierre para el sistema U-Pb en rutilo de aproximadamente 600 °C, el rutilo se convierte en un mineral muy útil para caracterizar el enfriamiento en rocas de alto grado metamórfico, complementario al zircón, para entonces poder descifrar procedencia y eventos de exhumación tanto en rocas sedimentarias siliciclásticas de cobertura, como en unidades del basamento. El termómetro de Zr en rutilo permite también un cálculo muy preciso de la temperatura de cristalización del rutilo. En el ejemplo presentado en este trabajo, que se enfoca al estudio de algunas unidades granulíticas colectadas en el Complejo Oaxaqueño del sur de México, la cristalización del rutilo ocurrió en el rango 808-873 °C. Las rocas metamórficas estudiadas, colectadas en diferentes localidades, permiten ver que el enfriamiento y exhumación del Complejo Oaxaqueño fue de hecho heterogéneo después del pico metamórfico Zapoteco (ca. 990 Ma) en las diferentes unidades tectónicas que lo componen. En la parte centro-norte (Nochixtlán-Oaxaca) el enfriamiento fue muy rápido, con una tasa de ca. 40 °C/Ma, mientras las porciones norte y sur enfriaron a una tasa mucho menor, de ca. 3 °C/Ma en el norte (al este de Coatepec, Pue.) y ca. 6 °C/Ma en el sur (entre Ejutla y Mihauatlán, Oax.). Esto puede deberse a una combinación de factores, como un colapso temprano de porciones del orógeno, un cambio de condiciones en la placa en subducción o, más en general, a un cambio repentino de las condiciones geodinámicas durante las primeras fases de la amalgamación de Rodinia.

Key words: rutilo; metamorfismo granulítico; geocronología U-Pb; composición química; Complejo Oaxaqueño; México

INTRODUCTION

The characterization of heavy minerals in siliciclastic rocks is a fundamental tool in modern geosciences, because it can unravel age, distributions, isotopic and chemical composition, facilitating the study of high-grade terrains, their tectonic evolution and exhumation history, or even yielding important constraints to provenance analysis (e.g., Morton and Hallsworth, 1999; Keppie et al., 2004; Schulze et al., 2004; Garzanti, 2016). The continuous improvement, and availability, of microanalytical techniques made the quantification of the composition of heavy minerals a standard procedure in many laboratories, taking advantage of reduced time and analytical costs. While zircon is considered the king of the heavy minerals used in provenance studies, thanks to its abundance as well as its chemical and mechanical strength (Park et al., 2010; Ortega-Flores et al., 2013; Spencer et al., 2014), in recent years other minerals, also quite abundant in detrital rocks, have been increasingly used as chemical and isotopic tracers. In particular, apatite and rutile have received a lot of interest in the geoscientific community (e.g., Luvizotto and Zack, 2008; Triebold et al., 2012; Abdullin et al., 2015; O’Sullivan et al., 2018; Rösel et al., 2019), because both serve as excellent petrologic witnesses, as well as geochronometers with closure temperatures lower than that of zircon in the U-Pb isotopic system (ca. 350-550 °C for apatite according to Cochrane et al., 2014, and ca. 550-600 °C for rutile, according to Cherniak, 2000 and Kooijman et al., 2012; compared with zircon, in excess of 900 °C, Cherniak and Watson, 2000). In high-grade metamorphic rocks, rutile provides excellent information on the post-peak exhumation history when U-Pb cooling ages are combined with thermometric calculations, based on the equilibrium between Ti and Zr (Watson et al., 2006; Tomkins et al., 2007).

In this paper we present new microanalytical data for rutile, obtained on separates from samples belonging to the Grenvillian, granulite-facies Oaxacan Complex in southern Mexico (Keppie et al., 2003; Solari et al., 2003, 2004b, 2014, 2020; Weber and Schulze, 2014), as an example of how the information extracted from this mineral can improve the understanding of high-grade rocks, their behavior during metamorphism, and exhumation history.

GEOLOGIC SETTING

The Oaxacan Complex (OC hereafter; Figure 1) constitutes the largest basement exposure of granulite facies, Grenvillian-age rocks in Mexico. Together with other minor known outcrops of similar rocks, it constitutes the framework of the Oaxaquia microcontinent, a Proterozoic piece that makes up the backbone of Mexico, as we know it today (Ortega-Gutiérrez et al., 1995; 2018). In more detail, the ca. 1.3-0.95 Ga OC is constituted by an assemblage of (mostly) metaplutonic and, in minor amount, metasedimentary rocks, all metamorphosed under granulite-facies conditions. High-grade metamorphism (ca. 750-800 °C and 7.5-8.5 kbar, Mora and Valley, 1985; Solari et al., 2004b; also see a review in Ortega-Gutiérrez et al., 2018) occurred during the Zapotecan Orogeny (Solari et al., 2003), at ca. 990 Ma. An igneous sequence of ca. 1.2 Ga (locally, up to 1.4 Ga), is either interpreted as indicative of an intra-oceanic arc or a rift (Ortega-Gutiérrez et al., 1995; Keppie et al., 2001 and 2003; Keppie and Dostal, 2007; Solari et al., 2014; Weber and Schulze, 2014; Ortega-Gutiérrez et al., 2018). A younger igneous pulse, dated at around 1.02 Ga, makes up an AMCG (anorthosite-mangerite-charnockite-granite) suite (Keppie et al., 2003). After the Zapotecan metamorphic peak, several post-tectonic, granitic to carbonated pegmatites intruded high-grade rocks of the OC. Their ages span from ca. 980 to ca. 965 Ma (Solari et al., 2003; Schepetilnikova et al., 2015), being related to magmatism generated by decompressional melting during or right after the Zapotecan event. The tectonic limits of the OC are: to the east, the Sierra de Juárez Mylonitic Complex as a tectonic limit with the Cuicateco terrane (Alaniz-Álvarez et al., 1994, 1996); to the west, the Caltepec fault zone with the Acatlán Complex (Elías-Herrera and Ortega-Gutiérrez, 2002; Elías-Herrera et al., 2007); to the south, the Chacalapa fault system against the Mesozoic-Tertiary Xolapa Complex (Tolson, 2005); to the north it is buried underneath Cenozoic continental deposits (Dávalos-Álvarez et al., 2007) and Quaternary continental deposits associated to the development of the Trans-Mexican Volcanic Belt.

Figure 1 a) General overview of the geology of Mexico, including main subdivision in tectonic domains. b) Schematic map of the Oaxacan Complex, southern Mexico, with sample localities. Modified after Solari et al. (2014)

Whereas the Mesoproterozoic to Neoproterozoic evolution of the OC has been tackled in several papers, the post-Zapotecan geologic history is known only locally. For instance, Solari et al. (2004a) studied some shear zones in the northern OC, describing their tectonic evolution and geologic significance during the Phanerozoic (mostly associating them with the geologic evolution from the Late Triassic to the Early Cretaceous, as a result of plate re-accommodation during the break up of Pangea and the opening of the Gulf of Mexico). Moreover, Keppie et al. (2004) reported some Tonian cooling ages (mostly by Ar-Ar), associating them to the switching between steep and flat-slab subduction and subsequent thermal relaxation of the lithosphere, and slow cooling regime during Cryogenian time that lead to the Oaxaquia breakup during the Ediacaran time (e.g., González-Guzmán et al., 2016). The construction of cooling curves would enormously benefit from more data obtained by different methods, and in minerals with different closure temperatures, collected at different tectonic levels, to correctly define the timing of exhumation and, eventually, the differences in the tectonic regime from one to another locality. The use of rutile thus helps to close one of those gaps, generating data that fit to characterize the post-Zapotecan exhumation of the Oaxacan Complex.

SAMPLE DESCRIPTION

After screening by petrography in thin section, as well as in heavy mineral concentrates, five samples were selected (geographic locations are illustrated in Figure 1).

The overall strike of the main OC structures bears, mostly, ENE-WSW, gently plunging toward NNW (Solari et al., 2003 and 2004a). At least for the northern OC, a stack of three tectonic sleeves was recognized and described by Solari et al. (2003). New field observations, not yet fully documented, suggest the continuation of piling up of tectonic sleeves along the whole OC. As an example for the existence of such tectonic sleeves that represent different depths in the crust, Schulze (2011) calculated for the southernmost OC sector of Pluma Hidalgo, higher metamorphic pressures of ca. 10 kb, possibly representative of a deeper tectonic slice when compared to the northern sectors studied by Solari et al. (2003). The samples selected for this work, which mostly scatter along the NW-SE length of the whole OC, are thus representative of different depths of those geologic units, now brought to the eroded topographic surface, but originally in an unknown, relative position.

Sample OC1101 is from a garnet-bearing granulite gneiss cropping out along the Federal Highway 190, along an EW Tertiary normal fault (Santa Lucia fault), in the uppermost tectonic slice described by Solari et al. (2003). It is made up of large (up to 5 mm) garnet, almost completely altered pyroxene, perthitic feldspar, quartz and plagioclase. In practice, it resembles an original, quartzo-feldsphatic sediment, such as a sandstone or a lithoarenite. Rutile, zircon and apatite are accessory minerals. Rutile is generally tabular to needle shaped, dark-brown colored, ranging in size up to 500 μm in length and up to 80 μm in width; some crystals show twinning (e.g., Figure 2a).

Figure 2 Scanning electron microscope (SEM) photomicrographs of some selected and mounted rutile crystals. The scale is indicated in each image. a) A twinned rutile crystal belonging to sample OC1101. b) A needle-shaped rutile crystal belonging to sample OC1006. c) Post-ablation, enlarged view of some mounted crystals belonging to sample OC1006. d) Pre-ablation image of some selected rutile crystals belonging to sample OC1016. 

Sample OC1006 is a banded and foliated granulitic gneiss, found south of Asunción Nochixtlán, constituting the host rock for charnockitic (ca. 1135±12 Ma, Solari et al., 2020) and syenitic (ca. 1221±8.2 Ma, Solari et al., 2020) intrusions. It bears quartz, perthitic feldspar, plagioclase (often sericitized), garnet, and sillimanite. The abundant aluminum-rich minerals are suggestive of a metasedimentary origin. Rutile is found either as euhedral crystals, or as inclusion needles in quartz. Their shapes are generally more elongated than in the previous sample, with an up to 1:13 elongation ratio (Figure 2b and 2c), ranging up to 1000 μm in length and up to 80 μm in width. Rutile is light colored, ranging from pale brown to honey brown.

Sample OC1014 is from a granulitic gneiss that occurs NW of Ejutla, whose paragenesis is indicative of a sedimentary protolith. It is made up of quartz, perthitic feldspar, plagioclase, abundant garnet, up to 7 mm in size, together with abundant graphite and secondary biotite. Rutile is an abundant accessory mineral, normally light to moderately brown colored, and up to 800 μm in length and up to 80 μm in width.

Sample OC1016 is from a granulitic orthogneiss cropping out East of Santa María Peñoles. It is faintly banded, and made up of quartz, micro-perthitic feldspar, and plagioclase. Some mafic minerals are completely replaced by an assemblage of oxide (magnetite and ilmenite), biotite and tremolite, tentatively indicating the original presence of pyroxene. Rutile is an accessory mineral together with zircon and minor apatite. It ranges from euhedral, short prismatic to stubby shapes, up to 800 μm in length, to subhedral shapes, in some cases less than 200 µm in length (Figure 2d). The maximum width in the observed crystals is 80 μm. It is also clearly darker than in the previous samples, ranging from dark brown to almost opaque, probably as a consequence of a higher iron content.

Sample OC1017 is a granulitic gneiss from an outcrop along the Federal Road 190, halfway from Ejutla to Mihauatlán. It contains abundant garnet, meso-perthitic feldspar, quartz, graphite, and sillimanite. Rutile and zircon are common accessory minerals. In this sample rutile is prismatic, sometimes with pyramidal terminations, ranging in size up to 1200 μm in length and with a maximum width of 80 μm, mostly reddish-brown colored; only few of the observed crystals are dark brown.

METHODOLOGY

Sample preparation

After the petrographic analysis aimed to confirm the presence of rutile, the selected samples were crushed and heavy minerals were separated using conventional techniques (e.g., Solari et al., 2007). The heavy minerals concentrate was further passed through a Frantz magnetic separator to eliminate magnetic minerals (with currents up to 2 A). The final concentrate of non-magnetic heavy minerals consisted of zircon, rutile, apatite, and eventually titanite. Among those phases, rutile is the only dark mineral, so hand picking was pretty straightforward under the binocular microscope. Approximately 120-150 rutile grains from each sample, representative of all morphologies, color and size, were selected and casted in a ring mount with epoxy resin, polished with sandpaper, and then imaged under a binocular microscope prior to analysis.

Rutile U-Pb analyses

U-Pb analyses were performed by laser ablation inductively-coupled plasma mass spectrometry (LA-ICPMS) employing a Thermo ICap Qc quadrupole ICPMS, coupled with a Resolution M050 excimer laser workstation, located at Laboratorio de Estudios Isotópicos, Centro de Geociencias, Universidad Nacional Autónoma de México (UNAM). A “squid” signal homogenizer was employed right after the ablation cell, before the ablated material entered the plasma. 350 ml/min of He were used as a carrier gas together with ca. 800 ml/min of Ar, mixed downstream with 4.5 ml/min of N2. A frequency of 5 Hz was utilized during this work, with a constant on-target fluence of 5 J/cm2, monitored at the beginning and end of each analytical session with an external energy monitor and systematically employing, throughout the entire duration of this study, an analytical spot of 60 μm. Background was acquired for 15 s, whereas the laser fired for 30 s, for a total sweep time of 0.530 ms, in which all isotopes of interest were acquired, simultaneously, for chemistry and geochronology (see Tables S1 and S2 in the supplementary material). We did not measure the depth of holes in rutile. However, comparative studies performed on zircons (Bernal et al., 2014) with the same instrument allow us to estimate the penetration depth for rutile in approximately 12-15 μm during the 30 s of analysis. Two rutile standards were employed for this work: R-10 as primary standard, whereas R-19 acted as secondary (control) standard. The accepted ages for the R-10 and R-19 standards are 1090±5 Ma and 493±5 Ma, respectively (Luvizotto et al., 2009; Zack et al., 2011).

To take into account the modeled down-hole fractionation observed in the primary rutile standard, the data reduction was performed using the Iolite 3.5 package (Paton et al., 2010; 2011), making use of the UComPbine data reduction scheme of Chew et al. (2014). The calculated age uncertainties corresponded to a two-standard error. Concordia diagrams for all the analyzed samples are plotted using IsoplotR (Vermeesch, 2018). Since rutile incorporates, at the moment of its crystallization, a variable amount of common Pb, this must be carefully evaluated to perform a proper correction of the isotope ratios. Fortunately, rutile does not incorporate virtually any Th in its crystalline network (e.g., Table S2). The common Pb can be thus calculated monitoring 208Pb, since 208Pb is formed by decay from 232Th but, in the absence of the radioactive isotope, all the contained 208Pb can be considered “common”. Therefore, the correction can be properly applied to the unknown samples (e.g., Figure 3). In general, we analyzed 50-60 rutile crystals from each sample. Finally, by monitoring the signal of particular major and trace elements, we were able to detect any anomaly, such as the presence and ablation of mineral inclusions that could affect the geochronological and chemical data. In such cases, the signal integration in Iolite was manually adjusted or the analysis was rejected. A further filtering was done on discordance, eliminating those analyses which were outside the -5 to 30 % range of concordance (uncolored ellipses in Figure 4).

Figure 3 Example of an unknown rutile sample, with data uncorrected (dashed ellipses) and corrected (green ellipses) for common Pb using the methodology explained in the text. Ellipses represent 2-sigma uncertainties. 

Figure 4 U-Pb concordia diagrams for rutile crystals belonging to the studied OC samples. Ellipses represent 2-sigma uncertainties. Insets: Kernel density estimation (KDE) of rutile ages (blue colored) compared with zircon age distributions (pink colored), obtained by Solari et al. (2014) for the same samples. KDE according to Vermeesch (2012). Valid number of analyses (n) for both rutile and zircon in each sample are indicated by pink and blue number, respectively. The beige and orange vertical stripes indicate the timing of the Zapotecan and Olmecan orogeny in the OC, respectively (cf.Solari et al., 2003; 2014). 

The secondary standard, measured over the course of these analytical sessions, yielded a concordia age of 495.9±5.9 Ma, in agreement with the aforementioned accepted age (Figure 5). The NIST-SRM 610 glass standard was also analyzed in each experiment (at the beginning, at the end and every 30 unknown analyses) to allow recalculation of trace and REE element concentrations in rutile crystals, employing 49Ti as internal standard. The reproducibility of NIST standard glasses by LA-ICPMS is very good, with relative standard deviations in general better than 1 % (e.g., Solari et al., 2010, 2015).

Figure 5 U-Pb concordia diagram for data of the secondary standard R-19 measured over the course of this study. Ellipses are at 2-sigma uncertainty. 

RESULTS

Rutile U-Pb Geochronology

Once corrected for common Pb, most of the rutile analyses yield concordant or slightly discordant results (Figure 4). All the samples have a unimodal distribution (Figure 4, insets), with apparent 206Pb/238U ages straddling the concordia in the 779-1031 Ma range (sample OC1101, Figure 4a), in the 917-1114 Ma range (sample OC1006, Figure 4b) or in the 911-1167 Ma range (sample OC1014, Figure 4c). The sample OC 1016 yields rutile analyses clustering in the 850-1009 Ma range, although few grains have younger ages, down to 779 Ma (Figure 4d). Finally, the sample OC1017 shows a main age cluster at 910-1026 Ma (Figure 4e).

Rutile trace-element composition and Zr-in thermometry

Common trace elements were measured in the studied rutile crystals on the same analyzed spot. Detrital rutile chemistry in younger, siliciclastic rocks, reflects its original formation environment. The characterization of the chemical composition in rutile crystallizing in basement rocks, such as the case of the OC, can thus draw important data to be used reference in the study of younger sediments (e.g., Zack et al., 2004; Triebold et al., 2012; Rösel, 2019; Martini et al., 2020). In general, the studied OC rutile grains fall in the metapelitic field of the Cr vs. Nb diagram proposed by Zack et al. (2004), characterized by high Nb and Cr values (compare OC samples in Figure 6a, with three examples of rutile belonging to mafic rocks from Mexico and Guatemala). The concentration of most of the chemical elements measured in OC rutile is variable (examples in Figure 6b to 6e, all analyzed elements in Table S2). The elemental distribution is, in several cases, similar from one sample to the other, with some elements enriched, especially in the OC-1016 sample (quartzofeldsphatic in composition), such as Mo, Zr, Hf, for instance. U is more or less equally distributed, with median concentrations ranging from 23 to 63 ppm (Figure 6e). Finally, it is important to underline the Zr-in thermometer behavior of the studied rutile crystals (according to the formulas of Tomkins et al., 2007): they define median crystallization temperatures in the 808-873 °C range, with the sample OC-1016 recording the highest values among the studied samples (Figure 6f).

Figure 6 Geochemical composition of analyzed rutile crystals from metapelitic samples of the Oaxacan Complex. a) Cr vs. Nb classification diagram of Zack et al. (2004). Rutile composition of three metabasic samples are included for comparison: ME13-21 is a Permian garnet-bearing amphibolite belonging to the Caltepec fault zone, Mexico (e.g.Elías-Herrera et al., 2007); Alex is a rutile sample belonging to late Precambrian-early Paleozoic anorthosite-gabbro sequences in Chiapas, Mexico (Cisneros de León et al., 2017); Gt-0327 is a Cretaceous mafic eclogite from the Sierra de Chuacús, Guatemala (Ortega-Gutiérrez et al., 2004). b) to e) box plots with the median value (number close to each box) of some of the elements analyzed in rutile for each OC sample. Concentrations are in ppm. f) Rutile temperature, in Celsius degrees, calculated using the Tomkins et al. (2007) Zr-in thermometer. The number close to each box corresponds to the calculated median temperature. 

DISCUSSION

Rutile closure temperature and Pb diffusion

As briefly stated in the introduction, rutile Pb closure temperatures are calculated according to different parameters and techniques. Recent studies have modified the initial findings of Mezger et al. (1989), who calculated a Pb closure temperature of ca. 400 °C, to higher values of ca. 600-630 °C (Cherniak, 2000; Kooijman et al., 2010 and 2012; Smye and Stockli, 2014). The diffusion studies performed by the former authors make also evident that Zr-in-rutile thermometry allows the calculation of near-metamorphic peak rutile crystallization temperatures, especially in granulite facies rocks (e.g., Kooijman et al., 2012). It is also clear that Pb diffusion in rutile is completely decoupled from Zr concentrations and distribution, making it impossible to correlate Zr-in-rutile temperatures with U-Pb ages. In this study we did not perform X-Y or depth profile analyses (cf.Smye and Stockli, 2014), so the age or concentration calculations are an average of the material sampled in the analyzed spot (60 μm in diameter and around 12-14 μm in depth, see also Bernal et al., 2014). Because the size of the studied rutile crystals (see Sample Description and Figure 2) and the granulite facies metamorphism that affected the OC rocks, whose temperature exceeded the Pb closure temperature reported above, we assume that most of the crystals were reset during the Zapotecan orogeny (ca. 990 Ma, Solari et al., 2003) and, thus, the calculated ages reflect the cooling through the Pb blocking temperature of ca. 600 °C. Although we do not discard the existence of eventual inherited rutile grains (for instance, some of the older rutile in sample OC1006, Figure 4b), for the present work in which we discuss their cooling after granulite metamorphism, we consider the maximum kernel density estimation (KDE hereafter) rutile distribution (insets in Figure 4) and the difference in time and temperature with the granulite metamorphic peak. The concordia lower intercept could give a more reliable age of the cooling event (see, for instance, Mezger and Krogstad, 1997). However, the age distribution of some of the studied samples (e.g., sample OC1006 and OC1016, Figures 4b and 4d, respectively) does not allow the calculation of a lower intercept in the concordia diagram. It is important to underline, however, that for those samples in which the lower intercept can be calculated (remainder samples in Figure 4), ages obtained are within error of the KDE maximum distribution. We are thus confident that the KDE maximum distribution is a reliable parameter that can be considered as the best approximation to the post-metamorphic peak cooling age for the studied samples.

Oaxacan Complex cooling curves

Combining the analytical results obtained in this research on rutile grains belonging to the OC, some interesting interpretations can be drawn. Chromium and Ni concentrations of the analyzed rutile crystals suggest a sedimentary protolith for the analyzed granulite rocks. This is important because it is sometimes complex to understand the protolith origin of granulite rocks (see, for instance, Solari et al., 2014, 2020).

The rutile KDE age distributions (Figure 4) allow a comparison with the detrital zircon grains of the same samples that were previously analyzed by Solari et al. (2014). Zircon, which is a mineral with a U-Pb closure temperature at about >900 °C, shows during magma crystallization or eventually, as is the case of OC, during the granulite-facies metamorphism, a variable age distribution in KDE diagrams that reflect either the protolith (age distributions older than ca. 1.0 Ga) or metamorphic crystallization age (age distribution of ca. 990 Ma). Arc-forming events (ca. 1.15-1.2 Ga, Keppie et al., 2001 and 2003; Weber et al., 2010; Weber and Schulze, 2014) are followed by two tectonothermal pulses, namely the Olmecan and Zapotecan Orogeny at ca. 1100 and ca. 990 Ma, respectively (Solari et al., 2003; Solari et al., 2014). While the Olmecan is mostly related to a yet-to-be understood, local migmatization event, the Zapotecan orogeny is the granulite-facies metamorphic event synchronous throughout Oaxaquia (see, for instance, Lawlor et al., 1999; Weber and Köhler, 1999; Solari et al., 2003 and 2014; Cameron et al., 2004; Weber and Schulze, 2014). While zircons yield multimodal distributions, witnessing the magmatic pulses in which they crystallized, or eventually, the granulite-facies metamorphism that allowed zircon recrystallization, rutile show unimodal distributions, reflecting their systematic U-Pb reset during granulite facies metamorphism at ca. 990 Ma. The age difference of the rutile main distribution from the Zapotecan event reflects the timing of cooling from 808-873 °C (rutile crystallization temperature according to the Zr-in thermometer, Figure 6F) to ca. 600 °C for the U-Pb closure temperature for rutile. The calculated cooling rates from the metamorphic peak temperature and age indicate, from N to S (being sample OC1101 the northernmost and sample OC1017 the southernmost ones, Figure 1), a range from ca 3 °C/Ma in the northernmost OC outcrops, to ca. 6 °C/Ma in the south. The sample OC1006 collected just south of Nochixtlán (Oax.) shows a maximum distribution at 985 Ma, that almost overlaps with analyzed zircon grains in the same sample (cf. Solari et al., 2014), indicative of a very fast cooling rate calculated with ca. 42 °C/Ma. Similarly, sample OC1014 shows a cooling rate of ca. 45 °C/Ma. Sample OC1016 presents an intermediate cooling rate of about 11 °C/Ma and two of the analyzed crystals show an even slower cooling rate of 1.5 °C/Ma (with ages of ca. 800 Ma, Figure 4D). It is possible that these grains were preserved as inclusions in garnet, where they are often observed in thin section, thus allowing such slow cooling rate.

Few published data are available to decipher the cooling history of the granulite basement in Mexico and, among those, rutile has not been employed yet. In the first study integrating cooling ages, Keppie et al. (2004) interpreted the whole OC as having cooled through one unique curve with an intermediate cooling rate of ca. 8 °C/Ma since the granulite metamorphism occurred at ca. 945 Ma, and then through a slower cooling rate of ca. 2 °C/Ma to expose the OC to the surface at about 710-760 Ma. Those events were associated to a tectonic switching from steep to flat slab subduction, and to thermal relaxation of the lithosphere, respectively. Our data are indicative of a more complex history, in which rutile plays an important role in identifying a heterogeneous cooling in the OC, which is indicative of the existence of several tectonic blocks, or slices, which behaved independently during the exhumation, right after the Zapotecan event. The shear zones previously characterized by Solari et al. (2003) in the northern OC are probably only some of the tectonic boundaries between those slices, mostly reactivated during the complex tectonics that affected the Mexican basement during the Phanerozoic (e.g., Solari et al., 2004a) and, based on the current rutile cooling ages, possibly acted already after the Zapotecan metamorphic peak, during the initial cooling and exhumation. Similar conclusions were recently drawn by Abdullin et al. (in press), whose data indicate that the heterogeneity characterized the OC cooling and exhumation even during the Mesozoic.

CONCLUSIONS

The application of rutile thermometry and U-Pb geochronology on crystals belonging to some metamorphic units of the Grenvillian Oaxacan Complex, southern Mexico, demonstrates that the exhumation after the Zapotecan granulite facies event (ca. 990 Ma) was heterogeneous among the different tectonic slices that compose the OC.

While this is a first approximation to U-Pb analysis in rutile, which further deserves a deeper study to better characterize the Pb diffusion during high-grade metamorphism, this example illustrates how combining the dating of different minerals with different closure temperatures (in this case, rutile and zircon, but the example holds true also when including apatite and micas, among others), collected in well-defined tectonic slices, is fundamental to construct meaningful cooling curves through geologic time in basement rocks.

Moreover, the chemical distribution of some of the analyzed rutile trace elements allow a further characterization, which can be useful in provenance analyses of younger sedimentary sequences in which rutile can be found as detrital heavy mineral (e.g., Ortega-Flores et al., 2018; Martini et al., 2020) due to its sensitivity to discriminate between protoliths.

SUPPLEMENTARY MATERIAL

Supporting supplementary Tables S1 y S2 can be found at the journal web site <http://rmcg.unam.mx/>, in the html version of this paper.

ACKNOWLEDGMENTS

This work is a direct contribution to the PAPIIT-DGAPA UNAM grant IN103417 to LAS. The first author thanks Centro de Geociencias, UNAM, for a scholarship received during his stay for the analytical work carried on for his Bachelor Thesis at Universidad Autónoma of San Luis Potosí, México, that finally resulted in the present work. Two anonymous journal reviewers are also thanked for their deep criticisms, which helped to better express some of the concepts reported here.

REFERENCES

Abdullin, F., Solé, J., Meneses-Rocha, J. J., Solari, L.A., Shchepetilnikova, V., Ortega-Obregón, C., 2015, LA-ICP-MS-based apatite fission track dating of the Todos Santos Formation sandstones from the Sierra de Chiapas (SE Mexico) and its tectonic significance: International Geology Review, 58(1), 32-48, doi:10.1080/00206814.2015.1055596 [ Links ]

Abdullin, F., Solari, L.A., Solé, J., Ortega-Obregón, C., 2021, Mesozoic exhumation history of the Grenvillian Oaxacan Complex, southern Mexico: Terra Nova. [ Links ]

Alaniz-Álvarez, S.A., Nieto-Samaniego, A.F., Ortega-Gutiérrez, F., 1994, Structural evolution of the Sierra Juárez Mylonitic Complex, state of Oaxaca, Mexico: Revista Mexicana de Ciencias Geológicas, 11(2), 147-156. [ Links ]

Alaniz-Álvarez, S.A., van der Heyden, P., Nieto-Samaniego, A.F., Ortega-Gutiérrez, F., 1996, Radiometric and kinematic evidence for Middle Jurassic strike-slip faulting in Southern Mexico related to the opening of the Gulf of Mexico: Geology, 24(5), 443-446. [ Links ]

Bernal, J.P., Solari, L.A., Gómez-Tuena, A., Ortega-Obregón, C., Mori, L., Vega-González, M., Espinosa-Arbeláez, D.G., 2014, In-situ 230Th/U dating of Quaternary zircons using LA-MCICPMS: Quaternary Geochronology, 23, 46-55, doi: 10.1016/j.quageo.2014.06.003 [ Links ]

Cameron, K., Lopez, R., Ortega-Gutiérrez, F., Solari, L., Keppie, J., Schulze, C., 2004, U-Pb geochronology and Pb isotopic compositions of leached feldspars: Constraints on the origin and evolution of Grenville rocks from eastern and southern Mexico: Geological Society of America Memoir, 197, 755-769. [ Links ]

Cherniak, D.J., 2000, Pb diffusion in rutile: Contributions to Mineralogy and Petrology, 139(2), 198-207. http://doi.org/10.1007/PL00007671 [ Links ]

Cherniak, D.J., Watson, E.B., 2000, Pb diffusion in zircon: Chemical Geology, 172(1), 5-24. [ Links ]

Chew, D.M., Petrus, J.A., Kamber, B.S., 2014, U-Pb LA-ICPMS dating using accessory mineral standards with variable common Pb: Chemical Geology, 363, 185-199, doi: 10.1016/j.chemgeo.2013.11.006 [ Links ]

Cisneros de León, A., Weber, B., Ortega-Gutiérrez, F., González-Guzmán, R., Maldonado, R., Solari, L., Schaaf, P., Manjárrez-Juárez, R., 2017, Grenvillian massif-type anorthosite suite in Chiapas, Mexico: Magmatic to polymetamorphic evolution of anorthosites and their Ti-Fe ores: Precambrian Research, 295, 203-226, doi:10.1016/j.precamres.2017.04.028 [ Links ]

Cochrane, R., Spikings, R. A., Chew, D., Wotzlaw, J.-F., Chiaradia, M., Tyrrell, S., Schaltegger, U, van der Lelij, R., 2014, High temperature (>350°C) thermochronology and mechanisms of Pb loss in apatite: Geochimica Et Cosmochimica Acta, 127(C), 39-56, doi:10.1016/j.gca.2013.11.028 [ Links ]

Dávalos-Álvarez, O.G., Nieto-Samaniego, A., Alaniz-Alvarez, S., Martínez-Hernández, E., Ramírez-Arriaga, E., 2007, Estratigrafía cenozoica de la región de Tehuacán y su relación con el sector norte de la falla de Oaxaca: Revista Mexicana de Ciencias Geológicas, 24(2), 197-215. [ Links ]

Elías-Herrera, M., Ortega-Gutiérrez, F., 2002, Caltepec fault zone: An Early Permian dextral transpressional boundary between the Proterozoic Oaxacan and Paleozoic Acatlán complexes, southern Mexico, and regional tectonic implications: Tectonics, 21(3) 4-1-4-18. [ Links ]

Elías-Herrera, M., Ortega-Gutiérrez, F., Sánchez-Zavala, J.L., Macías-Romo, C., Ortega-Rivera, A., Iriondo, A., 2007, The Caltepec fault zone: Exposed roots of a long-lived tectonic boundary between two continental terranes of southern México: Geological Society of America Special Papers, 422, 317-342, doi:10.1130/2007.2422(11) [ Links ]

Garzanti, E., 2016, From static to dynamic provenance analysis-Sedimentary petrology upgraded: Sedimentary Geology, 336, 3-13, doi:10.1016/j.sedgeo.2015.07.010 [ Links ]

González-Guzmán, R., Weber, B., Manjarrez-Juárez, R., Cisneros de León, A., Hecht, L., Herguera-García, J.C., 2016, Provenance, age constraints and metamorphism of Ediacaran metasedimentary rocks from the El Triunfo Complex (SE Chiapas, Mexico): evidence for Rodinia breakup and Iapetus active margin: International Geology Review 58, 2065-2091. [ Links ]

Keppie, J., Dostal, J., 2007, Rift-related basalts in the 1.2-1.3 Ga granulites of the northern Oaxacan Complex, southern Mexico: evidence for a rifted arc on the northwestern margin of Amazonia: Proceedings of the Geologists’ Association, 118(1), 63-74. [ Links ]

Keppie, J.D., Dostal, J., Ortega-Gutiérrez, F., Lopez, R., 2001, A Grenvillian arc on the margin of Amazonia: evidence from the southern Oaxacan Complex, southern Mexico: Precambrian Research 112, 165-181. [ Links ]

Keppie, J., Dostal, J., Cameron, K., Solari, L., 2003, Geochronology and geochemistry of Grenvillian igneous suites in the northern Oaxacan Complex, southern Mexico: tectonic implications: Precambrian Research, 120, 365-389. [ Links ]

Keppie, J.D., Solari, L.A., Ortega-Gutiérrez, F., Ortega-Rivera, A., Lee, J.K.W., Lopez, R., Hames, W.E., 2004, U-Pb and 40Ar/39Ar constraints on the cooling history of the northern Oaxacan Complex, southern Mexico: Tectonic implications, in Tollo, R.P., Corriveau, L., McLelland, J., Bartholomew, M.J. (eds.), Proterozoic tectonic evolution of the Grenville orogen in North America: Boulder, Colorado, Geological Society of America Memoir 197, 771-781. [ Links ]

Kooijman, E., Mezger, K., Berndt, J., 2010, Constraints on the U-Pb systematics of metamorphic rutile from in situ LA-ICP-MS analysis: Earth and Planetary Science Letters, 293, 321-330, doi: 10.1016/j.epsl.2010.02.047. [ Links ]

Kooijman, E., Smit, M.A., Mezger, K., Berndt, J., 2012, Trace element systematics in granulite facies rutile: implications for Zr geothermometry and provenance studies: Journal of Metamorphic Geology, 30, 397-412, doi: 10.1111/j.1525-1314.2012.00972.x [ Links ]

Lawlor, P., Ortega-Gutiérrez, F., Cameron, K., Ochoa-Camarillo, H., Lopez, R., Sampson, D., 1999, U-Pb geochronology, geochemistry, and provenance of the Grenvillian Huiznopala Gneiss of Eastern Mexico: Precambrian Research, v. 94, 73-99. [ Links ]

Luvizotto, G., Zack, T., 2008, Nb and Zr behavior in rutile during high-grade metamorphism and retrogression: An example from the Ivrea-Verbano Zone: Chemical Geology, 261, 303-317. [ Links ]

Luvizotto, G., Zack, T., Meyer, H., Ludwig, T., Triebold, S., Kronz, A., Münker, C., Stockli, D.F., Prowatke, S., Klemme, S., Jacob, D.E., von Eynatten, H., 2009 Rutile crystals as potential trace element and isotope mineral standards for microanalysis: Chemical Geology, 261, 346-369. [ Links ]

Martini, M., Solari, L., Peña-Guerrero, M., Zepeda-Martínez, M, Montomoli, C., 2020, Guidelines for assessing the provenance of Mesozoic and Cenozoic clastic successions sourced by pre-Jurassic basement complexes in Mexico: Journal of Sedimentary Research, 90, 513-532. [ Links ]

Mezger, K., Krogstad, E., 1997, Interpretation of discordant U-Pb zircon ages: an evaluation: Journal of Metamorphic Geology, 15, 127-140. [ Links ]

Mezger, K., Hanson, G., Bohlen, S., 1989, High-precision U-Pb ages of metamorphic rutile: application to the cooling history of high-grade terranes: Earth and Planetary Science Letters, 96, 106-118. [ Links ]

Mora, C., Valley, J., 1985, Ternary feldspar thermometry in granulites from the Oaxacan Complex, Mexico: Contributions to Mineralogy and Petrology, 89, 215-225. [ Links ]

Morton, A., Hallsworth, C., 1999, Processes controlling the composition of heavy mineral assemblages in sandstones: Sedimentary Geology, 124(1), 3-29. [ Links ]

Ortega-Flores, B., Solari, L., Lawton, T. F., Ortega-Obregón, C., 2013, Detrital-zircon record of major Middle Triassic-Early Cretaceous provenance shift, central Mexico: demise of Gondwanan continental fluvial systems and onset of back-arc volcanism and sedimentation: International Geology Review, 56(2), 237-261, doi: 10.1080/00206814.2013.844313 [ Links ]

Ortega-Flores, B., Martini, M., Guerrero-Moreno, S., Colás, V., Solari, L., Centeno-García, E., Silva-Romo, G., Grajales-Nishimura, M., 2018, Provenance analysis of Oligocene sandstone from the Cerro Pelón area, southern Gulf of Mexico: International Geology Review, 61(8), 915-925, doi: 10.1080/00206814.2018.1476922 [ Links ]

Ortega-Gutierrez, F., Ruiz, J., Centeno-Garcia, E., 1995, Oaxaquia, a Proterozoic microcontinent accreted to North America during the late Paleozoic: Geology, 23(12), 1127-1130. [ Links ]

Ortega-Gutiérrez, F., Solari, L., Solé, J., Martens, U., Gómez-Tuena, A., Morán-Icál, S., Reyes-Salas, M., Ortega-Obregón, C., 2004, Polyphase, High-Temperature Eclogite-Facies Metamorphism in the Chuacús Complex, Central Guatemala: Petrology, Geochronology, and tectonic Implications: International Geology Review, 46, 445-470. [ Links ]

Ortega-Gutiérrez, F., Elías-Herrera, M., Morán Zenteno, D., Solari, L., Weber, B., Luna-González, L., 2018, The pre-Mesozoic metamorphic basement of Mexico, 1.5 billion years of crustal evolution: Earth Science Reviews, 183, 2-37, doi:10.1016/j.earscirev.2018.03.006 [ Links ]

O’Sullivan, G.J., Chew, D.M., Morton, A.C., Mark, C., Henrichs, I. A., 2018, An Integrated Apatite Geochronology and Geochemistry Tool for Sedimentary Provenance Analysis: Geochemistry Geophysics Geosystems,19, 1309-1326, doi: 10.1002/2017GC007343 [ Links ]

Park, H., Barbeau, D. L., Jr, Rickenbaker, A., Bachmann Krug, D., Gehrels, G., 2010, Application of Foreland Basin Detrital‐Zircon Geochronology to the Reconstruction of the Southern and Central Appalachian Orogen: The Journal of Geology, 118(1), 23-44, doi:10.1086/648400 [ Links ]

Paton, C., Woodhead, J.D., Hellstrom, J.C., Hergt, J.M., Greig, A., Maas, R., 2010, Improved laser ablation U-Pb zircon geochronology through robust downhole fractionation correction: Geochemistry Geophysics Geosystems, 11, Q0AA06, doi: 10.1029/2009GC002618 [ Links ]

Paton, C., Hellstrom, J., Paul, B., Woodhead, J., Hergt, J., 2011, Iolite: Freeware for the visualisation and processing of mass spectrometric data: Journal of Analytical Atomic Spectrometry, 26, 2508, doi: 10.1039/c1ja10172b [ Links ]

Rösel, D., 2019, Interpretation and significance of combined trace element and U-Pb isotopic data of detrital rutile: a case study from late Ordovician sedimentary rocks of Saxo-Thuringia, Germany: International Journal of Earth Sciences, 108(1), 1-25, doi: 10.1007/s00531-018-1643-5 [ Links ]

Schulze, C., 2011, Petrología y geoquímica de las rocas del área de Pluma Hidalgo, Oaxaca e implicaciones tectónicas para el Proterozoico de Oaxaquia: México, Universidad Autónoma Nacional de México, PhD thesis, 311 pp. [ Links ]

Schulze, C., Keppie, J., Ortega-Rivera, A., Ortega-Gutiérrez, F., Lee, J., 2004, Mid-Tertiary cooling ages in the Precambrian Oaxacan Complex of southern Mexico: indication of exhumation and inland arc migration: Revista Mexicana de Ciencias Geológicas, 21(2), 203-211. [ Links ]

Shchepetilnikova, V., Solé, J., Solari, L., 2015, A chronological and chemical zircon study of some pegmatite dikes and lenses from the central part (Ayoquezco-Ejutla) of the Oaxacan Complex, southern Mexico: Revista Mexicana de Ciencias Geológicas, 32(1), 123-143. [ Links ]

Smye, A.J., Stockli, D.F., 2014, Rutile U-Pb age depth profiling: A continuous record of lithospheric thermal evolution: Earth and Planetary Science Letters, 408, 171-182, doi: 10.1016/j.epsl.2014.10.013 [ Links ]

Solari, L., Keppie, J.D., Ortega-Gutiérrez, F., Cameron, K., Lopez, R., Hames, W., 2003, 990 and 1100 Ma Grenvillian tectonothermal events in the northern Oaxacan Complex, southern Mexico: roots of an orogen: Tectonophysics, 365(1-4), 257-282, doi: 10.1016/S0040-1951(03)00025-8 [ Links ]

Solari, L., Keppie, J., Ortega-Gutiérrez, F., Ortega, A., Hames, W., Lee, J., 2004a, Phanerozoic structures in the Grenvillian Northern Oaxacan Complex, Southern Mexico: Result of thick-skinned tectonics: International Geology Review, 46, 614-628. [ Links ]

Solari, L., Keppie, J., Ortega-Gutiérrez, F., Cameron, K. L., López, R., 2004b, 990 Ma peak granulitic metamorphism and amalgamation of Oaxaquia, Mexico: U-Pb zircon geochronological and common Pb isotopic data: Revista Mexicana de Ciencias Geológicas, 21(2), 212-225. [ Links ]

Solari, L.A., Torres de León, R., Hernandez Pineda, G., Solé, J., Solís-Pichardo, G., Hernández-Treviño, T., 2007, Tectonic significance of Cretaceous Tertiary magmatic and structural evolution of the northern margin of the Xolapa Complex, Tierra Colorada area, southern Mexico: Bulletin of the Geological Society of America, 119, 1265-1279. [ Links ]

Solari, L., Gómez-Tuena, A., Bernal, J., Pérez-Arvizu, O., Tanner, M., 2010, U-Pb Zircon Geochronology with an Integrated LA-ICP-MS Microanalytical Workstation: Achievements in Precision and Accuracy: Geostandards and Geoanalytical Research, 34, 5-18. [ Links ]

Solari, L. A., Ortega-Gutiérrez, F., Elías-Herrera, M., Ortega-Obregón, C., Macías-Romo, C., Reyes-Salas, M., 2014, Detrital provenance of the Grenvillian Oaxacan Complex, southern Mexico: a zircon perspective: International Journal of Earth Sciences, http://doi.org/10.1007/s00531-013-0938-9 [ Links ]

Solari, L.A., Ortega-Obregón, C., Bernal, J.P., 2015, U-Pb zircon geochronology by LAICPMS combined with thermal annealing: Achievements in precision and accuracy on dating standard and unknown samples: Chemical Geology, 414, 109-123, doi: 10.1016/j.chemgeo.2015.09.008 [ Links ]

Solari, L.A., Ortega-Obregón, C., Ortega-Gutiérrez, F., Elías-Herrera, M., 2020, Origin and evolution of the Grenvillian Oaxacan Complex, southern Mexico: Hf isotopic and U-Pb geochronological constraints, in Martens, U., Molina-Garza, R.S. (eds.), Southern and central Mexico: basement framework, tectonic evolution, and provenance of Mesozoic-Cenozoic basins: GSA Special Paper 546, doi: 10.1130/SPE546 [ Links ]

Spencer, C.J., Prave, A.R., Cawood, P.A.,Roberts, N.M.W., 2014, Detrital zircon geochronology of the Grenville/Llano foreland and basal Sauk Sequence in west Texas, USA: Geological Society of America, Bulletin, 126(7-8), 1117-1128, doi: 10.1130/B30884.1 [ Links ]

Tomkins, H.S., Powell, R., Ellis, D.J., 2007. The pressure dependence of the zirconium-in-rutile thermometer: Journal of Metamorphic Geology, 25, 703-713. [ Links ]

Tolson, G., 2005, La falla Chacalapa en el sur de Oaxaca: Boletín de la Sociedad Geológica Mexicana, Volumen Conmemorativo del Centenario, Grandes Fronteras Tectónicas de México, Tomo LVII, 1, 111-122. [ Links ]

Triebold, S., Eynatten, von, H., Zack, T. 2012, A recipe for the use of rutile in sedimentary provenance analysis: Sedimentary Geology, 282, 268-275, doi: 10.1016/j.sedgeo.2012.09.008 [ Links ]

Vermeesch, P., 2012, On the visualization of detrital age distributions: Chemical Geology, 312-313, 190-194. [ Links ]

Vermeesch, P., 2018, IsoplotR: a free and open toolbox for geochronology: Geoscience Frontiers, 9, 1479-1493 [ Links ]

Watson, E.B., Wark, D.A., Thomas, J.B., 2006, Crystallization thermometers for zircon and rutile: Contributions to Mineralogy and Petrology, 151, 413-433. [ Links ]

Weber, B., Kohler, H., 1999, Sm-Nd, Rb-Sr and U-Pb geochronology of a Grenville Terrane in Southern Mexico: origin and geologic history of the Guichicovi Complex: Precambrian Research, 96, 245-262. [ Links ]

Weber, B., Scherer, E.E., Schulze, C., Valencia, V.A., Montecinos, P., Mezger, K., Ruiz, J., 2010, U-Pb and Lu-Hf isotope systematics of lower crust from central-southern Mexico - Geodynamic significance of Oaxaquia in a Rodinia Realm: Precambrian Research, 182, 149-162, doi: 10.1016/j.precamres.2010.07.007 [ Links ]

Weber, B., Schulze, C. 2014, Early Mesoproterozoic (>1.4 Ga) ages from granulite basement inliers of SE Mexico and their implications on the Oaxaquia concept - Evidence from U-Pb and Lu-Hf isotopes on zircon: Revista Mexicana de Ciencias Geológicas, 31(3), 377-394. [ Links ]

Zack, T., von Eynatten, H., Kronz, A., 2004, Rutile geochemistry and its potential use in quantitative provenance studies: Sedimentary Geology, 171, 37-58. [ Links ]

Zack, T., Stockli, D. F., Luvizotto, G. L., Barth, M. G., Belousova, E., Wolfe, M. R., Hinton, R. W., 2011, In situ U-Pb rutile dating by LA-ICP-MS: 208Pb correction and prospects for geological applications: Contributions to Mineralogy and Petrology, 162(3), 515-530, doi: 10.1007/s00410-011-0609-4 [ Links ]

Adame-Martínez, M.G., Solari, L.A., Ortega-Obregón, C., Abdullin, F., 2020, U-Pb geochronology of rutile: deciphering the cooling history of the Oaxacan Complex granulites, southern Mexico: Revista Mexicana de Ciencias Geológicas, v. 37, núm. 2, p. 135-145.

Received: November 26, 2019; Revised: June 02, 2020; Accepted: June 04, 2020

* E-mail: solari@unam.mx

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License