Rate Of Deglaciation Research Articles

Assessing the state, parameter-interlinkages and dynamic shift of glaciers in the western Himalaya

The present study remotely investigates multiple parameters of 25 well-distributed glaciers in the western Himalaya (WH) for ~1990–2015 period to understand their state. Glacier changes for 1990–2000 (pre-2000) and 2000–2015 (post-2000) were analysed separately to assess the shift in glacier dynamics. Results reveal a clear glacier terminus retreat (13.4 ± 2.1 m/y), deglaciation (0.2 ± 0.1%/y), surface velocity slowdown (21%) and increase in number of glacial lakes (28 to 38) as well their area (5%/y) during 1990–2015. Notably, average terminus retreat and deglaciation rates were reduced post-2000 (retreat: 12.5 ± 3.2 m/y; deglaciation: 0.1 ± 0.1%/y) compared to pre-2000 (retreat: 15.3 ± 6.8 m/y; deglaciation: 0.2 ± 0.2%/y). Further, the average surface velocity was found to be increasing (0.9%/y) over 1990–2000. This can be linked with suggested less negative or slightly balanced condition of the WH glaciers over 1990s by previous studies. Present study also finds low downwasting during pre-2000 (−0.24 ± 0.14 m/y). Post-2000, studied glaciers downwasted (−0.50 ± 0.06 m/y) and slowed down (2%/y) significantly showing a decadal time lag in geometrical adjustment. Supraglacial debris, largely influenced by slope characteristics, increased with a uniform rate of 0.9%/y during pre- and post-2000 periods. Systematic analysis of interlinkages among various glacier parameters reveals a large interdependency. Further, presence of debris cover and pro- and supra-glacial lakes modify the glacier response to a great extent. Overall, the present study shows that (a) studied glaciers were in overall depleting phase over 1990–2015, (b) balanced conditions during 1990s prevailed over the studied glaciers, followed by accelerated mass wastage, (c) topography and local surface conditions introduce the intraregional heterogeneity in glacier response to climate.

Paraglacial transformation and ice‐dammed lake dynamics in a high Arctic glacier foreland, Gåsbreen, Svalbard

AbstractThe structural change of the Svalbard landscape from glacially dominated to paraglacial, that has taken place since the termination of the Little Ice Age (ca. 1900) is expressed by the widespread retreat of glaciers and progressive exposure of glacial landforms. This study provides insights into the rate of post‐LIA deglaciation and associated paraglacial transformation in the foreland of Gåsbreen, one of the first ever investigated glacier systems in the Arctic. Glacier is situated in Sørkapp Land (Southern Spitsbergen), a region characterized by one of the fastest deglaciation rates in the entire Svalbard Archipelago. During the investigated period, 1938–2020, the Gåsbreen was in a recession that accelerated after 1990, leading to its foreland increasing from 2.2 to 5.8 km2. The dynamics of landscape change in the glacier foreland, exposed since the end of LIA, manifested in i.e. the formation of ice‐dammed lakes, degradation in the surface of ice‐cored moraines and the landforms that are underlain by dead‐ice. Mass movements and debris flow on ice‐cored moraines and fluvioglacial processes had a great influence on this transformation. Enhanced proglacial runoff intensified denudation, transport and accumulation of sediments, which resulted, in: an increase in the extent of sandurs and proglacial riverbeds, an increase in the area of glacial lakes, extending and changing of the course of rivers in the glacier foreland. At the same time, the major lake system in the area—Goësvatnet underwent several cycles of filling and draining, often through glacial outburst flooding events.

Half-a-century (1971–2020) of glacier shrinkage and climatic variability in the Bhaga basin, western Himalaya

Glacier shrinkage is a globally occurring phenomena. High-resolution change detection based on frequent mapping and monitoring of high-altitude glaciers is necessary to precisely evaluate future water availability and to understand glacier evolution under different climatic scenarios in the Hindukush-Karakoram-Himalayan (HKH) region. This also holds true for the Bhaga basin of the western Himalaya. This study investigates glacier and glacier lake changes in the Bhaga basin, over the last five decades based on satellite imagery including Corona KH4 (1971), Landsat 7 Enhanced Thematic Mapper Plus (ETM+; 2000), Linear Imaging Self-Scanning Sensor (LISS IV; 2013), and Sentinel 2 (2020). Regional temperature and precipitation trends were evaluated from gridded climatic datasets (1900–2020). In the Bhaga basin 306 glaciers (>0.2 km2) were mapped with a total area of 360.3 ± 4.0 km2, of which 55.7 ± 0.6 km2 was covered with debris in 2013. The total glacier covered area decreased by ∼8.2 ± 1.5 % (0.16 ± 0.03 % yr−1) during the entire observation period 1971–2020, with noticeable heterogeneity between tributary watersheds. In the past two decades (2000–2020), the deglaciation rate has increased significantly (0.25 % yr−1) compared to the previous decades (1971–2000; 0.12 % yr−1). Glacier lake area increased by 0.6 ± 0.1 km2 (0.012 km2 yr−1) between 1971 and 2020. The NCEP/NCAR climatic data reveals an increase of 0.63°C in temperature and a decrease of 6.39 mm in precipitation for the period 1948–2018. In comparison, APHRODITE data shows an increasing trend in temperature of 1.14°C between 1961 and 2015 and decreasing trend in precipitation of 31 mm between 1951 and 2007. Both NCEP/NCAR and APHRODITE data reveal significant temperature increase and precipitation decrease since the 1990s, which have probably augmented ice loss in the Bhaga basin during the early 21st century.

Application of first arrival seismic tomography in a glaciated basin: implications for paleo-ice stream development

AbstractIce streams are sites of ice-sheet drainage and together with other processes, such as calving, have an impact on deglaciation rates and ice-sheet mass balance. Proglacial lake deposits provide records of ice-sheet deglaciation and have the potential to supplement other paleoclimate records. Oneida Lake, northeastern USA, contains a thick proglacial lake sequence that buries evidence of ice streaming and a paleo-calving margin that developed during retreat of the Laurentide Ice Sheet. Previous high-resolution digital elevation models identified the Oneida Ice Stream from glacial landforms northwest of the lake. In this study, we utilize seismic refractions from a multichannel seismic (MCS) reflection dataset to estimate the thickness of glacial deposits using seismic tomography. With this method we constrain the depth to top of Paleozoic strata, especially in areas where the reflection data yielded poor outcomes and validate our reflection data in regions of good coverage. We demonstrate that where long offset seismic data are available, the first-arrival tomography method is useful in studies of formerly glaciated basins. Our study identifies a ~108 m thick sedimentary section and potentially long paleoclimate record in Oneida Lake, and identifies a paleotopographic low that likely encouraged formation of the Oneida Ice Stream.

Biogeochemical Dynamics of a Glaciated High‐Latitude Wetland

AbstractHigh‐latitude, coastal wetland biogeochemistry is dynamic in response to climate change, and yet we do not understand, and thus cannot fully predict, how crucial aspects of these systems will change in the future. Temperatures in the Northern Hemisphere have disproportionately increased 4°C in 30 years causing the rate of deglaciation to increase significantly in global high‐latitude river deltas. This will have a prolonged effect on local microbiome metabolism and biodiversity of the subsurface, influencing solid, liquid, and gaseous compounds in the system. Using sediment geochemical analyses, autonomous sampling techniques, and 16S rRNA gene sequencing, we have identified key processes that occur in the Copper River Delta, AK, a model system to study high‐latitude watersheds during rapid climate change. We calculated carbon accumulation rates upwards of 520 ± 60 g C m−2 yr−1 in outwash pond sediments nearest to the glaciers, which co‐occurred with pronounced suboxic peaks in Fe (III) and Mn (II). Sediment microbial communities across the outwash ponds are structured on the basis of total iron and manganese concentrations, proximity to glaciers, and organic matter content. Additionally, we revealed no methane accumulation in the ponds during ice‐cover, despite high organic matter content. High‐latitude wetland ecosystems are not only influenced by the changing climate, but also have the potential to impact carbon cycling considering high carbon burial rates. These findings show the importance of understanding changing biogeochemical processes in high‐latitude wetlands, as they have the potential to influence carbon cycling.

Assessment of snout analysis of Himalayan glaciers: impact studies on Pindari, Kafni, Sundardhunga, and Baljuri base camp glaciers.

There are several causes for the increasing rate of deglaciation, such as global warming, increase in the concentration of black carbon, and extensive use of fossil fuels which causes the change in the overall climate system and shifting glacier ecosystem. This study was conducted on Pindari valley glaciers part of lesser Himalaya in Uttarakhand. This study investigates to (1) monitor and map change in the frontal length or the snout region of a glacier that can be studied with the help of remote sensing techniques and (2) evaluate the decadal and annual retreat rate of the glacier from 1972 to 2018. The study applies both the maximum likelihood classifier and NDSI spectral indices based classification for extracting the glacier region for different periods. This study reveals a significant amount of retreats taking place in the selected glaciers, Pindari, Sundardhunga, Kafni, and Baljuri base camp glaciers, from 1972 to 2018 as 1719.95m, 1751.21m, 1057.01m, and 810.78m, respectively. The highest amount of change is noticed in Pindari and Sundardhunga glaciers, higher than ~ 1700m. The study analyses temporal variation of the annual and decadal retreat rate in the Pindari valley glaciers, which would be helpful for the further study of the other glaciers.

A 10Be-dated record of glacial retreat in Connemara, Ireland, following the Last Glacial Maximum and implications for regional climate

Late Pleistocene stadials are traditionally viewed as periods of severe, year-round cooling centred on the North Atlantic, where the abrupt onset and termination of stadial conditions is widely, though not unequivocally, thought to reflect the impact on Atlantic Meridional Circulation of meltwater fluxes from fringing ice masses. Widespread subaerial melting of land-based ice, however, implies a state of glacier-climate disequilibrium that, according to some conceptual models, had to persist for the duration of each stadial event rather than immediately prior to it. Yet, to date, the spatial extent of any such melting remains unclear, as do the magnitude and duration of any cognate atmospheric warming. To help address this ambiguity, we reconstructed the timing and nature of deglaciation of the Connemara ice centre, Ireland, which, being located immediately downwind of the North Atlantic, affords an ideal location for assessing terrestrial expressions of stadial and non-stadial conditions. We report 15 new cosmogenic beryllium-10 ages of glacial erratic boulders that reveal the timing and magnitude of fluctuations in this sector of the former Irish ice sheet during Heinrich Stadial 1 (HS1). Coupled with geomorphic mapping, these results indicate rapid, widespread deglaciation of the Connemara ice centre at approximately 17 ka, under full stadial conditions, and that the ice margin retreated back to regional source areas rather than eastward towards the centre of the island. Recognising the potential divergence from the traditional view of cold North Atlantic stadials, we propose that the rate and nature of deglaciation in Connemara during HS1 likely reflects enhanced subaerial melting during the summer ablation season. • Fifteen new 10 Be dates constrain the deglaciation of an ice centre in Connemara, Ireland to Heinrich Stadial 1. • Our chronology contrasts with the current paradigm and suggests summer melting occurred despite anomalous winter cooling. • The Connemara deglacial record supports a broader pattern of enhanced seasonality in the North Atlantic during stadials.

Climate change impact on cryosphere and streamflow in the Upper Jhelum River Basin (UJRB) of north-western Himalayas.

The critical significance of keeping the current information about the extent and dynamics of the cryosphere in the Himalayas cannot be understated. The climate of the Himalayas is vulnerable and interlinked with global-scale climate changes, and the hydrology of the region mainly depends on the cryosphere. This is the first study that has created glacier and glacier lake inventory that links the impact of cryosphere on streamflow to land system dynamic changes under the changing climate of the Upper Jhelum River Basin (UJRB) of the Kashmir Himalayan region. This study uses a series of satellite data (1980-2016) to assess the depletion of snow cover area (SCA), deglaciation, and dynamics of glacial lakes. Moreover, observational long-term hydrometeorological data were used to understand the variability in temperature, precipitation, and track changes of land system dynamics under depletion of streamflow. The results suggested an overall rise in temperature (TMax = 0.05 ºC a-1; TMin = 0.02 ºC a-1; Tavg = 0.06 ºC a-1) and a decrease in precipitation (2.9mm a-1) between 1980 and 2016 with a significant increase in annual average temperature and decrease in annual precipitation at stations located at higher altitudes. The SCA showed a significantly decreasing (p < 0.01) trend in the glacierized sub-basins with an annual rate of decrease of -0.78% a-1, -0.15% a-1, -0.03% a-1 -0.90% a-1 for Lidder, Sindh, Vishow, and Rambiara sub-basins, respectively. The findings of this study reveal the high occurrence of glacier disintegration and deglaciation. During the period 2010-2016, a rapid rate of deglaciation was observed (18.34 ± 0.14 km2), followed by 1992-2000 (15.61 ± 0.13 km2). The average rate of retreat was observed to be 6.81 ± 1.5m a-1 with a total retreat of 267 ± 80m during 1980-2016, which is higher than reported from surrounding mountain ranges in the Himalayas. The mapped 244 glacial and high-altitude lake inventory covers a total surface area of around 15 km2, with 5.87 km2 (40%) covered by 25 bedrock-dammed lakes. The glacial expansion and creation of new lakes are observed to be because of increasing glacier and snow melting between 1980 and 2016, which increases the risk of GLOF events in the future. The annual average discharge in UJRB significantly increased from 1991 to 1998 and was observed to be higher than the annual average of the respected gauge stations but shows significant depletion from 1998 onwards. The streamflow depletion under climate change is one of the reasons for land system dynamics in UJRB. The area under agriculture has decreased up to 63% with a massive expansion of built-up (399%), aquatic vegetation (523%), and plantation (765%) between 1992 and 2015.

Geophysical Challenges for Future Satellite Gravity Missions: Assessing the Impact of MOCASS Mission

The GRACE/GRACE-FO satellites have observed large scale mass changes, contributing to the mass budget calculation of the hydro-and cryosphere. The scale of the observable mass changes must be in the order of 300 km or bigger to be resolved. Smaller scale glaciers and hydrologic basins significantly contribute to the closure of the water mass balance, but are not detected with the present spatial resolution of the satellite. The challenge of future satellite gravity missions is to fill this gap, providing higher temporal and spatial resolution. We assess the impact of a geodetic satellite mission carrying on board a cold atom interferometric gradiometer (MOCASS: Mass Observation with Cold Atom Sensors in Space) on the resolution of simulated geophysical phenomena, considering mass changes in the hydrosphere and cryosphere. Moreover, we consider mass redistributions due to seamounts and tectonic movements, belonging to the solid earth processes. The MOCASS type satellite is able to recover 50% smaller deglaciation rates over a mountain range as the High Mountains of Asia compared to GRACE, and to detect the mass of 60% of the cumulative number of glaciers, an improvement respect to GRACE which detects less than 20% in the same area. For seamounts a significantly smaller mass eruption could be detected with respect to GRACE, reaching a level of mass detection of a submarine basalt eruption of 1.6 109 m3. This mass corresponds to the eruption of Mount Saint Helens. The simulations demonstrate that a MOCASS type mission would significantly improve the resolution of mass changes respect to existing geodetic satellite missions.

A post-glacial relative sea level curve for the central Douglas Channel area, British Columbia, Canada

We present a post-glacial relative sea level (RSL) curve for the Central Douglas Channel region on the northern Northwest Coast of North America spanning the last 14,500 years, constrained by 68 radiocarbon dated index and limiting points derived from multiple RSL proxies. We evaluate each proxy based on the reliability and specificity of the inference for indicating RSL position, allowing us to weight different methods of data collection. We determine that central Douglas Channel was ice-free following the Last Glacial Maximum by ∼14,500 BP and RSL was at least 90 m higher than today. Isostatic rebound caused RSL to fall to 21 m asl by 11,500 BP, though there may have been a glacial re-advance that would have paused RSL fall around the beginning of the Younger Dryas. RSL fell to 10–15 m asl by 10,000 BP, and continued to drop at a slower rate towards its current position, which it reached by ∼1800 years ago. This is the first well-constrained RSL reconstruction for a nearly 400 km stretch between the northern and central coasts of British Columbia, Canada. This reconstruction provides information about the timing and rate of deglaciation of western North America and exemplifies variability among sea level trends across the region. Long-term RSL reconstructions are integral to explorations of coastal change and archaeological inquiries into early coastal occupations.

Subsurface architecture of alpine icy debris fans: Integration of ground-penetrating radar and surface observations in Alaska and New Zealand

Icy debris fans (IDFs) are extremely dynamic supraglacial landforms at the mouths of bedrock catchments between valley glaciers and icecaps. Recent studies quantified the nature, pace, and volume of mass flow processes contributing ice and sediment to IDFs by integrating field observations, drone and time-lapse imagery, and terrestrial laser scanning. New geophysical data presented herein characterize the subsurface architecture of IDFs along the McCarthy Glacier in Alaska and the Douglas, La Perouse, and Mueller Glaciers in New Zealand. Ground Penetrating Radar (GPR) profiles and soundings from field surveys during 2013–2015 provide stratigraphic evidence of the following subsurface processes important in the delivery of ice/sediment through IDFs to glaciers: (1) transitional merging of fan IDF deposits with valley glacier ice in cirque settings, (2) progressive rotation (slumping) of IDF deposits preserved in the hanging walls of downfan-dipping crevasses/faults, (3) progradation of IDF deposits onto valley glaciers, and (4) buttressing of IDF deposits along the margins of valley glaciers. The GPR-based variations in subsurface architecture/processes are attributable to previously documented variations in ice/sediment supply from overlying icecaps/nevés and topographic interaction with valley glacier junctions.GPR profiles document major depositional units separated by prominent reflectors, interpreted as annual ablation lag surfaces where lithic sediments are relatively concentrated. Most IDFs contain less than five depositional units, which may reflect the transitory nature of ice/sediment in these dynamic landforms. The uppermost depositional package on IDFs revealed by GPR is interpreted as deposition from the year preceding the GPR survey because it appears unconsolidated and lacks diffractions common in ice imaged in older (lower) units. The thickness of the upper package is typically less than or equal to the thickness of new deposits recorded by time-lapse cameras spread evenly over fan surfaces, which is expected accounting for compaction and ablation.The new geophysical datasets, together with previous surface observations, bear directly on managing hazards and evaluating glacial budgets in alpine regions. In particular, documentation of progressive rotation, buttressing, and progradation in IDF subsurface stratigraphies indicates protracted episodes of landform instability during deglaciation, consistent with previous surface observations documenting exceptionally high rates of depositional resurfacing and slumping. GPR-based evidence for subsurface deformation suggests ice/sediment is transferred to valley glaciers through IDFs. This interpretation is supported by previous visual evidence, including surface sliding/glacial flow and exceptionally high rates of deposition on IDFs coeval with little change in fan size. Integrating the role of IDFs in glacier budget models may result in forecasting a lower rate of deglaciation than traditionally recognized for some glaciers decoupled from icecaps.

Paraglacial rock-slope deformations: sudden or delayed response? Insights from an integrated numerical modelling approach

Glacial and paraglacial processes have a major influence on rock slope stability in alpine environments. Slope deglaciation causes debuttressing, stress and hydro-mechanical perturbations that promote progressive slope failure and the development of slow rock slope deformation possibly evolving until catastrophic failure. Paraglacial rock slope failures can develop soon after or thousands of years after deglaciation, and can creep slowly accelerating until catastrophic failure or nucleate sudden rockslides. The roles of topography, rock properties and deglaciation processes in promoting the different styles of paraglacial rock slope failure are still elusive. Nevertheless, their comprehensive understanding is crucial to manage future geohazards in modern paraglacial settings affected by ongoing climate change. We simulate the different modes and timing of paraglacial slope failures in an integrated numerical modelling approach that couples realistic deglaciation histories derived by modelling of ice dynamics to 2D time-dependent simulations of progressive failure processes. We performed a parametric study to assess the effects of initial ice thickness, deglaciation rate, rock-slope strength and valley shape on the mechanisms and timing of slope response to deglaciation. Our results allow constraining the range of conditions in which rapid failures or delayed slow deformations occur, which we compare to natural Alpine case studies. The melting of thicker glaciers is linked to shallower rockslides daylighting at higher elevation, with a shorter response time. More pronounced glacial morphologies influences slope lifecycle and favour the development of shallower, suspended rockslides. Weaker slopes and faster deglaciations produce to faster slope responses. In a risk-reduction perspective, we expect rockslide differentiation in valleys showing a strong glacial imprint, buried below thick ice sheets during glaciation.

Proglacial Lakes Control Glacier Geometry and Behavior During Recession

AbstractIce‐contact proglacial lakes are generally absent from numerical model simulations of glacier evolution, and their effects on ice dynamics and on rates of deglaciation remain poorly quantified. Using the BISICLES ice flow model, we analyzed the effects of an ice‐contact lake on the Pukaki Glacier, New Zealand, during recession from the Last Glacial Maximum. The ice‐contact lake produced a maximum effect on grounding line recession &gt;4 times further and on ice velocities up to 8 times faster, compared to simulations of a land‐terminating glacier forced by the same climate. The lake contributed up to 82% of cumulative grounding line recession and 87% of ice velocity during the first 300 years of the simulations, but those values decreased to just 6% and 37%, respectively, after 5,000 years. Numerical models that ignore lake interactions will, therefore, misrepresent the rate of recession especially during the transition of a land‐terminating to a lake‐terminating environment.

Delayed and rapid deglaciation of alpine valleys in the Sawatch Range, southern Rocky Mountains, USA

Abstract. We quantify retreat rates for three alpine glaciers in the Sawatch Range of the southern Rocky Mountains following the Last Glacial Maximum using 10Be ages from ice-sculpted, valley-floor bedrock transects and statistical analysis via the BACON program in R. Glacier retreat in the Sawatch Range from at (100 %) or near (∼83 %) Last Glacial Maximum extents initiated between 16.0 and 15.6 ka and was complete by 14.2–13.7 ka at rates ranging between 35.6 and 6.8 m a−1. Deglaciation in the Sawatch Range commenced ∼2–3 kyr later than the onset of rising global CO2 and prior to rising temperatures observed in the North Atlantic region at the Heinrich Stadial 1–Bølling transition. However, deglaciation in the Sawatch Range approximately aligns with the timing of Great Basin pluvial lake lowering. Recent data–modeling comparison efforts highlight the influence of the large North American ice sheets on climate in the western United States, and we hypothesize that recession of the North American ice sheets may have influenced the timing and rate of deglaciation in the Sawatch Range. While we cannot definitively argue for exclusively North Atlantic forcing or North American ice sheet forcing, our data demonstrate the importance of regional forcing mechanisms for past climate records.

Processes controlling the development of talus slopes in SW Spitsbergen: The role of deglaciation and periglacial conditions

AbstractThe location of Svalbard at the interface between the warm Atlantic and cold Arctic oceans causes the terrestrial environment to be highly sensitive to contemporary climate warming. Talus slopes provide a component of glaciated areas that has been registering these changes on a scale of several thousand years. However, knowledge about their development during glacial recession is still limited. This paper fills this gap by providing unique data obtained by geophysical methods: electrical resistivity tomography (ERT) and ground‐penetrating radar (GPR), regarding the talus slopes in Revdalen (SW Spitsbergen), which was last glaciated in the Neoglacial period. The results indicate that the thickness of talus slopes depends first of all on the size of the sediment supply area and only secondarily on the stage of development. The initial content of buried glacial ice in the talus deposits is differential and depends on both the rate of deglaciation and the local intensity of rock wall denudation. Over time, as a result of creep, the presence of massive ice is limited to ever lower parts of the slope. Above, there is aggradation of pore ice in delivered debris material. At the end of this stage, the buried glacial ice can form, or co‐create together with pore ice, the core of subslope rock glaciers. The relatively long period since the beginning of the Revdalen deglaciation allowed a general model of the development of talus slopes in the polar environment to be prepared.

Deglaciation Rate of Selected Nunataks in Spitsbergen, Svalbard—Potential for Permafrost Expansion above the Glacial Environment

Spitsbergen has recently experienced a continuous deglaciation process, linked to both glacier front retreat and lowering of the glacier surface. This process is accompanied by permafrost aggradation from the top of the slopes down to the glacier. Here, the authors determine the rate of permafrost expansion in this type of vertical profile. To this end, seven nunataks across the island were analysed using Landsat satellite imagery, a high-resolution digital elevation model (ArcticDEM), and geoinformation software. Over the last 24–31 years, new nunataks gradually emerged from the ice cover at an average linear rate of 0.06 m a−1 per degree of increment of the slope of the terrain at an average altitude of approximately 640 m a.s.l. The analysis showed that the maximum rate of permafrost expansion down the slope was positively correlated with the average nunatak elevation, reaching a value of approximately 10,000 m2 a−1. In cold climates, with a mean annual air temperature (MAAT) below 0 °C, newly exposed land is occupied by active periglacial environments, causing permafrost aggradation. Therefore, both glacial and periglacial environments are changing over time concomitantly, with permafrost aggradation occurring along and around the glacier, wherever the MAAT is negative.

Modern Deglaciation of the Altai Mountains: Effects and Possible Causes

The analysis of satellite images revealed by 25% decrease of the Altai mountains’ glaciation area over last 50 years. In 2008–2017, deglaciation rate increased twice. This tendency is in good agreement with an observed increase in the Katun River flow by 9% in 2008–2017 as compared to 1940–1968 (under invariable total precipitation). The analysis of trends of main meteorological parameters based on weather station data and the ERA-Interim reanalysis demonstrated that statistically significant warming in the region occurs only during the warm season and does not exceed 0.5°С/10 years. For this reason, “atmosphere–glaciers” turbulent heat transfer has increased by 4 W/m2 in last two years, that caused an annual melting layer increase by 100 mm water equivalent (w.e.). However, the main reason for the Altai mountains’ deglaciation is an increase of downward solar radiation flux, which amounted to 5 W/m2 per decade and increased the melting layer by 365 mm w.e. per year. A positive trend in net radiation agrees well with a decrease in cloud amount, which is associated with an increase in the moisture divergence flux and geopotential height and with the weakening of zonal winds in the middle troposphere.

Last interglacial (MIS 5e) and Holocene paleohydrology and paleovegetation of midcontinental North America from Gulf of Mexico sediments

Marine Isotope Stage (MIS) 5e experienced slightly higher insolation and lower global ice volume compared to the pre-industrial climate, however, how climate and vegetation differ regionally during MIS 5e is poorly known. Here, we examined MIS 5e hydrologic and vegetation conditions in the Mississippi River Basin (MRB) in North America (NA) and compared this to the Holocene. To do this, we reconstructed paleohydrology and paleovegetation conditions during the last 150 kyr using plant waxes (n-alkanes) and their carbon (δ13Calk) and hydrogen (δ2Halk) isotopic compositions preserved in ODP 625B core from the northwestern Gulf of Mexico (GOM) near the MRB discharge area. In the early MIS 5e, the jet stream may have been located in its northern most position over the last 150 kyr indicated by the highest δ2Halk values resulting in high warm and moist air mass contribution from the GOM and tropical Atlantic to the MRB. This is likely caused by northern expansion of the Atlantic Warm Pool and a northward shift of the ITCZ with increased warming. The climate may have become drier towards late MIS 5e with precipitation concentrated in summer expanding the C4 grassland distribution, which persisted through the MIS 5b. We speculate that this may be due to prolonged warming and drying from early to mid MIS 5. There were notable increases in plant wax export from western dry grasslands during deglaciation (13.5 ka) and the mid-Holocene (5.3 ka). In contrast, plant wax delivery was low during the MIS 5e without significant source changes. Such contrast between the two interglacials is due to the earlier completion of the LIS retreat in late MIS 6 and an onset of a full interglacial climate during the MIS 5e compared to the Last Glacial and Holocene transition. This study presents timing and mechanisms of precipitation source and vegetation composition changes over the NA as a result of slight changes in insolation and rates of deglaciation. These results provide a baseline for understanding the sensitivity of regional hydrologic and vegetation conditions to slight changes in global climate.

Proglacial lakes and the retreat pattern of the southwest Laurentide Ice Sheet across Alberta, Canada

The extent of proglacial lakes following the initial separation of the southwest Laurentide Ice Sheet from the Cordilleran Ice Sheet and its eastward retreat from the Canadian Rocky Mountains has been reconstructed across regions of Alberta at a range of scales. However, to date, no studies have integrated all available geological information to produce a province-wide deglacial reconstruction that considers the evolution of proglacial lakes as components of the ice-marginal system.In this paper, we utilize a geologically constrained shoreline projection method with a high resolution digital elevation model to reconstruct the evolution of the ice-marginal system along the southwest LIS during the last deglaciation. This method provides new details on the configuration, volume, drainage history and routing of ∼240 proglacial lakes as they migrated across Alberta and establishes a succession of paleogeographic reconstructions that can place geological evidence of regional and local ice-flow reorganizations into a spatiotemporal context. These reconstructions demonstrate that although the evolution of proglacial lakes was largely driven by the topography of the emerging landscape and the configuration of the ice margin, positive feedbacks in the ice-marginal system, particularly where margins transitioned from terrestrial to subaqueous settings played a major role in deglacial ice dynamics. Narrow, ribbon-shaped lakes that paralleled the ice margin induced relatively minor changes in style and rate of deglaciation, whereas the evolution of progressively larger, longer-lived lakes extending obliquely to the ice margin promoted surging and subsequent rapid retreat.

Ice-flow patterns and precise timing of ice sheet retreat across a dissected fjord landscape in western Norway

We reconstruct patterns of ice flow and retreat of the southwestern Scandinavian Ice Sheet, from 2900 field observations of glacial striae and elevation measurements of 60 ice-marginal-deltas from a high-resolution LiDAR DEM. During the Last Glacial Maximum, ice flow was towards the west across the entire area, including across several-hundred meter deep north-south oriented fjords. During deglaciation, ice flow adjusted to topography and the dominant flow direction switched towards the south-west. We use a shoreline diagram constructed from relative sea-level curves to establish the age of each delta, which allows us to constrain the timing of retreat with almost decadal precision. Rapid ice sheet retreat commenced at the onset of the Holocene at 11,600 cal years BP. Retreat rates were 160 m a−1 in the deepest fjords, 60–80 m a−1 in shallower fjords, and even slower for land-terminating margins. The fastest retreat rates, 240 m a−1 and 340 m a−1, were experienced in the largest fjords, Hardangerfjorden and Sognefjorden, which border the study area to the south and north. Crosscutting glacial striae indicate that calving bays developed during retreat along the widest fjords. The combination of complex fjord topography with fast ice-margin retreat by iceberg calving, led to isolation of ice remnants on islands and peninsulas, a process that accelerated the overall rate of deglaciation. Ice-margin retreat paused between 11,300–11,100 cal years BP, probably due to cooling during the Preboreal Oscillation.