First Lu-Hf, δ18O and trace elements in zircon signatures from the Statherian Espinhaço anorogenic province (Eastern Brazil): geotectonic implications of a silicic large igneous province

Manuscript ID: 20180046. Received on: 04/17/2018. Approved on: 08/22/2018. ABSTRACT: The Espinhaço rift system encompasses taphrogenetic events from the Statherian to Tonian in the São Francisco-Congo (SFC) paleocontinent. The magmatism is represented mainly by metamorphosed anorogenic granites and rhyolites with subordinate amphibolites. Zircon U–Pb (LA-ICPMS and SHRIMP) ages from felsic (1748 ± 3 Ma and 1740 ± 8 Ma) and mafic (1725 ± 4 Ma) rock samples, coupled with previous studies suggest that the Espinhaço igneous province erupted from ca. 1.79 Ga to ca. 1.70 Ga. The felsic rocks show characteristics of A-type magmas. The negative εHf(t) data for meta-rhyolite zircons (-12.32 to -17.58), the moderate δ18O values (7.02 to 7.98) and the REE patterns suggest crustal melting related to an extensional environment. The mafic rock shows negative values of εHf(t) in zircons (-4.05 to -8.25) and moderate δ18O values (5.56 to 7.87). The results disclose a basaltic magmatism in continental intraplate setting whose parental magma could have been derived from the subcontinental lithospheric mantle with contamination of crustal material. These data coupled with coeval Espinhaço magmatism and mafic dyke swarms found to the south of the Espinhaço rift system reinforce the evidence of a long-lived Statherian silicic large igneous province (SLIP) on the SFC paleocontinental block.


INTRODUCTION
The outstanding development of chemical and isotopic analysis on zircon grains enhanced the understanding of nature, evolution and source of igneous rocks related to distinct tectonic environments and its geodynamic processes. Because zircon grains are able to record multiple geological events, through its morpho-chemical modifications, numerous investigations, as U-Pb and Lu-Hf systematics, as well as δ 18 O and REE abundances, have been focused on relating temporal evolution of the crust and the lithospheric mantle (Hoskin andSchaltegger 2003, Harley andKelly 2007). Even some representative zircon samples can disclose new findings and test previous hypothesis suggested by lithochemical and whole-rock isotopic methods. Such analytical methods of zircon spot analysis are here presented on samples of felsic and mafic rocks from the Espinhaço igneous province which rift-related anorogenic nature has been suggested since the early 1990's, based only on lithochemical data from felsic rocks (Chemale Jr. 1987, Dussin 1994, Fernandes 2001).
The São Francisco-Congo (SFC) paleocontinental block (Fig. 1A), amalgamated in the Rhyacian -Orosirian boundary (ca. 2.05 Ga), experienced a series of rifting events from the Statherian to Cryogenian (Pedrosa-Soares and Alkmim 2011, Chemale Jr. et al. 2012a, Danderfer et al. 2015, Heilbron et al. 2017. The oldest among those taphrogenic events took place in Early Statherian time, around 1.75-1.70 Ga, cutting across the São Francisco block in Central-Eastern Brazil, and culminating with the opening of the extensive N-Strending Espinhaço basin system associated with abundant acid magmatism (Fig. 1B).
To the east of the main exposure region of the Southern Espinhaço rift, Statherian felsic magmatism is also registered in the Guanhães Archean basement block (Fig. 1B), where large A-type granitic intrusions and rhyolitic volcanism have been referred to for decades (e.g., Dussin 1994, Fernandes 2001. These meta-rhyolites are found in close association with meta-mafic rocks (ortho-amphibolites), both lacking detailed analytical studies. Therefore, they are now the first bimodal suite related to the Espinhaço rift characterized with robust analytical data.
Besides new U-Pb (LA-MC-ICP-MS and SHRIMP) zircon ages, lithochemical and mineral chemistry data, we present the first Lu-Hf, δ 18 O and trace elements in zircon data for the mafic and felsic magmatism related to the Espinhaço rift system. Our robust dataset, together with a thorough data compilation from the literature, demonstrate the rift-related bimodal nature of the Espinhaço magmatic province. Furthermore, it supports the role of subcontinental lithospheric mantle on the petrogenesis of these anorogenic rock suite. Finally, we discuss paleotectonic correlations concerning the São Francisco paleocontinental block, envisaging the Statherian Espinhaço magmatic event as part of a silicic large igneous province (SLIP).  Alkmim (2004) and Pinto and Silva (2014).

GEOLOGICAL SETTING
The Espinhaço basin system extends for about 1.200 km along the N-S direction, cutting across the São Francisco paleocontinental block, and can be subdivided into five domains, namely: the Southern Espinhaço ridge, Guanhães block, Central Espinhaço ridge, Northern Espinhaço ridge, and the Chapada Diamantina (i.e., Diamantina plateau) domains (Fig. 1B). While the Northern Espinhaço ridge and Chapada Diamantina domains were mostly preserved from orogenic processes within the São Francisco craton, the Southern and Central Espinhaço ridge domains and Guanhães block were involved in the Brasiliano orogeny within the Neoproterozoic Araçuaí orogen, where they were being mostly metamorphosed in the greenschist facies but locally reaching upper amphibolite facies metamorphism , Cruz and Alkmim 2017. The Espinhaço Supergroup is mostly composed of siliciclastic sequences and magmatic rocks, with minor carbonate rocks, associated with three taphrogenic events that occurred during the Statherian (ca. 1.8-1.68 Ga), Calymmian-Ectasian (ca. 1.6-1.38 Ma) and Stenian (ca. 1.18 Ma) (Pedrosa-Soares and Alkmim 2011, Chemale Jr. et al. 2012a, Guadagnin et al. 2015, andreferences therein).
Located to the east of the Southern Espinhaço ridge domain (Fig. 1B), the Guanhães block stands out as a portion of the São Francisco paleocontinental region reworked within the Araçuaí Neoproterozoic orogen (Alkmim et al. 2006, Noce et al. 2007a, 2016. The Guanhães block encompasses Archean TTG migmatitic gneisses and associated granites, locally covered by Siderian to Neoproterozoic supracrustal successions, as well as a voluminous Statherian granitic intrusions and meta-ultramafic and meta-mafic rocks of unknown age (Noce et al. 2007a, Teixeira et al. 2017. Zircon U-Pb (SHRIMP) ages between 3150 Ma and 2710 Ma constrain the magmatic crystallization of both gneiss protoliths and granitoids that also show local metamorphic overprints around 2000 Ma and 570-500 Ma , 2016, Peixoto et al. 2015, Barrote 2016.
The Statherian magmatism, aim of this paper, mostly includes anorogenic meta-granites and meta-rhyolites, and minor meta-mafic rocks (Brito Neves et al. 1979, Dussin 1994. The meta-granitic rocks display a wide age spectrum, around 1790-1700 Ma and are represented by the Borrachudos, Catolé and Lagoa Real suites, located in the Guanhães block, Central and Northern Espinhaço domains, respectively (Turpin et al. 1988, Pimentel et al. 1994, Dussin et al. 2000, Fernandes 2001, Costa 2013. Meta-rhyolites dated around 1750 Ma, 1735 Ma and 1710 Ma suggest distinct pulses of felsic volcanism along the Espinhaço system (Machado et al. 1989, Babinski et al. 1994, Schobbenhaus et al. 1994, Danderfer et al. 2009, 2015. Statherian metamafic rocks dated between 1750 Ma and 1700 Ma occur in the Southern and Central Espinhaço ridge domains (Dussin 1994, Chemale et al. Jr. 2012a, Bezerra-Neto 2016, Silva 2016, Moreira 2017. Typical sedimentary rift deposits filled the Statherian Espinhaço basin system. They include sandstones, rudites and pelites related to alluvial fan, braided fluvial and lacustrine environments, with maximum depositional ages between 1.80 Ga and 1.68 Ga (Martins-Neto 2000, Chemale Jr. et al. 2012a, Santos et al. 2013, Danderfer et al. 2015. In the Guanhães block, where our detailed study area is located, the N-S-trending Alto Rio Guanhães unit comprises a metavolcano-sedimentary rock assemblage tectonically imbricated with the Archean basement, the Guanhães Group and the Statherian meta-granites (Fig. 2). Firstly recognized by Danderfer and Meireles (1987), the Alto Rio Guanhães unit records metamorphosed rocks from the greenschist to the amphibolite facies. The eastern portion of this unit is mostly composed by mafic schists and ortho-amphibolites, representing basic volcanic rocks and associated intrusive mafic bodies, with minor meta-rhyolites and metasedimentary rocks. The western portion mostly comprises quartzites, ferruginous quartzites, pelitic schists, calcsilicate rocks and iron formations, associated with meta-mafic rocks and meta-rhyolites (Fig. 2). Some meta-ultramafic bodies, mostly composed of tremolite-talc schists, are also found within the Alto Rio Guanhães unit.

RESULTS
Samples were collected from representative outcrops of the Alto Rio Guanhães unit and Borrachudos suite. The samples are free of weathering and hydrothermal alteration. We selected samples of mafic (sample J02: 18º38'13''S, 43º12'50''W) and rhyolitic (sample J03: 18º36'15''S, 43º12'52''W) rocks for LA-ICPMS and SIMS isotopic analysis on zircon grains. One sample of a typical syenogranite of the Borrachudos suite, collected from the Açucena pluton (LC40: 18º49'19''S, 42º09'14''W), was selected for zircon U-Pb SHRIMP analyses. Appendix A presents the descriptions of the analytical methods applied in the following sections.

Sample petrography and mineral chemistry
The samples selected for detailed studies are metamorphosed mafic and felsic rocks from the Alto Rio Guanhães unit, and meta-granites of the Borrachudos suite. The mafic rocks are fine-to coarse-grained ortho-amphibolites, representing volcanic/subvolcanic and plutonic rocks of basaltic composition. The felsic samples include volcanic to subvolcanic rocks of rhyolitic composition and related granites. Mineral chemistry data tables are found in the Supplementary Table A1.

Mafic rocks
The ortho-amphibolites are fine-to coarse-grained rocks, essentially composed of hornblende and plagioclase ( Fig. 3A and 3C). The main accessory minerals are quartz, titanite, epidote, apatite and ilmenite. Amphibole is largely dominant (up to 70% in modal content), forming granonematoblastic, stretched, rotated and/or sigmoidal grains along the regional foliation (Fig. 3B). It also exhibits a mottled texture outlined by fine-grained quartz inclusions with interlobate boundaries, suggesting statically recrystallization of exsolved silica. Decussated overgrown amphiboles imprint the anastomosed regional foliation and obliterate previous textures. Amphibole commonly shows zoned crystals, varying from  Mg-hornblende to tschermakite, with higher Mg/(Mg+Fe 2+ ) ratios and silica contents in cores than borders, which in turn are richer in Al VI , Al IV , Fe, and Ti (Fig. 4A).
The plagioclase grains (up to 40% in volume) form granoblastic aggregates. Except one sample with low-Ca andesine (An 31-44 ) -they are non-zoned grains ranging from labradorite (An 62 ) to anorthite (An 90 ) (Fig. 4B). Quartz, accompanied or not by epidote, mostly forms fine-grained xenoblastic inclusions within plagioclase and amphibole, as exsolved silica excess. The chemical variations in amphibole (from Mg-hornblende to more Al-Fe rich types) and plagioclase (with increasing Ca contents) are consistent with growth under prograde metamorphism.

Felsic rocks
The meta-rhyolites are porphyroclastic rocks formed by a well-foliated fine-grained matrix of microcline, quartz, plagioclase and biotite, enveloping milimetric eye-shaped porphyroclasts of quartz and microcline ( Fig. 3E and 3F). K-feldspar is the major mineral phase, containing mica and drop-like quartz inclusions, and ranging in composition between Or 93 Ab 07 and Or 96 Ab 04 (Fig. 4B). Quartz also forms symplectite intergrows with K-feldspar. Plagioclase crystals are oligoclase (An 14-20 ) in composition (Fig. 4B). They may show polysynthetic twinning, inclusions of quartz, mica and apatite, as well as saussurite. Brown biotite grains occur oriented according to the regional foliation and belong to the siderophyllite-annite series, showing a great variation in the Al IV contents, as well as high values and small variation in Fe/(Fe+Mg) ratios varying from 0.859 to 0.876. Accessory phases are essentially zircon and apatite.
The meta-granitoids are foliated to isotropic, mediumto coarse-grained, biotite-and/or hornblende-bearing syenogranite to monzogranite (Fig. 3G). Feldspars may show drop quartz inclusions, and vermicular intergrowth with quartz. Saussuritization of plagioclase is a common process. Brown biotite occurs oriented according to the regional foliation or forming spotted aggregates. Accessory minerals include apatite, zircon, titanite, magnetite and, occasionally, garnet.

Mafic rocks
Results from major and trace elements analysis for 11 basaltic samples from the Alto Rio Guanhães suite are given in Supplementary Table A2. The concentrations of most major and trace elements, as well as high-field-strength elements (HFSE; e.g., Nb, Zr, Y) and rare earth elements (REE), seem to represent the primary magmatic contents of the study samples. All samples have low LOI values (0.2 to 1.3 wt%).
The ortho-amphibolite is silica-oversaturated with normative quartz, hypersthene, albite and anorthite, which suggests that even after the crystallization of mafic minerals and feldspar there was still silica crystallizing as quartz. These samples are tholeiitic subalkaline basalts, with silica and alkali contents ranging from 45.1 to 49.3 wt% and 1.1 to 2.5 wt%, respectively. MgO values show a narrow variation ranging from 6.1 to 8.6 wt% with Mg# between 43.7 and 54.9.

Zircon U-Pb isotopic ages
An ortho-amphibolite (sample J02) and a meta-rhyolite (sample J03), both interlayered with metasedimentary rocks in the Alto Rio Guanhães unit, were selected for LA-ICPMS U-Pb, Lu-Hf and δ 18 O isotope and REE in zircon analysis. Representative CL images of zircons from these volcanic rocks are shown in Fig. 7, along with the locations of spots measured for the in situ isotopic analyses.

Ortho-amphibolite
Zircons of the ortho-amphibolite are subhedral and irregularly shaped grains with average lengths of 100-200 μm and length to width ratios around 2:1. The oscillatory concentric zoning is typical of magmatic origin (Fig. 7). They have moderate U content with variable Th/U (0.06-1.24) (Suppl. Data A3). Thirty one spot data (from a total of 63) provide a robust Concordia age of 1725.1 ± 3.9 Ma (MSWD = 0.025; Fig. 8A), which is interpreted to represent the crystallization age of the igneous protolith.

Meta-granite
Zircons from a hornblende-biotite syenogranite (Açucena pluton) of the Borrachudos suite were dated using the SHRIMP technique. The grains are homogeneous with medium luminescence in CL image and show oscillatory zoning. The analytical results of 12 spots in 10 crystals are shown in Supplementary  Table A3 and Fig. 8C. Eight spots belonging to the same population of crystals (MSWD = 0.69) yield a concordant grouping of 1740.5 ± 7.8 Ma, interpreted as the age of magmatic crystallization. Similar age is obtained by the 12 spots that define an upper intercept age of 1739.0 ± 8.8 Ma (MSWD = 0.73).

Hf isotope composition of zircons
Ten Hf isotopic spots analyses were conducted on zircons from the ortho-amphibolite (J02) and 12 analyses on those from the meta-rhyolite (J03). Analytical data are given in the Supplementary Table A4 and presented in Fig. 9A.
The zircons from the ortho-amphibolite show relatively uniform Hf isotopic compositions with radiogenic 176 Hf/ 177 Hf ratios of 0.281453 to 0.281571 for ages between 1717 and 1732 Ma. All analyzed zircon grains yielded consistently negative εHf (t), varying from -8.25 to -4.05, which suggests contamination during ascent through the crust. 176 Hf/ 177 Hf ratios for the analyzed spots of the meta-rhyolite zircons range from 0.28118 to 0.28132 with correspondent εHf (t), varying from -17.58 to -12.32 for ages between 1735 Ma and 1760 Ma, which suggests that this magma was essentially sourced from long-lived crustal rocks. The Hf (T DM ) model ages for the  (1990), Dussin (1994), Fernandes et al. (1994), Knauer and Grossi-Sad (1997), Oliveira (2002) and Silveira-Braga (2012).
ortho-amphibolite zircon grains range between 2.28 Ga and 2.44 Ga, while those for the meta-rhyolite vary from 2.62 Ga to 2.80 Ga, indicating the involvement of Archean and possibly Siderian to Orosirian rocks in magma genesis.

Zircon O isotopes
The oxygen isotopic compositions of zircon grains from the ortho-amphibolite and meta-rhyolite samples are listed in Supplementary Table A4. There are no consistent differences in δ 18 O between the rims and the cores of the analyzed grains in both samples.
The meta-rhyolite zircons have δ 18 O values varying from 6.8 to 8.0, with 53% of the analyzed zircons (n = 27) showing values between 7.5 and 8.0 (Fig. 9B). For the ortho-amphibolite zircons, δ 18 O values are variable from 5.6 to 7.9 (Fig. 9B), with 52% of the analyses between 7.0 and 7.5 (n = 16). Two grains have values around 5.6, which are similar to those of zircon crystallizing from uncontaminated mantle-derived magma (5.3 ± 0.3‰, Valley et al. 1998, Valley 2003. The meta-rhyolite zircons show a δ 18 O range typical of felsic rocks while the ortho-amphibolite zircon grains are consistent with contamination of mantle-derived magma that interacted with continental crust material (Valley et al. 2005, Bindeman 2008).

Trace element chemistry of zircons
Zircon grains from the ortho-amphibolite and meta-rhyolite have been analyzed for trace elements concentrations by LA-ICPMS. The measurements were performed in the zircon cores, in the same textural domains as the U-Pb and Lu-Hf analyses (Suppl. Data A4). The grains of both samples show very similar REE patterns, with HREE enrichment and LREE depletion relative to chondrite values ( Fig. 10A and 10B). The ortho-amphibolite zircons (Sm/La) N range from 48 to 431, (Lu/Gd) N values ranging between 11 to 24 and Th/U ratio from 0.3 to 0.8. The meta-rhyolite zircon grains display (Lu/Gd) N values between 11 to 19, and the Th/U ratio varies from 0.7 to 1. The REE patterns show positive Ce anomaly and negative Eu anomaly, features consistent with unaltered igneous zircon (Hoskin andSchaltegger 2003, Rubatto 2017). The negative Eu anomaly is typical of plagioclase fractionation while the positive Ce anomaly is related to oxidation processes or to an oxidizing environment (Hoskin and Schaltegger 2003).
The REE patterns are similar to populations of crustal zircon (Hoskin and Schaltegger 2003, Grimes et al. 2007, Harley and Kelly 2007 and overlap the REE concentration of zircons from oceanic crust and continental granitoids ( Fig. 10A and   10B). This correspondence between the REE concentration of zircons from different types of rocks can be explained by the REE compatibility in the zircon lattice (e.g., Grimes et al. 2007).
The zircon crystals from both samples have intermediate U/Yb ratio. The Hf and Y vs. U/Yb discriminant diagrams attest the continental nature for the meta-rhyolite zircon grains and indicate crustal input for the ortho-amphibolite zircons ( Fig. 10C and 10D).

DISCUSSION
Here we discuss the petrogenesis of the mafic and felsic magmatism related to the Espinhaço rift system, the correlations in view of the Statherian scenario in the São Francisco block, as well as the global inferences in relation to paleocontinental settings.

Magmatic crystallization
The mafic rocks from the Alto Rio Guanhães unit have compositions suggesting that fractional crystallization played an important role in the petrogenesis of the basaltic magmatism. These rocks are characterized by positive correlation of CaO, Na 2 O, Al 2 O 3 , Ni and Cr, in relation to Mg#, suggesting fractionation of olivine, pyroxene and plagioclase. The negative Sr anomaly of most samples corroborates plagioclase fractionation. The range between the ratios of the incompatible trace elements (Zr, Y, Nb) of the distinctly evolved samples is less than 12%. This and the absence of compositional gaps and the Pearson correlation coefficient values suggest that fractional crystallization process was involved in the generation of these rocks.
Considering that the ratios of incompatible elements that have similar geochemical behavior tend to remain constant during fractional crystallization, the plot in the Figure 11A shows the decreasing values of Ni accompanied by constancy of Zr/Nb ratio.

Assimilation and fractional crystallization (AFC)
In general, fractional crystallization is accompanied by assimilation of crustal material, mainly in rift-related environments. The negative values of εHf (t) in zircons between -4.05 to -8.25 for the ortho-amphibolite sample indicate a moderately juvenile magma significantly contaminated by continental crust material. The δ 18 O data registered for these zircon grains characterizes shifted magmatic δ 18 O values by crustal assimilation (item 3.3.3). Crustal input is also suggested by the U/Yb ratio of zircons (item 3.3.4).
Thus, in order to evaluate the influence of crustal contamination during the mafic magma ascent and/or residence in magma chamber, the DePaolo (1981) equation was used to verify the involvement of AFC processes. The results presented in Fig. 11B show that the values of Zr and Nb incompatible trace elements of the mafic samples, among others, can be explained by up to 30% of AFC process in the upper crust.

Magma mantle source
On the basis of the above discussion (items 4.1.1.1 and 4.1.1.2), fractional crystallization and crustal assimilation The AFC curve (red squares) shows ranges of 10% (between 0 and 90%). The arrow indicates that the tendency of the magmatic evolution of the analyzed mafic rocks (green diamond) is the same of the AFC curve, indicating values up to 30% AFC. The continental crust chemical data is from Taylor and McLennan (1985), and the mineral/melt partition coefficients for basaltic liquids are after the compilation of Rollinson (1993). (C) (La/Sm) N vs. (Tb/Yb) diagram to determine the source of the mafic magmas. The horizontal dashed line separates the fields for melting of garnet-bearing peridotite and of spinel-bearing peridotite Wang et al. (2002). Normalization values of the chondrite C1 after Sun and McDonough (1989). (D) Zr vs. Zr/Y discrimination plot (Pearce and Norry 1979). were important processes in the petrogenesis of the mafic magmatism in the Alto Guanhães unit. The similar REE and trace element patterns for theses rocks (Fig. 5A and 5B) suggest a common mantle source for them.
The amount of melting required to generate tholeiitic basaltic melts from the peridotitic mantle ranges from 10 to 20% (Jaques and Green 1980). Owing to the distinctive REE partition coefficients between garnet and spinel during melting of spinel peridotites vs. garnet peridotites, a spinel-bearing peridotite source is suggested for these rocks (Fig. 11C). The Alto Rio Guanhães mafic suite is characterized by low (Tb/Yb) N ratios (1.4-1.8) and (Gd/Yb) N ratios (≤ 2), typical of rocks generated by partial melting of a mantle source in the spinel-stable field and, thus, at depths less than 60-70 km (Wilson 2007). The magmatism shows variation of degrees of melting or source enrichment, according to the (La/Sm) N range (1.9-2.7). The limited range of (Tb/Yb) N ratios also suggests small variation in the depth of melting, relating to small variation in the crustal thickness during the progressive continental rifting.
La/Yb and La/Nb ratios can be used as indicative for discriminating fertile, enriched and depleted mantle sources. The meta-mafic rocks show (La/Yb) N > 1 (4.0-6.5) and (La/Nb) N > 1 (1.2-2.2) that are characteristics of an enriched source associated with the lithospheric mantle (Sun 1980, Humphris et al. 1985, Sun and McDonough 1989, as expected in continental rift systems. Despite the range of the major and trace elements, the meta-mafic samples from the Alto Rio Guanhães unit generally exhibit incompatible trace elements patterns more enriched than E-MORB, except for the values of Rb, Ba and Sr, which resemble N-MORB (Fig. 12). However, such signatures may be expected for evolved continental rift systems.

Tectonic setting
The mafic rocks of the Alto Rio Guanhães unit have never been dated in detail, as well as all other rocks of the Espinhaço system. The studied mafic rocks show patterns of incompatible trace elements comparable to those of intraplate basalts, such as high Ti and Zr, instead of basalts of volcanic arc and meso-oceanic ridge. Regarding the Pearce and Norry (1979) (Fig. 11D) and Wood (1980) diagrams, widely used for basaltic rock discrimination, all samples plot in the intraplate basaltic field, suggesting that the Alto Rio Guanhães mafic magmatism occurred in an attenuated continental lithosphere. Furthermore, the age of this magmatism fits very well with the age range of the felsic rocks found throughout the Espinhaço system. Thereby, this mafic magmatism can be related to the same continental taphrogenic process that developed during the Statherian in the São Francisco paleocontinental region.

Chemical proxies for anorogenic magmatism
The meta-rhyolites and intrusive granites in the Guanhães block show negative correlations between SiO 2 and TiO 2 , FeO (t) , CaO and Al 2 O 3 , suggesting fractionation of hornblende, biotite and Fe-Ti oxides during magma evolution. CaO and Al 2 O 3 contents coupled with compatible behavior of Sr, Eu and Ba with SiO 2 , point to fractionation of plagioclase. In the plots Sr vs. Rb (Fig. 13A), the variations in these contents seem to be mainly related to fractionation of plagioclase and K-feldspar.
The compatible/incompatible element ratios, such as Rb/Sr and K/Rb, are used to recognize the degree of fractionation of granitic magmas (e.g., Blevin 2003). In the SiO 2 vs.

A B
K/Rb and SiO 2 vs. Rb/Sr diagrams ( Fig. 13B and 13C), the granitic rock compositions in the study area plot in the moderately and strongly evolved fields, showing characteristics of differentiated magmas. This suggests an evolution for the granitic compositions by fractional crystallization processes. These felsic rocks have high SiO 2 and alkali contents, and are generally enriched in HFSE and REE. The Ga/Al*1000 ratios range from 2.8 to 7.5, which is typical of A-type granites, according to Whalen et al. (1987). These felsic magmatism also fits in the definition of A-type granitic rocks by Loiselle and Wones (1979) since they present, in addition to SiO 2 and alkalis, high levels of Nb, Zr, REE, negative Eu anomaly and low values of CaO and MgO. According to the geochemical subdivision of A-type granites proposed by Eby (1992), this felsic magmatism belongs to sub-type A2 (Fig. 13D).

Tectonic setting and magmatic evolution
A-type felsic magmas are usually generated in extension environments, where crustal stretching, mantle uplift and the heat transferred from the mantle can reach the shallowest levels of the crust (Loiselle and Wones 1979). However, the origin of A-type magmatism is quite controversial and the mechanics available in the literature describes both melting process of continental crust by underplating of melts derived from the mantle, and an evolution by fractional crystallization of basaltic magmas concomitant with crustal assimilation -AFC (e.g., Whalen et al. 1987, Eby 1992, Rollinson 1993, Peccerillo et al. 2003.
The studied meta-rhyolites and intrusive granites do not present evidence of being products of evolutionary processes involving mafic magmatism. The volume of felsic rocks is much larger than the mafic rocks and can not be explained by fractional crystallization processes of basaltic magmas. The surface area covered only by the Borrachudos granites is bigger than 6.000 km 2 . In addition, the felsic rocks show high ratios between incompatible trace elements, as Th/La, Zr/Ti and Th/Nb, in contrast to the values obtained for the mafic magmatism, which are up to 15 times smaller.
The felsic magmatism is characterized by highly negative values of εHf (t) (-12.32 to -17.58), indicating association with crustal material. The same is revealed by the positive δ 18 O values (7.02 to 7.98) and REE pattern of these zircon grains. The Hf T DM ages argue against a simple magmatic differentiation model, suggesting the generation of felsic magmas from an ancient continental crust. In the tectonic classification diagrams from Pearce et al. (1984), the samples plot in the field of intraplate granites ( Fig. 13E and 13F).

Statherian geochronological framework of rift-related magmatism of the São Francisco block
The Supplementary Table A5 summarizes the U-Pb and Pb-Pb ages of magmatic events related to the Statherian Espinhaço rift system recorded in the Guanhães block, Espinhaço Ridge and Chapada Diamantina, and also in correlated areas in the São Francisco block. These magmatic events are related to taphrogenetic processes that marked the São Francisco-Congo paleocontinent during the Late Paleoproterozoic (e.g., Dussin 1994, Danderfer et al. 2015, Cederberg et al. 2016, Dussin 2017.
The tholeiitic meta-mafic rocks of the Alto Rio Guanhães unit are the first detailed studied record of a Statherian mafic suite coexisting with the well-known felsic magmatism related to the Espinhaço rift system. The only geochronological record so far for a mafic rock in the Guanhães block was presented by Dussin et al. (2000), which yielded 1697 ± 10 Ma 207 Pb/ 206 Pb single zircon evaporation data as the crystallization age of an amphibolite intercalation on metasedimentary rocks from the Guanhães Group. In the Southern Espinhaço Range (Fig. 1), highly metamorphosed and altered meta-volcanic mafic rocks associated with the opening of the Espinhaço basin show ages between 1747 Ma and 1700 Ma (Dussin 1994, Chemale Jr. et al. 2012a, Bezerra-Neto 2016, Silva 2016. The precise in situ U-Pb dating of a meta-mafic rock from the Guanhães block presented in this study (1725 ± 4 Ma) is somewhat higher than the previously reported age by Dussin et al. (2000), but is within the age range of the meta-volcanic mafic rocks of the Southern Espinhaço. More recently, a meta-volcanic mafic rock from the Central Espinhaço ridge domain yielded an age of 1730 ± 8 Ma (Moreira 2017), which is virtually the same age obtained for the ortho-amphibolite of the Alto Guanhães unit. The 1748 ± 3 Ma age for the meta-rhyolite also constraints a deposition age for the sedimentary infill of the Alto Rio Guanhães basin and allows us to correlate it with the whole Statherian Espinhaço rift system. The age of the meta-rhyolite (1748 ± 3 Ma) of the Alto Rio Guanhães unit is somewhat older than the Conceição do Mato Dentro meta-rhyolite (1715 ± 5 Ma; Machado et al. 1989) of the Southern Espinhaço domain, which is one of the nearest occurrences of such rocks in relation to the study area. Except this case, the Alto Rio Guanhães meta-rhyolite age is very similar to the U-Pb ages obtained for most felsic volcanic and plutonic rocks found along the Espinhaço system, like the Planalto de Minas meta-rhyolite (1752 ± 2 Ma; Machado et al. 1989) and a volcaniclastic rock (1758 ± 18; Costa 2017), both of the Central Espinhaço domain; the Rio dos Remédios meta-rhyolites of the Chapada Diamantina domain (1748 ± 4 Ma and 1752 ± 4 Ma; Babinski et al. 1994, Schobbenhaus et al. 1994; the Sapiranga meta-volcanic rock (1740 ± 10 Ma; Danderfer et al. 2015) of the Northern Espinhaço domain; the Lagoa Real plutonic suite found between the Chapada Diamantina and Northern Espinhaço domains (1744 ± 2 Ma, Pimentel et al. 1994); and the here presented U-Pb SHRIMP age for the Borrachudos suite (1740 ± 8 Ma).
The metasedimentary rocks of the Statherian Espinhaço basin, as quartzites and meta-conglomerates, show maximum depositional ages between 1.80 and 1.68 Ga and are well exposed in all Espinhaço rift branches (Danderfer et al. 2009, Chemale Jr. et al. 2012a, Santos et al. 2013). In the eastern edge of the Southern Espinhaço ridge, near the Guanhães block, Rolim et al. (2016) found a maximum depositional age of ca. 1666 Ma for a quartzitic unit. The available geochronological data for the quartzites and banded iron formations of the Guanhães Group suggest a maximum deposition average age of ca. 2080 Ma, with a single Statherian zircon age of 1737 ± 19 Ma ( Barrote 2016). For the southern region of the Guanhães block, Statherian depositional ages (ca. 1668 Ma) were obtained to a meta-siliciclastic succession with iron formation intercalations (Carvalho et al. 2014, Silveira-Braga et al. 2015.
The U-Pb dataset for the whole Espinhaço system shows a wide age variation between 1790-1700 Ma, but a concentration in the time interval of ca. 1750-1710 Ma. The lower part of the rift system is characterized by large volume of acid magmatism, whereas the mafic magmatism is only registered in the upper supracrustal successions of the Statherian Espinhaço system, as well as in the Guanhães block. Contrasting with the mafic magmatism so far restricted in space, both the plutonic and volcanic felsic rocks are found along the whole Espinhaço rift system.
Intracontinental magmatism of Statherian age is also recorded in the south and west regions of the São Francisco paleocontinental block. To the south, the Pará de Minas mafic dyke swarm (Chaves 2001, Cederberg et al. 2016, which shows two Statherian dyke generations (ca. 1795 Ma and 1710 Ma; Cederberg et al. 2016), can be correlated in space and time to the Espinhaço rift system. Along the São Francisco northwestern region, it is well-documented the Statherian (1771-1768 Ma) meta-rhyolites and related A-type plutonic rocks of the Araí rift system (Pimentel et al. 1991).

Paleotectonic constrains for the São Francisco-Congo paleocontinent
The São Francisco and Congo paleocontinental blocks (Fig. 1A) amalgamated after an orogenic event that took place in the Rhyacian -Orosirian boundary. Since then, they became part of a paleocontinent that was only broken up by the opening of the South Atlantic Ocean in the Cretaceous. From the Rhyacian -Orosirian to the Neoproterozoic Brasiliano -Pan-African orogenic events only taphrogenic records have been reported with solid data in the focused region of the São Francisco -Congo paleocontinental block (Noce et al. 2007b, Pedrosa-Soares et al. 2008, Pedrosa-Soares and Alkmim 2011, Heilbron et al. 2017.
Distinct paleotectonic and paleogeographical reconstructions have been envisaged for the São Francisco-Congo paleocontinental block in Paleoproterozoic time, since those focusing the agglutination of Atlantica (ca. 2 Ga; Rogers 1996) and Columbia (1.90-1.70 Ga;e.g., Zhao et al. 2004, Hou et al. 2008, Cederberg et al. 2016 to the Central African block (Cordani et al. 2013, D'Agrella-Filho and. Based on geotectonic, geochronological and paleomagnetic evidence, D'Agrella-Filho and  suggested that the São Francisco-Congo, Kalahari, Borborema, Trans-Sahara, Rio de la Plata and smaller paleocontinental blocks assembled in the Central African block around 2 Ga, a paleocontinent diachronic but independent to the Columbia supercontinent (Fig. 14A).
Alternatively, there are models picturing the São Francisco-Congo block within the Columbia supercontinent (e.g., Rogers and Santosh 2002, Hou et al. 2008, Teixeira et al. 2017. Hou et al. (2008) envisaged a continuous subduction-related magmatic belt bordering Columbia that was amalgamated by orogenic systems developed from ca. 2.1 Ga to ca. 1.8 Ga. (Fig. 14B). In this configuration, the São Francisco -Congo craton, located in the border of Columbia, is juxtaposed to the Amazonian craton, linking the South America to Columbia through the Trans-Amazonian orogen. After this model, Danderfer et al. (2015) related the evolution of the Espinhaço rift system to a far-field continental extension induced by orogenic processes located in the Amazonian domain.
Recent paleotectonic models quoting the Statherian anorogenic magmatism suggest a paleocontinental link between the São Francisco and North China cratons (Peng 2015, Cederberg et al. 2016, Teixeira et al. 2017Xu et al. 2017). These authors point to similarities between the Pará de Minas dyke swarm, located in the southern tip of the São Francisco craton, and the magmatism related to the Xiong'er-Taihang LIP event, both interpreted as related to Statherian intracontinental setting (Fig. 14C).
In this framework, the magmatism related to the Espinhaço rift system, including the Alto Rio Guanhães mafic suite, can be correlated to the Pará de Minas mafic dyke swarm and, probably, to the Xiong'er-Taihang magmatism. This suggestion is supported by geochemical (e.g., trace element patterns; Fig. 15) and geochronological data from the Pará de Minas mafic dyke swarm (Chaves 2001), the Xiong'er volcanic rocks (Zhao et al. 2002, He et al. 2009, the Taihang dykes (Hou et al. 2001, Peng 2015, Wang et al. 2004, 2007 and the tholeiitic amphibolites of the Alto Rio Guanhães mafic suite (this work).
The younger dyke generation of the Pará de Minas dyke swarm (1717( -1702Cederberg et al. 2016) can be chrono-correlated to the metamorphosed volcanic and plutonic mafic rocks of the Southern Espinhaço rift (1747-1700 Ma; e.g., Silva 2016), as also to the tholeiitic mafic rocks of the Alto Rio Guanhães suite (1725 Ma According to the concept suggested by Bryan and Ernst (2008) and Ernst (2014), a Silicic Large Igneous Province (SLIP) presents the following five characteristics: ■ the extrusive volume is greater than 0.25 Mkm 3 (up to > 3 Mkm 3 ) and the total exposure area larger than 0.1 Mkm 2 ; ■ more than 80% of the province is represented by daciterhyolite, with transitional calc-alkaline I-type to A-type intraplate signature; ■ rhyolitic ignimbrite is the predominant rock; ■ the duration of magmatism is up to 40 Ma and shows a pulsed nature related to shorter intervals of 3 Ma to 10 Ma; ■ the province is exclusively continental. The generation of such voluminous felsic magma is related to large degrees of partial melting of a hydrated, calc-alkaline, amphibolitic to andesitic lower continental crust (Bryan et al. 2002, Bryan and.
In fact, the concept of Silicic Large Igneous Province (SLIP) was applied to the anorogenic magmatism related to the Statherian Espinhaço rift system by previous authors (Danderfer et al. 2015, Chaves et al. 2016, and it is now checked in the light of our data and compilation. The present-day exposed crustal levels along the whole Espinhaço rift domains show rocks metamorphosed from the low greenschist to upper amphibolite facies, implying in the erosion of significant amounts of the Statherian volcanic rocks. Therefore, evaluations of the total volume and surface area once occupied by the overall Statherian anorogenic rocks are certainly underestimated. In this scenario, the Statherian Espinhaço magmatism is mostly characterized by metamorphosed, A-type, alkali-calcic to calc-alkalic granites, rhyolites, dacites, ignimbrites and acid volcaniclastic rocks, with continental intraplate signature (see references quoted in previous items). They represent more than 80% in exposed area of the whole Statherian igneous rocks and their metamorphosed equivalents. The mafic magmatism is very restricted, much probably representing less than 10% in area (Fig. 1).
Although spreading along a total area of more than 250.000 km 2 , including the Espinhaço ridge, Chapada A B Figure 15. Trace element patterns of the Alto Rio Guanhães mafic suite (this study), Pará de Minas dyke swarm (Chaves 2001), the Xiong'er volcanic rocks (Zhao et al. 2002, He et al. 2009, the Taihang dykes (Hou et al. 2001, Peng 2015, Wang et al. 2004, 2007: (A) normalized to E-MORB; (B) normalized to Primitive Mantle. The normalization values are from Sun and McDonough (1989).
Diamantina and Guanhães block domains (Fig. 1), the occurrence areas of both felsic and mafic Statherian igneous rocks reach at least 70.000 km 2 , because unknown parts of them are hidden beneath younger rock layers and extensive Cenozoic covers (Fig. 1). Although this surface area (> 70.000 km 2 ) is apparently smaller than the SLIP definition requirement (0.1 Mkm 2 or 100.000 km 2 ), it only shows the erosional remains found in present-day exposed crustal levels, most of them reworked within the Neoproterozoic Araçuaí orogen. Allowing to such exposure conditions, any evaluation of the original volume of extrusive rocks may be unreliable, although it can be expected a very large volume (possibly > 0.25 Mkm 3 ) of them, because the minimum area (> 70.000 km 2 ) and distribution of the present-day exposure areas of plutonic and volcanic felsic rocks are along more than 250.000 km 2 of relatively deep exposed crustal levels.
The high-quality U-Pb ages on zircons from the Statherian igneous rocks of the Espinhaço rift system span from 1792 ± 7 Ma to 1703 ± 12 Ma, with most values ranging from ca. 1752 Ma to ca. 1710 Ma (see Suppl. Data A5). This time interval (ca. 1752-1710 Ma) comprises at least 80% age values (almost all of them related to felsic rocks) of the dated Espinhaço igneous rocks, and is in good agreement with the maximum duration (ca. 40 Ma) suggested for a SLIP (cf. Ernst 2008, Ernst 2014).
Therefore, bearing in mind the crustal level, exposed rocks and tectonic setting, the Espinhaço anorogenic province is, indeed, a good candidate to represent a SLIP, because it meets the main definition requirements, as the exclusive continental setting; the striking predominance of felsic (> 80% of A-type granites and rhyolites) over mafic (< 10%) rocks; a maximum duration of the felsic magmatism up to 40 Ma; and the very extensive, minimum total exposure area (despite the relatively deep erosional levels).
The data for the Alto Rio Guanhães mafic suite suggests a mantle upwelling beneath a stretched and thinned continental lithosphere, pointing to the involvement of the enriched lithospheric mantle in the generation of these rocks (Fig. 16). Although the trace elements pattern of this magmatism does not support a derivation from fertile sublithospheric composition, considering the life span of the Espinhaço anorogenic province, the contribution of fertile melts from a thermal anomaly (or mantle plume) can not be discarded.
As mentioned above (item 4.3.1), beyond the spatial-and chrono-correlation, the Alto Rio Guanhães mafic rocks display similar chemical characteristics with the Pará de Minas dyke swarm (Fig. 15) in the southern of the São Francisco craton. Therefore, these magmatic records can be either derived from a common parental magma or from a similar mantle source, suggesting they are part of the same igneous province. most likely generated by melting of enriched subcontinental lithospheric mantle in the spinel-stable field and display the involvement of crustal assimilation accompanied by fractional crystallization processes in the generation of this magmatism. The anorogenic felsic magmatism may be derived by partial melting of Archean lower crustal materials by underplating of melts derived from the mantle; ■ The Espinhaço rift system is characterized by a large volume of acid magmatism during its opening, whereas the mafic magmatism is only registered in the upper sequences of the Statherian rift system. Contrasting with the mafic rocks, so far restricted to the Southern and Central Espinhaço ranges and also the Guanhães block, the felsic rocks characterize a long-lived magmatism during the Statherian, with a ca. 1750-1710 Ma major magmatic event recorded throughout the whole Espinhaço rift system; ■ Contemporaneous Statherian magmatism is reported elsewhere on the Southern São Francisco craton and, together with the Pará de Minas dyke swarm, the Alto Rio Guanhães rocks represent a rifting event;

CONCLUSION
■ The main volumetric, compositional and tectonic parameters of the Espinhaço magmatism provide the evidence of a SLIP in the São Francisco and Congo paleocontinental blocks.

ACKNOWLEDGEMENT
The authors acknowledge the research financial support provided by the CNPq and CODEMIG. Joana Magalhães is grateful to the Brazilian Geological Survey (CPRM). Special thanks are given to Alexey Ulyanov for the assistance in the acquisition and reduction of trace elements in zircon data and to Katharina Marger for the assistance with δ 18 O data acquisition.

SUPPLEMENTARY DATA
Supplementary data associated with this article can be found in the online version: Suplementary Table A1,  Suplementary Table A2, Suplementary Table A3, Suplementary  Table A4, Suplementary Table A5.