Whiteschist genesis through metasomatism and metamorphism in the Monte Rosa nappe (Western Alps)

Whiteschists from the Monte Rosa Nappe were examined in the field as well as with petrographic, geochemical, and isotopic methods to constrain the controversial origin of these rocks in their Alpine metamorphic context. Whiteschists occur as ellipsoidal-shaped, decametric-sized bodies, within a Permian metagranite, and consist mainly of chloritoid, talc, phengite, and quartz. The transition from whiteschist to metagranite is marked by multiple sharp mineralogical boundaries defining concentric zones unrelated to Alpine deformation. The development of reaction zones, as well as the geometry of the whiteschist suggest a pervasive fluid infiltration, facilitated and canalized by reaction fingering. Whole-rock compositions of whiteschists and metagranites indicate an enrichment in MgO and H2O and depletion of Na2O, CaO, Ba, Sr, Pb, and Zn in the whiteschist relative to the metagranite. Trace- and rare-earth elements, together with all other major elements, notably K2O and SiO2, were within uncertainty not mobile. Such a K and Si saturated, Na undersaturated fluid is not compatible with previous interpretations of fluids derived from ultramafic rocks, evaporites, or Mg-enriched seawater due to mantle interactions. Together with the large variations in δD and δ18O values, this indicates large fluid fluxes during metasomatism. Calculated δD and δ18O values of fluids in equilibrium with the whiteschist support a magmatic–hydrothermal fluid source, as does the chemical alteration pattern. Bulk rock 87Sr/86Sr ratios in whiteschists confirm a pre-Alpine age of fluid infiltration. The 87Sr/86Sr ratios in whiteschists were subsequently partially homogenized in a closed system during Alpine metamorphism. In conclusion, the granite was locally affected by late magmatic–hydrothermal alteration, which resulted in sericite–chlorite alteration zones in the granite. The entire nappe underwent high-pressure metamorphism during the Alpine orogeny and the mineralogy of the whiteschist was produced during dehydration of the metasomatic assemblage under otherwise closed-system metamorphism. While each whiteschist locality needs to be studied in detail, this in-depth study suggests that many whiteschists found in granitic bodies in the Alps might be of similar origin.


Introduction
Whiteschists have first been described as talc and kyanite schists by Schreyer (1973) after the discovery of a number of occurrences of this special magnesium-aluminium-rich paragenesis in various geological settings over the world (Schreyer 1974). They are commonly found in continental units that underwent medium-to-high-pressure metamorphism and, depending on the metamorphic grade, various Mg-Al silicates are reported (Schreyer and Seifert 1969). The chemical composition of the rocks is relatively simple (Franz et al. 2013): K 2 O-FeO-MgO-Al 2 O 3 -SiO 2 -H 2 O (KFMASH). The rarity of these rocks and their simple composition has attracted the curiosity of petrologists, since there is no common protolith for this unusual rock.
Communicated by Steven Reddy.
Numerous protoliths have been proposed, but these are still debated.
In the Western Alps, whiteschists are found in the three Internal Crystalline Massifs: Monte Rosa, Gran Paradiso, and Dora Maira. Bearth (1952) and later Dal Piaz and Lombardo (1986) interpreted the whiteschists from the Monte Rosa as products of a metasomatic alteration of granites based on field evidence. Alternatively, a protolith consisting of meta-evaporites or Mg-rich pelites has been proposed by Chopin (1981) and Schreyer (1977). These two propositions have since been supplemented to include rocks of metasedimentary or metavolcanic origin in crustal domains (Chopin 1984;Schertl et al. 1991). The discrete and isolated character of the whiteschist occurrences within the Dora Maira orthogneisses later led Chopin et al. (1991) to favor Bearth's (1952) original interpretation though. The metasomatic origin of the Dora Maira whiteschist protolith was also supported by subsequent workers (Schertl et al. 1991;Gebauer et al. 1997;Compagnoni and Hirajima 2001). A metasomatized granite origin of the whiteschists in the Monte Rosa nappe was argued for by Pawlig and Baumgartner (2001).
The timing of alteration, as well as the source of the fluids is also a matter of debate. Ferrando (2012) argued that the whiteschists found in the Internal Crystalline Massifs of the Western Alps have many characteristics in common and, hence, have a common origin. Ferrando (2012) proposed that they are the result of Mg metasomatism, related to dehydration of ultramafic rocks, which interacted with ocean water. Based on geochemical and fluid-inclusion studies, she proposed that fluid-rock interaction is localized along highstrain zones. Different tectonic environments were suggested for the metasomatic event, including continental-oceanic transition zones, subduction, and continent-continent collision zones. Based on a stable isotope study, Sharp et al. (1993) proposed a metasomatic fluid originating from the dehydration of altered oceanic crust at great depth for the Dora Maira slice during the Alpine orogeny. In contrast, Pawlig and Baumgartner (2001) used stable isotope and whole-rock compositions to suggest that the alteration is the result of late-stage magmatic-hydrothermal fluid-rock interaction, which occurred during the Permian, thus prior to Alpine metamorphism. In a recent study, Chen et al. (2019) suggest that the metasomatic fluids of the Dora Maira whiteschist were reduced fluids, emphasizing the role of oxygen fugacity on the formation of whiteschist.
The present study contributes to the debate on the origin of whiteschists by presenting detailed field and geochemical data to clarify the nature and timing of the metasomatic alteration. The southwestern edge of the Monte Rosa nappe has been the focus of studies on whiteschists since Bearth's seminal work (Bearth 1952;Chopin and Monié 1984;Pawlig and Baumgartner 2001;Darbellay 2005). More recently, the whiteschists have been at the centre of highly controversial idea of tectonic over-pressure . A crucial part not addressed in the above-mentioned short paper is the origin of the chemical composition of the whiteschist, discussed in detail here. In the Val d'Ayas area, two whiteschists occur as 10-50 m large bodies within the Monte Rosa metagranite, and two others are located at the contact between the metagranite and the surrounding basement paragneisses. The two whiteschists occurring within the metagranite have been less affected by the late Alpine, retrograde overprint and maintained their Alpine peak metamorphic paragenesis of talc-chloritoid-phengite-quartz. Field and petrographic observations are complemented by an extensive whole rock and stable isotope geochemical study, to help identify the fluid source responsible for the metasomatic alteration leading to the whiteschist formation. Finally, rubidium-strontium whole-rock data will complement the field study to establish the timing of the fluid event.

Geological setting
The Monte Rosa nappe is a basement nappe belonging to the middle Penninic unit of the Western Alps (Bearth 1952). It forms a large backfold, structurally enclosed by the ultramafic and mafic rocks of the Zermatt-Saas unit above and the Antrona unit below (Steck et al. 2015). From a paleogeographic point of view, the Monte Rosa nappe represents the southernmost tip of the Briançonnais microcontinent that belongs to the European continental margin (Dal Piaz 2001;Steck et al. 2015).
The nappe consists of pre-Alpine polymetamorphic paragneisses ( Fig. 1), locally containing small volumes of amphibolite boudins and carbonate layers (Bearth 1952;Dal Piaz 1971), which underwent high-grade amphibolite facies regional metamorphism during the Variscan orogeny at conditions of 0.5-0.7 GPa and temperatures in excess of 650 °C, locally resulting in the partial melting of the basement rocks (Bearth 1952;Dal Piaz 1971;Ferrando et al. 2002). The basement was intruded by Permian granites, dated at 269 ± 4 Ma by SHRIMP U-Pb on zircon (Pawlig 2001). Chopin and Monié (1984) were the first to investigate the Alpine peak metamorphism of the Monte Rosa nappe using a whiteschist sample from the Upper Ayas valley (Chopin and Schreyer 1983;Chopin and Monié 1984). They suggested temperatures of 500 °C and pressures of 1.6 GPa. Later, Le Bayon et al. (2006) obtained peak conditions of 2.4 GPa at 500 °C. Luisier et al. (2019) present phase petrology arguments to suggest that the peak temperature and pressure in the whiteschists were 540-600 °C and 2.2 GPa, due to potential deviation from the lithostatic pressure, in a metagranite which records maximum pressures of 1.4 GPa. The age of the Alpine peak metamorphism was constrained by U-Pb on rutile at 42.6 ± 0.6 Ma (Lapen et al. 2007) and followed by a greenschist-facies retrogression dated at 41.0 to 38.5 Ma by Rb-Sr on metasediments from the Schistes Lustrés (Skora et al. 2015). Fission track zircon and apatite data from the metagranite around the whiteschists are 34 Ma and 14 Ma, respectively (Hurford et al. 1991), suggesting that this part of the Penninic domain was quickly exhumed after its peak metamorphism. The exhumation of the Monte Rosa nappe through extrusion and backfolding is related to major deformation structures that heterogeneously overprinted both the Alpine peak metamorphism in the metagranite and the polymetamorphic history of the paragneiss basement (Steck et al. 2015).

Analytical methods
Fresh samples were crushed using a hydraulic press and powdered in a tungsten carbide mill. Powders were first dried in the oven at 110 °C and calcined at 1050 °C for 3 h. LOI was calculated by weight difference before and after calcination. Calcined powders were mixed with Li 2 B 4 O 7 and melted to produce fused glass disks. Major elements were determined by X-ray fluorescence (XRF), using a PANalytical Axios-mAX spectrometer at the University of Lausanne, Switzerland. Ferrous iron was determined by a colorimetric method with a PerkinElmer Lambda 25 spectrophotometer at the University of Lausanne.
Glass disks were used for trace-element determination by LA-ICP-MS. A NewWave UP-193 ArF ablation system, combined with an Element XR at the University of Lausanne was employed. NIST SRM 612, was used as standard and Al 2 O 3 XRF compositions served as an internal standard.
The oxygen isotope composition of whole rocks was determined at the University of Lausanne using a CO 2 laser fluorination technique with F 2 as reagent (Lacroix and Vennemann 2015) coupled to a Thermo Finnigan MAT 253 gas source mass spectrometer. Whole-rock data are corrected to the session value of an in-house quartz standard, i.e., LS-1 Qtz with an accepted value of 18.1‰. The repeatability of the LS-1 determined as average of daily replicates is usually better than 0.2‰ (1 s). Data are given in conventional δ-notation, relative to Vienna Standard Mean Ocean Water (VSMOW).
Measurements of the hydrogen isotope compositions of minerals were made using high-temperature (1450 °C) reduction methods with He-carrier gas and a TC-EA linked to a Delta Plus XL mass spectrometer from Thermo-Finnigan on 2-4 mg-sized samples according to a method adapted after Bauer and Vennemann (2014). The precision of the in-house kaolinite and G1 biotite standards for hydrogen isotope analyses was better than ±2‰ for the method used; all values were normalized using a value of − 125‰ for the kaolinite standard and -65‰ for NBS-30. The data are given in the standard δ-notation, expressed relative to VSMOW in permil (‰).
Six samples from the Mezzalama whiteschist containing carbonates have been selected for measurements of the δ 13 C and δ 18 O values. All samples come from Zone C of the transition zone described below. One is a carbonate vein (14MR12) and the five other samples contain carbonate as a matrix phase. The carbonate powder from sample 14MR12 was obtained by micro-drilling and bulk rock powders were used for the other five samples. The C-and O-isotope composition of the carbonates were measured with a GasBench II connected to a Finnigan MAT DeltaPlus XL mass spectrometer, according to a method adapted after Spötl and Vennemann (2003). In-house standards of calcite are treated in the same way and run interspersed with the samples in the same sequence. Samples are normalized using an inhouse standard calibrated against δ 13 C and δ 18 O values of NBS-19 (+1.95 and -2.20‰, relative to VPDB). External reproducibility for the analyses estimated from replicate analyses of the in-house standard reacted at 70 °C (n=5) was ±0.04‰ for δ 13 C (expressed vs VPDB) and ±0.06‰ for δ 18 O (expressed vs VSMOW). 87 Sr/ 86 Sr ratios on bulk rock powders were measured at the Department of Mineralogy of the University of Geneva, Switzerland, after the method described in detail in Béguelin et al. (2015), using a Thermo Neptune PLUS multicollector ICP-MS. 87 Steck et al. 2015). The study area is outlined by the rectangle and whiteschist outcrops are marked by stars the mass of Sr as 87,6148999 g/mol, and a 86 Sr/ 88 Sr ratio of 0.1194 and the 84 Sr/ 88 Sr ratio of 0.006756. Strontium evolution diagram was calculated using lambda of 1.404 × 10 −11 .

Field observations and petrography
The study area is located in the Upper Ayas Valley, in the Italian Alps (Fig. 1). Here, the hinge of the Monte Rosa backfold crops out and is surrounded by the two flanks of the overlying Zermatt-Saas serpentinite and metabasic rocks. The western extent of the basement paragneiss is restricted and affected by the late deformation related to greenschistfacies metamorphism (Darbellay 2005;Steck et al. 2015) and the proximity to the thrust contacts between the Monte Rosa nappe and the Zermatt-Saas unit. The paragneiss in the NW of the Mezzalama refuge has a variably strong foliation related to the Alpine deformation, overprinting the pre-Alpine metamorphic schistosity. The intensity of the Alpine deformation is decreasing from the thrust contact in the NW towards the SE, over a distance of a few tens of meters. The granite intrusive contact on the NW of the Mezzalama refuge is affected by the deformation as well, as shown by aplite and pegmatite dykes in the paragneiss that are bent into the Alpine foliation. On the other side, the SE intrusive contact (visible along the Perazzispétz ridge) is well preserved and relatively undeformed. A more-detailed structural study can be found in Darbellay (2005) and Vaughan-Hammon et al. (in review). The metagranite has large regions of undeformed to slightly foliated metagranite, separated by high strain, anastomosing ductile shear zones. The undeformed metagranite is a K-feldspar porphyritic biotite metagranite, locally containing granodiorite enclaves or paragneiss xenoliths. Alignment of the K-feldspars is visible in the least deformed parts and is interpreted to be of magmatic origin. The three whiteschist outcrops are described in detail in the following sections and identified by arrows in Fig. 2.

Mezzalama whiteschist
The best-preserved whiteschist is the Mezzalama whiteschist, located to the N of the Mezzalama refuge (624991/084967). It has not been described in detail previously, but is the location of the study of Luisier et al. (2019) on pressure variations. It crops out as an ellipsoidal-shaped body of 40 × 20 m, elongated in the SW-NE direction. In this area, the whiteschist-to-metagranite transition occurs over approximately 2-6 m and we define this zone as transition zone in this study.
The deformation in this area is given by a weak, NW dipping schistosity (Darbellay 2005), which is consistent with the main foliation in the area. The foliation is related to the nappe exhumation and back-folding (Steck et al. 2015), associated with a greenschist-facies metamorphic overprint.
The foliation is weak in the porphyritic metagranite at the contact with the whiteschist as well as in the transition zone between the whiteschist and metagranite. The centre of the whiteschist remained largely unaffected by any deformation. The orientation of the whiteschist-to-metagranite transition zone is almost vertical on the NW and NE sides of the whiteschist. It dips steeply towards SW on the SW side and dips towards the SE on the SE side of the whiteschist. The NW-oriented foliation overprints and crosscuts the whiteschist-metagranite transition. This deformation event is younger than the metasomatic alteration and the whiteschist crystallisation. There is no shear zone or fault associated with this body (Fig. 3a, b). Moreover, straight pegmatite and aplite dykes crosscut the metagranite-whiteschist transition, without being deformed. Therefore, the whiteschist is a body, which is formed in-situ, surrounded by metasomatic alteration zones; it is not a xenolith, nor a tectonic slice or body mixed into the metagranite by any Alpine structural feature.
The transition between the whiteschist and the metagranite is characterized by two reaction fronts, delimited by sharp mineral boundaries. K-feldspar (Kfs) disappears abruptly in the transition zone ( Fig. 4a, b), while round, globular igneous quartz continuous straight into the whiteschist. This allows to visually track the geometry of the whiteschist-metagranite transition and confirms that the protolith of the whiteschist was indeed the Monte Rosa granite. The metasomatic zones can be straight (Fig. 4a) or lobate (Fig. 4b). The pegmatite and aplite dykes, which crosscut the metagranite-whiteschist boundary, are also turned into whiteschist (Fig. 4a), preserving their outlines in places due to their initially different grain size (Fig. 3b). These field observations are important evidence that excludes a xenolith origin of the whiteschist and confirm the previous interpretations that the metagranite is the protolith of the whiteschist (Pawlig and Baumgartner 2001;Marger et al. 2019).
The central, inner part of the whiteschist outcrop (hereafter referred to as zone A in ) consists of an assemblage of chloritoid, phengite, talc, and quartz, defining the peak Alpine paragenesis. Fine-grained kyanite, sericitic phengite, and chlorite are present as retrograde products in minor amounts (< 5 vol %), replacing chloritoid. Accessory minerals include zircon, tourmaline, monazite, florencite, allanite, apatite, rutile, pyrite, and iron oxide. The outer part of the whiteschist is labelled hereafter Zone A out . It becomes successively more retrogressed towards Zone B. Here, chlorite and white mica pseudomorphose chloritoid (Fig. 5b), chlorite flakes pseudomorphose talc, and fine-grained white mica and chlorite are found in the matrix, together with quartz. The subdivision of the Zone A into two subzones (inner, outer) is done to separate the well-preserved high-pressure paragenesis in the centre of the whiteschist outcrop from the greenschist-facies retrogressed rim of the whiteschist outcrop. The transition between the whiteschist in Zone A and Zone B is marked by the complete disappearance of chloritoid-even pseudomorph outlines are no longer visible-and the appearance of plagioclase neoblasts, leading to an assemblage of chlorite, phengite, quartz, and plagioclase (Fig. 5c). Zone C marks the disappearance of chlorite and appearance of calcite. Calcite is present in the matrix, or as distinct aggregates of crystals (Fig. 5d, e) or in veins. The mineral assemblage characteristic of Zone C is phengite, quartz, plagioclase, and calcite, with locally small amounts of biotite and garnet. K-feldspar appears a few tens of centimetres to few metres further away, in the porphyritic metagranite. The porphyritic metagranite preserved its igneous mineralogy of K-feldspar, plagioclase, biotite, and quartz (Fig. 5f). Metamorphic minerals are phengite, small garnets, and a fine-grained assemblage of zoisite/ clinozoisite, albite with white mica pseudomorphically replacing igneous-textured plagioclase .
A second small ellipsoidal body of around 1 square meter in size is found next to the main whiteschist described above and shares the same characteristics, except the absence of Zone C.
Overall, from structural and petrographic observations on the Mezzalama whiteschist, the whiteschist body is interpreted to have a tube-like geometry (Fig. 3d). Mineralogical zonation as described above at the whiteschist-to-metagranite transition is typical for a metasomatic alteration front (e.g., Korzhinskii 1968;Dipple et al. 1998).
The fact that the inner part of the whiteschist remained undeformed during Alpine metamorphism suggests that the crystallization of the peak pressure mineral assemblage in the whiteschist occurred under static conditions. The  Darbellay (2005). The three whiteschist outcrops (Mezzalama, Vera and Boundary) are labeled by arrows observation that the whiteschist is not related to any shear zone or fault, combined with the fact that the protolith of the whiteschist was not deformed before high-pressure metamorphism, suggests that the initial shape of the whiteschist protolith was not a planar structure, but rather a tube-like structure in 3D.

Boundary whiteschist
The second whiteschist occurrence is located at the boundary between the metagranite and the surrounding polymetamorphic paragneiss on the NW of the Mezzalama refuge (624934/085013), just a few hundred meters northwest of the main Mezzalama whiteschist. Due to intense deformation, the whiteschist shows a strongly recrystallized and foliated greenschist facies mineralogy in most places. The assemblage of chlorite, white mica, and quartz and the presence of relict chloritoid are similar to the A out of the Mezzalama whiteschists. Accelerated snow melting has revealed a further undeformed whiteschist just to the South of this location, along the contact, which due to its recent discovery, could not be studied here.  The picture is oriented NW -SE and corresponds to the red square in (b). Zone A consists of whiteschist, zone B is marked by the appearance of plagioclase and disappearance of talc and chloritoid or their pseudomorphs, zone C is marked by the appearance of calcite, ± bio-tite and garnet, and the disappearance of chlorite. The boundary to the porphyritic metagranite is marked by the disappearance of calcite, and the appearance of K-feldspar. The width of the transition zone is variable (0 to 6 meters). Details of the mineralogy can be found in Table 1. d 3D sketch illustrating the interpretation of the field geometry of the Mezzalama area. The outcropping shape is a 2D ellipsoid elongated in the NE -SW direction. Note that no shear zones are associated with the whiteschist. The steep boundaries around the outcrop implying a tube-like shape. The red line shows the the location of the samples profile discussed in Fig. 8a). The drawn lines illustrate the gradient in deformation intensity, decreasing from the NW to SE. Mineral abbreviations are: bi: biotite, cc: calcite, chl: chlorite, ctd: chloritoid, gt: garnet, Kfs: K-feldspar, plg: plagioclase, qz: quartz, tlc: talc, wm: white mica Page 7 of 24 10

Vera whiteschist
The third whiteschist outcrop is located to the SE of the Mezzalama refuge, on the left bank of the Piccolo Ghiacciaro di Vera, above a small lake (625150/ 084184). It was first described by Dal Piaz (1971) and studied in detail by Pawlig and Baumgartner (2001). This whiteschist is crosscut by a late Alpine top to the SW-dipping shear zone. The deformation overprints the whiteschist peak metamorphic assemblage, retrograding it locally into a greenschist facies assemblage. This top to the SW dipping shear zone is coherent with the backfolding, which occurred during greenschist facies metamorphism (Steck et al., 2015). The whiteschist is approximately 50 m long and 15 m wide, elongated in the NW-SE direction. The whiteschist body is slightly deformed by a schistosity that crosscuts the whiteschist-metagranite transition zones. Some blocks of whiteschist remained undeformed. The crosscutting relationship and associated greenschist overprint clearly shows that the deformation post-dates peak high-pressure metamorphism, as it is also the case for the Mezzalama whiteschist. The peak metamorphic assemblage consists of chloritoid, talc, phengite, quartz, and locally some kyanite is found. The deformed, retrogressed, whiteschist is similar to the Mezzalama whiteschist described from Zone A out . Zircon, monazite, allanite, tourmaline, rutile, iron oxide and pyrite are accessory minerals. The central part of the whiteschist was altered during pre-high-pressure metasomatism in the vicinity of a vein-like structure, which now contains sulfides, surrounded by a layer of garnet-rich whiteschist.

Whole-rock geochemistry and stable isotope compositions
The field observations described clearly show that the whiteschists in the Monte Rosa nappe are the result of metasomatic alteration of the granite protolith, in agreement with the interpretations of Pawlig and Baumgartner (2001). This conclusion was also reached by Marger et al. (2019), based on a comparative study of oxygen isotope and major-element compositions of tourmaline found both in whiteschists and metagranites. Below, we present whole-rock chemical and isotopic data on the granite protolith and its metasomatic alteration, the whiteschist. To estimate the compositional variations due to the alteration of the granite protolith, the ISOCON method was used (Grant 1986), combined with the least-squares approach of Baumgartner and Olsen (1995). It is worth noting here that any mass balance calculation requires additional constraints, as discussed in detail by Gresens (1967). He showed that enrichment in an element, i, can be obtained through a decrease of another element in the rock, even if element i was not at all transported; inversely, the concentration of element i can decrease in a metasomatic rock solely due to addition of another element j. Indeed, it is easy to show that, besides density and element concentrations of the altered Fig. 4 a Picture of the contact between zone C (top) and metagranite (bottom). The sharp boundary is marked by the disappearance of milky-white igneous K-feldspar (Kfs) grains form the granite protolith. Two dykes crosscutting the granite have also been turned into whiteschist. b Picture of a lobate transition between the granite and the whiteschist. Note the sharpness of the contact. Some larger quartz globules found in the granite can be traced into the whiteschist at this location and initial rocks, one additional equation is needed to mathematically obtain a unique solution (Gresens 1967;Grant 1986;Baumgartner and Olsen 1995;Ague and Haren 1996;Mukherjee and Gupta 2008). Typical assumptions leading to this additional equation might include constant volume or immobile elements such as Al or Zr (Dipple et al. 1990;Baumgartner and Olsen 1995;Ague 2003). The ISOCON approach by Baumgartner and Olsen (1995) uses a statistical approach. First, the variability of the protolith and metasomatic rock is established through multiple compositional analysis. The variance of the mean is then used to establish all possible element combinations, which are, within This approach is based on the principle that alteration of a rock will remove or add elements that are highly mobile, while other elements are considered to be immobile. The larger the group of elements is that behaves in a similar fashion, the more likely it is that they can be interpreted as immobile elements by this statistical approach. In addition, in the case of several possibilities, a set of immobile elements should be selected based on past studies that have suggested that certain elements are immobile (e.g., Al, Zr, Ti, etc.), adding geochemical knowledge to the statistical interpretation. For more detail, the reader is referred to Baumgartner and Olsen (1995). The results of the calculations can be graphically illustrated in an ISOCON diagram, in which the concentration of the elements in the unaltered rocks is plotted versus that of elements in the altered rock (Grant, 1986;Fig. 6). Here, we chose a logarithmic representation. Note that an ISOCON line, connecting the elements considered to be immobile, has a slope of one in this logarithmic diagram, while the y-axis value at log(c 0 ) = 0 is the mass change in the system (Baumgartner and Olsen 1995). Note that error bars are on a log scale and, therefore, are asymmetric or smaller than the marker size. The line defining the ISOCON passes through all immobile elements within the uncertainty, that is elements that are assumed not to change their concentration during alteration.
A total of 4 porphyritic metagranites, 3 samples from Zone B, 4 whiteschists from the Zone A out , and 8 whiteschists from Zone A in were selected, all from the Mezzalama area, for bulk rock major-and trace-element measurements (Tables 2, 3). 83 possible combinations of immobile elements are proposed by the ISOCON software (Baumgartner and Olsen 1995). The 25th first proposition of a combination of elements are listed as example in Online Supplemetary material S1. The chosen combination of immobile elements for the alteration of the metagranite-to-whiteschist zone A in is: Al 2 O 3 , FeO, Fe 2 O 3 , TiO 2 , and P 2 O 5 as major elements, and Sc, Ni, Y, Zr, Nb, Mo, La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, Hf, Th, and U as immobile trace elements. There is another solution to the ISOCON with the same number of elements, which results in Al 2 O 3 being slightly mobile, while SiO 2 becomes immobile. The results are very similar with respect to the mobile elements, but show a slightly different overall mass balance. The Alconserved scheme results in a 5 ±5% mass gain, while a silica immobile reference frame results in a 10 ±5% mass gain. Note that we do not use volume changes, since the volume of a rock depends on the mineralogy of the alteration process. Since the mineralogy changed during Alpine metamorphism, it is not possible to determine volume change without knowing the initial metasomatic mineralogy.
The results indicate an enrichment in MgO and H 2 O and a depletion of Na 2 O, CaO, Cu, Pb, Zn, Sr, and Ba to obtain the whiteschist from a granitic precursor. The results obtained agree with the previous ISOCON study from Vera whiteschist locality, which was studied by Pawlig and Baumgartner (2001), who found a mass gain of 17 ±6% for the granite-to-whiteschist alteration. Their study revealed the same trends, with MgO and H 2 O enrichment and CaO, Na 2 O depletion as well as immobile K 2 O. Hence, the two whiteschist outcrops were altered by fluids of a very similar chemical composition.
The relative enrichment in MgO and depletion in CaO and Na 2 O indicate that the fluid responsible for the metasomatic alteration of the granite was oversaturated in MgO, leading to MgO precipitation, and it was undersaturated in CaO and Na 2 O, resulting in leaching of these elements from the granite. As K 2 O and SiO 2 behave as immobile elements in the Vera whiteschist and the Mezzalama whiteschist, respectively, the fluid was saturated in these elements. A similar ISOCON plot was made for whiteschists belonging to Zone A out resulting in a similar compositional trend, only Fe 2 O 3 is somewhat lower in the metasomatic rock than in the granite protolith (Fig. 6c). This pattern reflects a similar fluid to be responsible for the composition of external parts of Zone A. Since zone A out and zone A in are very similar in original mineralogy, and the main difference is the greenschist facies overprint, this confirms that no major mass transport occurred during the late greenschist facies Alpine overprint.
ISOCON plots were made for rocks from Zone B (Fig. 7), which are intermediate between metagranites and whiteschists (Fig. 3c). Their composition was compared to the metagranite (Fig. 7A) and to the whiteschist (Fig. 7b). The results show that rocks from Zone B are closer in composition to the whiteschist (Fig. 6a) as they share the same compositional characteristics with an enrichment in MgO Thin section pictures illustrating the main petrographic characteristics of the different rock types found along the transition zone from the whiteschist to the porphyritic metagranite. Most of the rocks belonging to the Mezzalama whiteschist are undeformed or have been only slightly deformed during a post peak-pressure deformation event. a Zone A in (inner whiteschist); ctd, tlc, wm, chl and qz mark the well-preserved high pressure paragenesis (sample 14MR7). b Zone A out (outer whiteschist); ctd and tlc are partially to completely replaced by fine-grained chl and wm during the greenschist-facies, late Alpine retrogression. The outlines of ctd and tlc pseudomorphs are still visible. Chl is becoming more abundant in samples towards the outer parts of the whiteschist. c zone B; plagioclase porphyroblasts embed chl, wm and qz. d zone C; cc and qz vein, with small amount of wm; e zone C; greenschist-facies metagranite with interstitial matrix cc, plg, wm, qz, gt and bi; f Porphyritic metagranite, mostly undeformed, containing igneous Kfs, plg (pseudomorphosed by a fine-grained aggregate of zoisite, albite and white mica), bi, qz and fine-grained gt. Mineral abbreviations are: bi: biotite, cc: calcite, chl: chlorite, ctd: chloritoid, gt: garnet, Kfs: K-feldspar, plg: plagioclase, qz: quartz, tlc: talc, wm: white mica ◂ and H 2 O and depletion in CaO as well as Fe 2 O 3 . The similarity with the whiteschist from Zone A indicates that Zone B results from the alteration of the granite by the same fluid. However, rocks from Zone B have Na 2 O concentrations closer to the composition of the metagranite, and hence, the composition of the fluids becomes saturated in sodic plagioclase as the fluid reacted with the granite.
Whole-rock composition of 29 samples collected on a profile from the Mezzalama whiteschist to the porphyritic metagranite and the whiteschist at the contact with the paragneiss were measured (Table 4). The profile is shown in Fig. 3d by the red line. Based on the mass balance considerations above, the oxides MgO, Na 2 O, CaO, K 2 O as well as H 2 O are plotted. Figure 8a illustrates that the sodium front is located at the A/B zone boundary, while the MgO and CaO fronts are located at the B/C zone boundary. The stepwise changes along mineralogical fronts are characteristics of chromatographic alteration fronts in metasomatic environments described by the chromatographic theory (Korzhinskii 1968;Steefel and MacQuarrie 1996). δ 18 O and δD values were measured on 14 samples from this profile (Table 5). A distance versus δ 18 O (Fig. 8b, left vertical scale) and δD values (Fig. 8b, right vertical scale) plot shows that the isotopic compositions were changed throughout zones A out to C, but have constant values in the metagranite and in most of the samples in zone A in .  Compositions are in wt% for major elements. A mass change of 4.9 ± 5.0% is calculated using the following combination of immobile elements: Al 2 O 3 , FeO, Fe 2 O 3 , TiO 2 , and P 2 O 5 . b Result of the same calculation for trace elements, in ppm. Mobile trace elements are highlighted with filled symbols. c ISOCON diagram for major elements of the granite protolith ver-sus the external part of the whiteschist (Zone A out ). A mass change of -13.3 ± 7.5% is calculated. Note that all scales are logarithmic, as well as error bars. Error bars represent the standard deviation of the mean on the n analyses considered for both the granite and whiteschist (1SD). In both whiteschists (inner and outer), MgO and H 2 O are enriched, CaO and Na 2 O are leached. This implies a fluid oversaturated in MgO and undersaturated in CaO and Na 2 O The δ 18 O values of ca. 11.5‰ of the granite are typical for S-type granites (Harris et al. 1997). The stable isotope compositions decrease by about 3‰ for δ 18 O and increase by 20‰ for δD from metagranite to the whiteschist, suggesting that the system was at least partially fluid buffered. This is also suggested by the fact that whiteschists from both locations (Mezzalama and Boundary) have similar δ 18 O and δD values.

Carbonate stable isotope compositions
Calcite is found in both the Mezzalama whiteschist and the Vera whiteschist. In the Mezzalama area, they appear in Zone C of the transition zone, either as veins or as matrix carbonate. The carbonate is calcite in all cases. δ 13 C and δ 18 O values have been measured on six samples of carbonate veins and carbonate-bearing whiteschists from the Mezzalama area. Results are given in Table 6 together with those from 5 samples of the Vera whiteschist (Pawlig and Baumgartner 2001) for comparison. Mezzalama carbonates have a range from -7.6 to -6.7‰ in δ 13 C and from 9.3 to 10.8‰ in δ 18 O. Carbonates from Vera whiteschist from Pawlig and Baumgartner (2001) have a range from -9.8 to -9.6‰ in δ 13 C and 9.6 to 10.1‰ in δ 18 O. Carbonates from the Mezzalama whiteschist have slightly higher δ 13 C than the Vera whiteschist, for comparable δ 18 O. A plot of these isotopic compositions in a δ 13 C versus δ 18 O diagram (Fig. 9) supports a magmatic origin (Valley 1986;Hoefs 1987;Pawlig and Baumgartner 2001).

Rb-and Sr-isotope compositions
Several metagranites and whiteschists from the Mezzalama area have been selected for strontium (Sr) isotope analyses. Since the Rb/Sr ratios are much higher in the whiteschist than those of the granite, we expected to be able to determine the age of the hydrothermal alteration using a classical Rb/Sr isochrone. Whiteschist are from Zone A in (Table 1). The metagranites were collected from various parts of the study area and are undeformed metagranites. As shown in Fig. 6b, Sr was one of the most mobile elements during metasomatism of the granite. It was depleted along with Ca. In contrast, rubidium (Rb) behaves geochemically similar to potassium and was hence little affected by the metasomatism. Depending on the Rb/Sr of a rock, the isotopic ratios of 87 Rb/ 86 Sr and 87 Sr/ 86 Sr will have evolved differently. The larger the Rb/Sr ratio, the larger the change in 87 Sr/ 86 Sr with time. Ideally, one would obtain a Rb/Sr isochron if the system remained closed after metasomatism.
The 87 Sr/ 86 Sr ratios, in combination with Rb and Sr (ppm) content, were used to calculate 87 Rb/ 86 Sr ratios (Table 7; details in analytical methods) and plotted in a 87 Rb/ 86 Sr versus 87 Sr/ 86 Sr diagram (Fig. 10). Data from the Mezzalama whiteschists were compared with data from the Vera whiteschists (Pawlig 2001  Compositions are in wt% for major elements. Note that the scale is logarithmic, as well as error bars. Error bars represent the standard deviation of the mean on the n analyses, reported in Table 2a, considered for both the metagranite and whiteschist (1SD).
The alteration which produced zone B is similar to that of the whiteschist of zone A (Fig. 6a,  The resulting trends are clear, but unfortunately not linear, which precludes the use of the data for an isochron. It is well established that curved lines, as the one obtained for the present data set in Fig. 10, are typical for systems affected by multiple events; the isochron has been disturbed by later processes (e.g., Faure and Powell 1972).

Discussion
The field observations showed that: (i) the transition from whiteschist to metagranite is characterized by two zones, which have sharp mineralogical boundaries, (ii) the disappearance of K-feldspar from the igneous protolith towards the whiteschist is sharp and outlines a lobate boundary, while igneous quartz aggregates and crystals can be traced into the whiteschist zones; (iii) straight, undeformed pegmatite and aplite dykes crosscut the granite-whiteschist transition. The mass balance calculations performed above show that over 20 trace-and rare-earth elements, as well as Ti and Al can be considered to be immobile, when the whiteschist and granite are compared. This all confirms that the protolith of the whiteschist was indeed the granite. Pawlig and Baumgartner (2001) came to the same conclusion for the Vera whiteschist outcrop, where they suggested that two different granites were involved in the formation of slightly different whiteschist lithologies. Finally, Marger et al. (2019) demonstrated that tourmaline in the whiteschists has cores that correspond in zonation, major-element composition, and δ 18 O value to tourmaline in the adjacent granite for the Vera whiteschists. The field observations thus confirm that the whiteschists are metasomatic alteration zones within the granite. Hence, this further supports that both for the Vera and the Mezzalama whiteschists in the Monte Rosa nappe the protoliths are the granites.

Timing of the metasomatic event
Unfortunately, the Rb/Sr data do not define an isochron. The data form a curved array, suggesting that multiple metamorphic events influenced the Rb/Sr systematics. Hence, no precise age for the metasomatism can be defined. The strontium evolution diagram can, however, be used to obtain further information on the metamorphic evolution of the rocks concerned. Results of both this study as well as that of Pawlig (2001) are shown in Fig. 11. The 87 Sr/ 86 Sr ratios are homogeneous for the granite (0.7331-0.7358), indicated by the heavy line in Fig. 11. The 87 Sr/ 86 Sr evolution is also calculated back to 270 Ma, the age of the granite (Pawlig 2001), using today's Rb/Sr elemental ratio for the Monte Rosa granite. Unaltered granite evolved along this line with time. In a simple model, in which one metasomatic event increased the Rb/Sr ratio by decreasing the Sr concentration, as indicated by the mass balance calculations (without fractionation of the isotopes), one would expect one intersection for all whiteschist sample evolution lines with the granite evolution line; the intersection would define the age of the metasomatic event. Since the model whiteschist contains significantly more Rb, the radiogenic daughter, 87 Sr, would show a larger increase and plot above the granite evolution line, yielding the high 87 Sr/ 86 Sr measured today. Of course, the same age should be obtained, as defined by the isochron in a traditional 87 Sr/ 86 Sr versus 87 Rb/ 86 Sr diagram for a closed system. However, as this is not the case, the system was disturbed, resulting in multiple intersections with the granite for the various whiteschist samples. Nevertheless, all intersection points are older than the Alpine high-pressure metamorphism. They span about 200 Ma, from 80 Ma to 280 Ma. Alternatively, if a single, pervasive event homogenised the 87 Sr/ 86 Sr composition of the whiteschist after metasomatism-such as the prograde high-pressure metamorphic overprint-one would expect intersection of all 87 Sr/ 86 Sr evolution lines well above the granite evolution line, indicating the age of the homogenisation event. While multiple samples do intersect each other, no clear pattern evolves. Nevertheless, all whiteschist samples indicate that alteration occurred well before the Alpine metamorphism. The intersections in the age range of the Alpine metamorphism for several of the whiteschist samples could indicate a partial homogenization during devolatilization reactions during the prograde path of the Alpine high-pressure metamorphism.
A less likely, but possible alternative is that high 87 Sr/ 86 Sr could have been imposed from the alteration fluids; alternatively, a multiphase evolution, with a partial homogenization Fig. 8 a Whole rock composition of samples for the profile sampled at Mezzalama, from the whiteschist towards the contact of the metagranite with the country rock paragneisses, where a second, small whiteschist occurs (red line on Fig. 3d). The mobile elements identified in Fig. 6 and 7 (CaO, Na 2 O, MgO and H 2 O content (wt%)) are shown as a function of distance. Note that K 2 O does not change across the reaction front. b δD and δ 18 O versus distance. c δD versus δ 18 O for the same samples as (b) showing two distinct isotopic signatures, with high δD and low δ 18 O for whiteschists and low δD and high δ 18 O for metagranites. Note that the stable isotope fronts have penetrated further into the granite than the chemical alteration fronts. The whiteschist at the contact has the same isotopic composition as the main Mezzalama body ◂ could also yield these results. In this scenario, the rocks were infiltrated by fluids, carrying significant but variable amounts of radiogenic strontium. Typical sources of such radiogenic strontium would be an old crustal reservoir, such as the metapelites in the basement making up the host rocks of the Monte Rosa granite (Faure 1998).
We conclude that the strontium isotopic system was disturbed during the Alpine history, due to recrystallization and potential post-peak re-equilibration with fluids. Based on the strontium isotopic data, it is also clear that whiteschist alteration occurred prior to the Alpine metamorphism. The fact that all intersections are above the granite evolution line, e.g., the whiteschists have higher 87 Sr/ 86 Sr ratios, requires that a significant amount of time elapsed between alteration and partial homogenisation due to Alpine metamorphism, such that rocks with high Rb/Sr ratios could evolve towards high radiogenic 87 Sr. This alteration could have occurred during a late magmatic hydrothermal stage, shortly after emplacement of the granite, and during its cooling. This interpretation is in line with the results of Marger et al. (2019), and Pawlig and Baumgartner (2001), who showed, on the basis of microtextural relationships and tourmaline compositions, that the metasomatic alteration of the granite occurred before the onset of the Alpine metamorphism. However, due to the heavy Alpine metamorphic overprint, we cannot exclude that the metasomatism occurred later during the Mesozoic evolution of the Monte Rosa basement, yet still prior to the Alpine orogeny.

Nature of the metasomatic fluid altering the granite
Our new data based on field observations, radiogenic isotope compositions, and timing of the metasomatic alteration in the Monte Rosa nappe, as well as the stable isotope and whole-rock compositions require a pre-Alpine alteration. Considering the geologic history of the Monte Rosa nappe, we further consider two possible alternative scenarios: I) either seawater incursion during Mesozoic extension (e.g., Manatschal et al. 2000) or II) magmatic-hydrothermal fluids (Pawlig and Baumgartner 2001). Metasomatism by seawater fluids would occur during the Mesozoic extension, leading to the opening of the Piemont oceanic basin (Trümpy 1980). As part of the distal margin of the Briançonnais microcontinent (Steck et al. 2014), the Monte Rosa basement rocks are exposed along the ocean-continent transition. In this tectonic setting, seawater fluids could penetrate and locally interact with ultramafic rocks exhumed at the seafloor and Vera (Pawlig et al. 2001) Fig. 9 δ 13 C and δ 18 O isotopic compositions of calcite from Zone C in Mezzalama area (data in Table 6), indicating an igneous carbon source for both whiteschist outcrops. Data for the Vera outcrop come from Pawlig and Baumgartner (2001). Data for isotopic values of sedimentary carbonates and igneous carbonates are from Valley (1986) and Hoefs (1987)  interact on their way back to the surface with the Monte Rosa basement along large-scale detachment structures (Manatschal et al. 2000;Ferrando 2012). Metasomatism by late magmatic-hydrothermal fluids would occur during the cooling of the granite intrusion in the Permian. In such hydrothermal systems, surface-derived fluids mix with primary magmatic fluids in large-scale fluid convection cells, affecting the basement rocks as a result of the thermal perturbation associated with the igneous intrusion. An exact determination of the nature of the metasomatic fluids is hampered by the lack of the original alteration phases formed during alteration, since subsequent Alpine metamorphism crystallized new phases. Nevertheless, it is not straightforward to interpret the elemental changes. Indeed, the chemical changes introduced by fluid flowing through the rock are inversely proportional to the concentration changes observed in the fluid responsible for alteration along its flow direction (Phillips 1990); for more detail, see, for example, Baumgartner and Ferry (1991). In the literature, two endmember cases have been discussed for understanding and interpreting metasomatism: i) infiltration of a disequilibrium fluid driving mineral/rock composition changes (e.g., Korzhinskii 1968;Hofmann et al. 1974, Dipples andSteefel and MacQuarrie 1996); or ii) metasomatism due to changes in solubility as a consequence of pressure and temperature changes in flow direction (e.g., Phillips Table 2 Major element compositions in wt% of the metagranites, whiteschists (inner part), whiteschist (outer part) and rocks of the transition zone B from the Mezzalama area These data were used for mass balance calculations with the statistical ISOCON approach of Baumgartner and Olsen (1995)  1990; Baumgartner and Ferry 1991;Steefel and MacQuarrie 1996). Case (i) is readily understood as a fluid emanating for example from a typical ultramafic rock will be in contact with high-solubility Mg-bearing minerals and hence contain a high amount of total Mg. In contrast, the same fluid will be highly undersaturated with potassium, since serpentinites do not contain significant amounts of K-bearing minerals. If this ultramafic rock derived fluid infiltrates a granite it may produce a Mg metasomatism and, at the same time, leach potassium, for example. The second case is due to the fact that pressure and temperature change the solubility of minerals, and hence, flow in a P/T gradient results in concentration changes in the fluid, reflected by dissolution/precipitation of minerals in the rock. For example, a chloride-bearing fluid flowing down a temperature gradient-i.e., fluid is flowing towards lower temperatures-in a K-feldspar-albite-muscovite-quartz-bearing granite or pelitic rock will preferentially dissolve Na, and precipitate K, since the solubility of Na at equilibrium increases, while that of K decreases (e.g., Can 2002;Miron et al. 2016). A quantitative evaluation of the solubility of minerals is difficult, since the high-pressure metamorphism has changed the mineralogy of the protholith through closed system metamorphism, and hence, we propose a qualitative discussion here. The fact that the major elements-aluminium, silica, and even potassium-were neither added nor leached from the rock (Figs. 6, 7) suggests a fluid which has reached equilibrium with respect to these components in the granite, while sodium and calcium-mostly present in plagioclasewere leached and replaced by a Mg-rich silicate. A similar Table 4 Whole rock compositions (in wt%) of samples along the profile starting at the Mezzalama whiteschist (WS) towards the host rock paragneisses in the NW (red line on Fig. 2d alteration was observed by Baumgartner and Olsen (1995) for the porphyry copper Bingham (Utah, USA) sericite-biotite alteration zone, where mostly sodium and calcium were leached, together with minor iron, yielding a potassium-rich rock. We propose a similar alteration, but occurring at a slightly lower temperature, in the chlorite field. In that case, Mg would have been increasing as well, since chlorite in equilibrium with the same fluid as biotite is generally more Mg-rich (e.g., Spear 1995). A fluid emanating directly from serpentinites, as proposed by other workers for the whiteschists of the internal crystalline massifs (Sharp et al. 1993;Demény et al. 1997;Barnes et al. 2004;Ferrando et al. 2009), would have resulted in leaching of K 2 O and SiO 2 and an enrichment of Ni and Cr, which is not the case for the Monte Rosa whiteschists.
Evaporite-derived fluids, as suggested by Ferrando et al. (2009), could be a source for high Mg fluids, but they would have also high concentrations of sodium and calcium due to potential halite, dolomite, or calcite contents of the sediments. Since these elements were found to be leached from the granite protolith, such a source is unlikely.
Similarly, seawater-driven serpentinization of the mantle at the base of the continental crust during Mesozoic extension would be enriched in sodium and likely calcium from seawater in addition to magnesium. Hydration of the mantle would increase the salinity additionally, resulting in sodiumrich brines. A further argument against such a scenario is the missing of pervasive brittle structures, which would grant access to the interior part of the Monte Rosa rocks. Following Manatschal et al. (2000), we would expect fluids to be channelized along detachment faults, rather than undulating, reaction fingering type structures .
Finally, seawater fluids do not bear isotopic heavy strontium signatures (seawater 87 Sr/ 86 Sr = 0.709; Faure, 1972). Therefore, they would result in lower radiogenic Sr values during alteration, not seen in the measured high 87 Sr/ 86 Sr ratios of the whiteschists (Fig. 10). Consequently, field observations and whole rock as well as radiogenic isotopes argue against fluids during the Mesozoic extensional regime. This leaves, as most probable metasomatic event, the late magmatic-hydrothermal alteration of the granite.
The whiteschist δD and δ 18 O values can be used to estimate the composition of the fluid in equilibrium with their protolith, provided that the mineralogy of the metasomatic Oxygen isotope fractionation coefficients (used as ∆ min-H2O ) are from Zheng (1993) for chloritoid, talc, phengite, quartz, and from Cole (1985) for chlorite. Hydrogen isotope fractionation coefficients (∆ min-H2O ) are from Graham et al. (1987) for chlorite and from Suzuoki and Epstein (1976) for muscovite. The mole fractions of oxygen and hydrogen are calculated using estimated modal proportions to determine the whole-rock-H 2 O fractionation factor. The results are shown in Fig. 12. The metasomatic fluid compositions have δD and δ 18 O values that match a mix between primary magmatic fluids and meteoric fluids. Such fluids are consistent with a late magmatic-hydrothermal fluid, derived from the mixing of late magmatic and surface-derived fluids. Through circulation within the Monte Rosa paragneiss basement, fluids could have been enriched in radiogenic heavy strontium before reaching the granite, but highly radiogenic 87 Sr/ 86 Sr ratios could also have been obtained solely by Sr-leaching during metasomatism. Such fluids are highly reactive and susceptible to alter granites through reaction-induced fingering (e.g., Dewers and Ortoleva 1990;Omlin et al. 2017). Late magmatic-hydrothermal fluids succeed in explaining both field observations, major and trace elements, as well as stable and radiogenic isotopes. Moreover, they fit with the tectonic setting of the Monte Rosa basement during the Permian.

Conclusion
A comparison of the whole-rock chemical composition of the whiteschists and the granites confirms that the protolith of the whiteschists were produced in-situ from the enclosing granite. The strontium isotope compositions require a pre-Alpine age for the metasomatic alteration of the granite, in agreement with data presented in Marger et al. (2019). Hence, a subduction fluid origin as proposed by Sharp et al. (1993) and Ferrando et al. (2009) Table 7, for the Mezzalama whiteschist (full lines) and the Vera whiteschist (dashed lines; after Pawlig 2001). The dark grey band represents the duration of the Alpine orogenic event (Hurford et al. 1991;Lapen et al. 2007;Skora et al. 2015). The intersection between the lines for each sample indicate that the Rb/ Sr system was perturbed. We suggest that the Sr isotopic system is perturbed during the Alpine history by local, limited re-equilibration of the Sr isotopic system. This is likely due to dehydration reactions during the prograde reactions leading up to high-pressure metamorphism. Some additional re-equilibration might have occurred during the partial greenschist-facies overprint. On the right, two sketches illustrate an undisturbed isotopic system (above) and a partially reequilibrated system due to an isotopic perturbation (below) in order to clarify the interpretation compositions of the whiteschist and the chemical signature do not support such an origin.
The above-detailed discussion rather suggests that the major-element compositions of the whiteschists reflect late magmatic-hydrothermal alteration of the granite resulting in a sericite-chlorite-quartz alteration of the granite. The fluids were undersaturated with respect to Na and Ca, saturated with respect to K and Si, and oversaturated with respect to Mg. This is either related to their origin, that is their source rocks, or due to the changes in intensive variables (P, T) along the flow directions and intense fluid-rock exchange. Stable isotope compositions support fluids that have equilibrated with an igneous rock. A second, minor fluid event related to Alpine metamorphism is also suggested by the strontium isotope compositions. This leads us to propose the genetic model below.

Genetic model
The genetic model proposed involves four steps, as illustrated in Fig. 13: (i). the intrusion of the granitic magma in the upper crust during the Permian (Pawlig 2001) with an associated large-scale thermal perturbation in the country rocks.
A contact aureole was developed in the country rocks around the granite pluton (Fig. 13a). (ii). The thermal perturbation sets up a convection (e.g., Hanson 1995;Cook et al. 1997;Gerdes et al. 1998) leading to infiltration of surface-derived fluids. Fluids reaching the granite is focussed along the granite-paragneiss contact, and locally penetrated into the granite by reaction-induced fingering. Due to the large chemical disequilibrium between the granite and the fluid, pervasive alteration of the granite is induced along tube-like structures and/or in contact zones facilitating fluid access (Fig. 13b). At this point, the stable isotope composition is reset within the rocks, reflecting a typical hydrothermal mineralogy and composition (e.g., Baumgartner and Valley 2001). The rock produced is a phyllitic or argillitic alteration containing sericite, chlorite, quartz, or smectite-illite and kaolinite, depending on the prevailing temperature during alteration. In this kind of alteration zone, all plagioclase is dissolved, and K-feldspar altered to muscovite/ illite, and biotite is altered to chlorite by these slightly acidic hydrothermal solutions (e.g., Parry et al. 2002). The presence of sulphides described by Pawlig and Baumgartner (2001) further supports this interpretation. Strontium was, along with calcium, depleted due to plagioclase breakdown, while rubidium follows potassium into the muscovite/illite. Potentially, some radiogenic 87 Sr could have also been leached from the surrounding pre-Hercynian schists, though this cannot be demonstrated. During the following ca. 200 Ma years, the system is left undisturbed, accumulating radiogenic 87 Sr and thus increasing the 87 Sr/ 86 Sr ratios. How much radiogenic 87 Sr is produced will depend on the local Rb concentration. The overall chemical composition of the whiteschists including their bulk stable isotope compositions is not changed or perturbed during this period. (iii). Once the Monte Rosa basement was included into the Alpine orogenic wedge, prograde metamorphism lead to the Alpine high-pressure metamorphic assemblage of chloritoid-talc-phengite-quartz and the granites received a variable metamorphic overprint (Fig. 13c). During prograde metamorphism, the crystallisation of the whiteschist produced H 2 O, leading to local, partial homogenisation of the Rb/ Sr system. The small amounts of H 2 O lost during prograde metamorphism did not significantly change the chemistry or the stable isotope compositions. The partial homogenisation of the Rb/Sr system disturbed any potential isochrones, making an exact age determination impossible. (iv). During the exhumation and backfolding of the nappe, a foliation developed, associated with retrogression of the peak Alpine paragenesis under greenschist-facies conditions, which affected the entire area (Fig. 13d). Small amounts of fluid infiltration Fields for the primary magmatic fluids, the domain of ocean water, as well as the meteoric water line are shown for comparison. The fluids responsible for the metasomatism were likely a mix between surface waters and primary magmatic waters, typical for a late magmatic hydrothermal system occurred at this point, resulting in a local retrograde overprint of the high-pressure Alpine assemblage. Fluid circulation was localized on the metagranite-whiteschist contact. Fluid amounts reaching the interior of the whiteschists were too small to change the whiteschist isotopic compositions (see Luisier et al., in review).
The detailed field and geochemical study presented here confirms the interpretation of Luisier et al. (2019) in that the whiteschist body has indeed experienced the same metamorphic conditions as the metagranite. Hence, the differences in peak metamorphic pressures calculated for the metagranite (1.4 GPa) and the whiteschist (2.2 GPa) need to be explained. A proposed overpressure-either due to fluid overpressure or tectonic stress-as explained in detail in Luisier et al. (2019) may account for such pressure differences.
This genetic model accounts for both field observations and geochemical data of all whiteschist outcrops from the Monte Rosa nappe. As has been pointed out by many workers, whiteschist bodies are common in granitic bodies in the Western Alps (see Ferrando 2012 for a review). While the above discussion and data clearly show that no individual observation-neither field nor geochemistry-permits to uniquely demonstrate the origin of a whiteschist, it is, nevertheless, tempting to speculate that quite a few, if  The thermal perturbation drives a hydrothermal convective circulation system. The fluid circulation pattern is illustrated by arrows. Fluids either penetrate into the granite through focused paths by reaction fingering, leading to channels or tube-like structures, or they follow the contact between granite and paragneiss. We propose that the original alteration assemblage was a chlorite, sericite and quartz rock. c Closed-system Alpine high-pressure metamorphism, during which fluids were expelled, produce the peak paragenesis of the whiteschists, consisting of chloritoid, phengite, talc and quartz. During that stage, the granite becomes a metagranite. d The whiteschist area, and much of the granite, are only slightly affected by the strong foliation overprinting the Monte Rosa nappe during exhumation and backfolding. This lead locally to a greenschist-facies partial retrogression. The form of the crosses illustrates the deformation intensity gradient in the area, decreasing from NW to SE. The cross-section is modified after Darbellay (2005) not the majority of granite-hosted whiteschists might have a similar hydrothermal origin.