Contrasting tectono-metamorphic evolution of orogenic lower crust in the Bohemian Massif: A numerical model.

International audienc

Porous Melt Flow in Continental Crust—A Numerical Modeling Study

International audienceAbstract In continental crust, rapid melt flow through macroscopic conduits is usually envisaged as the most efficient form of melt transport. In contrast, there is growing evidence that in hot continental crust, grain‐scale to meso‐scale porous melt flow may operate over long distances and over millions of years. Here, we investigate the dynamics of such porous melt flow by means of two‐dimensional thermo‐mechanical numerical models using the code ASPECT. Our models are crustal‐scale and describe the network of pores through which the melt flows by permeability that depends on the spacing of the pores. Our results suggest that assuming realistic material properties, melt can slowly migrate in the hot and thick continental crust through pores with a characteristic spacing of 1 mm or larger. Despite its low velocity (millimeters to centimeters per year), over millions of years, such flow can create large partially molten zones in the middle‐lower crust and significantly affect its thermal state, deformation, and composition. We examined the role of the permeability, melt and solid viscosities, the slope of the melting curve and temperature conditions. We obtained contrasting styles of melt distribution, melt flow, and solid deformation, which can be categorized as melt‐enhanced convection, growth of partially molten diapirs and melt percolation in porosity waves. Our numerical experiments further indicate that grain‐scale porous flow is more likely in rocks where the melt productivity increases slowly with temperature, such as in metaigneous rocks

Re-evaluation of polyphase kinematic and 40Ar/39Ar cooling history of Moldanubian hot nappe at the eastern margin of the Bohemian Massif

International audienceA structural and geochronological 40Ar/39Ar study was performed in kilometre-scale middle and lower crustal lens-shaped domains dominated by a preserved subvertical foliation, surrounded by horizontally foliated migmatites. These domains occur within the Moldanubian nappe overlying the Brunia microcontinent at the eastern margin of the European Variscides. Three main deformation phases were recognized: subvertical S2 fabric trending NW–SE in lower crustal rocks and NE–SW in mid-crustal rocks. It is reworked by HT/MT horizontal fabric S3 along margins of crustal domains and in surrounding migmatites. S3 bears a prolate NE lineation parallel to the S2–S3 intersection in the lower crustal domain. In the middle crustal units, L3 is weak, connected to oblate strain and trends NE–SW parallel to the S2–S3 intersection. D4 non-coaxial shear deformation is mainly localized at the boundary between the Moldanubian nappe and Brunia and bears strong top to the NNE shear criteria. In order to constrain kinematics of the D3 deformation, strain modelling was performed to show that the Moldanubian hot nappe was frontally thrust over the Brunia indentor. The renewed D4 tangential movement only heterogeneously reactivates the horizontal S3. This evolution is recorded in 40Ar/39Ar amphibole cooling ages, which show two statistically significant Carboniferous peaks at ~342 and 332 Ma, which are also reflected by published detrital muscovite 40Ar/39Ar ages in the adjacent foreland basin. This geochronological record is correlated with progressive erosion of the topographically elevated upper crustal part of the Moldanubian nappe during D3 frontal thrusting, followed by greenschist facies D4 transpressive reactivation and subsequent erosion of high-grade parts of the nappe

European Variscan orogenic evolution as an analogue of Tibetan-Himalayan orogen: Insights from petrology and numerical modeling

International audienc

A community benchmark for viscoplastic thermal convection in a 2-D square box

Numerical simulations of thermal convection in the Earth’s mantle often employ a pseudoplastic rheology in order to mimic the plate-like behavior of the lithosphere. Yet the benchmark tests available in the literature are largely based on simple linear rheologies in which the viscosity is either assumed to be constant or weakly dependent on temperature. Here we present a suite of simple tests based on nonlinear rheologies featuring temperature, pressure, and strain rate-dependent viscosity. Eleven different codes based on the finite volume, finite element, or spectral methods have been used to run five benchmark cases leading to stagnant lid, mobile lid, and periodic convection in a 2-D square box. For two of these cases, we also show resolution tests from all contributing codes. In addition, we present a bifurcation analysis, describing the transition from a mobile lid regime to a periodic regime, and from a periodic regime to a stagnant lid regime, as a function of the yield stress. At a resolution of around 100 cells or elements in both vertical and horizontal directions, all codes reproduce the required diagnostic quantities with a discrepancy of at most $3% in the presence of both linear and nonlinear rheologies. Furthermore, they consistently predict the critical value of the yield stress at which the transition between different regimes occurs. As the most recent mantle convection codes can handle a number of different geometries within a single solution framework, this benchmark will also prove useful when validating viscoplastic thermal convection simula- tions in such geometries