The role of viscous heating in Barrovian metamorphism of collisional orogens: thermomechanical models and application to the Lepontine Dome in the Central Alps

Thermal models for Barrovian metamorphism driven by doubling the thickness of the radiogenic crust typically meet difficulty in accounting for the observed peak metamorphic temperature conditions. This difficulty suggests that there is an additional component in the thermal budget of many collisional orogens. Theoretical and geological considerations suggest that viscous heating is a cumulative process that may explain the heat deficit in collision orogens. The results of 2D numerical modelling of continental collision involving subduction of the lithospheric mantle demonstrate that geologically plausible stresses and strain rates may result in orogen‐scale viscous heat production of 0.1 to >1 μW m^?3, which is comparable to or even exceeds bulk radiogenic heat production within the crust. Thermally induced buoyancy is responsible for crustal upwelling in large domes with metamorphic temperatures up to 200 °C higher than regional background temperatures. Heat is mostly generated within the uppermost mantle, because of large stresses in the highly viscous rocks deforming there. This thermal energy may be transferred to the overlying crust either in the form of enhanced heat flow, or through magmatism that brings heat into the crust advectively. The amplitude of orogenic heating varies with time, with both the amplitude and time‐span depending strongly on the coupling between heat production, viscosity and collision strain rate. It is argued that geologically relevant figures are applicable to metamorphic domes such as the Lepontine Dome in the Central Alps. We conclude that deformation‐generated viscous dissipation is an important heat source during collisional orogeny and that high metamorphic temperatures as in Barrovian type metamorphism are inherent to deforming crustal regions.     

Intrusion of ultramafic magmatic bodies into the continental crust: Numerical simulation

Intrusions of ultramafic bodies into the lower density continental crust are documented for a large variety of tectonic settings spanning continental shields, rift systems, collision orogens and magmatic arcs. The intriguing point is that these intrusive bodies have a density higher by 300-500 kg m⁻³ than host rocks. Resolving this paradox requires an understanding of the emplacement mechanism. We have employed finite differences and marker-in-cell techniques to carry out a 2D modeling study of intrusion of partly crystallized ultramafic magma from sublithospheric depth to the crust through a pre-existing magmatic channel. By systematically varying the model parameters we document variations in intrusion dynamics and geometry that range from funnel- and finger-shaped bodies (pipes, dikes) to deep seated balloon-shaped intrusions and flattened shallow magmatic sills. Emplacement of ultramafic bodies in the crust lasts from a few kyr to several hundreds kyr depending mainly on the viscosity of the intruding, partly crystallized magma. The positive buoyancy of the sublithospheric magma compared to the overriding, colder mantle lithosphere drives intrusion while the crustal rheology controls the final location and the shape of the ultramafic body. Relatively cold elasto-plastic crust (T_Moho = 400 °C) promotes a strong upward propagation of magma due to the significant decrease of plastic strength of the crust with decreasing confining pressure. Emplacement in this case is controlled by crustal faulting and subsequent block displacements. Warmer crust (T_Moho = 600 °C) triggers lateral spreading of magma above the Moho, with emplacement being accommodated by coeval viscous deformation of the lower crust and fault tectonics in the upper crust. Strong effects of magma emplacement on surface topography are also documented. Emplacement of high-density, ultramafic magma into low-density rocks is a stable mechanism for a wide range of model parameters that match geological settings in which partially molten mafic-ultramafic rocks are generated below the lithosphere. We expect this process to be particularly active beneath subduction-related magmatic arcs where huge volumes of partially molten rocks produced from hydrous cold plume activity accumulate below the overriding lithosphere.     

Cl/Br ratios and stable chlorine isotope analysis of magmatic–hydrothermal fluid inclusions from Butte,Montana and Bingham Canyon,Utah

A bulk geochemical study has been carried out on fluid inclusion leachates extracted from quartz veins from porphyry Cu deposits in Butte, Montana, USA and Bingham Canyon, Utah, USA. The leachates mostly represent low-salinity magmatic–hydrothermal fluid inclusions. Their halogen ratios (Br/Cl) of fluid inclusion leachates were determined by ion chromatography, and δ³⁷Cl values of the leachates were measured by continuous-flow isotope ratio mass spectrometry. Br/Cl ratios from early pre-Main stage and later Main stage veins at Butte range from 0.60 to 1.88 × 10⁻³ M. Ratios are similar in pre-Main stage veins with sericite bearing selvages and Main stage samples ranging from 0.81 to 1.08 × 10⁻³ and from 0.92 to 1.88 × 10⁻³ M, respectively, clustering below seawater (1.54 × 10⁻³ M) and overlapping mantle values (~1–2 × 10⁻³ M). Two samples associated with early pre-Main stage potassic alteration yield distinctly lower Br/Cl ratios of 0.60 and 0.64 × 10⁻³ M. Butte δ³⁷Cl values range from −0.8‰ to −2.3‰ with no significant difference between pre-Main stage and Main stage samples. Br/Cl ratios for quartz veins from Bingham Canyon range from 0.18 to 3.68 × 10⁻³ M. Br/Cl ratios from Bingham range above and below previously reported for porphyry copper deposits. In contrast to Butte, δ³⁷Cl values for Bingham are lower, ranging from −0.9‰ to −4.1‰. In the absence of any processes which can significantly fractionate chlorine isotopes at high temperatures, we suggest that the porphyry system at Bingham, and to a lesser extent at Butte, have inherited negative chlorine isotopic signatures from the subducting slab generated at low temperatures.     

Exhumation mechanisms of melt‐bearing ultrahigh pressure crustal rocks during collision of spontaneously moving plates

A series of 2D petrological–thermomechanical numerical experiments was conducted to: (i) characterize the variability of exhumation mechanisms of ultrahigh pressure metamorphic (UHPM) rocks during collision of spontaneously moving plates and (ii) study the possible geodynamic effects of melting at ultrahigh pressure conditions for the exhumation of high‐temperature–ultrahigh pressure metamorphic (HT–UHPM) rocks. To this end, the models include fluid‐ and melt‐induced weakening of rocks. Five distinct modes of exhumation of (U)HPM rocks associated with changes in several parameters in the models of plate collision and continent subduction are identified as follows: vertical crustal extrusion, large‐scale crustal stacking, shallow crustal delamination, trans‐lithospheric diapirism, and channel flow. The variation in exhumation mechanisms for (U)HPM rocks in numerical models of collision driven by spontaneously moving plates contrasts with the domination of the channel flow mode of exhumation in a majority of the published results from numerical models of collision that used a prescribed plate convergence velocity and/or did not include fluid‐ and melt‐induced weakening of rocks. This difference in the range of exhumation mechanisms suggests that the prescribed convergence velocity condition and the neglect of fluid‐ and melt‐related weakening effects in the earlier models may inhibit development of several important collisional processes found in our experiments, such as slab breakoff, vertical crustal extrusion, large‐scale stacking, shallow crustal delamination and relamination, and eduction of the continental plate. Consequently, the significance of channel flow for the exhumation of UHPM rocks may have been overstated based on the results of the earlier numerical experiments. In addition, the results from this study extend over a larger proportion of the high‐temperature range of P–T conditions documented from UHPM rocks, including those retrieved from HT–UHPM rocks, than the results of experiments from previous numerical models. In particular, the highest peak metamorphic temperatures (up to 1000 °C) are recorded in the case of the vertical crustal extrusion model in which subducted continental crust is subjected to a period of prolonged heating by asthenospheric mantle abutting the continental side of the vertically hanging slab. Nonetheless, some extreme temperature conditions which have been suggested for the Kokchetav and Bohemian massifs, perhaps up to 1100–1200 °C, are still to be achieved in experiments using numerical models.     

The Fe–C–O–H–N system at 6.3–7.8 GPa and 1200–1400 °C: implications for deep carbon and nitrogen cycles

Interactions in a Fe–C–O–H–N system that controls the mobility of siderophile nitrogen and carbon in the Fe⁰-saturated upper mantle are investigated in experiments at 6.3–7.8 GPa and 1200–1400 °C. The results show that the γ-Fe and metal melt phases equilibrated with the fluid in a system unsaturated with carbon and nitrogen are stable at 1300 °C. The interactions of Fe₃C with an N-rich fluid in a graphite-saturated system produce the ε-Fe₃N phase (space group P6₃/mmc or P6₃22) at subsolidus conditions of 1200–1300 °C, while N-rich melts form at 1400 °C. At IW- and MMO-buffered hydrogen fugacity (fH₂), fluids vary from NH₃- to H₂O-rich compositions (NH₃/N₂?>?1 in all cases) with relatively high contents of alkanes. The fluid derived from N-poor samples contains less H₂O and more carbon which mainly reside in oxygenated hydrocarbons, i.e., alcohols and esters at MMO-buffered fH₂ and carboxylic acids at unbuffered fH₂ conditions. In unbuffered conditions, N₂ is the principal nitrogen host (NH₃/N₂?≤?0.1) in the fluid equilibrated with the metal phase. Relatively C- and N-rich fluids in equilibrium with the metal phase (γ-Fe, melt, or Fe₃N) are stable at the upper mantle pressures and temperatures. According to our estimates, the metal/fluid partition coefficient of nitrogen is higher than that of carbon. Thus, nitrogen has a greater affinity for iron than carbon. The general inference is that reduced fluids can successfully transport volatiles from the metal-saturated mantle to metal-free shallow mantle domains. However, nitrogen has a higher affinity for iron and selectively accumulates in the metal phase, while highly mobile carbon resides in the fluid phase. This may be a controlling mechanism of the deep carbon and nitrogen cycles.     

Batagayite,CaZn2(Zn,Cu)6(PO4)4(PO3OH)3·12H2O,a new phosphate mineral from Këster tin deposit (Yakutia,Russia): occurrence and crystal structure

Mineralogy and Petrology - Batagayite, CaZn2(Zn,Cu)6(PO4)4(PO3OH)3·12H2O, is a new secondary phosphate mineral from the Këster deposit, Arga-Ynnykh-Khai massif, NE Yakutia, Russia. It is...     

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.     

Metapelites of the Kanskiy Granulite Complex (Eastern Siberia): Kinked P-T Paths and Geodynamic Model

The Southern Yenisey Range (Eastern Siberia) consists of thegranulite-facies Kanskiy complex bordered by the lower-gradeYeniseyskiy and Yukseevskiy complexes. Samples of metapeliteof the Kanskiy complex typically show characteristic garnet-formingreaction textures and near-isobaric cooling P–T paths.An important new result of this study concerns the differencein shape of the P–T paths from different parts of theKanskiy granulite complex: metapelites collected 8 km from theboundary with the Yeniseyskiy complex followed a linear pathwith dP/dT 0·006 kbar/°C; metapelites collected3 km from this boundary reveal a kinked P–T path withan interval of burial cooling (dP/dT –0·006 kbar/°C).The difference in the shape of the P–T paths is supportedby the chemical zoning of garnet studied in the second groupof samples. A mechanism of buoyant exhumation of granulite issuggested by comparison with the results of numerical modelling,which indicate that such a diversity of P–T paths mayresult from a transient disturbance of the thermal structureby rapid differential movement of material from different crustallevels. To arrive at a correct tectonic interpretation, thewhole assemblage of interrelated P–T paths of metamorphicrocks collected from different localities within the same complexmust be studied. KEY WORDS: crustal diapirism; exhumation; granulites; numerical modelling; P–T path     

Inherent gravitational instability of thickened continental crust with regionally developed low- to medium-pressure granulite facies metamorphism

Petrological arguments show that regionally developed low- to medium-pressure, high-temperature granulite facies metamorphism may critically enhance the lowering of crustal density with depth. This leads to gravitational instability of homogeneously thickened continental crust, mainly due to changes in mineral assemblages and the thermal expansion of minerals in conjunction with the exponential lowering of the effective viscosity of rocks with increasing temperature. It is argued that crustal processes of gravitational redistribution (crustal diapirism) contributing to the exhumation of granulite facies rocks may be activated in this way.