From orthogneiss to migmatite: Geochemical assessment of the melt infiltration model in the Gföhl Unit (Moldanubian Zone, Bohemian Massif)

The Gföhl Unit is the largest migmatite terrain of the Variscan orogenic root domain in Europe. Its genesis has been until now attributed to variable degrees of in situ partial melting. In the Rokytná Complex (Gföhl Unit, Czech Republic) there is a well-preserved sequence documenting the entire migmatitization process on both outcrop and regional scales. The sequence starts with (i) banded orthogneiss with distinctly separated monomineralic layers, continuing through (ii) migmatitic mylonitic gneiss, (iii) schlieren migmatite characterised by disappearance of monomineralic layering and finally to (iv) felsic nebulitic migmatite with no relics of the original banding.

While each type of migmatite shows a distinct whole-rock geochemical and Sr–Nd isotopic fingerprint, the whole sequence evolves along regular, more or less smooth trends for most of the elements. Possible mechanisms which could account for such a variation are that the individual migmatite types (i) are genetically unrelated, (ii) originated by equilibrium melting of a single protolith, (iii) formed by disequilibrium melting (with or without a small-scale melt movement) or (iv) were generated by melt infiltration from external source. The first scenario is not in agreement with the field observations and chemistry of the orthogneisses/migmatites. Neither of the remaining hypotheses can be ruled out convincingly solely on whole-rock geochemical grounds. However in light of previously obtained structural, petrologic and microstructural data, this sequence can be interpreted as a result of a process in which the banded orthogneiss was pervasively, along grain boundaries, penetrated by felsic melt derived from an external source.

In terms of this melt infiltration model the individual migmatites can be explained by different degrees of equilibration between the bulk rock and the passing melt. The melt infiltration can be modelled as an open-system process, characterised by changes of the total mass/volume and accompanied by gains/losses in many of the major- and trace elements. The modelling of the mass balance resulted in identification of a component added by a heterogeneous nucleation of feldspars, quartz and apatite from the passing melt. This is in line with the observed presence of new albitic plagioclase, K-feldspar and quartz coatings as well as resorption of relict feldspars. At the most advanced stages (schlieren and nebulitic migmatites) the whole-rock trace-element geochemical variations document an increasing role for fractional crystallization of the K-feldspar and minor plagioclase, with accessory amounts of monazite, zircon and apatite.

The penetrating melt was probably (leuco-) granitic, poor in mafic components, Rb rich, with low Sr, Ba, LREE, Zr, U and Th contents. It probably originated by partial melting of micaceous quartzo-feldspathic rocks.

If true and the studied migmatites indeed originated by a progressive melt infiltration into a single protolith resembling the banded orthogneiss, this until now underappreciated process would have profound implications regarding rheology and chemical development of anatectic regions in collisional orogens.     

Thermal and mechanical behaviour of the orogenic middle crust during the syn‐ to late‐orogenic evolution of the Variscan root zone,Bohemian Massif

We combine structural observations, petrological data and ⁴⁰Ar–³⁹Ar ages for a stack of amphibolite facies metasedimentary units that rims high‐P (HP) granulite facies felsic bodies exposed in the southern Bohemian Massif. The partly migmatitic Varied and Monotonous units, and the underlying Kaplice unit, show a continuity of structures that are also observed in the adjacent Blanský les HP granulite body. They all exhibit an earlier NE?SW striking and steeply NW‐dipping foliation (S3), which is transposed into a moderately NW‐dipping foliation (S4). In both the Varied and Monotonous units, the S3 and S4 foliations are characterized by a Sil–Bt–Pl–Kfs–Qtz–Ilm±Grt assemblage, with occurrences of post‐D4 andalusite, cordierite and muscovite. In the Monotonous unit, minute inclusions of garnet, kyanite, sillimanite and biotite are additionally found in plagioclase from a probable leucosome parallel to S3. The Kaplice unit shows rare staurolite and kyanite relicts, a Sil–Ms–Bt–Pl–Qtz±Grt assemblage associated with S3, retrogressed garnet?staurolite aggregates during the development of S4, and post‐D4 andalusite, cordierite and secondary muscovite. Mineral equilibria modelling for representative samples indicates that the Varied unit records conditions higher than ~7 kbar at 725 °C during the transition from S3 to S4, followed by a P?T decrease from ~5.5 kbar/750 °C to ~4.5 kbar/700 °C. The Monotonous unit shows evidence of partial melting in the S3 fabric at P?T above ~8 kbar at 740–830 °C and a subsequent P?T decrease to 4.5–5 kbar/700 °C. The Kaplice unit preserves an initial medium‐P prograde path associated with the development of S3 reaching peak P?T of ~6.5 kbar/640 °C. The subsequent retrograde path records 4.5 kbar/660 °C during the development of S4. ⁴⁰Ar–³⁹Ar geochronology shows that amphibole and biotite ages cluster at c. 340 Ma close to the HP granulite, whereas adjacent metasedimentary rocks preserve c. 340 Ma amphibole ages, but biotite and muscovite ages range between c. 318 and c. 300 Ma. The P?T conditions associated with S3 imply an overturned section of the orogenic middle crust. The shared structural evolution indicates that all mid‐crustal units are involved in the large‐scale folding cored by HP granulites. The retrograde P–T paths associated with S4 are interpreted as a result of a ductile thinning of the orogenic crust at a mid‐crustal level. The ⁴⁰Ar–³⁹Ar ages overlap with U–Pb zircon ages in and around the HP granulite bodies, suggesting a short duration for the ductile thinning event. The post‐ductile thinning late‐orogenic emplacement of the South Bohemian plutonic complex is responsible for a re‐heating of the stacked units, reopening of argon system in mica and a tilting of the S4 foliation to its present‐day orientation.     

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.     

Juxtaposition of Barrovian and migmatite domains in the Chinese Altai: a result of crustal thickening followed by doming of partially molten lower crust

Ordovician metasedimentary rocks are the oldest and most extensive sedimentary sequence in the Chinese Altai. They experienced two major episodes of deformation (D1 and D2) resulting in the formation of juxtaposed Barrovian‐type and migmatite domains. D1 is characterized by a penetrative sub‐horizontal fabric (S1), and D2 is marked by upright folds (F2) with NW–SE‐trending axial planes in shallow crustal levels and by sub‐vertical transposition foliations (S2) in the high‐grade cores of large‐scale F2 antiforms. In the Barrovian‐type domain, successive growth of biotite, garnet and staurolite is observed in the S1 fabric. Kyanite included in garnet and plagioclase in the migmatite domain is interpreted to have formed also in S1. In the biotite and garnet zones, the spaced S2 cleavage is marked by biotite and muscovite, and in the staurolite and kyanite zones, the penetrative S2 fabric is characterized by sillimanite, locally with late cordierite. Phase equilibria modelling indicates that the S1 fabric was associated with an increase in pressure and temperature under Barrovian‐type conditions in both domains. The S2 fabric was related to decompression, in which rocks in the biotite and garnet zones well preserve the peak assemblage, and the higher grade rocks in the staurolite and kyanite zones re‐equilibrated to different degrees under high‐temperature/low‐pressure (HT/LP) conditions. The D1 metamorphic history is attributed to the progressive burial related to Early–Middle Palaeozoic crustal thickening and the metamorphism associated with D2 is interpreted to result from exhumation by vertical extrusion. The extrusion of hot rocks was contemporaneous with the formation of gneiss domes accompanied by the intrusion of juvenile magmas at middle crustal levels during the Middle Palaeozoic. Consequently, there is a genetic link between the Barrovian‐type and migmatite domains related to continuous transition of the Barrovian‐type fabric into the HT/LP one during development of domal structures in the southern Altai orogenic belt. This study has a broad impact on the understanding of the thermo‐mechanical behaviour of accretionary orogenic systems worldwide. The lower crustal flow and doming of hot crust, so far reported only in continental collisional settings, seems to be also an integral mechanism responsible for both horizontal and vertical redistribution of accreted material prior to continental collision.     

The Variscan belt: correlations and plate dynamics

International Journal of Earth Sciences -     

Metamorphic inheritance of Rheic passive margin evolution and its early‐Variscan overprint in the Teplá‐Barrandian Unit,Bohemian Massif

Rift‐related regional metamorphism of passive margins is usually difficult to observe on the surface, mainly due to its strong metamorphic overprint during the subsequent orogenic processes that cause its exposure. However, recognition of such a pre‐orogenic evolution is achievable by careful characterization of the polyphase tectono‐metamorphic record of the orogenic upper plate. A multidisciplinary approach, involving metamorphic petrology, P–T modelling, structural geology and in situ U‐Pb monazite geochronology using laser‐ablation split‐stream inductively coupled plasma mass spectrometry, was applied to unravel the polyphase tectono‐metamorphic record of metapelites at the western margin of the Teplá‐Barrandian domain in the Bohemian Massif. The study resulted in discovery of three tectono‐metamorphic events. The oldest event M1 is LP–HT regional metamorphism with a geothermal gradient between 30 and 50 °C km^?1, peak temperatures up to 650 °C and of Cambro‐Ordovician age (c. 485 Ma). The M1 event was followed by M2‐D2, which is characterized by a Barrovian sequence of minerals from biotite to kyanite and a geothermal gradient of 20–25 °C km^?1. D2‐M2 is associated with a vertical fabric S2 and was dated as Devonian (c. 375 Ma). Finally, the vertical fabric S2 was overprinted by a D3‐M3 event that formed sillimanite to chlorite bearing gently inclined fabric S3 also of Devonian age. The high geothermal gradient of the M1 event can be explained as the result of an extensional, rift‐related tectonic setting. In addition, restoration of the deep architecture and polarity of the extended domain before the Devonian history – together with the supracrustal sedimentary and magmatic record – lead us to propose a model for formation of an Ordovician passive continental margin. The subsequent Devonian evolution is interpreted as horizontal shortening of the passive margin at the beginning of Variscan convergence, followed by detachment‐accommodated exhumation of lower‐crustal rocks. Both Devonian shortening and detachment occurred in the upper plate of a Devonian subduction zone. The tectonic evolution presented in this article modifies previous models of the tectonic history of the western margin of the Teplá‐Barrandian domain, and also put constraints on the evolution of the southern margin of the Rheic ocean from the passive margin formation to the early phases of Variscan orogeny.     

P–T–t–D records of Early Palaeozoic Andean-type shortening of a hot active margin: The Dunhuang block in NW China

High-pressure (HP) granulites form either in the domain of the subducted plate during continental collision or in supra-subduction systems where the thermally softened upper plate is shortened and thickened. Such a discrepancy in tectonic setting can be evaluated by metamorphic pressure–temperature–time-deformation (P–T–t–D) paths. In the current study, P–T–t–D paths of Early Palaeozoic HP granulite facies rocks, in the form of metabasic lenses enclosed in migmatitic metapelite, from the Dunhuang block, NW China, are investigated in order to constrain the nature of the HP rocks and shed light on the geodynamic evolution of a modern hot orogenic system in an active margin setting. The rocks show a polyphase evolution characterized by (1) relics of horizontal or gently dipping fabric (S1) preserved in cores of granulite lenses and in garnet porphyroblasts, (2) a N-S trending sub-vertical fabric (S2) preserved in low-strain domains and (3) upright folds (F3) associated with a ubiquitous steep E-W striking axial planar foliation (S3). Garnet in the granulites preserves relics of a prograde mineral assemblage M1a equilibrated at ~11.5 kbar and ~770–780°C, whereas the matrix granulite assemblage (M1b) from the S1 fabric attained peak pressure at ~13.5 kbar and ~850°C. The granulites were overprinted at ~8–11 kbar and ~850–900°C during crustal melting (M2) followed by partial re-equilibration (M3) at ~8 kbar and ~625°C. A garnet Lu–Hf age of 421.6 ± 1.2 Ma dates metamorphism M1, while a garnet Sm–Nd age of 385.3 ± 4.0 Ma reflects M3 cooling of the granulites. The mineral assemblage, M1, of the host migmatitic metapelite formed at ~9–12.5 kbar and ~760–810°C, partial melting and migmatization (M2) occurred at ~7 kbar and ~760°C and re-equilibration (M3) at ~5–6 kbar and ~675°C. A garnet Lu–Hf age of 409.7 ± 2.3 Ma dates thermal climax (M2) and a garnet Sm–Nd age of 356 ± 11 Ma constrains M3 for the migmatitic metapelites. The timing of this late phase is also bracketed by an emplacement age of syntectonic granite dated at c. 360 Ma. Decoupling of M1 and M2 P–T evolutions between the mafic granulites and migmatitic metapelites indicates their different positions in the crustal column, while the shared pressure–temperature (P–T) evolution M3 suggests formation of a mélange-like association during the late stages of orogeny. The high-pressure event D1-M1 is interpreted as a result of Late Silurian–Early Devonian moderate crustal thickening of a thermally softened and thinned pre-orogenic crust. The high-temperature (HT) re-equilibration D2-M2 is interpreted as a result of Mid-Devonian shortening of the previously thickened crust, possibly due to ‘Andean-type’ underthrusting. The D3-M3 event reflects Late Devonian supra-subduction shortening and continuous erosion of the sub-crustal lithosphere. This tectono-metamorphic sequence of events is explained by polyphased Andean-type deformation of a ‘Cascadia-type’ active margin, which corresponds to a supra-subduction tectonic switching paradigm.     

Coupling of P–T–t–D histories of eclogite and metagreywacke—Insights to late Ordovician–Silurian crustal folding events recorded in the Beishan Orogen (NW China)

The Beishan complex is composed of orthogneiss and metagreywacke that both enclose bodies of eclogite and serves as a unique example for comparative petrological study of all these lithologies. The rocks show the earliest regional steep N-S striking fabric (S2) preserved in low strain domains that are reworked by ubiquitous steep N-NE dipping cleavage (S3). The eclogite shows an almost isotropic fabric defined by an M1 assemblage of Grt–Cpx–Amp–Qz–Rt–Ilm that is locally retrogressed to M2-3 amphibolite facies mineral assemblages, with P–T peak at 20–21 kbar and 750–775°C and retrogression to 2–3kbar and 530–550°C. The typical mineral assemblage of the host metagreywackes is Bt–Ms–Pl–Qz−Chl–Ilm±Grt. Rare Al-rich metagreywacke layers are composed of Grt–Ky–St±Sil−And–Bt–Ms–Pl–Qz±Chl±Rt–Ilm giving a P–T path with peak at 8–8.5kbar and ~670°C correlated with the S2 fabric and retrogression to ~2.5kbar and 525–550°C correlated with the S3 foliation. In two eclogite samples, the garnet-whole rock-clinopyroxene Lu–Hf isochrons give ages of 461.9±1.6 Ma and 462.0±6.2 Ma interpreted as reflecting average age of garnet formation, and Sm–Nd isochrons give ages of 453.6±2.7 Ma and 452.8±3.0 Ma interpreted as dating near-peak metamorphism. In metagreywacke, in-situ U–Pb dating of monazite gives two groups of ages of 445–440 Ma (Mnz cores) and 436–429 Ma (Mnz rims), interpreted as reflecting the metamorphic peak and retrogression. Our results show that eclogite was formed during Ordovician by subduction of a continental crust (D1). Eclogite and metagreywacke underwent partly decoupled P–T–t–D paths until their juxtaposition at mid-crustal levels during a first late Ordovician–early Silurian D2 shortening. Coupling of their P–T–t–D paths occurred during exhumation in the Silurian and a second and orthogonal D3 shortening event. The data from the Beishan Orogen are consistent with a collisional intra-Gondwanan orogen located south of the Central Asian Orogenic Belt.     

Microstructure mechanism map of dynamically recrystallized marble

The relative nucleus density (RND) model of dynamically recrystallized grain size [Sakai, T., Jonas, J.J. 1984. Dynamic recrystallization: mechanical and microctructutal consideration. Acta metallurgica, 32, 198–209] was applied to experimentally and to naturally deformed marbles that have undergone dynamic recrystallization. The model shows that a relationship between initial grain size (D₀) and stable dynamically recrystallized grain size (D_S) for a given value of temperature-corrected strain-rate (Z) controls grain size evolution during dynamic recrystallization. New microstructural mechanism maps (MM-maps) for experimentally and naturally deformed marbles (based on previously published data) were defined in log grain size–log Z space and show two distinct regions of grain reduction and grain coarsening. The boundary between these two regions corresponds to an equation relating dynamically recrystallized grain size and temperature corrected strain rate, as proposed in this work. The new MM-map was used to trace semi-quantitatively microstructural and grain size evolution in naturally deformed marbles that underwent dynamic recrystallization at different thermal conditions. The boundary between grain coarsening and grain reduction does not necessarily coincide with the boundary between rotation and migration recrystallization mechanisms. Assessment of available natural data shows that the boundary condition D₀ = 2D_S between grain-coarsening and grain-reduction introduced by Sakai and Jonas [Sakai, T., Jonas, J.J. 1984. Dynamic recrystallization: mechanical and microctructutal consideration. Acta metallurgica, 32, 198–209] is not required for naturally deformed marble.