The hydrochemistry of Lake Vostok and the potential for life in Antarctic subglacial lakes

Our understanding of Lake Vostok, the huge subglacial lake beneath the East Antarctic Ice Sheet, has improved recently through the identification of key physical and chemical interactions between the ice sheet and the lake. The north of the lake, where the overlying ice sheet is thickest, is characterized by subglacial melting, whereas freezing of lake water occurs in the south, resulting in ～210 m of ice accretion to the underside of the ice sheet. The accreted ice contains lower concentrations of the impurities normally found in glacier ice, suggesting a net transfer of material from meltwater into the lake. The small numbers of microbes found so far within the accreted ice have DNA profiles similar to those of contemporary surface microbes. Microbiologists expect, however, that Lake Vostok, and other subglacial lakes, will harbour unique species, particularly within the deeper waters and associated sediments. The extreme environments of subglacial lakes are characterized by high pressures, low temperatures, permanent darkness, limited nutrient availability, and oxygen concentrations that are derived from the ice that provides the meltwater. Copyright © 2002 John Wiley & Sons, Ltd.     

Chemical depletion of sediment under the Greenland Ice Sheet

The chemical composition of sediment sampled from a network of ice boreholes contacting the bed of the ablation zone of the Greenland Ice Sheet is compared to the composition of intact bedrock clasts. This sediment is enriched in silica and depleted in reactive cations compared to the underlying bedrock. In order to test whether these differences result from grain size biases either in sample collection or due to natural sorting, sediment samples were analyzed by grain size aliquots. Enrichment of silicon and depletion of cations is consistent across grain size classes and the compositions of bedrock and subglacial sediment are statistically separate. The difference in composition between subglacial sediment and rock aligns closely with the composition of dissolved solutes in waters sampled from the same field sites and is dissimilar to the composition of the sites' suspended sediment. This implies that chemical weathering rather than disproportionate physical removal of friable minerals is responsible for the compositional differences between rock and sediment. Mass balance analysis implies 3–10% of the sediment's mass is lost to solute dissolution, with approximately double that amount precipitated as clay minerals (a large portion of which may have been physically expelled). This result implies that temperate ice sheet subglacial environments may be more chemically active than previously realized. Copyright © 2016 John Wiley & Sons, Ltd.     

Evidence for seasonal subglacial outburst events at a polythermal glacier,Finsterwalderbreen, Svalbard

Bulk runoff and meteorological data suggest the occurrence of two meltwater outburst events at Finsterwalderbreen, Svalbard, during the 1995 and 1999 melt seasons. Increased bulk meltwater concentrations of Cl^? during the outbursts indicate the release of snowmelt from storage. Bulk meltwater hydrochemical data and suspended sediment concentrations suggest that this snowmelt accessed a chemical weathering environment characterized by high rock:water ratios and long rock–water contact times. This is consistent with a subglacial origin. The trigger for both the 1995 and 1999 outbursts is believed to be high rates of surface meltwater production and the oversupply of meltwater to areas of the glacier bed that were at the pressure melting point, but which were unconnected to the main subglacial drainage network. An increase in subglacial water pressure to above the overburden pressure lead to the forcing of a hydrological connection between the expanding subglacial reservoir and the ice‐marginal channelized system. The purging of ice blocks from the glacier during the outbursts may indicate the breach of an ice dam during connection. Although subglacial meltwater issued continually from the glacier terminus via a subglacial upwelling during both melt seasons, field observations showed outburst meltwaters were released solely via an ice‐marginal channel. It is possible that outburst events are a seasonal phenomenon at this glacier and reflect the periodic drainage of meltwaters from the same subglacial reservoir from year to year. However, the location of this reservoir is uncertain. A 100 m high bedrock ridge traverses the glacier 6·5 km from its terminus. The overdeepened area up‐glacier from this is the most probable site for subglacial meltwater accumulation. Copyright © 2001 John Wiley & Sons, Ltd.     

The hydrochemistry of runoff from a ‘cold-based’ glacier in the High Arctic (Scott Turnerbreen,Svalbard)

There are still relatively few hydrochemical studies of glacial runoff and meltwater routing from the high latitudes, where non-temperate glacier ice is frequently encountered. Representative samples of glacier meltwater were obtained from Scott Turnerbreen, a ‘cold-based’ glacier at 78° N in the Norwegian high Arctic archipelago of Svalbard, during the 1993 melt season and analysed for major ion chemistry. Laboratory dissolution experiments were also conducted, using suspended sediment from the runoff. Significant concentrations of crustal weathering derived SO²⁻₄ are present in the runoff, which is characterized by high ratios of SO²⁻₄: (SO²⁻₄+HCO⁻₃) and high p(CO₂). Meltwater is not routed subglacially, but flows to the glacier terminus through subaerial, ice marginal channels, and partly flows through a proglacial icing, containing highly concentrated interstitial waters, immediately afront the terminus. The hydrochemistry of the runoff is controlled by: (1) seasonal variations in the input of solutes from snow- and icemelt; (2) proglacial solute acquisition from the icing; and (3) subaerial chemical weathering within saturated, ice-cored lateral moraine adjoining drainage channels at the glacier margins, sediment and concentrated pore water from which is entrained by flowing meltwater. Diurnal variations in solute concentration arise from the net effects of variable sediment pore water entrainment and dilution in the ice marginal streams. Explanation of the hydrochemistry of Scott Turnerbreen requires only one major subaerial flow path, the ice marginal channel system, in which seasonally varying inputs of concentrated snowmelt and dilute icemelt are modified by seepage or entrainment of concentrated pore waters from sediment in lateral moraine, and by concentrated interstitial waters from the proglacial icing, supplied by leaching, slow drainage at grain intersections or simple melting of the icing itself. The ice marginal channels are analogous neither to dilute supra/englacial nor to concentrated subglacial flow components. © 1998 John Wiley & Sons, Ltd.     

Quaternary ice thinning of David Glacier in the Terra Nova Bay region,Antarctica

Understanding the history of Antarctic glaciation is important for interpreting paleoclimatic changes and estimating the changes in climate, sea level, and ice volume in the future. Ice core studies of the East Antarctic Ice Sheet (EAIS) and marine sediment cores from the entire Ross Sea have employed numerous proxies to reconstruct the glacial history of the Antarctic region. However, the ice and marine core records can be biased because of their specific locations, such as the uppermost accumulation zone or the terminus of the ablation zone, thereby introducing significant uncertainties in ice modeling. In this study, we analyzed 34 new ¹⁰Be and ²⁶Al samples from four benches that were glaciated in the past by David glacier and incorporate the present ice-free flat surfaces. We suggest that the David glacier experienced monotonic and stepwise vertical lowering along the flanks of Mt. Priestley since the early Pleistocene. The uppermost bedrock benches on Mt. Priestley were exposed at 1.77 ± 0.32 Ma, with no evidence of subsequent overriding by readvancing ice. At Mt. Priestley, the David glacier has been characterized by a cold-based regime since 1.77 Ma, with a denudation rate of only ∼16 cm/Ma, corresponding to the regional transition from warm to cold-based glaciation at 3.5 Ma. Simple exposure ages from two lower benches date to Marine Isotope Stage (MIS) 7 (234.1 ± 13.1 ka; 545 m asl) and MIS 4 (64.8 ± 13.7 ka; 222 m asl), suggesting that, since MIS 8, the overall lowering of glaciers has remained monotonic. The upper bench marks the lower limit of the MIS 8 glacial period and the upper limit of Penultimate Glacial Maximum (MIS 6), while the lower landform defines the upper limit of the last glacial period (MIS 4–2). The magnitude of Quaternary ice thinning at the David glacier was the highest (∼990 m) in the present terminal area (i.e., the most sensitive ablation zone), in contrast to the other outlet glaciers draining into the Terra Nova Bay, which experienced less ice lowering. The combination of the terrestrial (in situ ¹⁰Be and ²⁶Al) and previous marine (authigenic ¹⁰Be) cosmogenic data used in our study document the history of lowering of the David glacier driven by climatic changes during the Pleistocene. Both deglaciation and glaciation were limited during the mid-Pleistocene transition (MPT) and prior to the mid-Bruhnes event (MBE), due to the prevailing cold and arid climate, whereas deglaciation was dominant during other warm periods.     

Landscape preservation under Fennoscandian ice sheets determined from in situ produced Be and Al

Some areas within ice sheet boundaries retain pre-existing landforms and thus either remained as ice free islands (nunataks) during glaciation, or were preserved under ice. Differentiating between these alternatives has significant implications for paleoenvironment, ice sheet surface elevation, and ice volume reconstructions. In the northern Swedish mountains, in situ cosmogenic ¹⁰Be and ²⁶Al concentrations from glacial erratics on relict surfaces as well as glacially eroded bedrock adjacent to these surfaces, provide consistent last deglaciation exposure ages (∼8-13 kyr), confirming ice sheet overriding as opposed to ice free conditions. However, these ages contrast with exposure ages of 34-61 kyr on bedrock surfaces in these same relict areas, demonstrating that relict areas were preserved with little erosion through multiple glacial cycles. Based on the difference in radioactive decay between ²⁶Al and ¹⁰Be, the measured nuclide concentration in one of these bedrock surfaces suggests that it remained largely unmodified for a minimum period of 845₋₄₁₈⁺⁴⁶¹ kyr. These results indicate that relict areas need to be accounted for as frozen bed patches in basal boundary conditions for ice sheet models, and in landscape development models. Subglacial preservation also implies that source areas for glacial sediments in ocean cores are considerably smaller than the total area covered by ice sheets. These relict areas also have significance as potential long-term subglacial biologic refugia.     

“Difference Dating”: A novel approach towards dating alpine glacial moraines

The East Antarctic Ice Sheet responds sluggishly to shifts in climate. To capture subtle changes in Antarctic climate, researchers have focused instead on smaller alpine and cirque glaciers that fringe the ice sheet throughout the McMurdo Dry Valleys. The exposure ages of glacial moraine boulders scatter widely and often incorporate large amounts of inheritance, prohibiting the construction of a regional deglaciation chronology. We present a new sampling technique that takes advantage of ubiquitous desert pavements and allows for the detection of inheritance in overlying glacial moraine boulders. Our approach requires a large sample set, but offers increased confidence in modeling moraine age, an acceptable trade-off considering the need for more refined Antarctic paleoclimate reconstructions. Using the beryllium-10 system in sandstone quartz, we show that single exposure ages collected from moraine boulder tops are frequently inaccurate, and consistently over- and underestimate moraine age. Difference Dating offers a new approach to dating alpine glacial moraines independent from traditional boulder exposure age dating.     

Subglacial Preservation of Valley Morphology at Amundsenisen,Western Dronning Maud Land,Antarctica

On the high altitude polar plateau of Amundsenisen, western Dronning Maud Land, East Antarctica, a subglacial valley, with a broad horizontal valley floor interpreted as a sediment floodplain or valley delta, was studied by radio echo sounding. In addition, a small, probably glacial, valley was mapped within the same subglacial massif. Basal ice temperatures were calculated using field data on precipitation, air temperature and ice sheet thickness. Discoveries of old landforms which have been preserved more or less intact beneath the former Fennoscandian and Laurentide ice sheets have received increasing attention during the last decade. The aim of this study is to investigate whether preservation of landforms occurs under the East Antarctic Ice Sheet, and to discuss under that climatological and glaciological circumstances preservation may take place. The results show that the ice sheet covering the investigated localities is frozen to bed, and therefore has an insignificant erosional capability. The observations suggest that a large-scale subglacial sediment deposit and a small valley formed by glacial erosion have survived beneath a cold-based ice sheet marginal zone for a long time period. The process of glacial preservation, recognized for bedrock features and tentatively observed for sediment accumulations, should act on similar large-scale landforms under any cold-based ice sheet, present or past. On the basis of existing studies of the age and stability of the East Antarctic Ice Sheet, a Middle Pliocene age is suggested for the preserved landforms. The presence of the presumed sediment-filled valley further indicates that no prolonged periods of basal melting have occurred at the Amundsenisen study area during the ice sheet history, which includes the Quaternary glaciation periods. Finally, calculations of basal temperature for localities at different altitudes within the same subglacial massif were used to demonstrate local altitudinal control of glacial preservation. © 1997 by John Wiley & Sons, Ltd.     

Distribution of ice thickness and subglacial topography of the “Chinese Wall” around Kunlun Station,East Antarctica

As fundamental parameters of the Antarctic Ice Sheet, ice thickness and subglacial topography are critical factors for studying the basal conditions and mass balance in Antarctica. During CHINARE 24 (the 24th Chinese National Antarctic Research Expedition, 2007/08), the research team used a deep ice-penetrating radar system to measure the ice thickness and subglacial topography of the “Chinese Wall” around Kunlun Station, East Antarctica. Preliminary results show that the ice thickness varies mostly from 1600 m to 2800 m along the “Chinese Wall”, with the thickest ice being 3444 m, and the thinnest ice 1255 m. The average bedrock elevation is 1722 m, while the minimum is just 604 m. Compared with the northern side of the ice divide, the ice thickness is a little greater and the subglacial topography lower on the southern side, which is also characterized by four deep valleys. We found no basal freeze-on ice in the Gamburtsev Subglacial Mountains area, subglacial lakes, or water bodies along the “Chinese Wall”. Ice thickness and subglacial topography data extracted from the Bedmap 2 database along the “Chinese Wall” are consistent with our results, but their resolution and accuracy are very limited in areas where the bedrock fluctuates intensely. The distribution of ice thickness and subglacial topography detected by ice-penetrating radar clarifies the features of the ice sheet in this “inaccessible” region. These results will help to advance the study of ice sheet dynamics and the determination of future locations of the GSM’s geological and deep ice core drilling sites in the Dome A region.     

Preserved landscapes underneath the Antarctic Ice Sheet reveal the geomorphological history of Jutulstraumen Basin

The landscape of Antarctica, hidden beneath kilometre-thick ice in most places, has been shaped by the interactions between tectonic and erosional processes. The flow dynamics of the thick ice cover deepened pre-formed topographic depressions by glacial erosion, but also preserved the subglacial landscapes in regions with moderate to slow ice flow. Mapping the spatial variability of these structures provides the basis for reconstruction of the evolution of subglacial morphology. This study focuses on the Jutulstraumen Glacier drainage system in Dronning Maud Land, East Antarctica. The Jutulstraumen Glacier reaches the ocean via the Jutulstraumen Graben, which is the only significant passage for draining the East Antarctic Ice Sheet through the western part of the Dronning Maud Land mountain chain. We acquired new bed topography data during an airborne radar campaign in the region upstream of the Jutulstraumen Graben to characterise the source area of the glacier. The new data show a deep relief to be generally under-represented in available bed topography compilations. Our analysis of the bed topography, valley characteristics and bed roughness leads to the conclusion that much more of the alpine landscape that would have formed prior to the Antarctic Ice Sheet is preserved than previously anticipated. We identify an active and deeply eroded U-shaped valley network next to largely preserved passive fluvial and glacial modified landscapes. Based on the landscape classification, we reconstruct the temporal sequence by which ice flow modified the topography since the beginning of the glaciation of Antarctica.     

The 1996 eruption at Gjálp,Vatnajökull ice cap,Iceland: efficiency of heat transfer,ice deformation and subglacial water pressure

The 13-day-long Gjálp eruption within the Vatnajökull ice cap in October 1996 provided important data on ice–volcano interaction in a thick temperate glacier. The eruption produced 0.8 km³ of mainly volcanic glass with a basaltic icelandite composition (equivalent to 0.45 km³ of magma). Ice thickness above the 6-km-long volcanic fissure was initially 550–750 m. The eruption was mainly subglacial forming a 150–500 m high ridge; only 2–4% of the volcanic material was erupted subaerially. Monitoring of the formation of ice cauldrons above the vents provided data on ice melting, heat flux and indirectly on eruption rate. The heat flux was 5–6×10⁵ W m^-2 in the first 4 days. This high heat flux can only be explained by fragmentation of magma into volcanic glass. The pattern of ice melting during and after the eruption indicates that the efficiency of instantaneous heat exchange between magma and ice at the eruption site was 50–60%. If this is characteristic for magma fragmentation in subglacial eruptions, volcanic material and meltwater will in most cases take up more space than the ice melted in the eruption. Water accumulation would therefore cause buildup of basal water pressure and lead to rapid release of the meltwater. Continuous drainage of meltwater is therefore the most likely scenario in subglacial eruptions under temperate glaciers. Deformation and fracturing of ice played a significant role in the eruption and modified the subglacial water pressure. It is found that water pressure at a vent under a subsiding cauldron is substantially less than it would be during static loading by the overlying ice, since the load is partly compensated for by shear forces in the rapidly deforming ice. In addition to intensive crevassing due to subsidence at Gjálp, a long and straight crevasse formed over the southernmost part of the volcanic fissure on the first day of the eruption. It is suggested that the feeder dyke may have overshot the bedrock–ice interface, caused high deformation rates and fractured the ice up to the surface. The crevasse later modified the flow of meltwater, explaining surface flow of water past the highest part of the edifice. The dominance of magma fragmentation in the Gjálp eruption suggests that initial ice thickness greater than 600–700 m is required if effusive eruption of pillow lava is to be the main style of activity, at least in similar eruptions of high initial magma discharge.Editorial responsibility: J. Donnelly-Nolan     

Sedimentary processes at ice sheet grounding-zone wedges revealed by outcrops,Washington State (USA)

Grounding-zone wedges (GZWs) mark the grounding terminus of flowing marine-based ice streams and, in the presence of an ice shelf, the transition from grounded ice to floating ice. The morphology and stratigraphy of GZWs is predominantly constrained by seafloor bathymetry, seismic data, and sediment cores from deglaciated continental shelves; however, due to minimal constraints on GZW sedimentation processes, there remains a general lack of knowledge concerning the production of these landforms. Herein, outcrop observations are provided of GZWs from Whidbey Island in the Puget Lowlands (Washington State, USA). These features are characterized by prograded diamictons bounded by glacial unconformities, whereby the lower unconformity indicates glacial advance of the southern Cordilleran Ice Sheet and the upper unconformity indicates locally restricted ice advance during GZW growth; the consistent presence of an upper unconformity supports the hypothesis that GZWs facilitate ice advance during landform construction. Based on outcrop stratigraphy, GZW construction is dominated by sediment transport of deformation till and melt-out of entrained basal debris at the grounding line. This material may be subsequently remobilized by debris flows. Additionally, there is evidence for subglacial meltwater discharge at the grounding line, as well as rhythmically bedded silt and sand, indicating possible tidal pumping at the grounding line. A series of GZWs on Whidbey Island provides evidence of punctuated ice sheet movement during retreat, rather than a rapid ice sheet lift-off. The distance between adjacent GZWs of 10²–10³ m and the consistency in their size relative to modern ice stream grounding lines suggests that individual wedges formed over decades to centuries. © 2018 John Wiley & Sons, Ltd.     

Chemical weathering since the Little Ice Age recorded in lake sediments: a high‐resolution proxy of past climate

Because of the different geochemical behaviour of rubidium and strontium in earth surface processes, variations of the Rb/Sr ratios in lake sediments were used as a geochemical proxy of chemical weathering and past climate in a single watershed. Low magnetic susceptibility, low CaCO₃, low Sr concentration and, hence, high Rb/Sr ratio in the lake sediments indicate weak chemical weathering under a cold but wet climate during the Little Ice Age (LIA) in the closed Daihai Lake watershed. The concordant change in both Sr and CaCO₃ concentrations with δ¹⁸O values in the Dunde ice core suggests that weak chemical weathering during the wet LIA was controlled by air temperature. After the LIA, however, precipitation played a dominant role in chemical weathering. Copyright © 2001 John Wiley & Sons, Ltd.     

The hydrochemistry of meltwaters draining a polythermal-based,high Arctic glacier,south Svalbard: I. The ablation season

Solute and runoff time-series at Finsterwalderbreen, Svalbard, provide evidence for considerable basal routing of water and the existence of at least two contrasting subglacial chemical weathering environments. The hydrochemistry of a subglacial upwelling provides evidence for a snowmelt-fed subglacial reservoir that dominates bulk runoff during recession flow. High concentrations of Cl⁻ and crustal ions, high pCO₂ and ratios of [*SO²⁻₄/(*SO²⁻₄+HCO⁻₃)] close to 0·5 indicate the passage of snowmelt through a subglacial weathering environment characterized by high rock:water ratios, prolonged residence times and restricted access to the atmosphere. At higher discharges, bulk runoff becomes dominated by icemelt from the lower part of the glacier that is conveyed through a chemical weathering environment characterized by low rock:water ratios, short residence times and free contact with atmospheric gases. These observations suggest that icemelt is routed via a hydrological system composed of basal/ice-marginal, englacial and supraglacial components and is directed to the glacier margins by the ice surface slope. Upwelling water flows relatively independently of icemelt to the terminus via a subglacial drainage system, possibly constituting flow through a sediment layer. Cold basal ice at the terminus forces it to take a subterranean routing in its latter stages. The existence of spatially discrete flow paths conveying icemelt and subglacial snowmelt to the terminus may be the norm for polythermal-based glaciers on Svalbard. Proglacial mixing of these components to form the bulk meltwaters gives rise to hydrochemical trends that resemble those of warm-based glaciers. These hydrochemical characteristics of bulk runoff have not been documented on any other glacier on Svalbard to date and have significance for understanding interactions between thermal regime and glacier hydrology. © 1998 John Wiley & Sons, Ltd.     

An erosional origin for drumlins of NW Poland

Drumlins are landforms essential to understanding of ice sheet movement over soft beds, sediment transport along the ice/bed interface, and the formation of a wide range of glacial deposits. Although investigated more than any other glacial landform, the origin of drumlins remains contentious. Using high-resolution LiDAR imagery and field data, we investigate the geomorphology and internal composition of one of the biggest drumlin fields in the North European Lowland. The Stargard drumlin field consists of over 1300 drumlins and related streamlined subglacial bedforms in a terminal part of a major Weichselian palaeo-ice stream of the southern Scandinavian Ice Sheet. The drumlins are typically 600-800 m long, 200-250 m wide, 3-6 m high and have axial elongation ratios ~2 but in some cases exceeding 15. Several subzones inferred from drumlin morphometry exist reflecting different ice flow dynamics. The most elongated drumlins occur in areas where ice moved down-slope and where thick fine-grained deposits of low hydraulic conductivity occur in the substratum. The largest portion of land occupied by drumlins and the greatest frequency density of drumlins occur where the ice moved up-slope. Stargard drumlins are composed of a wide variety of glacial deposits including various types of tills and meltwater sediments, which range from undisturbed to heavily deformed. There is no correlation between the deposits in the drumlins and the drumlin forms indicating that the deposits pre-date the drumlinizing process. It is suggested that the drumlin field was generated by a combination of direct glacial erosion and subglacial meltwater erosion by removing antecedent material from the inter-drumlin areas and streamlining the resultant bumps. Our data support the search for a unifying theory of drumlin formation and suggest erosion as the most plausible single mechanism generating drumlin landscapes. © 2019 John Wiley & Sons, Ltd.     

Paired 10Be sampling of polished bedrock and erratic boulders to improve dating of glacial landforms: an example from the Western Alps

Cosmogenic nuclide dating of glacial landforms may lead to ambiguous results for ice retreat histories. The persistence of significant cosmogenic concentrations inherited from previous exposure may increase the apparent exposure ages for polished bedrocks affected by limited erosion under ice and for erratic boulders transported by glaciers and previously exposed in high-altitude rock walls. In contrast, transient burying by moraines, sediments and snow decreases the apparent exposure age. We propose a new sampling strategy, applied to four sites distributed in the Arc and Arve valleys in the Western Alps, to better constrain the factors that can bias exposure ages associated with glacial processes. We used the terrestrial cosmogenic nuclide ¹⁰Be (TCN) to estimate the exposure time from paired sampling of depth profiles in polished bedrock and on overlying erratic boulders. For a given sampling site, the exposure ages for both the polished bedrock and boulder are expected to be the same. However, in six cases out of seven, boulders had significantly higher ¹⁰Be surface concentrations than those of the associated polished surfaces. In present and past glacial processes, the ¹⁰Be distribution with depth for boulders and bedrocks implies the presence of an inheritance concentration of ¹⁰Be. Our study suggests that ¹⁰Be concentrations in erratic boulders and in polished bedrocks provide maximum and minimum exposure ages of the glacial retreat, respectively. © 2019 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

A GIS-based multi-proxy analysis of the evolution of subglacial dynamics of the Quebec–Labrador ice dome,northeastern Quebec,Canada

The glacial records of the inner-core regions of the Laurentide Ice Sheet (LIS) document complex yet coherent patterns reflecting ice-sheet change (e.g. ice-divide migration), providing unique insights into past glacial conditions. This study develops a conceptual model of subglacial dynamics evolution within a major ice-dispersal centre of the LIS in northeastern Quebec, Canada using a GIS-based analysis of the surficial geologic record. Multiple proxies of subglacial conditions (subglacial streamlined landforms, lake density and lake area over thin drift/bedrock) were analysed through grid-overlay techniques and then classified based on different proxy variables ranging from highly mobile warm-based to immobile cold-based conditions. An additional proxy (till blanket) was used to identify areas of thick till deposition, but with few proxies (few lake or landform metrics). Based on local ice-flow reconstructions, the most ‘relict’ glacial terrain zone (GTZ1) has warm-based conditions over 66% of its area and is remarkably well preserved, suggesting laterally extensive warm-based conditions during the oldest identified ice-flow phase. This relict glacial terrain is partially overprinted by two subsequent ice-flow phases in spatially restricted zones in the northeast (73% warm-based), east-central (41% warm-based), and northwest (33% warm-based) of the study area. A zone of more sluggish conditions (only 3% warm-based) was identified in the highlands at the centre of the study area, characterized by thin till cover, few landforms, yet with large patches of relatively abundant small lakes, indicative of areal scouring. No clear evidence of sustained cold-based conditions (i.e. high chemical index of alteration values or high ¹⁰Be abundances) was found in the study area. These results suggest that warm-based conditions (active erosion and/or deposition) were uniformly widespread during the earliest ice-flow phase, later becoming more spatially restricted with broader sluggish ice conditions. These spatially restricted regions of warm-based subglacial regimes were likely controlled by surrounding and down-flow ice streaming. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Meteoritic and bedrock constraints on the glacial history of Frontier Mountain in northern Victoria Land, Antarctica

In 2001, a small H4 chondrite, Frontier Mountain (FRO) 01149, was found on a glacially eroded surface near the top of Frontier Mountain, Antarctica, about 600 m above the present ice level. The metal and sulphides are almost completely oxidized due to terrestrial weathering. We used a chemical leaching procedure to remove weathering products, which contained atmospheric ¹⁰Be and ³⁶Cl in a ratio similar to that found in Antarctic ice. The FRO 01149 meteorite has a terrestrial age of 3.0 ± 0.3 Myr based on the concentrations of the cosmogenic radionuclides ¹⁰Be, ²⁶Al and ³⁶Cl. This age implies that FRO 01149 is the oldest stony meteorite (fossil meteorites excluded) discovered on Earth. The noble gas cosmic ray exposure age of FRO 01149 is ~ 30 Myr. The meteorite thus belongs to the 33 Myr exposure age peak of H-chondrites.The bedrock surface on which FRO 01149 was found has wet-based glacial erosional features recording a former high-stand of the East Antarctic ice sheet. This ice sheet evidently overrode the highest peaks (> 2800 m a.s.l.) of the inland sector of the Transantarctic Mountains in northern Victoria Land. We argue that FRO 01149 was a local fall and that its survival on a glacially eroded bedrock surface constrains the age of the last overriding event to be older than ~ 3 Myr. The concentrations of in-situ produced cosmogenic ¹⁰Be, ²⁶Al and ²¹Ne in a glacially eroded bedrock sample taken from near the summit of Frontier Mountain yield a surface exposure age of 4.4 Myr and indicate that the bedrock was covered by several meters of snow. The exposure age is also consistent with bedrock exposure ages of other summit plateaus in northern Victoria Land.     

Subglacial abrasion rates at Goldbergkees,Hohe Tauern,Austria, determined from cosmogenic 10Be and 36Cl concentrations

We report concentrations of cosmogenic ¹⁰Be and ³⁶Cl used to determine erosion depths in the recently deglaciated bedrock at Goldbergkees in the Eastern Alps. The glacier covered the sampling sites during the Little Ice Age (LIA) until c. 1940. The youngest ages calculated from these concentrations match the known exposure time after the post‐LIA exposure of <100 years. The apparent age (no cover, no erosion) of most samples, however, is significantly older. We show that the measured nuclide concentrations represent subglacial erosion depths, rather than exposure times. In particular, erosion depths calculated using ¹⁰Be and ³⁶Cl concentrations of individual samples match well, whereas apparent ³⁶Cl ages are consistently older than ¹⁰Be ages. The bedrock at the ‘youngest’ surfaces was deeply eroded (≥ 297 cm) by the Goldbergkees during the late Holocene. In contrast, bedrock at the margin of the LIA ice extent was eroded ≤35 cm. These values convert to subglacial erosion rates on the order of 0.1 mm/a to >5 mm/a. While modeled erosion rates depend on the duration of glacial cover and erosion intrinsic to the different exposure scenarios used for calculation (700–3300 years), modeled total erosion depths are insensitive (5–20% change). Analysis of erosion depths on the transverse valley profile shows a general trend of greatest erosion part way up the valley side and less erosion under thin ice at the lateral margin. A second profile along the valley axis indicates depth of erosion is greatest where the ice abuts the foot of the investigated bedrock riegel and at its lee side just beyond the crest. Copyright © 2016 John Wiley & Sons, Ltd.     

Weak chemical weathering during the Little Ice Age recorded by lake sediments

Low magnetic susceptibility, low Sr content and hence high Rb/Sr ratio in the lake sediment sequence indicate a weak chemical weathering process under arid and cold climate of the Little Ice Age in a single closed lake watershed. According to different geochemical behavior between rubidium and strontium in earth surface processes, variation of Rb/Sr ratios in the lake sediment sequence can be used as an effective geochemical proxy with definite climatic significance of chemical weathering in watershed. Unlike chemical weathering process in tropic zone and modern temperate-humid climate, concordant changes in both Sr content and magnetic susceptibility with δ¹⁸O values of Dunde ice core suggest that the weak chemical weathering was controlled by air temperature during the Little Ice Age maximum. After the Little Ice Age, chemical weathering intensity was controlled also gradually by precipitation with increasing in temperature.