Relationship between fluid-bearing and fluid-absent invariant points and a petrogenetic grid for a greenschist facies assemblage in the system CaO-MgO-Al2O3-SiO2-CO2-H2O

In the 6 component system CaO-MgO-Al₂O₃-SiO₂-CO₂-H₂ with 9 solid phases (quartz, plagioclase, epidote, tremolite, talc, chlorite, magnesite, calcite, dolomite) and a fluid phase, all 17 possible fluid-absent reactions have been set up and balanced. Using molar entropy and volume data for the solid phases, these reactions are arranged in P-T space about the 8 possible fluid-absent invariant points after the method of Schreinemakers. Field observations in Ordovician greenschist facies basic volcanics at Sofala N.S.W., indicate that neither talc+epidote nor magnesite+calcite are stable under the conditions of metamorphism. Assuming these conditions to apply to the theoretical study here, the fluid-absent invariant points are arranged in a relative fashion with fluid-absent reactions subdividing P-T space into smaller areas.A scheme which permits a fluid of composition (i.e. a fluid containing CO₂ and H₂O together with other components), is modeled by treating H₂O as a mobile component independent of CO₂, and by allowing values that lie off the locus of binary H₂O-CO₂. Taking into account that neither talc+epidote nor magnesite +calcite is to be permitted, the fluid scheme is used to set up and balance all 39 possible fluid-bearing reactions. These are then arranged about 20 valid fluid-bearing invariant points in space after the method of Korzhinskii and Sehreinemakers.A characteristic solid phase assemblage is defined for each P-T area using chemographic relations inherent from the fluid-absent boundary reactions. The fluid-bearing invariant points that have a solid assemblage compatible with the characteristic assemblage in a particular P-T area are stable within the P-T regime of that area. When these stable fluidbearing invariant points are arranged in a relative fashion in space, they outline a fluid grid which can be used to study the possible effects of local variation in X _fluid over the particular P-T regime.Symbols Used U chemical potential - S entropy - V molar volume - n coefficient of a phase in a reaction - X mole fraction - T temperature - P pressure - F number of degrees of freedom - C number of components - p number of phases - s solid - slope of reaction - 1 quartz - 2 plagioclase - 3 epidote - 4 tremolite - 5 talc - 6 chlorite - 7 dolomite - 8 magnesite - 9 calcite     

Petrology of metamorphic rocks from the Atbashy complex,Southern Tien-Shan,Kyrgyzstan

The high-to ultrahigh-pressure metamorphic rocks of the Atbashy complex were petrologically investigated. The eclogites of the Choloktor Formation show a prograde evolution from epidote-blueschist facies(P = 17-21 kbar and T = 450-515 ℃) to peak eclogite-UHP conditions(P = 26-29 kbar and T = 545-615 ℃) with a subsequent epidote-amphibolite and greenschist facies overprint. The micaschists of the Choloktor Formation also show a clockwise P-T path from blueschist/epidote-blueschist facies conditions through peak eclogite facies conditions(P = 21-23 kbar and T = 530-580 ℃) to retrograde epidote-amphibolite and greenschist facies stages. A comparison of the P-T paths in the eclogites and mica-schists of Choloktor Formation reveal that they may have shared their P-T history from peak to retrograde stages. The mica-schists of the Atbashy Formation record peak metamorphism of P = 10-12 kbar and T = 515-565 ℃, which indicates that the highest grade of regional metamorphism in the Atbashy Ridge was epidote-amphibolite facies.The newly obtained P-T conditions for the mica-schists of Choloktor Formation indicate that sheets of sedimentary rocks were brought to great depths along the subduction zone and they metamorphosed under eclogite facies HP conditions. The eclogite blocks were amalgamated with mica-schists of Choloktor Formation in the eclogite facies HP conditions and together they experienced isothermal decompression to ～40 km. During this path, the eclogites and mica-schists of Choloktor Formation docked with mica-schists of Atbashy Formation at 10-12 kbar and 515-565 ℃, and from this depth(～40 km) the whole sequence was exhumed together. These new results improve our understanding of high-pressure metamorphism in subduction-related accretionary prism zones and the exhumation processes of deeply-seated rocks in the Atbashy HP-UHP complex.     

DasT-P-Feld der Diagenese und niedrigtemperierten Metamorphose aufgrund von Mineralreaktionen

It is attempted to determineT, P values for the lowgrade metamorphic facies. From geologic-petrographic observations and from experiments mineralreactions are outlined which characterize the zeolitic-, the glaucophane- and the greenschistfacies, respectively. Only such reactions are considered which are univariant if P_f=P_s; thus aP-T grid can be arrived at. Experimental data on the equilibria of the relevant reactions is taken from the literature and from own experiments. Experimental results are always checked against known field-observations. Contrary to current opinium, we arrive at rather higher temperatures for the beginning of the lowgrade metamorphic facies: Diagenesis. From sedimentation up to slightly below 300° C. Zeolitic facies. From slightly below 300° C up to about 400° C. Greenschistfacies. From about 400° C up to about 550° C. Glaucophaneschistfacies. Under pressures of at least 6–7 kb this facies begins somewhat above 300° C, probably at about 330° C, grading into high-pressure greenschistfacies.     

A thermodynamic model for multicomponent melts,with application to the system CaO-Al2O3-SiO2

A thermodynamic model is proposed for calculation of liquidus relations in multicomponent systems of geologic interest. In this formulation of mineral-melt equilibria, reactions are written in terms of the liquid oxide components, and balanced on the stoichiometry of liquidus phases. In order to account for non-ideality in the liquid, a ‘Margules solution’ is derived in a generalized form which can be extended to systems of any number of components and for polynomials of any degree. Equations are presented for calculation of both the excess Gibbs free energy of a solution and the component activity coefficients.Application to the system CaO-Al₂O₃-SiO₂ at one atmosphere pressure is achieved using linear programming. Thermodynamic properties of liquidus minerals and the melt are determined which are consistent with adopted error brackets for available calorimetric and phase equilibrium data. Constraints are derived from liquidus relations, the CaO-SiO₂ binary liquid immiscibility gap, solid-solid $\text{P-T}$ reactions, and measured standard state entropies, enthalpies, and volumes of minerals in this system.Binary and ternary liquidus diagrams are recalculated by computer programs which trace cotectic boundaries and isothermal sections while checking each point on a curve for metastability. The maximum differences between calculated and experimentally determined invariant points involving stoichiometric minerals are 17°C and 1.5 oxide weight per cent. Because no solid solution models have been incorporated, deviations are larger for invariant points which involve non-stoichiometric minerals.Calculated heats of fusion, silica activities in the melt, and heats of mixing of liquids compare favorably with experimental data, and suggest that this model can be used to supplement the limited amount of available data on melt properties.     

Zur Stabilität von Margarit im System CaO-Al2O3-SiO2-H2O

An increasing number of occurrences of margarite have been reported in the last years. However, previous experimental investigations in the system CaO-Al₂O₃-SiO₂-H₂O are limited to the synthesis of margarite and to the upper stability limit according to the reaction (1) 1 margarite?1 anorthite +1 corundum +1 H₂O (Chatterjee, 1971; Velde, 1971). Since margarite often occurs together with quartz, the upper stability limit of margarite in the presence of quartz is of special interest. Therefore, the reactions (5) 1 margarite +1 quartz ?1anorthite +1 kyanite/andalusite +1 H₂O and (6) 4 margarite+3 quartz ? 2 zoisite+5 kyanite+3 H₂O were investigated experimentally using mixtures of natural margarite (from Chester, Mass., USA), quartz, kyanite, andalusite, zoisite, and synthetic anorthite. The indicated equilibrium temperatures at water pressures equal to total pressure are: 515± 25°C at 4 kb, 545 ±15°C at 5 kb, 590±10°C at 7 kb, and 650±10°C at 9 kb for reaction (5), and 651±11°C at 10 kb, 648 ± 8°C at 12.5kb, and 643±13°C at 15kb for reaction (6), respectively. Besides this, additional brackets for equilibrium temperatures were determined for the above cited reaction (1): 520±10°C at 3 kb, 580±10°C at 5 kb, and 640± 20°C at 7 kb. On the basis of these experimentally determined reactions (1), (5), and (6) and of the reactions (3) 2 zoisite +1 kyanite? 4 anorthite +1 corundum +1 H₂O (7) 2 zoisite +1 kyanite +1 quartz ? 4 anorthite +1 H₂O and (10) 1 pyrophyllite ? 1 andalusite/kyanite+3 quartz+1 H₂O for which experimental or, in the case of reaction (3), calculated data were already available, a pressure-temperature diagram with 3 invariant points and 11 univariant reactions was developed using the method of Schreinemakers. This diagram, summarizing both experimental and phase relation studies, allows conclusions about the conditions under which margarite has been formed in nature. Margarite is limited to low grade metamorphism at water pressures up to approximately 3.5 kb; in the presence of quartz, margarite is even limited to low grade metamorphism at water pressures up to 5.5 kb. Only at water pressures higher than the values stated before margarite, and margarite+quartz, respectively, can occur in medium grade metamorphism (as defined by Winkler, 1970 and 1973). For the combined occurrence of margarite+quartz and staurolite as reported by Harder (1956) and Frey (personal communication, 1973) it may be estimated that water pressure has been greater than approximately 5.5 kb, wheras temperature has been in the range from 550 to 650°C. Furthermore, the present study shows that the assemblage zoisite+kyanite (+ H₂O) is an indicator of both pressure [P _{H 2 O}> approximately 9kb]and temperature [T> approximately 640 to 650° Cat water Pressures up to 15 kb].     

The thermodynamic properties and phase relations of some minerals in the system CaO-AI2O3-SiO2-H2O

The heat capacities of lawsonite, margante, prehnite and zoisite have been measured from 5 to 350 K with an adiabatic-shield calorimeter and from 320 to 999.9 K with a differential-scanning calorimeter. At 298.15 K, their heat capacities, corrected to end-member compositions, are 66.35, 77.30, 79.13 and 83.84 cal K^?1 mol^?1; their entropies are 54.98, 63.01, 69.97 and 70.71 cal K^?1 mol^?1, respectively. Their high-temperature heat capacities are described by the following equations (in calories, K, mol): Lawsonite (298–600 K): Cp° = 66.28 + 55.95 × 10^?3T ? 15.27 × 10⁵T^?2 Margarite (298–1000 K): Cp° = 101.83 + 24.17 × 10^?3T ? 30.24 × 10⁵T^?2 Prehnite (298–800 K): Cp° = 97.04 + 29.99 × 10^?3T ? 25.02 × 10⁵T^?2 Zoisite (298–730 K): Cp° = 98.92 + 36.36 × 10^?3T ? 24.08 × 10⁵T^?2 Calculated Clapeyron slopes for univariant equilibria in the CaO-Al₂O₃-SiO₂-H₂O system compare well with experimental results in most cases. However, the reaction zoisite + quartz = anorthite + grossular + H₂O and some reactions involving prehnite or margarite show disagreements between the experimentally determined and the calculated slopes which may possibly be due to disorder in experimental run products. A phase diagram, calculated from the measured thermodynamic values in conjunction with selected experimental results places strict limits on the stabilities of prehnite and assemblages such as prehnite + aragonite, grossular + lawsonite, grossular + quartz, zoisite + quartz, and zoisite + kyanite + quartz. The presence of this last assemblage in eclogites indicates that they were formed at moderate to high water pressure.     

A petrogenetic grid for pelitic schists in the system SiO2-Al2O3-FeO-MgO-K2O-H2O

A quantitative petrogenetic grid for pelitic schists in the system KFMASH that includes the phases garnet, chlorite, biotite, chloritoid, cordierite, staurolite, talc, kyanite, andalusite, sillimanite, and pyrophyllite (with quartz, H₂O and muscovite or K-feldspar in excess) is presented. The grid is based on thermodynamic data of Berman et al. (1985) and Berman (1988) for endmember KFASH and KMASH equilibria and natural Fe-Mg partitioning for the KFMASH system. Calculation of P-T slopes and the change in Fe/(Fe+Mg) along reactions in the KFMASH system were made using the Gibbs method. In addition, the effect on the grid of MnO and CaO is evaluated quantitatively. The resulting grid is consistent with typical Buchan and Barrovian parageneses at medium to high grades. At low grades, the grid predicts an extensive stability field for the paragenesis chloritoid+biotite which arises because of the unusual facing of the reaction chloritoid+biotite + quartz+H₂O = garnet+chlorite+muscovite, which proceeds to the right with increasing T in the KFMASH system. However, the reaction proceeds to the left with increasing T in the MnKFASH system so the assemblage chloritoid + biotite is restricted to bulk compositions with high Fe/(Fe+Mg+Mn). Typical metapelites will therefore contain garnet+chlorite at low grades rather than chloritoid + biotite.     

The stability of sapphirine-quartz and hypersthene-sillimanite-quartz assemblages: an experimental investigation in the system FeO−MgO−Al2O3−SiO2 under H2O and CO2 conditions

Phase relations and mineral chemistry involving the phases garnet (Gt), spinel (Sp), hypersthene (Hy), sapphirine (Sa), cordierite (Cd), sillimanite (Sil) and quartz (Qz) have been experimentally determined in the system FMAS (FeO−MgO−Al₂O₂−SiO₂) under low fO₂ and for various H₂O/CO₂ conditions. Several compositions were studied with 100 (Mg/Mg+Fe) ratio ranging from 64 to 87 with excess quartz and sillimanite. Our data do not show any differences in Gt−Cd stability and composition as a function of H₂O, CO₂ and H₂O−CO₂ (±CH₄) content, in good agreement with a previous experimental study at lower temperature (Aranovich and Podlesskii 1983). At 1,000° C and 11 kbar, under CO₂-saturated conditions, cordierite grew from a crystalline mix unseeded with cordierite. Thus, under water-absent conditions, cordierite will have a high-P stability field in the presence of CO₂. If water has a pressure stabilizing effect on cordierite, then our results would indicate that the effects of H₂O and CO₂ are of the same magnitude at high temperature. Our data support the theoretical P-T grid proposed by Hensen (1986) for high-T metapelites and are largely consistent with the high-temperature experimental data of Hensen and Green (1973). The univariant boundary Gt+Cd=Hy+Sil+Qz, which marks the disappearance of Hy−Sil−Qz assemblages, has a negative dP/dT slope above 1,000° C and a positive one below this temperature. Extrapolation of our data to iron-free systems shows that the high-P breakdown limit of Mg-cordierite has a negative slope in the range 1,025–1,300° C and probably positive below 1,000° C. This indicates a maximum of stability for Mg-cordierite at around 1,000° C and 13 kbar. Because of the curvature of the univariant reactions En+Sil=Py+Qz, Mg−Cd=En+Sil+Qz and Gt+Cd=Hy+Sil+Qz, the iron-free invariant point involving the phases Py, En, Cd, Sil and Qz probably does not exist. Sapphirine—Qz-bearing assemblages are stable only at temperatures above 1,050° C. At 1,075° C, the joint Gt−Sa is stable up to 11 kbar. At higher pressure, garnet, sapphirine and quartz react according to the reaction Gt+Sa+Qz=Hy+Sil. Reequilibrated sapphirines are more aluminous than the theoretical endmember Mg₂Al₄SiO₁₀ due to AlAl=MgSi substitutions [100(Al₂O₃/Al₂O₃+FeO+MgO) in experimental sapphirines ranges from 50.5 to 52.2]. Sapphirine in the assemblage Sa−Cd−Sil−Qz shows a decrease in Al content with decreasing temperature and pressure, such that the alumina isopleths for sapphirine have a slight negative dP/dT slope. A similar decrease in Al content of sapphirine with temperature is also observed in Sa−Sil−Qz assemblages.     

Fe-Ti-oxides in metamorphic basites from the Eastern Alps,Austria: a contribution to the formation of solid solutions of natural Fe-Ti-oxide assemblages

The development of Fe-Ti oxide assemblages in basic rocks from the Penninic series of the southern Venediger rea, Austria, during polyphase Alpine metamorphism has been studied. Textural and compositional relations indicate thorough reequilibration of the opaque mineral assemblages during late Barrovian metamorphism at essentially static conditions of lower amphibolite to greenschist facies. In contrast, silicate mineralogy of the preceeding blueschist to eclogite facies metamorphism might still be preserved to a large extent. Chemical adjustment of the Fe-Ti oxide minerals to decreasing P-T conditions is characterized by (1) formation of complex intergrowths of ilmenite and hematite solid solutions (<550° C), (2) the decomposition of hemo-ilmenite ₁ to ferrianilmenite₂+magnetite+rutile and of ilmeno-hematite₁ to titanhematite₂+rutile±magnetite (<450° C), and (3) low-grade oxidation of ferrianilmenite₂ to magnetite+hematite-rutile intergrowths or hematite +rutile and of titanhematite₂ to hematite-rutile intergrowths (≦400° C). Chemical equilibrium is suggested by the regular partitioning of Cr, V, Mg and Mn between coexisting hemo-ilmenite, ilmeno-hematite, and magnetite. The hematite-ilmenite miscibility gap has been delimited on the basis of the bulk compositions of the exsolved phases and the temperature estimates obtained from Fe-Ti oxide thermometry.     

Amphibole zonation in metabasites as a guide to the evolution of metamorphic conditions

In low grade metabasites the amphibole components tremolite, glaucophane, edenite and tschermakite have their activities controlled by interactions with the excess components albite, clinozoisite, chlorite, quartz and H₂O vapor. Three types of reaction are involved, (i) Those in which only components of condensed phases take part: isopleths of equilibrium constant are straight lines in the P-T plane. (ii) Dehydration reactions in which entropy change due to change in Al coordination is of the same sign as that due to dehydration: isopleths of constant K are positive at low pressure and negative at high pressure. (iii) Dehydration reactions in which entropy change due to Al coordination change is opposite in sign to that of dehydration: isopleths of constant K loop in the P-T plane with positive slopes at low and at high pressure. Zonation in naturally occurring amphiboles records the evolution of metamorphic conditions in particular rocks. In an example from the eastern Alps (Austria) early conditions calculated as 15 kb, 200 ° C evolve upgrade to 6 kb, 525 ° C implying concurrent heating and erosion. The record of evolving conditions may span some 30 Ma of geological history.     

P--T Stabilities of Laumontite, Wairakite, Lawsonite, and Related Minerals in the System CaAl2Si2O8-SiO2-H2O

Hydrothermal investigation of the bulk composition CaO.Al₂O₃.4SiO₂+excessH₂O has been conducted using conventional techniques over thetemperature ranges 200–450 °C and 500–6000 barsP_fluid. A number of reactions have been studied by employingmineral mixtures consisting of reactants and products in about9: 1 and 1: 9 ratios. The phase relations were deduced fromrelatively long experiments by observing which seeded assemblagedisappeared or decreased markedly in one of the paired run charges. Laumontite was synthesized in the laboratory, probably for thefirst time. Laumontite was grown from seeded wairakite to over99 per cent using a weak NaCl solution. The refractive indicesof the synthetic material are about = 1.504 and = 1.514. Theaverage unit cell dimensions are a₀ = 14.761±0.005 Å;b₀ = 13.077±0.005 Å; c₀ = 7.561±0.003 Å;and ß = 112.02°±0.04°. Within the errorof measurement, the optical properties and cell parameters arein good agreement with those of natural laumontite. The equilibriumdehydration of laumontite involves two reactions: (1) laumontite= wairakite+2H₂O, passing through about 230 °C at 0.5 kb,255±5 °C at 1 kb, 282±5 °C at 2 kb, 297±5°C at 3 kb and 325±5 °C at 6 kb; and (2) laumontite= lawsonite+2 quartz+2H₂O, taking place at about 210 °Cat 3 kb and 275 °C at 3.2 kb. Above 300 °C, the equilibriumcurve for the solid-solid reaction (3) lawsonite+2 quartz =wairakite passes through 305 °C, 3.4 kb and 390 °C,4.4 kb. Equilibrium has been demonstrated unambiguously forthe above three reactions. The hydrothermal decomposition ofnatural laumontite above its own stability limit appears tobe a very slow process. Combined with previously published equilibria determined hydrothermallyfor wairakite, the phase relations are further investigatedby chemographic analysis interrelating the phases, laumontite,wairakite, lawsonite, anorthite, prehnite+kaolinite, and 2 pumpellyite+kaolinitein the system CaAl₂Si₂O₈-SiO₂-H₂O. This synthesis allowed theconstruction of a semiquantitative petrogenetic grid applicableto natural parageneses and the delineation of the physical conditionsfor the various low-grade metamorphic facies in low µCO₂environments. The similar stratigraphic zonations, consistentlyfound in a variety of environments, are recognized to be a functionof burial depth, geothermal gradient, and mineralogical andchemical composition of the parental rocks. Departures fromthe normal sequences are believed to be due to the combinationsof mineralogical variations, availability of H₂O, differencesin the ratio µCO₂/µH₂O, and the rate of reaction.The possible P-T boundaries for diagenesis, the zeolite facies,the lawsonite-albite facies, the prehnite-pumpellyite facies,and the adjacent metamorphic facies are illustrated diagrammatically.     

Experimental phase relations of hydrous peridotites modelled in the system H2O-CaO-MgO-Al2O3,-SiO2

The hydration of peridotites modelled by the system H₂O-CaO-MgO-Al₂O₃-SiO₂ has been treated theoretically after the method of Schreinemakers, and has been investigated experimentally in the temperature range 700°–900° C and in the pressure range of 8–14 kbar. In the presence of excess forsterite and water, the garnet- to spinel-peridotite transition boundary intersects the chlorite dehydration boundary at an invariant point situated at 865±5° C and 15.2±0.3 kbar. At lower pressures, a model spinel lherzolite hydrates to both chlorite- and amphibole-bearing assemblages at an invariant point located at 825±10° C and 9.3±0.5 kbar. At even lower pressures the spinel-to plagioclase-peridotite transition boundary intersects the dehydration curve for amphibole+forsterite at an invariant point estimated to lie at 855±10° C and 6.5±0.5 kbar.Both chlorite and amphibole were characterized along their respective dehydration curves. Chlorite was found to shift continuously from clinochlore, with increasing temperature, to more aluminous compositions. Amphibole was found to be tremolitic with a maximmum of 6 wt.% Al₂O₃.The experimentally determined curves in this study were combined with the determined or estimated stability curves for hydrous melting, plagioclase, talc, anthophyllite, and antigorite to obtain a petrogenetic grid applicable to peridotites, modelled by the system H₂O-CaO-MgO-Al₂O₃-SiO₂, that covers a wide range of geological conditions. Direct applications of this grid, although quite limited, can be made for ultramafic assemblages that have been extensively re-equilibrated at greenschist to amphibolite facies metamorphism and for some highgrade ultramafic assemblages that display clear signs of retrogressive metamorphism.     

Phase relations in the greenschist-blueschist-amphibolite-eclogite facies in the system Na2O-CaO-FeO-MgO-Al2O3-SiO2-H2O (NCFMASH), with application to metamorphic rocks from Samos, Greece

Calculated phase equilibria among the minerals sodic amphibole, calcic amphibole, garnet, chloritoid, talc, chlorite, paragonite, margarite, omphacite, plagioclase, carpholite, zoisite/clinozoisite, lawsonite, pyrophyllite, kyanite, sillimanite, quartz and H₂O are presented for the model system Na₂O-CaO-FeO-MgO-Al₂O₃-SiO₂-H₂O (NCFMASH), which is relevant for many greenschist, blueschist, amphibolite and eclogite facies rocks. Using the activity-composition relationships for multicomponent amphiboles constrained by Will and Powell (1992), equilibria containing coexisting calcic and sodic amphiboles could be determined. The blueschist–greenschist transition reaction in the NCFMASH system, for example, is defined by the univariant reaction sodic amphibole + zoisite = calcic amphibole + chlorite + paragonite + plagioclase (+ quartz + H₂O) occurring between approximately 420 and 450 °C at 9.5 to 10 kbar. The calculated petrogenetic grid is a valuable tool for reconstructing the PT-evolution of metabasic rocks. This is shown for rocks from the island of Samos, Greece. On the basis of mineral and whole rock analyses, PT-pseudosections were calculated and, together with the observed mineral assemblages and reaction textures, are used to reconstruct PT-paths. For rocks from northern Samos, pseudomorphs after lawsonite preserved in garnet, the assemblage sodic amphibole-garnet-paragonite-chlorite-zoisite-quartz and the retrograde appearance of albitic plagioclase and the formation of calcic amphibole around sodic amphibole constrain a clockwise PT-path that reaches its thermal maximum at some 520 °C and 19 kbar. The derived PT-trajectory indicates cooling during exhumation of the rocks and is similar to paths for rocks from the western part of the Attic-Cycladic crystalline complex. Rocks from eastern Samos indicate lower pressures and are probably related to high-pressure rocks from the Menderes Massif in western Turkey. Received: 8 July 1997 / Accepted: 11 February 1998     

Calculated greenschist facies mineral equilibria in the system CaO−FeO−MgO−Al2O3−SiO2−CO2−H2O

The kinetic problems associated with the experimental determination of reactions among complex solidsolution phases at low temperatures have hindered our understanding of the phase relations in greenschist facies rocks. In the absence of reliable experimental data, we have used the new, expanded internally-consistent thermodynamic dataset of Holland and Powell (1990), to present calculated phase equilibria for the system CaO–FeO–MgO–Al₂O₃–SiO₂–H₂O–CO₂ (CaFMASCH) with quartz in eccess, in the range 400°–500°C at low to intermediate pressures, involving the minerals amphibole, chlorite, anorthite, clinozoisite, dolomite, chloritoid, garnet, margarite, andalusite, and calcite. By solving independent sets of non-linear equations formed from equilibrium relationships, we calculate not only the loci of reactions in pressuretemperature-x(CO₂) space, but also the compositions of coexisting minerals in terms of the substitutions, FeMg_-1 and (Fe,Mg)SiAl_-1Al_-1. Invariant, univariant and divariant equilibria are calculated and discussed in relation to naturally-occurring greenschist facies metabasic and siliceous dolomitic mineral assemblages. We thus avoid the use of activity-corrected curves so commonly presented in the literature as a substitute for genuine univariant phase diagram boundaries.     

Thermodynamics and behavior of γ-Mg2SiO4 at high pressure: Implications for Mg2SiO4 phase equilibrium

Raman spectra of γ-Mg₂SiO₄ taken to 200 kbar were used to calculate entropy and heat capacity at various P-T conditions. These new thermodynamic data on γ-MgSiO₄, similar data on MgSiO₃ perovskite (pv), previous data on β-MgSiO₄ and MgO (mw), and previous volumetric data of all phases were used to calculate the phase boundaries in the Mg₂SiO₄ phase diagram. Our resulting slope for the β→γ transition (50±4 bar K^-1) is in excellent agreement with recent multi-anvil studies. The slopes for the β→pv+MgO and γ→pv+MgO are-7±3 and -25±4 bar K^-1, respectively, and are consistent with our CO₂ laser heated diamond anvil studies. These slopes result in a β-γ-MgO+pv triple point at approximately 229 kbar and 2260 K for the iron free system.     

Coupling thermodynamic modeling and high-resolution in situ LA-ICP-MS monazite geochronology: evidence for Barrovian metamorphism late in the Grenvillian history of southeastern Ontario

The Flinton Group is a greenschist to upper amphibolite facies package of metasediments in southeastern Ontario that was metamorphosed during the Ottawan Orogeny. Thermodynamic modeling of metapelitic mineral assemblages suggests an increase in peak conditions of metamorphism across the 40 km wide study area from 3.5 to 7.9 kbar and 540 to 715 °C. Garnet isopleth thermobarometry applied to the cores of compositionally zoned porphyroblasts reveals remarkably similar P-T conditions of initial crystallization at approximately 3.7–4.0 kbar and 512–520 °C, corresponding to a relatively high geothermal gradient of ca. 34–45 °C km^?1. It is inferred from modeling and reaction textures that metamorphism was along Barrovian P-T paths. Major and trace element zoning in garnet from one sample records a complex growth history as evidenced by major and trace element zoning and the distribution of xenotime, allanite and monazite inclusions. High-resolution (6 μm) LA-ICP-MS U-Pb geochronology performed on monazite in the rock matrix and included in the outer 150 μm of garnet rim-ward of a Y annulus revealed an age of 976?±?4 Ma. The age is interpreted to reflect monazite growth at the expense of allanite and apatite late in garnet’s growth history over the P-T interval 4.5–6.8 kbar and 540–640 °C. This new age estimate for near peak metamorphism fits well into the regional framework but is significantly younger than previously reported ages for Ottawan metamorphism. Based on microstructures this new age suggests that compressional tectonics were operating much later in the history of the Grenville of southeastern Ontario than previously thought.     

The heat capacity of a natural monticellite and phase equilibria in the system CaO-MgO-SiO2-CO2

The heat capacity of a natural monticellite (Ca_1.00Mg_.09Fe_.91Mn_.01Si_0.99O_3.99) measured between 9.6 and 343 K using intermittent-heating, adiabatic calorimetry yields C_p⁰(298) and S₂₉₈⁰ of 123.64 ± 0.18 and 109.44 ± 0.16 J · mol⁻¹K⁻¹ respectively. Extrapolation of this entropy value to end-member monticellite results in an S⁰₂₉₈ = 108.1 ± 0.2 J · mol⁻¹K⁻¹. High-temperature heat-capacity data were measured between 340–1000 K with a differential scanning calorimeter. The high-temperature data were combined with the 290–350 K adiabatic values, extrapolated to 1700 K, and integrated to yield the following entropy equation for end-member monticellite (298–1700 K): S_T⁰(J · mol⁻¹K⁻¹) = S₂₉₈⁰ + 164.79 In T + 15.337 · 10⁻³T + 22.791 · 10⁵T⁻² − 968.94. Phase equilibria in the CaO-MgO-SiO₂ system were calculated from 973 to 1673 K and 0 to 12 kbar with these new data combined with existing data for akermanite (Ak), diopside (Di), forsterite (Fo), merwinite (Me) and wollastonite (Wo). The location of the calculated reactions involving the phases Mo and Fo is affected by their mutual solid solution. A best fit of the thermodynamically generated curves to all experiments is made when the S⁰₂₉₈ of Me is 250.2 J · mol⁻¹ K⁻¹ less than the measured value of 253.2 J · mol⁻¹ K⁻¹.A best fit to the reversals for the solid-solid and decarbonation reactions in the CaO-MgO-SiO₂-CO₂ system was obtained with the ΔG⁰₂₉₈ (kJ · mole⁻¹) for the phases Ak(−3667), Di(−3025), Fo(−2051), Me(−4317) and Mo(−2133). The two invariant points − Wo and −Fo for the solid-solid reactions are located at 1008 ± 5 K and 6.3 ± 0.1 kbar, and 1361 ± 10 K and 10.2 ± 0.2 kbar respectively. The location of the thermodynamically generated curves is in excellent agreement with most experimental data on decarbonation equilibria involving these phases.     

Petrology and metamorphism of glaucophane eclogites in Changning-Menglian suture zone,Bangbing area,southeast Tibetan Plateau: An evidence for Paleo-Tethyan subduction

High/ultrahigh-pressure (HP/UHP) metamorphic complexes, such as eclogite and blueschist, are generally regarded as significant signature of paleo-subduction zones and paleo-suture zones. Glaucophane eclogites have been recently identified within the Lancang Group characterized by accretionary mélange in the Changning-Menglian suture zone, at Bangbing in the Shuangjiang area of southeastern Tibetan Plateau. The authors report the result of petrological, mineralogical and metamorphism investigations of these rocks, and discuss their tectonic implications. The eclogites are located within the Suyi blueschist belt and occur as tectonic lenses in coarse-grained garnet muscovite schists. The major mineral assemblage of the eclogites includes garnet, omphacite, glaucophane, phengite, clinozoisite and rutile. Eclogitic garnet contains numerous inclusions, such as omphacite, glaucophane, rutile, and quartz with radial cracks around. Glaucophane and clinozoisite in the matrix have apparent optical and compositional zonation. Four stages of metamorphic evolution can be determined: The prograde blueschist facies (M₁), the peak eclogite facies (M₂), the decompression blueschist facies (M₃) and retrograde greenschist facies (M₄). Using the Grt-Omp-Phn geothermobarometer, a peak eclogite facies metamorphic P-T condition of 3000–3270 MPa and 617–658°C was determined, which is typical of low-temperature ultrahigh-pressure metamorphism. The comparison of the geological characteristics of the Bangbing glaucophane eclogites and the Mengku lawsonite-bearing retrograde eclogites indicates that two suites of eclogites may have formed from significantly different depths or localities to create the tectonic mélange in a subduction channel during subduction of the Triassic Changning-Menglian Ocean. The discovery of the Bangbing glaucophane eclogites may represent a new oceanic HP/UHP metamorphic belt in the Changning-Menglian suture zone.©2021 China Geology Editorial Office.     

Experimental study of subsolidus phase relations and mixing properties of pyroxene in the system CaO-Al2O3-SiO2

Subsolidus phase relations in the system CaO-Al₂O₃-SiO₂ (CAS) were experimentally determined with tight reversals of several univariant curves and with 14 equilibration experiments containing the assemblage pyroxene + anorthite, where pyroxene is a binary solid solution of Ca-Tschermak (CaTs-CaAl₂SiO₆) and Ca-Eskola (CaEs-Ca_0.5AlSi₂O₆) endmembers.Reversals were obtained on the following reactions (bar, °C): 3An = Gr + 2Ky + Q (P = 22T ? 700), 3An + Cor = Gr + 3Ky (P = 21.8T ? 950), 3CaTs= Gr + 2Cor(P = 55T ? 53900), and 6CaTs_{(1 ? x)}CaEs_x = 2(1 ? 2x)Gr + 4(1 ? 2x)Cor + 9xAn. Observed slopes indicate 9.8 J/mol · K of Al-Si disorder in Ca-Tschermak pyroxene and 5.3 J/mol·K of Al-Si disorder in anorthite, at 1300°C. It is suggested that Al-Si disorder in anorthite increases by 1.9 J/mol · K from 700°C to 1300°C.Compositions of CaTs-CaEs pyroxene in equilibrium with anorthite and PbO-rich liquid were experimentally determined at 1400–1430°C and 22.7–30.8 kbar. Microprobe measurements gave compositions which are consistent with an ideal pyroxene solution and the following parameters for the reaction 3An = 2CaTs + 2CaEs (J, bar, K): 2RTln(X_CaTs · X_CaEs) + 60200 + 86.4T ? (5.06 + 13 × 10^?7P)P = 0, resulting in ΔH⁰_j = ?39.8 kJ/mol and S⁰ = 461.8 J/mol · K for the Ca-Eskola endmember at 1300°C. The obtained properties of the Ca-Eskola component are necessary for thermobarometry based on pyroxene bearing assemblages containing plagioclase, quartz, or kyanite.     

Template plotting

The template transformation allows rapid straight line plotting of fluid-bearing reactions for conditions exceeding the two phase field in water-carbon dioxide. The system CaO-MgO-SiO₂-H₂O-CO₂ with the six solid phase tale, tremolite, quartz, magnesite, calcite, and dolomite is treated at 2Kb. These phases, in T-NH₂O space at constant pressure, participate in two invariant points and eight reactions. Published T-NH₂O data for three of these reactions at 2Kb are transformed using the template into LnK-1/T space to yield the parameters of the best (by eye) straight lines. The parameters of the other five reaction lines are derived by linear combination of these parameters. The reactions belonging to one invariant point in T-NH₂O space will be straight lines that intersect at a single point in ln K-1/T space. This point is not the transformation of the invariant point, and it occurs because of the linear relationship between the reactions. Errors need to be considered in transforming the system back into T-NH₂O space, especially when reactions lines with shallow slopes are involved.