Exhumation of high-pressure rocks of the Kokchetav massif: facts and models

The exhumation of ultrahigh-pressure (UHP) metamorphic units from depths more than 100-120 km is one of the most intriguing questions in modern petrology and geodynamics. We use the diamondiferous Kumdy-Kol domain in the Kokchetav Massif to show that exhumation models should take into consideration initially high uplift velocities (from 20 down to 6 cm/year) and the absence of the deformation of UHP assemblages. The high rate of exhumation are indicated by ion microprobe (SHRIMP) dating of zircons from diamondiferous rocks and supported by the low degree of nitrogen aggregation in metamorphic diamonds.Diamondiferous rocks in the Kumdy-Kol domain occur as steeply dipping (60°-80°) thin slices (few hundred metres) within granite-gneiss. Using geological, petrological and isotopic-geochemical data, we show that partial melting of diamondiferous metamorphic rocks occurred; a very important factor which has not been taken into account in previous models.Deformation of diamondiferous rocks at Kumdy-Kol is insignificant; diamond inclusions in garnet are often intergrown with mica crystals carrying no traces of deformation. All these facts could be explained by partial melting of metapelites and granitic rocks in the Kumdy-Kol domain. The presence of melt is responsible for an essential reduction of viscosity and a density difference (Δρ) between crustal rocks and mantle material and reduced friction between the upwelling crustal block, the subducting and overriding plates. Besides Δρ, the exhumation rate seems to depend on internal pressure in the subducting continental crustal block which can be regarded as a viscous layer between subducting continental lithosphere and surrounding mantle.We construct different models for the three stages of exhumation: a model similar to “corner flow” for the first superfast exhumation stage, an intermediate stage of extension (most important from structural point of view) and a very low rate of exhumation in final diapir+erosional uplift.     

Exhumation Paths of High-Pressure-Low-Temperature Metamorphic Rocks from the Lycian Nappes and the Menderes Massif (SW Turkey): a Multi-Equilibrium Approach

The Menderes Massif and the overlying Lycian Nappes occupy anextensive area of SW Turkey where high-pressure–low-temperaturemetamorphic rocks occur. Precise retrograde P–T pathsreflecting the tectonic mechanisms responsible for the exhumationof these high-pressure–low-temperature rocks can be constrainedwith multi-equilibrium P–T estimates relying on localequilibria. Whereas a simple isothermal decompression is documentedfor the exhumation of high-pressure parageneses from the southernMenderes Massif, various P–T paths are observed in theoverlying Karaova Formation of the Lycian Nappes. In the uppermostlevels of this unit, far from the contact with the MenderesMassif, all P–T estimates depict cooling decompressionpaths. These high-pressure cooling paths are associated withtop-to-the-NNE movements related to the Akçakaya shearzone, located at the top of the Karaova Formation. This zoneof strain localization is a local intra-nappe contact that wasactive in the early stages of exhumation of the high-pressurerocks. In contrast, at the base of the Karaova Formation, alongthe contact with the Menderes Massif, P–T calculationsshow decompressional heating exhumation paths. These paths areassociated with severe deformation characterized by top-to-the-eastshearing related to a major shear zone (the Gerit shear zone)that reflects late exhumation of high-pressure parageneses underwarmer conditions. KEY WORDS: exhumation; high-pressure–low-temperature metamorphism; multi-equilibrium P–T estimates; Lycian Nappes; Menderes Massif     

Provenance and Magmatic-Metamorphic Evolution of a Variscan Island-Arc Complex: Constraints from U-Pb Dating, Petrology, and Geospeedometry of the Kyffhauser Crystalline Complex, Central Germany

The Kyffhäuser Crystalline Complex, Central Germany, formspart of the Mid-German Crystalline Rise, which is assumed torepresent the Variscan collision zone between the East Avalonianterrane and the Armorican terrane assemblage. High-precisionU–Pb zircon and monazite dating indicates that sedimentaryrocks of the Kyffhäuser Crystalline Complex are youngerthan c. 470 Ma and were intruded by gabbros and diorites between345 ± 4 and 340 ± 1 Ma. These intrusions had magmatictemperatures between 850 and 900°C, and caused a contactmetamorphic overprint of the sediments at P–T conditionsof 690–750°C and 5–7 kbar, corresponding toan intrusion depth of 19–25 km. At 337 ± 1 Ma themagmatic–metamorphic suite was intruded by granites, syenitesand diorites at a shallow crustal level of some 7–11 km.This is inferred from a diorite, and conforms to P–T pathsobtained from the metasediments, indicating a nearly isothermaldecompression from 5–7 to 2–4 kbar at 690–750°C.Subsequently, the metamorphic–magmatic sequence underwentaccelerated cooling to below 400°C, as constrained by garnetgeospeedometry and a previously published K–Ar muscoviteage of 333 ± 7 Ma. With respect to P–T–D–tdata from surrounding units, rapid exhumation of the KCC canbe interpreted to result from NW-directed crustal shorteningduring the Viséan. KEY WORDS: contact metamorphism; U–Pb dating; hornblende; garnet; Mid-German Crystalline Rise; P–T pseudosection     

Petrological evidence for shock‐induced high‐P metamorphism in a gabbro

Microlites (minute spherulitic, dendritic, skeletal, acicular and poikilitic crystals) diagnostic of crystallization in quenched melt or glass in fault rocks have been used to infer fossil earthquakes. High‐P microlites and crystallites are described here in a variably eclogitized gabbro, the wallrock to the coesite‐bearing eclogite breccia at Yangkou in the Chinese Su‐Lu high‐P metamorphic belt. The studied hand specimens are free of discernible shear deformation, although microfractures are not uncommon under the microscope. In the least eclogitized gabbro, the metagabbro, stellate growths of high‐P minerals on the relict igneous minerals are common. Dendritic garnet crystals (<1?5 μm) grew around rutile and/or phengite replacing ilmenite and biotite, respectively. Skeletal garnet also rims broken flakes of igneous biotite and mechanically twinned augite. Radial intergrowths of omphacite and quartz developed around relict igneous orthopyroxene and are rimmed by skeletal or poikilitic garnet where a Ti‐bearing mineral relict is present. Acicular epidote, kyanite and phengite crystallites are randomly distributed in a matrix of Na‐rich plagioclase, forming the pseudomorphs after igneous plagioclase. In the more eclogitized gabbro, the coronitic eclogite located closer to the eclogite breccia, all the igneous minerals broke down into high‐P assemblages. Thick coronas of poikilitic garnet grew between the pseudomorphs after igneous plagioclase and ferromagnesian minerals. The igneous plagioclase is replaced by omphacite crystallites, with minor amounts of phengite and kyanite. Thermodynamic modelling of the plagioclase pseudomorphs shows an increase in P–T in the wallrock from the metagabbro to the coronitic eclogite, and the P–T variation is unrelated to H₂O content. The fluid‐poor pressure overstepping scenario is unsupported both by phase diagram modelling and by whole‐rock chemical data, which show that the various types of eclogitized gabbro are all fairly dry. A large pressure difference of >2 GPa between the metagabbro and the coesite‐bearing eclogites ~20 m apart cannot be explained by the subduction hypothesis because this would require a depth difference of >60 km. The microlites and crystallites are evidence for dynamic crystallization due to rapid cooling because constitutional supercooling was unlikely for the plagioclase pseudomorphs. The lack of annealing of the broken biotite and augite overgrown by strain free skeletal garnet is consistent with a transient high‐P–T event at a low ambient temperature (<300 °C), probably in the crust. Therefore, the eclogitization of the wallrock to the eclogite breccia was also coseismic, as proposed earlier for the eclogite facies fault rocks. The outcrop‐scale P–T variation and the transient nature of the high‐P–T event are inconsistent with the other existing tectonic models for high‐P metamorphism. The fact that the less refractory but denser biotite is largely preserved while the more refractory but less dense plagioclase broke down completely into high‐P microlite assemblages in the metagabbro indicates a significant rise in pressure rather than temperature. Given that the metamorphic temperatures are far below the melting temperatures of most of the gabbroic minerals under fluid‐absent conditions, stress‐induced amorphization appears to be the more likely mechanism of the coseismic high‐P metamorphism.     

An Eocene/Oligocene blueschist‐/greenschist facies P–T loop from the Cycladic Blueschist Unit on Naxos Island,Greece: Deformation‐related re‐equilibration vs. thermal relaxation

Geothermobarometric and geochronological work indicates a complete Eocene/early Oligocene blueschist/greenschist facies metamorphic cycle of the Cycladic Blueschist Unit on Naxos Island in the Aegean Sea region. Using the average pressure–temperature (P–T) method of thermocalc coupled with detailed textural work, we separate an early blueschist facies event at 576 ± 16 to 619 ± 32°C and 15.5 ± 0.5 to 16.3 ± 0.9 kbar from a subsequent greenschist facies overprint at 384 ± 30°C and 3.8 ± 1.1 kbar. Multi‐mineral Rb–Sr isochron dating yields crystallization ages for near peak‐pressure blueschist facies assemblages between 40.5 ± 1.0 and 38.3 ± 0.5 Ma. The greenschist facies overprint commonly did not result in complete resetting of age signatures. Maximum ages for the end of greenschist facies reworking, obtained from disequilibrium patterns, cluster near c. 32 Ma, with one sample showing rejuvenation at c. 27 Ma. We conclude that the high‐P rocks from south Naxos were exhumed to upper mid‐crustal levels in the late Eocene and early Oligocene at rates of 7.4 ± 4.6 km/Ma, completing a full blueschist‐/greenschist facies metamorphic cycle soon after subduction within c. 8 Ma. The greenschist facies overprint of the blueschist facies rocks from south Naxos resulted from rapid exhumation and associated deformation/fluid‐controlled metamorphic re‐equilibration, and is unrelated to the strong high‐T metamorphism associated with the Miocene formation of the Naxos migmatite dome. It follows that the Miocene thermal overprint had no impact on rock textures or Sr isotopic signatures, and that the rocks of south Naxos underwent three metamorphic events, one more than hitherto envisaged.     

Polymetamorphism in garnet micaschists of the Saualpe Eclogite Unit (Eastern Alps,Austria), resolved by automated SEM methods and EMP–Th–U–Pb monazite dating

Polymetamorphic garnet micaschists from the Austroalpine Saualpe Eclogite Unit (Kärnten, Austria, Eastern Alps) display complex microstructural and mineral–chemical relationships. Automated scanning electron microscopy routines with energy dispersive X‐ray (EDX) spectral mapping were applied for monazite detection and garnet mineral–chemical characterization. When the Fe, Mg, Mn and Ca element wt% compositions are used as generic labels for garnet EDX spectra, complex zonations and porphyroblast generations can be resolved in complete thin sections for selective electron‐microprobe analyses. Two garnet porphyroblast generations and diverse monazite age populations have been revealed in low‐Ca and high‐Al‐metapelites. Garnet 1 has decreasing Mn, constant Ca and significantly increasing Mg from cores to rims. Geothermobarometry of garnet 1 assemblages signals a crystallization along a M1 prograde metamorphism at ~650 °C/6–8 kbar. Sporadic monazite 1 crystallization started at c. 320 Ma. Subsequent pervasive 300–250 Ma high‐Y and high‐Gd monazite 1 formation during decompression coincided with the intrusion of Permian and Early Triassic pegmatites. Monazite 1 crystallized along the margin of garnet 1. Coronas of apatite and allanite around the large 320–250 Ma monazite signal a retrogressive stage. These microstructures suggest a Carboniferous‐to‐Early‐Permian age for the prograde M1 event with garnet 1. Such a M1 event at an intermediate‐P/T gradient has not yet been described from the Saualpe, and preceded a Permo‐Triassic low‐P stage. The M2 event with garnet 2 postdates the corona formation around Permian monazite. Garnet 2 displays first increasing X_Ca at decreasing X_Mg, then increasing X_Ca and X_Mg, and finally decreasing X_Ca with increasing X_Mg, always at high Ca and Mg, and low Mn. This records a P–T evolution which passed through eclogite facies conditions and reached maximum temperatures at ~750 °C/14 kbar during decompression‐heating. A monazite 2 population (94–86 Ma) with lower Y and Gd contents crystallized at decreasing pressure during the Cretaceous (Eo‐Alpine) metamorphism M2 at a high‐P/T gradient. The Saualpe Eclogite Unit underwent two distinct clockwise metamorphic cycles at different P–T conditions, related to continental collisions under different thermal regimes. This led to a characteristic distribution pattern of monazite ages in this unit which is different from other Austroalpine basement areas.     

Diamond in metasedimentary crustal rocks from Pohorje,Eastern Alps: a window to deep continental subduction

We report the first finding of diamond and moissanite in metasedimentary crustal rocks of Pohorje Mountains (Slovenia) in the Austroalpine ultrahigh‐pressure (UHP) metamorphic terrane of the Eastern Alps. Microscopic observations and Raman spectroscopy show that diamond occurs in situ as inclusions in garnet, being heterogeneously distributed. Under the optical microscope, diamond‐bearing inclusions are of cuboidal to rounded shape and of pinkish, yellow to brownish colour. The Raman spectra of the investigated diamond show a sharp, first order peak of sp³‐bonded carbon, in most cases centred between 1332 and 1330 cm^?1, with a full width at half maximum between 3 and 5 cm^?1. Several spectra show Raman bands typical for disordered graphitic (sp²‐bonded) carbon. Detailed observations show that diamond occurs either as a monomineralic, single‐crystal inclusion or it is associated with SiC (moissanite), CO₂ and CH₄ in polyphase inclusions_. This rare record of diamond occurring with moissanite as fluid‐inclusion daughter minerals implies the crystallization of diamond and moissanite from a supercritical fluid at reducing conditions. Thermodynamic modelling suggests that diamond‐bearing gneisses attained P–T conditions of ≥3.5 GPa and 800–850 °C, similar to eclogites and garnet peridotites. We argue that diamond formed when carbonaceous sediment underwent UHP metamorphism at mantle depth exceeding 100 km during continental subduction in the Late Cretaceous (c. 95–92 Ma). The finding of diamond confirms UHP metamorphism in the Pohorje Mountains, the most deeply subducted part of Austroalpine units.     

Experimental determination of REE partition coefficients between zircon,garnet and melt: a key to understanding high‐T crustal processes

The partitioning of rare earth elements (REE) between zircon, garnet and silicate melt was determined using synthetic compositions designed to represent partial melts formed in the lower crust during anatexis. The experiments, performed using internally heated gas pressure vessels at 7 kbar and 900–1000 °C, represent equilibrium partitioning of the middle to heavy REE between zircon and garnet during high‐grade metamorphism in the mid to lower crust. The D_REE (zircon/garnet) values show a clear partitioning signature close to unity from Gd to Lu. Because the light REE have low concentrations in both minerals, values are calculated from strain modelling of the middle to heavy REE experimental data; these results show that zircon is favoured over garnet by up to two orders of magnitude. The resulting general concave‐up shape to the partitioning pattern across the REE reflects the preferential incorporation of middle REE into garnet, with D_Gd (zircon/garnet) ranging from 0.7 to 1.1, D_Ho (zircon/garnet) from 0.4 to 0.7 and D_Lu (zircon/garnet) from 0.6 to 1.3. There is no significant temperature dependence in the zircon–garnet REE partitioning at 7 kbar and 900–1000 °C, suggesting that these values can be applied to the interpretation of zircon–garnet equilibrium and timing relationships in the ultrahigh‐T metamorphism of low‐Ca pelitic and aluminous granulites.     

Low‐P and high‐T metamorphism of basalts: Insights from the Sudbury impact melt sheet aureole and thermodynamic modelling

Low‐pressure and high‐temperature (LP–HT) metamorphism of basaltic rocks, which occurs globally and throughout geological time, is rarely constrained by forward phase equilibrium modelling, yet such calculations provide valuable supplementary thermometric information and constraints on anatexis that are not possible to obtain from conventional thermometry. Metabasalts along the southern margin of the Sudbury Igneous Complex (SIC) record evidence of high‐grade contact metamorphism involving partial melting and melt segregation. Peak metamorphic temperatures reached at least ~925°C at ~1–3 kbar near the SIC contact. Preservation of the peak mineral assemblage indicates that most of the generated melt escaped from these rocks leaving a residuum characterized by a plagioclase–orthopyroxene–clinopyroxene–ilmenite‐magnetite±melt assemblage. Peak temperatures reached ~875°C up to 500 m from the SIC lower contact, which marks the transition to metabasalts that only experienced incipient partial melting without melt loss. Metabasalts ~500 to 750 m from the SIC contact are characterized by a similar two‐pyroxene mineral assemblage, but typically contain abundant hornblende that overgrew clino‐ and orthopyroxene along an isobaric cooling path. Metabasalts ~750 to 1,000 m from the SIC contact are characterized by a hornblende–plagioclase–quartz–ilmenite assemblage indicating temperatures up to ~680°C. Mass balance and phase equilibria calculations indicate that anatexis resulted in 10–20% melt generation in the inner ~500 m of the aureole, with even higher degrees of melting towards the contact. Comparison of multiple models, experiments, and natural samples indicates that modelling in the Na₂O–CaO–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O₂ (NCFMASHTO) system results in the most reliable predictions for the temperature of the solidus. Incorporation of K₂O in the most recent amphibole solution model now successfully predicts dehydration melting by the coexistence of high‐Ca amphibole and silicate melt at relatively low pressures (~1.5 kbar). However, inclusion of K₂O as a system component results in prediction of the solidus at too low a temperature. Although there are discrepancies between modelling predictions and experimental results, this study demonstrates that the pseudosection approach to mafic rocks is an invaluable tool to constrain metamorphic processes at LP–HT conditions.