Maria-Gema Llorens

Converging flow and anisotropy cause large-scale folding in Greenland's ice sheet

The increasing catalogue of high-quality ice-penetrating radar data provides a unique insight in the internal layering architecture of the Greenland ice sheet. The stratigraphy, an indicator of past deformation, highlights irregularities in ice flow and reveals large perturbations without obvious links to bedrock shape. In this work, to establish a new conceptual model for the formation process, we analysed the radar data at the onset of the Petermann Glacier, North Greenland, and created a three-dimensional model of several distinct stratigraphic layers. We demonstrate that the dominant structures are cylindrical folds sub-parallel to the ice flow. By numerical modelling, we show that these folds can be formed by lateral compression of mechanically anisotropic ice, while a general viscosity contrast between layers would not lead to folding for the same boundary conditions. We conclude that the folds primarily form by converging flow as the mechanically anisotropic ice is channelled towards the glacier.

When do folds unfold during progressive shear

Folded layers in rocks can be stretched again, potentially unfolding the folds back to straight layers. Little is known, however, about how to recognize partly or even entirely unfolded layers. When folded layers can unfold, what determines their mechanical behavior, and how can we recognize them in the field? In order to address these questions, we present a series of numerical simulations of the stretching of previously folded single layers and multilayers in simple shear. Layers do not completely unfold when they undergo softening before or during the stretching process, or when adjacent competent layers prevent them from unfolding. Intrafolial folds and cusp-like folds adjacent to straight layers as well as variations in fold amplitudes and limb lengths of irregular folds are indicative of stretching of a fold train.

Subgrain Rotation Recrystallization During Shearing: Insights From Full-Field Numerical Simulations of Halite Polycrystals

We present, for the first time, results of full-field numerical simulations of subgrain rotation recrystallization of halite polycrystals during simple shear deformation. The series of simulations show how microstructures are controlled by the competition between (i) grain size reduction by creep by dislocation glide and (ii) intracrystalline recovery encompassing subgrain coarsening (SGC) by coalescence through rotation and alignment of the lattices of neighboring subgrains. A strong grain size reduction develops in models without intracrystalline recovery, as a result of the formation of high-angle grain boundaries when local misorientations exceed 15°. The activation of subgrain coarsening associated with recovery decreases the stored strain energy and results in grains with low intracrystalline heterogeneities. However, this type of recrystallization does not significantly modify crystal preferred orientations. Lattice orientation and grain boundary maps reveal that this full-field modelling approach is able to successfully reproduce the evolution of dry halite microstructures from laboratory deformation experiments, thus opening new opportunities in this field of research. We demonstrate how the mean subgrain boundary misorientations can be used to estimate the strain accommodated by dislocation glide using a universal scaling exponent of about 2/3, as predicted by theoretical models. In addition, this strain gauge can be potentially applied to estimate the intensity of intracrystalline recovery, associated with temperature, using quantitative crystallographic analyses in areas with strain gradients.

Mechanics of Fold Development in Pure- and Simple Shear

Folds not only help to determine the orientation and magnitude of shortening, but contain much more information on deformation mechanics, kinematics and rheology. The use of folds for strain analysis requires a first-order quantification of the relationships between parameters that determine fold geometry. We present series of 2-dimensional simulations of single-layer folding under pure and simple shear using the software packages ELLE and BASIL. The orientation of the competent layer, the competence contrast between layer and matrix, the stress exponent of the power-law viscous material as well as the boundary conditions are systematically varied to observe the influence of these parameters on fold development. Analyses of stress heterogeneities show that stress localizations are not regularly distributed at the initial stages of folding. Folding of a competent layer requires lower energy (work) in simple shear than in pure shear conditions.