Phase equilibria modelling of kyanite‐bearing anatectic paragneisses from the central Grenville Province

Kyanite‐bearing paragneisses from the Manicouagan Imbricate Zone and its footwall (high‐P belt of the central Grenville Province) preserve evidence of partial melting with development of metamorphic textures involving biotite–garnet ± kyanite ± plagioclase ± K‐feldspar–quartz. Garnet in these rocks displays a variety of zoning patterns with respect to Ca. Pseudosection modelling in the Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O (NCKFMASHTO) system using measured bulk rock compositions accounts for the textural evolution of two aluminous and two sub‐aluminous samples from the presumed thermal peak to conditions at which retained melt solidified. The prograde features are best explained by pseudosections calculated with compositions to account for melt loss. The intersection of isopleths of grossular content and Fe/(Fe + Mg) relating to large porphyroblasts of garnet provide constraints on the P–T conditions of the metamorphic peak. These P–T estimates are considered to be minima because of the potential for diffusional modification of the composition of garnet at high‐T and during the early stages of cooling. However, they are consistent with textural observations and pseudosection topology, with peak assemblages best preserved in rocks for which the calculated pseudosections predict only small changes in mineral proportions in the P–T interval, in which retrograde reactions are inferred to have occurred between the thermal peak and the solidus. Maximum P–T conditions (14.5–15.5 kbar and 840–890 °C) and steep retrograde P–T paths inferred for rocks from the Manicouagan Imbricate Zone are comparable with those determined for mafic rocks from the same area. In contrast, maximum P–T conditions of 12.5–13 kbar and 815–830 °C and flatter P–T paths are inferred for the rocks of the footwall to the Manicouagan Imbricate Zone. The general consistency between textures, mineral compositions and the topologies of the calculated pseudosections suggests that the pseudosection approach is an appropriate tool for inferring the P–T evolution of high‐P anatectic quartzo‐feldspathic rocks.     

An automated method for the calculation of P–T paths from garnet zoning,with application to metapelitic schist from the Kootenay Arc,British Columbia,Canada

An automated method for the calculation of P–T paths based on garnet zoning is presented and used to interpret zoning in metapelitic schist from the southern Canadian Cordillera. The approach adopted to reconstruct the P–T path is to match garnet compositions along a radial transect with predictions from thermodynamic forward models, while iteratively modifying the composition to account for fractional crystallization. The method is applied to a representative sample of garnet‐ and staurolite‐bearing schist from an amphibolite facies Barrovian belt in the southern Canadian Omineca belt. Garnet zoning in these schists is concentric and largely continuous from core to rim. Three zones are present, the first two of which coincide with sector‐zoned cores of garnet crystals. Similar zoning is developed in rocks that contain or lack staurolite, respectively, suggesting garnet growth was restricted to the initial part of the prograde P–T path prior to the development of staurolite. Growth zoning in large garnet crystals has not been significantly modified by diffusion. This interpretation is based on zoning characteristics of garnet crystals and is further supported by results of a forward model incorporating the effects of simultaneous fractional crystallization and intracrystalline diffusion. The P–T path calculated for this rock includes an initial, linear stage with a high dP/dT, and a later stage dominated by heating. The approach adopted in this study may have application to other garnet‐bearing rocks in which growth zoning is preserved.     

Constraining the conditions of Barrovian metamorphism in Sikkim,India: P–T–t paths of garnet crystallization in the Lesser Himalayan Belt

Garnet crystallization in metapelites from the Barrovian garnet and staurolite zones of the Lesser Himalayan Belt in Sikkim is modelled utilizing Gibbs free energy minimization, multi‐component diffusion theory and a simple nucleation and growth algorithm. The predicted mineral assemblages and garnet‐growth zoning match observations remarkably well for relatively tight, clockwise metamorphic P–T paths that are characterized by prograde gradients of ～30 °C kbar^?1 for garnet‐zone rocks and ～20 °C kbar^?1 for rocks from the staurolite zone. Estimates for peak metamorphic temperature increase up‐structure toward the Main Central Thrust. According to our calculations, garnet stopped growing at peak pressures, and protracted heating after peak pressure was absent or insignificant. Almost identical P–T paths for the samples studied and the metamorphic continuity of the Lesser Himalayan Belt support thermo‐mechanical models that favour tectonic inversion of a coherent package of Barrovian metamorphic rocks. Time‐scales associated with the metamorphism were too short for chemical diffusion to substantially modify garnet‐growth zoning in rocks from the garnet and staurolite zones. In general, the pressure of initial garnet growth decreases, and the temperature required for initial garnet growth was reached earlier, for rocks buried closer toward the MCT. Deviations from this overall trend can be explained by variations in bulk‐rock chemistry.     

Age and P–T conditions of the Gridino‐type eclogite in the Belomorian Province,Russia

Although eclogites in the Belomorian Province have been regarded as Archean in age and among the oldest in the world, there are also multiple studies that have proposed a Paleoproterozoic age. Here, we present new data for the Gridino‐type eclogites, which occur as boudins and metamorphosed dykes within tonalite–trondhjemite–granodiorite gneisses. Zircon from these eclogites has core and rim structures. The cores display high Th/U ratios (0.18–0.45), negative Eu anomalies and strong enrichment in HREE, and have Neoarchean U–Pb ages of c. 2.70 Ga; they are interpreted to be magmatic in origin. Zircon cores have δ¹⁸O of 5.64–6.07‰ suggesting the possibility of crystallization from evolved mantle‐derived magmas. In contrast, the rims, which include the eclogite facies minerals omphacite and garnet, are characterized by low Th/U ratios (<0.035) and flat HREE patterns, and yield U–Pb ages of c. 1.90 Ga; they are interpreted to be metamorphic in origin. Zircon rims have elevated δ¹⁸O of 6.23–6.80‰, which was acquired during eclogite facies metamorphism. Based on petrography and phase equilibria modelling, we recognize a prograde epidote amphibolite facies mineral assemblage, the peak eclogite facies mineral assemblage and a retrograde high‐P amphibolite facies mineral assemblage. The peak metamorphic conditions of 695–755°C at >18 kbar for the Gridino‐type eclogites suggest an apparent thermal gradient of <39–42°C/kbar for the Lapland–Kola collisional orogeny.     

Influence of tectonic overpressure on P–T paths of HP–UHP rocks in continental collision zones: thermomechanical modelling

The principle of lithostatic pressure is habitually used in metamorphic geology to calculate burial/exhumation depth from pressure given by geobarometry. However, pressure deviation from lithostatic, i.e. tectonic overpressure/underpressure due to deviatoric stress and deformation, is an intrinsic property of flow and fracture in all materials, including rocks under geological conditions. In order to investigate the influences of tectonic overpressure on metamorphic P–T paths, 2D numerical simulations of continental subduction/collision zones were conducted with variable brittle and ductile rheologies of the crust and mantle. The experiments suggest that several regions of significant tectonic overpressure and underpressure may develop inside the slab, in the subduction channel and within the overriding plate during continental collision. The main overpressure region that may influence the P–T paths of HP–UHP rocks is located in the bottom corner of the wedge‐like confined channel with the characteristic magnitude of pressure deviation on the order of ～0.3 GPa and 10–20% from the lithostatic values. The degree of confinement of the subduction channel is the key factor controlling this magnitude. Our models also suggest that subducted crustal rocks, which may not necessarily be exhumed, can be classified into three different groups: (i) UHP‐rocks subjected to significant (≥0.3 GPa) overpressure at intermediate subduction depth (50–70 km, P = 1.5–2.5 GPa) then underpressured at depth ≥100 km (P ≥ 3 GPa); (ii) HP‐rocks subjected to ≥0.3 GPa overpressure at peak P–T conditions reached at 50–70 km depth in the bottom corner of the wedge‐like confined subduction channel (P = 1.5–2.5 GPa); (iii) lower‐pressure rocks formed at shallower depths (≤40 km depth, P ≤ 1 GPa), which are not subjected to significant overpressure and/or underpressure.     

The significance of Mn‐rich ilmenite and the determination of P–T paths from zoned garnet in metasedimentary rocks from the western Cape Breton Highlands,Nova Scotia

The Jumping Brook Metamorphic Suite in the western Cape Breton Highlands of Nova Scotia is part of an inverted Barrovian sequence that formed during a Late Silurian–Early Devonian promontory–promontory collision in the Canadian Appalachians. In this study, systematic discrepancies between geochemical observations and thermodynamic model predictions led to the discovery of a systematic relationship linking the style of garnet core isopleth intersection (GCII) to the pyrophanite (MnTiO₃) component of co‐existing ilmenite. Samples that yielded tight GCIIs at or near the garnet‐in curve were found to contain ilmenite with negligible pyrophanite components, whereas samples yielding GCIIs far removed (up to 105°C) from the garnet‐in curve were found to contain ilmenite with significant pyrophanite and/or ecandrewsite (ZnTiO₃) components. Based on petrographic and geochemical observations, Mn(±Zn)‐rich ilmenite are interpreted to have sequestered Mn throughout prograde metamorphism due to sluggish intracrystalline diffusion. The amount of reactive Mn input into the thermodynamic models from whole‐rock analyses were, in some cases, overestimated, resulting in garnet‐in curve topologies that extend to erroneously low P–T conditions. Modifications to the whole‐rock chemistry that account for Mn sequestration into ilmenite, however, yielded robust model results. Our results show that, in addition to uncertainties in thermodynamic data sets and phenomenon related to reaction kinetics, Mn‐rich ilmenite may superimpose additional complexities related to the interpretation of predicted equilibria involving garnet. Numerical simulations of garnet crystallization were used to infer P–T paths of metamorphism for one sample from the garnet zone (Mn corrected) and two samples from the staurolite zone (Mn uncorrected) of the inverted sequence. Model results are remarkably similar among the three samples and indicate that garnet crystallization occurred along relatively steep (31–37°C/km) clockwise P–T paths. The peak conditions of garnet crystallization and metamorphism (560–590°C, 7.4–8.0 kbar) are interpreted to have been attained approximately simultaneously, such that the paths are characterized by tight prograde‐to‐retrograde transitions. The hairpin nature of the P–T paths is interpreted to represent the onset of thrust‐related exhumation and isograd inversion along ductile shear zones, consistent with available field and geochronological constraints.     

Metamorphic P–T profile and P–T path discontinuity across the far‐eastern Nepal Himalaya: investigation of channel flow models

Thermobarometric data and compositional zoning of garnet show the discontinuities of both metamorphic pressure conditions at peak‐T and P–T paths across the Main Central Thrust (MCT), which juxtaposes the high‐grade Higher Himalayan Crystalline Sequences (HHCS) over the low‐grade Lesser Himalaya Sequences (LHS) in far‐eastern Nepal. Maximum recorded pressure conditions occur just above the MCT (～11 kbar), and decrease southward to ～6 kbar in the garnet zone and northward to ～7 kbar in the kyanite ± staurolite zone. The inferred nearly isothermal loading path for the LHS in the staurolite zone may have resulted from the underthrusting of the LHS beneath the HHCS. In contrast, the increasing temperature path during both loading and decompression (i.e. clockwise path) from the lowermost HHCS in the staurolite to kyanite ± staurolite transitional zone indicates that the rocks were fairly rapidly buried and exhumed. Exhumation of the lowermost HHCS from deeper crustal depths than the flanking regions, recording a high field pressure gradient (～1.2–1.6 kbar km^?1) near the MCT, is perhaps caused by ductile extrusion along the MCT, not the emplacement along a single thrust, resulting in the P–T path discontinuities. These observations are consistent with the overall scheme of the model of channel flow, in which the outward flowing ‘HHCS’ and inward flowing ‘LHS’ are juxtaposed against each other and are rapidly extruded together along the ‘MCT’. A rapid exhumation by channel flow in this area is also suggested by a nearly isothermal decompression path inferred from cordierite corona surrounding garnet in gneiss of the upper HHCS. However, peak metamorphic temperatures show a progressive increase of temperature structurally upward (～570–740 °C) near the MCT and roughly isothermal conditions (～710–810 °C) in the upper structural levels of the HHCS. The observed field temperature gradient is much lower than those predicted in channel flow models. However, the discrepancy could be resolved by taking into account heat advection by melt and/or fluid migration, as these can produce low or nearly no field temperature gradient in the exhumed midcrust, as observed in nature.     

A curious case of agreement between conventional thermobarometry and phase equilibria modelling in granulites: New constraints on P–T estimates in the Antarctica segment of the Musgrave–Albany–Fraser–Wilkes Orogen

The Windmill Islands region in Wilkes Land, east Antarctica, preserves granulite facies metamorphic mineral assemblages that yield seemingly comparable P–T estimates from conventional thermobarometry and mineral equilibria modelling. This is uncommon in granulite facies terranes, where conventional thermobarometry and phase equilibria modelling generally produce conflicting P–T estimates because peak mineral compositions tend to be modified by retrograde diffusion processes. In situ U–Pb monazite geochronology and calculated metamorphic phase diagrams show that the Windmill Islands experienced two phases of high thermal gradient metamorphism during the Mesoproterozoic. The first phase of metamorphism is recorded by monazite ages in two widely separated samples and occurred at c. 1,305 Ma. This event was regional in extent, involved crustally derived magmatism and reached conditions of ~3.2–5 kbar and 690–770°C corresponding to very high thermal gradients of >150°C/kbar. The elevated thermal regime is interpreted to reflect a period of extension or increased extension in a back‐arc setting that existed prior to c. 1,330 Ma. The first metamorphic event was overprinted by granulite facies metamorphism at c. 1,180 Ma that was coeval with the intrusion of charnockite. This event involved peak temperatures of ~840–850°C and pressures of ~4–5 kbar. A phase of granitic magmatism at c. 1,250–1,210 Ma, prior to the intrusion of the charnockite, is interpreted to reflect a phase of compression within an overall back‐arc setting. Existing conventional thermobarometry suggests conditions of ~4 kbar and 750°C for M₁ and 4–7 kbar and 750–900°C for M₂. The apparent similarities between the phase equilibria modelling and existing conventional thermobarometry may suggest either that the terrane cooled relatively quickly, or that the P–T ranges obtained from conventional thermobarometry are sufficiently imprecise that they cover the range of P–T conditions obtained in this study. However, without phase equilibria modelling, the veracity of existing conventional P–T estimates cannot be evaluated. The calculated phase diagrams from this study allow the direct comparison of P–T conditions in the Windmill Islands with phase equilibria models from other regions in the Musgrave–Albany–Fraser–Wilkes Orogen. This shows that the metamorphic evolution of the Wilkes Land region is very similar to that of the eastern Albany–Fraser Orogen and Musgrave Province in Australia, and further demonstrates the remarkable consistency in the timing of metamorphism and the thermal gradients along the ~5,000 km strike length of this system.     

P–T evolution of Precambrian eclogite in the Sveconorwegian orogen,SW Sweden

Conditions of the prograde, peak‐pressure and part of the decompressional P–T path of two Precambrian eclogites in the eastern Sveconorwegian orogen have been determined using the pseudosection approach. Cores of garnet from a Fe–Ti‐rich eclogite record a first prograde and syn‐deformational stage along a Barrovian gradient from ~670 °C and 7 kbar to 710 °C and 8.5 kbar. Garnet rims grew during further burial to 16.5–19 kbar at ~850–900 °C, along a steep dP/dT gradient. The pseudosection model of a kyanite‐bearing eclogite sample of more magnesian bulk composition confirms the peak conditions. Matrix reequilibration associated with subsequent near‐isothermal decompression and partial exhumation produced plagioclase‐bearing symplectites replacing kyanite and clinopyroxene at an estimated 850–870 °C and 10–11 kbar. The validity of the pseudosections is discussed in detail. It is shown that in pseudosection modelling the fractionation of FeO in accessory sulphides may cause a significant shift of field boundaries (here displaced by up to 1.5 kbar and 70 °C) and must not be neglected. Fast burial, exhumation and subsequent cooling are supported by the steepness of both the prograde and the decompressional P–T paths as well as the preservation of garnet growth zoning and the symplectitic reaction textures. These features are compatible with deep tectonic burial of the eclogite‐bearing continental crust as part of the underthrusting plate (Eastern Segment, continent Baltica) in a collisional setting that led to an effectively doubled crustal thickness and subsequent exhumation of the eclogites through tectonic extrusion. Our results are in accordance with regional structural and petrological relationships, which demonstrate foreland‐vergent partial exhumation of the eclogite‐bearing nappe along a basal thrust zone and support a major collisional stage at c. 1 Ga. We argue that the similarities between Sveconorwegian and Himalayan eclogite occurrences emphasize the modern style of Grenvillian‐aged tectonics.     

Tectonometamorphic evolution of an intracontinental orogeny inferred from P–T–t–d paths of the metapelites from the Rehamna massif (Morocco)

New petrographic and microstructural observations, mineral equilibria modelling and U/Pb (monazite) geochronological studies were carried out to investigate the relationships between deformation and metamorphism across the Rehamna massif (Moroccan Variscan belt). In this area, typical Barrovian (muscovite to staurolite) zones developed in Cambrian to Carboniferous metasedimentary rocks that are distributed around a dome‐like structure. First assemblages are characterized by the presence of locally preserved andalusite, followed by prograde evolution culminating at 6 kbar and 620 °C in the structurally deepest staurolite zone rocks. This Barrovian sequence was subsequently uplifted to supracrustal levels, heterogeneously reworked at greenschist facies conditions, which was followed locally by static growth of andalusite, indicating heating to 2.5–4 kbar and 530–570 °C. The ²⁰⁶Pb/²³⁸U monazite age of 298.3 ± 4.1 Ma is interpreted as minimum age of peak metamorphic conditions, whereas the ages of 275.8 ± 1.7 Ma and 277.0 ± 1.1 Ma date decompression and heating at low pressure, in agreement with previous dating of Permian granitoids intruding the Rehamna massif. The prograde metamorphism occurred during thickening and associated horizontal flow in the deeper crust (S1 horizontal schistosity). The horizontally disposed metamorphic zones were subsequently uplifted by a regional scale antiform during ongoing N–S compression. The re‐heating of the massif follows late massive E–W shortening, refolding and retrograde shearing of all previous fabrics coevally with regionally important intrusions of Permian granitoids. We argue that metamorphic evolution of the Rehamna massif occurred several hundred kilometres from the convergent plate boundaries in the interior of continental Gondwanan plate. The tectonometamorphic history of the Rehamna massif is put into Palaeozoic plate tectonic perspective and Late Carboniferous reactivation of (Devonian)–Early Carboniferous basins formed during stretching of the north Gondwana margin and formation of the Palaeotethys Ocean. The inherited heat budget of these magma‐rich basins plays a role in the preferential location of this intracontinental orogen. It is shown that rapid transition from lithospheric stretching to compression is characterized by specific HT type of Barrovian metamorphism, which markedly differs from similar Barrovian sequences along Palaeozoic plate boundaries reported from Variscan Europe.     

Osumilite–melt interactions in ultrahigh temperature granulites: phase equilibria modelling and implications for the P–T–t evolution of the Eastern Ghats Province,India

The exposed residual crust in the Eastern Ghats Province records ultrahigh temperature (UHT) metamorphic conditions involving extensive crustal anatexis and melt loss. However, there is disagreement about the tectonic evolution of this late Mesoproterozoic–early Neoproterozoic orogen due to conflicting petrological, structural and geochronological interpretations. One of the petrological disputes in residual high Mg–Al granulites concerns the origin of fine‐grained mineral intergrowths comprising cordierite + K‐feldspar ± quartz ± biotite ± sillimanite ± plagioclase. These intergrowths wrap around porphyroblast phases and are interpreted to have formed by the breakdown of primary osumilite in the presence of melt trapped in the equilibration volume by the melt percolation threshold. The pressure (P)–temperature (T) evolution of four samples from three localities across the central Eastern Ghats Province is constrained using phase equilibria modelling in the chemical system Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–Fe₂O₃ (NCKFMASHTO). Results of the modelling are integrated with published geochronological results for these samples to show that the central Eastern Ghats Province followed a common P–T–t history. This history is characterized by peak UHT metamorphic conditions of 945–955 °C and 7.8–8.2 kbar followed by a slight increase in pressure and close‐to‐isobaric cooling to the conditions of the elevated solidus at 940–900 °C and 8.5–8.3 kbar. In common with other localities from the Eastern Ghats Province, the early development of cordierite before osumilite and the peak to immediate post‐peak retrograde reaction between osumilite and melt to produce the intergrowth features requires that the prograde evolution was one of contemporaneous increasing pressure with increasing temperature. This counter‐clockwise (CCW) evolution is evaluated for one sample using inverse phase equilibria modelling along a schematic P–T path of 150 °C kbar^?1 starting from the low P–T end of the prograde P–T path as constrained by the phase equilibria modelling. The inverse modelling is executed by step‐wise down temperature reintegration of sufficient melt into the residual bulk chemical composition at the P–T point of the 1 mol.% melt isopleth at each step, representing the melt remaining on grain boundaries after each prograde drainage event, to reach the melt connectivity transition (MCT) of 7 mol.%. The procedure is repeated until a plausible protolith composition is recovered. The result demonstrates that clastic sedimentary rocks that followed a CCW P–T evolution could have produced the observed mineral assemblages and microstructures preserved in the central Eastern Ghats Province. This study also highlights the role of melt during UHT metamorphism, particularly its importance to both chemical and physical processes along the prograde and retrograde segments of the P–T path. These processes include: (i) an increase in diffusive length scales during the late prograde to peak evolution, creating equilibration volumes larger than a standard thin section; (ii) the development of retrograde mineral assemblages, which is facilitated if some melt is retained post‐peak; (iii) the presence of melt as a weakening mechanism and the advection of heat by melt, allowing the crust to thicken; and (iv) the effect of melt loss, which makes the deep crust both denser and stronger, and reduces heat production at depth, limiting crustal thickening and facilitating the transition to close‐to‐isobaric cooling.     

Crustal architecture above the high-pressure belt of the Grenville Province in the Manicouagan area: new structural, petrologic and U–Pb age constraints

Recent U–Pb age determinations and P–T estimates allow us to characterize the different levels of a formerly thickened crust, and provide further constraints on the make up and tectono-thermal evolution of the Grenville Province in the Manicouagan area. An important tectonic element, the Manicouagan Imbricate zone (MIZ), consists of mainly 1.65, 1.48 and 1.17 Ga igneous rocks metamorphosed under 1400–1800 MPa and 800–900 °C at 1.05–1.03 Ga, during the Ottawan episode of the Grenvillian orogenic cycle, coevally with intrusion of gabbro dykes in shear zones. The MIZ has been interpreted as representing thermally weakened deep levels of thickened crust extruded towards the NW over a parautochthonous crustal-scale ramp. Mantle-derived melts are considered as in part responsible for the high metamorphic temperatures that were registered.New data show that mid-crustal levels structurally above the MIZ are represented by the Gabriel Complex of the Berthé terrane, that consists of migmatite with boudins of 1136±15 Ma gabbro and rafts of anatectic metapelite with an inherited monazite age at 1478±30 Ma. These rocks were metamorphosed at about the same time as the MIZ (metamorphic zircon in gabbro: 1046±2 Ma; single grains of monazite in anatectic metapelite: 1053±2 Ma) and under the same T range (800–900 °C) but at lower P conditions (1000–1100 MPa). They are mainly exposed in an antiformal culmination above a high-strain zone, which has tectonic lenses of high P–T rocks from the MIZ and is intruded by synmetamorphic gabbroic rocks. This zone is interpreted as part of the hangingwall of the MIZ during extrusion. A gap of 400 MPa in metamorphic pressures between the tectonic lenses and the country rocks, together with the broad similarity in metamorphic ages, are consistent with rapid tectonic transport of the high P–T rocks over a ramp prior to the incorporation of the mafic lenses in the hangingwall.Between the antiformal culmination of the Gabriel Complex and the MIZ 1.48 Ga old granulites of the Hart Jaune terrane are exposed. They are intruded by unmetamorphosed 1228±3 Ma gabbro sills and 1166±1 Ma anorthosite. Hart Jaune Terrane represents relatively high crustal levels that truncate the MIZ-Gabriel Complex contact and are preserved in a synformal structure.Farther south, the Gabriel Complex is overlain by the Banded Complex, a composite unit including 1403+32/−25 Ma granodiorite and 1238+16/−13−1202+40/−25 Ma granite. This unit has been metamorphosed under relatively low-P (800 MPa) granulite-facies conditions. Metamorphic U–Pb data, limited to zircon lower intercept ages (971±38 Ma and 996±27 Ma) and a titanite (990±5 Ma) age, are interpreted to postdate the metamorphic peak.The general configuration of units along the section is consistent with extrusion of the MIZ during shortening and, finally, normal displacement along discrete shear zones.     

Subduction initiation in the Neo‐Tethys: constraints from counterclockwise P–T paths in amphibolite rocks of the Nagaland Ophiolite Complex,India

High‐P metamorphic rocks that are formed at the onset of oceanic subduction usually record a single cycle of subduction and exhumation along counterclockwise (CCW) P–T paths. Conceptual and thermo‐mechanical models, however, predict multiple burial–exhumation cycles, but direct observations of these from natural rocks are rare. In this study, we provide a new insight into this complexity of subduction channel dynamics from a fragment of Middle‐Late Jurassic Neo‐Tethys in the Nagaland Ophiolite Complex, northeastern India. Based on integrated textural, mineral compositional, metamorphic reaction history and geothermobarometric studies of a medium‐grade amphibolite tectonic unit within a serpentinite mélange, we establish two overprinting metamorphic cycles (M₁–M₂). These cycles with CCW P–T trajectories are part of a single tectonothermal event. We relate the M₁ metamorphic sequence to prograde burial and heating through greenschist and epidote blueschist facies to peak metamorphism, transitional between amphibolite and hornblende‐eclogite facies at 13.8 ± 2.6 kbar, 625 ± 45 °C (error 2σ values) and subsequent cooling and partial exhumation to greenschist facies. The M₂ metamorphic cycle reflects epidote blueschist facies prograde re‐burial of the partially exhumed M₁ cycle rocks to peak metamorphism at 14.4 ± 2 kbar, 540 ± 35 °C and their final exhumation to greenschist facies along a relatively cooler exhumation path. We interpret the M₁ metamorphism as the first evidence for initiation of subduction of the Neo‐Tethys from the eastern segment of the Indus‐Tsangpo suture zone. Reburial and final exhumation during M₂ are explained in terms of material transport in a large‐scale convective circulation system in the subduction channel as the latter evolves from a warm nascent to a cold and more mature stage of subduction. This Neo‐Tethys example suggests that multiple burial and exhumation cycles involving the first subducted oceanic crust may be more common than presently known.     

Linking thermodynamic modelling,Lu–Hf geochronology and trace elements in garnet: new P–T–t paths from the Sevier hinterland

Major element, trace element and Lu–Hf geochronological data from amphibolite facies pelitic schist in the Raft River and Albion Mountains of northwest Utah and southern Idaho indicate that garnet grew during increasing pressure, interpreted to be the result of tectonic burial and crustal thickening during Sevier orogenesis. Garnet growth was interrupted by hiatuses interpreted from discontinuities in major element zonation. Pressure–temperature paths were determined from the pre‐hiatus portions of the garnet chemical zoning profiles and indicate an increase of ~2 kbar and ~50 °C in the western Raft River Mountains. Garnet Lu–Hf dates of 150 ± 1 Ma in the western Raft River Mountains and 138.7 ± 0.7 Ma and 132 ± 5 Ma in the southern Albion Mountains indicate the timing of garnet growth. Lutetium garnet zoning profiles indicate that the Lu–Hf ages are biased towards the post‐hiatus or outer pre‐hiatus segments, indicating that the determined ages likely post‐date the recorded P–T path history or date the tail end of the paths. Crustal thickening associated with Sevier orogenesis in the western Raft River Mountains thus began slightly before 150 ± 1 Ma, in the Late Jurassic. This study shows that integrating P–T paths determined from garnet growth zoning with Lu–Hf garnet geochronology and in situ garnet trace element analyses is an effective approach for interpreting and dating deformation events in orogenic belts.     

Metamorphic textures and P–T evolution of high-P granulites from the Lelukuau terrane, NE Grenville Province

High‐pressure (HP) granulites and eclogitized metagabbro are exposed along an orogen‐parallel high‐P belt that was developed at c. 1050–1020 Ma in the NE Grenville Province. Among these rocks, mafic granulites derived from a Labradorian anorthosite suite of the Lelukuau terrane contain garnet, Al‐Na diopside, and, depending on bulk composition, plagioclase and kyanite. Moreover, the distribution of phases is influenced by the original igneous texture. For instance, in high X_MgO leucocratic varieties, garnet porphyroblasts nucleated together with kyanite in An‐rich cores of plagioclase domains whereas in low X_MgO rocks garnet occurs together with clinopyroxene within formerly igneous ferromagnesian domains and kyanite is missing. In contrast, garnet pseudomorphs after igneous plagioclase in melanocratic varieties display evidence of earlier corona development. Metamorphic textures are consistent with a two stage evolution: (a) development of garnet and Al‐Na‐diopside (Cpx₁) under high‐P metamorphic conditions, concomitant with elimination of plagioclase in the mesocratic to melanocratic varieties; and (b) partial loss of Al‐Na from Cpx₁ resulting in production of new andesitic plagioclase, and growth of new clinopyroxene (Cpx₂) after garnet and quartz in leucocratic to mesocratic rocks consistent with decompression. Widespread equilibrium textures between garnet‐Pl₂‐Cpx₂ and/or reset Cpx₁ are consistent with development at the thermal peak. Estimated P–T conditions for the presumed thermal peak fall in the range 1500–1800 MPa and 800–900 °C and are comparable to those recorded by eclogitized gabbros from other parts of the high‐P belt of the NE Grenville province. Low jadeite content of clinopyroxene from the HP granulites is attributed to the low bulk Na₂O/(Na₂O + CaO) of these rocks relative to common basaltic compositions. Scarcity of apparent retrograde textural overprint in both the HP granulites and the eclogites suggests fast subsequent cooling, consistent with extrusion of the high‐P belt towards the foreland shortly after the metamorphic peak.     

Anticlockwise pressure–temperature paths record Variscan upper‐plate exhumation: Example from micaschists of the Porto Vecchio region,Corsica

To better understand the evolution of deep‐seated crust of the Variscan orogen in the Sardinia‐Corsica region, we studied garnet‐bearing micaschists which were sampled 3 km east and 15 km northeast of Porto Vecchio, south‐eastern Corsica. After a careful investigation of the textural relations and compositions of minerals, especially of zoned garnet, a P–T path was reconstructed using contoured P–T pseudosections. U–Th–Pb dating of monazite in the micaschists was undertaken with the electron microprobe. The micaschists from both localities were formed along similar anticlockwise P–T paths. The prograde branch of these paths starts at 3 kbar close to 600°C in the P–T field of sillimanite and reaches peak conditions at 7 kbar and 600 (15 km NE of Porto Vecchio) to 630°C (3 km E of Porto Vecchio). The metamorphism at peak P–T conditions happened c. 340 Ma based on low‐Y (<0.65 wt% Y₂O₃) monazite. Ages of monazite with high‐Y contents (>2 wt% Y₂O₃), which probably have formed before garnet, scatter around 362 Ma. The retrograde branch of the P–T paths passes through 4 kbar at ~550°C. We conclude that the micaschists belong to a common metasedimentary sequence, which extends over the Porto Vecchio region and is separated from other metamorphic rock sequences in the north and the south by major tectonic boundaries. This sequence had experienced peak pressures which are lower than those determined for metamorphic rocks, such as micaschist and gneiss, from north‐eastern Sardinia. At present, we favour a continent–continent collisional scenario with the studied metasedimentary sequence buried during the collisional event as part of the upper plate. The contemporaneous high‐P metamorphic rocks from NE Sardinia were part of the upper portion of the lower plate. The addressed rocks from both plates were exhumed in an exhumation channel.     

P–T–t evolution of granulite facies metamorphism and partial melting in the Winding Stair Gap,Central Blue Ridge,North Carolina,USA

The Winding Stair Gap in the Central Blue Ridge province exposes granulite facies schists, gneisses, granofelses and migmatites characterized by the mineral assemblages: garnet–biotite–sillimanite–plagioclase–quartz, garnet–hornblende–biotite–plagioclase–quartz ± orthopyroxene ± clinopyroxene and orthopyroxene–biotite–quartz. Multiple textural populations of biotite, kyanite and sillimanite in pelitic schists support a polymetamorphic history characterized by an early clockwise P–T path in which dehydration melting of muscovite took place in the stability field of kyanite. Continued heating led to dehydration melting of biotite until peak conditions of 850 ± 30 °C, 9 ± 1 kbar were reached. After equilibrating at peak temperatures, the rocks underwent a stage of near isobaric cooling during which hydrous melt ± K‐feldspar were replaced by muscovite, and garnet by sillimanite + biotite + plagioclase. Most monazite crystals from a pelitic schist display patchy zoning for Th, Y and U, with some matrix crystals having as many as five compositional zones. A few monazite inclusions in garnet, as well as Y‐rich cores of some monazite matrix crystals, yield the oldest dates of c. 500 Ma, whereas a few homogeneous matrix monazites that grew in the main foliation plane yield dates of 370–330 Ma. Culling and analysis of individual spot dates for eight monazite grains yields three age populations of 509 ± 14 Ma, 438 ± 5 Ma and 360 ± 5 Ma. These data suggest that peak‐temperature metamorphism and partial melting in the central Blue Ridge occurred during the Salinic or Taconic orogeny. Following near isobaric cooling, a second weaker thermal pulse possibly related to intrusion of nearby igneous bodies resulted in growth of monazite c. 360 Ma, coinciding with the Neoacadian orogeny.     

Constraints from monazite and xenotime growth modelling in the MnCKFMASH‐PYCe system on the P–T path of a metapelite from Shillong‐Meghalaya Plateau: implications for the Indian shield assembly

The Palaeo‐Mesoproterozoic metapelite granulites from northern Garo Hills, western Shillong‐Meghalaya Gneissic Complex (SMGC), northeast India, consist of resorbed garnet, cordierite and K‐feldspar porphyroblasts in a matrix comprising shape‐preferred aggregates of biotite±sillimanite+quartz that define the penetrative gneissic fabric. An earlier assemblage including biotite and sillimanite occurs as inclusions within the garnet and cordierite porphyroblasts. Staurolite within cordierite in samples without matrix sillimanite is interpreted to have formed by a reaction between the sillimanite inclusion and the host cordierite during retrogression. Accessory monazite occurs as inclusions within garnet as well as in the matrix, whereas accessory xenotime occurs only in the matrix. The monazite inclusions in garnet contain higher Ca, and lower Y and Th/U than the matrix monazite outside resorbed garnet rims. On the other hand, matrix monazite away from garnet contains low Ca and Y, and shows very high Th/U ratios. The low Th/U ratios (<10) of the Y‐poor garnet‐hosted monazite indicate subsolidus formation during an early stage of prograde metamorphism. A calculated P–T pseudosection in the MnCKFMASH‐PYCe system indicates that the garnet‐hosted monazite formed at <3 kbar/600 °C (Stage A). These P–T estimates extend backward the previously inferred prograde P–T path from peak anatectic conditions of 7–8 kbar/850 °C based on major mineral equilibria. Furthermore, the calculated P–T pseudosections indicate that cordierite–staurolite equilibrated at ~5.5 kbar/630 °C during retrograde metamorphism. Thus, the P–T path was counterclockwise. The Y‐rich matrix monazite outside garnet rims formed between ~3.2 kbar/650 °C and ~5 kbar/775 °C (Stage B) during prograde metamorphism. If the effect of bulk composition change due to open system behaviour during anatexis is considered, the P–T conditions may be lower for Stage A (<2 kbar/525 °C) and Stage B (~3 kbar/600 °C to ~3.5 kbar/660 °C). Prograde garnet growth occurred over the entire temperature range (550–850 °C), and Stage‐B monazite was perhaps initially entrapped in garnet. During post‐peak cooling, the Stage‐B monazite grains were released in the matrix by garnet dissolution. Furthermore, new matrix monazite (low Y and very high Th/U ≤80, ~8 kbar/850–800 °C, Stage C), some monazite outside garnet rims (high Y and intermediate Th/U ≤30, ~8 kbar/800–785 °C, Stage D), and matrix xenotime (<785 °C) formed through post‐peak crystallization of melt. Regardless of textural setting, all monazite populations show identical chemical ages (1630–1578 Ma, ±43 Ma). The lithological association (metapelite and mafic granulites), and metamorphic age and P–T path of the northern Garo Hills metapelites and those from the southern domain of the Central Indian Tectonic Zone (CITZ) are similar. The SMGC was initially aligned with the southern parts of CITZ and Chotanagpur Gneissic Complex of central/eastern India in an ENE direction, but was displaced ~350 km northward by sinistral movement along the north‐trending Eastern Indian Tectonic Zone in Neoproterozoic. The southern CITZ metapelites supposedly originated in a back‐arc associated with subducting oceanic lithosphere below the Southern Indian Block at c. 1.6 Ga during the initial stage of Indian shield assembly. It is inferred that the SMGC metapelites may also have originated contemporaneously with the southern CITZ metapelites in a similar back‐arc setting.     

P–T evolution of the Precambrian Metamorphic Complex,NW Iran: a study of metapelitic rocks

The Mahneshan Metamorphic Complex (MMC) is one of the Precambrian terrains exposed in the northwest of Iran. The MMC underwent two main phases of deformation (D₁ and D₂) and at least two metamorphic events (M₁ and M₂). Critical metamorphic mineral assemblages in the metapelitic rocks testify to regional metamorphism under amphibolite‐facies conditions. The dominant metamorphic mineral assemblage in metapelitic rocks (M₁) is muscovite, biotite I, Garnet I, staurolite, Andalusite I and sillimanite. Peak metamorphism took place at 600–620°C and ∼7 kbar, corresponding to a depth of ca. 24 km. This was followed by decompression during exhumation of the crustal rocks up to the surface. The decrease of temperature and pressure during exhumation produced retrograde metamorphic assemblages (M₂). Secondary phases such as garnet II biotite II, Andalusite II constrain the temperature and pressure of M₂ retrograde metamorphism to 520–560°C and 2.5–3.5 kbar, respectively. The geothermal gradient obtained for the peak of metamorphism is 33°C km⁻¹, which indicates that peak metamorphism was of Barrovian type and occurred under medium‐pressure conditions. The MMC followed a ‘clockwise’ P–T path during metamorphism, consistent with thermal relaxation following tectonic thickening. The bulk chemistry of the MMC metapelites shows that their protoliths were deposited at an active continental margin. Together with the presence of palaeo‐suture zones and ophiolitic rocks around the high‐grade metamorphic rocks of the MMC, these features suggest that the Iranian Precambrian basement formed by an island‐arc type cratonization. Copyright © 2010 John Wiley & Sons, Ltd.     

Thermomechanical evolution of the high‐grade core in the nappe zone of Variscan Sardinia,Italy: the role of shear deformation and granite emplacement

In the nappe zone of the Sardinian Variscan chain, the deformation and metamorphic grade increase throughout the tectonic nappe stack from lower greenschist to upper amphibolite facies conditions in the deepest nappe, the Monte Grighini Unit. A synthesis of petrological, structural and radiometric data is presented that allows us to constrain the thermal and mechanical evolution of this unit. Carboniferous subduction under a low geothermal gradient (~490–570 °C GPa^?1) was followed by exhumation accompanied by heating and Late Carboniferous magma emplacement at a high apparent geothermal gradient (~1200–1450 °C GPa^?1). Exhumation coeval with nappe stacking was closely followed by activity on a ductile strike‐slip shear zone that accommodated magma intrusion and enabled the final exhumation of the Monte Grighini Unit to upper crustal levels. The reconstructed thermo‐mechanical evolution allows a more complete understanding of the Variscan orogenic wedge in central Sardinia. As a result we are able to confirm a diachronous evolution of metamorphic and tectonic events from the inner axial zone to the outer nappe zone, with the Late Variscan low‐P/high‐T metamorphism and crustal anatexis as a common feature across the Sardinian portion of the Variscan orogen.