No abstract

No abstract

<p><em>The processes that govern rock (trans)formation (deposition, deformation, segregation, metamorphism) can result in the development of layering and rock fabrics. Rocks can thus exhibit extrinsic or intrinsic anisotropy at various spatial scales. Anisotropy has important mechanical consequences, in particular, for strain localisation in the lithosphere. This effect is typically not included in geodynamic models. Mechanical anisotropy can be modelled by explicitly modelled by numerically resolving layers of different strengths. Due to the expensive computational cost, this approach is not suitable for large scale geodynamic models. The latter may rather benefit from an upscaling approach that involves anisotropic constitutive laws.  To model the evolution of such material Mühhlaus, (2002) proposed the use of the director vector which corresponds to a single orientation that is changing throughout the process of deformation. We have implemented visco-elasto-plastic anisotropic constitutive laws and the director vector approach in the geodynamic simulation tool MDoodz7.0. Here we present  the rheological implementation, we show some simple simulations involving anisotropic flow and discuss the potential role of anisotropy for large-scale geodynamic processes.</em></p>

<p>The mechanical characteristics of a shell, having a double curvature, are fundamentally different to the characteristics of a plate, having no curvature in its undeformed state. Geometrically, the Earth’s lithosphere is a shell rather than a plate. However, most geodynamic numerical models applied to study the deformation of the lithosphere do not consider this curvature. It is currently unclear whether the shell-type geometry of the lithosphere has a significant impact on lithosphere deformation on the scale of few 1000 kilometers. This study investigates the importance of considering lithospheric shells and compares numerical results of a shortening shell-type and plate-type lithosphere. We apply the two-dimensional state-of-the-art thermo-mechanical code MDoodz (Duretz et al. 2021). We consider a shortening lithosphere in an initially curved and in an initially rectangular geometry and calculate the spatio-temporal stress distribution inside the deforming lithosphere. We further present preliminary results on the effects and relative importance of various softening mechanism, leading to strain localization and subduction initiation, such as thermal softening, grain size reduction, or anisotropy generation due to fabric development.</p> <p> </p> <p>REFERENCES</p> <p>Duretz T., R. de Borst and P. Yamato (2021), Modeling Lithospheric Deformation Using a Compressible Visco-Elasto-Viscoplastic Rheology and the Effective Viscosity Approach, Geochemistry, Geophysics, Geosystems, Vol. 22 (8), e2021GC009675</p>

<p>Strain localization and associated softening mechanisms in a deforming lithosphere are important for subduction initiation or the generation of tectonic nappes during orogeny. Many strain localization and softening mechanisms have been proposed as being important during the viscous, creeping, deformation of the lithosphere, such as thermal softening, grain size reduction, reaction-induced softening or anisotropy development. However, which localization mechanism is the controlling one and under which deformation conditions is still contentious. In this contribution, we focus on strain localization in viscous material due to the generation of anisotropy, for example due to the development of a foliation. We numerically model the generation and evolution of anisotropy during two-dimensional viscous simple shear in order to quantify the impact of anisotropy development on strain localization and on the effective softening. We calculate the finite strain ellipse during viscous deformation. The aspect ratio of the finite strain ellipse serves as proxy for the magnitude and evolution of anisotropy, which determines the ratio of normal to tangential viscosity. To track the orientation of the anisotropy during deformation we apply a director method. We benchmark our implementation of anisotropy by comparing results of two independently developed numerical algorithms based on the finite difference method: one algorithm employs a direct solver and the other a pseudo-transient iterative solver. We will present results of our numerical simulations and discuss their application to natural shear zones.</p>