Landform Record Research Articles

Net evolution of subglacial sediment transport in the Quebec–Labrador sector of the Laurentide Ice Sheet

The Laurentide Ice Sheet's (LIS) interior had a dynamic polythermal base, but the spatiotemporal variations of subglacial processes related to ice divide migration and other transient changes remain largely unknown, limiting our understanding of regional glacial dynamics. Previous studies focused on the regional glacial landform record, while ice sheet models lacked detailed parameterization within these regions, leading to an overestimation of cold-based subglacial conditions' extent and duration. In this study, glacial sediment dispersal patterns as identified by heavy minerals, clasts, and multivariate statistics of till matrix geochemistry were used to assess ice sheet dynamics within the Quebec–Labrador sector of the LIS. The earliest ice-flow phase produced and transported till across the study area (>175 km). However, major oxide data from till matrix geochemistry show a correlation with underlying bedrock, and this relationship is relatively common in areas of thin till cover and resistant bedrock lithologies. These results suggest a switch from an early phase of widespread erosion and long, sustained sediment transport to one of more limited erosion, perhaps abrasion dominant and shorter transport. Till compositional data and related dispersal patterns add supporting evidence to earlier ice sheet reconstructions based on ice-flow indicators and 10Be data together suggesting a transition from widespread uniform warm-based conditions during the earliest ice flow, followed by the development of an ice divide, its migration, and more sporadic warm-based conditions. Consequently, a thorough understanding of ice-flow history is essential for ice sheet modelling and future mineral exploration programs in inner ice sheet regions of the LIS.

Footprint of the Baltic Ice Stream: geomorphic evidence for shifting ice stream pathways

The Baltic Ice Stream, a large fast‐flowing sector of the last Fennoscandian Ice Sheet that occupied the present‐day Baltic Sea basin, was first conceptualized in the earliest days of glacial geological research in Scandinavia. Landform and sedimentological evidence from the terrestrial margins support the concept and numerical ice‐sheet models demonstrate its existence and possible evolution. However, with evidence for the Baltic Ice Stream thus far limited to the terrestrial periphery, its true form, scale, function, and role in deglaciation have proven enigmatic. Here we present geomorphological evidence directly from the Baltic seabed that confirms the existence of and sheds light on the behaviour of the Baltic Ice Stream. Based on an extensive collection of high‐, moderate‐ and low‐resolution bathymetric terrain models covering a large proportion of the Baltic Sea floor, and complemented by LiDAR‐data for the Baltic islands, we have identified and mapped >20 000 individual subglacial bedforms, meltwater landforms and grounding line landforms. We reconstruct a six‐stage sequence of ice flow and retreat, finding that streaming was persistent in the Baltic but that pathways were variable in extent, timing and duration: different sectors of the Baltic exhibit asynchronous streaming and out‐of‐phase grounding line changes. During deglaciation, grounding line re‐advances occurred in both the southwestern and the northern Baltic Proper, and, while abundant iceberg ploughmarks attest to calving as a significant ice loss mechanism, lobate margins suggest supply to the Baltic catchment was consistently high. Our reconstruction is limited by a fragmentary geomorphic record. Here we put forward a first hypothesis for how the Baltic Ice Stream evolved, and hope it stimulates new geomorphic, stratigraphical and core data collection to extend the landform record, provide insights into subglacial and grounding line processes, and constrain the chronology for Baltic Ice Stream flow and retreat.

The surficial geology record of ice stream catchment dynamics and ice-divide migration in the Quebec-Labrador sector of the Laurentide Ice Sheet

The surficial geology of an approximately 30 000 km2 region in northeastern Quebec and western Labrador was mapped to elucidate the glacial dynamics of a region that experienced ice divide migrations and determine the net effect of this evolution on the glacial landscape. This region is located in one of the major inner-ice dispersal centres of the Laurentide Ice Sheet, the Quebec-Labrador Sector, and therefore provides important information on the overall dynamics of the interior of an ice sheet throughout a glaciation. Along with detailed surficial mapping, landform analysis was conducted to assess changes in ice stream catchments and related ice divide migration and understand their impact on the overall landform record. This work is a synthesis of a four-year mapping program that resulted in the completion of eight 1:100 000 scale surficial geology maps that included three years of field-based investigations. Mapping has revealed a complex landscape with a diverse suite of surficial deposits and a fragmented and overprinted mosaic of landforms, which provides evidence for polythermal subglacial conditions that evolved throughout glaciation.

Chronology and erosion rate of the Pinedale glaciation, Colorado Front Range (USA), inferred from the sedimentary record of glacial Lake Devlin

Glacial and periglacial sediments and landforms record the chronology of glaciation and amount of Pleistocene erosion during colder periods that added substantially to global sediment budgets and contributed to the global CO2 cycle. The now-drained glacial Lake Devlin, dammed in a Front Range tributary valley by a glacier in the North Branch of Boulder Creek (Colorado, USA) preserves an important sedimentary archive of the ca. 32−14 ka Pinedale glaciation, recording both paleoclimate information and an integrated measure of glacial and periglacial erosion rates over a full glacial cycle. Despite rapid erosion of fine-grained deposits after the lake drained, most sediment generated during Pinedale time remains as legacy deposits in the catchment. Geomorphic evidence and dating of glaciolacustrine sediment from surface exposures demonstrate that the ca. 30 ka Pinedale glacial advance was nearly as extensive as the local Late Glacial Maximum at ca. 20 ka. Sedimentary archives dated by 14C, optically stimulated luminescence, and cosmogenic nuclides extend earlier studies (Madole et al., 1973) of pollen and magnetic susceptibility (MS) in cores from the glaciolacustrine deposits of Lake Devlin and of Pinedale climate. Records suggest short-term warming and biotic change at ca. 15 ka after ∼14 kyr of cold, dry conditions punctuated by MS peaks at ca. 26.5 ka, 20 ka, and 16.5 ka. Lake Devlin drained catastrophically after ca. 14 ka, millennia after ice had retreated upvalley from the lateral moraine that dammed the lake. Sediment production during the Pinedale was equivalent to a periglacial and glacial erosion rate of ∼70 mm kyr−1, several times higher than long-term rates in the adjacent Front Range, but much lower than rates measured where modern glaciers are eroding weak bedrock in zones of rapid rock uplift, such as SSE Alaska, USA. Data from the Lake Devlin basin contribute to contemporary discussions of how glacial erosion influences the global CO2 cycle.

Assessing ice sheet models against the landform record: The Likelihood of Accordant Lineations Analysis (LALA) tool

AbstractPalaeo‐ice sheets leave behind a rich database regarding their past behaviour, recorded in the landscape in the form of glacial geomorphology. The most numerous landform created by these ice sheets are subglacial lineations, which generate snapshots of the direction of ice flow at fixed (yet typically unknown) points in time. Despite their relative density within the landform record, the information provided by subglacial lineations is currently underutilised in tests of numerical ice sheet models. To some extent, this is a consequence of ongoing debate regarding lineation formation, but predominantly, it reflects the lack of rigorous model‐data comparison techniques that would enable lineation information to be properly integrated. Here, we present the Likelihood of Accordant Lineations Analysis (LALA) tool. LALA provides a statistically rigorous measure of the log‐likelihood of a supplied ice sheet simulation through comparison of simulation output with both the location and direction of observed lineations. Given an ensemble of ice sheet simulations, LALA provides a formal, and statistically underpinned, quantitative assessment of each simulation's quality‐of‐fit to mapped lineations. This enables a comparison of each simulation's relative plausibility, including identification of the most likely ice sheet simulations amongst the ensemble. This is achieved by modelling lineation formation as a marked Poisson point process and comparison of observed to modelled flow directions using the von Mises distribution. LALA is flexible—users can adapt parameters to account for differing assumptions regarding lineation formation, and for variations in the level of precision required for differing model‐data comparison experiments. We provide guidelines and rationale for assigning parameter values, including an assessment of the variability between users when mapping lineations. Finally, we demonstrate the utility of LALA through application to an ensemble of simulations of the last British‐Irish Ice Sheet. This comparison highlights the benefits of LALA over previous tools and demonstrates some of the considerations of experimental design required when identifying the fit between ice sheet model simulations and the landform record.

MIS 3 age of the Veiki moraine in N Sweden – Dating the landform record of an intermediate-sized ice sheet in Scandinavia

ABSTRACT The Veiki moraine in northern Sweden, a geomorphologically distinct landscape of ice-walled lake plains, has been interpreted to represent the former margin of an intermediate-sized pre–Last Glacial Maximum (LGM) Fennoscandian ice sheet, but its age is debated as either marine isotope stage (MIS) 5c or MIS 3. We have applied optically stimulated luminescence (OSL) and radiocarbon dating to four sites within the northern part of the Veiki moraine to establish its chronology. The radiocarbon ages provide only minimum ages and most OSL ages have low precision due to poor luminescence characteristics and problems with incomplete bleaching, leading to two alternative ages. In either case, the OSL dating places the Veiki moraine formation in MIS 3 (best estimate 56–39 ka). Sedimentation continued in the low-lying centers of some plateaus (ice-walled lake plains) during MIS 3 and during the Holocene, with a break during the Last Glacial Maximum when the area was ice covered. We speculatively constrain the broad timing further by relating the sequence of events to other climate records. We suggest that ice margin retreat to the west of the Veiki area took place during Greenland Interstadial (GI) 16.1 (58.0–56.5 ka) and that limited ice advances, which led to debris-covered ice margins in the Veiki zone, occurred during the following stadials GS-16.1 to 15.1 (56.5–54.2 ka). The GI-14 interstadial, which began 54.2 ka and lasted ~5.9 ka, could then be the period when the ice within the dead-ice landscape melted, first leading to ice-walled lakes and later to the inversed topography characteristic of the Veiki landscape.

Reconstructing the advance and retreat dynamics of the central sector of the last Cordilleran Ice Sheet

The advance of the Cordilleran Ice Sheet (CIS) towards its Last Glacial Maximum (LGM) configuration and its subsequent retreat remain poorly understood. We use the glacial landform record to determine ice dynamics for the central sector of the CIS in northern British Columbia, Canada, beneath the LGM ice divide. We classify seventy ice-flow indicator flowsets based on morphology, elevation, orientation and cross-cutting relationships into one of three stages, whereby stage 1 is oldest and stage 3 youngest. Combined with ice-contact geomorphology, our reconstruction highlights complex changes in ice flow over time as a result of ice divide migrations through the LGM and deglacial phases. The orientation and distribution of landforms indicates active post-LGM ice retreat westward through the Cassiar and Omineca mountains. We map the regional distribution of independent mountain glaciers, ice caps, and ice fields that regrew during a cooling event in the Late Glacial and show that some of these readvance glaciers were subsequently overrun by advancing outlet glaciers of the CIS. We use the cross-cutting relationship between readvance glaciers and CIS outlet glaciers and available chronological data to reconstruct the eastern CIS margin during the Late Glacial for the first time.

The role of glacial dynamics in the development of ice divides and the Horseshoe Intersection Zone of the northeastern Labrador Sector of the Laurentide Ice Sheet

The Labrador Sector of the Laurentide Ice Sheet is characterized by a complex network of ice divides and an extensive landform record outlining two broad and opposing ice flows that are separated by a narrow Horseshoe Intersection Zone (HIZ). This geomorphic system gave rise to contrasting reconstructions, which reflect uncertainties on the temporal evolution of ice divides, the nature and age of the main landform systems and the overall pattern of ice retreat during the last deglaciation. Here, we address these issues through systematic mapping of glacial landforms and ice-flow indicators in a large area of northeastern Quebec and Labrador. The application of cosmogenic (10Be and 26Al) dating to esker-fed glaciomarine deltas and different rock surfaces brings new constraints on the chronological framework and insights on the subglacial thermal regime. Our results outline four main landform assemblages, two of which relating to a major system of opposite ice flows that delineate the position of the eastern Ancestral Labrador ice divide. Our reconstruction shows that the HIZ occurs several tens of km to the west of this divide and that these two features are genetically distinct. The HIZ relates to the late-glacial development of ice streams showing a massive convergent ice flow into Ungava Bay. The upstream extent of this ice flow system shows a strong coupling with the outline of Ungava Bay drainage divide, suggesting that the configuration of the HIZ could be the expression of a topographic control on the late-glacial ice sheet dynamics. Landforms also indicate a significant shift in the deglaciation mode, which evolved from an overall warm-based to an areally-confined cold-based ice retreat near the center of the former ice mass – a change in dynamics that also played a role in the configuration of the HIZ. Overall, areas indicative of cold-based ice conditions are limited to highly elevated terrains, thus ruling out models arguing for an extensive cold-based ice cover over Ungava Bay throughout the last deglaciation.

Quaternary landscape evolution of patagonia at the Chilean Triple Junction: Climate and tectonic forcings

Glacial and fluvial landforms record the recent history of the Earth surface and have the potential to provide information about uplift and climate changes during Quaternary times. Patagonia hosts uniquely well-preserved moraine deposits, along with outwash and fluvial terraces. This offers the possibility of determining the timing of glacier advances since the Great Patagonian Glaciation (ca. 1.2 Ma) and allows us to investigate landscape response to climate and tectonic forcing.In this study, we present twenty-eight new in situ 10Be exposure ages for glacial deposits, outwash sediments, and fluvial terraces at the latitude of the Chile Triple Junction (∼46°30′S). We combine this database with a detailed mapping and a morphological analysis of the Deseado River terrace system. The glacial and outwash deposits document at least five major glacial advances in the Lago Buenos Aires valley before the deposition of the Moreno system during Marine Isotope Stage (MIS) 8. Key advances occurred during MIS 12 and MIS 16. Further ice expansions that resulted in Telken V and VI/VII outwash surfaces date to 780 ka and ca. 1.2 Ma respectively. Overall, our data show the difficulty of dating pre-LGM glacial deposits using surface exposure methods and reinforce the idea that outwash surfaces are better targets to determine the age of glacial sequences in the region.The Deseado River terrace system yields stratigraphically consistent ages that range from 1167 to 447 ka, which indicate that these terraces formed during Pleistocene times. Our analysis suggests that terraces T3, T7 and T8 may have been formed during cold stages, close to the glacial-interglacial transitions. We also calculated incision rates that reveal a strong initial pulse, from ca. 1200 to 900 ka, which was followed by significantly weaker river incision pulses.Finally, our geomorphological analysis confirms that the oldest fluvial terraces of the Deseado River system exhibit a NE-SW tilting through a broad area, next to the Atlantic Ocean. This was previously interpreted as the surface response to a long-wavelength dynamic uplift, a phenomenon that should be reassigned to the early Pleistocene.

Using fluvial terraces as distributed deformation offset markers: Implications for deformation kinematics of the North Qilian Shan Fault

Knowledge of the kinematics of thrusts is the key to understanding mountain building processes in compressive ranges; however, relatively little attention has been paid to distributed deformation (on-fault and off-fault such as rotation and warping) across a thrust fault. Distributed deformation is widespread along the northern front of the Qilian Shan, in the NE Tibetan Plateau. In the Daciyao River, of the central North Qilian Shan, the East Yumu Shan Fault (EYF) splays into four parallel faults at the piedmont. The Daciyao River flows across the deformation zone and its terrace landform record provides a valuable geomorphic marker for reconstructing the cumulative deformation. Terrace deformation across the fault splays was estimated with measured longitudinal profiles using differential GPS, and the fault slip was reconstructed by applying the trishear model, a widely applied model used to explain the geometry of basement-related structures. At least 10 strath terraces are formed with bedrock bases developed into Neogene red beds along the mountain front. The age of one of the terraces (T8) was constrained by Optically Stimulated Luminescence (OSL) dating of the overlying aeolian loess, which provide a minimum age (68 ka) for the age of terrace abandonment. The survey data and trishear model results suggest that the fault splays are bending-moment faults (BMF), originating from a fault propagation fold. Based on the fault and fold geometry across the EYF, the total vertical slip rate is estimated as 0.9–1.2 mm/yr. Slip rates on the EYF are consistent along-strike during the Late Pleistocene but they are twice as high during the Holocene, indicating a high potential for large earthquakes along this fault. This work also indicates that to calculate accurate total slip rate, terrace heights must be measured outside of the distributed broad zone.

Proposing a new conceptual model for the reconstruction of ice dynamics in the SW sector of the Scandinavian Ice Sheet (SIS) based on the reinterpretation of published data and new evidence from optically stimulated luminescence (OSL) dating

Abstract. We propose a new concept of the Weichselian ice dynamics in the south-western sector of the Baltic Sea depression. The review of existing geochronological data from Germany, Denmark and southernmost Sweden in combination with new optically stimulated luminescence (OSL) data from the German Oder Lobe area is the basis for a reassessment and an improvement of previous ice dynamic models. Factors like the pre-existing topography, glaciotectonic features and the occurrence of till beds and inter-till deposits of varying origin are also taken into consideration for our process-based reconstruction of the sedimentary environments close to the ice margin and hence the ice dynamics of the Scandinavian Ice Sheet (SIS). During the early MIS 3 (marine isotope stage), the late MIS 3 and MIS 2, the SIS advanced into present-day terrestrial areas around the south-western Baltic Sea Basin. The first ice advance during the warming phase in early MIS 3 is poorly documented as the Ellund–Warnow Advance in Germany but may be correlated with the numerically dated Ristinge Advance in Denmark and Sweden. The late MIS 3 advance in contrast is reliably documented. It shaped the landforms of the Brandenburg Advance and the maximum Weichselian ice extent in the Oder Lobe area in north-eastern Germany and occurred contemporaneously with the Klintholm Advance in southern Sweden and Denmark. The lack of a corresponding till in various cliff profiles along the Baltic Sea coastline between southern Schleswig-Holstein and the island of Rügen can be explained by the distinct lobate structure of this ice advance, which was strongly guided by the pre-existing low-lying topography. We propose the horst of Bornholm, Denmark, acting as an ice divide, with ice-dammed lakes existing on the lee side between two glacier lobes. This lobate structure had not been considered in previous conceptual models, which led to seemingly conflicting chronological and stratigraphical interpretations. Our introduction of the lobate structure for the first time resolves these contradictions and integrates the data in a coherent model. The dynamics of the MIS 2 readvance to the Last Glacial Maximum (LGM) extent were clearly different to the previous advance and were most likely characterized by a more uniformly advancing ice front with a less lobate structure which also overrode the horst of Bornholm and the island of Rügen. This advance reached the maximum Weichselian ice extent in some parts of the south-western SIS, but, in the Oder Lobe area, it is proven to have terminated at a lesser extent than the early MIS 3 advance, but it did shape the most prominent morphological landform record of the last glacial cycle. In order to advance the reconstruction of Weichselian ice dynamics in the future, we strongly suggest using both an MIS-based terminology and a process-based approach in the interpretation of geochronological data to live up to the dynamic nature of continental ice sheets.

Late Cenozoic evolution of the Ariège River valley (Pyrenees) constrained by cosmogenic 26Al/10Be and 10Be/21Ne dating of cave sediments

Despite increasing knowledge on the orogenic phases of the Pyrenees, the Neogene evolution of the range remains poorly constrained. The central Pyrenees, particularly the Ariège River valley and its terrace systems and glacial extensions, are key to reconstructing Pyrenean evolution during the Neogene. However, few terrace relics remain on the piedmont edges. To overcome this limitation and temporally extend the dataset obtained from terraces, we focus on alluvium-filled horizontal epiphreatic passages developed in limestone karstic networks. These landforms record the transient position of former local base levels during the process of valley deepening, similar to fluvial terraces. The alluvium fills of the studied caves in the Tarascon-sur-Ariège area, in the transition zone between the upper range and the piedmont, therefore enable the reconstruction of the geodynamic evolution of the Ariège River valley. All studied caves are developed on at least eight levels. Based on burial durations determined by 26Al/10Be and 10Be/21Ne terrestrial cosmogenic nuclide dating of samples from epiphreatic levels, our results indicate Langhian to Messinian (~13–5 Ma) incision rates averaging 48 m Myr−1. However, the obtained record does not allow us to accurately retrace the Pliocene evolution of the area due to the lack of known caves at corresponding levels. Moreover, raised local base levels during glacial phases both make the record more complex and call into question the methodological approach in terms of potential internal sediment remobilization and mixing related to implied re-flooding periods.

The changing extent of marine-terminating glaciers and ice caps in northeastern Svalbard since the ‘Little Ice Age’ from marine-geophysical records

Climate warming in Svalbard since the end of the ‘Little Ice Age’ early in the 20th century has reduced glacier extent in the archipelago. Previous attempts to reconstruct ‘Little Ice Age’ glacier limits have encountered problems in specifying the area of tidewater glacier advances because it is difficult to estimate the past positions of their marine termini. Multibeam echo-sounding data are needed to map past glacier extent offshore, especially in open-marine settings where subaerial lateral moraines cannot be used due to the absence of fjord walls. We use the submarine glacial landform record to measure the recent limits of advance of over 30 marine-terminating northeastern Svalbard glaciers and ice caps. Our results demonstrate that previous work has underestimated the ice-covered area relative to today by about 40% for northeastern Svalbard (excluding southeast Austfonna) because marine-geophysical evidence in the form of submarine terminal moraines was not included. We show that the recent ice extent was 1753 km2larger than today over our full area of multibeam data coverage; about 5% of the total modern ice cover of Svalbard. It has often been assumed that moraine ridges located within a few kilometres of modern ice fronts in Svalbard represent either a ‘Little Ice Age’ maximum or relate to surge activity over the past century or so. In the marine environment of northeastern Svalbard, this timing can often be confirmed by reference to early historical maps and aerial photographs. Assemblages of submarine glacial landforms inshore of recently deposited terminal moraines suggest whether a recent advance may be a result of surging or ‘Little Ice Age’ climatic cooling relative to today. However, older terminal moraines do exist in the archipelago, as shown by radiocarbon and10Be dating of Holocene moraine ridges.

Evolution of the Norwegian plateau icefield Hardangerjøkulen since the ‘Little Ice Age’

The maximum ‘Little Ice Age’ (LIA) glacier extent provides a significant baseline to assess long-term glacier change and to place currently observed rates of glacier recession in a broader temporal context. To that end, we examine the evolution of the plateau icefield Hardangerjøkulen since the LIA. First, we reconstruct Hardangerjøkulen’s maximum LIA extent (~AD 1750) and subsequent recession based on the glacial landform record and aided by historical map interpretation. Ice-marginal moraines, glacial drift limits, trimlines, and identifiable erosion and weathering boundaries provide evidence of a LIA icefield with an area of 110 km2. Existing LIA model simulations of Hardangerjøkulen are not yet fully able to reproduce our reconstructed extent. Second, we compile a set of remotely sensed icefield outlines from successive time points in the 20th and 21st century to calculate icefield area and length change since the LIA. This reveals a substantial reduction in icefield size, with a total area loss of 41 km2 (37%; 2% 10 a–1) by 2010 and a cumulative frontal retreat averaging 1.3 km (29%; 5 m a–1) by 2013. Icefield recession has been greatest since the end of the 20th century, when rates of areal shrinkage increased to 6.5–10% 10 a–1 in 1995–2010, and the rate of average terminus retreat accelerated to 17 m a–1 in 2003–2010. Third, we present a relative dating approach, based on the known age of the different icefield outlines, that allows bracketing ages to be assigned to all ice-marginal landforms between any two outlines. This approach shows that episodes of moraine formation vary temporally between individual outlet glaciers of Hardangerjøkulen, suggesting that the moraine record of a single outlet glacier alone may not be sufficient to derive an icefield-wide picture of past ice advances, and thereby climate fluctuations.

High‐resolution bathymetric mapping reveals subaqueous glacial landforms in the Arctic alpine lake Tarfala, Sweden

ABSTRACTIn Arctic alpine regions, glacio‐lacustrine environments respond sensitively to variations in climate conditions, impacting, for example,glacier extent and rendering former ice‐contact lakes into ice distal lakes and vice versa. Lakefloors may hold morphological records of past glacier extent, but remoteness and long periods of ice cover on such lakes make acquisition of high‐resolution bathymetric datasets challenging. Lake Tarfala and Kebnepakte Glacier, located in the Kebnekaise mountains, northern Sweden, comprise a small, dynamic glacio‐lacustrine system holding a climate archive that is not well studied. Using an autonomous surface vessel, a high‐resolution bathymetric dataset for Lake Tarfala was acquired in 2016, from which previously undiscovered end moraines and a potential grounding line feature were identified. For Kebnepakte Glacier, structure‐from‐motion photogrammetry was used to reconstruct its shape from photographs taken in 1910 and 1945. Combining these methods connects the glacial landform record identified at the lakefloor with the centennial‐scale dynamic behaviour of Kebnepakte Glacier. During its maximum 20th century extent, attained c. 1910, Kebnepakte Glacier reached far into Lake Tarfala, but had retreated onto land by 1945, at an average of 7.9 m year–1. Copyright © 2019 John Wiley & Sons, Ltd.

The Tekapo Glacier, New Zealand, during the Last Glacial Maximum: An active temperate glacier influenced by intermittent surge activity

Quaternary glaciations have created impressive landform assemblages that can be used to understand palaeo-glacier extent, character and behaviour, and hence past global and local glacier forcings. However, in the southern hemisphere and especially in New Zealand, the Quaternary glacial landform record is relatively poorly investigated with regard to glaciological properties. In this study, a 1 m digital elevation model (DEM) was generated from airborne LiDAR data and supplemented with aerial imagery and field observations to analyse the exceptionally well-preserved glacial geomorphology surrounding Lake Tekapo, New Zealand. We describe a rich suite of Last Glacial Maximum (LGM) and recessional ice-marginal, subglacial, supraglacial, glaciofluvial and glaciolacustrine landform assemblages. These represent two landsystems comprising i) fluted till surfaces with low-relief push moraine ridges; and ii) crevasse-squeeze ridges, ‘zig-zag’ eskers and attenuated lineations. The former landsystem records the behaviour of an active temperate glacier and the latter landsystem, which is superimposed upon and inset within the former, strongly suggests intermittent surge phases. The two landsystem signatures indicate a sequential change in ice-marginal dynamics during recession that was likely to have been partially non-climatically driven. Overall, we present the first evidence of surge-type glacier behaviour in New Zealand.

Submarine geomorphology of the northeastern Baffin Island fiords and cross-shelf troughs

ABSTRACT This paper presents a series of 50 geomorphological maps of the seabed of northeastern Baffin Island fiords and cross-shelf troughs, in eastern Arctic Canada. The mapping was produced using swath bathymetry imagery from multiple bathymetry datasets. A total of 24 types of landform were systematically mapped (>55,000 landforms) and reflect processes linked to subglacial, ice-marginal or paraglacial to postglacial environments. The landforms record the transition of the seabed from being covered by a marine-terminating ice sheet to the establishment of postglacial conditions. The landform assemblages allow the distinction of ice-flows orientation and ice-stream pathways along fiords and cross-shelf troughs. The multiple moraines and grounding-zone wedges indicate that the ice margin stabilized during retreat and that the overall deglaciation occurred by steps. These maps provide a framework for future investigations in northeastern Baffin Island fiords and shelf and can also provide a template for future seabed geomorphological studies in Arctic Canada.

The Glacial Geomorphology of the Ice Cap Piedmont Lobe Landsystem of East Mýrdalsjökull, Iceland

A surficial geology and geomorphology map of the forelands of the Sandfellsjökull and Oldufellsjökull piedmont lobes of the east Mýrdalsjökull ice cap is used to characterise the historical and modern landscape imprint in a glacial landsystems context. This serves as a modern analogue for palaeoglaciological reconstructions of ice cap systems that operated outlet lobes of contrasting dynamics, but the subtle variability in process-form regimes is encoded in the geomorphology. The landsystems of the two piedmont lobes reflect significantly different process-form regimes, and hence contrasting historical glacier dynamics, despite the fact that they are nourished by the same ice cap. The Sandfellsjökull landsystem displays the diagnostic criteria for active temperate glacier operation, including arcuate assemblages of inset minor push moraines and associated flutings, kame terrace and ice-dammed lake deposits, linear sandar directed by overridden moraine arcs, and since 1945, features, such as ice-cored, pitted, and glacially pushed outwash fans that are linked to englacial esker networks representative of recession into an overdeepening. Moraine plan forms have also changed from weakly crenulated and discontinuous curvilinear ridges to sawtooth features and crevasse-squeeze ridges and till eskers in response to changing proglacial drainage conditions. The Oldufellsjökull landsystem displays subtle signatures of jökulhlaup-driven surges, including sparse and widely spaced moraine clusters that are separated by exceptionally long flutings. The subtlety of the surge imprint at Oldufellsjökull was recognised only by comparison with nearby Sandfellsjökull, suggesting that palaeo-surging has likely been under-estimated in the ancient landform record. Hence, the simple imprint of sparse and widely spaced moraine clusters that are separated by exceptionally long flutings should be included as possible surge-diagnostic criteria.

Timing and dynamics of glaciation in the Ikh Turgen Mountains, Altai region, High Asia

Abstract Spanning the northern sector of High Asia, the Altai region contains a rich landform record of glaciation. We report the extent, chronologies, and dynamics of two paleoglaciers on opposite flanks of the Ikh Turgen mountains (In Russian: Chikhacheva Range), straddling the border between Russia and Mongolia, using a combination of remote sensing-based glacial geomorphological mapping, 10Be surface exposure dating, and geomorphometric analysis. On the eastern side (Mongolia), the Turgen-Asgat paleoglacier, with its potential for developing a large accumulation area (∼257 km2), expanded 40 km down valley, and mean ages from a latero-frontal moraine indicate deglaciation during marine oxygen isotope stage (MIS) 3 (45.1 ± 1.8 ka, n = 4) and MIS 2 (22.8 ± 3.3 ka, n = 5). These minimum age constraints are consistent with other 10Be glacial chronologies and paleoclimate records from the region, which indicates glacier culmination during cold and wet conditions coinciding with MIS 3 (piedmont-style glaciation; inferred for a few sites across the region) and glacier culmination during cold and dry conditions coinciding with MIS 2 (mainly valley-style glaciation; inferred from several sites across the region). On the western side (Russia), the Boguty paleoglacier had a smaller accumulation area (∼222 km2), and advanced 30 km down valley across a low gradient forefield. Surface exposure ages from two moraine complexes on this side of the mountains exhibit wide scatter (∼14–53 ka, n = 8), making paleoclimate inferences and comparison to other proxies difficult. Ice surface profile reconstructions imply that the two paleoglaciers likely shared an ice divide.

BRITICE Glacial Map, version 2: a map and GIS database of glacial landforms of the last British–Irish Ice Sheet

During the last glaciation, most of the British Isles and the surrounding continental shelf were covered by the British–Irish Ice Sheet (BIIS). An earlier compilation from the existing literature (BRITICE version 1) assembled the relevant glacial geomorphological evidence into a freely available GIS geodatabase and map (Clark et al. 2004: Boreas 33, 359). New high‐resolution digital elevation models, of the land and seabed, have become available casting the glacial landform record of the British Isles in a new light and highlighting the shortcomings of the V.1 BRITICE compilation. Here we present a wholesale revision of the evidence, onshore and offshore, to produce BRITICE version 2, which now also includes Ireland. All published geomorphological evidence pertinent to the behaviour of the ice sheet is included, up to the census date of December 2015. The revised GIS database contains over 170 000 geospatially referenced and attributed elements – an eightfold increase in information from the previous version. The compiled data include: drumlins, ribbed moraine, crag‐and‐tails, mega‐scale glacial lineations, glacially streamlined bedrock (grooves, roches moutonnées, whalebacks), glacial erratics, eskers, meltwater channels (subglacial, lateral, proglacial and tunnel valleys), moraines, trimlines, cirques, trough‐mouth fans and evidence defining ice‐dammed lakes. The increased volume of features necessitates different map/database products with varying levels of data generalization, namely: (i) an unfiltered GIS database containing all mapping; (ii) a filtered GIS database, resolving data conflicts and with edits to improve geo‐locational accuracy (available as GIS data and PDF maps); and (iii) a cartographically generalized map to provide an overview of the distribution and types of features at the ice‐sheet scale that can be printed at A0 paper size at a 1:1 250 000 scale. All GIS data, the maps (as PDFs) and a bibliography of all published sources are available for download from: https://www.sheffield.ac.uk/geography/staff/clark_chris/britice.