Ice‐marginal lakes can affect glacier dynamics but are ignored in studies of the evolution of the Greenland ice sheet (GrIS) and of peripheral mountain glaciers and ice caps (PGICs). Here we show ...that lakes occupy 10% of the GrIS ice margin and occur on 5% of PGICs. Ice velocity at the GrIS margin is enhanced by ∼ 25% at lakes versus on land. Mean ice discharge into lakes is ∼4.9 Gt.yr, which is ∼1% of ice discharged through marine termini. We locate thousands of subglacial overdeepenings within which 7,404 km2 of future lakes could form, all of which will be ice‐marginal at some time. Future lakes in the west and east will be restricted to the margin of the GrIS and within alpine valleys, respectively. This status and possible future leads us to contend that lakes should be incorporated into projections of Greenland ice loss.
Plain Language Summary
Lakes forming on the edges of glaciers can affect ice dynamics but are not included within numerical models of the Greenland ice sheet (GrIS) or of peripheral mountain glaciers and ice caps (PGICs). We show that lakes presently occupy 10% of the GrIS ice margin and occur on 5% of PGICs. We determine that ice is typically moving 25% faster in the vicinity of lakes than at land margins. We estimate that ice discharge into lakes is ∼4.9 Gt.yr, which is 1% of that through marine glaciers. We identify thousands of depressions in the landscape that could hold future lakes as glaciers retreat from them. Overall, understanding of lake evolution should be incorporated into projections of Greenland ice loss and runoff.
Key Points
Lakes occupy 10% of the ice margin and vary between 2% and 26% by region
Ice velocity is ∼25% greater into lakes compared to onto land and total ice discharge into lakes is ∼4.9 Gt.yr−1
Thousands of overdeepenings have potential to host more and larger lakes in the future but not all will be ice‐marginal simultaneously
Climate change poses a considerable threat to the biodiversity of high altitude ecosystems worldwide, including cold‐water river systems that are responding rapidly to a shrinking cryosphere. Most ...recent research has demonstrated the severe vulnerability of river invertebrates to glacier retreat but effects upon other aquatic groups remain poorly quantified. Using new data sets from the European Alps, we show significant responses to declining glacier cover for diatoms, which play a critical functional role as freshwater primary producers. Specifically, diatom α‐diversity and density in rivers presently fed by glaciers will increase with future deglaciation, yet β‐diversity within and between sites will reduce because declining glacier influence will lower the spatiotemporal variability of glacier cover and its associated habitat heterogeneity. Changes in diatom assemblage composition as glacier cover declined were associated strongly with increasing riverbed stability and water temperature. At the species level, diatoms showed a gradation of responses; for example, Eunotia trinacria, found exclusively at river sites with high (≥52%) catchment glacier cover, may be affected negatively by ice loss. Conversely, seven taxa confined to sites with no glacier cover, including Gomphonema calcareum, stand to benefit. Nineteen (22%) taxa were noted as threatened, endangered, rare or decreasing on the Red List of Algae for Germany, with most at sites ≤26% glacier cover, meaning further ice loss may benefit these diatoms. However, six taxa found only in rivers ≥28% glacier cover may require reclassification of their Red List conservation status, as this habitat is threatened by deglaciation. Our identification of clear links between decreasing glacier cover and river diatom biodiversity suggests there could be significant reorganization of river ecosystems with deglaciation, for example, through alterations to primary production, biogeochemical cycles, and the shifting resource base of alpine freshwater food webs which lack significant allochthonous energy inputs.
Glacier retreat is known to threaten the biodiversity of freshwaters worldwide, but little is understood of the effects on algal communities that play key roles in primary production, biogeochemical cycling and the resource base of alpine freshwater food webs. This research investigated the biodiversity of benthic diatom assemblages in rivers of the European Alps to determine their response to reducing catchment glacier cover. Diatom α‐diversity and density increased but β‐diversity was reduced. Six conservation Red‐List taxa may require reclassification as they were found exclusively ≥28% glacier cover and are therefore under threat from ongoing glacier loss.
•We explore sediment transport and morphological changes during multiple floods.•Flood sequences lead to cumulative effects on changes in river bed geometry.•In-channel scour dramatically increased ...the conveyance capacity of a channel.•Cumulative channel adjustments greatly affect local and downstream flood hazards.
Effects of flood-induced bed elevation and channel geometry changes on flood hazards are largely unexplored, especially in the case of multiple floods from the same site. This study quantified the evolution of river channel and floodplain geometry during a repeated series of hypothetical extreme floods using a 2D full hydro-morphodynamic model (LHMM). These experiments were designed to examine the consequences of channel geometry changes on channel conveyance capacity and subsequent flood dynamics. Our results revealed that extreme floods play an important role in adjusting a river channel to become more efficient for subsequent propagation of floods, and that in-channel scour and sediment re-distribution can greatly improve the conveyance capacity of a channel for subsequent floods. In our hypothetical sequence of floods the response of bed elevation was of net degradation, and sediment transport successively weakened even with floods of the same magnitude. Changes in river channel geometry led to significant impact on flood hydraulics and thereby flood hazards. We found that flood-induced in-channel erosion can disconnect the channel from its floodplain resulting in a reduction of floodwater storage. Thus, the frequency and extent of subsequent overbank flows and floodplain inundation decreased, which reduced downstream flood attenuation and increased downstream flood hazard. In combination and in summary, these results suggest that changes in channel capacity due to extreme floods may drive changes in flood hazard. The assumption of unchanging of river morphology during inundation modelling should therefore be open to question for flood risk management.
The Jostedalsbreen ice cap is mainland Europe's largest ice cap and accommodates 20% (458 km2 in 2019) of the total glacier area of mainland Norway. Jostedalsbreen and its meltwater contribute to ...global sea-level rise and to local water management, hydropower and tourism economies and livelihoods. In this study, we construct a digital terrain model (DTM) of the ice cap from 1966 aerial photographs, which by comparing to an airborne LiDAR DTM from 2020, we compute changes in surface elevation and geodetic mass balances. The area mapped in both surveys cover about 3/4 of the ice cap area and 49 of 82 glaciers. The measured glacier area has decreased from 363.4 km2 in 1966 to 332.9 km2 in 2019, i.e. a change of −30 km2 or −8.4% (−0.16% a−1), which is in line with the percentage reduction in area for Jostedalsbreen as a whole. The mean geodetic mass balance over the 49 glaciers was −0.15 ± 0.01 m w.e. a−1, however, large variability is evident between glaciers, e.g. Nigardsbreen (−0.05 m w.e. a−1), Austdalsbreen (−0.28 m w.e. a−1) and Tunsbergdalsbreen (−0.36 m w.e. a−1) confirming differences also found by the glaciological records for Nigardsbreen and Austdalsbreen.
Ice thickness distribution and volume are both key parameters for glaciological and hydrological applications. This study presents VOLTA (Volume and Topography Automation), which is a Python script ...tool for ArcGISTM that requires just a digital elevation model (DEM) and glacier outline(s) to model distributed ice thickness, volume and bed topography. Ice thickness is initially estimated at points along an automatically generated centreline network based on the perfect-plasticity rheology assumption, taking into account a valley side drag component of the force balance equation. Distributed ice thickness is subsequently interpolated using a glaciologically correct algorithm. For five glaciers with independent field-measured bed topography, VOLTA modelled volumes were between 26.5% (underestimate) and 16.6% (overestimate) of that derived from field observations. Greatest differences were where an asymmetric valley cross section shape was present or where significant valley infill had occurred. Compared with other methods of modelling ice thickness and volume, key advantages of VOLTA are: a fully automated approach and a user friendly graphical user interface (GUI), GIS consistent geometry, fully automated centreline generation, inclusion of a side drag component in the force balance equation, estimation of glacier basal shear stress for each individual glacier, fully distributed ice thickness output and the ability to process multiple glaciers rapidly. VOLTA is capable of regional scale ice volume assessment, which is a key parameter for exploring glacier response to climate change. VOLTA also permits subtraction of modelled ice thickness from the input surface elevation to produce an ice-free DEM, which is a key input for reconstruction of former glaciers. VOLTA could assist with prediction of future glacier geometry changes and hence in projection of future meltwater fluxes.
•The VOLTA model is presented to estimate ice thickness distribution and volume.•The model requires just glaciers outlines and a DEM as inputs.•VOLTA can rapidly estimate volume and bed topography for multiple glaciers.•The model was tested against glaciers with independent field based observations.•Modelled volumes were between 26.5% (underestimate) and 16.6% (overestimate).
Accelerated glacier and ice sheet retreat and thinning in recent decades has profound consequences for catchment sediment supply with attendant repercussions for nutrient cycling, carbon fluxes and ...natural resource management. This paper evaluates the impacts of deglaciation on sediment yields from glaciated, deglaciating and recently-deglaciated catchments. It summarises the key characteristics of sediment yields from glaciated catchments to be that they span five orders of magnitude, vary with latitude and are greatest in high-relief and tectonically-active regions. We review the available quantitative data on sediment yields from glaciated catchments and we comment extensively on spatio-temporal variability to understand global to local and inter- and intra-catchment controls. Significant gaps in the available sediment yield data and also in our knowledge of sediment sources, pathways and sinks are identified. We constrain a set of novel approaches by which these gaps could be addressed. In particular, we suggest that the opportunities presented by emerging datasets and analytical methods enabling landcover changes, Digital Elevation Model (DEM) change detection, analyses of connectivity and analyses of sediment plumes are exciting and these approaches should become practical tools for understanding intra- and inter-catchment sediment yields from deglaciating landscapes. We showcase preliminary studies utilising these datasets and they are used to formulate hypotheses designed to stimulate further research.
•Explanation of importance of understanding sediment yields from deglaciating catchments.•Concise summary of the key characteristics of sediment yields from glaciated catchments.•Review of spatio-temporal trends in global data compilation.•Identification of gaps in the available data and therefore our knowledge.•Suggestions for capitalising on emerging datasets and methods to address these gaps.
There is a lack of detailed information on the Holocene evolution of the west Greenland ice margin, not least because it was farther inland than at present and thus at present is covered by ice. ...Suggestions have been put forward, both of relatively quick ice-margin retreat and of relatively stable ice-margin positions. This study presents the first exploitation of sediments from an ice-dammed basin. Sediment that is rich in organic material records a period of time sufficient for vegetation to thrive; thus, without a lake and with a distal and diminished ice mass relative to the present. In contrast, sediment composed of suspension-settling deposits and with drop stones records an ice-dammed lake and a proximal calving ice-margin. Overall, we evidence relatively coarse-grained glacifluvial sedimentation predominantly from an ice-marginal delta and/or a proglacial, braided river soon after the early Holocene deglaciation of this area. Subsequent mid-Holocene aeolian activity deposited leaves, some vegetation (roots) developed in the basin, and ponding of water formed organic-rich "gyttja" sediments and thin layers of peat. The lake then became ice marginal, and ice advances are interpreted to have occurred at approximately 4000 cal. yr BP, and between 2776 ± 26 cal. yr BP and 2440 ± 45 cal. yr BP. The Little Ice Age ice-margin advance at Russell Glacier apparently reached its maximum extent after 147 ± 39 cal. yr BP, whereas at Isunnguata Sermia the maximum could have been as early as 245 ± 26 cal. yr BP. Given that ice-marginal lakes are becoming globally ubiquitous, improved resolution and new geological data on lake and ice-margin dynamics and interactions, such as ice-dammed lake-basin sediments, will be important for numerical models for assessing past and future ice-mass evolution.
Abstract
Rapid changes observed today in mountain glaciers need to be put into a longer-term context to understand global sea-level contributions, regional climate-glacier systems and local landscape ...evolution. In this study we determined volume changes for 400 mountain glaciers across the Southern Alps, New Zealand for three time periods; pre-industrial “Little Ice Age (LIA)” to 1978, 1978 to 2009 and 2009 to 2019. At least 60 km
3
± 12 km
3
or between 41 and 62% of the LIA total ice volume has been lost. The rate of mass loss has nearly doubled from − 0.4 m w.e year
−1
during 1,600 to 1978 to − 0.7 m w.e year
−1
at present. In comparison Patagonia has lost just 11% of it’s LIA volume. Glacier ice in the Southern Alps has become restricted to higher elevations and to large debris-covered ablation tongues terminating in lakes. The accelerating rate of ice loss reflects regional-specific climate conditions and suggests that peak glacial meltwater production is imminent if not already passed, which has profound implications for water resources and riverine habitats.