Search

Two glacial deposits in the Gesso valley ( Maritime, Alps) have been 10 Be-dated at 20 140±1080 (weighted mean± SD) and 16 590±970 years, respectively, thus constraining the Last Glacial Maximum ( LGM) and Gschnitz stadials in the southwestern part of the Alps. The LGM age is chronologically coherent with MIS 2 and synchronous with most other LGM moraines in the Alps. The Gschnitz stadial also appears to be in agreement with the ages obtained from other Alpine sites and with Heinrich Event I. This suggests that the Alpine glaciers reacted simultaneously and essentially synchronously with the climate change associated with Heinrich Event 1. The Equilibrium Line Altitudes ( ELAs) of the LGM and Gschnitz reconstructed palaeoglaciers are 1850 and 1910 m a.s.l., respectively. The ELA comparison across the Alps indicates that the palaeoclimate of the Maritime Alps during the LGM was rather different from that of other Alpine sectors. However, the similar Gschnitz ELA value between the Gesso valley and other sites across the mountain chain indicates that Alpine glaciers responded with the same intensity to the climate change associated with Heinrich Event I. Overall, these results suggest that the interaction between the atmospheric circulation of air masses and local Alpine orography was more complex than has previously been argued. [ABSTRACT FROM AUTHOR]

Four boulder samples from the Piano del Praiet frontal moraine in the Gesso della Barra Valley (Maritime Alps) have been 10Be dated. The results give a weighted mean age of 11 340±370 (870) yr, constraining the frontal moraine to the Egesen glacial stadial, during the Younger Dryas cold phase. By applying the same 10Be production rate to other Egesen moraines previously dated in the Alps, we obtain similar ages for all of them. This suggests a synchroneity of the Egesen deglaciation in the European Alps at the end of the Younger Dryas. From the palaeoshape of the Egesen glacier, reconstructed by means of geomorphological mapping, an Equilibrium Line Altitude depression (δELA) of −520 to −530 m, with respect to the present-day ELA, and of −260 to −320 m, with respect to the Little Ice Age ELA, has been calculated. Comparison with other Alpine sector δELAs indicates that the Maritime Alps experienced humid climatic conditions during the Younger Dryas. [ABSTRACT FROM AUTHOR]

We map glacial cirques, and analyse spatial variability in their altitude and aspect to derive a long-term, time-integrated, perspective on climate patterns during former periods of mountain glaciation (likely spanning multiple Quaternary glaciations) in Britain and Ireland. The data reveal that, although air temperatures were important, exposure to moisture-bearing air masses was the key factor in regulating sites of former mountain glacier formation, and indicate that during such periods, moisture supply was largely controlled by North Atlantic westerlies, with notable inland precipitation gradients (precipitation decreasing inland), similar to present day. In places, trends in cirque altitude may also reflect regional differences in the extent of cirque deepening, controlled by the dimensions and dynamics of the glaciers that came to occupy them. Specifically, comparatively deep cirques in coastal locations may reflect the former presence of dynamic (fed by moisture from the North Atlantic), but comparatively small, glaciers (largely confined to their cirques). By contrast, decreasing cirque depth further inland, may reflect the former presence of larger and/or less dynamic ice masses, occupying comparatively continental climatic conditions. [ABSTRACT FROM AUTHOR]

Glacier reconstructions are widely used in palaeoclimatic studies and this paper presents a new semi-automated method for generating glacier reconstructions: GlaRe, is a toolbox coded in Python and operating in ArcGIS. This toolbox provides tools to generate the ice thickness from the bed topography along a palaeoglacier flowline applying the standard flow law for ice, and generates the 3D surface of the palaeoglacier using multiple interpolation methods. The toolbox performance has been evaluated using two extant glaciers, an icefield and a cirque/valley glacier from which the subglacial topography is known, using the basic reconstruction routine in GlaRe. Results in terms of ice surface, ice extent and equilibrium line altitude show excellent agreement that confirms the robustness of this procedure in the reconstruction of palaeoglaciers from glacial landforms such as frontal moraines. [ABSTRACT FROM AUTHOR]

Multi-temporal inventories of glacierised regions provide an improved understanding of water resource availability. In this study, we present a Landsat-based multi-temporal inventory of glaciers in four Upper Indus sub-basins and three internal drainage basins in the Ladakh region for the years 1977, 1994, 2009 and 2019. The study records data on 2257 glaciers (of individual size >0.5 km 2) covering an area of ∼7923±106 km 2 which is equivalent to ∼30 % of the total glacier population and ∼89 % of the total glacierised area of the region. Glacier area ranged between 0.5±0.02 and 862±16 km 2 , while glacier length ranged between 0.4±0.02 and 73±0.54 km. Shayok Basin has the largest glacierised area and glacier population, while Tsokar has the least. Results show that the highest concentration of glaciers is found in the higher elevation zones, between 5000 and 6000 m a.s.l., with most of the glaciers facing towards the NW–NE quadrant. The error assessment shows that the uncertainty, based on the buffer-based approach, ranges between 2.6 % and 5.1 % for glacier area, and 1.5 % and 2.6 % for glacier length with a mean uncertainty of 3.2 % and 1.8 %, respectively. This multitemporal inventory is in good agreement with previous studies undertaken in parts of the Ladakh region. The new glacier database for the Ladakh region will be valuable for policy-making bodies, and future glaciological and hydrological studies. The data can be viewed and downloaded from PANGAEA, 10.1594/PANGAEA.940994 (Soheb et al., 2022). [ABSTRACT FROM AUTHOR]

Understanding the evolution of the ice-bed interface is fundamentally important for gaining insight into the dynamics of ice masses and how subglacial landforms are created. However, the formation of the suite of landforms generated at this boundary — subglacial bedforms — is a contentious issue that is yet to be fully resolved. Bedforms formed in aeolian, fluvial, and marine environments either belong to separate morphological populations or are thought to represent a continuum of forms generated by the same governing processes. For subglacial bedforms, a size and shape continuum has been hypothesised, yet it has not been fully tested. Here we analyse the largest data set of subglacial bedform size and shape measurements ever collated (96,900 bedforms). Our results show that flutes form a distinct population of narrow bedforms. However, no clear distinction was found between drumlins and megascale glacial lineations (MSGLs), which form a continuum of subglacial lineations. A continuum of subglacial ribs also exists, with no clear size or shape distinctions indicating separate populations. Furthermore, an underreported class of bedform with no clear orientation to ice flow (quasi-circular bedforms) overlaps with the ribbed and lineation continua and typically occurs in spatial transition zones between the two, potentially merging these three bedform types into a larger continuum. [ABSTRACT FROM AUTHOR]

A toolbox for the automated calculation of glacier equilibrium-line altitudes (ELAs) using the Accumulation Area Ratio, Area-Altitude Balance Ratio, Area-Altitude and Kurowski methods is presented. These are the most commonly-used methods of ELA calculation in palaeo-glacier reconstructions. The toolbox has been coded in Python and runs in ArcGIS requiring only the reconstructed surface of the palaeo-glacier (a DEM) as input. Through fast and automatic calculation this toolbox simplifies the process of ELA determination and can successfully work both for a single glacier and for large datasets of multiple glaciers. [ABSTRACT FROM AUTHOR]

Ice low-switching, which can involve changes in ice flow velocity and direction, is crucial to a full understanding of ice masses and their response to climate change. A topographically controlled ice flow switch near a glacier margin was recently documented at Breiðamerkurjökull, southeast Iceland, where the central flow unit migrated eastward in response to variations in subglacial topography and the influence of Jökulsárlón glacial lagoon. This site provides an opportunity to study the geomorphic response to ice-margin reconfiguration. Investigating contemporary processes can offer valuable insights into analogous landforms created during the deglaciation of palaeo-ice sheets. The landform assemblage and topographic setting of our Icelandic study site is compared to a palaeo-example from Alberta, Canada, which was once covered by the Laurentide ice sheet. Uncrewed aerial vehicle-(UAV) derived data was used to assess the geomorphic response to this switching and related processes across a 1.5 km2 area of the central flow unit which deglaciated between 2010 and 2023. From 2010 to 2017, the landscape featured streamlined subglacial material, a stable subglacial esker system and proglacial lakes (Landsystem A), shifting to a spillway-dominated system between 2018 and 2023 (Landsystem B). Since 2018 this section of Breiðamerkurjökull has been retreating across a reverse slope bed, resulting in the formation of quasi-annual ice-marginal spillways. Meltwater impoundment at the ice margin, formed ice-contact lakes which eventually initiated ice-margin parallel spillways draining proglacial meltwater along the local land-surface gradient, towards Jökulsárlón. As the ice retreats, an ice-contact lake forms again at the new margin and initiates the erosion of the next ice-marginal spillway. The geomorphological signature demonstrates how subglacial topography and ice-flow switching can significantly influence ice and meltwater dynamics. Since the glacier flow-switch, part of the central unit is now lake-terminating with areas of the margin evolving into a stagnant system, as it is now cut off from the accumulation centre. Therefore, Landsystem B could be analogous to regions of ice stream shut down and where ice masses retreated across reverse slope beds. For example, the Pakowki Lake region of Southeastern Alberta displays a similar landform assemblage and is presented as a palaeo-example in this work. Such insights are important for assessing the efficacy of numerical models in reconstructing the finer scale dynamics of past ice sheets during retreat. • Landform mapping at Breiðamerkurjökull (Iceland) from UAV-derived data reveals two distinct landform assemblages. • Landsystem A (2010-2017): flat subglacial bed, subglacial eskers and proglacial lakes. • Landsystem B (2018-2023): reverse slope bed and ice-marginal spillways. • Landsystem B is analogous to SE Alberta, providing support for existing spillway reconstructions in the region. • The geomorphic response highlights the significance of subglacial topography on ice and meltwater dynamics. [ABSTRACT FROM AUTHOR]

Cirques are glacially eroded, bowl-shaped depressions, characterised by steep headwalls and flat or overdeepened floors. Given their association with past glaciers, cirques are sometimes used as proxies for palaeoclimate. However, cirques are shaped over multiple glacial cycles, and their usefulness as palaeoclimate indicators therefore remains open to question. In this paper, we map 3984 glacier-free cirques across the Scandinavian Peninsula and analyse variations in cirque floor altitude (CFA). We explore the relationships between CFAs and cirque aspect, latitude, longitude, and distance to the coast. We test the validity of using CFAs as indicators of palaeoclimate through comparison with the equilibrium-line altitudes (ELAs) of 513 modern cirque glaciers. Results indicate that both CFAs and modern cirque-glacier ELAs decrease with latitude and vary with aspect, being generally lowest on east-facing slopes. However, the clearest and strongest trend in both CFAs and modern cirque glacier ELAs is an increase in elevation with distance from the modern coast (i.e., distance 'inland'). This likely indicates that similar climatic gradients, particularly an inland reduction in precipitation, acted to regulate former sites of glacier initiation (reflected by CFAs) and modern glacier ELAs. This would imply that CFAs are a useful proxy for palaeoclimate. However, we note that both CFAs and modern ELAs reflect the general topography of this region (with increasing elevations moving inland), and the glacial history of the area (indirectly linked to palaeoclimate) may have played a role in regulating where cirques have formed. For these reasons, we suggest that palaeoclimatic interpretations derived from CFAs should be treated with caution. • We analyse cirque floor altitudes as a proxy for palaeo ELAs across the Scandinavian Peninsula. • Results indicate that precipitation and temperature play a key role in regulating past and present ELAs. • The dominant trend in cirque floor altitudes shows increasing values with distance inland from the coast. [ABSTRACT FROM AUTHOR]

We present new chronological data and Equilibrium Line Altitude (ELA) information for palaeoglaciers in the Maritime Alps during the Last Glacial Maximum (LGM) and the early deglaciation. Three relatively small catchments were investigated to test if the response of small (1–10 km2) glaciers to LGM climatic forcing was distinguishable from that of larger glacial systems. Palaeoglacier geometries and ELAs were reconstructed using geomorphological evidence and toolboxes in a geographic information system. Chronological control was provided through 10Be surface exposure dating of erratic boulders located on frontal and lateral moraine ridges. Our data indicate two phases of glacier advance or stabilisation in the Maritime Alps, the first occurring at ca. 25 to 24 ka and the second at ca. 20.5 to 19 ka. This is consistent with ages that have been reported from larger outlet lobes along the south-western Alpine fringe, where the second LGM advance is usually preserved in the form of a recessional moraine. Within the uncertainty of the dating technique, there are no distinguishable differences between the responses of small and large Alpine glaciers to the LGM climate. The reconstructed ELAs for palaeoglaciers in the Maritime Alps indicate a considerable variability, which appears to be linked to differences in received solar radiation, from ca. 1600 m a.s.l. for north-facing cirque glaciers to almost 2000 m a.s.l. for glaciers with a southerly aspect. We argue that such influence of solar radiation needs to be considered when using ELAs for palaeoclimatic inferences at regional scales. The calculated ELAs are up to 450 m higher than those in the northern Apennines or the Corsican mountains, indicating a relatively dry LGM climate on the Italian side of the Maritime Alps. • Glaciers in the Maritime Alps advanced at 25 to 24 ka and at 20.5 to 19 ka. • Synchronous advances of small and large glaciers across the south-western Alps. • Catchment aspect and solar radiation had a crucial influence on palaeoglacier ELAs. • Regional glacier ELA of 1800 m a.s.l. indicates a relatively dry LGM climate. [ABSTRACT FROM AUTHOR]

Glaciogenic sediments are present in many hydrocarbon-producing basins across the globe but their complex nature makes it difficult to characterise the reservoir-quality sedimentary units. Despite this, Ordovician glacial deposits in North Africa, and Carboniferous-Permian glaciogenic sequences in the Middle East, have been proven to host significant, economical, hydrocarbon accumulations. Additionally, discoveries have been made in the shallow (<1000 m below seabed), glacial, Pleistocene sedimentary succession of the North Sea (e.g. Peon and Aviat). This paper provides a predictive exploration framework in the form of a conceptual model of glaciogenic sediment-landform distributions. The model is based on the extensive onshore glacial sedimentary record integrated with available offshore data. It synthesises the published knowledge, drawing heavily on glacial landsystem models, glacial geomorphology and sedimentology of glaciogenic deposits to provide a novel conceptual model allowing for the efficient description and interpretation of glacial sediments and landforms in the subsurface. Subsequently, land-terminating and water-terminating ice sheet depositional systems are described and discussed, with respect to ice advance and retreat cycles. This detailed description focuses on the macro-scale stratigraphic organisation of glacial sediments with relation to the ice margin, aiding the prediction of glaciogenic sediment distributions, and their likely geometry, architecture and connectivity as reservoirs. • A novel conceptual model for the distribution of glaciogenic reservoirs. • First-pass interpretation tool for complex glaciogenic sequences. • Ice margin position controlling reservoir quality. • Identification of landforms and sediments crucial to reservoir predictability. [ABSTRACT FROM AUTHOR]

This study is focused on the geomorphology of surface sedimentary environments contained in the fluvial, tidal and wave-dominated areas of the Niger Delta. GIS techniques applied to high-quality Landsat and SPOT images were used to identify and map landforms and processes, as well as quantitatively characterise their metrics and spatial distribution. A detailed analysis shows that the extent of the Niger Delta is 70,000 km2, contrary to previous publications that reported 75,000 km2. The delta has also been remapped and classified into mega sedimentary environments: the upper deltaic plain is 69% of the total extent, while the lower deltaic plain and delta front are 25% and 6% respectively. Other subunits were distinguished and mapped within the upper deltaic plain such as fluvial channels, point bars, braid bars, oxbow lakes, and other lakes. Tidal channels, beaches and spits were identified in the lower deltaic plain and delta front. The geometries of these sedimentary bodies (landforms) appear to be relatively scale-invariant, thus meaning that from the measurement of some dimensional parameters (e.g. the length) it is possible to estimate all others (e.g. the width and overall extent). This is an aspect of great relevance in subsurface analyses where it is often difficult to predict the extent of sedimentary bodies. The correlation between channel width and length (r2 = 0.9), is the strongest relationship found between sedimentary body metrics. • This work presents the geomorphology of Niger Delta surface sedimentary environments. • GIS is applied to map landforms and quantify their metrics and spatial distribution. • Niger Delta is classified as tide-dominated, wave-influenced and fluvial-effect-Tw f. • Present-day Niger Delta geomorphology is relevant to subsurface reservoir studies. [ABSTRACT FROM AUTHOR]

Successful adaptation to climate change at regional scales can often depend on understanding the nature of geomorphological responses to climate change at those scales. Here we use evidence from landscapes which are known to be environmentally sensitive to show that geomorphological change in response to shifts in climate can be highly nonlinear. Our study sites are two mountain massifs on the western coast of Ireland. Both sites have similar geological and Pleistocene glacial histories and are similar topographically, geomorphologically and in their climate histories. We show that despite these similarities their response to late Holocene, climate change has differed. Both massifs have responded to short-term climate changes over the last 4500 years that are considered to have been uniform across the region, but these climate changes have resulted in highly differentiated and nonlinear landscape responses. We argue this reflects nonlinearity in the forcing–response processes at such scales and suggests that current approaches to modelling the response of such systems to future climate change using numerical climate models may not accurately capture the landscape response. We end by discussing some of the implications for obtaining decision-relevant predictions of landscape responses to climatic forcing and for climate change adaptation and planning, using regional climate models. [ABSTRACT FROM AUTHOR]

Volcanic activity can have a notable impact on glacier behaviour (dimensions and dynamics). This is evident from the palaeo-record, but is often difficult to observe for modern glaciers. However, documenting and, if possible, quantifying volcanic impacts on modern glaciers is important if we are to predict their future behaviour (including crucial ice masses such as the West Antarctic Ice Sheet) and to monitor and mitigate glacio-volcanic hazards such as floods (including jökulhlaups) and lahars. This review provides an assessment of volcanic impacts on the behaviour of modern glaciers (since AD 1800) by presenting and summarising a global dataset of documented examples. The study reveals that shorter-term (days-to-months) impacts are typically destructive, while longer-term (years-to-decades) are more likely protective (e.g., limiting climatically driven ice loss). However, because these events are difficult to observe, particularly before the widespread availability of global satellite data, their frequency and importance are likely underestimated. The study also highlights that because the frequency and nature of volcano-glacier interactions may change with time (e.g., glacier retreat may lead to an increase in explosive volcanic activity), predicting their future importance is difficult. Fortunately, over coming years, continued improvements in remotely sensed data will increase the frequency, and enhance the quality, of observations of volcanic impacts on glaciers, allowing an improved understanding of their past and future operation. [ABSTRACT FROM AUTHOR]

Some surge-type glaciers on the High-Arctic archipelago of Svalbard have large glacitectonic composite ridge systems at their terrestrial margins. These have formed by rapid glacier advance into proglacial sediments during the active surge phase, creating multicrested moraine complexes. Such complexes can be formed during single surge advances or multiple surges to successively less-extensive positions. The few existing studies of composite ridge systems have largely relied on detailed information on internal structure and sedimentology to reconstruct their formation and links to surge processes. However, natural exposures of internal structure are commonly unavailable, and the creation of artificial exposures is often problematic in fragile Arctic environments. To compensate for these issues, we investigate the potential for reconstructing composite ridge system formation based on geomorphological evidence alone, focusing on clear morphostratigraphic relationships between ridges within the moraine complex and relict meltwater channels/outwash fans. Based on mapping at the margins of Finsterwalderbreen (in Van Keulenfjorden) and Grønfjordbreen (in Grønfjorden), we show that relict meltwater channels that breach outer parts of the composite ridge systems are in most cases truncated upstream within the ridge complex by an inner pushed ridge or ridges at their ice-proximal extents. Our interpretation of this relationship is that the entire composite ridge system is unlikely to have formed during the same glacier advance but is instead the product of multiple advances to successively less-extensive positions, whereby younger ridges are emplaced on the ice-proximal side of older ridges. This indicates that the Finsterwalderbreen composite ridge system has been formed by multiple separate advances, consistent with the cyclicity of surges. Being able to identify the frequency and magnitude of former surges is important as it provides insight into the past behaviour of surge-type glaciers and, if absolute dating is possible, allows for the assessment of surge-type glacier response to climate change on decadal to centennial timescales. Although further investigations into the internal structure of these deposits should be sought where possible, our study demonstrates that geomorphology could be an invaluable tool for reconstructing the formation of composite ridge systems. [ABSTRACT FROM AUTHOR]

A glacial geomorphological analysis of three valleys on Mt. Mavrovouni (North Pindus Mountains, Greece) is presented alongside a pilot study using cosmogenic 36Cl to obtain surface exposure ages from iron-rich ophiolite glacial and periglacial boulders. At least three distinct morphostratigraphic units of glacial (moraines) and periglacial (relict pronival ramparts) origin have been identified. Four 36Cl surface exposure ages were obtained from the stratigraphically youngest glacial and periglacial deposits. Although this limited dataset with relatively large uncertainties cannot support a robust geochronology, the ages are consistent with the 36Cl-based chronologies of limestone-derived moraines on Mt. Tymphi (NW Greece) and Mt. Chelmos (S Greece), confirming that the last glaciers on this massif formed during the Last Glacial Maximum as also indicated by other studies in the Pindus mountains. At the same time it provides confidence in the suitability of 36Cl dating for iron-rich samples, such as ophiolites, using an updated 36Cl model that incorporates improved production rates for iron spallation. The presented preliminary chronology of moraines and pronival ramparts is based on those ages as well as on local and regional morphostratigraphic correlations. The stabilisation of the most extensive Late Pleistocene glaciers took place during the Last Glacial Maximum, at 27.0 ± 6.5 ka whereas the presence of pronival ramparts dated at 20.2 ± 4.8 ka suggests persisting cold and arid conditions. Older, still undated glacial deposits exist lower in the valleys which can be attributed to the Middle Pleistocene major glaciation phases (MIS 12/MIS 6), based on their relative morphostratigraphic position within the glacial sedimentary sequence. • Benefits of ophiolite for dating glacial landforms compared with limestone • Pilot study of 36Cl surface exposure dating on ophiolites with high iron content • Preliminary glacial chronology using 36Cl ages and morphostratigraphic correlations • Two distinct glaciations ascribed to the Late and Middle Pleistocene respectively. • Last glaciation on Mt. Mavrovouni in NW Greece occurred around the global LGM. [ABSTRACT FROM AUTHOR]

Trough Mouth Fans (TMF) are sedimentary depocenters located at glaciated continental margins and consist predominantly of glacigenic debris flow deposits. The Foula wedge is a Pleistocene TMF accumulated offshore West of Shetland over the Northeast Atlantic margin. This study presents an analysis of a 3D seismic reflection dataset imaging the distal Foula wedge basin fan deposits between 1010 and 1100 m water depth, directly downslope from a gully system which was active untill the end of the last deglaciation. Results reveal, in unprecedented detail, the basal surface of this fan system and its internal complex architecture. Features typical of both debris flow deposits and turbidites are identified, including a basin channel network with linear and diverging erosional features forming distinctive terminal lobes, stacked and backstepping. The study links the seafloor morphology of the basin fan with its subsurface geomorphology, showing connection with the downslope gully system to the east. It presents evidence for a complex distal depositional system on glaciated margins, characterised by heterogeneous sediment delivery processes and deposits. A conceptual evolution model is proposed, with a glacigenic debris flow-dominated TMF at the LGM, subsequently influenced by meltwater discharges, with deposition occurring as a function of the shelf margin and slope paleo-morphology, slope substrate composition, interaction of downslope and along slope processes and ice-margin dynamics. • 3D architecture of a glacigenic basin fan system revealed in unprecedented detail. • Subsurface seismic geomorphology highlights a more heterogenous basin depositional system than previously interpreted. • Articulated basin channel network characterised by features typical of both debris flow deposits and turbidites. • Sediment delivery processes controlled by paleo-topography, substrate and meltwater in response to the ice-margin dynamics. [ABSTRACT FROM AUTHOR]

Abstract: Despite their importance in understanding glaciological processes and constraining large-scale flow patterns in palaeo-glaciology, there is little consensus as to how drumlins are formed. Attempts to solve the ‘drumlin problem’ often fail to address how they are created from an initially flat surface in the absence of obvious cores or obstacles. This is a key strength of the instability theory, which has been described in a suite of physically-based mathematical models and proposes that the coupled flow of ice and till causes spontaneous formation of relief in the till surface. Encouragingly, model predictions of bedform height and length are consistent with observations and, furthermore, the theory has been applied to a range of subglacial bedforms and not just drumlins. However, it has yet to confront the myriad observations relating to the composition and internal structure of drumlins and this could be seen as a major deficiency. This paper is a first attempt to assess whether the instability theory is compatible with the incredible diversity of sediments and structures found within drumlins. We summarise the underlying principles of the theory and then describe and attempt to explain the main types of drumlin composition (e.g. bedrock, till, glaciofluvial sediments, and combinations thereof). Contrary to a view which suggests that the presence of some sedimentary sequences (e.g. horizontally stratified cores) is inconsistent with the theory, we suggest that one would actually expect a diverse range of constituents depending on the inheritance of sediments that pre-date drumlin formation, the duration and variability of ice flow, and the balance between erosion and deposition (till continuity) at the ice–bed interface. We conclude that the instability theory is compatible with (and potentially strengthened by) what is known about drumlin composition and, as such, offers the most complete and promising solution to the drumlin problem to date. [Copyright &y& Elsevier]

The Antarctic continent reached its current polar location ~83 Ma and became shrouded by ice sheets ~34 Ma, coincident with dramatic global cooling at the Eocene-Oligocene boundary. However, it is not known whether the first Antarctic glaciers formed immediately prior to this or were present significantly earlier. Here we show that mountain glaciers were likely present in the Transantarctic Mountains during the Late Palaeocene (~60-56 Ma) and middle Eocene (~48-40 Ma). Temperate (warm-based) glaciers were prevalent during the Late Eocene (~40-34 Ma) and, in reduced numbers, during the Oligocene (~34-23 Ma), before larger, likely cold-based, ice masses (including ice sheets) dominated. Some temperate mountain glaciers were present during the Miocene Climatic Optimum (~15 Ma), before a widespread switch to cold-based glaciation. Our findings highlight the longevity of glaciation in Antarctica and suggest that glaciers were present even during the Early-Cenozoic greenhouse world., (© 2022. The Author(s).)