Meteoric 10Be concentrations from saprolite and till in northern Sweden: Implications for glacial erosion and age

We examine ¹⁰Be concentration in two pit profiles in the Parkajoki area at ∼67°N on the northern Fennoscandian shield in northern Sweden. Due to repeated cover by cold-based, non-erosive ice sheets, the area retains many relict non-glacial features, including tors and saprolites. In the examined pit profiles, gruss-type saprolite developed from weathering of intermediate igneous rocks is overlain unconformably by Weichselian till.Our results show that ¹⁰Be concentrations found in the till greatly exceed the levels of ¹⁰Be that can have accumulated since deglaciation at ∼11 ka and are comparable to those reported from Pliocene and Early Pleistocene tills in North America. Old tills with grussified boulders at depth were excavated in the Parkajoki area and correlations with neighbouring parts of Finland indicate a Middle Pleistocene or older age. Evidence from pit excavations and geochemistry shows that the underlying saprolites have been truncated by glacial erosion and that previously weathered material has been incorporated into the till sequence. Hence, ¹⁰Be inventories in the tills are dominated by material recycled from Middle Pleistocene or older soils, near-surface sediments and saprolite, and cannot be used to date the periods of till deposition. The retention of relict ¹⁰Be in the tills nonetheless confirms minimal glacial erosion.Concentrations of meteoric ¹⁰Be in the saprolites are lower than any reported saprolite concentrations measured in other settings. Uncertainty in the pre-glaciation ¹⁰Be concentrations in the saprolites makes age determinations difficult. One possibility is that that the saprolite had higher ¹⁰Be concentrations in the past but that saprolite formation ended after glaciation and burial by till and that the ¹⁰Be has substantially decayed. Modelling of the meteoric ¹⁰Be depth profiles in this case suggests that the saprolites in the Parkajoki area were formed at a minimum of 2 Ma. Erosion of the saprolite allows an older age of up to ∼5 Ma, with up to 250 cm of material removed and incorporated into later tills. A second possibility is that concentrations of meteoric ¹⁰Be in the saprolite were originally lower, with formation of the saprolite in a period or periods of ice- and permafrost-free conditions before 0.8 Ma.     

Meteoric 10Be speciation in subglacial sediments of East Antarctica

A sequential chemical extraction procedure was developed and tested to investigate the utility of meteoric ¹⁰Be as a tracer for authigenic mineral formation beneath the East Antarctic Ice Sheet. Subglacial meltwater is widely available under the Antarctic Ice Sheet and dissolved gases within it have the potential to drive chemical weathering processes in the subglacial environment. Meteoric ¹⁰Be is a cosmogenic nuclide with a half-life of 1.39·10⁶ years that is incorporated into glacier ice, therefore its abundance in the subglacial environment in Antarctica is meltwater dependent. It is known to adsorb to fine-grained particles in aqueous solution, precipitate with amorphous oxides/hydroxides, and/or be incorporated into authigenic clay minerals during chemical weathering. The presence of ¹⁰Be in chemical weathering products derived from beneath the ice therefore indicates chemical weathering processes in the subglacial environment. Freshly emerging subglacial sediments from the Mt. Achernar blue ice moraine were subject to chemical extractions where these weathering phases were isolated and ¹⁰Be concentrations therein quantified. Optimization of the phase isolation was developed by examining the effects of each extraction on the sample mineralogy and chemical composition. Experiments on ¹⁰Be desorption revealed that pH 3.2–3.5 was optimal for the extraction of adsorbed ¹⁰Be. Vigorous disaggregation of the samples before grain size separations and acid extractions is crucial due to the incorporation of the nuclide in clay minerals and its preferential absorption to clay-sized particles. ¹⁰Be concentrations of 2–22·10⁷ atoms·g⁻¹ measured in oxides and clay minerals in freshly emerging sediments strongly indicate subglacial chemical weathering in the catchment of the Mt. Achernar moraine. Based on total ¹⁰Be sample concentrations, local basal melt rates, and ¹⁰Be ice concentrations, sediment-meltwater contact in the subglacial environment is on the order of thousands of years per gram of underlying fine sediment. Strong correlation (R = 0.97) between ¹⁰Be and smectite abundance in the sediments supports authigenic clay formation in the subglacial environment. This suggests meteoric ¹⁰Be is a useful tool to characterize subglacial geochemical weathering processes under the Antarctic Ice Sheet.     

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.     

Evidence for glaciation in the Northern Hemisphere back to 44 Ma from ice-rafted debris in the Greenland Sea

The widely accepted age estimate for the onset of glaciation in the Northern Hemisphere ranges between 2 and 15 million years ago (Ma). However, recent studies indicate the date for glacial onset may be significantly older. We report the presence of ice-rafted debris (IRD) in ～ 44 to 30 Ma sediments from the Greenland Sea, evidence for glaciation in the North Atlantic during the Middle Eocene to Early Oligocene. Detailed sedimentological evidence indicates that glaciers extended to sea level in the region, allowing icebergs to be produced. IRD may have been sourced from tidewater glaciers, small ice caps, and/or a continental ice sheet.     

Geochemical analyses of air from an ancient debris-covered glacier,Antarctica

We examined air trapped in ancient ice from three shallow cores (<35 m deep) recovered from stagnant portions of the Mullins glacier, an 8 km long debris-covered alpine glacier in the McMurdo Dry Valleys that is overlain by several in-situ volcanic ash-fall deposits. Previously reported ⁴⁰Ar/³⁹Ar dates on ash-fall in the vicinity of the core sites average 4.0 Ma, and underlying ice is presumably as old in some areas. We analyzed the elemental and isotopic composition of O₂, N₂, and Ar and total air content of the glacial ice. We also dated the trapped air directly to an uncertainty of ±220 kyr (1σ) by measuring its ⁴⁰Ar/³⁶Ar and ³⁸Ar/³⁶Ar ratios. Our results suggest that the air analyzed is likely a mixture of ancient atmosphere trapped at the time of ice formation and more recent air introduced via cracks in the ice that penetrate to at least 33 m. The isotopic signatures of gases have been complicated by gas loss, as well as a mixture of thermal and gravitational fractionation. The oldest age estimated for the trapped air dates to 1.6 Ma, indicating that the original air is at least as old as 1.6 ± 0.2 Ma. A convergence to older ice ages with increasing depth in the deepest core analyzed (33 m) hints at the possibility that pristine air might be recovered at greater depths. Minor interstitial debris present in the glacial ice (<1%), along with geochemical evidence for in-situ microbial respiration, prohibit direct analysis of CO₂. We measured the triple isotopic composition of O₂ as a proxy for CO₂ and infer that, in the air represented in our ice samples, CO₂ concentrations are within the range observed over the last 800 ka.     

Chronology of Early to Mid-Pleistocene sediments in the northern North Sea: New evidence from amino acid and strontium isotope analyses

Sediments deposited during glacial-interglacial cycles through the Early to Mid-Pleistocene in the North Sea are chronologically poorly constrained. To contribute to the chronology of these units, amino acid racemization (AAR) and strontium (Sr) isotope analyses have been performed on samples from four shallow borings and one oil well along a transect in the northern North Sea. D/L Asp (aspartic acid) values obtained through reverse-phase liquid chromatography in the benthic foraminiferal species Elphidium excavatum is focused on because of consistent results and a good stratigraphic distribution of this benthic species. For the Early Pleistocene, an age model for the well 16/1–8, from the central part of the northern North Sea based on Sr ages allows for dating of the prograding wedges filling the pre-Quaternary central basin. A regional calibration curve for the racemization of Asp in Elphidium excavatum is developed using published ages of radiocarbon-dated samples and samples associated with the previously identified Bruhnes/Matuyama (B/M) paleomagnetic boundary and a Sr age from this study. Based on all the available geochronological evidence, samples were assigned to marine oxygen isotope stages (MIS) with uncertainties on the order of 10–70 ka.Sr ages suggest a hiatus of <2 million years (Ma) possibly due to non-deposition or low sedimentation between the Utsira Formation (Pliocene) and the Early Pleistocene. An increase in sedimentation rates around 1.5 ± 0.07 Ma (∼MIS 51) may partly be due to sediment supply from rivers from the south-east and partly due to the extension of ice sheet around 1.36 ± 0.07 Ma from the Norwegian coast to the central North Sea. A possible basin-wide glaciation occurred around 1.1 Ma (∼ MIS 32) (upper regional unconformity/top of unit Q4 in this study), resulting in erosion and regional unconformity. Two interglacials in the Norwegian Channel have been dated: the Radøy Interglacial to 1.07 ± 0.01 Ma (possibly MIS 31, the ‘super interglacial’), and the Norwegian Trench Interglacial to 0.50 ± 0.02 Ma (possibly MIS 13). A massive till unit identified at the same stratigraphic level in all shallow borings may partly represent an extensive MIS 12 glaciation. This study shows that the combined use of amino acid racemization data and Sr isotope chronology can refine the chronological ambiguities of Quaternary North Sea sediments related partly to the impact of glacial processes.     

Chronology of Lateglacial ice flow reorganization and deglaciation in the Gotthard Pass area,Central Swiss Alps,based on cosmogenic 10Be and in situ 14C

We reconstruct the timing of ice flow reconfiguration and deglaciation of the Central Alpine Gotthard Pass, Switzerland, using cosmogenic ¹⁰Be and in situ ¹⁴C surface exposure dating. Combined with mapping of glacial erosional markers, exposure ages of bedrock surfaces reveal progressive glacier downwasting from the maximum LGM ice volume and a gradual reorganization of the paleoflow pattern with a southward migration of the ice divide. Exposure ages of ∼16–14 ka (snow corrected) give evidence for continuous early Lateglacial ice cover and indicate that the first deglaciation was contemporaneous with the decay of the large Gschnitz glacier system. In agreement with published ages from other Alpine passes, these data support the concept of large transection glaciers that persisted in the high Alps after the breakdown of the LGM ice masses in the foreland and possibly decayed as late as the onset of the Bølling warming. A younger group of ages around ∼12–13 ka records the timing of deglaciation following local glacier readvance during the Egesen stadial. Glacial erosional features and the distribution of exposure ages consistently imply that Egesen glaciers were of comparatively small volume and were following a topographically controlled paleoflow pattern. Dating of a boulder close to the pass elevation gives a minimum age of 11.1 ± 0.4 ka for final deglaciation by the end of the Younger Dryas. In situ ¹⁴C data are overall in good agreement with the ¹⁰Be ages and confirm continuous exposure throughout the Holocene. However, in situ ¹⁴C demonstrates that partial surface shielding, e.g. by snow, has to be incorporated in the exposure age calculations and the model of deglaciation.     

10Be surface-exposure age dating of the Last Glacial Maximum in the northern Pamir (Tajikistan)

Knowledge of the spatial and temporal variations in Alpine glaciations is essential for reconstructing the regional and global timing of ice ages. This study investigates glacial deposits at the mouth of the Muksu catchment in the northern Pamir using ¹⁰Be surface-exposure age dating. We sampled boulders from the furthest downstream recessional moraine (20 samples) and five lateral moraines (41 samples) near the former terminus of the Fedchenko Glacier, the longest (∼72 km) present-day Alpine glacier of the Pamir. After the identification of outliers, the boulder population of the recessional moraine yielded a mean exposure age of 17.5 ± 1.9 ka. The maximum exposure age of the lateral moraines, collected ∼5 km up-valley of the recessional moraine, is 18.2 ± 1.7 ka. The boulder ages reflect glacial deposition during the Last Glacial Maximum (Marine Isotope Stage 2) in the region; they are in accordance with published glacial deposition ages in the western Tian Shan.     

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.     

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.     

P–PINI: A cosmogenic nuclide burial dating method for landscapes undergoing non-steady erosion

Existing methods of cosmogenic nuclide burial dating perform well provided that sediment sources undergo steady rates of erosion and the samples experience continuous exposure to cosmic rays. These premises exert important limitations on the applicability of the methods. And yet, high mountain sediment sources are rife with transient processes, such as non-steady erosion by glacial quarrying and/or landsliding, or temporary cosmic-ray shielding beneath glaciers and/or sediment. As well as breaching the premises of existing burial dating methods, such processes yield samples with low nuclide abundances and variable ²⁶Al/¹⁰Be ratios that may foil both isochron and simple burial-age solutions. P–PINI (Particle-Pathway Inversion of Nuclide Inventories) is a new dating tool designed for dating the burial of sediments sourced from landscapes characterized by abrupt, non-steady erosion, discontinuous exposure, and catchments with elevation-dependent ²⁶Al/¹⁰Be production ratios. P–PINI merges a Monte Carlo simulator with established cosmogenic nuclide production equations to simulate millions of samples (¹⁰Be–²⁶Al inventories). The simulated samples are compared statistically with ¹⁰Be–²⁶Al measured in field samples to define the most probable burial age. Here, we target three published ¹⁰Be–²⁶Al datasets to demonstrate the versatility of the P–PINI model for dating fluvial and glacial sediments. (1) The first case serves as a robust validation of P–PINI. For the Pulu fluvial gravels (China), we obtain a burial age of 1.27 ± 0.10 Ma (1σ), which accords with the isochron burial age and two independent chronometers reported in Zhao et al. (2016) Quaternary Geochronology 34, 75–80. The second and third cases, however, reveal marked divergence between P–PINI and isochron-derived ages. (2) For the fluvial Nenana Gravel (USA), we obtain a minimum-limiting burial age of 4.5 ± 0.7 Ma (1σ), which is compatible with unroofing of the Alaska Range starting ∼ 6 Ma, while calling into question the Early Pleistocene isochron burial age presented in Sortor et al. (2021) Geology 49, 1473–1477. (3) For the Bünten Till (Switzerland), we obtain a limiting burial age of <204 ka (95th percentile range), which conforms with the classical notion of the most extensive glaciation in the northern Alpine Foreland assigned to the Riss glaciation (sensu marine isotope stage 6) contrary to the isochron burial age presented in Dieleman et al. (2022) Geosciences, 12, 39. Discrepancies between P–PINI and the isochron ages are rooted in the challenges posed by the diverse pre-burial ²⁶Al/¹⁰Be ratios produced under conditions characteristic of high mountain landscapes; i.e., non-steady erosion, discontinuous cosmic-ray exposure, and elevation-dependent ²⁶Al/¹⁰Be production ratios in the source region, which are incompatible with the isochron method, but easily accommodated by the stochastic design of P–PINI.     

Formation,mechanism and significance of alluvial-dammed lakes in Golmud River catchment,north-eastern Qinghai-Tibetan Plateau

This study investigated a series of dammed lakes and downstream-adjacent alluvial fans in the upstream to middle reaches of the Golmud River in the eastern Kunlun Mountain, on the north-eastern Qinghai-Tibetan Plateau (QTP). An optically stimulated luminescence (OSL) chronology shows the sediments of five dammed lakes developed from c. 45–40, 30–25, 18–14, and 12–8 ka, corresponding to MIS 3b, late MIS 3a, Last Deglaciation, and early Holocene, respectively. The remote sensing data show these dammed lakes have a total area of 109.4 km², with the lake volume of more than 4.0 km³. Symmetric alluvial fans from north–south tributary valleys produced OSL ages of c. 61–52, 42–31, 26–20, and 16–10 ka, corresponding to glaciation periods: the MIS 3c and MIS 3a, MIS 2, and the Last Deglaciation. This suggests that glacial activity is responsible for the alluvial fan development, where dammed rivers occurred first, but lake formation did not take place synchronously until later periods of strong hydrologic activity, resulting from northward intrusions of the Indian summer monsoon (ISM) or glacier melt. Thus, the blocking pattern is that river valleys were dammed during periods of glacial activity and lakes formed during wet periods. The lake formation and subsequent drainage may have resulted in: (i) impeded headwater incision and strengthening of downstream dissection; (ii) enriched the halite and potash in the distal Qarhan Salt Lake through hydrologic and hydrochemical processes of abundant water input, the salt lake expansion, salt redissolution from playa and final resedimentation during later dry periods. The alluvial-dammed lake pattern in the mountain-basin systems of eastern Kunlun Mountain offers a model for assessing the linkages between monsoon dynamics, geomorphic processes and distal salt lake evolutions in other arid regions.     

Ice cap erosion patterns from bedrock 10Be and 26Al,southeastern Tibetan Plateau

Quantifying glacial erosion contributes to our understanding of landscape evolution and topographic relief production in high altitude and high latitude areas. Combining in situ ¹⁰Be and ²⁶Al analysis of bedrock, boulder, and river sand samples, geomorphological mapping, and field investigations, we examine glacial erosion patterns of former ice caps in the Shaluli Shan of the southeastern Tibetan Plateau. The general landform pattern shows a zonal pattern of landscape modification produced by ice caps of up to 4000 km² during pre-LGM (Last Glacial Maximum) glaciations, while the dating results and landforms on the plateau surface imply that the LGM ice cap further modified the scoured terrain into different zones. Modeled glacial erosion depth of 0–0.38 m per 100 ka bedrock sample located close to the western margin of the LGM ice cap, indicates limited erosion prior to LGM and Late Glacial moraine deposition. A strong erosion zone exists proximal to the LGM ice cap marginal zone, indicated by modeled glacial erosion depth >2.23 m per 100 ka from bedrock samples. Modeled glacial erosion depths of 0–1.77 m per 100 ka from samples collected along the edge of a central upland, confirm the presence of a zone of intermediate erosion in-between the central upland and the strong erosion zone. Significant nuclide inheritance in river sand samples from basins on the scoured plateau surface also indicate restricted glacial erosion during the last glaciation. Our study, for the first time, shows clear evidence for preservation of glacial landforms formed during previous glaciations under non-erosive ice on the Tibetan Plateau. As patterns of glacial erosion intensity are largely driven by the basal thermal regime, our results confirm earlier inferences from geomorphology for a concentric basal thermal pattern for the Haizishan ice cap during the LGM. © 2018 John Wiley & Sons, Ltd.     

A locally calibrated,late glacial 10Be production rate from a low-latitude,high-altitude site in the Peruvian Andes

Well-dated records of tropical glacier fluctuations are essential for developing hypotheses and testing proposed mechanisms for past climate changes. Since organic material for radiocarbon dating is typically scarce in low-latitude, high-altitude environments, surface exposure-age dating, based on the measurement of in situ produced cosmogenic nuclides, provides much of the chronologic information on tropical glacier moraines. Here, we present a locally calibrated ¹⁰Be production rate for a low-latitude, high-altitude site near Quelccaya Ice Cap (∼13.95°S, 70.89°W, 4857 m asl) in the southeastern Peruvian Andes. Using an independent age (12.35 +0.2, −0.02 ka) of the late glacial Huancané IIa moraines based on thirty-four bracketing radiocarbon ages and twelve ¹⁰Be concentrations of boulders on the moraines, we determine a local production rate of 43.28 ± 2.69 atoms gram⁻¹ year⁻¹ (at g⁻¹ yr⁻¹). Reference ¹⁰Be production rates (i.e., production rates by neutron spallation appropriate for sea-level, high-latitude sites) range from 3.97 ± 0.09 to 3.78 ± 0.09 at g⁻¹ yr⁻¹, determined using scaling after Lal (1991) and Stone (2000) and depending on our assumed boulder surface erosion rate. Since our boulder surface erosion rate estimate is a minimum value, these reference production rates are also minimum values. A secondary control site on the Huancané IIIb moraines suggests that the ¹⁰Be production rates are at least as low as, or possibly lower than, those derived from the Huancané IIa moraines. These sea-level, high-latitude production rates are at least 11–15% lower than values derived using the traditional global calibration dataset, and they are also lower than those derived from the late glacial Breque moraine in the Cordillera Blanca of Peru. However, our sea-level, high-latitude production rates agree well with recently published, locally calibrated production rates from the Arctic, New Zealand, and Patagonia. The production rates presented here should be used to calculate ¹⁰Be exposure ages in low-latitude, high-altitude locations, particularly in the tropical Andes, and should improve the ability to compare the results of studies using ¹⁰Be exposure-age dating with other chronological data.     

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.     

Large-volume silicic volcanism in Kamchatka: Ar–Ar and U–Pb ages,isotopic, and geochemical characteristics of major pre-Holocene caldera-forming eruptions

The Kamchatka Peninsula in far eastern Russia represents the most volcanically active arc in the world in terms of magma production and the number of explosive eruptions. We investigate large-scale silicic volcanism in the past several million years and present new geochronologic results from major ignimbrite sheets exposed in Kamchatka. These ignimbrites are found in the vicinity of morphologically-preserved rims of partially eroded source calderas with diameters from ～ 2 to ～ 30 km and with estimated volumes of eruptions ranging from 10 to several hundred cubic kilometers of magma. We also identify and date two of the largest ignimbrites: Golygin Ignimbrite in southern Kamchatka (0.45 Ma), and Karymshina River Ignimbrites (1.78 Ma) in south-central Kamchatka. We present whole-rock geochemical analyses that can be used to correlate ignimbrites laterally. These large-volume ignimbrites sample a significant proportion of remelted Kamchatkan crust as constrained by the oxygen isotopes. Oxygen isotope analyses of minerals and matrix span a 3‰ range with a significant proportion of moderately low-δ¹⁸O values. This suggests that the source for these ignimbrites involved a hydrothermally-altered shallow crust, while participation of the Cretaceous siliceous basement is also evidenced by moderately elevated δ¹⁸O and Sr isotopes and xenocryst contamination in two volcanoes. The majority of dates obtained for caldera-forming eruptions coincide with glacial stages in accordance with the sediment record in the NW Pacific, suggesting an increase in explosive volcanic activity since the onset of the last glaciation 2.6 Ma. Rapid changes in ice volume during glacial times and the resulting fluctuation of glacial loading/unloading could have caused volatile saturation in shallow magma chambers and, in combination with availability of low-δ¹⁸O glacial meltwaters, increased the proportion of explosive vs effusive eruptions. The presented results provide new constraints on Pliocene–Pleistocene volcanic activity in Kamchatka, and thus constrain an important component of the Pacific Ring of Fire.     

Streamflow variation due to glacier melting and climate change in upstream Heihe River Basin,Northwest China

Streamflow simulation is often challenging in mountainous watersheds because of incomplete hydrological models, irregular topography, immeasurable snowpack or glacier, and low data resolution. In this study, a semi-distributed conceptual hydrological model (SWAT-Soil Water Assessment Tool) coupled with a glacier melting algorithm was applied to investigate the sensitivity of streamflow to climatic and glacial changes in the upstream Heihe River Basin. The glacier mass balance was calculated at daily time-step using a distributed temperature-index melting and accumulation algorithm embedded in the SWAT model. Specifically, the model was calibrated and validated using daily streamflow data measured at Yingluoxia Hydrological Station and decadal ice volume changes derived from survey maps and remote sensing images between 1960 and 2010. This study highlights the effects of glacier melting on streamflow and their future changes in the mountainous watersheds. We simulate the contribution of glacier melting to streamflow change under different scenarios of climate changes in terms of temperature and precipitation dynamics. The rising temperature positively contributed to streamflow due to the increase of snowmelt and glacier melting. The rising precipitation directly contributes to streamflow and it contributed more to streamflow than the rising temperature. The results show that glacial meltwater has contributed about 3.25 billion m³ to streamflow during 1960–2010. However, the depth of runoff within the watershed increased by about 2.3 mm due to the release of water from glacial storage to supply the intensified evapotranspiration and infiltration. The simulation results indicate that the glacier made about 8.9% contribution to streamflow in 2010. The research approach used in this study is feasible to estimate the glacial contribution to streamflow in other similar mountainous watersheds elsewhere.     

Rates of rockwall slope erosion in the upper Bhagirathi catchment,Garhwal Himalaya

Rockwall slope erosion is defined for the upper Bhagirathi catchment using cosmogenic Beryllium-10 (¹⁰Be) concentrations in sediment from medial moraines on Gangotri glacier. Beryllium-10 concentrations range from 1.1 ± 0.2 to 2.7 ± 0.3 × 10⁴ at/g SiO₂, yielding rockwall slope erosion rates from 2.4 ± 0.4 to 6.9 ± 1.9 mm/a. Slope erosion rates are likely to have varied over space and time and responded to shifts in climate, geomorphic and/or tectonic regime throughout the late Quaternary. Geomorphic and sedimentological analyses confirm that the moraines are predominately composed of rockfall and avalanche debris mobilized from steep relief rockwall slopes via periglacial weathering processes. The glacial rockwall slope erosion affects sediment flux and storage of snow and ice at the catchment head on diurnal to millennial timescales, and more broadly influences catchment configuration and relief, glacier dynamics and microclimates. The slope erosion rates exceed the averaged catchment-wide and exhumation rates of Bhagirathi and the Garhwal region on geomorphic timescales (10³−10⁵ years), supporting the view that erosion at the headwaters can outpace the wider catchment. The ¹⁰Be concentrations of medial moraine sediment for the upper Bhagirathi catchment and the catchments of Chhota Shigri in Lahul, northern India and Baltoro glacier in Central Karakoram, Pakistan show a tentative relationship between ¹⁰Be concentration and precipitation. As such there is more rapid glacial rockwall slope erosion in the monsoon-influenced Lesser and Greater Himalaya compared to the semi-arid interior of the orogen. Rockwall slope erosion in the three study areas, and more broadly across the northwest Himalaya is likely governed by individual catchment dynamics that vary across space and time. © 2019 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons, Ltd.     

Quantifying glacial erosion on a limestone bed and the relevance for landscape development in the Alps

Glacial erosion is the basic process that has shaped the landscapes of the Alps. Despite intense research over centuries, and the use of various techniques, determination of glacial erosion rates remains challenging. This is not only because the location where the process occurs is almost inaccessible, but also because it is dependent on many different factors, including ice thickness and velocity, glacier thermal regime and lithology. Reported glacial erosion rates range over several orders of magnitude (0.01 to >10 mm a⁻¹). Most studies focus on crystalline bedrock, whereas few researchers have investigated glacial erosion on limestone. Here we analyse glacially polished bedrock surfaces at the recently deglaciated forefield of the Tsanfleuron glacier, Swiss Alps. The nearly horizontally bedded limestone hosts a well-developed karst system. Meltwater from the glacier drains into the subsurface within a few metres of the ice margin. By combining geomorphological mapping, measurement of cosmogenic ³⁶Cl concentrations of glacially eroded bedrock surfaces and a numerical model (MECED), we quantify at each sample location the amount of rock removed during glacier occupation. The glacial erosion rates calculated from these values range from 0 to 0.08 mm a⁻¹. These are orders of magnitude lower than values measured at comparable sites on crystalline bedrock. The high ³⁶Cl concentrations we measured show that the Tsanfleuron glacier was unable to effectively erode the gently dipping, strongly karstified limestone. We suggest that this effect may play a key role in formation and preservation over many glacial cycles of high-elevation, low-relief limestone plateaus in the Alps. © 2020 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