Folding in power-law viscous multi-layers

We study high-amplitude folding in layered rocks with two-dimensional numerical simulations. We employ the finite-element method to model shortening of an incompressible multi-layer with power-law viscous rheology. The Lagrangian numerical mesh is deformed and re-meshed to accurately follow the layer interfaces. Three settings are considered: (i) pure shearing of a confined multi-layer, (ii) simple shearing of a multi-layer above a detachment, and (iii) slump folding owing to gravity sliding. In our pure shear simulations, finite-amplitude folds always develop despite confinement and thin weak interlayers. The fold shapes can be significantly irregular, resulting from initial geometrical heterogeneities that are perturbations of the layer interfaces and differences in layer thickness. The bulk normal viscosity of the multi-layer decreases significantly with progressive folding. This structural softening decreases the bulk normal viscosities by a factor of 2-20. For simple shear, the multi-layer does not develop asymmetric fold shapes significantly. Fold axial planes in the multi-layer are mostly curved and not parallel. For slump folding, fold shapes can be significantly asymmetric exhibiting strongly curved fold axial planes and overturned fold limbs. The rheology of the competent layers has a major impact on the fold shapes for gravity-driven multi-layer folding.