Discrete-element modelling of detachment folding

A two‐dimensional, discrete‐element modelling technique is used to investigate the initiation and growth of detachment folds in sedimentary rocks above a weak décollement level. The model depicts the sedimentary rocks as an assemblage of spheres that obey Newton's equations of motion and that interact with elastic forces under the influence of gravity. Faulting or fracturing between neighbouring elements is represented by a transition from repulsive–attractive forces to solely repulsive forces. The sedimentary sequence is mechanically heterogeneous, consisting of intercalated layers of markedly different strengths and thicknesses. The interlayering of weak and strong layers within the sedimentary rocks promotes the localization of flexural flow deformation within the weak layers. Even with simple displacement boundary conditions, and straightforward interlayering of weak and strong layers, the structural geometries that develop are complex, with a combination of box, lift‐off and disharmonic detachment fold styles forming above the décollement. In detail, it is found that the modelled folds grow by both limb rotation and limb lengthening. The combination of these two mechanisms results in uplift patterns above the folds that are difficult, or misleading, to interpret in terms of simple kinematic models. Comparison of modelling results with natural examples and with kinematic models highlights the complexities of structural interpretation in such settings.