Long length scales of element transport during reaction texture development in orthoamphibole-cordierite gneiss: Thor-Odin dome,British Columbia,Canada

First-order factors controlling the textural and chemical evolution of metamorphic rocks are bulk composition and pressure–temperature–time (P–T–t) path. Although it is common to assume that major element bulk composition does not change during regional metamorphism, rocks with reaction textures such as corona structures record evidence for major changes in effective bulk composition (EBC) and therefore provide significant insight into the scale, pathways, and mechanisms of element transport during metamorphism. Quantifying changes in EBC is essential for petrologic applications such as calculation of phase diagrams (pseudosections). The progressive growth of complex corona structures on garnet and Al₂SiO₅ porphyroblasts in orthoamphibole-cordierite gneiss Thor-Odin dome (British Columbia, Canada) reduced the EBC volume of the rock during metamorphism and therefore had a dramatic effect on the evolution of the stable mineral assemblage. These rocks contain a chemical and textural record of metamorphic reactions and preserve 3D networks (reaction pathways) connecting corona structures. These coronal networks record long (>cm) length scales of localized element transport during metamorphism. P–T, T–X, and P–X pseudosections are used to investigate the control of effective bulk composition on phase assemblage evolution. Despite textural complexity and evidence for disequilibrium, mineral assemblages and compositions were successfully modeled and peak metamorphic conditions estimated at 750°C and 9 kbar. These results illustrate how textural and chemical changes during metamorphism can be evaluated using an integrated petrographic and pseudosection approach, highlight the importance of effective bulk composition choice for application of phase equilibria methods in metamorphic rocks, and show how corona structures can be used to understand the scale of compositional change and element transport during 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     

Fractionation of bulk rock composition due to porphyroblast growth: effects on eclogite facies mineral equilibria, Pam Peninsula, New Caledonia

Chemically zoned porphyroblasts in metamorphic rocks indicate that diffusional processes could not maintain equilibrium conditions on a grain scale during porphyroblast growth or establish it afterwards. An effect of this inability to maintain equilibrium is the progressive removal of elements forming garnet cores from any metamorphic reaction that occurs at the porphyroblast boundaries or in the matrix of the rock. To examine this effect on mineral assemblages, the Bence–Albee matrix correction was applied to X‐ray intensity maps collected using eclogite samples from northern New Caledonia in order to determine the chemical composition of all parts of the sample. The manipulation of these element maps allows a quantitative analysis of the fractionation of the bulk rock composition between garnet cores and the matrix. A series of calculated equilibrium‐volume compositions represents the change in matrix chemistry with progressive elemental fractionation as a consequence of prograde garnet growth under high‐P conditions. Pressure–temperature pseudosections are calculated for these compositions, in the CaO–Na₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O system. Assemblages, modal proportions and mineral textures observed in the New Caledonian eclogites can be closely modelled by progressively ‘removing’ elements forming garnet cores from the bulk rock composition. The pseudosections demonstrate how chemical fractionation effects the peak metamorphic assemblage, prograde textures and the development of retrograde assemblages.     

Pressure-temperature paths from P-T pseudosections and zoned garnets: potential, limitations and examples from the Zanskar Himalaya, NW India

The extraction of P-T histories from metamorphic rocks provides a valuable dataset for the elucidation of the tectonic mechanisms for orogeny. While continued re-equilibration frequently obliterates early information, garnet zonation and inclusion assemblages can often surmount this problem. The task is more difficult in high variance assemblages or if inclusions are not preserved, but one approach is to use pseudosections that are specific to the bulk composition of a given rock. In the latter case, the compositions and abundances of all the minerals are fixed at a given P-T point such that, if the effective bulk composition is known, the garnet composition alone can be used to reconstruct the history. Here, we explore this approach using examples from the Zanskar Himalaya, NW India. Pseudosections have been calculated for four pelitic to semipelitic rocks from the Zanskar Himalaya and have been contoured for garnet composition. The calculations adequately model the mineral assemblages in the rocks and predict the presence of chlorite in the early assemblage where chlorite is found as inclusions within garnet. Moreover, the pseudosections successfully model the garnet core compositions, with all three independent compositional contours overlapping at a single pressure and temperature. This occurs at ∼550 °C and at pressures varying from 3–7 kbar for the four rocks studied. We have been less successful, however, at modelling garnet compositions beyond the cores because fractionation of the effective bulk composition is caused by garnet growth itself. However, in this case, a combination of the␣pseudosection and conventional thermobarometry using␣Fe-Ti inclusions and matrix phases allows us to reconstruct␣the entire P-T history. The resulting P-T paths record burial of 3–5 kbar without significant temperature increase followed by isobaric heating of 50–100 °C. This evolution is consistent with Himalayan collision in the early Tertiary but a combination of the P-T data presented here and published geochronological data suggests renewed thrusting south of the suture zone in the Oligocene. In addition, the data demonstrate that no extra heat source is required to cause melting of the Himalayan crust in the Miocene. While melting could have occurred both by dehydration during decompression or in the presence of a fluid, the lack of garnet resorption does suggest decompression was rapid and followed quickly by cooling. This scenario favours melting by decompression. Received: 17 July 1997 / Accepted: 6 April 1998     

Orthopyroxene–sillimanite–quartz assemblages: distribution, petrology, quantitative P–T–X constraints and P–T paths

Granulite facies magnesian metapelites commonly preserve a wide array of mineral assemblages and reaction textures that are useful for deciphering the metamorphic evolution of a terrane. Quantitative pressure, temperature and bulk composition constraints on the development and preservation of characteristic peak granulite facies mineral assemblages such as orthopyroxene + sillimanite + quartz are assessed with reference to calculated phase diagrams. In NCKFMASH and its chemical subsystems, peak assemblages form mainly in high‐variance fields, and most mineral assemblage changes reflect multivariant equilibria. The rarity of orthopyroxene–sillimanite–quartz‐bearing assemblages in granulite facies rocks reflects the need for bulk rock X_Mg of greater than approximately 0.60–0.65, with pressures and temperatures exceeding c. 8 kbar and 850 °C, respectively. Cordierite coronas mantling peak minerals such as orthopyroxene, sillimanite and quartz have historically been used to infer isothermal decompression P–T paths in ultrahigh‐temperature granulite facies terranes. However, a potentially wide range of P–T paths from a given peak metamorphic condition facilitate retrograde cordierite growth after orthopyroxene + sillimanite + quartz, indicating that an individual mineral reaction texture is unable to uniquely define a P–T vector. Therefore, the interpretation of P–T paths in high‐grade rocks as isothermal decompression or isobaric cooling may be overly simplistic. Integration of quantitative data from different mineral reaction textures in rocks with varying bulk composition will provide the strongest constraints on a P–T path, and in turn on tectonic models derived from these paths.     

有效全岩成分的计算方法及其在变质相平衡模拟中的应用

变质相平衡模拟是变质岩领域近几十年最重要的进展之一，它已经成为确定变质作用P-T-t轨迹和探索变质演化过程的有力工具。变质岩的矿物组合不但与其形成的温度（T）和压力（P）条件有关，而且受控于岩石的全岩成分（X）。但是变质岩通常是不均匀的并且往往保留两期以上的矿物组合，因此计算不同成分域或不同变质演化期次的有效全岩成分是模拟P-T视剖面图的核心问题之一。在中-低温变质岩中，石榴石变斑晶的生长会不断地将其核部成分"冻结"而不参与后续变质反应，这导致根据实测全岩成分计算的P-T视剖面图无法有效地模拟石榴石幔部或边部生长阶段的变质演化过程。"瑞利分馏法"和"球体积法"利用电子探针实测的石榴石成分环带可以模拟计算石榴石各个生长阶段所对应的有效全岩成分，本文推荐使用这两个方法来处理石榴石变斑晶的分馏效应问题。相比较而言，石榴石在高温变质岩中通常无法保留生长阶段的成分环带特征，这是因为石榴石成分在高温条件下会发生扩散再平衡，并同时与多数基质矿物达到热力学平衡，这时一般不需要考虑石榴石的分馏效应。但是高温变质岩通常会发生部分熔融并伴随熔体的迁移，进而改变岩石的有效全岩成分。因此，通过P-T视剖面图模拟熔体迁移前后的变质演化过程需要使用"相平衡法"计算迁移的熔体成分以及熔体迁移前后岩石的有效全岩成分。此外，后成合晶与反应边是变质岩中最常见的退变质反应结构，但是后成合晶或反应边中的矿物之间并未达到热力学平衡。这种情况需要结合岩相学观察和矿物成分，利用最小二乘法确定后成合晶或反应边中发生的平衡反应方程式，进而获取变质反应发生时的有效全岩成分并通过计算P-T视剖面图来估算退变质的温压条件。除此之外，岩石体系中三价铁（Fe₂O₃）和H₂O含量的估算一直以来都是相平衡模拟研究中的难点，本文推荐使用P/T-X（Fe³⁺/Fe^tot，MH₂O）视剖面图来确定这两个组分的含量，这是因为P/T-X图可以估算各个变质演化阶段或特定矿物组合的Fe₂O₃或H₂O含量。     

An investigation of the relationship between bulk composition, inferred reaction progress and fluid-flow parameters for layered micaceous carbonates from Maine, USA

Reaction progress exhibited by multivariant assemblages in micaceous limestones can provide an excellent record of metamorphic fluid flow. However, it is necessary to understand the sensitivity of these assemblages to bulk‐composition parameters. Here, analysis of bulk composition on different scales and pseudosection construction are used to draw conclusions on relationships between bulk composition, fluid flow and reaction progress. Issues addressed include the effects of bulk composition on the mineralogical evolution of micaceous carbonates, the sensitivity of bulk composition to bulk‐composition sampling methods, the magnitude of cross‐layer fluid‐composition gradients, the potential for metasomatism to drive reaction progress, and the relative timing of reaction in adjacent layers. Pseudosections successfully represent observed mineral assemblages, constrain the position of reactions in T–X(CO₂) space, and allow assessment of the sensitivity of reaction position, inferred reaction progress and calculated fluid fluxes to uncertainties in bulk composition. The scale of bulk‐composition sampling affects bulk compositions, calculated modes, predicted mineral assemblages and calculated fluid compositions. Larger samples record an average of different lithological subdomains, while point‐count‐derived bulk compositions are subject to uncertainties related to the small number of sample points. The optimum bulk composition for pseudosection purposes probably lies between measured bulk compositions. Results suggest that reaction progress in some extensively reacted layers was driven by infiltration of H₂O‐rich fluid which flowed or diffused parallel to layering, perpendicular to layering in response to fluid‐composition gradients, and out of veins. Small variations in fluid composition across layering (ΔX(CO₂) < 0.02) were maintained by internal buffering by the mineral assemblages. Internal buffering must also have driven samples up a sequence of narrow low‐variance fields in T–X(CO₂) space, and so reaction in adjacent layers must have close to simultaneous. Metasomatic effects on reaction progress are likely to have been small, so long as the porosity was low.     

Calculated phase equilibria in K2O-FeO-MgO-Al2O3-SiO2-H2O for sapphirine-quartz-bearing mineral assemblages

Sapphirine, coexisting with quartz, is an indicator mineral for ultrahigh‐temperature metamorphism in aluminous rock compositions. Here a new activity‐composition model for sapphirine is combined with the internally consistent thermodynamic dataset used by THERMOCALC, for calculations primarily in K₂O‐FeO‐MgO‐Al₂O₃‐SiO₂‐H₂O (KFMASH). A discrepancy between published experimentally derived FMAS grids and our calculations is understood with reference to H₂O. Published FMAS grids effectively represent constant a_H2O sections, thereby limiting their detailed use for the interpretation of mineral reaction textures in compositions with differing H₂O. For the calculated KFMASH univariant reaction grid, sapphirine + quartz assemblages occur at P–T in excess of 6–7 kbar and 1005 °C. Sapphirine compositions and composition ranges are consistent with natural examples. However, as many univariant equilibria are typically not ‘seen’ by a specific bulk composition, the univariant reaction grid may reveal little about the detailed topology of multi‐variant equilibria, and therefore is of limited use for interpreting the P–T evolution of mineral assemblages and reaction sequences. Calculated pseudosections, which quantify bulk composition and multi‐variant equilibria, predict experimentally determined KFMASH mineral assemblages with consistent topology, and also indicate that sapphirine stabilizes at increasingly higher pressure and temperature as X_Mg increases. Although coexisting sapphirine and quartz can occur in relatively iron‐rich rocks if the bulk chemistry is sufficiently aluminous, the P–T window of stability shrinks with decreasing X_Mg. An array of mineral assemblages and mineral reaction sequences from natural sapphirine + quartz and other rocks from Enderby Land, Antarctica, are reproducible with calculated pseudosections. That consistent phase diagram calculations involving sapphirine can be performed allows for a more thorough assessment of the metamorphic evolution of high‐temperature granulite facies terranes than was previously possible. The establishment of a a‐x model for sapphirine provides the basis for expansion to larger, more geologically realistic chemical systems (e.g. involving Fe³⁺).     

Challenges in high-pressure granulite metamorphism in the era of pseudosections: reaction textures, compositional zoning and tectonic interpretation with examples from the Bohemian Massif

High‐pressure granulites are an important record of geodynamic processes in overthickened or subducted continental crust. Orthopyroxene‐free assemblages in granitic (ternary feldspar(s) + quartz + garnet + kyanite + rutile), intermediate (ternary feldspar(s) + quartz + garnet + clinopyroxene ± kyanite ±rutile ± titanite) and basic (garnet + clinopyroxene + plagioclase ± quartz + rutile) compositions indicate formation conditions at mantle depths. Clinopyroxene compositions in Variscan high‐pressure granulites are unusual in that they include omphacite (in plagioclase‐bearing rocks thus not eclogite) and Al‐rich diopside (i.e. indicating high Ca‐Tschermak content), with both yielding temperatures above 900 °C. Problems such as compositional zoning, multiple generations of key phases in reaction domains and unmixing of high‐temperature solid‐solution phases during cooling (ternary feldspars, omphacite) clearly indicate disequilibrium and require very careful interpretation as to which phases and compositions possibly represent a former equilibrium association. Pressure–temperature (P–T) determination by the pseudosection method, although allowing prediction of mineral assemblages, compositions and molar proportions for a fixed bulk composition for modelled P–T conditions, still requires reliable activity–composition information for the key phases feldspar and clinopyroxene as well as an interpretation of former equilibrium compositions in the investigated samples, i.e. the same restrictions applying to conventional thermobarometry. The interpretations of some recently determined pseudosections for the composition of Variscan clinopyroxene‐bearing high‐pressure granulites contradict numerous published P–T paths. However, quantitative information from thermobarometry or pseudosections must be integrated with key petrographic observations. In the case of the Variscan example, it is argued that petrographic observations and published P–T paths are consistent with mineral assemblages predicted in pseudosections and support existing tectonometamorphic models.     

Retrograde replacement of andalusite by Ca–Na mica in chloritoid-bearing metapelites. PTX modelling of rocks with different Al content in the MnNCKFMASH system

Andalusite porphyroblasts are totally pseudomorphosed by margarite–paragonite aggregates in aluminous pelites containing the peak mineral assemblage andalusite, chlorite, chloritoid, margarite, paragonite, quartz ± garnet, in a NW Iberia contact area. Equilibria at low P–T are investigated using new KFMASH and (mainly) MnCNKFMASH grids constructed with Thermocalc 3.21. P–T and T–X pseudosections with phase modal volume isopleths are constructed for compositions relatively richer and poorer in andalusite to model the assemblages in an andalusite‐bearing rock that contains a thin andalusite‐rich band (ARB) during retrogression. Their compositions, prior to retrogression, are used in the modelling, and have been retrieved by restoring the pseudomorph‐forming elements into the current‐depleted matrix, except for Al₂O₃ which is assumed to be immobile. Compositional differences between the thin band and the rest of the rock have not resulted in differences in andalusite porphyroblast retrogression. The absence of chloritoid resorbtion implies either a pressure increase at constant reacting‐system composition, or that its composition changed during retrogression at constant pressure, by becoming enriched in the progressively replaced andalusite porphyroblasts. T–X pseudosections at 1 kbar model this latter process using as end‐members in X, first, the restored original rock and ARB compositions, and, then the same process, taking into account the change in composition of both as retrogression proceeded. The MnNCKFMASH pseudosections of rocks with different Al contents facilitate making further deductions on the rock‐composition control of the resulting assemblages upon retrogression. Andalusite eventually disappears in relatively Al‐poor rocks, resulting, as in this study, in a rock formed by chloritoid–chlorite as the only FM minerals, plus margarite–paragonite pseudomorphs of andalusite. In rocks richer in Al, chlorite would progressively disappear and a kyanite/andalusite–chloritoid assemblage would eventually be stable at retrograde conditions. The Al‐silicate, stable during retrogression in Al‐rich rocks, indicates pressure conditions and hence the tectonic context under which retrogression took place.     

Calculated phase equilibria for MORB compositions: a reappraisal of the metamorphic evolution of lawsonite eclogite

Pseudosections calculated with thermocalc predict that lawsonite‐bearing assemblages, including lawsonite eclogite, will be common for subducted oceanic crust that experiences cool, fluid‐saturated conditions. For glaucophane–lawsonite eclogite facies conditions (500–600 °C and 18–28 kbar), MORB compositions are predicted in the NCKMnFMASHO system to contain glaucophane, garnet, omphacite, lawsonite, phengite and quartz, with chlorite at lower temperature and talc at higher temperature. In these assemblages, the pyrope content in garnet is mostly controlled by variations in temperature, and grossular content is strongly controlled by pressure. The silica content in phengite increases linearly with pressure. As the P–T conditions for these given isopleths are only subtly affected by common variations in bulk‐rock compositions, the P–T pseudosections potentially present a robust geothermobarometric method for natural glaucophane‐bearing eclogites. Thermobarometric results recovered both by isopleth and conventional approaches indicate that most natural glaucophane–lawsonite eclogites (Type‐L) and glaucophane–epidote eclogites (Type‐E) record similar peak P–T conditions within the lawsonite stability field. Decompression from conditions appropriate for lawsonite stability should result in epidote‐bearing assemblages through dehydration reactions controlled by lawsonite + omphacite = glaucophane + epidote + H₂O. Lawsonite and omphacite breakdown will be accompanied by the release of a large amount of bound fluid, such that eclogite assemblages are variably recrystallized to glaucophane‐rich blueschist. Calculated pseudosections indicate that eclogite assemblages form most readily in Ca‐rich rocks and blueschist assemblages most readily in Ca‐poor rocks. This distinction in bulk‐rock composition can account for the co‐existence of low‐T eclogite and blueschist in high‐pressure terranes.     

Petrology and isotopic evolution of granulites from central Madurai Block (southern India): reference to Ediacaran crustal evolution

The Madurai Block, constituting part of the southern granulite terrain in southern India, has contributed significantly towards understanding the UHT (ultrahigh-temperature) granulites that serve as a window into the mid-lower continental crust. The dominant rock types are charnockites, sapphirine-bearing granulites, garnet cordierite gneisses, and quartzites. Significant textural relations reveal multiphase reactions responsible for the formation of diverse mineral parageneses during prolonged metamorphic history of the area. Prograde reaction is evident from the textural relationship where biotite/sillimanite relics are seen as inclusion in garnet/orthopyroxene, suggesting dehydration reactions. The symplectitic assemblages that formed during isothermal decompression involve a series of cordierite-forming reactions, followed by retrogression and cooling. Variety of mineral assemblages present in the rocks of this area offer a wide spectrum of P–T sensors that provide details on the physical conditions of metamorphism. For the rigorous interpretation of the P–T path in the Perumalmalai area, quantitative phase diagrams (P–T pseudosections) have been constructed and contoured for the compositional as well as modal isopleths of involved mineral phases. The rocks of Perumalmalai area document a clockwise decompression P–T trajectory, consistent with crustal thickening followed by extensional collapse. SHRIMP U–Pb ages from zircon associated with sapphirine-bearing granulite facies rocks of Perumalmalai area suggest a widespread Ediacaran tectonothermal event. The occurrence of Ediacaran UHT metamorphism followed by isothermal decompression in the Madurai Block is consistent with the timing and physical conditions associated with the formation of East African Orogen during the amalgamation of Gondwana.     

Shooting at a moving target: phase equilibria modelling of high‐temperature metamorphism

Thermodynamic modelling and calculation of P–T pseudosections are commonly employed for quantifying the P–T evolution of metamorphic rocks. A key assumption involved in interpreting a P–T pseudosection is that the bulk‐rock composition used is representative of the effective bulk composition (EBC) from which apparently equilibrated mineral assemblages grew. Choosing an EBC can be difficult in cases where the rock has evolved significantly throughout the P–T history and has become domainal for whatever reason (e.g. loss of fluid and/or melt), particularly at suprasolidus conditions. During partial melting, melt migration may not only change the bulk composition by melt loss but also may generate local variations due to the variable consumption/loss of melt from domain to domain to create volumes of rock that were once internally equilibrated in the presence of a grain boundary melt, but which departed from equilibrium as inter‐granular mobility was slowed by local reductions in melt volume. As well as careful consideration of an EBC, the results of thermodynamic modelling are highly dependent on the specific thermodynamic data set and solution models used, as updates to these data sets may lead to substantially different calculated phase equilibria. This contribution addresses: (1) how consideration of evolving EBCs at multiple scales of observation can be used to resolve the history of complex high‐grade rocks, and (2) how use of different thermodynamic data sets and a–x models (i.e. thermocalc ds5.5 v. ds6) can result in different interpretations of metamorphic evolution. This study investigates the evolution of a mineralogically heterogeneous and texturally complex hand sample of granulite from the Gruf Complex (Central Alps). At the hand‐specimen scale, an EBC can be identified and used to constrain the P–T conditions at which the ‘whole rock’ was last in mutual equilibrium, in the presence of intergranular melt that has subsequently been lost or consumed. Smaller macrodomains (~cm scale) and microdomains (~mm scale) can be identified that represent subsequent evolution during and after melt channelization and loss, and P–T pseudosections can be calculated for the compositions of these domains. Using this approach reveals that the sample experienced a clockwise P–T path marked by near‐isothermal decompression following attainment of peak UHT conditions (~960 °C, 8.5 kbar). The approach enables construction of a P–T history of a rock for which P–T pseudosections are otherwise difficult to interpret. Thermodynamic modelling using ds6 yields similar results to those stated above, but suggests: (1) near‐isothermal decompression occurred over a wider pressure range (~0.5 kbar v. 1.5 kbar), and (2) that not all microdomains record this part of the P–T evolution.     

胶北地块高压与低压泥质麻粒岩的相平衡关系与p-T演化轨迹

利用最新的内洽性热力学数据库和THERMOCALC3.21程序对胶北地块高压与低压泥质麻粒岩的相平衡关系进行了定量分析。计算了胶北地块高压泥质麻粒岩、低压泥质麻粒岩和夕线石榴黑云片岩等代表性富铝岩石KFMASH(K2O-FeO-MgO-Al2O3-SiO2-H2O)体系的p-T视剖面图,再现了这些岩石随温压条件变化可能出现的各种矿物组合与矿物成分变化,发现原岩成分不同的变质岩石,尽管变质演化过程有所差异,但在麻粒岩相变质条件下所形成的矿物组合一致。通过计算泥质岩石在高压(p=1.0GPa)和低压(p=0.5GPa)条件下的T-X视剖面图,发现极度富铁、贫镁的岩石,在高压麻粒岩相条件下并不会生成含蓝晶石的特征矿物组合,在低压麻粒岩相条件下也不会生成含堇青石的特征矿物组合。将样品实际观测结果与p-T视剖面图的计算结果对比,确定胶北地块高压泥质麻粒岩变质峰期的温压条件为830~860℃,1.25~1.4GPa,峰期后呈现顺时针样式的p-T演化轨迹,反映陆壳先碰撞增厚、后又快速减薄的地质动力学过程;确定胶北地块低压泥质麻粒岩变质峰期的温压条件为790~820℃,0.62~0.68GPa,峰期后呈现近等压冷却的p-T演化轨迹。     

A sapphirine–phlogopite–cordierite paragenesis in a low-P amphibolite facies terrain, Arunta Inlier, Australia

Summary Silica-undersaturated phlogopite schists from the Cackleberry Metamorphics, Arunta Inlier, central Australia, preserve relatively low-temperature sapphirine-bearing parageneses that developed during low-pressure upper amphibolite facies metamorphism. Peak metamorphic phlogopite–cordierite–sapphirine assemblages are interpreted to have formed during the same event recorded in nearby metapelites, at c.3 kbar and 650–700 °C. Initial cooling of the terrain resulted in the breakdown of sapphirine to corundum–chlorite–phlogopite and corundum–spinel–chlorite assemblages. Further retrogression at greenschist facies conditions resulted in the replacement of sapphirine by diaspore–chlorite intergrowths. The reaction textures are consistent with a near-isobaric heating-cooling path at low-pressure, and provide evidence for the stability of sapphirine at c.700 °C at low pressures in rocks of an appropriate Mg- and Fe³⁺-rich bulk composition. Received August 15, 2001 accepted December 27, 2001     

Quantitative determination of metamorphic reaction history: mass balance relations between groundmass and mineral inclusion assemblages in metamorphic rocks

Qualitative and quantitative information about metamorphic reaction history and PT paths may be obtained from mineral inclusions in garnet by comparing the mineralogy, distribution, and compositions of paragenetically-related inclusions with minerals in the groundmass assemblage. Using the algebraic technique of singular value decomposition (SVD), we document mass balance relations between inclusion and groundmass assemblages in metapelitic rocks from two metamorphic terranes that experienced different peak metamorphic conditions, and whose transition from inclusion to groundmass assemblage records different PT path segments relative to peak conditions. We calculate mass balances relating an inclusion assemblage consisting in part of armored relics of chloritoid to groundmass mineral assemblages in a kyanite-staurolite mica schist from the Solitude Range, British Columbia, and an inclusion assemblage of kyanite, staurolite, and rutile to groundmass minerals in a sillimanite-cordierite gneiss from the Skagit Gneiss, North Cascade Range, Washington. Mass balances for each rock are consistent with reaction histories inferred from petrographic observations. In the Solitude Range schist, the results of mass balance calculations are consistent with the growth of staurolite and garnet at the expense of chloritoid during prograde metamorphism and suggest that chlorite, although not preserved as an inclusion, was involved in initial staurolite growth. In the Skagit sillimanite gneiss, mass balance relations exist between the inclusion suite, which formed during high pressure metamorphism, and the associated groundmass assemblage, which equilibrated at high temperature but much lower pressure. Mass balance does not exist between the groundmass of the Skagit sillimanite gneiss and the groundmass of a nearby kyanite-staurolite schist that has been proposed as a possible lower-grade equivalent of the sillimanite-bearing rocks. These results indicate that, although compositional modification and selective preservation of minerals must be taken into account, mineral inclusion suites may nevertheless preserve enough compositional information to allow reconstruction of complete or nearly complete pre-existing assemblages. This information may not be retrievable from any other source if no lower-grade equivalents of the rocks of interest are exposed.     

P–T–t path of the Hercynian low-pressure rocks from the Mandatoriccio complex (Sila Massif, Calabria, Italy): new insights for crustal evolution

The tectono‐metamorphic evolution of the Hercynian intermediate–upper crust outcropping in eastern Sila (Calabria, Italy) has been reconstructed, integrating microstructural analysis, P–T pseudosections, mineral isopleths and geochronological data. The studied rocks belong to a nearly complete crustal section that comprises granulite facies metamorphic rocks at the base and granitoids in the intermediate levels. Clockwise P–T paths have been constrained for metapelites of the basal level of the intermediate–upper crust (Umbriatico area). These rocks show noticeable porphyroblastic textures documenting the progressive change from medium‐P metamorphic assemblages (garnet‐ and staurolite‐bearing assemblages) towards low‐P/high‐T metamorphic assemblages (fibrolite‐ and cordierite‐bearing assemblages). Peak‐metamorphic conditions of ～590 °C and 0.35 GPa are estimated by integrating microstructural observations with P–T pseudosections calculated for bulk‐rock and reaction‐domain compositions. The top level of the intermediate–upper crust (Campana area) recorded only the major heating phase at low‐P (～550 °C and 0.25 GPa), as documented by the static growth of biotite spots and of cordierite and andalusite porphyroblasts in metapelites. In situ U–Th–Pb dating of monazite from schists containing low‐P/high‐T metamorphic assemblages gave a weighted mean U–Pb concordia age of 299 ± 3 Ma, which has been interpreted as the timing of peak metamorphism. In the framework of the whole Hercynian crustal section the peak of low‐P/high‐T metamorphism in the intermediate‐to‐upper crust took place concurrently with granulite facies metamorphism in the lower crust and with emplacement of the granitoids in the intermediate levels. In addition, decompression is a distinctive trait of the P–T evolution both in the lower and upper crust. It is proposed that post–collisional extension, together with exhumation, is the most suitable tectonic setting in which magmatic and metamorphic processes can be active simultaneously in different levels of the continental crust.     

榴辉岩中传统地质温压计新解：来自PT视剖面图的证据

石榴石-单斜辉石(GC)温度计和石榴石-单斜辉石-多硅白云母(GCP)压力计是确定榴辉岩形成温压条件的最常用方法,二者主要依据石榴石、绿辉石和多硅白云母中相组分之间的交换和转换变质反应.依据MORB成分计算的PT视剖面图表明,在不同榴辉岩矿物组合中,控制3个矿物相成分变化的相组分之间的变质反应不同.在低温含绿泥石、滑石和蓝闪石榴辉岩组合中,石榴石和绿辉石的镁含量主要受到含水矿物脱水反应的控制,并都随温度升高而升高,二者之间的铁镁交换反应并不起主要作用.因此,在自然界含有蓝闪石等含水矿物的低温榴辉岩中,由于绿辉石相对富镁而常常导致GC温度计结果偏低.在含有硬柱石的高压-超高压榴辉岩中,石榴石中的钙含量受到硬柱石的控制,随着压力升高或温度降低,硬柱石含量增加,使石榴石中钙降低,此时石榴石-绿辉石-多硅白云母之间的转换反应对石榴石成分的影响会很微弱,由于石榴石相对贫钙而导致GCP压力计结果偏低.在含有蓝晶石的中温高压-超高压榴辉岩中,矿物成分的变化受到石榴石-绿辉石之间的铁镁交换反应和石榴石．绿辉石．多硅白云母-蓝晶石-石英/柯石英之间的一系列转换反应控制,因此,GC和GCP温压计都能给出相对合理的结果.在低压普通角闪石榴辉岩中,石榴石和绿辉石中的镁含量主要反应压力变化,有时并不指示变质作用温度.在含有蓝闪石等含水矿物的低温榴辉岩中,Thermocalc程序中的平均温压(avPT)方法可以给出比较合适的温度,但压力结果与GCP压力计一样也会偏低一些.在蓝闪石和绿帘石等含水矿物消失后的中温蓝晶石榴辉岩中,avPT方法难以给出合理的PT信息.相对来说,视剖面图方法能够给出最多的PT信息,是目前确定变质岩PT条件的最好方法.     

The formation of eclogite facies metatroctolites and a general petrogenetic grid in Na2O–CaO–FeO–MgO–Al2O3–SiO2–H2O (NCFMASH)

Eclogite facies metatroctolites from a variety of Western Alps localities (Voltri, Monviso, Lanzo, Allalin, Zermat–Saas, etc.) that preserve textural evidence of their original form as bimineralic olivine‐plagioclase rocks are considered in terms of calculated mineral equilibria in the system Na₂O‐CaO‐FeO‐MgO‐Al₂O₃‐SiO₂‐H₂O (NCFMASH). Pseudosections, based on a new petrogenetic grid for NCFMASH presented here, are used to unravel the metamorphic history of the metatroctolites, considering the rocks to consist of different composition microdomains corresponding to the original olivine and plagioclase grains. On the basis that the preservation of the mineral assemblage in each microdomain will tend to be from where on a rock's P–T path the metamorphic fluid phase is used up via rehydration reactions, P–T pseudosections contoured for water content, and P–T path‐M_H2O (amount of water) pseudosections, are used to examine fluid behaviour in each microdomain. We show that the different microdomains are likely to preserve their mineral assemblages from different places on the P–T path. For the olivine microdomain, the diagnostic mineral assemblage is chloritoid + talc (+ garnet + omphacite). The preservation of this assemblage, in the light of the closed system P–T path‐M_H2O relationships, implies that the microdomain loses its metamorphic fluid as it starts to decompress, and, in the absence of subsequent hydration, the high pressure mineral assemblage is then preserved. In the plagioclase microdomain, the diagnostic assemblage is epidote (or zoisite) + kyanite + quartz suggesting a lower pressure (of about 2 GPa) than for the olivine microdomain. In the light of P–T path‐M_H2O relationships, development of this assemblage implies breakdown of lawsonite across the lawsonite breakdown reaction, regardless of the maximum pressure reached. It is likely that the plagioclase microdomain was mainly fluid‐absent prior to lawsonite breakdown, only becoming fluid‐present across the reaction, then immediately becoming fluid‐absent again.