Lateral plucking as a mechanism for elongate erosional glacial bedforms: explaining megagrooves in Britain and Canada

Megagrooves are kilometre‐scale linear topographic lows carved in bedrock, separated by ridges, typically in areas of largely devoid of till. They have been reported from several areas covered by Pleistocene glaciations, such as Canadian Northwest (NW) Territories, Michigan and NW Scotland. Here we report two previously undocumented megagroove fields from Ungava, Canada, and northern England, and present new analyses of the megagrooves from NW Scotland. This paper seeks to determine the nature of the lithological and structural controls on the occurrence and formation of megagrooves. Analysis of both geomorphological and bedrock properties shows that megagrooves are generally:

1. confined to well stratified or layered bedrock, such as (meta)sedimentary rocks with closely spaced joints, and tend not to occur on massive rocks such as gneiss or granite, or thick‐bedded sedimentary rocks;
2. subparallel to palaeo‐ice flow and the strike of the strata; and tend not to occur where palaeo‐ice flow is at high angles to the strike of strata;
3. produced by significant glacial erosion by sustained unidirectional ice flow.

Detailed analysis of megagrooves in NW Scotland shows that neither glacio‐fluvial erosion, nor differential abrasion was the dominant mechanism of formation. A mechanism, here termed ‘lateral plucking’, is suggested that involves block plucking on rock steps parallel to ice flow. Removal of joint‐bounded blocks from such rock steps involves a component of rotation along a vertical axis. Block removal may be enhanced by a direct component of shear stress onto the vertical stoss sides. The lateral plucking mechanism results in horizontal erosion at right angles to the ice flow, and enhances the groove/ridge topography. Megagrooves are potentially useful as palaeo‐ice flow indicators in areas devoid of till, and can thus complement the palaeo‐ice stream datasets which are presently largely based on soft‐sediment landform studies. British Geological Survey © NERC 2011     

Retreating ice: research in pro‐glacial areas matters

Glacier forefields are landscapes in transition from glacial to non‐glacial conditions; this implies intense geomorphic, hydrological and ecological dynamics with important on‐ and off‐site effects. This special issue collects 13 papers covering recent research in both (sub‐)polar and alpine pro‐glacial environments that focus on (i) pro‐glacial sediment sources, (ii) pro‐glacial rivers, (iii) pro‐glacial lakes, (iv) ground water and ice, and (v) the development of soil and vegetation in its interplay with morphodynamics. Advances in mapping, surveying and geophysical techniques form the basis for research perspectives related to the historical evolution of pro‐glacial areas, the understanding of complex interactions of multiple processes, and the effects of continued glacier recession. Copyright © 2015 John Wiley & Sons, Ltd.     

Size,shape and spatial arrangement of mega‐scale glacial lineations from a large and diverse dataset

Mega‐scale glacial lineations (MSGLs) are a characteristic landform on ice stream beds. Solving the puzzle of their formation is key to understanding how ice interacts with its bed and how this, in turn, influences the dynamics of ice streams. However, a comprehensive and detailed characterization of this landform's size, shape and spatial arrangement, which might serve to test and refine formational theories, is largely lacking. This paper presents a detailed morphometric analysis and comparison of 4043 MSGLs from eight palaeo‐ice stream settings: three offshore (Norway and Antarctica), four onshore (Canada), and one from under a modern ice stream in West Antarctica. The length of MSGLs is lower than previously suggested (mode 1000–2000 m; median 2892 m), and they initiate and terminate at various locations on an ice stream bed. Their spatial arrangement reveals a pattern that is characterized by an exceptional parallel conformity (80% of all mapped MSGLs have an azimuth within 5° from the mean values), and a fairly constant lateral spacing (mode 200–300 m; median 330 m), which we interpret as an indication that MSGLs are a spatially self‐organized phenomenon. Results show that size, shape and spatial arrangement of MSGLs are consistent both within and also generally between different ice stream beds. We suggest this results from a common mechanism of formation, which is largely insensitive to local factors. Although the elongation of MSGLs (mode 6–8; median 12.2) is typically higher than features described as drumlins, these values and those of their width (mode 100–200 m; median 268 m) overlap, which suggests the two landforms are part of a morphological continuum and may share a similar origin. We compare their morphometry with explicit predictions made by the groove‐ploughing and rilling instability theories of MSGL formation. Although the latter was most compatible, neither is fully supported by observations. © 2014 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.     

Spatio‐temporal analysis of glacial ice area distribution of Hunza River Basin,Karakoram region of Pakistan

There is a high degree of uncertainty about the state and fate of Pakistan's Karakoram glaciers due to data scarcity in high altitude regions. They are thought to be less vulnerable to climatic change because they behave differently as compared with eastern Himalayas. This study measures the decadal temporal changes in the glacial ice area of Karakoram's Hunza River Basin, one of the eight subbasins of Upper Indus Basin. An attempt has been made to investigate the relationship between glacial ice area changes and calculated values of precipitation, temperature and run‐off. A combination of satellite and field‐based approach is applied. Output includes maps of glacial ice hypsometries of eight glacial ice subregions of Hunza River Basin for 3 years (i.e., 1989, 2002, and 2010). The results show a decreasing trend in the glacial ice‐covered area signifying a reduction of 20.47% with the largest reduction being in the lower elevation bands. There is presently no conclusive answer as to why glacial ice in the Karakoram is acting differently from the near‐global indication of glacial ice changes. Climate data from high altitudes are needed to find answer for this anomalous behaviour.     

In situ cosmogenic 10Be dating of the Quaternary glaciations in the southern Shaluli Mountain on the Southeastern Tibetan Plateau

It is generally considered that four-times ice age happened during the Quaternary epoch on the Tibetan Plateau. However, the research on the chronology of the four-times ice age is far from enough. The Shaluli Mountain on the Southeastern Tibetan Plateau is an ideal place for plaeo-glacier study, because there are abundant Quaternary glacial remains there. This paper discusses the ages of the Quaternary glaciations, based on the exposure dating of roche moutonnée, moraines and glacial erosion surfaces using in situ cosmogenic isotopes ¹⁰Be. It is found that the exposure age of the roche moutonnée at Tuershan is 15 ka, corresponding to Stage 2 of the deep-sea oxygen isotope, suggesting that the roche moutonnée at Tuershan is formed in the last glacial maximum. The exposure age of glacial erosion surface at Laolinkou is 130–160 ka, corresponding to Stage 6 of the deep-sea oxygen isotope. The oldest end moraine at Kuzhaori may form at 421–766 kaBP, corresponding to Stages 12–18 of the deep-sea oxygen isotope. In accordance with the climate characteristic of stages 12, 14, 16 and 18 reflected by the deep-sea oxygen isotope, polar ice cores and loess sequence, the oldest end moraine at Kuzhaori may form at stage 12 or stage 16, the latter is more possible.     

Digital landscapes of deglaciation: identifying Late Quaternary glacial lake outburst floods using LiDAR

High resolution DEMs obtained from LiDAR topographic data have led to improved landform inventories (e.g. landslides and fault scarps) and understanding of geomorphic event frequency. Here we use airborne LiDAR mapping to investigate meltwater pathways associated with the Tweed Valley palaeo ice‐stream (UK). In particular we focus on a gorge downstream of Palaeolake Milfield, previously mapped as a sub‐glacial meltwater channel, where the identification of abandoned headcut channels, run‐up bars, rock‐cut terrace surfaces and eddy flow features attest to formation by a sub‐aerial glacial lake outburst flood (GLOF) caused by breaching of a sediment dam, likely an esker ridge. Mapping of these landforms combined with analysis of the gorge rim elevations and cross‐section variability revealed a two phase event with another breach site downstream following flow blockage by higher elevation drumlin topography. We estimate the magnitude of peak flow to be 1–3 × 10³ m³/s, duration of the event to range from 16–155 days, and a specific sediment yield of 10⁷–10⁹ m³/km²/yr. We identified other outburst pathways in the lower Tweed basin that help delineate an ice margin position of the retreating Tweed Valley ice stream. The results suggest that low magnitude outburst floods are under‐represented in Quaternary geomorphological maps. We therefore recommend regional LiDAR mapping of meltwater pathways to identify other GLOFs in order to better quantify the pattern of freshwater and sediment fluxes from melting ice sheets to oceans. Despite the relatively low magnitude of the Till outburst event, it had a significant impact on the landscape development of the lower Tweed Valley through the creation of a new tributary pathway and triggering of rapid knickpoint retreat encouraging new regional models of post‐glacial fluvial landscape response. Copyright © 2015 John Wiley & Sons, Ltd.     

High resolution mapping of supra‐glacial drainage pathways reveals link between micro‐channel drainage density,surface roughness and surface reflectance

This paper reports on the use of a small unmanned aerial vehicle (sUAV) carrying a standard compact camera, to construct a high resolution orthomosaic (OM) and digital elevation model (DEM) over the lower reaches of the glacier Midtre Lovénbreen, Svalbard. Structure from Motion (SfM) techniques were used to build the OM and DEM, and together these reveal insights into the nature of supra‐glacial drainage. Major meandering supra‐glacial drainage pathways show clear dynamism, via meander cutoffs and abandoned channels. In addition, the imagery reveals a very extensive network of smaller channels that may well carry substantial amounts of water. This network of channels is in part controlled by the structure of the glacier, but in turn, these channels have a significant impact on the ice surface. Roughness of the ice surface is higher where channels are most extensive. In addition, we find a relationship between channel density and surface reflectance, such that greater channel density is associated with lower reflectance values. Given the role of surface reflectance and roughness in the energy balance of glaciers, it is therefore apparent that extensive networks of small supra‐glacial channels across such glaciers have the potential to have an important impact on energy exchanges between the atmosphere and the ice surface. Copyright © 2015 John Wiley & Sons, Ltd.     

CLustre: semi‐automated lineament clustering for palaeo‐glacial reconstruction

Datasets containing large numbers (>10 000) of glacial lineaments are increasingly being mapped from remotely sensed data in order to develop a palaeo‐glacial reconstruction or ‘inversion’. The palimpsest landscape presents a complex record of past ice flow and deconstructing this information into a logical history is an involved task. One stage in this process requires the identification of sets of genetically linked lineaments that can form the basis of a reconstruction. This paper presents a semi‐automated algorithm, CLustre, for lineament clustering that uses a locally adaptive, region growing, methodology. After outlining the algorithm, it is tested on synthetic datasets that simulate parallel and orthogonal cross‐cutting lineaments, encompassing 1500 separate classifications. Results show robust classification in most scenarios, although parallel overlap of lineaments can cause false positive classification unless there are differences in lineament length. Case studies for Dubawnt Lake and Victoria Island, Canada, are presented and compared with existing datasets. For Dubawnt Lake 9 out of 14 classifications directly match incorporating 89% of lineaments. For Victoria Island 57 out of 58 classifications directly match incorporating 95% of lineaments. Differences are related to small numbers of unclassified lineaments and parallel cross‐cutting lineaments that are of a similar length. CLustre enables the automated, repeatable, assignment of lineaments to flow sets using defined user criteria. This is important as qualitative visual interpretation may introduce bias, potentially weakening the testability of palaeo‐glacial reconstructions. In addition, once classified, summary statistics of lineament clusters can be calculated and subsequently used during the reconstruction process. Copyright © 2015 John Wiley & Sons, Ltd.     

Geomorphic evolution of small river–lake‐systems in northeast Germany during the Late Quaternary

This study investigates the post‐glacial development of four small river–lake systems in the Weichselian belt of northern central Europe. The valleys investigated are part of an immature drainage system characterized by frequent and abrupt changes in flow direction and the presence of numerous stagnant‐ice depressions in the valley course. The depressions contain thick sedimentary sequences which provide excellent archives for the reconstruction of the post‐glacial valley development. Study results indicate that the valleys reuse segments of former subglacial meltwater channels. During the Late Pleniglacial these channels carried meltwater streams. Stagnant‐ice melting occurred in stages from the Oldest Dryas to the early Holocene and was often followed by the formation of lakes in the valley course. Flow reversals occurred during the Late‐glacial–Holocene transition and were in response to general base‐level lowering caused by stagnant‐ice melting, headwater erosion and lake overspills. Lacustrine deposition typically started during the early Late‐glacial comprising mainly silicate gyttjas, whereas organic gyttjas and peats accumulated during the Allerød. The Younger Dryas is associated with a marked increase in fluvial and aeolian sedimentation, and lake‐level high stands. This was followed by early Holocene lake‐level low stands and a subsequent stabilization phase with decreasing silicate input and increasing organic lacustrine deposition. In general, dramatic changes in Late Pleniglacial to early Holocene sedimentation suggest that small‐scale catastrophic events played a more important role in triggering geomorphic changes then previously recognized. Infilling continued until peat accumulation and terrestrialization of lake basins became widespread during the mid‐ to late Holocene. Beginning in the late Holocene anthropogenic influences become important mainly involving an increase in sediment supply due to forest clearing and land use, followed by mill stowage, river course correction and anthropogenic lake‐level manipulations. Copyright © 2007 John Wiley & Sons, Ltd.     

A large landslide event in a post‐glacial landscape: rethinking glacial legacy

Threlkeld Knotts (c. 500 m above sea level) in the English Lake District has hitherto been considered to be a glacially‐modified intrusion of microgranite. However, its surface features are incompatible with glacial modification; neither can these nor the subsurface structures revealed by ground‐penetrating radar (GPR) be explained by post‐glacial subaerial processes acting on a glacially‐modified microgranite intrusion. Here we re‐interpret Threlkeld Knotts as a very large post‐glacial landslide involving the microgranite, with an estimated volume of about 4 × 10⁷ m³. This interpretation is tested against published and recent information on the geology of the site, the glacial geomorphic history of the area and newly‐acquired GPR data. More than 60 large post‐Last Glacial Maximum (LGM) rock–slope failures have significantly modified the glaciated landscape of the Lake District; this is one of the largest. Recognition of this major landslide deposit in such a well‐studied environment highlights the need to continuously re‐examine landscapes in the light of increasing knowledge of geomorphic processes and with available technology in currently active or de‐glaciating environments. Copyright © 2012 John Wiley & Sons, Ltd.     

Mega-scale glacial lineations and cross-cutting ice-flow landforms

Landsat images reveal a previously unsuspected large-scale pattern of streamlining within drift that is assumed to reflect former phases of ice flow. Such a glacial grain can be regarded as a landform assemblage comprised of a number of components. Drumlins and megaflutes form part of the pattern, but in addition there are two previously undocumented ice-moulded landform elements: streamlined lineations of much greater proportions, referred to as mega-scale glacial lineations, and a distinctive cross-cutting topology within the grain. The ice-moulded landform assemblage is described and illustrated with reference to examples from Canada. Possible modes of genesis of such landforms are discussed and their glaciological implications outlined. The discovery of this pattern indicates the pervasive nature of subglacial deformation of sediment, and demands a radical re-interpretation of ice sheet dynamics.     

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.     

Pressure‐driven,shoreline currents in a perennially ice‐covered,pro‐glacial lake in Antarctica,identified from a LiCl tracer injected into a pro‐glacial stream

The distribution of streamwater within ice‐covered lakes influences sub‐ice currents, biological activity and shoreline morphology. Perennially ice‐covered lakes in the McMurdo Dry Valleys, Antarctica, provide an excellent natural laboratory to study hydrologic–limnologic interactions under ice cover. For a 2 h period on 17 December 2012, we injected a lithium chloride tracer into Andersen Creek, a pro‐glacial stream flowing into Lake Hoare. Over 4 h, we collected 182 water samples from five stream sites and 15 ice boreholes. Geochemical data showed that interflow travelled West of the stream mouth along the shoreline and did not flow towards the lake interior. The chemistry of water from Andersen Creek was similar to the chemistry of water below shoreline ice. Additional evidence for Westward flow included the morphology of channels on the ice surface, the orientation of ripple marks in lake sediments at the stream mouth and equivalent temperatures between Andersen Creek and water below shoreline ice. Streamwater deflected to the right of the mouth of the stream, in the opposite direction predicted by the Coriolis force. Deflection of interflow was probably caused by the diurnal addition of glacial runoff and stream discharge to the Eastern edge of the lake, which created a strong pressure gradient sloping to the West. This flow directed stream momentum away from the lake interior, minimizing the impact of stream momentum on sub‐ice currents. It also transported dissolved nutrients and suspended sediments to the shoreline region instead of the lake interior, potentially affecting biological productivity and bedform development. Copyright © 2014 John Wiley & Sons, Ltd.     

Geomorphology under ice streams: Moving from form to process

Ice streams are integral components of an ice sheet's mass balance and directly impact on sea level. Their flow is governed by processes at the ice‐bed interface which create landforms that, in turn, modulate ice stream dynamics through their influence on bed topography and basal shear stresses. Thus, ice stream geomorphology is critical to understanding and modelling ice streams and ice sheet dynamics. This paper reviews developments in our understanding of ice stream geomorphology from a historical perspective, with a focus on the extent to which studies of modern and palaeo‐ice streams have converged to take us from a position of near‐complete ignorance to a detailed understanding of their bed morphology. During the 1970s and 1980s, our knowledge was limited and largely gleaned from geophysical investigations of modern ice stream beds in Antarctica. Very few palaeo‐ice streams had been identified with any confidence. During the 1990s, however, glacial geomorphologists began to recognise their distinctive geomorphology, which included distinct patterns of highly elongated mega‐scale glacial lineations, ice stream shear margin moraines, and major sedimentary depocentres. However, studying relict features could say little about the time‐scales over which this geomorphology evolved and under what glaciological conditions. This began to be addressed in the early 2000s, through continued efforts to scrutinise modern ice stream beds at higher resolution, but our current understanding of how landforms relate to processes remains subject to large uncertainties, particularly in relation to the mechanisms and time‐scales of sediment erosion, transport and deposition, and how these lead to the growth and decay of subglacial bedforms. This represents the next key challenge and will require even closer cooperation between glaciology, glacial geomorphology, sedimentology, and numerical modelling, together with more sophisticated methods to quantify and analyse the anticipated growth of geomorphological data from beneath active ice streams. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

On the morphological characteristics of overdeepenings in high‐mountain glacier beds

Overdeepenings, i.e. closed topographic depressions with adverse slopes in the direction of flow, are characteristic for glacier beds and glacially sculpted landscapes. Quantitative information about their morphological characteristics, however, has so far hardly been available. The present study provides such information by combining the analysis of (a) numerous bed overdeepenings below still existing glaciers of the Swiss Alps and the Himalaya‐Karakoram region modelled with a robust shear stress approximation and (b) detailed bathymetries from recently exposed lakes in the Peruvian Andes. The investigated overdeepenings exist where glacier surface slopes are low (< 5°–10°), occur in bedrock or morainic material and are most commonly a fraction of a kilometre squared in surface area, hundreds of metres long, about half the length in width and tens of metres deep. They form under conditions of low to high basal shear stresses, at cirque, confluence, trunk valley and terminus positions. The most striking phenomenon, however, is the high variability of their geometries: Depths, surface areas, lengths and widths of the overdeepenings vary over orders of magnitude and are only weakly – if at all – interrelated. Inclinations of adverse slopes do not differ significantly from those of forward slopes and are in many cases higher than so far assumed theoretical limits for supercooling of ascending water and corresponding closure of sub‐glacial channels. Such steep adverse slopes are a robust observation and in support of recently developed new concepts concerning the question about where supercooling of sub‐glacial water and closure of ice channels can or must occur. However, the question of when and under what climatic, topographic and ice conditions the overdeepenings had formed remains unanswered. This open question constitutes a key problem concerning the interpretation of observed overdeepenings, the understanding of the involved glacio‐hydraulic processes and the possibility of realistic predictive modelling of overdeepening formation. Copyright © 2016 John Wiley & Sons, Ltd.     

Limited modification of mid-latitude landscapes by ice sheets: The case of northeast Scotland

The paper uses a case study in Scotland to examine the amount and processes of landscape modification by Quaternary ice sheets. There is an inverse correlation between the distribution of landforms of glacial erosion and pre-glacial landscape remnants in northeast Scotland. The implication is that in places ice sheets can preserve a pre-glacial landscape unscathed, while elsewhere they remove the pre-glacial weathered rock. The location of glacial protection or erosion is strongly influenced by the topography and its influence on former ice sheet flow and basal thermal regime. The classic glacially eroded landscape of areal scouring can be produced by the removal of only 10–50 m of weathered rock. Furthermore rock basins, often regarded as the hallmark of glacial erosion, may be directly inherited from the pre-glacial pattern of deep weathering.     

A well‐preserved Eemian incised‐valley fill in the southern North Sea Basin,Belgian Continental Shelf ‐ Coastal Plain: Implications for northwest European landscape evolution

This paper demonstrates that the Belgian Continental Shelf and coastal plain occupy a key position between the depositional North Sea Basin and the erosional area of the Dover Strait as it is an area where erosional landforms and fragmented sedimentary sequences provide new evidence on northwest European landscape evolution. The study area hosts 20–30 m thick penultimate to last glacial sand‐dominated sequences that are preserved within the buried palaeo‐Scheldt Valley. Here, we build on the results of previous seismo‐ and lithostratigraphical studies, and present new evidence from biostratigraphical analysis, OSL dating and depth‐converted structure maps, together revealing a complex history of deposition and landscape evolution controlled by climate change, sea‐level fluctuations and glacio‐isostasy. This study presents strong new supportive evidence on the development of the incised palaeo‐Scheldt Valley landform that became established towards the end of the penultimate glacial period (MIS 6; Saalian) as a result of glacio‐isostatic forebulge updoming, proglacial lake drainage and subsequent collapse of a forebulge between East Anglia and Belgium following ice‐sheet growth, disintegration and retreat in areas to the north. The majority of the incised‐valley fill is of estuarine to shallow marine depositional context deposited during the transgression and high‐stand of the last interglacial (MIS 5e: Eemian). A thin upper part of the valley fill consists of last glacial (MIS 5d‐2: Weichselian) fluvial sediments that show a gradual decrease and retreat of fluvial activity to inland, upstream reaches of the valley system until finally the valley ceases to exist as the combined result of climate‐driven aeolian activity and possibly also glacio‐isostatic adjustment. Thus, strong contrasts exist between the palaeo‐Scheldt Valley and estuary systems of the penultimate glacial maximum to Last Interglacial (Saalian, Eemian), the beginning of the Last Glacial (Weichselian Early Glacial and Early‐Middle Pleniglacial), and the Last Glacial Maximum to Holocene. Copyright © 2018 John Wiley & Sons, Ltd.     

Pro‐glacial soil variability and geomorphic activity – the case of three Swiss valleys

Soils in pro‐glacial areas are often approached from a chronosequence viewpoint. In the chronosequence approach, the objective is to derive rates of soil formation from differences in properties between soils of different age. For this reason, in chronosequence studies, soils are sampled in locations that are assumed geomorphically stable and that have different age. As a result, these studies do not necessarily yield a complete view of soil variability in pro‐glacial areas, and may miss important relations between geomorphology and soil development. In this contribution, we present new soil observations from three closely related pro‐glacial areas in Switzerland. These observations were intended to get closer to a complete view of soil variability, and to assess impacts from factors other than time on soil development. About 40 soils were visited in each pro‐glacial valley in a combined design‐convenience sampling scheme and described in the field. Linear modelling was used to assess effects of time and topographic factors on soil properties. The time since glacial retreat turned out to rarely explain more than half of the variation in soil properties, and a linear model combining effects of time and topographic variables explained typically about half of the variation in each pro‐glacial valley. Models differed and were not transferable between valleys. Apparently, time and the present‐day shape of the landscape combined are insufficient information to accurately predict soil properties. Field evidence points to the importance of the geomorphic history and regime of the valleys as a reason for this. Copyright © 2014 John Wiley & Sons, Ltd.     

A comparison of solute fluxes and sources from glacial and non‐glacial catchments over contrasting melt seasons

Solute and runoff fluxes from two adjacent alpine streams (one glacial and one non‐glacial) were investigated to determine how the inorganic solute chemistry of runoff responded to seasonal and interannual changes in runoff sources and volume, and to differences in physical catchment properties. Intercatchment differences in solute composition were primarily controlled by differences in catchment geology and the presence of soils, whereas differences in total solute fluxes were largely dependent on specific discharge. The glacial stream catchment had higher chemical denudation rates due to the high rates of flushing (higher specific discharge). The non‐glacial Bow River had higher overall concentrations of solutes despite the greater prevalence of more resistant lithologies in this catchment. This is likely the result of both longer average water–rock contact times, and a greater supply of protons from organic soils and/or pyrite oxidation. Increases in snowpack depth/snowmelt runoff reduced the retention of nitrate in the Bow River catchment (i.e. increased nitrate export), probably by reducing net biological uptake, or by reducing the proportion of runoff that had contact with biologically active soil horizons that tend to remove nitrate. The two streams exhibited opposite solute flux responses to climate perturbations over three melt seasons (1998, 1999, and 2000). The 1998 El Niño event resulted in an unusually thin winter snowpack, and increased runoff and solute fluxes from the glacial catchment, but decreased fluxes from the Bow River catchment. Solute fluxes in the Bow River increased proportionally to discharge, indicating that increased snowmelt runoff in this catchment resulted in a proportional increase in weathering rates. In contrast, the proportional variation in solute flux in the glacial stream was only ∼70–80% of the variation in water flux. This suggests that increased ablation of glacier ice and the development of subglacial channels during the 1998 El Niño year apparently reduced the average water–rock contact time in the glacial catchment relative to seasons when the subglacial drainage system was primarily distributed in character. Copyright © 2005 John Wiley & Sons, Ltd.