Numerical modelling of subglacial ribs, drumlins, herringbones, and mega‐scale glacial lineations reveals their developmental trajectories and transitions

• Published in: Earth Surface Processes and Landforms
• Main Authors: Ely, Jeremy C., Stevens, David, Clark, Chris D., Butcher, Frances E. G.
• Other Authors: H2020 European Research Council, Natural Environment Research Council
• Format: Article in Journal/Newspaper
• Language: English
• Published: Wiley 2023
• Subjects: Ice Sheet
• Online Access: http://dx.doi.org/10.1002/esp.5529; https://onlinelibrary.wiley.com/doi/pdf/10.1002/esp.5529; https://onlinelibrary.wiley.com/doi/full-xml/10.1002/esp.5529

Abstract Initially a matter of intellectual curiosity, but now important for understanding ice‐sheet dynamics, the formation of subglacial bedforms has been a subject of scientific enquiry for over a century. Here, we use a numerical model of the coupled flow of ice, water, and subglacial sediment to explore the formation of subglacial ribs (i.e., ribbed moraine), drumlins and mega‐scale glacial lineations (MSGLs). The model produces instabilities at the ice–bed interface, which result in landforms resembling subglacial ribs and drumlins. We find that a behavioural trajectory is present. Initially subglacial ribs form, which can either develop into fields of organized drumlins, or herringbone‐type structures misaligned with ice flow. We present potential examples of these misaligned bedforms in deglaciated landscapes, the presence of which means caution should be taken when interpreting cross‐cutting bedforms to reconstruct ice flow directions. Under unvarying ice flow parameters, MSGLs failed to appear in our experiments. However, drumlin fields can elongate into MSGLs in our model if low ice–bed coupling conditions are imposed. The conditions under which drumlins elongate into MSGLs are analogous to those found beneath contemporary ice streams, providing the first mechanism, rather than just an association, for linking MSGLs with ice stream flow. We conclude that the instability theory, as realized in this numerical model, is sufficient to explain the fundamental mechanics and process‐interactions that lead to the initiation of subglacial bedforms, the development of the distinctive types of bedform patterns, and their evolutionary trajectories. We therefore suggest that the first part of the longstanding ‘drumlin problem’ – how and why they come into existence – is now solved. However, much remains to be discovered regarding the exact sedimentary and hydrological processes involved.