Xisheng Xu
email: xsxu@nju.edu.cn

Saskia Erdmann
email: saskia.erdmann@cnrs-orleans.fr

Kai Zhao

T A Mather

Thermodynamic modeling for an incrementally fractionated granite magma system: Implications for the origin of igneous charnockite

Keywords:

Thermodynamic modeling for an incrementally fractionated granite magma system: Implications for the origin of igneous charnockite

Introduction

Magma fractionation plays an important role in the forma-tion and differentiation of the continental crust (e.g. [START_REF] Keller | Volcanicplutonic parity and the differentiation of the continental crust[END_REF], and in the enrichment of elements of economic interest (e.g. [START_REF] Mustard | Gold and metal enrich-ment in natural granitic melts during fractional crystallization[END_REF]. Many studies on low-viscosity, mafic layered magma series have suggested that cumulate formation by gravitational settling is a key process of magma fractionation (e.g. [START_REF] Wager | Types of igneous cumulates[END_REF]. Other studies have advocated that cumulates form by in situ crystallization in solidification fronts and subse-quent removal of evolved melt (e.g. [START_REF] Latypov | Plagioclase compositions give evidence for in situ crystallization under horizontal flow conditions in mafic sills[END_REF][START_REF] Holness | Assessing the role of compaction in the formation of adcumulates: a microstructural perspective[END_REF]. Silicic magmas have been considered to be largely stagnant during much of their lifetime in crustal reser-voirs, given their high viscosities (e.g. [START_REF] Huber | Homogenization processes in silicic magma chambers by stirring and mushification (latent heat buffering)[END_REF][START_REF] Glazner | Magmatic life at low Reynolds number[END_REF], but vigorous magma-flow-driven crystalliquid segregation of crystal-rich magma is at least locally recorded [START_REF] Paterson | Interpreting magmatic fabric patterns in plutons[END_REF][START_REF] Zak | Four magmatic fabrics in the Tuolumne batholith, central Sierra Nevada, California (USA): implications for interpreting fabric patterns in plutons and evolution of magma chambers in the upper crust[END_REF]. The formation of crystal-poor rhyo-lites or other evolved melts (e.g. aplites) requires efficient seg-regation of melts from crystal mush through hindered settling, compaction and/or filter pressing [START_REF] Bachmann | On the origin of crystal-poor rhyolites: ex-tracted from batholithic crystal mushes[END_REF][START_REF] Gelman | Identifying the crystal graveyards remaining after large silicic eruptions[END_REF]. Direct microstructural evidence for melt ex-traction and crystal accumulation, however, is rarely examined in intrusive magma bodies [START_REF] Fiedrich | Mineralogical, geochemical, and textural indicators of crystal accumulation in the Adamello Batholith (northern Italy)[END_REF], and mechanisms of melt-crystal segregation thus commonly remain hypothetical [START_REF] Holness | Melt segregation from silicic crystal mushes: a critical ap-praisal of possible mechanisms and their microstructural record[END_REF].

Thermodynamic modeling (e.g. the using of the rhyolite-MELTS, Magma Chamber Simulator, or Perple_X program packages) is a powerful tool for constraining phase relationships, intensive crys-tallization parameters, and the compositional evolution of liquids and solids (e.g. [START_REF] Connolly | Computation of phase equilibria by linear programming: a tool for geodynamic modeling and its application to subduction zone decarbon-ation[END_REF][START_REF] Gualda | Rhyolite-MELTS: a modified calibration of MELTS optimized for silica-rich, fluid-bearing magmatic systems[END_REF][START_REF] Spera | Energy-constrained open-system magmatic pro-cesses I: general model and energy-constrained assimilation and fractional crys-tallization (EC-AFC) formulation[END_REF]Zhao et al., 2017a). The basic requirement for the application of most thermodynamic models is that the systems considered have approached total equilibrium (i.e. equilibrium for all chemical components and phases in the system) or at least lo-cal equilibrium (i.e. equilibrium for crystal rims and host liquids) in a compositionally closed system. Mass and energy exchange as a result of open-system processes, e.g. magma recharge, assimilation, and fractional crystallization (leading to melt and/or volatile loss), may cause partial disequilibrium on the mineral scale, e.g. crystal zoning [START_REF] Humphreys | Chemical evolution of intercumulus liquid, as recorded in plagioclase overgrowth rims from the Skaergaard intrusion[END_REF], and deviation of whole-rock compositions from those of primary melts (e.g. [START_REF] Montel | Partial melting of metagreywackes, part II. Compo-sitions of minerals and melts[END_REF]. The open-system evolution of magma systems through the addi-tion of material, e.g. magma recharge or country-rock assimilation [START_REF] Spera | Energy-constrained open-system magmatic pro-cesses I: general model and energy-constrained assimilation and fractional crys-tallization (EC-AFC) formulation[END_REF], has been numerically constrained for many systems. Modeling the compositional and the mineralogi-cal effects resulting from the loss of liquids from evolving magma systems is still rare (with few exceptions, such as the studies of [START_REF] Gelman | Identifying the crystal graveyards remaining after large silicic eruptions[END_REF]Lee and Morton, 2015), yet crucial for quantifying parental magma composition and intensive system variables.

This study investigates such open-system evolution, using the peraluminous Jiuzhou pluton, South China, as an example. The pluton comprises orthopyroxene-bearing granite (charnockite) and orthopyroxenefree granite with assemblages that permit the char-acterization of crystallization pressure, temperature, and melt H2O content. Previous work has shown that the system likely experi-enced partial melt extraction, which was crucial for the preser-vation of orthopyroxene in solidification of granite magma (Zhao et al., 2017a). To further assess possible melt extraction, and to quantify the magmatic system evolution, we have (1) carried out geochemical modeling to evaluate the potential magmatic pro-cesses including magma mixing, restite entrainment and fractional crystallization, and (2) characterized the pluton-and mineral-scale processes of magmatic evolution using whole-rock compositions, mineral textures and compositions, crystal size distribution, mi-crostructure, and thermodynamic modeling utilizing the Perple_X program package of [START_REF] Connolly | Computation of phase equilibria by linear programming: a tool for geodynamic modeling and its application to subduction zone decarbon-ation[END_REF]. Our approach can be applied to constrain the role of melt extraction and cumulate formation in H2O-poor to moderately H2O-rich granitic magma systems, which may contribute new insights into the processes and the conditions of granite fractionation that lead to compositional zoning and de-viation of whole-rock composition from primary melt composition.

Geological background

The Qinzhou Bay Granitic Complex (QBGC), South China, com-prises several strongly peraluminous (ASI > 1.1) granitic intru-sions including the orthopyroxene-bearing to orthopyroxene-free Jiuzhou pluton (Fig. 1). The granitic rocks are classified as S-type granites derived from partial melting of a crustal source with minor or no input of mantle-derived magma [START_REF] Qi | Origin of the Darongshan-Shiwandashan S-type granitoid belt from southeastern Guangxi: geochemical and Sr-Nd-Hf isotopic constraints[END_REF][START_REF] Zhao | Origin of the granulite enclaves in Indo-Sinian peraluminous granites, South China and its implication for crustal anatexis[END_REF][START_REF] Jiao | Metasedimentary melting in the formation of charnockite: Petrological and zircon U-Pb-Hf-O isotope evidence from the Darongshan S-type granitic complex in southern China[END_REF]. Restitic granulite enclaves hosted by the Jiuzhou pluton are inferred to record source partial melting at ∼0.7-0.8 GPa and ∼900-1000 • C and granite emplace-ment at ∼0.2-0.3 GPa (Zhao et al., 2017b).

Mineralogy and petrology of the Jiuzhou pluton

Field relations, microtexture and microstructure

The Jiuzhou pluton is exposed over ∼200 by ∼10 km extend-ing in NE-SW direction (Fig. 1). The charnockites crop out in the southwest part of the intrusion, while the northeast part of the in-trusion exposes orthopyroxenefree granites. The charnockite sam-ples moreover crop out at lower altitude than the orthopyroxene-free granite samples, i.e. at <150 m versus >150 m above sea level (Table S1 in supplementary file 1).

The charnockites are massive, and layering or a preferred fab-ric are not evident on the outcrop scale, while the samples show a weak preferred fabric on the hand specimen-and the micro-scale. Plagioclase crystals are euhedral or subhedral with a low apparent aspect ratio of ∼1.0:1.2:1.8 (R 2 = 0.9). They locally form monomin-eralic clusters (Fig. 2a) and chains (Fig. 3a). Some of the crystals show undulose extinction, mechanical twins (Fig. 2c-e), and/or subgrains (Fig. 2d). Rarely, the crystals are also slightly bent (by <2 • , Fig. 2c). The plagioclase crystal size distribution (CSD) curves show upward concave deflections for large grain sizes (>3.8 mm), except for sample JZ09 with an approximately linear CSD, while they have decreasing grain numbers at crystal sizes of <1 mm (Fig. 4). Quartz crystals are subhedral to anhedral and slightly elongated (Fig. 3a). They commonly show undulose extinction or chessboard subgrain patterns (Fig. 2d). Alkali feldspar crystals have subhedral to anhedral shape with tendrils locally extending into the matrix (Fig. 3a). Some of the crystals show a preferred ori-entation with their long axis sub-parallel to the main plagioclase chains (Fig. 3a). Few of the crystals also show subgrains and/or undulose extinction. Biotite crystals are present as either single crystals with euhedral to subhedral shape or as anhedral crys-tals surrounding orthopyroxene. Some of the euhedral crystals are significantly bent (by up to ∼20 • ) or kinked (Fig. 2f and Fig. S3 in supplementary file 1). Cordierite crystals are locally present as subhedral-euhedral crystals or as inclusions in orthopyroxene, plagioclase, and alkali feldspar. Deformation features in cordierite have not been observed, but they may have been obliterated by its significant sericitization. Orthopyroxene crystals are present as eu-hedral crystals and as subhedral to anhedral crystals with irregular grain boundary outlines against partial or complete rims of biotite (Fig. 2c). Orthopyroxene shows no deformation features.

The orthopyroxene-free granites are massive, showing no layer-ing or a preferred fabric on outcrop-, hand specimen-, or micro-scale. Plagioclase crystals are euhedral to subhedral with a low apparent aspect ratio of ∼1.0:1.3:2.0 (R 2 = 0.8). They show a ran-dom to slightly clustered distribution (Fig. 3b,c). Contrary to the charnockites, plagioclase in the orthopyroxene-free granites shows no mechanical twinning or subgrains. The plagioclase CSDs are ap-proximately linear at crystal sizes of >1 mm, while they show de-creasing grain numbers at crystal sizes of <1 mm (Fig. 4). Quartz crystals are subhedral to anhedral. They show no significant elon-gation or chessboard subgrain patterns, while undulose extinction is common. Alkali feldspar crystals are subhedral to poikilitic. De-formation features have not been observed. Biotite crystals are euhedral to subhedral. Deformation of biotite is rare, but weakly bent biotite with undulose extinction has been observed. Cordierite crystals are subhedral to euhedral, which are partly altered to finegrained sericite and chlorite. Garnet crystals (supplementary table 2) are occasionally present, which have subhedral to anhedral shape and biotite rims. Both cordierite and garnet show no evi-dence for deformation.

Plagioclase and alkali feldspar compositions and zoning

Plagioclase and alkali feldspar core-rim compositions were ac-quired using a JEOL JXA-8800 electron microprobe at Nanjing Uni-versity, China. Details of the method are described in supplemen-tary file 1. Plagioclase in the charnockite (sample JZ01-1) forms anorthite-rich crystals with weak oscillatory and gradual normal zoning, of which some show marginal step-zoning to anorthite-poor overgrowths. The anorthite-rich crystals, including the cores of the step-zoned crystals, have ∼An 60-45 composition (on average ∼An50) (Fig. 5). The anorthite-poor overgrowths have ∼An42-30 and on average ∼An36 composition (Fig. 5). The plagioclase crys-tals thus show a significantly bimodal composition with a compo-sitional gap at ∼An45-42 (Fig. 5c). The anorthite-poor overgrowths are asymmetrical and typically <100 μm wide (Fig. 5a). They are present on one or more crystal surfaces, forming partial rims. They are absent from plagioclaseplagioclase grain boundaries in pla-gioclase clusters and chains (Fig. 2a), but occur against all other mineral types. The plagioclase-plagioclase grain boundaries are characterized by long, mostly straight and at times slightly curved or serrated 2D outlines (Fig. 2a). The anorthite-rich cores occasionally show deformation twins, while deformation twins have not been observed in the overgrowths (Fig. 2e). Crystals with and without the anorthite-poor overgrowths have closely comparable shape and size. The plagioclase crystals with the anorthite-poor overgrowths are locally concentrated in the thin sections, where ∼30% of the crystals show stepzoning (Fig. 3a). Overall, the over-growths make up ∼2 vol% of the charnockite samples. Alkali feldspar also shows core-rim step zoning, where cores and over-growths have ∼Or80 and ∼Or85 composition, respectively.

The overgrowths are partially developed. They make up a larger pro-portion of the samples (∼4 vol%) than the plagioclase overgrowths. and they are locally concentrated around the large phenocrysts (Fig. 3a). The overgrowths are devoid of undulose extinction or subgrains.

Plagioclase in the orthopyroxene-free granite (JZ10) shows sim-ilar zoning patterns as plagioclase in the charnockite (Fig. 5d,e). The anorthitepoor overgrowths are, however, more common and wider than those on plagioclase in the charnockite. They are developed on ∼50-90% of the crystals, typically forming their outer ∼100-300 μm wide rims, amounting to ∼8-10 vol% of the samples. They are predominantly developed as complete rims in 2D, they also occur on plagioclase-plagioclase grain bound-aries (Fig. 2b), and they occur throughout the samples (inset of Fig. 3b,c). The plagioclase crystal cores have ∼An 57-42 (on aver-age ∼An 47 ) composition, and they thus have lower anorthite con-tent than the plagioclase cores in the charnockite. The step-zoned, anorthite-poor overgrowths have ∼An42-30 and on average ∼An35. Cores and overgrowths of the plagioclase crystals do not show a compositional gap as they do in the charnockite, but ∼An39-42 compositions are rare (Fig. 5f). Few of the alkali feldspar crystals also show core-rim step-zoning from ∼Or85 cores to ∼Or90 over-growths.

The cores are completely mantled by the overgrowths, while the majority of the crystals have a composition equivalent to that of the overgrowths.

Overgrowths and unzoned, orthoclase-rich crystals make up ∼10-12 vol% of the orthopyroxene-free granite samples.

Constraints on the geochemical variation of the Jiuzhou pluton

The Jiuzhou rocks define a continuous compositional trend for all major elements from low-SiO2 (63-67 wt%) charnockites to high-SiO2 (68-70 wt%) orthopyroxene-free granites (Fig. 6a-f). The compositions of the charnockites are significantly outside the compositional fields of experimentally produced, metapelite-and metagraywacke-derived partial melts (Fig. 6a-f). The orthopyrox-ene-free granites also differ from, but approach or partly over-lap with the composition of experimental melts (Fig. 6a-f). This suggests that most (or all) of the rocks exposed in the Jiuzhou pluton represent magmas that deviate from primary melt composi-tions. In the following sections, mass balance calculations are used to constrain the plutonscale processes that account for the Ji-uzhou whole-rock compositional variation. The phase compositions required in the calculation were determined by thermodynamic modeling using the Perple_X program [START_REF] Connolly | Computation of phase equilibria by linear programming: a tool for geodynamic modeling and its application to subduction zone decarbon-ation[END_REF]version 6.7.4). Details of the method and modeling results are described in the supplementary file 1.

Possible restite entrainment

Restite enclaves consisting of mainly garnet, cordierite, spinel, biotite, alkali feldspar, and minor orthopyroxene and plagioclase are locally present in the Jiuzhou pluton (Zhao et al., 2017b). The minerals in the enclaves form inclusion-rich porphyroblasts or subhedral-anhedral matrix crystals. Garnet crystals, which only locally occur in the orthopyroxene-free granite, have subhedral to anhedral shape and biotite rims, and they may thus repre-sent restite. The entrainment and bulk assimilation of the garnet-and cordieriterich source-derived restite would have significantly increased whole-rock Al2O3 and FeO and decreased whole-rock K2O and SiO2 contents (Fig. 6af). The addition of 20-30 wt% restite to a magma with composition equivalent to the evolved orthopyroxene-free granites could have produced magmas with SiO2-Al2O3-MgO-FeO-K2O-Na2O compositions equivalent to the charnockites (Fig. 6a-b, d-f), but cannot account for the signifi-cantly higher CaO contents of the charnockites compared with the orthopyroxenefree granites (Fig. 6c).

Possible fractionation in solidification fronts

Fractionation of silicic magmas is thought to be inefficient at low crystallinities (i.e. at crystallinities of <∼40 vol%), but may become efficient when crystals form a framework and magma con-vection ceases (e.g., at crystallinities of >∼40-50 vol%; [START_REF] Bachmann | On the origin of crystal-poor rhyolites: ex-tracted from batholithic crystal mushes[END_REF]. The continuous compositional trends and ver-tical zoning from relative primitive charnockite to more evolved, orthopyroxene-free granite may suggest incremental fractiona-tion, which means that magma may have crystallized in a so-lidification front, which incrementally advanced from the first cooled magma chamber margin into its later cooled center (i.e. in situ fractional crystallization, see box-model and Table S1 in supplementary file 1; [START_REF] Langmuir | Geochemical consequences of in situ crystallization[END_REF][START_REF] Nishimura | In situ fractional crystallization observed in the Os-umi granodiorite batholith[END_REF][START_REF] Hertogen | SilMush: a procedure for modeling of the geochemical evolution of silicic magmas and granitic rocks[END_REF]. Some early solidification fronts may still contain melt when most of the original magma chamber was solid-rich given that significant near-solidus crys-tallization (cotectic, eutectic, or minimum crystallization) will buffer temperature across most of the system (cf. [START_REF] Huber | Homogenization processes in silicic magma chambers by stirring and mushification (latent heat buffering)[END_REF].

To model such incremental fractionation, mass balance equa-tions were solved step-wise fashion. Free parameters used in the modeling are the solidification zone mass fraction (0.1-0.2, cf. Hertogen and Mareels, 2016), the critical crystallinity for melt ex-traction (0.4-0.5, cf. [START_REF] Bachmann | On the origin of crystal-poor rhyolites: ex-tracted from batholithic crystal mushes[END_REF], and the poros-ity of the cumulate pile. The porosity values were determined using thermodynamic modeling and plagioclase compositions (see section 5.2). The other input parameters were determined by ther-modynamic modeling. We highlight that different input values for the three parameters do not, in any case, change the direction of the modeling trends and our general conclusions, but they affect the absolute composition of each modeled increment/solidification zone. The higher the critical crystallinity and porosity values, the more felsic is the composition of the modeled cumulate pile, while the choice of solidification zone mass fraction only determines the coarseness of the modeling (e.g. see also [START_REF] Hertogen | SilMush: a procedure for modeling of the geochemical evolution of silicic magmas and granitic rocks[END_REF]. Further details of the modeling are provided in supplemen-tary file 1, while the results are given in supplementary file 3 and summarized in Fig. 5a-f. The modeling predicts that the cumu-late pile (cumulus solid plus unextracted melt) has ∼62 wt% to 71 wt% SiO2. The evolved magma in the residual chamber compris-ing ∼10% of the mass of the initial magma chamber is modeled to have reached ∼72 wt% SiO2. The most felsic sample of the plu-ton (Fig. 6; [START_REF] Jiao | Metasedimentary melting in the formation of charnockite: Petrological and zircon U-Pb-Hf-O isotope evidence from the Darongshan S-type granitic complex in southern China[END_REF] with ∼72 wt% SiO2 may represent or at least approach the finial liquid-rich differentiates, which have a lower proportion (<∼10 vol%) than the solid-rich cumulates. The modeled compositional variation thus matches the range of the Ji-uzhou whole-rock SiO2 contents, but also the corresponding FeO, MgO, CaO, and K2O contents (Fig. 6b, c andd). The modeled Al2O3 contents are slightly higher (by <2 wt%) than the natural whole-rock contents, while the modeled Na2O contents are slightly lower (by <0.5 wt%) than the natural whole-rock contents (Fig. 6a,f). From the first solidification front to the residual magma chamber, the concentrations of H2O are predicted to have increased from ∼2.0 to 3.5 wt% in the evolving liquids and from ∼0.8 to 3.0 wt% in the bulk cumulates, respectively (Fig. 6g). Samples in the gray fields contain orthopyroxene. (g) The calculated bulk H2 O content of the cumulate piles (cumulus solids plus trapped melt, symbolized by diamonds) and the evolved magma in the residual magma chamber by crosses) for each modeled increment/solidification front.

Discussion

Magma processes recorded in silicic magma systems: examples from the Jiuzhou pluton

Many granitic intrusions worldwide have bulk compositions that deviate from pure liquid compositions (e.g. [START_REF] Montel | Partial melting of metagreywackes, part II. Compo-sitions of minerals and melts[END_REF]. The Jiuzhou pluton is a case example, exposing granitic rocks with ∼63-72 wt% SiO2 that partly differ from the compositions of experimentally produced melts (Fig. 6a-f). Pro-cesses that may have modified the composition of granitic magmas in general include i) magma mixing with mantle-derived melts; ii) entrainment of restitic material; and/or iii) fractional crystal-lization and cumulate formation. We show that these processes can be quantified by combining established geochemical and pet-rographic analysis with thermodynamic modeling.

The role of mantle-derived melt input in the Jiuzhou magma system

Peraluminous granites have mostly sedimentary sources, but many are now known to have experienced mixing with mantle-derived melt (e.g. [START_REF] Stevens | Selective peritectic garnet entrainment as the origin of geochemical diversity in S-type granites[END_REF]. Previous work has shown that mantle-derived melts have made no significant contribution to the strongly peraluminous Jiuzhou pluton with its enriched Sr-Nd-Hf isotopic composition [START_REF] Qi | Origin of the Darongshan-Shiwandashan S-type granitoid belt from southeastern Guangxi: geochemical and Sr-Nd-Hf isotopic constraints[END_REF]Zhao et al., 2017b). The absence of microgranitoid enclaves, the normal zoning of pla- [START_REF] White | Ultrametamorphism and granitoid genesis[END_REF]. f-g , parental magma compositions are calculated based on mass fraction of interstitial melt and whole-rock composition of sample JZ01-1 and JZ10. The composition of interstitial melt is approximate to that of minimum melts. H2 O contents were estimated by our thermodynamic modeling (see supplementary file 3). Fe 2 O3 *: all Fe reported as Fe2O3 .

gioclase, and the lack of patchy zoning in plagioclase (Fig. 5a,d) or any of the other minerals, also preclude significant recharge by mantle-derived melts. The small-amplitude anorthite oscillations in plagioclase cores are likely a record of crystallization during early-stage convective stirring of the magma chamber and/or a reflection of melt boundary layer crystallization (e.g. [START_REF] Humphreys | Chemical evolution of intercumulus liquid, as recorded in plagioclase overgrowth rims from the Skaergaard intrusion[END_REF]. The absence of mantle-derived melt input makes the Ji-uzhou pluton an ideal end-member for quantification with no complication of intermittent variation in composition or other in-tensive parameters.

The role of restite entrainment in the Jiuzhou magma system

Peraluminous granites may contain variable amounts of restite (from a few inherited zircon cores to possibly several wt%; e.g. [START_REF] Stevens | Selective peritectic garnet entrainment as the origin of geochemical diversity in S-type granites[END_REF]. Restitic material is present and locally abun-dant in the Jiuzhou pluton, forming centimeter-size enclaves and micro-restite (Zhao et al., 2017b). The whole-rock geochemical trends, however, cannot be explained by restite entrainment as the restite is too plagioclase-depleted to reproduce the composi-tional variation (Fig. 6c). The volume of the observed micro-restite is moreover small (<3 vol%), while restite entrainment would have had to be on the order of 20-30 wt% to account for the variations in all oxide components (Fig. 6a-f). Samples that do not contain visible restite may comprise largely assimilated restite components (e.g. that were assimilated as a result of decompression), but apart from minor amounts of xenocrysts (e.g. garnet) no evidence for such restite components exists even on the microscale. Orthopy-roxene crystals, including the biotite-rimmed crystals, represent primary magmatic crystals which were partially resorbed by their evolving host melt (Zhao et al., 2017a). We further note that our restitebearing samples (JZ04, JZ03-4) do not have the lowest SiO2 content, while our sample with the lowest SiO2 content (JZ01-1) does not contain obvious micro-restite. We thus posit that restite entrainment has not exerted a major control on the compositional variation of the Jiuzhou pluton. For other granitic systems, it is im-perative to assess the amount and the type of restite present: If restite occurs in the form of xenoliths or xenocrysts, then all work has to focus on samples or sample fractions that are free of these components. If restite was assimilated, then it is crucial when and how much assimilation took place.

The role of fractionation and cumulate formation in the Jiuzhou magma system

Incremental fractionation with cumulate formation and melt extraction can explain the whole-rock compositional variation of the Jiuzhou pluton (as modeled in section 4.2 and shown in Fig. 6) and potentially that of other granitic plutons globally. In the fol-lowing, we discuss the record of mineralscale features. Part of the plagioclase crystals, i.e. ∼30% of the crystals in the charnockites and ∼50-90% of the crystals in the orthopyroxene-free granites, show characteristic step-zoning with subhedral to subround cores and partial to complete, low-anorthite overgrowths. Low-anorthite overgrowths may form by late-or post-magmatic fluid alteration, but this can be ruled out for the Jiuzhou system, as such alter-ation should have produced albite-rich compositions and zones rich in micro pores (e.g. [START_REF] Larsson | Hydrothermal al-teration of plagioclase and growth of secondary feldspar in the Hengill Volcanic Centre, SW Iceland[END_REF], which are not ob-served.

Plagioclase with magmatic step-zoning is, however, also rel-atively common in some mafic igneous cumulates [START_REF] Humphreys | Chemical evolution of intercumulus liquid, as recorded in plagioclase overgrowth rims from the Skaergaard intrusion[END_REF][START_REF] Namur | Efficiency of compaction and compositional convection during mafic crystal mush solidification: the Sept Iles layered intrusion[END_REF][START_REF] Namur | Crystallization of interstitial liquid and latent heat buffering in solidifying gabbros: Skaergaard intrusion[END_REF]. In these cumulates, part of the plagioclase crystals have largely unzoned cores, which are surrounded by a mantle with decreasing anor-thite content, followed by an overgrowth of constant low-anorthite composition with widths of up to several hundred micrometers. The plagioclase cores are mostly invoked to have formed at the top of a crystal mush and their mantles during burial in this mush (e.g. [START_REF] Namur | Crystallization of interstitial liquid and latent heat buffering in solidifying gabbros: Skaergaard intrusion[END_REF]. The unzoned overgrowths may reflect postcumulus growth of interstitial melt during continuous convec-tive melt exchange with a large, evolved melt reservoir [START_REF] Namur | Efficiency of compaction and compositional convection during mafic crystal mush solidification: the Sept Iles layered intrusion[END_REF] or crystallization at decreasing undercooling buffered by latent heat release [START_REF] Namur | Crystallization of interstitial liquid and latent heat buffering in solidifying gabbros: Skaergaard intrusion[END_REF].

Convective melt exchange with a large melt reservoir may be possible in the case of mafic system, but not in the case of highly viscous silicic system at high crystal fraction [START_REF] Huber | Homogenization processes in silicic magma chambers by stirring and mushification (latent heat buffering)[END_REF][START_REF] Glazner | Magmatic life at low Reynolds number[END_REF]. Significant latent heat release at near-solidus conditions may have decreased undercooling in the solidi-fication fronts of the Jiuzhou system and buffered the composition of plagioclase and alkali feldspar overgrowths in the Jiuzhou sys-tem (cf. [START_REF] Namur | Crystallization of interstitial liquid and latent heat buffering in solidifying gabbros: Skaergaard intrusion[END_REF]. The largest amount of latent heat in granitic system is released by near-solidus crystallization [START_REF] Huber | Homogenization processes in silicic magma chambers by stirring and mushification (latent heat buffering)[END_REF]. The calculated composition of the interstitial melt in the Jiuzhou system is close to a minimum composition (Table 1; Fig. 6), and we therefore posit that thermally and compositionally buffered crystallization of this melt at a ternary minimum pro-duced the compositionally unzoned overgrowths. The plagioclase and alkali feldspar overgrowths are partially developed on crystals in the charnockites (on ∼2%) and orthopyroxene-free granites (on ∼4-6%), which suggests partial (or significant) melt expulsion. The expulsion of this minimum melt upward in the crystal mush is in-ferred to have partially resorbed cumulus crystals (which typically have subhedral to subround cores) and then resumed crystalliza-tion at the ternary minimum. To our knowledge, quantitative stud-ies of mineral core-rim zoning from stratigraphically controlled granite samples that may record melt extraction are not available for comparison, but will be essential for further assessing the im-portance and for quantifying melt extraction processes in granitic magma system.

The driving forces of fractionation

Mafic igneous cumulates, such as those of the Skaergaard intru-sion, are generally layered on the hand specimen-to outcrop-scale, characterized by contrasting mineral assemblages, modes, compo-sitions and/or textures [START_REF] Holness | Assessing the role of compaction in the formation of adcumulates: a microstructural perspective[END_REF]. In the Jiuzhou plu-ton and in many granite intrusions worldwide, layering on the outcrop-scale is not developed, but a weak fabric may be evi-dent in hand specimen and thin section. For the Jiuzhou rocks, this fabric is defined by the preferred orientation of non-equant crystals and subparallel chains of plagioclase and other minerals (Fig. 2a and Fig. 3). It may have formed by magma flow [START_REF] Holness | Assessing the role of compaction in the formation of adcumulates: a microstructural perspective[END_REF], crystal settling [START_REF] Jerram | Quantifying the building blocks of igneous rocks: are clustered crystal frameworks the foundation?[END_REF], and/or crystal-mush compaction [START_REF] Philpotts | Crystal-mush compaction and the origin of pegmatitic segregation sheets in a thick flood-basalt flow in the Mesozoic Hartford basin[END_REF]. Viscous deformation by dislocation creep is locally recorded by (1) undulose extinction, subgrains, and deformation twins in plagioclase and rare bending of crystals (Fig. 2c-e); (2) rare undulose extinction and subgrains in alkali feldspar; (3) undulose extinction and chessboard subgrain patterns in quartz (Fig. 2d); and (4) bending and kinking of bi-otite (Fig. 2f). The deformation features observed in plagioclase, alkali feldspar, quartz, and biotite of the charnockites record lim-ited deformation by dislocation creep at super-solidus and/or at sub-solidus conditions, where the quartz chessboard subgrain patterns indicate viscous deformation at >600 • C [START_REF] Law | Deformation thermometry based on quartz c-axis fabrics and re-crystallization microstructures: a review[END_REF]. The relative scarcity of deformation features in weak phases such as biotite, their absence in the low-anorthite plagioclase overgrowths (Fig. 2e) and in the alkali feldspar rims and tendrils, are consis-tent with super-solidus deformation that ceased prior to interstitial melt crystallization, although minor subsolidus deformation may have affected biotite and quartz.

The contact boundaries of plagioclase crystals in crystal clus-ters and chains of the charnockites, which are characterized by long, mostly straight and at times slightly curved or serrated 2D outlines (Fig. 2a), may record heterogeneous nucleation, syn-neusis, or dissolutionreprecipitation/diffusion creep (e.g. [START_REF] Holness | The influence of interfacial energies on igneous microstructures[END_REF]. Heterogeneous nucleation should produce a random texture and not the observed preferred orientation of the rock-forming crystals. Synneusis would prefer-entially arrange crystals in a subparallel orientation along the lowenergy crystal faces (e.g. ( 010) interfaces of plagioclase, [START_REF] Vernon | A Practical Guide to Rock Microstructures[END_REF], which is not the case. Long contact lengths of crystals may also form by "welding", i.e. by dissolution-reprecipitation/diffu-sion creep under gravitational compaction [START_REF] Schmidt | Settling and com-paction of olivine in basaltic magmas: an experimental study on the time scales of cumulate formation[END_REF]. Deformation-driven dissolution-reprecipitation initially produces serrated contacts, but they rapidly anneal at high temperature (e.g. [START_REF] Mukai | Textural evolution of plagio-clase feldspar across a shear zone: implications for deformation mechanism and rock strength[END_REF]. The plagioclase and alkali-feldspar overgrowths are, moreover, asymmetrical on crystals in the charnockites, and the largest volumes of these overgrowths occur in apparent pres-sure shadow zones of the phenocrysts, which supports the interpretation that deformation-driven dissolution-reprecipitation was important (cf. [START_REF] Holness | Assessing the role of compaction in the formation of adcumulates: a microstructural perspective[END_REF].

Viscous deformation should have been negligible in the case of fabric formation by magma flow and crystal settling, but im-portant in the case of viscous compaction [START_REF] Holness | Assessing the role of compaction in the formation of adcumulates: a microstructural perspective[END_REF][START_REF] Fiedrich | Mineralogical, geochemical, and textural indicators of crystal accumulation in the Adamello Batholith (northern Italy)[END_REF]. We therefore suggest that the fabrics and the deformation of plagioclase, alkali feldspar, biotite and quartz in the charnockite are developed by crystal mush compaction with evi-dence for significant diffusion and limited dislocation creep, while a contribution of magmatic flow to the rock fabric cannot be ruled out. Below, we evaluate (1) if the magma mush in the Jiuzhou plu-ton was thick enough to self-compact; (2) if its compaction may have exerted enough strength to induce diffusion and/or disloca-tion creep.

Compaction-driven fractionation

In the absence of regional deformation, gravity-driven com-paction is commonly invoked to drive melt segregation in granitic magma system (e.g. [START_REF] Fiedrich | Mineralogical, geochemical, and textural indicators of crystal accumulation in the Adamello Batholith (northern Italy)[END_REF]. Gravity-driven com-paction, however, requires that the overlying crystal mush is thick and dense enough to provide sufficient strength, which can be as-sessed by the critical/characteristic compaction length [START_REF] Mckenzie | The generation and compaction of partially molten rock[END_REF]. This length scale depends on the rheological properties of the crystal framework and the interstitial melt. The bulk and shear viscosities of the crystal framework have not been experimentally constrained, but can be conservatively estimated using the rhe-ological model of [START_REF] Caricchi | Non-Newtonian rheology of crystal-bearing magmas and implications for magma ascent dynamics[END_REF] and [START_REF] Costa | A model for the rheology of particle-bearing suspensions and partially molten rocks[END_REF] (details are provided in supplementary file 1). For strain rates of ∼10 -12 -10 -9 s -1 (cf. [START_REF] Nicolas | Flow mechanism and viscosity in basaltic magma chambers[END_REF] Using these values, we then calculate a characteristic compaction length of 0.2-10 m for crystal mush at porosities of ≤0.1-0.5 (as detailed in section 5 of supplementary file 1), which is close to the compaction length of <1 m that has been estimated by Lee et al. (2015) for an evolved granite intrusion.

Individual outcrops in the Jiuzhou pluton expose thicknesses of up to ≥100 m without significant layering or obvious internal contacts, which could represent individual sill or crystal mush layer (see details in section 5 in supplementary file 1). On the basis of our calculations, field ob-servations, and the systematic, vertical compositional pluton-scale zoning (Fig. 6 and Table S1), we thus conclude that compaction of the Jiuzhou crystal mush was physically feasible.

Evidence for viscous deformation (by diffusion and dislocation creep) increases from the stratigraphically lower, orthopyroxene-bearing units (charnockite sample JZ01, Fig. 2a and3a), over the middle units (orthopyroxene-free granite sample JZ10, Fig. 3b), to the higher units (orthopyroxene-free granite sample JZ06, Fig. 2b and3c). The CSD patterns determined for plagioclase are con-sistent with increasing crystal accumulation towards the lower stratigraphic levels of the pluton (Fig. 4). Two of our charnock-ite samples have concave CSD patterns, where curves are deflected upward for grain sizes of >∼4 mm (Fig. 4). Such grain size distri-butions and CSD curves may reflect the modification of initially linear CSDs by crystal accumulation [START_REF] Higgins | Origin of megacrysts in granitoids by textural coarsening: a crystal size distribution (CSD) study of microcline in the Cathedral Peak Gran-odiorite, Sierra Nevada, California[END_REF]. The small number of grains <1mm in size may reflect pressure dissolution by compaction [START_REF] Schmidt | Settling and com-paction of olivine in basaltic magmas: an experimental study on the time scales of cumulate formation[END_REF], or limited nucleation of small crystals in the already crystal-rich cumulates (Holness and Ver-non, 2015). Our most evolved charnockite sample (JZ09) and the orthopyroxene-free granites, in contrast, have approximately linear CSDs (Fig. 4), which suggests that the accumulation of crystals was less significant [START_REF] Higgins | Origin of megacrysts in granitoids by textural coarsening: a crystal size distribution (CSD) study of microcline in the Cathedral Peak Gran-odiorite, Sierra Nevada, California[END_REF].

The strength required for deformation by diffusion and dis-location creep is estimated to have been ∼0.01-40 kPa and ∼200-5000 kPa, respectively, for strain rates of ∼10 -12 -10 -9 s -1 (cf. [START_REF] Kirby | Rheology of the lithosphere: selected topics[END_REF][START_REF] Nicolas | Flow mechanism and viscosity in basaltic magma chambers[END_REF]; see calculations in supplementary file 1). For an estimated gravita-tional pressure gradient of ∼3 kPa/m (cf. [START_REF] Schmidt | Settling and com-paction of olivine in basaltic magmas: an experimental study on the time scales of cumulate formation[END_REF], a crystal mush layer >0.01-10 m thick exerts then enough strength to induce diffusion creep, while a crystal mush layer >70-1000 m thick is required to induce dislocation creep. Dislocation creep would thus have operated in a thicker mush layer and at a later stage of crystallization and compaction than diffusion creep, i.e. significant deformation by diffusion creep would have been fol-lowed by minor dislocation creep.

The undeformed, unzoned plagioclase and alkali feldspar over-growths are interpreted to have crystallized from trapped melt subsequent to mush compaction. In the charnockites and ortho-pyroxene-free granites, these overgrowths comprise ∼6 vol% and ∼18-22 vol%, respectively. A significant proportion of quartz (higher than that of plagioclase and alkali feldspar, cf. [START_REF] Wilke | Experimental investigation of the effect of Ca, Fe and Ti on cotectic compositions of the rhyolitic system[END_REF][START_REF] Wilke | The effect of anorthite content and water on quartz-feldspar cotectic compositions in the rhyolitic system and implications for geobarometry[END_REF] must have also crystallized from the near-solidus interstitial melt, but the distribution of this late-stage quartz can-not be estimated on the basis of our petrographic observations and imaging. We can nevertheless roughly estimate the fraction of the quartz based on the modal proportion of low-anorthite and alkali feldspar overgrowths in the Qz-Ab-Or-An model system [START_REF] Wilke | Experimental investigation of the effect of Ca, Fe and Ti on cotectic compositions of the rhyolitic system[END_REF], i.e. ∼4 vol% for the charnockite and ∼11-13 vol% for the orthopyroxene-free granite. We thus suggest that the proportion of trapped melt is significantly lower in the charnockites (∼10 vol%) relative to the orthopyroxene-free granites (∼29-35 vol%), which is consistent with their stratigraphically lower position and more primitive composition (Table S1 in supplementary file 1).

The recorded amount of viscous deformation by dislocation creep in the Jiuzhou plutonand in many other granite plu-tons worldwideis significantly weaker than in many mafic layered series, e.g. in the Skaergaard intrusions (where disloca-tion creep may not have occurred at super-solidus conditions and for which viscous compaction remains controversial, [START_REF] Holness | Assessing the role of compaction in the formation of adcumulates: a microstructural perspective[END_REF]. However, as we have detailed above, deformation by diffusion creep (dissolution-reprecipitation) may be more impor-tant in silicic systems than deformation by dislocation creep. Ev-idence for diffusion creep may be elusive, as obvious markers for this process (e.g. truncated oscillatory zoning patterns indicat-ing partial dissolution) are rare. Further systematic and detailed microstructural studies for granite plutons are clearly needed to critically assess if (1) the recorded viscous deformation is caused by gravitationally-driven viscous compaction or other subsolidus processes; and

(2) how important viscous compaction is as a fractionation driving force in general (e.g. [START_REF] Fiedrich | Mineralogical, geochemical, and textural indicators of crystal accumulation in the Adamello Batholith (northern Italy)[END_REF][START_REF] Holness | Melt segregation from silicic crystal mushes: a critical ap-praisal of possible mechanisms and their microstructural record[END_REF].

Thermodynamic modeling for cumulate systems

Equilibrium thermodynamic modeling for partial disequilibrium cumulate systems is per se problematic and an oversimplifica-tion. For example, solidus plagioclase is predicted by thermody-namic modeling to crystallize with ∼An48 composition from a magma with the composition of charnockite sample JZ01-1, and with ∼An40 composition for a magma with the composition of the orthopyroxene-free granite sample JZ10 at system H2 O con-tents of 1 and 2 wt% (see details of the modeling in supplemen-tary file 1). These modeled compositions of solidus plagioclase are lower than the average compositions of the analyzed plagio-clase crystals (excluding the mass and composition of the low-anorthite overgrowths), i.e. ∼An50 and ∼An47 in samples JZ01-1 and JZ10, respectively (Fig. 5c,f). If the cumulus phases of the magmas reached total equilibrium, then the thermodynamic mod-eling using sample JZ01-1 and JZ10 compositions should predict plagioclase compositions equivalent to those characteristic for the natural samples. If ∼10 and ∼30 wt% of an evolved, interstitial melt with an invariant minimum composition (cf. [START_REF] Nishimura | In situ fractional crystallization observed in the Os-umi granodiorite batholith[END_REF], i.e. with ∼75 wt% SiO2 (Table 1), is subtracted from the bulk composition of samples JZ01-1 and JZ10, then plagioclase with ∼An50 and ∼An47 composition is predicted to crystallize at the solidus, matching the average composition of the natural crys-tals (Fig. 5c). The fractions of trapped melt are thus inferred to have been ∼10±5 wt% (∼12±5% in volume) and ∼30±5 wt% (∼33±5% in volume) in the most mafic charnockite (JZ01-1) and in the most mafic orthopyroxenefree granite (JZ10) samples, re-spectively. These estimates are consistent with the estimated melt fractions based on our petrological observations (see section 5.1.4). They are lower than the porosities of loosely packed cumulates (∼50-60 vol%, [START_REF] Schmidt | Settling and com-paction of olivine in basaltic magmas: an experimental study on the time scales of cumulate formation[END_REF], but are consistent with cu- mulate porosities that can be achieved by predominantly mechan-ical compaction in the case of the orthopyroxene-free granites (cf. [START_REF] Higgins | The origin of laminated and massive anorthosite, Sept Iles lay-ered intrusion, Québec, Canada[END_REF].

The inferred trapped melt fraction of ∼10±5 wt% in the charnockites is consistent with mechanical compaction fol-lowed by significant viscous compaction (cf. [START_REF] Schmidt | Settling and com-paction of olivine in basaltic magmas: an experimental study on the time scales of cumulate formation[END_REF].

The fraction of melt lost from the solidification fronts during compaction can be further estimated to have been ∼40 wt% and ∼20 wt% for the charnockite and the orthopyroxene-free magmas, respectively, if we assume that melt extraction took place at a crys-tallinity of ∼50 wt%. The composition of the parental magma (i.e. the magma from which the cumulus minerals crystallized) of the most primitive solidification front represented by the most mafic charnockite sample (JZ01-1) can thus be reconstructed to have had ∼68 wt% SiO2 (Table 1), which is close to the average compo-sition of all our samples, while the composition of the parental magma for the most mafic orthopyroxene-free sample (JZ10) is more evolved with ∼70 wt% SiO2

(Table 1). Thermodynamic mod-eling using these reconstructed parental magma compositions will derive the most realistic constraints for opensystem fractiona-tion, while the uncertainties of thermodynamic modeling using cumulate compositions may be significant (Fig. 7; see details in section 5.3 in the supplementary file 1). We posit that thermody-namic modeling predicts comparable phase relations for relatively dry granitic system with moderately evolved composition at near solidus condition within ≤30-50 • C for fractionated systems with ≤30-45 wt% melt loss compared with the model results of the parental magma (with 68-70 wt% SiO2 ). The modeling, however, overestimates the phase boundaries at near liquidus condition and thus the liquidus temperature by 50-80 • C. These uncertainties add to those inherent in the modeling approach (∼20-60 • C; Zhao et al., 2017a). The combined uncertainty may be small enough to be within the uncertainty of independent estimates (e.g. from con-ventional mineral thermometers) as in the case of the previously modeled phase relations for the Jiuzhou charnockite (Zhao et al., 2017a). The combined uncertainties may, however, also become significant (≥100 • C), limiting a meaningful interpretation. Dise-quilibrium features, such as the step-zoned plagioclase crystals are thus important to be considered, even if their proportion appears relatively minor (as in our charnockite sample JZ01-1), and a re-construction of the parental magma composition is required.

A conceptual model for fractionation and cumulate formation in silicic magma systems

Fractionation of silicic magma systems is commonly consid-ered to involve the segregation of interstitial melt from a ho-mogeneous crystal mush (e.g. [START_REF] Bachmann | On the origin of crystal-poor rhyolites: ex-tracted from batholithic crystal mushes[END_REF][START_REF] Gelman | Identifying the crystal graveyards remaining after large silicic eruptions[END_REF]. The solid-rich fractionates of the Jiuzhou pluton are, however, not homogeneous, ranging from charnockite to orthopyroxene-free granite, advocating incremental fractiona-tion (e.g. as in the concept of [START_REF] Langmuir | Geochemical consequences of in situ crystallization[END_REF]. Below, we build a conceptual model (Fig. 8a-g) for the magmatic evolution of the Ji-uzhou pluton using the most mafic charnockite sample JZ01-1 and the orthopyroxene-free granite sample JZ10 as end-members.

Thermodynamic modeling using the reconstructed parental magma composition predicts that plagioclase, orthopyroxene, quartz, alkali feldspar and biotite saturate at >800-790 • C and at crystallinities of >40-50 wt% before melt extraction (Fig. 7a). The orthopyroxene-out boundary is predicted at ∼780 • C and at a crystallinity of ∼55 wt%; orthopyroxene is thus predicted to have been stable prior and subsequent to melt extraction (Fig. 7a). Melt extraction increases the survival chances of orthopyroxene, as this process leaves a relatively small proportion of melt for reaction with orthopyroxene. Many other granite plutons worldwide, specif-ically many garnet-and/or cordierite-bearing granites, may equally record such a fractionation process, as garnet and cordierite are commonly non-solidus phases similar with orthopyroxene in the Jiuzhou magmas.

In the more evolved orthopyroxene-free granites of the Ji-uzhou pluton, plagioclase, cordierite, quartz, and biotite saturate at >770-750 • C and 40-50 wt% crystallinities, and thus before melt extraction (Fig. 7c). Orthopyroxene should have formed initially, but is totally resorbed at ∼780 • C and at ∼40 wt% crystallinity (Fig. 7c). For magma more evolved than the orthopyroxene-free sample JZ10, the orthopyroxene stability field shrinks, i.e. the orthopyroxene-out boundary shifts from a crystallinity of ∼40 wt% for a magma with ∼70 wt% SiO2 to a crystallinity of ∼20 wt% for a magma with ∼72 wt% SiO2 (Fig. 9). The change of orthopy-roxene stability field with magma composition was, in addition to significant melt extraction, a key factor that permitted the local preservation of orthopyroxene in the marginal, relatively primitive solidification fronts of the Jiuzhou intrusion (Fig. 8g), yet not in the more evolved interior of the intrusion.

Conclusions and implications

Many volcanic and plutonic magma systems with intermediate to silicic compositions are invoked to have fractionated by crystal-liquid segregation, but microstructural evidence for this process is rarely reported. Our detailed microstructural observations re-vealed the crystal accumulation and melt extraction in the intru-sive magma body, and further combination of petrographic study with thermodynamic modeling quantitatively constrained the pro-portions of cumulate phases, trapped melt fractions and the composition of parental magma(s). For the studied Jiuzhou pluton, gravitational compaction has reduced the trapped melt fraction in crystal mush from ∼30 wt% at the upper level to ∼10 wt% at the lower level of the pluton, which results in the deviation of the whole-rock compositions from the metasediment-sourced experimental melt compositions. The studied rocks show continu-ous compositional zoning consistent with incremental fractionation from parental magmas with ∼68 wt% SiO2 composition. Using the reconstructed parental magma as compared to whole-rock compo-sitions is crucial for thermodynamic modeling of phase relations and intensive parameters, which otherwise produce erroneous es-timates. Many granite plutons worldwide also have compositions deviating from those of primary melts, which may have experi-enced fractionation driven by gravitational compaction or by other processes such as hindered settling and gas-driven filter pressing. Detailed microstructural studies combining thermodynamic mod-eling for such systems will permit a critical appraisal of their frac-tionation processes with implications for understanding plutonic-volcanic connections, formation of mineral deposits and crustal differentiation.

Fig. 1 .

 1 Fig. 1. (a) Tectonic schematic map of South China and its surrounding plates (Zhao et al., 2017b and references therein). (b) Geological map showing the area in which the orthopyroxene-bearing to orthopyroxene-free Jiuzhou pluton is exposed. Filled and open stars within the Jiuzhou pluton mark our sample locations for charnockite and orthopyroxene-free granite, respectively.

Fig. 2 .

 2 Fig. 2. Microphotographs of the charnockites and the orthopyroxene-free granites from the Jiuzhou pluton, South China, in cold-cathodoluminescence (CL) imaging mode and under crossed polarized light. (a, b) Cold-CL images showing plagioclase clusters in the charnockite (a) and in the orthopyroxene-free granite (b). The low-anorthite overgrowths (dark green) on anorthite-rich cores (bright green) in the charnockite are developed on some of the external plagioclase grain boundaries, but they are absent from plagioclase-plagioclase grain boundaries within the cluster. The overgrowths on plagioclase in the orthopyroxene-free granite, however, form complete, albeit asym-metrical rims. The white dashed lines denote the boundaries between partial overgrowths and cores. The red solid lines indicate the contact boundaries between cores and cores. (c) Tapering deformation twins in plagioclase (marked by arrows), which are slightly bent at an angle of ∼2 • . (d) Deformation-twinned plagioclase crystal with un-dulose extinction, and quartz crystal with chessboard subgrain pattern and undulose extinction (marked by arrows). (e) Plagioclase with a low-anorthite overgrowth. Note that the deformation twins taper at the boundary of the overgrowth and core. (f) Euhedral biotite crystal between plagioclase crystals, which is significantly kinked. Mineral abbreviations follow Whitney and Evans (2010). (For interpretation of the colors in the figure(s), the reader is referred to the web version of this article.)

Fig. 3 .

 3 Fig. 3. Digitized images based on cold-CL, thin-section scale maps showing the mineral assemblages and core-rim step-zoning in plagioclase and alkali feldspar (a) for charnockite and (b and c) for orthopyroxene-free granites. (a) The charnockite comprises highly clustered plagioclase, which link to form chains (see the spatial distribution pattern of plagioclase in the bottom-right inset). The black (and red) dashed lines mark the transition from anorthite-rich plagioclase (and orthoclase-poor alkali feldspar) cores to their anorthite-poor (and orthoclase-rich) overgrowths. (b and c) The orthopyroxene-free granite is characterized by a random to slightly clustered distribution of plagioclase and the other mineral types. The corresponding cross-polarized photographs and cold-CL images are provided in supplementary file 1 (Figs. S1 and S2). Acc in the legend represent accessories or epoxy.

Fig. 4 .

 4 Fig. 4. Crystal size distribution (CSD) curves for plagioclase from the charnockite and the orthopyroxene-free granite of the Jiuzhou pluton. The curves are labeled with their sample numbers.

Fig. 5 .

 5 Fig. 5. Backscattered electron images showing plagioclase crystals with step-zoning in (a) charnockite sample JZ01-1 and (b) orthopyroxene-free granite sample JZ10. (c, d) Compositional profiles indicated in (a) and (b). (e and f) Histograms showing the anorthite content of anorthite-rich plagioclase cores and anorthite-poor overgrowths. The solid and dashed arrows denote the average composition of the natural crystals excluding the composition of the overgrowths and the solidus plagioclase composition predicted by thermodynamic modeling (detailed in Fig. S2 in supplementary file 1).

Fig. 6 .

 6 Fig. 6. (a-f) Harker diagrams showing the composition of the Jiuzhou rocks and the modeled compositional evolution for (1) restite entrainment (solid black lines with triangles), (2) incremental fractionation (solid blue lines with diamonds). The solid and open red stars represent whole-rock compositional data from Zhao et al. (2017b); the gray circles represent whole-rock compositional data from Qi et al. (2007), Zhao et al. (2012), and Jiao et al. (2015). The open black stars denote the average composition of our samples. The dashed black lines with green squares indicate the evolution of the interstitial melt. The yellow fields highlight the compositional range of experimental liquids produced by partial melting of metasediments (see compiled dataset in supplementary file 2).

  at a porosity of ∼0.1. At higher porosi-ties, the shear and the bulk viscosities of the crystal framework decrease. The shear viscosity of the interstitial melt with ∼75 wt% SiO2 and ∼3.0 wt% H2 O at ∼800 • C and with ∼5.6 wt% H2 O at ∼750 • C can be estimated to be ∼10 5.4 -10 6.2 Pa s using the model of Giordano et al. (2008) (see details in supplementary file 1).

Fig. 7 .

 7 Fig. 7. Modeled P-T phase diagrams and modal abundance of crystals formed at 0.2 GPa (fraction in wt%) for the reconstructed parental magma compositions (Ta-ble 2) of sample JZ01-1 (a, b) and sample JZ10 (c, d). The dashed lines in (c) and (d) show the modeling resultsfor sample JZ01-1 for an initial melt content of 2 wt% and for sample JZ10 for an initial 2.7 wt% melt H 2 O content (see detailed model-ing results in supplementary file 1). Solid lines (a, and b) mark the phase stability fields and dashed lines labeled with numbers denote the plagioclase anorthite iso-pleths (in mol%). Phases are labeled inside their stability field. Solid and dotted blue lines with squares in (b, and d) denote the modeled crystallization sequence using the constructed parental magma compositions and bulk compositions of samples as a function of temperature and crystallinity. Note that the modeled alkali feldspar-, quartz-and cordierite-in boundaries for both compositions are comparable within 30-50• C, while the temperatures of biotite-in, orthopyroxene-out and cordierite-out boundaries do not change. Solid and dotted red lines in (b, d, f and h) denote the melt fraction as a function of temperature.

Fig. 8 .

 8 Fig. 8. Schematic cartoon illustrating incremental fractionation (in situ fractional crystallization) proposed for the Jiuzhou pluton. (a-c) Cartoons showing the proposed crys-tallization history with decreasing temperature for the Jiuzhou charnockite magmas. (d-f) Cartoons showing the proposed crystallization history for the orthopyroxene-free granite magmas. (g) Conceptual snapshot summarizing the inferred zonation and paths of melt migration at crystallinities of ∼40-50 wt%. The dashed green arrows denote the direction of interstitial melt expulsion. The numbers in the legend indicate that most minerals are texturally and/or chemicall bimodal: 1 and 2 for orthopyroxene repre-sent euhedral and resorbed anhedral crystals; 1 and 2 for biotite represent anhedral biotite surrounding orthopyroxene and euhedral biotite; 1 and 2 for plagioclase represent crystals with weak zoning or step-zoned crystals; 1 and 2 for quartz and alkali feldspar represent cumulus and postcumulus crystals (see details in text).

Fig. 9 .

 9 Fig.9. The crystallinity at the orthopyroxene-out boundary as a function of evolving magma composition in the fractionating magma chamber in our model system (de-tails are reported in supplementary file 3). The gray area highlights the extraction window, i.e. the crystallinities for which melt extraction is inferred. Orthopyroxene was preserved only in solidification zones in which it formed early at crystallinities of ≥40%.

Table 1

 1 Whole-rock compositions and calculated model compositions used in our numerical modeling.

		JZ01-1	a	JZ10	b	Average c	Minimum melt d	Minimum melt e	Parental magma	f for JZ01-1	Parental magma	g for JZ10
	SiO2	62.95		68.63	67.72	74.74	75.73	68.25		70.48
	TiO2	1.11		0.74		0.80	-	0.28	0.63		0.54
	Al2O3	15.63		14.05	14.48	14.70	12.96	13.86		13.34
	Fe2O3*	7.06		5.15		5.40	1.91	1.71	4.56		4.04
	MnO	0.09		0.07		0.07	-	0.03	0.07		0.06
	MgO	2.16		1.64		1.80	0.31	0.31	1.30		1.23
	CaO	3.64		2.46		2.51	0.49	1.12	2.21		1.86
	Na2O	2.29		2.23		2.09	1.93	3.11	2.23		2.21
	K2O	3.43		3.56		3.63	6.52	4.85	4.52		4.24
	H2O	-		-		-	3.5	-	2.0		2.7
	Note: crystallinity of 50%. a-b , Data are after Zhao et al. (2017b); e , minimum melt composition reported by c , the average composition of our samples.	d , minimum melt composition estimated from thermodynamic modeling using the average composition at a

Acknowledgements

We would like to sincerely thank Editor Tamsin Mather, re-viewer Chris Huber, and an anonymous reviewer for their many constructive comments and suggestions, especially those related to the role of deformation and compaction. We would also like to thank Guan C.G. and Zhou W.M. for access to their facility and support with the cathodoluminescence imaging. This work was fi-nancially supported by the Natural Science Foundation of China ( 41430208) and the program B for Outstanding PhD candidates at Nanjing University (201702B071).

Appendix A. Supplementary material

Supplementary material related to this article can be found on-line at https://doi.org/10.1016/j.epsl.2018.07.039.