Kyanite Sillimanite Research Articles

Bulk Compositional Influence on Diverse Metapelitic Mineral Assemblages in the Whetstone Lake Area, Ontario

Abstract Understanding the interplay between bulk composition and metamorphic grade underpins our interpretations of metamorphism in orogenic belts. The focus of this study is the regional garnet–staurolite–kyanite–sillimanite metamorphic sequence of the Whetstone Lake area, southeastern Ontario. In the kyanite and lower sillimanite zones of this area, there is exceptional diversity in metapelitic mineral assemblages that cannot be accounted for by differences in metamorphic grade. We present a data set of petrographic observations, phase proportions, whole-rock geochemical compositions, and mineral compositions, from thirty-two samples that encapsulate the range of assemblages found in these zones. Differences in bulk composition are the primary control on mineral assemblage development. Whole-rock XMg = molar MgO/(MgO + FeO) and ${\textrm{X}}_{\textrm{Fe}^{3+}}=\kern0.5em \textrm{molar}\ 2\times{\textrm{Fe}}_2{\textrm{O}}_3/\left(2\times{\textrm{Fe}}_2{\textrm{O}}_3+\textrm{FeO}\right)$ exert the greatest control on the observed mineral assemblages, whilst variation in MnO, K2O, and Al2O3 have a secondary influence. We use a set of quality factors ( Duesterhoeft & Lanari, 2020) to test the ability of thermodynamic models to reproduce the observed mineral assemblages, modal abundances, and mineral compositions in the diverse bulk compositions at Whetstone Lake. Eight samples were selected for phase equilibrium modelling, for which two bulk compositions were calculated for each sample: (1) a whole-rock bulk composition based on an X-ray fluorescence analysis and (2) a carefully considered local bulk composition based on combining mineral proportions with representative mineral compositions, as obtained from a single thin section. Our modelling uses thermodynamic data set 6.2 ( Holland & Powell, 2011) and the solution models of White et al. (2014a, 2014b) that incorporate several Fe3+ end members needed to model the natural data. Modelling in both types of bulk composition broadly predicted mineral assemblages that match those observed. In addition, predicted mineral assemblage fields overlap within uncertainty between 620°C and 675°C and between 6.5 and 7.5 kbar, consistent with the limited range of grade represented by the natural rocks. Predicted modal abundances better match those observed when phase diagrams are constructed using local bulk compositions compared to whole-rock bulk compositions. Despite the acceptable agreement between predicted and observed mineral assemblages, consistent discrepancies are found between predicted and observed mineral compositions. These include overestimation of ${\textrm{X}}_{\textrm{Mg}}^{\ast }$ in garnet, staurolite, and cordierite, overestimation of Ti in staurolite and biotite, underestimation of Si in biotite, and overestimation of Al and underestimation of Fe3+, Fe2+, and Mg in muscovite. The Whetstone Lake suite of this study will be useful to test the predictive capability of future thermodynamic models.

Advantages and Limitations of Combined Diffusion-Phase Equilibrium Modelling for Pressure–Temperature–Time History of Metamorphic Rocks

ABSTRACT This paper presents and discusses the results of phase diagram (Perple_X) and diffusion modelling (CZGM, or Compositional Zoning and its Modification by diffusion) to constrain the P–T path of metamorphism. The approach is based on the best fits between the zoning profile in measured garnet and that obtained by the intersections of garnet isopleths calculated by phase diagram modelling using whole rock bulk composition. The model was applied to garnets in natural rocks of various metamorphic grades, which were formed within different geotectonic environments. To compare the sequence of compositional change during Barrovian-type metamorphism, well-studied pelitic rocks from garnet–staurolite, kyanite–sillimanite, and sillimanite-K-feldspar metamorphic zones were selected. Garnets with two-stepped core and rim profiles that were formed during two different metamorphic stages or events were used for pressure–temperature (P–T) path constraint of each stage or event. For high-grade rocks, in which the original zoning profile in garnet was severely modified, the diffusion of the initial zoning profile was quantified to estimate the timescale of the metamorphic event. These rocks include high- to ultra-high-pressure rocks, which were subjected to thermal overprinting during collisional orogenesis. The results of the application of this approach allow for deciphering the reason why the calculated profile by phase diagram modelling does not fit with that of the measured garnet from low-grade rocks, in which garnet has preserved the original compositional zoning. This includes garnets whose nucleation was shifted from the garnet-in boundary to higher temperatures and pressures, as well as garnet crystallised during different metamorphic stages or events. Finally, the P–T paths in high-grade rocks were constrained after the multicomponent diffusion in garnet was quantified, and this was used for further P–T-time path constraint of metamorphism in the rocks.

The Maz Metasedimentary Series (Western Sierras Pampeanas, Argentina). A relict basin of the Columbia supercontinent?

AbstractThe Maz Metasedimentary Series is part of the Maz Complex that crops out in the sierras of Maz and Espinal (Western Sierras Pampeanas) and in the Sierra de Umango (Andean Frontal Cordillera), northwestern Argentina. The Maz Complex is found within a thrust stack of Silurian age, which later underwent open folding. The Maz Metasedimentary Series mainly consists of medium-grade garnet–staurolite–kyanite–sillimanite schists and quartzites, with minor amounts of marble and calc-silicate rocks. Transposed metadacite dykes have been recognized along with amphibolites, metagabbros, metadiorites and orthogneisses. Schist, quartzite and metadacite samples were analysed for SHRIMP U–Pb zircon dating. The Maz Metasedimentary Series is polymetamorphic and records probably three metamorphic events during the Grenvillian orogeny, at c. 1235, 1155 and 1035 Ma, and a younger metamorphism at c. 440–420 Ma resulting from reactivation during the Famatinian orogeny. The sedimentary protoliths were deposited between 1.86 and 1.33–1.26 Ga (the age of the Andean-type Grenvillian magmatism recorded in the Maz Complex), and probably before 1.75 Ga. The main source areas correspond to Palaeoproterozoic and, to a lesser magnitude, Meso-Neoarchaean rocks. The probable depositional age and the detrital zircon age pattern suggest that the Maz Metasedimentary Series was laid down in a basin of the Columbia supercontinent, mainly accreted between 2.1 and 1.8 Ga. The sedimentary sources were diverse, and we hypothesize that deposition took place before Columbia broke up. The Rio Apa block, and the Río de la Plata, Amazonia and proto-Kalahari cratons, which have nearby locations in the palaeogeographic reconstructions, were probably the main blocks that supplied sediments to this basin.

Geodynamic environments of the origin of poly- and monometamorphic complexes in the Southern Altai metamorphic belt, Central Asian orogenic belt

Tectonic sheets of various size along the southern slope of the Mongolian and Chinese Altai ranges and in eastern Kazakhstan include high-grade metamorphic rocks, which are collectively referred to as the Southern Altai Metamorphic Belt. Rocks of the sheets show traces of amphibolite-facies elevated-pressure metamorphism of the kyanite–sillimanite type M2. Some of the tectonic sheets display evidence of polymetamorphism: the rocks preserve textures and mineral assemblages of an earlier metamorphic episode (of elevated temperature and relatively low pressure) of the andalusite–sillimanite facies series M1. The earlier metamorphic episode occurred at 390–385 Ma, and the later one, at ~370–356 Ma. The protoliths of the high-grade metamorphic rocks were mostly Early Paleozoic terrigenous rocks and subordinate amounts of volcanic rocks analogous to the weakly metamorphosed or unmetamorphosed rocks in their northern surroundings. Typical rocks of the tectonic sheets are mafic dikes and massifs of the Gashun Nuur Complex, which were emplaced between metamorphic episodes M1 and M2. According to their geochemistry and Nd isotopic parameters, most of the metabasites are similar to enriched basalts of mid-oceanic ridges and oceans plateaus. The quantitatively subordinate group of the layered mafic bodies displays geochemical characteristics of subduction-related rocks. Correlations between the metamorphic events and magmatism in the continental (Mongolian and Chinese Altai) and paleoceanic (Trans-Altai Gobi and eastern Junggar) regions led us to suggest a geodynamic model for the development of the Southern Altai Metamorphic Belt. The volcano-terrigenous rocks, which were later metamorphosed, were accumulated mostly in the Early Paleozoic as an accretion wedge on an active continental margin. The earlier episode of high-temperature metamorphic M1 and coeval large-scale calc–alkaline magmatism occurred at the same active continental margin after the magmatic front shifted southward (in modern coordinates). The emplacement of the swarms of mafic bodies of the Gashun Nuur Complex and simultaneous rifting in the southern Chinese Altai were triggered by the subduction of an spreading ridge of an oceanic or backarc basin beneath the active margin. The second metamorphic episode (elevated-pressure metamorphism) M2 and overthrusting in the structures of the Altai are correlated with deformations at low angles and the transition from oceanic to continental volcanism in the Trans-Altai Gobi and Junggar. These tectonic processes were induced by the accretion of a system of mid-Paleozoic ensimatic island arcs of the Trans-Altai Gobi and Junggar to the Altai margin of the Siberian paleocontinent.

Geodynamic Environments of the Origin of Poly- and Monometamorphic Complexes in the Southern Altai Metamorphic Belt, Central Asian Orogenic Belt

Tectonic sheets of various size along the southern slope of the Mongolian and Chinese Altai ranges and in eastern Kazakhstan include high-grade metamorphic rocks, which are collectively referred to as the Southern Altai Metamorphic Belt. Rocks of the sheets show traces of amphibolite-facies elevated-pressure metamorphism of the kyanite–sillimanite type M2. Some of the tectonic sheets display evidence of polymetamorphism: the rocks preserve textures and mineral assemblages of an earlier metamorphic episode (of elevated temperature and relatively low pressure) of the andalusite–sillimanite facies series M1. The earlier metamorphic episode occurred at 390–385 Ma, and the later one, at ~370–356 Ma. The protoliths of the high-grade metamorphic rocks were mostly Early Paleozoic terrigenous rocks and subordinate amounts of volcanic rocks analogous to the weakly metamorphosed or unmetamorphosed rocks in their northern surroundings. Typical rocks of the tectonic sheets are mafic dikes and massifs of the Gashun Nuur Complex, which were emplaced between metamorphic episodes M1 and M2. According to their geochemistry and Nd isotopic parameters, most of the metabasites are similar to enriched basalts of mid-oceanic ridges and oceans plateaus. The quantitatively subordinate group of the layered mafic bodies displays geochemical characteristics of subduction-related rocks. Correlations between the metamorphic events and magmatism in the continental (Mongolian and Chinese Altai) and paleoceanic (Trans-Altai Gobi and eastern Junggar) regions led us to suggest a geodynamic model for the development of the Southern Altai Metamorphic Belt. The volcano-terrigenous rocks, which were later metamorphosed, were accumulated mostly in the Early Paleozoic as an accretion wedge on an active continental margin. The earlier episode of high-temperature metamorphism M1 and coeval large-scale calc–alkaline magmatism occurred at the same active continental margin after the magmatic front shifted southward (in modern coordinates). The emplacement of the swarms of mafic bodies of the Gashun Nuur Complex and simultaneous rifting in the southern Chinese Altai were triggered by the subduction of an spreading ridge of an oceanic or backarc basin beneath the active margin. The second metamorphic episode (elevated-pressure metamorphism) M2 and overthrusting in the structures of the Altai are correlated with deformations at low angles and the transition from oceanic to continental volcanism in the Trans-Altai Gobi and Junggar. These tectonic processes were induced by the accretion of a system of mid-Paleozoic ensimatic island arcs of the Trans-Altai Gobi and Junggar to the Altai margin of the Siberian paleocontinent.

Kinetic control of staurolite–Al2SiO5 mineral assemblages: Implications for Barrovian and Buchan metamorphism

AbstractThe distribution and textural features of staurolite–Al2SiO5 mineral assemblages do not agree with predictions of current equilibrium phase diagrams. In contrast to abundant examples of Barrovian staurolite–kyanite–sillimanite sequences and Buchan‐type staurolite–andalusite–sillimanite sequences, there are few examples of staurolite–sillimanite sequences with neither kyanite nor andalusite anywhere in the sequence, despite the wide (~2.5 kbar) pressure interval in which they are predicted. Textural features of staurolite–kyanite or staurolite–andalusite mineral assemblages commonly imply no reaction relationship between the two minerals, at odds with the predicted first development (in a prograde sense) of kyanite or andalusite at the expense of staurolite in current phase diagrams. In a number of prograde sequences, the incoming of staurolite and either kyanite, in Barrovian sequences, or andalusite, in Buchan‐type sequences, is coincident or nearly so, rather than kyanite or andalusite developing upgrade of a significant staurolite zone as predicted. The width of zones of coexisting staurolite and either kyanite, in Barrovian sequences, or andalusite, in Buchan‐type sequences, is much wider than predicted in equilibrium phase diagrams, and staurolite commonly persists upgrade until its demise in the sillimanite zone. We argue that disequilibrium processes provide the best explanation for these mismatches. We suggest that kyanite (or andalusite) may develop independently and approximately contemporaneously with staurolite by metastable chlorite‐consuming reactions that occur at lower P–T conditions than the thermodynamically predicted staurolite‐to‐kyanite/andalusite reaction, a process that involves only modest overstepping (<15°C) of the stable chlorite‐to‐staurolite reaction and which is favoured, in the case of kyanite, by advantageous nucleation kinetics. If so, the pressure difference between Barrovian kyanite‐bearing sequences and Buchan andalusite‐bearing sequences could be ~1 kbar or less, in better agreement with the natural record. The unusual width of coexistence of staurolite and Al2SiO5 minerals, in particular kyanite and andalusite, can be accounted for by a combination of lack of thermodynamic driving force for conversion of staurolite to kyanite or andalusite, sluggish dissolution of staurolite, and possibly the absence of a fluid phase to catalyse reaction. This study represents an example of how kinetic controls on metamorphic mineral assemblage development have to be considered in regional as well as contact metamorphism.

Mesoarchean convergent margin processes and crustal evolution: Petrologic, geochemical and zircon U–Pb and Lu–Hf data from the Mercara Suture Zone, southern India

The Mercara Shear Zone is sandwiched between the Western Dharwar Craton and the Coorg Block in the Southern Granulite Terrain of India, and is marked by steep gravity gradients interpreted to suggest the presence of underplated high-density material in the lower crust. Here we present geological, petrological and geochemical data, together with zircon U–Pb ages and Lu–Hf isotopes from a suite of metaigneous (TTG-related gneisses, charnockite, metagabbro, mafic granulite) and metasedimentary (quartz mica schist, khondalite, garnet biotite gneiss, kyanite–sillimanite bearing metapelite) rocks from this zone. Geochemical data on the magmatic suite suggests formation through subduction-related arc magmatism, whereas the metasediments represent volcano-sedimentary trench sequences. Phase equilibrium modeling of mafic granulites from the Mercara Shear Zone suggests P–T range of 10–12kbar at 700°C to 900°C. The zircon data yield weighted mean 207Pb/206Pb ages of 3229±80Ma for metagabbro, 3168±25Ma for the charnockite, and 3181±20Ma for the mafic granulite. Ages ranging from 3248±28Ma to 3506±26Ma were obtained from zircons in the kyanite/sillimanite bearing metapelite, 3335±44Ma from khondalite, 3135±14Ma from garnet biotite gneiss, 3145±17Ma to 3292±57Ma from quartz mica schist and 3153±15Ma to 3252±36 from TTG gneiss. The tightly defined ages of 3.1 to 3.2Ga from igneous zircons in the magmatic suite suggest prominent Mesoarchean convergent margin magmatism. The timing of high grade metamorphism as constrained from metamorphic overgrowths in zircons is ca. 3.0Ga which might mark the collisional event between the Western Dharwar Craton and the Coorg Block. Hf isotope features suggest magma derivation mostly from juvenile sources and the Lu–Hf model ages indicate that the crust building might have also involved partial recycling of basement rocks as old as ca. 3.8Ga. Our study defines the Mercara Shear Zone as a terrane boundary, and possible Mesoarchean suture along which the Coorg Block was accreted to the Western Dharwar Craton. The accretion of these continental fragments might have coincided with the birth of the oldest supercontinent “Ur”.

Reconstruction of the metamorphic P–T path from the garnet zoning in aluminous schists from the Tsogt Block, Mongolian Altai

The paper presents original authors’ data on aluminous schists in the Tsogt tectonic plate in the Southern Altai Metamorphic Belt. The nappe includes a medium-temperature/medium-pressure zonal metamorphic complex, whose metamorphic grade varies from the greenschist to epidote-amphibolite facies. The garnet and garnet–staurolite schists contain three garnet generations of different composition and morphology. The P–T metamorphic parameters estimated by mineralogical geothermometers and geobarometers and by numerical modeling with the PERPLEX 668 software provide evidence of two successive metamorphic episodes: high-gradient (of the andalusite–sillimanite type, geothermal gradient approximately 40–50°/km) and low-gradient (kyanite–sillimanite type, geothermal gradient approximately 27°/km). The P-T parameters of the older episode are T = 545–575°C and P = 3.1–3.7 kbar. Metamorphism during the younger episode was zonal, and its peak parameters were T = 560–565°C, P = 6.4–7.2 kbar for the garnet zone and T = 585–615°C, P = 7.1–7.8 kbar for the staurolite zone. The metamorphism evolved according to a clockwise P–T path: the pressure increased during the first episode at a practically constant temperature, and then during the second episode, the temperature increased at a nearly constant pressure. Such trends are typical of metamorphism related to collisional tectonic settings and may be explained by crustal thickening due to overthrusting. The regional crustal thickening reached at least 15–18 km.

Mesoproterozoic suturing of Archean crustal blocks in western peninsular India: Implications for India–Madagascar correlations

The Kumta and Mercara suture zones welding together Archean crustal blocks in western peninsular India offer critical insights into Precambrian continental juxtapositions and the crustal evolution of eastern Gondwana. Here we present the results from an integrated study of the structure, geology, petrology, mineral chemistry, metamorphic P–T conditions, zircon U–Pb ages and Lu–Hf isotopes of metasedimentary rocks from the two sutures. The dominant rocks in the Kumta suture are greenschist- to amphibolite-facies quartz–phengite schist, garnet–biotite schist, chlorite schist, fuchsite schist and marble. The textural relations, mineral chemistry and thermodynamic modelling of garnet–biotite schist from the Kumta suture indicate peak metamorphic P–T conditions of ca. 11kbar at 790°C, with detrital SHRIMP U–Pb zircon ages ranging from 3420 to 2547Ma, εHf (t) values from −9.2 to 5.6, and TDMc model ages from 3747 to 2792Ma. The K–Ar age of phengite from quartz–phengite schist is ca. 1326Ma and that of biotite from garnet–biotite schist is ca. 1385Ma, which are interpreted to broadly constrain the timing of metamorphism related to the suturing event. The Mercara suture contains amphibolite- to granulite-facies mylonitic quartzo-feldspathic gneiss, garnet–kyanite–sillimanite gneiss, garnet–biotite–kyanite–gedrite–cordierite gneiss, garnet–biotite–hornblende gneiss, calc-silicate granulite and metagabbro. The textural relations, mineral chemistry and thermodynamic modelling of garnet–biotite–kyanite–gedrite–cordierite gneiss from the Mercara suture indicate peak metamorphic P–T conditions of ca. 13kbar at 825°C, followed by isothermal decompression and cooling. For pelitic gneisses from the Mercara suture, LA-ICP-MS U–Pb zircon ages vary from 3249 to 3045Ma, εHf (t) values range from −18.9 to 4.2, and TDMc model ages vary from 4094 to 3314Ma. The lower intercept age of detrital zircons in the pelitic gneisses from the Mercara suture ranges from 1464 to 1106Ma, indicating the approximate timing of a major lead-loss event, possibly corresponding to metamorphism, and is broadly coeval with events in the Kumta suture. Synthesis of the above results indicates that the Kumta and Mercara suture zones incorporated sediments from Palaeoarchean to Mesoproterozoic sources and underwent high-pressure metamorphism in the late Mesoproterozoic. The protolith sediments were derived from regions containing juvenile Palaeoarchean crust, together with detritus from the recycling of older continental crust. Integration of the above results with published data suggests that the Mesoproterozoic (1460–1100Ma) Kumta and Mercara suture zones separate the Archean (3400–2500Ma) Karwar–Coorg block and Dharwar Craton in western peninsular India. Based on regional structural and other geological data we interpret the Kumta and Mercara suture zones as extensions of the Betsimisaraka suture of eastern Madagascar into western India.

Evolution of high-pressure mafic granulites and pelitic gneisses from NE Madagascar: Tectonic implications

The occurrence of high-pressure mafic–ultramafic bodies within major shear zones is one of the indicators of paleo-subduction. In mafic granulites of the Andriamena complex (north-eastern Madagascar) we document unusual textures including garnet–clinopyroxene–quartz coronas that formed after the breakdown of orthopyroxene–plagioclase–ilmenite. Textural evidence and isochemical phase diagram calculations in the Na2O–CaO–K2O–FeO–MgO–Al2O3–SiO2–H2O–TiO2 system indicate a pressure–temperature (P–T) evolution from an isothermal (780°C) pressure up to c. 24kbar to decompression and cooling. Such a P–T trajectory is typically attained in a subduction zone setting where a gabbroic/ultramafic complex is subducted and later exhumed to the present crustal level during oceanic closure and final continental collision. The present results suggest that the presence of such deeply subducted rocks of the Andriamena complex is related to formation of the Betsimisaraka suture. LA-ICPMS U–Pb zircon dating of pelitic gneisses from the Betsimisaraka suture yields low Th/U ratios and protolith ages ranging from 2535 to 2625Ma. A granitic gneiss from the Alaotra complex yields a zircon crystallization age of ca. 818Ma and Th/U ratios vary from 1.08 to 2.09. K–Ar dating of muscovite and biotite from biotite–kyanite–sillimanite gneiss and garnet–biotite gneiss yields age of 486±9Ma and 459±9Ma respectively. We have estimated regional crustal thicknesses in NE Madagascar using a flexural inversion technique, which indicates the presence of an anomalously thick crust (c. 43km) beneath the Antananarivo block. This result is consistent with the present concept that subduction beneath the Antananarivo block resulted in a more competent and thicker crust. The textural data, thermodynamic model, and geophysical evidence together provide a new insight to the subduction history, crustal thickening and evolution of the high-pressure Andriamena complex and its link to the terminal formation of the Betsimisaraka suture in north-eastern Madagascar.

Monazite Dating of Prograde and Retrograde P–T–d paths in the Barrovian terrane of the Thaya window, Bohemian Massif

Monazite laser ablation–split-stream inductively coupled plasma–mass spectrometry (LASS) was used to date monazite in situ in Barrovian-type micaschists of the Moravian zone in the Thaya window, Bohemian Massif. Petrography and garnet zoning combined with pseudosection modelling show that rocks from staurolite–chlorite, staurolite, kyanite and kyanite–sillimanite zones record burial in the S1 fabric under a moderate geothermal gradient from 4–4·5 kbar and ∼530–540°C to 5 kbar and 570°C, 6–7 kbar and 600–640°C, 7·5–8 kbar and 630–650°C, and 8 kbar and 650°C, respectively. In the kyanite and kyanite–sillimanite zones, garnet rim chemistry and local syntectonic replacement of garnet by sillimanite–biotite aggregates point to re-equilibration at 5·5–6 kbar and 630–650°C in the S2 fabric. Heterogeneously developed retrograde shear zones (S3) are marked by widespread chloritization, but minor chlorite is present in the studied samples. Monazite abundance and size increase with metamorphic grade from 5 µm in the staurolite–chlorite zone to >100 µm in the kyanite and kyanite–sillimanite zones. Irrespective of the monazite-forming reaction, this is interpreted as the onset of limited prograde monazite growth at staurolite grade, and continued prograde monazite growth after the kyanite-in reaction, compatible with conditions of about 5·5 kbar and 570°C and 7·5 kbar and 630°C from pseudosection modelling. Monazite is zoned, showing embayments and sharp boundaries between zones, with low Y in the staurolite zone, high-Y cores and low-Y rims in the kyanite zone, and high-Y cores, a low-Y mantle and a high-Y rim in the sillimanite zone. The 207Pb-corrected 238U/206Pb ages from three samples range from 344 ± 7 to 330 ± 7 Ma, irrespective of metamorphic grade. The dates from monazite inclusions are interpreted as the ages of the staurolite- and kyanite-in reactions along the prograde path at 340 and 337 ± 7 Ma, respectively. The monazite in the matrix (and some inclusions) is interpreted as dating the prograde crystallization at (340–337) ± 7 Ma within the S1 fabric, and then being affected by recrystallization at or down to 332 ± 7 Ma in the S2 and S3 fabrics. The two groups of data, for 340–337 and 332 Ma, are significantly different when only their in-run uncertainties (±1–3 Myr) are compared and indicate a 9 ± 3 Myr period of monazite (re)crystallization. A systematic increase in heavy rare earth element (HREE) content with decreasing monazite age from 344 to 335 Ma is correlated with growth on the prograde P–T path; the drop in HREE of monazite at 335–328 Ma is assigned to recrystallization. The presence of chlorite even in the least retrogressed samples witnesses limited external fluid availability on the retrograde P–T path. Migration of this fluid was probably responsible for heterogeneous fluid-assisted recrystallization and resetting of original prograde monazite, even where included in garnet, staurolite or kyanite. It is suggested that the rocks passed the chlorite-in reaction on the retrograde path at 332 ± 7 Ma. The timing of burial in the Thaya window, a deformed part of the underthrust Brunia microcontinent, was coeval with exhumation of granulites and migmatites of the Moldanubian orogenic root at c. 340 Ma.

A new perspective on metamorphism and metamorphic belts

The discovery of ultrahigh-pressure rocks from collision-type orogenic belts has revolutionized the classic interpretation of (1) progressive and retrogressive metamorphism recorded on surface exposures of regional metamorphic belts, (2) geochronology of the various stages of metamorphism, (3) origin of metamorphic textures, (4) P– T– t path, (5) metamorphic facies series, (6) exhumation model, and (7) role of fluids during regional metamorphism. Based mainly on our recent studies of the Kokchetav, Dabie Shan, Indonesia, Franciscan and Sanbagawa belts, we suggest the following revolutionary paradigm shifts. The so-called mineral isograds defined on the maps of regional metamorphic belts were a misunderstanding of the progressive dehydration reaction during subduction because extensive late-stage hydration has mostly obliterated the progressive minerals in pelitic–psammitic and metabasic rocks. Progressively zoned garnet has survived as the sole progressive mineral that was unstable with the majority of matrix-forming minerals. The classic Barrovian isograds should therefore be carefully re-examined. The well-documented SHRIMP chronology of spot-dating zoned zircons with index minerals from low- P in the core, through HP–UHP in the mantle to low- P on the rim clearly shows that the slow exhumation speed of 23–40 My from mantle depth to mid-crustal level was followed by mountain building with doming at latest stage. Extensive hydration of the UHP–HP unit occurred due to fluid infiltration underneath, when the UHP–HP unit intruded the low-grade to unmetamorphosed unit at a mid-crustal level. Most deformation textures such as mineral lineations, porphyroblasts, pull-apart, or boudinaged amphiboles, formed during extensive hydration at the late stage do not constrain the progressive stress regime. The P– T–time path determined by thermobarometry using mineral inclusions in garnet and forward modeling of garnet zoning, independent of the matrix minerals, indicates an anticlockwise trend in the P– T space, and follows an independent P– T change in the metamorphic facies series. This is consistent with the numerically calculated geotherm along the Wadati–Benioff plane. Collision-type orogenic belts have long been regarded as being characterized by the intermediate-pressure type metamorphic facies series. The kyanite–sillimanite is an apparent type facies series formed by the late-stage extensive hydration. In contrast, the original high- P to ultrahigh- P type facies series with an anticlockwise kink-point at around 10 kb is a progressive type. A collision-type regional metamorphic belt crops out as a very thin unit sandwiched between overlying and underlying low- P or weakly metamorphosed units. The metamorphic belt has an aspect ratio (thickness vs width) of 1:100, and it extends for several hundreds to a thousand km. It resembles a thin mylonitic intrusion from the mantle extending from 100 to 200 km depth into the crustal rock unit. The underlying unit is thermally metamorphosed in the andalusite–sillimanite type facies series. The major reason for the misunderstanding of the progressive metamorphism in collision-type orogenic belts is the underestimation of the role of fluids derived from the underlying low-grade metamorphic unit, when juxtaposed at a mid-crustal level. The circulation of fluids along the plate boundaries is more important than a P– T change.

Kyanite–sillimanite metamorphism of the Precambrian complexes, Transangarian region of the Yenisei Ridge

Four Precambrian metamorphic complexes in the vicinity of regional faults in the Transangarian region of the Yenisei Ridge were examined. Based on geothermobarometry and P-T path calculations, our geological and petrological studies showed that the Neoproterozoic medium-pressure metamorphism of the kyanite–sillimanite type overprinted regionally metamorphosed low-pressure andalusite-bearing rocks at about 850 Ma. A positive correlation between rock ages and P-T estimates for the kyanite-sillimanite metamorphism provide evidence of the regional structural and tectonic heterogeneity. The medium-pressure metamorphism was characterized by (1) the development of deformational structures and textures, and kyanite-bearing blastocataclasites (blastomylonites) with sillimanite, garnet, and staurolite after andalusite-bearing regional metamorphic rocks; (2) insignificant apparent thickness of the zone of medium-pressure zonal metamorphism (from 2.5 to 7 km), which was localized in the vicinity of the overthrusts; (3) a low metamorphic field gradient during metamorphism (from 1–7 to 12 °C/km); and (4) a gradual increase in lithostatic pressure towards the thrust faults. These specific features are typical of collisional metamorphism during overthrusting of continental blocks and are evidence for near-isothermal loading. This event was justified within the framework of the crustal tectonic thickening model via rapid overthrusting and subsequent rapid uplifting and erosion. The results obtained allowed us to consider medium-pressure kyanite-bearing metapelites as a product of collision metamorphism, formed either by unidirectional thrusting of rock blocks from Siberian craton onto the Yenisei Ridge in the zones of regional faults (Angara, Mayakon, and Chapa areas) or by opposite movements in the zone of splay faults of higher orders (Garevka area).

GEOMETRY OF ISOGRADIC, ISOTHERMAL, AND ISOBARIC SURFACES: INTERPRETATION AND APPLICATION

An isograd represents a line on a map resulting from the intersection of an isogradic surface with the topography. It is inferred to represent a metamorphic reaction and is inherently diachronous because of the time required to move heat, fluid, and reactions through the rocks. Relief of 500 m or more or reconstruction of post-metamorphic faulting and folding will be essential to estimate the orientation of an isogradic surface. Ductile rocks require that isobaric surfaces dip gently during metamorphism, and they must do so in restored sections. Isotherms can be constructed into cross-sections because the angle between the isograds and isotherms are a function of the P–T slope of univariant reaction curves. Near Mica Creek, British Columbia, we can infer an 1100-bar difference in pressure along a kyanite–sillimanite (ky–sil) isograd, based upon structural relief. The experimentally determined ky–sil P–T curve suggests a P/T slope of ~20 bars/°C. We calculate the apparent variation in P and T along the isograd and reconstruct the isotherm geometry. Dehydration-reaction curves, at relatively high P, have steep P–T slopes, and the associated isogradic surfaces can approximate isothermal surfaces. Some dehydration isograds are “smeared out”, but for regional scales, these latter isogradic surfaces at high P still approximate the isothermal surface geometry. Data on isogradic surfaces can be used to set limits on parts of the P–T grids for pelitic rocks. Locally, bathozone-defining isobaric surfaces can be reconstructed; we suggest a ΔP of 2 kbar for the garnet – biotite – kyanite bathozone.

Polymetamorphism in the NE Shackleton Range, Antarctica: Constraints from Petrology and U-Pb, Sm-Nd, Rb-Sr TIMS and in situ U-Pb LA-PIMMS Dating

Metapelitic rock samples from the NE Shackleton Range, Antarctica, include garnet with contrasting zonation patterns and two age spectra. Garnet porphyroblasts in K-rich kyanite–sillimanite– staurolite–garnet–muscovite–biotite schists from Lord Nunatak show prograde growth zonation, and give Sm–Nd garnet, U–Pb monazite and Rb–Sr muscovite ages of 518 ± 5, 514 ± 1 and 499 ± 12 Ma, respectively. Geothermobarometry and P–T pseudo-section calculations in the model system CaO–Na2O–K2O– TiO2–MnO–FeO–MgO–Al2O3–SiO2–H2O are consistent with garnet growth during prograde heating from 540°C/7 kbar to 650°C/7·5 kbar, and partial resorption during a subsequent P–T decrease to <650°C at <6 kbar. All data indicate that rocks from Lord Nunatak were affected by a single orogenic cycle. In contrast, garnet porphyroblasts in K-poor kyanite–sillimanite– staurolite–garnet–cordierite–biotite-schists from Meade Nunatak show two growth stages and diffusion-controlled zonation. Two distinct age groups were obtained. Laser ablation plasma ionization multicollector mass spectrometry in situ analyses of monazite, completely enclosed by a first garnet generation, yield ages of c. 1700 Ma, whereas monazite grains in open garnet fractures and in most matrix domains give c. 500 Ma. Both age groups are also obtained by U–Pb thermal ionization mass spectrometry analyses of matrix monazite and zircon, which fall on a discordia with lower and upper intercepts at 502 ± 1 and 1686 ± 2 Ma, respectively. Sm–Nd garnet dating yields an age of 1571 ± 40 Ma and Rb–Sr biotite analyses give an age of 504 ± 1 Ma. Integrated geochronological and petrological data provide evidence that rocks from Meade Nunatak underwent a polymetamorphic Barrovian-type metamorphism: (1) garnet 1 growth and subsequent diffusive garnet annealing between 1700 and 1570 Ma; (2) garnet 2 growth during the Ross Orogeny at c. 500 Ma. During the final orogenic event the rocks experienced peak P–T conditions of about 650°C/7·0 kbar and a retrograde stage at c. 575°C/4·0 kbar.

Optimization of some key geothermobarometers for pelitic metamorphic rocks

AbstractI will consider mainly geothermobarometry in medium-grade pelitic rocks, including the garnet-biotite (GB) geothermometer, the Grossular-Al silicate-plagioclase (GASP) geobarometer, and the muscovite-almandine-biotite-sillimanite (MABS) geobarometer. For GB (Holdaway, 2000) experimental data and estimated biotite ΔWTi were used to optimize two exchange parameters and four biotite Margules parameters. Using stepwise linear regression, experimental vs. calculated T were constrained to lie on a line with slope of one and intercept of zero, maximizing r2. The best model involves experiments by Ferry and Spear (1978) and Perchuk and Lavrent’eva (1983), suggesting minimal viAl in the Ferry and Spear product biotite. For GASP (Holdaway, 2001), end-member experimental data do not adequately constrain the equilibrium. I used the GB model above, and allowed the end-member curve to rotate about the best-constrained part of the GASP end-member data. The end-member curve was further constrained with the kyanite-sillimanite (K-S) boundary using published chemical data on 76 pelitic schist samples from 11 localities, rejecting Low-Grs and low-An samples. The Fuhrman and Lindsley (1988) plagioclase model gives the best results. For MABS, work in progress involves 61 samples from the 11 localities which have muscovite analyses. Biotite Margules parameters were based on the GB model and McMullin et al. (1991). The MABS end-member curve was calibrated by comparison of P values determined using trial MABS data and GASP results. The P values for the 61 samples agree well with the K-S boundary, and sillimanite-bearing rocks of west-central Maine all fall in the sillimanite field. Preliminary biotite values are: GAnn = –5149198 – 412.05 T, WAlFe = –14023 + 28.14 T, WAlMg = –259582 + 308.44 T, WTiFe = 124842 – 98.67 T, WTiMg = –186148 + 271.72 T. For geobarometry, the Berman (1988, revised 1992) database was used with adjustable H and S of grossular for GASP and H and S of annite for MABS. The accuracy of currently available databases, activity models and mole fraction models is not adequate for good geothermobarometry, without further refinement. Adjustable parameters tend to compensate for error in activity models, mole fraction models and databases.

Pro‐ and retrogradeP–Tevolution of granulites of the Beit Bridge Complex (Limpopo Belt, South Africa): constraints from quantitative phase diagrams and geotectonic implications

AbstractInterpretations based on quantitative phase diagrams in the system CaO–Na2O–K2O–TiO2–MnO–FeO–MgO–Al2O3–SiO2–H2O indicate that mineral assemblages, zonations and microstructures observed in migmatitic rocks from the Beit Bridge Complex (Messina area, Limpopo Belt) formed along a clockwiseP–Tpath. That path displays a progradeP–Tincrease from 600 °C/7.0 kbar to 780 °C/9–10 kbar (pressure peak) and 820 °C/8 kbar (thermal peak), followed by aP–Tdecrease to 600 °C/4 kbar. The data used to construct theP–Tpath were derived from three samples of migmatitic gneiss from a restricted area, each of which has a distinct bulk composition: (1) a K, Al‐rich garnet–biotite–cordierite–sillimanite–K‐feldspar–plagioclase–quartz–graphite gneiss (2) a K‐poor, Al‐rich garnet–biotite–staurolite–cordierite–kyanite–sillimanite–plagioclase–quartz–rutile gneiss, and (3) a K, Al‐poor, Fe‐rich garnet–orthopyroxene–biotite–chlorite–plagioclase–quartz–rutile–ilmenite gneiss. Preservation of continuous prograde garnet growth zonation demonstrates that the pro‐ and retrogradeP–Tevolution of the gneisses must have been rapid, occurring during a single orogenic cycle. These petrological findings in combination with existing geochronological and structural data show that granulite facies metamorphism of the Beit Bridge metasedimentary rocks resulted from an orogenic event during the Palaeoproterozoic (c. 2.0 Ga), caused by oblique collision between the Kaapvaal and Zimbabwe Cratons. Abbreviations followKretz (1983).

Petrogenesis of andalusite–kyanite–sillimanite veins and host rocks, Sanandaj‐Sirjan metamorphic belt, Hamadan, Iran

AbstractQuartz‐rich veins in metapelitic schists of the Sanandaj‐Sirjan belt, Hamadan region, Iran, commonly contain two Al2SiO5 polymorphs, and, more rarely, three coexisting Al2SiO5 polymorphs. In most andalusite and sillimanite schists, the types of polymorphs in veins correlate with Al2SiO5 polymorph(s) in the host rocks, although vein polymorphs are texturally and compositionally distinct from those in adjacent host rocks; e.g. vein andalusite is enriched in Fe2O3 relative to host rock andalusite. Low‐grade rocks contain andalusite + quartz veins, medium‐grade rocks contain andalusite + sillimanite + quartz ± plagioclase veins, and high‐grade rocks contain sillimanite + quartz + plagioclase veins/leucosomes. Although most andalusite and sillimanite‐bearing veins occur in host rocks that also contain Al2SiO5, kyanite‐quartz veins crosscut rocks that lack Al2SiO5 (e.g. staurolite schist, granite). A quartz vein containing andalusite + kyanite + sillimanite + staurolite + muscovite occurs in andalusite–sillimanite host rocks. Textural relationships in this vein indicate the crystallization sequence andalusite to kyanite to sillimanite. This crystallization sequence conflicts with the observation that kyanite‐quartz veins post‐date andalusite–sillimanite veins and at least one intrusive phase of a granite that produced a low‐pressure–high‐temperature contact aureole; these relationships imply a sequence of andalusite to sillimanite to kyanite. Varying crystallization sequences for rocks in a largely coherent metamorphic belt can be explained by P–T paths of different rocks passing near (slightly above, slightly below) the Al2SiO5 triple point, and by overprinting of multiple metamorphic events in a terrane that evolved from a continental arc to a collisional orogen.

Oxygen isotope thermometry of quartz–Al2SiO5veins in high-grade metamorphic rocks on Naxos island (Greece)

Diffusion models predict that peak metamorphic temperatures are best recorded by the oxygen isotope fractionation between minerals in a bi-mineralic rock in which a refractory accessory mineral with slow oxygen diffusion rate is modally minor to a mineral with a faster diffusion rate. This premise is demonstrated for high-grade metamorphism on the island of Naxos, Greece, where quartz–kyanite oxygen isotope thermometry from veins in high-grade metamorphic pelites gives temperatures of 635–690 °C. These temperatures are in excellent agreement with independent thermometry for the regional M2 peak metamorphic conditions and show that the vein minerals isotopically equilibrated at the peak of metamorphism. Quartz–sillimanite fractionations in the same veins give similar temperatures (680±35 °C) and suggest that the veins grew near to the kyanite–sillimanite boundary, corresponding to pressures of 6.5 to 7.5 kbar for temperatures of 635–685 °C. By contrast, quartz–kyanite and quartz–biotite pairs in the host rocks yield lower temperature estimates than the veins (590–600 and 350–550 °C, respectively). These lower apparent temperatures are also predicted from calculations of diffusional resetting in the polyphase host-rock system. The data demonstrate that bimineralic vein assemblages can be used as accurate thermometers in high-temperature rocks whereas retrograde exchange remains a major problem in many polymineralic rocks.

U–Pb isotopic constraints on diachronous metamorphism in the northern Monashee complex, southern Canadian Cordillera

U–Pb isotopic data from the northern Monashee complex, one of the deepest structural exposures in the southern Canadian Cordillera, indicate that the age of metamorphism varies according to structural position in a 6 km thick section. This metamorphism resulted in an unusual sequence in which rocks with the lowest‐grade mineral assemblage (kyanite–sillimanite–staurolite–muscovite) are underlain and overlain by higher‐grade rocks. Xenotime and monazite U–Pb dates vary progressively from 64 Ma in the structurally highest rocks to 49 Ma in the deepest rocks. Discordant U–Pb ages from Proterozoic and Cretaceous monazite and titanite are used to interpret the thermal significance of the early Tertiary dates. The discordant analyses define linear arrays with lower intercepts that broadly overlap with early Tertiary, and the amount of discordance varies with structural level; it is least in the deeper rocks and greatest in higher rocks. Electron microprobe work showed that the monazite discordance in the deeper rocks resulted from Tertiary mineral overgrowth and recrystallization rather than Pb diffusion. We use previous studies of Pb diffusion and the fact that Proterozoic monazite and titanite suffered only negligible to moderate amounts of diffusive Pb loss to contend that elevated temperatures (c. 600–650 °C are inferred from pelitic mineral assemblages) existed in the deeper rocks for a short duration, perhaps a few million years. The downwards younging 64–49 Ma U–Pb dates are interpreted as closely reflecting xenotime and monazite growth ages rather than cooling ages or substantially reset ages based on the lack of Pb diffusion in monazite and the previously obtained 40Ar/39Ar data which suggest that rapid cooling occurred immediately after the U–Pb dates. In addition, growth ages are interpreted as thermal peak ages based on U–Pb dates from coeval kyanite‐bearing leucosomes, the consistent nature of the U–Pb dates throughout the study area, and petrographic relationships which suggest that monazite grew before or during development of the syn‐metamorphic foliation. These interpretations lead us to conclude that metamorphism was diachronous according to structural level, with higher rocks attaining peak temperatures and cooling rapidly while deeper rocks were heating towards a thermal peak that was attained a few million years later. This thermal scenario requires that higher rocks cannot have been the heat source for the deeper metamorphism, as was previously proposed.