Uncoupled U/Pb and REE response in zircon during the transformation of eclogite to mafic and intermediate granulite (Blanský les,Bohemian Massif)

An eclogite–mafic granulite occurs as a rare boudin within a felsic kyanite–K‐feldspar granulite in a low‐strain zone. Its boundary is marked by significant metasomatism–diffusional gain of potassium at the centimetre‐scale, and probable infiltration of felsic melt on a larger scale. This converted the eclogite–mafic granulite into an intermediate‐composition, ternary‐feldspar‐bearing granulite. Based on inclusions in garnet, the peak P–T conditions of the original eclogite are 18 kbar at 850–950 °C, with later matrix re‐equilibration at 12 kbar and 950 °C. Four samples from the transition of the eclogite–mafic granulite through to the intermediate granulite were studied. In the eclogite, REE patterns in the garnet core show no Eu anomaly, compatible with crystallization in the absence of plagioclase and consistent with eclogite facies conditions. Towards the rim of garnet, LREE decrease, and a weak negative Eu anomaly appears, reflecting passage into HP granulite facies conditions with plagioclase present. The rims of garnet next to ternary feldspar in the intermediate granulite show the lowest LREE and deepest Eu anomalies. Zircon from the four samples was analysed by LASS (laser ablation–split‐stream inductively coupled plasma–mass spectrometry). It shows U–Pb ages from 404 ± 4.0 to 331 ± 3.3 Ma, with a peak at 340 ± 4.0 Ma corresponding to the likely exhumation of the rocks to 12 kbar. Older ages from zircon with steep HREE patterns indicate the minimum age of the protolith, and ages <360 ± 4.0 Ma are interpreted to correspond to the eclogite facies metamorphism. Only some zircon grains ≤350 ± 4.0 Ma have flat HREE patterns, suggesting that these are primarily modified protolith grains, rather than new zircon crystallized in the eclogite‐ or granulite facies. The metasomatic processes that converted the eclogite–mafic granulite to an intermediate granulite may have facilitated zircon modification as zircon in the intermediate granulite has flat HREE and ages of 340 ± 4.0 Ma. The difference between the oldest and youngest ages with flat REE patterns indicates a 16 ± 5.6 Ma period of zircon modification in the presence of garnet.     

The tectono‐metamorphic evolution of the very low‐grade hangingwall constrains two‐stage gneiss dome formation in the Montagne Noire (Southern France)

The Montagne Noire in the southernmost French Massif Central is made of an ENE‐elongated gneiss dome flanked by Palaeozoic sedimentary rocks. The tectonic evolution of the gneiss dome has generated controversy for more than half a century. As a result, a multitude of models have been proposed that invoke various tectonic regimes and exhumation mechanisms. Most of these models are based on data from the gneiss dome itself. Here, new constraints on the dome evolution are provided based on a combination of very low‐grade petrology, K–Ar geochronology, field mapping and structural analysis of the Palaeozoic western Mont Peyroux and Faugères units, which constitute part of the southern hangingwall of the dome. It is shown that southward‐directed Variscan nappe‐thrusting (D₁) and a related medium‐P metamorphism (M₁) are only preserved in the area furthest away from the gneiss dome. The regionally dominant pervasive tectono‐metamorphic event D₂/M₂ largely transposes D₁ structures, comprises a higher metamorphic thermal gradient than M₁ (transition low‐P and medium‐P metamorphic facies series) and affected the rocks between c. 309 and 300 Ma, post‐dating D₁/M₁ by more than 20 Ma. D₂‐related fabrics are refolded by D₃, which in its turn, is followed by dextral‐normal shearing along the basal shear zone of both units at c. 297 Ma. In the western Mont Peyroux and Faugères units, D₂/M₂ is largely synchronous with shearing along the southern dome margin between c. 311 and 303 Ma, facilitating the emplacement of the gneiss dome into the upper crust. D₂/M₂ also overlaps in time with granitic magmatism and migmatization in the Zone Axiale between c. 314 and 306 Ma, and a related low‐P/high‐T metamorphism at c. 308 Ma. The shearing that accompanied the exhumation of the dome therefore was synchronous with a peak in temperature expressed by migmatization and intrusion of melts within the dome, and also with the peak of metamorphism in the hangingwall. Both, the intensity of D₂ fabrics and the M₂ metamorphic grade within the hangingwall, decrease away from the gneiss dome, with grades ranging from the anchizone–epizone boundary to the diagenetic zone. The related zonation of the pre‐D₃ metamorphic field gradients paralleled the dome. These observations indicate that D₂/M₂ is controlled by the exhumation of the Zone Axiale, and suggest a coherent kinematic between the different crustal levels at some time during D₂/M₂. Based on integration of these findings with regional geological constraints, a two‐stage exhumation of the gneiss dome is proposed: during a first stage between c. 316 and 300 Ma dome emplacement into the upper crust was controlled by dextral shear zones arranged in a pull‐apart‐like geometry. The second stage from 300 Ma onwards was characterized by northeast to northward extension, with exhumation accommodated by north‐dipping detachments and hangingwall basin formation along the northeastern dome margin.     

Metamorphic petrology of a high‐T/low‐P granulite terrane (Damara belt,Namibia) – Constraints from pseudosection modelling and high‐precision Lu–Hf garnet‐whole rock dating

Migmatites comprise a minor volume of the high‐grade part of the Damara orogen of Namibia that is dominated by granite complexes and intercalated metasedimentary units. Migmatites of the Southern Central Zone of the Damara orogen consist of melanosomes with garnet+cordierite+biotite+K‐feldspar, and leucosomes, which are sometimes garnet‐ and cordierite‐bearing. Field evidence, petrographic observations, and pseudosection modelling suggest that, in contrast to other areas where intrusion of granitic magmas is more important, in situ partial melting of metasedimentary units was the main migmatite generation processes. Pseudosection modelling and thermobarometric calculations consistently indicate that the peak‐metamorphic grade throughout the area is in the granulite facies (~5 kbar at ~800°C). Cordierite coronas around garnet suggest some decompression from peak‐metamorphic conditions and rare andalusite records late, near‐isobaric cooling to <650°C at low pressures of ~3 kbar. The inferred clockwise P–T path is consistent with minor crustal thickening through continent–continent collision followed by limited post‐collisional exhumation and suggests that the granulite facies terrane of the Southern Central Zone of the Damara orogen formed initially in a metamorphic field gradient of ~35–40°C/km at medium pressures. New high‐precision Lu–Hf garnet‐whole rock dates are 530 ± 13 Ma, 522.0 ± 0.8 Ma, 520.8 ± 3.6 Ma, and 500.3 ± 4.3 Ma for the migmatites that record temperatures of ~800°C. This indicates that high‐grade metamorphism lasted for c. 20–30 Ma, which is compatible with previous estimates using Sm–Nd garnet‐whole rock systematics. In previous studies on Damara orogen migmatites where both Sm–Nd and Lu–Hf chronometers have been applied, the dates (c. 520–510 Ma) agree within their small uncertainties (0.6–0.8% for Sm–Nd and 0.1–0.2% for Lu–Hf). This implies rapid cooling after high‐grade conditions and, by implication, rapid exhumation at that time. The cause of the high geothermal gradient inferred from the metamorphic conditions is unknown but likely requires some extra heat that was probably added by intrusion of magmas from the lithospheric mantle, i.e., syenites that have been recently re‐dated at c. 545 Ma. Some granites derived from the lower crust at c. 545 Ma are the outcome rather than the cause of high‐T metamorphism. In addition, high contents of heat‐producing elements K, Th, and U may have raised peak temperatures by 150–200°C at the base of the crust, resulting in the widespread melting of fertile crustal rocks. The continuous gradation from centimetre‐scale leucosomes to decametre‐scale leucogranite sheets within the high‐grade metamorphic zone suggests that leucosome lenses coalesced to form larger bodies of anatectic leucogranites, thereby documenting a link between high‐grade regional metamorphism and Pan‐African magmatism. In view of the close association of the studied high‐T migmatites with hundreds of synmetamorphic high‐T granites that invaded the terrane as metre‐ to decametre‐wide sills and dykes, we postulate that crystallization of felsic lower crustal magma is, at least partly, responsible for heat supply. Late‐stage isobaric cooling of these granites may explain the occurrence of andalusite in some samples.     

Tectonometamorphic evolution of the Atbashi high‐P units (Kyrgyz CAOB,Tien Shan): Implications for the closure of the Turkestan Ocean and continental subduction–exhumation of the South Kazakh continental margin

The South Tien Shan (STS) belt results from the last collision event in the western Central Asian Orogenic Belt (CAOB). Understanding its formation is of prime importance in the general framework of the CAOB. The Atbashi Range preserves high‐P (HP) rocks along the STS suture, but still, its global metamorphic evolution remains poorly constrained. Several HP units have been identified: (a) a HP tectonic mélange including boudins of mafic eclogites in a sedimentary matrix, (b) a large (>100 km long) high‐P metasedimentary unit (HPMU) and (c) a lower blueschist facies accretionary prism. Raman Spectroscopy on carbonaceous material combined with phengite and chlorite multiequilibria and isochemical phase diagram modelling indicates that the HPMU recorded homogeneous P–T conditions of 23–25 kbar and 560–570°C along the whole unit. ⁴⁰Ar/³⁹Ar dating on phengite from the HPMU ranges between 328 and 319 Ma at regional scale. These ages are interpreted as (re‐) crystallization ages of phengite during T_max conditions at a pressure range of 20–25 kbar. Thermobarometry on samples from the HP tectonic mélange provides similar metamorphic peak conditions. Thermobarometry on the blueschist to lower greenschist facies accretionary prism indicates that it underwent P–T conditions of 5–6 kbar and 290–340°C, highlighting a 17–20 kbar pressure gap between the HPMU‐tectonic mélange units and the accretionary prism. Comparison with available geochronological data suggests a very short time span between the prograde path (340 Ma), HP metamorphic peak (330 Ma), the T_max (328–319 Ma) and the final exhumation of the HPMU (303–295 Ma). Extrusion of the HPMU, accommodated by a basal thrust and an upper detachment, was driven by buoyant forces from 70–75 km up to 60 km depth, which directly followed continental subduction and detachment of the HPMU. At crustal depths, extrusion was controlled by collisional tectonics up to shallow levels. Lithological homogeneity of the HPMU and its continental‐derived character from the North Tien Shan suggest this unit corresponds to the hyper‐extended continental margin of the Kazakh continent, subducted southward below the north continental active margin of the Tarim craton. Integration of the available geological data allows us to propose a general geodynamic scenario for Tien Shan during the Carboniferous with a combination of (a) N‐dipping subduction below the Kazakh margin of Middle Tien Shan until 390–340 Ma and (b) S‐dipping subduction of remaining Turkestan marginal basins between 340 and 320 Ma.     

Protracted garnet growth in high‐P eclogite: constraints from multiple geochronology and P–T pseudosection

Understanding convergent margin processes requires determination of the onset and the termination of subduction, the duration of subduction‐zone metamorphism, and the subduction zone polarity. Garnet growth and intracrystalline zonation can be used to constrain the timing, duration and kinetics of tectonometamorphic processes. An eclogite from the Huwan shear zone in the Hong'an orogen was investigated with combined pseudosection analysis and multiple geochronologies. The pseudosection analysis illustrates that garnet growth is continuous and along an early near‐isothermal trajectory followed by a near‐isobaric heating path from 1.9 GPa/500 °C to 2.4 GPa/575 °C and subsequent near‐isothermal decompression. ⁴⁰Ar/³⁹Ar dating of an amphibole inclusion in garnet from the eclogite yielded an age of 310 ± 5 Ma, which is consistent with a U–Pb age of 305 ± 3 Ma for the metamorphic zircon within uncertainty. Garnet core and rim material produced Lu–Hf ages of 296.9 ± 3.8 and 256.9 ± 3.9 Ma respectively; the latter is consistent with its Sm–Nd age of 254.3 ± 4.6 Ma for the same aliquots. Similarly, limited zircon U–Pb ages of c. 257 Ma were obtained in zircon rims with garnet inclusions. These ages were interpreted to bracket the period of garnet growth and the difference of up to c. 40 Ma is best explained by protracted garnet growth. We propose that the rocks represent detachment of part of the downgoing slab and remained free of significant compression/decompression or heating/cooling close to the subduction channel, most likely underplating the mantle wedge, for a long time. These rocks were incorporated into the following subduction channel due to the successive entry of the buoyant materials, and exhumed at some time later than c. 254 Ma. The increasing observations of protracted garnet growth and long‐lived subduction in various orogens worldwide demand more sophisticated geodynamic models.     

Neoproterozoic evolution and Cambrian reworking of ultrahigh temperature granulites in the Eastern Ghats Province,India

The time‐scales and P–T conditions recorded by granulite facies metamorphic rocks permit inferences about the geodynamic regime in which they formed. Two compositionally heterogeneous cordierite–spinel‐bearing granulites from Vizianagaram, Eastern Ghats Province (EGP), India, were investigated to provide P–T–time constraints using petrography, phase equilibrium modelling, U–Pb geochronology, the rare earth element composition of zircon and monazite, and Ti‐in‐zircon thermometry. These ultrahigh temperature (UHT) granulites preserve discrete compositional layering in which different inferred peak assemblages are developed, including layers bearing garnet–sillimanite–spinel, and others bearing orthopyroxene–sillimanite–spinel. These mineral associations cannot be reproduced by phase equilibrium modelling of whole‐rock compositions, indicating that the samples became domainal on a scale less than that of a thin section, even at UHT conditions. Calculation of the P–T stability fields for six compositional domains within which the main rock‐forming minerals are considered to have attained equilibrium suggests peak metamorphic conditions of ~6.8–8.3 kbar at ~1,000°C. In most of these domains, the subsequent evolution resulted in the growth of cordierite and final crystallization of melt at an elevated (residual) H₂O‐undersaturated solidus, consistent with <1 kbar of decompression. Concordant U–Pb ages obtained by SHRIMP from zircon (spread 1,050–800 Ma) and monazite (spread 950–800 Ma) demonstrate that crystallization of these minerals occurred during an interval of c. 250 Ma. By combining LA‐ICP‐MS U–Pb zircon ages with Ti‐in‐zircon temperatures from the same analysis sites, we show that the crust may have remained above 900°C for a minimum of c. 120 Ma between c. 1,000 and c. 880 Ma. Overall, the results suggest that, in the interval 1,050 to 800 Ma, the evolution of the Vizianagaram granulites culminated with UHT conditions from c. 1,000 Ma to c. 880 Ma, associated with minor decompression, before further zircon crystallization at c. 880–800 Ma during cooling to the solidus. However, these rocks are adjacent to the Paderu–Anantagiri–Salur crustal block to the NW that experienced counterclockwise P–T–t paths, and records similar UHT peak metamorphic conditions (7–8 kbar, ~950°C) followed by near‐isobaric cooling, and has a similar chronology during the Neoproterozoic. The limited decompression inferred at Vizianagaram may be explained by partial exhumation due to thrusting of this crustal block over the adjacent Paderu–Anantagiri–Salur crustal block. The residual granulites in both blocks have high concentrations of heat‐producing elements and likely remained hot at mid‐crustal depths throughout a period of relative tectonic quiescence in the interval 800–550 Ma. During the Cambrian Period, the EGP was located in the hinterland of the Denman–Pinjarra–Prydz orogen. A later concordant population of zircon dated at 511 ± 6 Ma records crystallization at temperatures of ~810°C. This age may record a low‐degree of melting due to limited influx of fluid into hot, weak crust in response to convergence of the Crohn craton with a composite orogenic hinterland comprising the Rayner terrane, EGP, and cratonic India.     

Post‐Variscan exhumation of the Central Anti‐Atlas (Morocco) constrained by zircon and apatite fission‐track thermochronology

The origin of the Anti‐Atlas relief is one of the currently debated issues of Moroccan geology. To constrain the post‐Variscan evolution of the Central Anti‐Atlas, we collected nine samples from the Precambrian basement of the Bou Azzer‐El Graara inlier for zircon and apatite fission‐track thermochronology. Zircon ages cluster between 340 ± 20 and 306 ± 20 Ma, whereas apatite ages range from 171 ± 7 Ma to 133 ± 5 Ma. Zircon ages reflect the thermal effect of the Variscan orogeny (tectonic thickening of the ca. 7 km‐thick Paleozoic series), likely enhanced by fluid advection. Apatite ages record a complex Mesozoic–Cenozoic exhumation history. Track length modelling yields evidence that, (i) the Precambrian basement was still buried at ca. 5 km depth by Permian times, (ii) the Central Anti‐Atlas was subjected to (erosional) exhumation during the Triassic‐Early Cretaceous, then buried beneath ca. 1.5 km‐thick Cretaceous‐Paleogene deposits, (iii) final exhumation took place during the Neogene, contemporaneously with that of the High Atlas.     

Preconditioned IDR(s) iterative solver for non‐symmetric linear system associated with FEM analysis of shallow foundation

Non‐associated flow rule is essential when the popular Mohr–Coulomb model is used to model nonlinear behavior of soil. The global tangent stiffness matrix in nonlinear finite element analysis becomes non‐symmetric when this non‐associated flow rule is applied. Efficient solution of this large‐scale non‐symmetric linear system is of practical importance. The standard Krylov solver for a non‐symmetric solver is Bi‐CGSTAB. The Induced Dimension Reduction [IDR(s)] solver was proposed in the scientific computing literature relatively recently. Numerical studies of a drained strip footing problem on homogenous soil layer show that IDR(s = 6) is more efficient than Bi‐CGSTAB when the preconditioner is the incomplete factorization with zero fill‐in of global stiffness matrix K_ep (ILU(0)‐K_ep). Iteration time is reduced by 40% by using IDR(s = 6) with ILU(0)‐K_ep. To further reduce computational cost, the global stiffness matrix K_ep is divided into two parts. The first part is the linear elastic stiffness matrix K_e, which is formed only once at the beginning of solution step. The second part is a low‐rank matrix Δ, which is re‐formed at each Newton–Raphson iteration. Numerical studies show that IDR(s = 6) with this ILU(0)‐K_e preconditioner is more time effective than IDR(s = 6) with ILU(0)‐K_ep when the percentage of yielded Gauss points in the mesh is less than 15%. The total computation time is reduced by 60% when all the recommended optimizing methods are used. Copyright © 2013 John Wiley & Sons, Ltd.     

Stable staurolite–cordierite assemblages in K‐poor metapelitic schists in Aston and Hospitalet gneiss domes of the central Pyrenees (France,Andorra)

Staurolite–cordierite assemblages are common in mica schists of the Aston and Hospitalet gneiss domes of the central Axial Zone, Pyrenees (France, Andorra). Within a 200 m wide zone, staurolite, cordierite and andalusite porphyroblasts contain inclusion trails that preserve the same stage of development of a crenulation cleavage, strongly suggesting that all three phases are contemporaneous. Their syntectonic growth occurred during a short period at the beginning of the formation of the dominant schistosity (S₂) of the domes. Staurolite and cordierite touching each other further indicates an equilibrium relationship. Whole‐rock analyses show that some staurolite–cordierite schists are depleted in K₂O compared to post‐Archean shales (PAAS) and amphibolite facies pelites. Analysis of the st‐crd paragenesis in K‐poor schists without muscovite using KFMASH and MnNCKFMASH petrogentic grids, pseudosections and AFM compatibility diagrams predicts stable conditions at pressures of ~3.5 kbar at 575 °C. For metapelites with intermediate X_Mg values (0.7 >  X_Mg >0.48) a ‘muscovite‐out window’ exists from 550–650 °C at 3.5 kbar in the KFMASH system. Conventional thermobarometry (GB‐GASP, AvT‐AvP) and petrogenetic grids show an isobaric P–T path to peak temperatures of ~650 °C, supported by the presence of sillimanite‐K‐feldspar gneiss and migmatites. LP‐HT metamorphism in the Aston dome is related to early Carboniferous (c. 339 Ma) granitic intrusions into the dome core. As metamorphism is directly linked with the formation of the main S₂ schistosity, the temporal relations demonstrated in this study conflict with previous studies which constrained LP‐HT metamorphism and the development of flat‐lying schistosity to the late Carboniferous (315–305 Ma) – at least in the eastern Axial Zone.     

Monazite in a Variscan mylonitic paragneiss from the Münchberg Metamorphic Complex (NE Bavaria) records Cadomian protolith ages

Small oval‐shaped, unshielded monazite grains found in a Variscan garnet–muscovite‐bearing mylonitic paragneiss from the Liegendserie unit of the Münchberg Metamorphic Complex in the northwestern Bohemian Massif, central Europe, yield only pre‐Variscan ages. These ages, determined with the electron microprobe, have maxima at c. 545, 520 and 495 Ma and two side‐maxima at 455 and 575 Ma, and are comparable with previously determined ages of detrital zircon reported from paragneisses elsewhere in the NW Bohemian Massif. The pressure (P)–temperature (T) history of this mylonitic paragneiss, determined from contoured P–T pseudosections, involved an initial stage at 6 kbar/600 °C, reaching peak P–T conditions of 12.5 kbar/670 °C with partial melting, followed by mylonitization and retrogression to 9 kbar/610 °C. The monazite, representing detrital grains derived from igneous rocks of a Cadomian provenance between 575 and 455 Ma, has survived these Variscan metamorphic/deformational events unchanged because this mineral has probably never been outside its P–T stability field during metamorphism.     

Geochronological and petrological constraints from the evolution in the Saxon Granulite Massif,Germany, on the Variscan continental collision orogeny

Controversy over the plate tectonic affinity and evolution of the Saxon granulites in a two‐ or multi‐plate setting during inter‐ or intracontinental collision makes the Saxon Granulite Massif a key area for the understanding of the Palaeozoic Variscan orogeny. The massif is a large dome structure in which tectonic slivers of metapelite and metaophiolite units occur along a shear zone separating a diapir‐like body of high‐P granulite below from low‐P metasedimentary rocks above. Each of the upper structural units records a different metamorphic evolution until its assembly with the exhuming granulite body. New age and petrologic data suggest that the metaophiolites developed from early Cambrian protoliths during high‐P amphibolite facies metamorphism in the mid‐ to late‐Devonian and thermal overprinting by the exhuming hot granulite body in the early Carboniferous. A correlation of new Ar–Ar biotite ages with published P–T–t data for the granulites implies that exhumation and cooling of the granulite body occurred at average rates of ~8 mm/year and ~80°C/Ma, with a drop in exhumation rate from ~20 to ~2.5 mm/year and a slight rise in cooling rate between early and late stages of exhumation. A time lag of c. 2 Ma between cooling through the closure temperatures for argon diffusion in hornblende and biotite indicates a cooling rate of 90°C/Ma when all units had assembled into the massif. A two‐plate model of the Variscan orogeny in which the above evolution is related to a short‐lived intra‐Gondwana subduction zone conflicts with the oceanic affinity of the metaophiolites and the timescale of c. 50 Ma for the metamorphism. Alternative models focusing on the internal Variscan belt assume distinctly different material paths through the lower or upper crust for strikingly similar granulite massifs. An earlier proposed model of bilateral subduction below the internal Variscan belt may solve this problem.     

Influence of dissolution/reprecipitation reactions on metamorphic greenschist to amphibolite facies mica 40Ar/39Ar ages in the Longmen Shan (eastern Tibet)

Linking ages to metamorphic stages in rocks that have experienced low‐ to medium‐grade metamorphism can be particularly tricky due to the rarity of index minerals and the preservation of mineral or compositional relicts. The timing of metamorphism and the Mesozoic exhumation of the metasedimentary units and crystalline basement that form the internal part of the Longmen Shan (eastern Tibet, Sichuan, China), are, for these reasons, still largely unconstrained, but crucial for understanding the regional tectonic evolution of eastern Tibet. In situ core‐rim ⁴⁰Ar/³⁹Ar biotite and U–Th/Pb allanite data show that amphibolite facies conditions (~10–11 kbar, 530°C to 6–7 kbar, 580°C) were reached at 210–180 Ma and that biotite records crystallization, rather than cooling, ages. These conditions are mainly recorded in the metasedimentary cover. The ⁴⁰Ar/³⁹Ar ages obtained from matrix muscovite that partially re‐equilibrated during the post peak‐P metamorphic history comprise a mixture of ages between that of early prograde muscovite relicts and the timing of late muscovite recrystallization at c. 140–120 Ma. This event marks a previously poorly documented greenschist facies metamorphic overprint. This latest stage is also recorded in the crystalline basement, and defines the timing of the greenschist overprint (7 ± 1 kbar, 370 ± 35°C). Numerical models of Ar diffusion show that the difference between ⁴⁰Ar/³⁹Ar biotite and muscovite ages cannot be explained by a slow and protracted cooling in an open system. The model and petrological results rather suggest that biotite and muscovite experienced different Ar retention and resetting histories. The Ar record in mica of the studied low‐ to medium‐grade rocks seems to be mainly controlled by dissolution–reprecipitation processes rather than by diffusive loss, and by different microstructural positions in the sample. Together, our data show that the metasedimentary cover was thickened and cooled independently from the basement prior to c. 140 Ma (with a relatively fast cooling at 4.5 ± 0.5°C/Ma between 185 and 140 Ma). Since the Lower Cretaceous, the metasedimentary cover and the crystalline basement experienced a coherent history during which both were partially exhumed. The Mesozoic history of the Eastern border of the Tibetan plateau is therefore complex and polyphase, and the basement was actively involved at least since the Early Cretaceous, changing our perspective on the contribution of the Cenozoic geology.     

Protracted thrusting followed by late rapid cooling of the Greater Himalayan Sequence,Annapurna Himalaya,Central Nepal: Insights from titanite petrochronology

The Greater Himalayan Sequence (GHS) has commonly been treated as a large coherently deforming high‐grade tectonic package, exhumed primarily by simultaneous thrust‐ and normal‐sense shearing on its bounding structures and erosion along its frontal exposure. A new paradigm, developed over the past decade, suggests that the GHS is not a single high‐grade lithotectonic unit, but consists of in‐sequence thrust sheets. In this study, we examine this concept in central Nepal by integrating temperature–time (T–t) paths, based on coupled Zr‐in‐titanite thermometry and U–Pb geochronology for upper GHS calcsilicates, with traditional thermobarometry, textural relationships and field mapping. Peak Zr‐in‐titanite temperatures are 760–850°C at 10–13 kbar, and U–Pb ages of titanite range from c. 30 to c. 15 Ma. Sector zoning of Zr and distribution of U–Pb ages within titanite suggest that diffusion rates of Zr and Pb are slower than experimentally determined rates, and these systems remain unaffected into the lower granulite facies. Two types of T–t paths occur across the Chame Shear Zone (CSZ). Between c. 25 and 17–16 Ma, hangingwall rocks cool at rates of 1–10°C/Ma, while footwall rocks heat at rates of 1–10°C/Ma. Over the same interval, temperatures increase structurally upwards through the hangingwall, but by 17–16 Ma temperatures converge. In contrast, temperatures decrease upwards in footwall rocks at all times. While the footwall is interpreted as an intact, structurally upright section, the thermometric inversion within the hangingwall suggests thrusting of hotter rocks over colder from c. 25 to c. 17–16 Ma. Retrograde hydration that is restricted to the hangingwall, and a lithological repetition of orthogneiss are consistent with thrust‐sense shear on the CSZ. The CSZ is structurally higher than previously identified intra‐GHS thrusts in central Nepal, and thrusting duration was 3–6 Ma longer than proposed for other intra‐GHS thrusts in this region. Cooling rates for both the hangingwall and footwall of the CSZ are comparable to or faster than rates for other intra‐GHS thrust sheets in Nepal. The overlap in high‐T titanite U–Pb ages and previously published muscovite ⁴⁰Ar/³⁹Ar cooling ages imply cooling rates for the hangingwall of ≥200°C/Ma after thrusting. Causes of rapid cooling include passive exhumation driven by a combination of duplexing in the Lesser Himalayan Sequence, and juxtaposition of cooler rocks on top of the GHS by the STDS. Normal‐sense displacement does not appear to affect T–t paths for rocks immediately below the STDS prior to 17–16 Ma.     

Provenance of the HP–HT subducted margin in the Variscan belt (Cabo Ortegal Complex,NW Iberian Massif)

The Variscan Upper Allochthon is a continental‐affinity terrane that recorded a Cambrian–Ediacaran magmatic arc generation, a subsequent transition to a passive margin, and a collision‐related high‐P metamorphism during the Devonian–Carboniferous amalgamation of Pangea. The objective of this article is to decipher which continental margin subducted in the Devonian high‐P–high‐T (HP–HT) event. To do so, a provenance study is presented using combined U–Pb (n = 613) and Lu–Hf (n = 463) isotopic LA–ICP–MS zircon analyses and Sm–Nd whole–rock (n = 5) determinations. These analyses have been performed on five samples of the Banded Gneisses (Cabo Ortegal Complex, NW Iberia), which forms a part of the HP–HT bottom member of the Upper Allochthon. Palaeozoic–Neoproterozoic zircon ages (34.7%) have a maximum abundance at 522–512 Ma, peaks at 575, 561, 545 Ma and minor abundance peaks between 780 and 590 Ma, and show from their Lu–Hf compositions a volcanic arc mixing pattern. This arc was probably related to the Cadomian arc system. The Mesoproterozoic population is scarce and scattered (2.8%), and due to its Lu–Hf pattern, it is proposed that this population is also West Africa Craton derived. The Paleoproterozoic population (39.6%) is concentrated at 2.07 Ga and it is linked to the Eburnean Orogeny, where depleted mantle derived magmas intruded an Archean craton margin. This craton is represented by the Archean population (22.8%), which is grouped at 3.0, 2.68‐2.61 and 2.52‐2.48 Ga, and shows long‐term reworking processes and at least two juvenile magma intrusions. These data show that the Variscan Upper Allochthon has a West African provenance and therefore, it strongly suggests that the NW Iberian allochthonous complexes and their correlative European terranes are also West Africa derived. These results allow us to finally clarify that the first high‐P event, recorded during the eo‐Variscan amalgamation of Pangea, was attained by the subduction of the margin of Gondwana under Laurussia.     

Migmatites record multiple episodes of crustal anatexis and geochemical differentiation in the Sulu ultrahigh‐pressure metamorphic zone,eastern China

Partial melting of ultrahigh‐pressure (UHP) metamorphic rocks is common during collisional orogenesis and post‐collisional reworking, indicating that determining the timing and processes involved in this partial melting can provide insights into the tectonic evolution of collisional orogens. This study presents the results of a combined whole‐rock geochemical and zirconological study of migmatites from the Sulu orogen in eastern China. These data provide evidence of multiple episodes of crustal anatexis and geochemical differentiation within the UHP metamorphic rocks. The leucosomes contain higher concentrations of Ba and K and lower concentrations of the rare earth elements (REE), Th and Y, than associated melanosomes and granitic gneisses. The leucosomes also have homogenous Sr–Nd–O isotopic compositions that are similar to proximal (i.e. within the same outcrop) melanosomes, suggesting that the anatectic melts were generated by the partial melting of source rocks that are located within individual outcrops. The migmatites contain zircons with six different types of domains that can be categorized using differences in structures, trace element compositions, and U–Pb ages. Group I domains are relict magmatic zircons that yield middle Neoproterozoic U–Pb ages and contain high REE concentrations. Group II domains represent newly grown metamorphic zircons that formed at 230 ± 1 Ma during the collisional orogenesis. Groups III, IV, V, and VI zircons are newly grown anatectic zircons that formed at 222 ± 2 Ma, 215 ± 1 Ma, 177 ± 2 Ma, and 152 ± 2 Ma, respectively. The metamorphic zircons have higher Th/U and lower (Yb/Gd)_N values, flat heavy REE (HREE) patterns with no significantly negative Eu anomalies relative to the anatectic zircons, which are characterized by low Th/U ratios, steep HREE patterns, and negative Eu anomalies. The first two episodes of crustal anatexis occurred during the Late Triassic at c. 222 Ma and c. 215 Ma as a result of phengite breakdown. The other two episodes of anatexis occurred during the Jurassic period at c. 177 Ma and c. 152 Ma and were associated with extensional collapse of the collision‐thickened orogen. The majority of Triassic anatectic zircons and all of the Jurassic zircons are located within the leucosomes, whereas the melanosomes are dominated by Triassic metamorphic zircons, suggesting that the leucosomes within the migmatites record more episodes of crustal anatexis. Both metamorphic and anatectic zircons have elevated ε_Hf(t) values compared with relict magmatic zircon cores, suggesting that these zircons contain non‐zircon Hf derived from material with more radiogenic Hf isotope compositions. Therefore, the Sulu and Dabie orogens experienced different episodes of reworking during the exhumation and post‐collisional stages.     

Zircon (re)crystallization during short‐lived,high‐P granulite facies metamorphism (Eger Complex,NW Bohemian Massif)

The Eger Complex in the northwestern Bohemian Massif consists mainly of amphibolite facies granitic gneisses containing a subordinate volume of felsic granulites. Microstructural changes and modelling of metamorphic conditions for both rock types suggest a short‐lived static heating from ~760 to ~850 °C at a constant pressure of ~16 kbar, which led to the partial granulitization of the granitoid rocks. Detailed study of the protolith zircon modifications and modelling of the Zr re‐distribution during the transition from amphibolite to granulite facies suggests that the development of c. 340 Ma old zircon rims in the granulite facies sample is the result of recrystallization of older (c. 475 Ma) protolith zircon. This study suggests that the partial granulitization is a result of a short exposure of the Eger Complex metagranitoids to a temperature of ~850 °C at the base of an arc/fore‐arc domain and their subsequent rapid exhumation during the Lower Carboniferous collision along the western margin of the Bohemian Massif.     

Extruding the Borborema Province (NE‐Brazil): a two‐stage Neoproterozoic collision process

We propose that the development of the Borborema Province from 620 to 570 Ma resulted from two discrete collisional events. Collision I, along the West Gondwana Orogen on the west side of the Province, took place at c. 620–610 Ma as the result of collision between the Parnaíba Block, as the forefront of the much larger Amazonian‐West Africa Craton, and the old basement of the Borborema Province. The suture zone related to this collision was reactivated by a dextral transform zone (the Transbrasiliano Lineament), allowing the Borborema Province to approach and collide against the São Francisco Craton in the south at c. 590–580 Ma marking collision II along the Sergipano Orogen. The combined stresses related to eastward push from collision I and northward push from the cratonic indentation into a thickened lithosphere gave rise to an extensive network of strike‐slip shear zones across the Province, forcing its northeastward extrusion.     

Reply to the Discussion by Rose of Phelps et al. (2014) ‘Oceanographic and eustatic control of carbonate platform evolution and sequence stratigraphy on the Cretaceous (Valanginian–Campanian) passive margin,northern Gulf of Mexico’, Sedimentology, 61, 461–496

Discussion points raised by Rose ( 2016 ) concentrate on late Albian stratigraphic relationships between formations of the East Texas Basin and the San Marcos Arch of the Comanche Platform in the northern Gulf of Mexico. Criticisms of Phelps et al. (2014) regarding stratigraphic nomenclature, palaeogeography and regional lithostratigraphic correlations generally focus on interpretive aspects of the study or do not account for the full scope of published information. Revisions to the top Aptian–Albian Supersequence boundary by Rose are incompatible with the relative location of a subaerial unconformity, as well as deepening lithofacies trends and retrogradational stratigraphic patterns below the interpreted boundary. Rose's placement of the top Aptian–Albian Supersequence boundary precisely at the Albian–Cenomanian stage boundary also implies ca 1·4 Ma of diachroneity in second order sea‐level patterns between the northern Gulf of Mexico and other documented global sedimentary basins.     

Phase equilibria modelling and LASS monazite petrochronology: P–T–t constraints on the evolution of the Priest River core complex,northern Idaho

Phase equilibria modelling, laser‐ablation split‐stream (LASS)‐ICP‐MS petrochronology and garnet trace‐element geochemistry are integrated to constrain the P–T–t history of the footwall of the Priest River metamorphic core complex, northern Idaho. Metapelitic, migmatitic gneisses of the Hauser Lake Gneiss contain the peak assemblage garnet + sillimanite + biotite ± muscovite + plagioclase + K‐feldspar ± rutile ± ilmenite + quartz. Interpreted P–T paths predict maximum pressures and peak metamorphic temperatures of ~9.6–10.3 kbar and ~785–790 °C. Monazite and xenotime ²⁰⁸Pb/²³²Th dates from porphyroblast inclusions indicate that metamorphism occurred at c. 74–54 Ma. Dates from HREE‐depleted monazite formed during prograde growth constrain peak metamorphism at c. 64 Ma near the centre of the complex, while dates from HREE‐enriched monazite constrain the timing of garnet breakdown during near‐isothermal decompression at c. 60–57 Ma. Near‐isothermal decompression to ~5.0–4.4 kbar was followed by cooling and further decompression. The youngest, HREE‐enriched monazite records leucosome crystallization at mid‐crustal levels c. 54–44 Ma. The northernmost sample records regional metamorphism during the emplacement of the Selkirk igneous complex (c. 94–81 Ma), Cretaceous–Tertiary metamorphism and limited Eocene exhumation. Similarities between the Priest River complex and other complexes of the northern North American Cordillera suggest shared regional metamorphic and exhumation histories; however, in contrast to complexes to the north, the Priest River contains less partial melt and no evidence for diapiric exhumation. Improved constraints on metamorphism, deformation, anatexis and exhumation provide greater insight into the initiation and evolution of metamorphic core complexes in the northern Cordillera, and in similar tectonic settings elsewhere.     

Thermal evolution of the Tseel terrane,SW Mongolia and its relation to granitoid intrusions in the Central Asian Orogenic Belt

The timing and thermal effects of granitoid intrusions into accreted sedimentary rocks are important for understanding the growth process of continental crust. In this study, the petrology and geochronology of pelitic gneisses in the Tseel area of the Tseel terrane, SW Mongolia, are examined to understand the relationship between igneous activity and metamorphism during crustal evolution in the Central Asian Orogenic Belt (CAOB). Four mineral zones are recognized on the basis of progressive changes in the mineral assemblages in the pelitic gneisses, namely: the garnet, staurolite, sillimanite and cordierite zones. The gneisses with high metamorphic grades (i.e. sillimanite and cordierite zones) occur in the central part of the Tseel area, where granitoids are abundant. To the north and south of these granitoids, the metamorphic grade shows a gradual decrease. The composition of garnet in the pelitic gneisses varies systematically across the mineral zones, from grossular‐rich garnet in the garnet zone to zoned garnet with grossular‐rich cores and pyrope‐rich rims in the staurolite zone, and pyrope‐rich garnet in the sillimanite and cordierite zones. Thermobarometric analyses of individual garnet crystals reveal two main stages of metamorphism: (i) a high‐P and low‐T stage (as recorded by garnet in the garnet zone and garnet cores in the staurolite zone) at 520–580 °C and 4.5–7 kbar in the kyanite stability field and (ii) a low‐P and high‐T stage (garnet rims in the staurolite zone and garnet in the sillimanite and cordierite zones) at 570–680 °C and 3.0–6.0 kbar in the sillimanite stability field. The earlier high‐P metamorphism resulted in the growth of kyanite in quartz veins within the staurolite and sillimanite zones. The U–Pb zircon ages of pelitic gneisses and granitoids reveal that (i) the protolith (igneous) age of the pelitic gneisses is c. 510 Ma; (ii) the low‐P and high‐T metamorphism occurred at 377 ± 30 Ma; and (iii) this metamorphic stage was coeval with granitoid intrusion at 385 ± 7 Ma. The age of the earlier low‐T and high‐P metamorphism is not clearly recorded in the zircon, but probably corresponds to small age peaks at 450–400 Ma. The low‐P and high‐T metamorphism continued for c. 100 Ma, which is longer than the active period of a single granitoid body. These findings indicate that an elevation of geotherm and a transition from high‐P and low‐T to low‐P and high‐T metamorphism occurred, associated with continuous emplacement of several granitoids, during the crustal evolution in the Devonian CAOB.