Regional Equilibrium Line Altitude Research Articles

Reconstruction of glaciers in the western boundary of the Altiplano (18.5°-19°S): Singularities and insights on potential drivers of past advances

Today, glaciers in the western Altiplano are very scarce, even on peaks exceeding 6000 m. In this mostly ice-free landscape, however, moraines and other glacial deposits are commonly found attesting to quite different climate conditions that favored the advance of glaciers in the past. Although other areas of the Altiplano have been commonly accounted for regarding paleoclimate reconstructions, western Altiplano has been often overlooked. Here we present a detailed map of glacial landforms and a paleoglacier reconstruction from a test area in the western Altiplano between 18.5° and 19°S. We reconstructed regional equilibrium line altitudes (ELA) for several moraine stages, representing extensive past glacier advances in the region. During a prominent and ubiquitous ‘Principal Moraine (PM)’ stage glaciers from most peaks and all orientations advanced to altitudes of up to 4000 m asl. Reconstructed PM ELAs along 90 valleys range from 4400 to 5000 m asl. ELA distribution reveals a strong aspect-dependency at the western boundary of the Altiplano, with ELAs 300 m lower on west-facing glaciers than east-facing glaciers. The coincidence of such steep gradient with a prominent NW-SE ridge accounts for topographic control on the precipitation and thus in the advance of glaciers along the western boundary of the Altiplano. To the east, ELAs from the peaks towering the Altiplano are comparable to ELAs from east-facing glaciers at the western Altiplano boundary but show little to no aspect-dependency. Since these patterns cannot be solely explained by increased moisture advection from the Amazon basin, we suspect that western-sourced moisture associated with increased frequency of cold fronts and cut-off events played an important role in the glacier dynamics at this latitude. Nevertheless, further investigations are required to evaluate the relative role of both precipitation regimes on the glacier dynamics from the westernmost Altiplano.

The impact of local topoclimatic factors on marginal Pleistocene glaciation in the Northern Romanian Carpathians

In the Northern Romanian Carpathians (NRC) small Pleistocene cirque glaciers have formed in several isolated mountain massifs exceeding 1800 m asl. This paper brings forward new geomorphological evidence of marginal glaciation in the Călimani, Suhard, and Gurghiu Mountains which are the southernmost glaciated areas of the NRC.We reconstructed the extent and ice-surface geometry for 12 palaeo-glaciers in the study area during the Maximal Ice Extent (MIE) which is attributed to LGM. The equilibrium line altitudes (ELAs) calculated using the area-altitude-balance-ratio (AABR) 1.6. were between 1740 m and 1870 m which is the highest value in the Eastern Carpathians and one of the highest in the entire Carpathians. The specific features of the Călimani and Gurghiu Mountains are the extensive, gentle outer slopes of volcanic calderas that play an important role as additional areas for snowblow accumulation. Inclusion of the potential snow contribution area (snowblow and avalanche accumulation) in the ELA calculation resulted in an additional upward ELA shift up to 66 m (up to 40% of the glacier elevation range).Our data indicate that small cirque glaciation in the NRC could develop when mountain ranges are at least 100 m above the extrapolated regional ELA. The general ELA pattern in the NRC shows a positive trend (∼3 m/km) which could reflect both the general rise in temperature and the starvation of precipitation toward the southeast. However, the observed eastward ELA rise on W-E transects in the Eastern Carpathians shows a pure precipitation effect and indicates the influence of zonal atmospheric circulation in the far interior of Europe during the LGM.

DISTRIBUTION AND EVOLUTION OF ICE APRONS IN A CHANGING CLIMATE IN THE MONT-BLANC MASSIF (WESTERN EUROPEAN ALPS)

Abstract. Ice Apron (IA) is a poorly studied ice feature, commonly existing in all the world’s major mountain regions. This study aims to map the locations of the IAs in the Mont Blanc massif (MBM), making use of the very high-resolution optical satellite images from 2001, 2012 and 2019. 423 IAs were identified and accurately delineated in the MBM on the images from 2019, and their topographic characteristics were studied. We generated our own Digital Elevation Model (DEM) at 4 m resolution since the freely available products predominantly suffer from significant inconsistencies, especially in steep mountain areas. Results show that most IAs exist at elevations above the regional Equilibrium Line Altitude (ELA), on steep slopes, on concave surfaces, on northern and southern aspects and on the most rugged terrains. They are also commonly associated with steep slope glaciers as 85% of them occur on these glaciers’ headwalls. A comparison between 2001 and 2019 shows that IAs have lost around 29% of their area over a period of 18 years. This is significant and the rate of area loss is very alarming in comparison with the larger glacier bodies. We also studied the effect of topographic parameters on the area loss. We found that topographic factors like slope, aspect, curvature, elevation and Terrain Ruggedness Index (TRI) strongly influence the rate of area loss of IAs.

Holocene glacier and ice cap fluctuations in southwest Greenland inferred from two lake records

Glaciers and ice caps (GICs) respond rapidly to changes in temperature and precipitation. Thus, records of their past fluctuations yield valuable information on past climate. However, relatively little is known about the long-term, Holocene history of Greenland’s local GICs, peripheral to the Greenland Ice Sheet. Here we report sediment records of Holocene glacier activity from two distally fed glacial lakes near Buksefjord, southwest Greenland. The two lakes’ watersheds host modern GICs of contrasting size. The Pers Lake (informal name) watershed drains part (3 km2) of a single small ice cap. In contrast, nearby Lake T3’s (informal name) watershed drains numerous GICs totaling 100 km2. At the time it emerged from the sea ∼8.6 ka BP, Pers Lake was receiving no glacial meltwater input. Sediment physical and geochemical properties indicate persistent meltwater input and regrowth of the ice cap within the Pers Lake catchment beginning at ∼1.4 ka BP, after almost 3000 years of sporadic meltwater input beginning ∼4.3 ka BP. The ice cap above Pers Lake reached a maximum late Holocene extent during the final phase of the Little Ice Age (LIA), ∼0.1 ka BP. The complementary Lake T3 sediment record suggests continued meltwater input from the larger suite of upstream GICs from the time of the lake’s isolation from the sea ∼8.4–7.5 ka BP through to the present. This indicates that some GICs here probably survived the Holocene Thermal Maximum (HTM), although were significantly reduced in size for an extended period (of unknown age and duration). Combined with evidence from Pers Lake and prior studies that show GICs at low and intermediate elevations in this region melted away completely during the HTM, and evidence for GIC presence at Lake T3, we provide lower and upper bounds on regional HTM equilibrium-line altitudes (ELAs). We estimate that regional ELAs were between ∼1370 and 1470 m above sea level in the early-to-middle Holocene. From the middle to late Holocene, our results, along with other regional GIC studies, indicate progressive lowering of regional glacier ELAs in response to Neoglacial summer cooling of ∼2.7 °C, assuming no change in precipitation.

Morphometry and spatial distribution of glacial cirques in the Western Fuegian Andes of Argentina, southernmost South America

Glacial cirques comprise a distinctive element of relief in alpine landscapes, and their morphometry, distribution and spatial arrangement are generally linked to morphoclimatic factors. This study provides an analysis of morphometric characteristics, floor altitude, aspect and distribution of 251 first-order and 99 "inner" cirques with regard to environmental agents in a selected alpine-type area of the southwestern Fuegian Andes, in Southern Argentina. Parametric variables were extracted and quantified from satellite images and a Digital Elevation Model, and integrated to a Geographical Information System for their systematization and processing. Significant associations were identified between the different geometric and spatial parameters evaluated. Variability in cirque morphometry along the study area suggests that their evolution was regulated to a large degree by glacial dynamics related to Beagle palaeoglacier activity, and by the W-E topographic gradient of the mountain section. The record indicates that most cirques in the Fuegian Andes have undergone allometric enlargement in the three dimensions (i.e., lesser vertical erosion compared to walls' horizontal expansion), with local variations in size and shape controlled by ice residence time, aspect and altitude, whereas geological structure is supposed to have affected cirque formation along the entire mountain chain. Cirque azimuth potentially responds to the combined influence of climate (solar radiation and westerly winds), structure (pervasive faulting and rock bedding) and topography in their development. In particular, dominance of southeasterly aspects is interpreted as a consequence of a relatively cloud-free atmosphere and aridity during past glacial periods in the region, consistently with published evidence. The absence of aspect-related variation in cirque floor altitude is considered to indicate incidence of extensive glaciations in the area (i.e., mostly ice-sheet conditions, with mountain peaks far above the regional equilibrium line altitude or ELA), as well as formation of younger, lower cirques during successive cold periods of different intensity. Occurrence of "inner" cirques inside larger (preexisting) cirques is thought to reflect a temporary positioning of the ELA just above the floor of the container cirque, permitting short-lived, small cirque-type glaciers to develop along structural bedrock weaknesses.

Late Pleistocene glaciers to present-day snowpatches: a review and research recommendations for the Marrakech High Atlas

There are no glaciers today in the High Atlas, Morocco. However, there is evidence that niche glaciers and late-lying snowpatches in the High Atlas were present as recently as the last century and there are at least four sites where snowpatches appear to survive some summer seasons today. Many other sites also support non-perennial late-lying snow below steep shaded north and northeast-facing cliffs at altitudes > 3100 m. Coarse sediment ridges interpreted as moraines or pronival ramparts enclose most of these snowpatches. These landforms most likely record the positions of former niche glaciers and late-lying snowpatches in the Little Ice Age. The niche glaciers and late-lying snowpatches survived below the regional equilibrium line altitude because of strong local topoclimatic controls. In addition to strong shading, many of the current late-lying snowpatches are fed by long deep gullies which funnel avalanching snow from the cirque backwalls. The disappearance of many perennial snowpatches in the last few decades coincides with a strong trend towards warmer summer air temperatures since the 1970s (> 2 °C). However, inter-annual changes in snowpack mass balance are affected by very large variations (> 400% variability) in winter precipitation. A new research programme is underway investigating the history of late-lying snow and cirque glaciers in the High Atlas. A particular focus of this research is to utilise geomorphological and geochronological evidence to understand fluctuations in snow and ice through the Holocene and link this to continuous records of environmental change in the High Atlas region.

Glacial Geomorphology and Preliminary Glacier Reconstruction in the Jablanica Mountain, Macedonia, Central Balkan Peninsula

Although glacial landforms on the Balkan Peninsula have been studied since the 19th century, only scarce data are available about the extent of the former glaciations in the Central Balkan Peninsula, the transition zone between the Mediterranean and Central Europe. Glacial features of the Jablanica Mt. were mapped, described and classified into morphostratigraphic units. A revised glacio-geomorphological map was produced and glacial landforms were assigned to six morphostratigraphic units. Ten primary and two secondary cirques were identified in the upper parts of the studied valleys, while downstream the valleys were steep and glacially shaped with several glacial steps and thresholds. Cirque and valley morphology indicate that subglacial deepening was limited within the cirques and was more intensive in the valley sections during more extensive glacial phases. The largest reconstructed glaciers were 4.6–7 km long, while the last cirque glaciers were only a few hundred meters long. Using morphostratigraphic data, a glacier reconstruction was carried out for the largest mapped glacial extent. On the basis of glacial geomorphology, a former equilibrium-line altitude (ELA) of ~1800 m and glacier cover of 22.6 km2 were estimated during this stage. The local ELA values were compared to the regional ELA record and enabled to tentatively attribute a MIS 6 age for the reconstructed maximum ice extent in the study area.

Reconstruction and paleoclimatic significance of late Quaternary glaciers in the Tararua Range, North Island, New Zealand

Reconstructed mountain glaciers are routinely used as a proxy for climate in the search for evidence of interhemispheric climate fluctuations during the Quaternary. In New Zealand, valley glaciers at the Last Glacial Maximum (LGM) extended from an ice sheet centred on New Zealand's Southern Alps to below present-day sea level. In contrast, evidence of LGM glacial activity on the North Island is rare. Here, a glacioclimatic reconstruction is presented of two former glaciers in the Tararua Range (41°S) in the southern North Island. At Mt. Aston, an isolated cirque basin contains landform evidence of a marginal niche glacier. At Park Valley, the lateral moraine of a larger cirque glacier has yielded published cosmogenic isotope ages. The paleoglacier reconstruction shows that paleoequilibrium line altitudes increased northwards across New Zealand during the local LGM. Hence, at this latitude only topography >1200 m above present day sea-level was of sufficient elevation to allow small glaciers to form. The Mt Aston glacier covered only 0.18 km2 with an equilibrium line altitude (ELA) of c. 1287 m above present sea level. A mean ELA glacier thickness of c. 25 m gives a basal shear stress at the ELA of c. 100 kPa−1, with a mean summer (December, January, February, DJF) temperature at the ELA of no lower than 5.5 °C below present, below which precipitation would have been insufficient to support the reconstructed glacier. Implied LGM paleo-temperatures from both the ELA reconstruction and the glaciological reconstruction broadly accord with other paleo-temperature proxies from the North Island. Park Valley glacier covered c. 0.45 km2 with an ELA of c. 1270 m and a mean ELA basal shear stress of 65 kPa. Its balance discharge was 9 × greater than at Mt Aston. It appears to have been glaciologically viable across a wide range of paleotemperatures: thus, the more marginal glacier is a more useful paleoclimatic indicator because it places a maximum limit on LGM temperature depression, which the larger glacier does not. ELAs of both glaciers closely approximate the regional LGM ELA trend surface. The paleo-glacier reconstructions imply that together with temperature driving LGM paleo-ELA depression, changes in south-westerly airflow over New Zealand, bringing moisture-laden but cool air, maximized snowfall and minimised winter melt. The corollary is that patterns of Quaternary glacier fluctuations may be interpreted as responses, at least in-part, to precipitation-driven changes, and secondly, North Island glaciation was probably more extensive than previously assumed.

Hypsometric control on glacier mass balance sensitivity in Alaska and northwest Canada

AbstractGlacier hypsometry provides a first‐order approach for assessing a glacier's response to climate forcings. We couple the Randolph Glacier Inventory to a suite of in situ observations and climate model output to examine potential change for the ∼27,000 glaciers in Alaska and northwest Canada through the end of the 21st century. By 2100, based on Representative Concentration Pathways (RCPs) 4.5–8.5 forcings, summer temperatures are predicted to increase between +2.1 and +4.6°C, while solid precipitation (snow) is predicted to decrease by −6 to −11%, despite a +9 to +21% increase in total precipitation. Snow is predicted to undergo a pronounced decrease in the fall, shifting the start of the accumulation season back by ∼1 month. In response to these forcings, the regional equilibrium line altitude (ELA) may increase by +105 to +225 m by 2100. The mass balance sensitivity to this increase is highly variable, with the most substantive impact for glaciers with either limited elevation ranges (often small (<1 km2) glaciers, which account for 80% of glaciers in the region) or those with top‐heavy geometries, like icefields. For more than 20% of glaciers, future ELAs, given RCP 6.0 forcings, will exceed the maximum elevation of the glacier, resulting in their eventual demise, while for others, accumulation area ratios will decrease by >60%. Our results highlight the first‐order control of hypsometry on individual glacier response to climate change, and the variability that hypsometry introduces to a regional response to a coherent climate perturbation.

The gravel fans of upper Glen Roy, Lochaber, Scotland: their importance for understanding glacial, proglacial and glaciolacustrine dynamics during the Younger Dryas cold period in an Atlantic margin setting

Over the past 200 years significant research effort has gone into explaining the origin of the obvious former shorelines in Glen Roy (the so-called “Parallel Roads”). The large gravel deposit at the mouth of Glen Turret has attracted similar interest, but a solution to its origin and age remains contested: the same applies to the associated gravel fans in upper Glen Roy. This paper presents the results of systematic mapping and instrumental levelling of these features and new evidence from two previously unrecorded gravel fans in the nearby Allt Chonnal valley. Interpretation of altitudinal and lithofacies data indicates that all the fans (including the one at the mouth of Glen Turret) were deposited in a series of ice-dammed lakes during the Loch Lomond Stadial (Younger Dryas). The largest gravel fans were deposited in the shallow lake heads of the 260m, 325m and 350m lakes in upper Glen Roy, infilling these areas to the extent that deposition was in part subaerial. The absence of foreset bedding from the deposits is explained by the relatively shallow depth of lake waters, which inhibited development of classic Gilbert-type deltas and encouraged Hjulström-type fans. The previous assumption that gravel deposition into the 325m and 350m lakes was relatively limited is shown to be erroneous. The Allt Chonnal gravel fans, deposited into these lakes, have an estimated combined volume of 6,000,000m3 deposited in about 200 years. Significant gravel deposition into these former lakes continued at least until glaciers started to retreat in Glen Roy. When glacier retreat began, gravel deposition was reduced by over 85%. This research also concludes that the glacier which deposited the Turret fan emanated from Glen Gloy, supplied not only by ice from the south but also from ice that spilled into Glen Gloy from the adjacent Great Glen. The Loch Lomond Stadial (LLS) ice cap mapped in the Monadhliath Mountains was connected to the west Highland ice cap, which at the time filled the Great Glen. Regional equilibrium line altitude (ELA) estimates, based solely on the Monadhliath LLS ice cap evidence, may therefore be altitudinally too high.

Holocene glacier variability and Neoglacial hydroclimate at Ålfotbreen, western Norway

Glaciers and small ice caps respond rapidly to climate perturbations (mainly winter precipitation, and summer temperature), and the mass-balance of glaciers located in western Norway is governed mainly by winter precipitation (Pw). Records of past Pw can offer important insight into long-term changes in atmospheric circulation, but few proxies are able to accurately capture winter climate variations in Scandinavia. Reconstructions of equilibrium-line-altitude (ELA) variations from glaciers that are sensitive to changes in Pw therefore provide a unique opportunity to quantify past winter climate in this region. Here we present a new, Holocene glacier activity reconstruction for the maritime ice cap Ålfotbreen in western Norway, based on investigations of distal glacier-fed lake sediments and modern mass balance measurements (1963–2010). Several lake sediment cores have been subject to a suite of laboratory analyses, including measurements of physical parameters such as dry bulk density (DBD) and loss-on-ignition (LOI), geochemistry (XRF), surface magnetic susceptibility (MS), and grain size distribution, to identify glacial sedimentation in the lake. Both radiocarbon (AMS 14C) and 210Pb dating were applied to establish age-depth relationships in the sediment cores. A novel approach was used to calibrate the sedimentary record against a simple ELA model, which allowed reconstruction of continuous ELA changes for Ålfotbreen during the Neoglacial (when Ålfotbreen was present, i.e. the last ∼1400 years). Furthermore, the resulting ELA variations were combined with an independent summer temperature record to calculate Neoglacial Pw using the ‘Liestøl equation’. The resulting Pw record is of higher resolution than previous reconstructions from glaciers in Norway and shows the potential of glacier records to provide high-resolution data reflecting past variations in hydroclimate. Complete deglaciation of the Ålfotbreen occurred ∼9700 cal yr BP, and the ice cap was subsequently absent or very small until a short-lived glacier event is seen in the lake sediments ∼8200 cal yr BP. The ice cap was most likely completely melted until a new glacier event occurred around ∼5300 cal yr BP, coeval with the onset of the Neoglacial at several other glaciers in southwestern Norway. Ålfotbreen was thereafter absent (or very small) until the onset of the Neoglacial period ∼1400 cal yr BP. The ‘Little Ice Age’ (LIA) ∼650–50 cal yr BP was the largest glacier advance of Ålfotbreen since deglaciation, with a maximum extent at ∼400–200 cal yr BP, when the ELA was lowered approximately 200 m relative to today. The late onset of the Neoglacial at Ålfotbreen is suggested to be a result of its low altitude relative to the regional ELA. A synthesis of Neoglacial ELA fluctuations along the coast of Norway indicates a time-transgressive trend in the maximum extent of the LIA, which apparently seems to have occurred progressively later as we move northwards. We suggest that this trend is likely due to regional winter precipitation differences along the coast of Norway.

A new model for global glacier change and sea-level rise

The anticipated retreat of glaciers around the globe will pose far-reaching challenges to the management of fresh water resources and significantly contribute to sea-level rise within the coming decades. Here, we present a new model for calculating the 21st century mass changes of all glaciers on Earth outside the ice sheets. The Global Glacier Evolution Model (GloGEM) includes mass loss due to frontal ablation at marine-terminating glacier fronts and accounts for glacier advance/retreat and surface Elevation changes. Simulations are driven with monthly near-surface air temperature and precipitation from 14 Global Circulation Models forced by the RCP2.6, RCP4.5 and RCP8.5 emission scenarios. Depending on the scenario, the model yields a global glacier volume loss of 25-48% between 2010 and 2100. For calculating glacier contribution to sea-level rise, we account for ice located below sea-level presently displacing ocean water. This effect reduces glacier contribution by 11-14%, so that our model predicts a sea-level equivalent (multi-model mean +-1 standard deviation) of 79+-24 mm (RCP2.6), 108+-28 mm (RCP4.5) and 157+-31 mm (RCP8.5). Mass losses by frontal ablation account for 10% of total ablation globally, and up to 30% regionally. Regional equilibrium line altitudes are projected to rise by 100-800 m until 2100, but the effect on ice wastage depends on initial glacier hypsometries.

Reconstruction and regional significance of the coire breac palaeoglacier, glen esk, eastern grampian highlands, scotland

We present a glaciological and climatic reconstruction of a former glacier in Coire Breac, an isolated cirque within the Eastern Grampian plateau of Scotland, 5 km from the Highland edge. Published glacier reconstructions of presumed Younger Dryas‐age glaciers in this area show that equilibrium line altitudes decreased steeply towards the east coast, implying a arctic maritime glacial environment. Extrapolation of the ELA trend surface implies that glaciers should have existed in suitable locations on the plateau, a landscape little modified by glaciation. In Coire Breac, a 0.35 km2 cirque glacier existed with an equilibrium line altitude of 487 ± 15 m above present sea level. The equilibrium line altitude matches closely the extrapolated regional equilibrium line altitude trend surface for Younger Dryas Stadial glaciers. The mean glacier thickness of 24 m gives an ice volume of 7.8 × 106 m3, and a maximum basal shear stress of c. 100 kPa−1. Ablation gradient was c. –0.0055 m m−1, with a mean July temperature at the equilibrium line altitude of c. 5.1°C. The reconstruction implies an arctic maritime climate of low precipitation with local accumulation enhanced by blown snow, which may explain the absence of other contemporary glaciers nearby. Reconstructed ice flow lines show zones of flow concentration around the lower ice margin which help to explain the distribution of depositional facies associated with a former debris cover which may have delayed eventual glacier retreat. No moraines in the area have been dated, so palaeoclimatic interpretations remain provisional, and a pre‐Lateglacial Interstadial age cannot be ruled out.

Reconstruction of Late Pleistocene glaciers on Mount Lovćen, Montenegro

The study of Pleistocene glaciations in the Dinaric Alps began at the end of the 19th century. Since then, the interpretation of glaciation in this area has not seen much study until the beginning of the 21st century, when a few modern geomorphological studies with numerical dating and equilibrium line altitude (ELA) reconstructions were carried out in the SE part of the Dinaric Alps (where Mount Lovćen is situated). The morphology of Mount Lovćen consists of a karst landscape modified by glacial processes. This paper presents geomorphological evidence of glaciations on Mount Lovćen, along with a reconstruction of paleoglacier dimensions and estimations of the paleo-equilibrium line altitude. Mapping of glacial landforms on Mount Lovćen revealed at least two different glacial events. Only limited patches of older glacial deposits were identified, and therefore the reconstruction of glacier dimensions and ELA were only performed for the younger glaciation. This glacier extended over an area of at least 39.9 km2 and reached a maximum thickness of about 290 m. The regional ELA for the Mount Lovćen glaciers was calculated as 1259 m.

Was Mount Taylor glaciated in the Late Pleistocene? An analysis based on field evidence and regional equilibrium line altitudes

Was Mount Taylor glaciated in the Late Pleistocene? An analysis based on field evidence and regional equilibrium line altitudes

The Younger Dryas in the English Lake District: reconciling geomorphological evidence with numerical model outputs

The geomorphology of the south‐western and central Lake District, England is used to reconstruct the mountain palaeoglaciology pertaining to the Lateglacial and Younger Dryas. Limitations to previous ice‐mass reconstructions and consequent palaeoclimatic inferences include: (i) the use of static (steady‐state) glacier reconstructions, (ii) the assumption of a single‐stage Younger Dryas advance, (iii) greatly varying ice‐volume estimates, (iv) inexplicable spatial variations in ELA (Equilibrium Line Altitude), and (v) a lack of robust extent chronology. Here we present geomorphological mapping based on aerial photography and the NextMap Britain Digital Elevation Model, checked by ground survey. Former glacier extents were inferred and ELAs were calculated using the Balance Ratio method of Osmaston. Independently, a time‐dependant 2‐D ice‐flow model was forced by a regional ELA history that was scaled to the GRIP record. This provided a dynamic reconstruction of a mountain ice field that allowed for non‐steady‐state glacier evolution. Fluctuations in climate during the Younger Dryas resulted in multiple glacial advance positions that show agreement with the location of mapped moraines, and may further explain some of the ELA variations found in previous local and static reconstructions. Modelling based on the GRIP record predicts three phases: an initial maximum extent, a middle minor advance or stillstand, and a pronounced but less extensive final advance. The comparisons find that the reconstructions derived from geomorphological evidence are effective representations of steady‐state glacier geometries, but we do propose different extents for some glaciers and, in particular, a large former glacier in Upper Eskdale.

An approach to derive regional snow lines and glacier mass change from MODIS imagery, western North America

Abstract. We describe a method to calculate regional snow line elevations and annual equilibrium line altitudes (ELAs) from daily MODIS imagery (MOD02QKM) on large glaciers and icefields in western North America. An automated cluster analysis of the cloud-masked visible and near-infrared bands at 250 m resolution is used to delineate glacier facies (snow and ice) for ten glacierized regions between 2000–2011. For each region and season, the maximum observed value of the 20th percentile of snow-covered pixels (ZS(20)) is used to define a regional ELA proxy (ELAest). Our results indicate significant increases in the regional ELA proxy at two continental sites (Peyto Glacier and Gulkana Glacier) over the period of observation, though no statistically significant trends are identified at other sites. To evaluate the utility of regional ELA proxies derived from MOD02QKM imagery, we compare standard geodetic estimates of glacier mass change with estimates derived from historical mass balance gradients and observations of ZS(20) at three large icefields. Our approach yields estimates of mass change that more negative than traditional geodetic approaches, though MODIS-derived estimates are within the margins of error at all three sites. Both estimates of glacier mass change corroborate the continued mass loss of glaciers in western North America. Between 2000 and 2009, the geodetic change approach yields mean annual rates of surface elevation change for the Columbia, Lillooet, and Sittakanay icefields of −0.29 ± 0.05, −0.26 ± 0.05, and −0.63 ± 0.17 m a−1, respectively. This study provides a new technique for glacier facies detection at daily timescales, and contributes to the development of regional estimates of glacier mass change, both of which are critical for studies of glacier contributions to streamflow and global sea level rise.

Palaeoglacial and palaeoclimatic conditions in the NW Pacific, as revealed by a morphometric analysis of cirques upon the Kamchatka Peninsula

The distribution of glacial cirques upon the Kamchatka peninsula, Far Eastern Russia, is systematically mapped from satellite images and digital elevation model data. A total of 3758 cirques are identified, 238 of which are occupied by active glaciers. The morphometry of the remaining 3520 cirques is analysed. These cirques are found to show a very strong N bias in their azimuth (orientation), likely resulting from aspect-related variations in insolation. The strength of this N bias is considered to indicate that former glaciation upon the peninsula was often ‘marginal’, and mainly of cirque-type, with peaks extending little above regional equilibrium-line altitudes. This is supported by the fact that S and SE-facing cirques are the highest in the dataset, suggesting that glacier-cover was rarely sufficient to allow S and SE-facing glaciers to develop at low altitudes. The strength of these azimuth-related variations in cirque altitude is thought to reflect comparatively cloud-free conditions during former periods of glaciation. It is suggested that these characteristics, of marginal glaciation and comparatively cloud-free conditions, reflect the region's former aridity, which was likely intensified at the global Last Glacial Maximum, and during earlier periods of ice advance, as a result of the development of negative pressure anomalies over the North Pacific (driven by the growth of the Laurentide Ice Sheet), combined with other factors, including an increase in the extent and duration of sea ice, a reduction in global sea levels, cooler sea surface temperatures, and the localised growth of mountain glaciers. There is published evidence to suggest extensive glaciation of the Kamchatka Peninsula at times during the Late Quaternary, yet the data presented here appear to suggest that such phases were comparatively short-lived, and that smaller cirque-type glaciers were generally more characteristic of the period.

The plateau glacier in the Sierra de Béjar (Iberian Central System) during its maximum extent. Reconstruction and chronology

A detailed study of the glacial morphology in the Sierra de Béjar (Iberian Central System) provided a data set of geomorphic indicators to reconstruct the paleoglaciers developed in this mountain area during the last glacial cycle (Late Pleistocene). Applying a physical-based glacier model and using the geomorphic indicators, a three dimensional reconstruction of the ice mass during the maximum extent of the glaciers has been carried out. We used this reconstruction to project hypsometric curves over the former glaciers and to estimate the ELAs (Equilibrium Line Altitudes) of the paleoglaciers for their stage of maximum extent. At this stage the Sierra de Béjar hosted a plateau glacier, considered as a dome-shaped icecap around 57km2 in area. According to our estimations, the maximum thickness of the ice was 211m, the minimum elevation of paleoglaciers 1210m asl, and the regional ELA was at 2010m asl. During later stages, reduction in ice mass due to deglaciation caused the icecap to evolve into an icefield, and finally the main glacier was fragmented in valley and cirque glaciers. The geochronological data obtained with 10Be provides an age of ~27ka for the maximum extent of the glaciers (GM), whereas the global Last Glacial Maximum (LGM) represents a younger stage of the Sierra de Béjar glacier evolution. Finally, the new data obtained in the Sierra de Béjar allow evaluating the influence of some factors such as the continentality and latitudinal location, in the development of glacial processes in these areas of the Iberian Central System.

Cirques, peaks, and precipitation patterns in the Swiss Alps: Connections among climate, glacial erosion, and topography

Glacial erosion, a process influenced by climate, has been implicated in limiting the relief of mountain ranges. However, climate itself is sensitive to large-scale topography, suggesting that climate, topography, and glacial processes form a coupled system. Large spatial gradients in precipitation exist in the southern Swiss Alps, allowing us to study this coupling in a region of climate variability. More than 500 cirques were identified, and neighboring peaks were found to co-vary in elevation with cirque floors ( R 2 0.64). Cirque headwall relief does not vary with precipitation or cirque floor altitude. These relationships confirm the hypothesis that cirque formation restricts peak altitudes via slope processes that limit the relief of cirque headwalls. We compared the position of the regional equilibrium line altitude (ELA) estimated from modern climate to a surface defined by the cirque floors. The modern ELA and cirque floor surfaces are similar in shape, illustrating the impact of spatial variability in precipitation on glacial processes and topography. Precipitation variability in this region is partially dictated by large-scale topography. Therefore, precipitation patterns, glacial process domains, and topography must evolve together.