Open Access, thanks to the authors’ funding. Article available to all readers via GSW (http://ammin.geoscienceworld.org) and the MSA web site.

Abstract

This review presents a compositional database of primary anatectic granitoid magmas, entirely based on melt inclusions (MI) in high-grade metamorphic rocks. Although MI are well known to igneous petrologists and have been extensively studied in intrusive and extrusive rocks, MI in crustal rocks that have undergone anatexis (migmatites and granulites) are a novel subject of research. They are generally trapped along the heating path by peritectic phases produced by incongruent melting reactions. Primary MI in high-grade metamorphic rocks are small, commonly 5–10 μm in diameter, and their most common mineral host is peritectic garnet. In most cases inclusions have crystallized into a cryptocrystalline aggregate and contain a granitoid phase assemblage (nanogranitoid inclusions) with quartz, K-feldspar, plagioclase, and one or two mica depending on the particular circumstances. After their experimental remelting under high-confining pressure, nanogranitoid MI can be analyzed combining several techniques (EMP, LA-ICP-MS, NanoSIMS, Raman). The trapped melt is granitic and metaluminous to peraluminous, and sometimes granodioritic, tonalitic, and trondhjemitic in composition, in agreement with the different P-T-aH2O conditions of melting and protolith composition, and overlap the composition of experimental glasses produced at similar conditions. Being trapped along the up-temperature trajectory—as opposed to classic MI in igneous rocks formed during down-temperature magma crystallization—fundamental information provided by nanogranitoid MI is the pristine composition of the natural primary anatectic melt for the specific rock under investigation. So far ~600 nanogranitoid MI, coming from several occurrences from different geologic and geodynamic settings and ages, have been characterized. Although the compiled MI database should be expanded to other potential sources of crustal magmas, MI data collected so far can be already used as natural “starting-point” compositions to track the processes involved in formation and evolution of granitoid magmas.

Introduction

Granitoid rocks and their extrusive counterparts represent the largest products of magmatic activity on continental crust. In zones of crustal thickening associated with collisional orogens, there are several examples of crystalline basements, which have experienced repeated episodes of high-temperature metamorphism, partial melting, and extraction of granitoid melts (Sawyer et al. 2011; Brown 2013). On the other hand, in continental margin arc systems there is much field and geochemical evidence of melting of crustal rocks and hybridization between crustal- and mantle-derived magmas, resulting in the formation of batholiths (Barnes et al. 2002; Annen et al. 2006; Kemp et al. 2007; Gray and Kemp 2009). Therefore, anatexis of crustal protoliths and formation of anatectic granitoid magmas are key processes for the formation of new crust and its geochemical differentiation (Brown and Rushmer 2006; Hacker et al. 2011; Sawyer et al. 2011). Field geology, experimental petrology, thermodynamic calculations, numerical modeling, whole-rock geochemistry and Hf-O isotopes in zircon represent the major tools for understanding the origin and evolution of anatectic granitoid magmas (see Annen et al. 2006; Brown 2013; Clemens 2006; Kemp et al. 2007; Sawyer et al. 2011; White et al. 2011). Despite the vigorous debate on the petrogenesis of granitoid magmas that began more than 60 years ago (Read 1948) conflicting views still exist, in particular concerning the melting processes at the source and the differentiation processes that modify granitoid chemistry (Clemens and Watkins 2001; Clemens and Stevens 2012, 2015; Brown 2013; Weinberg and Hasalovà 2015a, 2015b and discussion therein).

Here we present a novel and alternative method to investigate crustal magmas, providing a review of the current research on melt inclusions (MI) in high-grade metamorphic rocks, the source where crustal magmas originate. After discussing the origin of MI by incongruent melting in migmatites, granulites, and anatectic enclaves and presenting the data collected mainly by the research group of the authors, along with the adopted experimental and analytical strategies, we highlight how such tiny inclusions of melt can constrain the origin and evolution of granitoid magmas produced by anatexis of crustal rocks.

Changing the Viewpoint: Entrapment of melt inclusions upon heating

MI are small volumes of silicate melt enclosed in minerals (Roedder 1984). First recognized in the 19th century by Sorby (1858), MI hosted in minerals of igneous rocks have been long studied and widely used in the last 40 years (see data set reported in Kesler et al. 2013), becoming a worldwide-accepted technique to investigate melt evolution in magmatic systems (Clocchiatti 1975; Roedder 1979; Sobolev 1996; Frezzotti 2001; Schiano 2003; Bodnar and Student 2006; Thomas and Davidson 2012; Audétat and Lowenstern 2014; and references therein). These MI (hereafter “classic MI”) are generally trapped during crystallization of magmas, i.e., along the cooling path where the host crystal is crystallizing from the melt that is being entrapped (Fig. 1). MI-bearing phenocrysts in lavas are a typical example of this mode of occurrence. Being trapped along the down-temperature trajectory (Fig. 1), such classic MI represent therefore snapshots on the liquid line of descent of silicate liquids (i.e., they represent variably evolved melts) and, particularly in the case of felsic systems, they have represented a key tool to investigate and understand the genesis of various ore deposits (Audétat and Lowenstern 2014 and references therein). Even in the case of entrapment by the first liquidus phase (for example olivine in some mantle-derived magmas), classic MI often approximate the composition of the parental melt (i.e., the most primitive melt found in an area) and not that of the primary one [i.e., the melt derived directly by partial melting of the source; see Fig. 1 in Bartoli et al. (2014)].

To represent primary melt compositions, MI should be entrapped in the source region, before magmas undergo segregation and differentiation processes such as mixing and/or mingling with other magma batches, double-diffusive convection, Soret diffusion, cumulus phenomena, fractional crystallization, volatile exsolution and assimilation of residual or exotic material (e.g., Burnham 1967; DePaolo 1981; Lesher et al. 1982; Barbarin 1988; Candela 1997; Weinberg et al. 2001; Stevens et al. 2007; Brown 2013). Indeed, all these processes produce final products—i.e., most of the extrusive and intrusive rocks we collect in the field—which are markedly different from the original melt produced at the source.

The continental crust (i.e., the source region of the anatectic crustal magmas) mostly melts through incongruent melting reactions (Weinberg and Hasalovà 2015a and references therein). Along with magma crystallization, such reactions fulfill the fundamental requirement to form MI, i.e., the growth of a mineral in the presence of a melt phase. Indeed during incongruent melting, both a silicate liquid and a solid phase—the peritectic mineral—are produced contemporaneously along the up-temperature path (Fig. 1). In such a scenario, the peritectic phase can trap droplets of the melt with which it is growing. Shielded by the peritectic host, these MI can provide us with the pristine geochemical signatures (“the starting point” as stated by Sawyer et al. 2011) of magmas formed by crustal anatexis. Although the possibility of trapping droplets of melt during incongruent melting may appear intuitive, MI in high-grade metamorphic rocks have received only little consideration in the past decades. After the first detailed microstructural and microchemical description of MI in regionally metamorphosed granulites of southern India (Cesare et al. 2009), many other occurrences have been reported worldwide in the last six years (e.g., Cesare et al. 2011, 2015; Ferrero et al. 2012, 2014, 2015; Groppo et al. 2012; Mosca et al. 2012; Bartoli et al. 2013a, 2014; Darling 2013; Kawakami et al. 2013; Barich et al. 2014; Carosi et al. 2015; Massonne 2014), demonstrating that MI in high-grade metamorphic rocks are more common than expected. More rarely, MI have also been described in enclaves of partially melted continental crust brought to the Earth’s surface in lavas. This is the case of the anatectic metasedimentary enclaves from the Neogene Volcanic Province (NVP) of southern Spain (Cesare et al. 1997—also known as “erupted migmatites,” Zeck 1970). From all the above considerations it is clear that petrologists can utilize a novel tool, i.e., MI trapped by peritectic minerals, to gain previously unavailable information and to better understand the origins and evolution of crustal magmas.

Microstructural characterization of melt inclusions

In most migmatites, granulites, and anatectic enclaves investigated so far, MI have been mostly found in garnets of the melt-depleted melanosome portions of the rock (Table 1), because this mineral is the most common peritectic phase associated with crustal melting reactions in a wide PTX range (see Cesare et al. 2011; Ferrero et al. 2012 and references therein). Other MI-bearing phases are plagioclase and zircon, and less commonly cordierite, andalusite, ilmenite, hercynite, and monazite (Table 1; Cesare et al. 2015). Crustal enclaves contain abundant MI both in reactants and peritectic phases, e.g., plagioclase and garnet (Cesare and Maineri 1999; Cesare et al. 2003; Acosta-Vigil et al. 2007, 2010). Indeed, during melting with some degree of reaction overstepping, some reactant phases in the rock matrix may recrystallize, equilibrate with the melt and trap MI (Cesare et al. 2015).

From a microstructural point of view, MI may occur isolated, distributed in zonal arrangements, or in clusters with random distribution within the host mineral (Fig. 2; Ferrero et al. 2012; Barich et al. 2014). These modes of occurrence unequivocally demonstrate the primary nature of MI (Roedder 1979). MI are generally small (<20 μm) and often show negative crystal shape (Figs. 3a–3d). Locally, larger MI have been described [up to 150–200 μm, Barich et al. (2014); Ferrero et al. (2015)]. Variable degrees of crystallization have been observed and MI range from glassy to fully crystallized (Figs. 2e and 3a–3e). The most common daughter minerals found within crystallized MI are quartz, alkali feldspar, plagioclase, biotite, and muscovite (Figs. 3a–3c). In MI from near ultrahigh-pressure conditions, metastable polymorphs of albite, K-feldspar, and quartz have also been identified (Fig. 3b; see also Ferrero et al. 2016). A large number of minerals trapped during inclusion formation has been found within MI, the most common being graphite, zircon, rutile, and kyanite. Owing to the very fine grain-size (often <1 μm) and the granitic phase assemblage observed in totally crystallized MI, they have been referred to as “nanogranites” (Cesare et al. 2009). Because recent microchemical investigations have demonstrated the existence of granodiorite, trondhjemite, and tonalite compositions (Carosi et al. 2015; this study), the totally crystallized polycrystalline MI will be referred, hereafter, to as “nanogranitoid inclusions.” Nanogranitoids can display a diffuse micro- to nano-porosity, which contains liquid H2O, as evidenced by micro-Raman mapping (Bartoli et al. 2013a). Nanogranitoids may also coexist with primary C-bearing fluid inclusions (Fig. 3f). As both inclusion types are primary, they indicate conditions of fluid-melt immiscibility, thus suggesting COH fluid-present anatexis (Cesare et al. 2007; Ferrero et al. 2011). Preserved glassy MI are abundant in anatectic enclaves owing to the fast cooling (quench) related to the magma uprise and extrusion (Fig. 3d). On the other hand, nanogranitoids and partially crystallized MI are typically found in slowly cooled, regionally metamorphosed migmatite and granulite terranes (Figs. 3a–3c). Surprisingly, some preserved glassy MI coexisting with nanogranitoids have been recognized also in migmatitic terranes (Fig. 3e; see also Cesare et al. 2009; Bartoli et al. 2013a, 2013b; Barich et al. 2014). This geological oddity (i.e., the glass preservation in slowly cooled rocks) has been ascribed to either pore size effect (Cesare et al. 2009), the heterogeneous distribution of nucleation sites (Ferrero et al. 2012), or the high viscosity of the trapped melt (Bartoli et al. 2015).

Chemical characterization of melt inclusions

Re-melting nanogranitoid inclusions

The experimental re-melting of nanogranitoid inclusions is a prerequisite to retrieve the original composition (in terms of major and volatile contents) of the trapped melt. Because anatexis takes place in rocks of the middle to lower continental crust, at the bottom of thickened continental crust, and even in continental materials subducted to mantle depths, re-melting experiments of nanogranitoids must be performed under high confining pressure, using a piston-cylinder apparatus, to prevent the inclusion from decrepitating, volatile loss, and melt contamination. Indeed, as a consequence of the H2O loss, crystallized inclusions can be completely re-melted only at temperatures higher than the trapping—i.e., H2O loss results in an increased solidus temperature—leading to an inevitable melt-host interaction (see details in Bartoli et al. 2013c). The first attempts to remelt nanogranites using the routine technique in igneous petrology, namely the high-temperature heating stage at ambient pressure (Esposito et al. 2012), produced the decrepitation of the nanogranitoid inclusions, making them poorly suited for a geochemical study (Cesare et al. 2009).

The P-T conditions to be used during piston-cylinder remelting should be initially established based on phase-equilibria modeling or classic thermobarometry, which can provide a fair estimate of the (presumed) P-T conditions of MI entrapment (Bartoli et al. 2013a, 2013c; Ferrero et al. 2015). Conditions of further experiments are then fine-tuned based on the results of the first runs. Despite the time consuming pre- and post-run sample preparation (compared to the heating stage method), and the trial-and-error nature of the experimental approach, this method can successfully re-melt the nanogranitoids under high-confining pressure (see Bartoli et al. 2013c; Ferrero et al. 2015).

Building up the compositional database

Working with MI in partially melted rocks (i.e., their identification, re-homogenization, and chemical characterization) is challenging due to their small diameter, often <20 μm (Figs. 2 and 3). Despite these difficulties, originally glassy and successfully re-melted MI can be analyzed for the concentrations of most elements in the periodic table. Major element contents are routinely measured by electron microprobe (EMP), analyzing secondary standards in the same analytical session to monitor the alkali migration process and eventually calculate correction factors (Ferrero et al. 2012). Indeed, the small size of MI often demands the use of a beam diameter of ~1 μm (Bartoli et al. 2013a), resulting inevitably in Na loss during EMP analysis of hydrous felsic glasses (Morgan and London 1996, 2005a). The laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) has been used to measure trace element contents of MI, drilling the glassy MI exposed on the host surface or the nanogranitoid inclusions entirely enclosed in the host. In the latter case, one of the major elements measured by EMP in coexisting glassy or re-melted MI can be used as internal standard for the deconvolution of the mixed signal (Halter et al. 2002). H2O, which is the most important volatile dissolved in crustal magmas (Burnham 1975), can be quantified by different methods: deficiency from 100% in the EMP analysis (Acosta-Vigil et al. 2007), Raman spectroscopy (Bartoli et al. 2013a; Ferrero et al. 2015), and secondary ion mass spectrometry (SIMS/NanoSIMS: Bartoli et al. 2013c, 2014). Using NanoSIMS and Raman spectroscopy, the H2O concentration of glass can be determined with a high spatial resolution of 1–2 μm. Bartoli et al. (2014) applied different methods to measure H2O in MI from amphibolite-facies migmatites and observed that, compared to NanoSIMS concentrations, the difference of EMP totals from 100% generally yields slightly higher-H2O contents (up to ~15% relative), whereas Raman spectroscopy underestimates H2O (up to approximately 20% relative) when compared to NanoSIMS values. Ferrero et al. (2015) reported lower differences between H2O values determined by Raman and EMP difference method (~9% relative on the average value).

By applying all these techniques, ~600 MI—including remelted nanogranitoids in migmatites and granulites and glassy MI in anatectic enclaves—have been characterized to date (Fig. 4), leading to the construction of an extensive geochemical database of MI in high-grade metamorphic rocks. In detail, the compositional database consists of ~600 and ~200 major and trace elements analyses, respectively, and ~140 measurements of H2O by either SIMS, NanoSIMS, or Raman. The interested reader may refer to the Supporting Information (Supplementary Tables1 1 and 2) and to Cesare et al. (2015) for tables where all the presented and discussed analyses are reported.

The analyzed MI come from three continents (Europe, Asia, and Africa) and are related to different geologic and geodynamic settings (Table 1) such as: (1) the late orogenic, low-to-medium P extensional setting represented by the Neogene Volcanic Province of southern Spain; (2) the low-P dynamothermal aureole of Ojén (Betic Cordillera, S Spain); (4) the low- and medium-P collision-related magmatism of the Maghrebian (La Galite, Tunisia) and Himalayan (Kali Gandaki, Nepal) mountain belts, respectively; (4) the low-to-medium P and ultrahigh-T regional metamorphism of the Kerala Kondalite Belt (India); and (5) the near ultrahigh-P metamorphism of the Bohemian Massif (Orlica-Śnieżnik Dome, central Europe). MI have been found in peraluminous metasedimentary rocks and metaluminous metaigneous rocks, and are inferred to have formed at conditions varying from 670 to 950 °C and 4 to 27 kbar, from Proterozoic to Quaternary/Pleistocene (Table 1). A complete list of published occurrences of MI in migmatites and granulites, along with a description of the host rock and conditions of melting, is reported in Cesare et al. (2015). For the sake of completeness, it should be noted that the research group of the writers has recently found MI also in Proterozoic and Archean rocks from the Americas, as well as, from Antarctica and Africa, but their microstructural and microchemical characterization is still in progress.

Reliability of the MI as records of melt compositions

During and after entrapment of melt its composition may be modified by several processes, making the MI unrepresentative of the bulk melt in the rock (Roedder 1984). Many processes such as the crystallization of the host on MI walls, the diffusional exchange of cations between the crystalline phases within nanogranitoids and the host mineral, the nucleation of bubbles by volume contraction and the retrograde fluid infiltration have been widely discussed by Acosta-Vigil et al. (2010, 2012), Bartoli et al. (2013c, 2014), and Cesare et al. (2015). Here, we consider two processes, with which igneous petrologists working with classic MI are more familiar: the boundary layer effect and diffusive H2O re-equilibration.

Boundary layers may develop at the crystal-melt boundary, when the rate of crystal growth is faster than the rate of cation diffusion (Bacon 1989). This results in a layer immediately adjacent to the growing crystals that is generally depleted in compatible elements and enriched in incompatible elements (Kent 2008). If these boundary layer melts are trapped as MI, their chemical compositions are obviously not representative of the bulk melt. It should be noted, however, that: (1) examples of boundary layers have been mainly observed in experimental runs, whereas many studies on natural classic MI do not report compositions affected by boundary layer phenomena—see review papers of Kent (2008) and Audétat and Lowenstern (2014), and (2) experimental boundary layers commonly formed at high degrees of undercooling (~100–200 °C) around skeletal and dendritic crystals (e.g., Faure and Schiano 2005; Morgan and London 2005b). In particular, Faure and Schiano (2005) have demonstrated that the melt trapped in polyhedral crystals is representative of the parental melt, whereas the composition of MI hosted in skeletal or dendritic forsterite crystals lies away from the predicted liquid line of descent. It is important to stress again that MI in high-grade metamorphic rocks are trapped along the heating path in polyhedral peritectic crystals (as opposed to the MI formation during cooling/undercooling in crystallizing magmas). Results from a thorough study of primary MI hosted in garnets of migmatites from the Ojén unit (Betic Cordillera, S Spain) and formed at increasing temperatures (from ~700 to ~820 °C), show that MI compositions are mostly controlled by melting conditions, and not by boundary layer systematics. Thus, Fe (compatible with respect to the host garnet) and K (incompatible) are enriched in MI formed at higher temperatures, whereas low-temperature MI have higher H and Na (both incompatible) (Bartoli et al. 2014, 2015). Therefore, the boundary layer effect does not seem to have played an important role in shaping the geochemistry of these trapped melts. Likewise, mass balances show that boundary layer phenomena do not control the differences in major elements compositions among MI hosted in different minerals of anatectic enclaves of SE Spain (Cesare et al. 2015). Regarding the trace element concentrations, MI in several plagioclase and garnet crystals show remarkably similar concentrations in incompatible elements (e.g., Li, Rb, Cs, B, Be, Zn, As, Zr, Th, U), indicating that boundary layer phenomena are negligible for these elements. Conversely, the trace elements compatible or very slightly incompatible in the host (Sr, Ba, and Eu for plagioclase, and Y and HREE for garnet) are depleted in the MI due to interaction with the host phase (see Acosta-Vigil et al. 2012). In conclusion, the available data seem to indicate that MI are representative of the bulk melt for the major and incompatible trace elements.

Recently, the reliability of olivine-hosted classic MI in preserving the pre-eruptive H2O contents of degassed lavas has been challenged by experimental studies, which have induced huge H2O variations in MI within olivine (e.g., from ~3.8 wt% H2O to ~800 ppm; Hauri 2002; Massare et al. 2002; Portnyagin et al. 2008; Gaetani et al. 2012). On the other hand, Audétat and Lowenstern (2014) have noted that most of the above experiments were done under unrealistic conditions and that many studies on natural MI in olivine demonstrate that the original volatile content of the trapped melt can be preserved after entrapment. Without entering in this “magmatic” dispute, we note that the two contrasting mechanisms of primary MI entrapment (i.e., on the liquidus during crystallization of magmas, and on or close to the solidus during incongruent melting of crustal rocks) along with the contrasting P-T histories may result in different processes affecting the H2O content of the trapped melts. Compared with crystalline rocks, which have experienced prograde anatexis, olivine-hosted MI within magma batches represent the most suitable system for diffusive H2O re-equilibration. Indeed, during ascent of magmas, important gradients of H2O concentration and pressure can be continuously formed between the MI and the external magma. On the contrary, in the partially melted crust MI-bearing peritectic phases coexist with discrete fractions of melt scattered in a mostly solid matrix. As discussed by Bartoli et al. (2014), the large differences in physical parameters (e.g., temperature and, in turn, hydrogen diffusivities, pressure, and size of MI-bearing crystals) can result in times for the diffusive H2O re-equilibration that are >10 orders of magnitude greater in migmatitic and granulitic terranes than in the plumbing of magmatic systems.

It might be argued that if crustal rocks undergo a long thermal history, reaching temperatures much higher than the trapping temperatures and producing higher amounts of melt, H2O contents in MI could re-equilibrate with the external drier melt. In this case, three important observations must be taken into account. First, there is much field, microstructural, and geochemical evidence indicating that most of the crustal melts generated in the source area have been extracted (Brown 2013). Second, melt accumulation and extraction from anatectic rocks is expected to occur very early after crossing the solidus and to continue during the prograde path, at very low melt fractions (<0.07; Rosenberg and Handy 2005; Jamieson et al. 2011), leaving residual granulites enriched in peritectic minerals (Brown 2013). Third, very rapid rates of melt extraction have been reported in the literature. For example, Harris et al. (2000) and Villaros et al. (2009) calculated residence times of S-type granites at the source as short as 50 and 500 yr, respectively. Therefore, high-temperature conditions do not implicitly mean that MI in residual migmatites and granulites experienced diffusive H2O re-equilibration. The successful re-homogenization of nanogranitoids at temperatures similar to those inferred by phase equilibria calculations or thermobarometry, as reported by Bartoli et al. (2013a, 2013c) and Ferrero et al. (2015), strongly supports the idea that MI can behave as a closed system. Rather, MI compositions could be compromised by fluid infiltration, hydration, and retrogression occurring during the subsolidus low-temperature evolution in crystalline basements. However, these processes leave clear microstructural evidence such as the presence of low-temperature minerals within the nanogranitoid inclusions like chlorite (replacing biotite) and pyrite [see Fig. 19 in Cesare et al. (2015)]. All the above observations strongly suggest that the MI database corresponds to that of primary anatectic melts.

Major and trace elements, and H2O contents

In general, the compiled compositional database for MI in high-grade metamorphic rocks reports compositions rich in SiO2 (~64–79 wt%; on a hydrous basis), with variable Al2O3 (~10–18 wt%), Na2O (~1–6 wt%), and K2O contents (~0–8 wt%) (Fig. 4). Most MI have very low to low-moderate FeOt+MgO+TiO2 concentrations (~0.25–2.5 wt%), whereas a small fraction of the data set (~10%) show values up to ~4.0–5.0 wt%. CaO is low (<1.5 wt%), except for MI from Kali Gandaki where CaO content can reach values as high as 4 wt% (Fig. 4). H2O shows extremely variable concentrations (0–12 wt%) that are related to different P-T conditions and fluid regimes during crustal melting.

Regarding the trace elements, to date the data set is limited to two case studies from S Spain, i.e., crustal enclaves from El Hoyazo and migmatites from Ojén. Despite the different tectonic settings, anatectic history, and host rocks and crystals, the trace elements patterns show remarkably similar characteristics when compared with the composition of the upper continental crust: MI are generally enriched in elements (Cs, Rb, Li, Be, and B) partially to totally controlled by muscovite and biotite, and depleted in elements (Ba, Th, REE, Sr, Zr, and Hf) hosted by feldspars and accessory phases (zircon and monazite) [Fig. 5; see also Fig. 13 of Acosta-Vigil et al. (2012)]. Overall, MI from Ojen migmatites are characterized by higher contents of Cs, Rb, B, U, Nb, Zr, Sm, V, and Zn, lower contents of Ba and Pb, and variable concentrations of Th, La, and Ce with respect to MI in NVP enclaves (Fig. 5). This may reflect differences in bulk host rock compositions and/or mechanisms of anatexis [i.e., different melting reactions, different extent of equilibration between melt and residuum; see Acosta-Vigil et al. (2010, 2012)].

One composition, different types of granitoids

Approximately 20 different classification schemes for granitoid rocks have been proposed in the last 40 yr (Barbarin 1999; Frost et al. 2001 and references therein). To describe the composition of crustal magmas trapped in peritectic phases of high-grade metamorphic rocks, we have followed some of the most commonly used classification schemes. Based on their normative compositions and considering the compositional fields of O’Connor (1965), the majority of MI analyzed up to date correspond to granites, whereas nanogranitoids from Kali Gandaki (Nepal) plot in the granodiorite, trondhjemite, and tonalite fields (Fig. 6). MI are generally slightly [ASI = 1.0–1.1; ASI = mol. Al2O3/(CaO+Na2O+K2O)] to strongly peraluminous, with ASI up to 1.5 (Fig. 7). However, some proportion of the MI enclosed in orthogneisses of the Bohemian Massif and tonalites of La Galite, are metaluminous or located at the boundary between the metaluminous and peralkaline fields (ASI = 0.8–1.0, alkalinity index up to −0.02) (Fig. 7). Chappell and White (1974) proposed the first, modern geochemical scheme for the classification of granitic (s.l.) rocks, and recognized two distinct granitoid types: the S-type (supracrustal type) derived from metasedimentary rocks and the I-type (infracrustal type) inferred to have formed from metaigneous sources (see also Chappell 1984, 1999; Chappell and White 1992, 2001). Considering their ASI, most analyzed MI corresponds to S-type granitoids, with only 8% of the MI plotting in the I-type field (Fig. 7). Remarkably, almost the totality of MI produced from anatexis of metaigneous rocks correspond to I-type granitoid compositions (Fig. 7).

Frost et al. (2001) introduced a geochemical classification for granitoid rocks that has achieved wide use (see also Frost and Frost 2008). Using the compositional parameters suggested by these authors, MI are mostly ferroan granitoids, from alkali-calcic to calc-alkalic, even though some spread can be observed (Fig. 8). MI from Kali Gandaki rocks, instead, depart markedly from the other crustal melts and plot in the calcic field (Fig. 8). As stated by Frost and Frost (2008), calcic compositions are typical of melts produced by fluid-present melting.

Insights on the origins and evolution of granitoid magmas

Melt inclusions and the model (haplo)granite system

The origin of granitoid rocks has been strongly debated in the middle half of the past century by N.L. Bowen and H.H. Read (see Bowen 1948; Gilluly 1948; Read 1948). This controversy ended a decade later, when the magmatic origin of granite (s.l.) was definitely demonstrated via experimental petrology by Tuttle and Bowen (1958). These authors showed that melting of a ternary mixture of quartz, albite, and orthoclase—i.e., a simple granitic system called haplogranite (Qtz-Ab-Or), starts at relatively low temperatures, typical of the continental crust, provided that the system is H2O-saturated. After the pioneering work of Tuttle and Bowen (1958), an impressive number of studies have been performed to understand the effect of pressure and H2O, and to a lesser extent of the other minor components, on the phase relations of the haplogranitic system (Johannes and Holtz 1996; see below). Also, whole-rock data have been often compared to the eutectic/cotectic compositions of haplogranite system to identify primary compositions of crustal melts (e.g., Olsen and Grant 1991; Symmes and Ferry 1995; Slagstad et al. 2005) or to infer the depth of magma-storage (e.g., Blundy and Cashman 2001; Almeev et al. 2012; Gualda and Ghiorso 2013).

Because the haplogranitic system is a relevant reference frame to discuss granite petrogenesis, and because MI in high-grade metamorphic terranes seems to provide the compositions of the primary melts produced by melting of natural crustal rocks, it is important to see how MI data plot in the pseudoternary normative Qtz-Ab-Or diagram. Here, granitoid MI from each locality show some spread and a distinctive composition with respect to each other (Fig. 9a). They are commonly located toward the center of the diagram, often in the Qtz field and above the corresponding cotectic lines at the pressure of melting, between the H2O-saturated eutectic point and the Qtz-Or join (Fig. 9a). Considering the relatively wide compositional spectrum of the host rocks, it is clear from Figure 9a that the concept of haplogranitic “minimum melt” composition finds little applicability in natural multi-component systems, as even low-temperature melts may significantly depart from eutectic compositions. Indeed, natural rocks are more complex than a model system such as the haplogranite, and additional components such as Ca, Ti, Fe, and Mg play an important role in modifying the liquidus and solidus phase relationships. For example, Wilke et al. (2015) noted that the effect of An on the phase relationships of the haplogranite system is much more pronounced if small amounts of Fe and Ti (~1 and ~0.2 wt%, respectively) are present.

In Figure 9b we have summarized the effects of P, H2O, and other components (e.g., Al and Ca) on the position of the minimum and eutectic points and of the cotectic lines in the ternary system Qtz-Ab-Or, as reported in the literature. MI often cluster above the corresponding cotectic lines, which, in accordance with experimental studies, suggests that excess Al and Ca play an important role in moving the eutectic compositions and the position of cotectic lines toward the Qtz apex (Figs. 9a and 9b). Most MI show some spread in the Qtz/feldspar ratio compared to their small variation in Ab/Or, forming trends at relatively high angle with the cotectic line (Fig. 9a). This spread of MI data parallel to the Qtz-feldspar direction suggests a control by diffusion in the melt for the major element compositions at the time of MI entrapment (Fig. 9b). Indeed, it has been experimentally shown that the sluggish diffusion of Si and Al vs. rapid diffusion of alkalis in granite melts produced by partial melting of quartzo-feldspathic protoliths result in linear trends in the Qtz-Ab-Or diagram [see Fig. 10 in Acosta-Vigil et al. (2006)].

The decrease of aH2O moves the eutectic point toward the Qtz-Or join without affecting the Qtz content (Fig. 9b). The majority of MI are located between the composition of the H2O-saturated eutectic and the Qtz-Or side. Although this shift toward higher Or/Ab ratios can be partly related to the presence of Ca in natural systems (Fig. 9b), the collected data support the widespread idea that anatexis and formation of crustal melts generally occur under H2O-undersaturated conditions (e.g., Thompson 1990; Clemens and Watkins 2001). It is important to note, however, that the absence or low proportion of H2O-saturated melts does not represent a proof for fluid-absent melting, as even H2O-fluxed melting is generally expected to produce final melts undersaturated in H2O (Weinberg and Hasalovà 2015b).

The reduction of the aH2O during anatexis can be related to two different scenarios: (1) melting along a prograde path, in which all the free H2O present in the rock below the H2O-saturated solidus is dissolved in the melt at the solidus and melting then continues through fluid-absent reactions, and (2) presence of other components in excess fluids, like the carbonic species CO2 and CH4, that reduce the activity of H2O. A clear example of scenario 1 is represented by the Ojén MI. The peraluminous metagraywacke of Ojén underwent prograde anatexis and MI in metatexites represent the first melt produced immediately beyond the fluid-saturated solidus at ~660–700 °C, ~4.5–5 kbar (Bartoli et al. 2013c), whereas MI in diatexites reflect the composition of melt formed at higher temperature, under H2O-undersaturated conditions (~820 °C; Bartoli et al. 2015). In agreement with this evolution and with the bulk-rock composition, MI in diatexites show higher Or/Ab compared to MI in metatexites (Fig. 9a; see also Bartoli 2012), and the H2O content of MI decreases from ~7 to ~3 wt% with increasing temperature (Bartoli et al. 2015).

In scenario 2, the presence in the fluid of carbonic species characterized by low solubility in granitoid melts can result in the phenomenon of immiscibility. As a matter of fact, cordierite, garnet, and plagioclase crystals from enclaves of the NVP (SE Spain) preserve spectacular microstructural evidence of immiscible trapping (Fig. 3f), where granitic MI coexist with primary fluid inclusions that are CO2-dominated (>85 mol%), with minor amounts of N2 and CH4, and traces of CO and H2 (see Cesare et al. 2007; Ferrero et al. 2011). Fluid-melt immiscibility is also documented in metaluminous granitoids at La Galite, Tunisia (Ferrero et al. 2014). The origin of CO2-rich fluids in the deep crust is strongly debated (Santosh and Omori 2008; Huizenga and Touret 2012). In the graphite-bearing enclaves from NVP the occurrence of CO2 has been interpreted as internally generated by Fe3+ reduction and graphite oxidation during incongruent melting of biotite (Cesare et al. 2005). Conversely, infiltration of CO2 from the mantle or from other crustal sources has been postulated for rocks from La Galite (Ferrero et al. 2014).

MI from Kerala Khondalite Belt (India) plot very far from haplogranitic minimum melts, close to the Qtz-Or sideline and point to melting temperatures well above minimum or eutectic temperatures, in agreement with the ultrahigh-temperature (>900 °C) origin of the rocks and the very low-H2O content (~0.7 wt%; Cesare et al. 2009) of the MI. Similar compositions have been obtained by experimental melting of pelitic protoliths at 5–7 kbar, 900–950 °C (Fig. 9a). On the other hand, nanogranitoids from Kali Gandaki (Nepal) are clearly displaced toward the Qtz-Ab sideline and reflect the compositions of crustal melts produced by H2O-fluxed melting (Carosi et al. 2015). The latter inclusions show the highest H2O contents (up to 12 wt% by EMP differences) and CaO, and the lowest K2O (Fig. 4), partly overlapping the composition of experimental glasses produced in the presence of excess H2O at variable P-T conditions (6–14 kbar, 675–775 °C; Fig. 9a). Indeed, the presence of free H2O decreases the plagioclase + quartz solidus more strongly than it depresses the stability of micas, producing tonalite/granodiorite/trondhjemite melts (Weinberg and Hasalovà 2015a and references therein). Melt inclusions from the orthogneiss of the Bohemian Massif, re-homogenized at 875 °C and 27 kbar, show the highest content of normative Ab (Fig. 9a) explained by melting at near ultrahigh-P conditions (increasing pressure moves haplogranitic eutectics closer the Ab apex; Fig. 9b). Experimental melts produced at 24–30 kbar, 870–950 °C by melting of different crustal protoliths plot in the same area (Fig. 9a; Ferrero et al. 2015).

Characterizing the melting processes at the source

Granitoid MI in high-grade metamorphic rocks are providing more and more robust and consistent information on the chemistry of natural crustal melts (see above). In addition, they represent a new tool to constrain mechanisms (i.e., melting reactions, melting mode, extent of equilibration with the solid residue) and conditions (P-T-aH2O) of formation of crustal magmas, and hence open a new window into the anatectic history of the partially melted continental crust (Fig. 1; Acosta-Vigil et al. 2010).

An example of this novel approach is the study of enclaves from El Hoyazo (Neogene Volcanic Province, SE Spain; Acosta-Vigil et al. 2007, 2010, 2012). These rocks represent fragments of continental crust that underwent partial melting and extraction of about 30–60 wt% of granitic melt (Cesare et al. 1997; Cesare 2008). Here, melt (undevitrified glass owing to rapid cooling of the host dacitic lava) has been chemically investigated both within plagioclase and garnet (see compilation of Fig. 5) and in the rock matrix where it often occurs along foliation planes (i.e., the equivalent of leucosomes in regional migmatites and granulites). Melt inclusions show the highest concentrations of Li, Cs, and B, whereas the matrix glass is characterized by higher contents of first row transition elements (FRTE), Y, Zr, Th, REE, and higher Rb/Cs and Rb/Li (Fig. 10a; see Acosta-Vigil et al. 2010). The concentrations of Zr and LREE were used to calculate the zircon and monazite saturation temperatures (Watson and Harrison 1984; Montel 1993), to link MI and potential melt-producing reactions. This approach represents an important advance in granitoid petrology, as the initial temperature of magmas at the source has been long considered to be an inaccessible parameter for allochthonous granites (see Miller et al. 2003). Melt inclusions in plagioclase yield the lowest temperatures (~665–715 °C), whereas MI in garnet (~685–750 °C) and matrix glass (~695–800 °C) show higher temperatures (Fig. 10b). Because the analyzed melts seem to be undersaturated to some extent in the accessory phases (Acosta-Vigil et al. 2012), the obtained temperatures have to be considered as minimum estimates (Miller et al. 2003; Acosta-Vigil et al. 2010). In summary, based on the composition of MI within several hosts and matrix glasses, calculated accessory mineral saturation temperatures, and microstructural relationships, it was concluded that the deep crust beneath the Neogene Volcanic Province of SE Spain underwent rapid heating and melting: MI in plagioclase reflect the earliest granitic melts produced by fluid-present to fluid-absent muscovite melting, whereas MI in garnet were produced simultaneously to slightly later via fluid-absent melting of muscovite (Fig. 10c). In contrast, the generation of granitic melts characterized by higher contents of FRTE and REE (the matrix melt) occurred at higher temperatures by biotite dehydration-melting reactions (Fig. 10c). In all cases melt and residual minerals were at, or close to, equilibrium with respect to most trace elements, except for garnet, and some undersaturation in accessory minerals (Acosta-Vigil et al. 2012).

The “starting-point” composition

To understand which granitoid rocks reflect primary crustal melts and, subsequently, what changes crustal magmas have undergone either near the source region or during their journey toward Earth’s surface, “starting-point” compositions are required (Sawyer et al. 2011). Because leucosomes commonly do not represent primary melts (Sawyer 2008 and discussion therein), a widespread approach among crustal petrologists is to assume the composition of experimental glasses as representative of crustal melt composition at the source, and to use such composition in mass-balance calculations to track the processes involved in magma evolution and crustal differentiation (e.g., Milord et al. 2001; Solar and Brown 2001; Barnes et al. 2002; Guernina and Sawyer 2003; Stevens et al. 2007; Sawyer 2008; Hacker et al. 2011).

The fundamental advance provided by the study of nanogranitoid MI in high-grade metamorphic rocks is that they make accessible the composition of the natural primary anatectic melt for the specific rock under investigation at its specific P-T-XH2O conditions (Bartoli et al. 2013a; Ferrero et al. 2015; this work). Despite the number of MI analyses is already significant (see above), the compiled MI database should be expanded to cover other potential sources of crustal magmas (e.g., intermediate to mafic protoliths) and a wider spectrum of P-T-aH2O conditions under which continental crust may melt (i.e., to date the number of MI case studies is still quite low compared to the number of experimental studies). However, the already significant MI data set collected so far can be already useful to start discussing some of the previously proposed inferences and models regarding crustal melting and the formation and evolution of granitoid magmas. Clearly, to make inferences on petrogenetic processes affecting anatectic magmas, and in the absence of specific primary melt compositions from MI for their particular study area, authors should consider using the MI compositions in the published data set for the conditions ( P-T-XH2O-bulk rock composition) closest to the particular rocks that they are investigating.

Figures 11a–11d provides the compositional variations of MI in terms of mol% of FeTot+Mg+Ti, Na+Ca, K, and Si+Al. The major crystalline phases are reported along with MI data. These types of diagrams have been previously used to model crustal reworking by anatexis (see Solar and Brown 2001; Barnes et al. 2002; Korhonen et al. 2010; Farina et al. 2015; Yakymchuk et al. 2015). Host rock compositions form different clusters according to their nature (for example, metapelites show the highest FeTot+Mg+Ti values, whereas metagraywackes are enriched in Si and K); MI plot away from the FeTot+Mg+Ti apex, into the bottom half of the diagrams at variable (Na+Ca)/K and (Na+Ca)/(Si+Al) ratios, providing a slightly variable though distinctive melt composition (Figs. 11a–11d). This is in agreement with the fact that the composition of the source and conditions of anatexis play a major control on melt chemistry.

The Himalayan Leucogranites (HL) are commonly considered to represent nearly pure melts produced by low-temperature (~700 °C) melting of metapelites and metagraywackes (Le Fort et al. 1987). As a matter of fact, most HL overlap with the low-temperature MI (i.e., those hosted in garnet from the Ronda metatexites and in garnet and plagioclase of NVP enclaves), supporting their primary nature (Figs. 12e and 12f). Some HL compositions, however, do not resemble MI compositions and plot in the proximity of the lower portion of the field occupied by most MI data (Figs. 11e and 11f). It is important to note that the difference in mafic components (Fe, Mg, and Ti) between HL and MI cannot be related to contamination of MI by the host phase during analysis, as MI hosted in Fe-, Mg-, and Ti-free hosts such as plagioclase show similar and sometimes higher Fe+Mg+Ti contents with respect to MI hosted in garnet (Figs. 11b and 11d). On one hand, the observed discrepancy can be related to an incomplete MI data set (i.e., incomplete sampling of MI in high-grade metamorphic rocks). On the other hand, it could indicate a non-primary nature of some HL. For instance, Scailett et al. (1990) stated that some geochemical features of Gangotri leucogranites (for which compositions sometimes depart from those of MI; Fig. 11e) can be explained by fractionation of the early crystallizing phases. Their compositions often depart from the compositions of primary HL in the opposite direction of the Bt and Ms fields (Figs. 11e and 11f). As a matter of fact, some Gangotri leucogranites have lower Co, V, and B contents than MI in Ronda metatexites and lower Li content than MI in garnet of NVP, in agreement with a mica fractionation process. It should be noted that petrographic and experimental evidence indicates that biotite was likely the first mineral to crystallize in many HL (Scaillet et al. 1995). Therefore, the available MI data seem to support a primary nature for the majority of HL, even though fractional crystallization processes might have occurred locally, in agreement with results from previous studies (i.e., Le Fort et al. 1987; Scaillet et al. 1990). It is important to note that MI from Kali Gandaki (Nepal) represent Ca-rich melts produced during Eocene prograde metamorphism (Carosi et al. 2015; Iaccarino et al. 2015) and, therefore, they are not related to the Miocene HL formed during exhumation stage (Harris and Massey 1994).

The Layos Granite (Toledo Complex, central Spain) represents a S-type, strongly peraluminous Crd-bearing granitoid suite, varying from quartz-rich tonalite to melamonzogranites and formed at about 800–850 °C by biotite melting of metapelites (Barbero and Villaseca 1992). Most of these rocks plot away from MI data sets, toward the FeTot+Mg+Ti apex (Figs. 11e and 11f). The authors explained the observed compositional variability of the suite by a “restite unmixing” model (Barbero and Villaseca 1992). This model was proposed by Chappell et al. (1987) to explain the geochemical variability of granitoid rocks and states that the majority of granites represent mixtures between a low-T melt and the restite, i.e., solid materials that were residual from the source. Therefore, granites were thought to image their source in a simple way: compositional variations observed within granitoid suites reflect varying degrees of restite unmixing with a minimum or near-minimum temperature (haplogranitic) melt (Chappell et al. 1987). However, this model has received several criticisms from the geochemical, petrographic, and physical point of view (see Wall et al. 1987; Clemens and Mawer 1992; Clemens 2003; Vernon 2010; Clemens and Stevens 2012 for additional details). The fact that the composition of the primary melt produced in the source may differ markedly from that of minimum melt, as clearly demonstrated by MI study (see discussion above), demonstrates that the assumption that the melt component has an eutectic composition is not applicable in natural systems. Recently, Stevens et al. (2007) and Clemens and Stevens (2012) proposed a revised version of the restite unmixing model, called “peritectic assemblage entrainment” or “selective peritectic phase entrainment.” Following this model, the more mafic granitoid rocks result from the entrainment of the solid peritectic products of the melting reaction to the melt (although how this would occur is not fully explained), rather than from the entrainment of a proportion of the total solid mineral assemblage remaining after the melting reaction. Regarding the Layos Granite, only one data point plots in the MI field (overlapping the composition of MI trapped in Crd from metapelites and formed at ~850 °C; Figs. 11e and 11f) and, in turn, seems to reflect a primary origin. The rest of data plots between this primary melt composition and the FeTot+Mg+Ti apex and does not form a trend toward the (Si+Al)/10 apex. Similarly, in the (FeTot+Mg+Ti)–(Si)–(Al) diagram (not shown), Layos granites do not show compositional trends toward the Si and Al apexes. According to the model proposed by Stevens et al. (2007), this chemical signature can be explained by differential entrainment of peritectic cordierite (up to 30% of Crd has been documented in these rocks), without the involvement of residual Qtz and Als, i.e., the restite unmixing model does not seem to be applicable in this specific case.

The Malani rhyolites (western India) are considered to be the product of high-T crustal anatexis (Maheshwari et al. 1996), and the fact that the majority of lava compositions match those of MI supports their primary origin (Fig. 11c). Strongly peraluminous, K-rich Malani rhyolites resemble the glassy MI from Kherala Khondalite Belt, in agreement with the inference that the rhyolites were likely produced by melting of metapelites at temperatures ≥850 °C (Maheshwari et al. 1996). At the same time, the comparison with MI would suggest a more complex scenario for the genesis of these anatectic lavas. Malani rhyolites have variable (Na+Ca)/K ratio and some lava compositions overlap the compositions of MI found in NVP metapelitic enclaves, Ronda metatexites and La Galite granodiorite and formed at variable temperatures (~700–820 °C) (Fig. 11c). These lavas range from peraluminous, through metaluminous to peralkaline. Notably, almost the totality of metaluminous/peralkaline lavas approach the metaluminous/peralkaline compositions of MI from La Galite granodiorite. The observed compositional variability can reflect the presence of distinct magma batches in the Malani rhyolite complex, formed and extracted during prograde heating of a compositionally heterogeneous crustal source. In this sense, other work has verified the occurrence of a heterogeneous crystalline basement (composed of both granodiorites, orthogneisses, and metasedimentary rocks) beneath Malani rhyolites (Pandit et al. 1999; Sharma 2004). Instead a single composition seems to reflect the entrainment of mafic minerals (Fig. 11c).

In Figure 12 we have compared the composition of MI with a data set of ~520 experimental glasses reported in literature and produced at variable P-T-aH2O conditions and using different starting materials, i.e., metasedimentary and felsic (meta)granitoid protoliths. Remarkably, MI and glasses from experiments seem to be characterized by similar compositions and compositional ranges, even though experimental melts may show higher FeTot+Mg+Ti values than MI (Fig. 12). Koester et al. (2002) reported glasses, produced at 700 °C and 5–15 kbar, characterized by anomalously high FeO, TiO2, and MgO contents (e.g., FeO ~9.5 wt%; see red dots in Fig. 12), in contrast with their low-temperature origin. The authors argued that the lack of chemical equilibrium was responsible for the heterogeneous chemical compositions of melt observed in near-solidus experiments (Koester et al. 2002; see also Gardien et al. 1995). Indeed, as any other petrologic tool, partial melting experiments may suffer from potential pitfalls, in particular: impossibility of bulk compositional changes (in experiments there is no H2O and melt loss), slow kinetics of reactions, loss of iron from silicate phases to experimental capsules, lack of equilibration, overstepping of melting reactions and difficulty of obtaining reliable glass analyses in near-solidus experiments (Ratajeski and Sisson 1999; Gardien et al. 2000; White et al. 2011; Webb et al. 2015). Nonetheless, when the low-temperature (≤720 °C) experiments are considered, with the exception of the analyses reported in Koester et al. (2002), the produced melts show low and variable FeTot+Mg+Ti contents, overlapping the compositions of low-temperature MI (Fig. 12).

From all the above considerations, it is clear that MI are an useful tool with which to discuss the petrogenesis of granitoid rocks, and can complement experimental glasses to constrain differentiation processes occurring in anatectic magmas. Researchers working on petrogenesis of crustal granitoids can now look for and study MI in migmatites and granulites considered to be the source region of the studied magmas, or in peritectic garnet entrained in the granitoid rocks to fix the real primary composition for their specific case study. In particular, as the number of experimental studies presently reporting trace element composition of melt is limited, MI are an important contribution that provides the trace element compositions of natural anatectic melts.

Implications

Nanogranitoid melt inclusions in high-grade metamorphic rocks are a novel subject of research and are becoming a matter of considerable scientific interest, as they allow the analysis in situ of the natural primary melts produced by crustal melting. The size of nanogranitoid MI (commonly 10 μm in diameter) and of daughter minerals crystallized in them (often ≤1 μm) pose analytical challenges, and the accurate re-homogenization and compositional characterization of these small data repositories is far from being routine. However, we have shown here that, albeit slowed down by many experimental and analytical difficulties, the study of these MI opens new horizons for metamorphic and igneous petrology. Technological improvement in the coming years will make the study of MI easier, and we believe that their chemical characterization, in particular of volatiles, trace elements, and radiogenic isotopes, will provide exciting new results in many topics, such as granitoid petrogenesis, UHT and UHP metamorphism, secular variations of felsic magmatism, volatile recycling, crust formation, and silicate Earth differentiation.

Acknowledgments

Many thanks to Edward Sawyer and Chris Yakymchuck for their reviews of the manuscript and to Cal Barnes for the editorial handling. This research benefitted from funding from the Italian Ministry of Education, University, Research (Grant SIR RB-SI14Y7PF to O.B. and Grant PRIN 2010 TT22SC to B.C.), from Padova University (Progetti per Giovani Studiosi 2013 to O.B., Progetto di Ateneo CPDA107188/10 to B.C. and Piscopia—Marie Curie Fellowship GA No. 600376 to A.A.V.), from the Ministerio de Ciencia e Innovación of Spain (grant CGL2007-62992 to A.A.V.), and from the Alexander von Humboldt Foundation, the German Federal Ministry for Education and Research and the Deutsche Forschungsgemeinschaft (Project FE 1527/2-1 to S.F.). The research leading to these results has received funding from the European Commission, Seventh Framework Programme, under Grant Agreement no. 600376.

1
Deposit item AM-16-75541, Supplemental Tables. Deposit items are free to all readers and found on the MSA web site, via the specific issue’s Table of Contents (go to http://www.minsocam.org/MSA/AmMin/TOC/).
Manuscript handled by Calvin Barnes