Search

High-resolution topographic modeling has become more accessible due to the development of structure from motion (SfM)-image-matching algorithms in digital photogrammetry. Large archival databases of historical aerial photographs are available in university, public, and government libraries, commonly as paper copies. The photographs can be in poor condition (i.e., deformed by humidity, scratched, or annotated). In addition, the negatives, as well as metadata, may be missing. Processing such photographs using classic stereo-photogrammetry is difficult and, in many instances, impossible. SfM can be applied to these photosets to access the valuable archive of geomorphic changes over the past century. In this paper, we illustrate the utility of the SfM technique using 568 digitized vertical aerial photographs of Mount Meager volcano, located in southwestern British Columbia, Canada. We use the aerial photographs, which span the period from 1947 to 2006, to track glaciers and glacier–landslide interactions on the volcano. Over this period, glaciers have thinned and retreated, interrupted by minor advances in the 1960s and 1970s. Landslides are frequent on the volcano and contribute to debris cover on the glaciers affecting the ablation process. SfM processing of the aerial photographs allowed us to unlock geomorphic information and reconstruct landscape change that would otherwise have been impossible. The results from SfM provide a visually effective way of presenting landscape change to a broad public audience, as a form of virtual geoheritage. The approach can thus be broadly applied in scientific and professional practices for improving land planning and hazard management.

Several global datasets of glacier thickness exist, but the number of observations from western Canada are sparse and spatially biased. To supplement these limited observations, we measured ice thickness with ice penetrating radar on five glaciers in the Columbia Mountains, Canada. Our radar surveys, when combined with previous surveys for two glaciers in the Rocky Mountains, total 182 km of transects that represent 34 672 point measurements of ice thickness. Our measurements are, on average, 38% thicker than previous surface inversion model estimates of glacier thickness. Using our measurements within a cross-validation scheme, we model ice thickness with the Open Global Glacier Model (OGGM) driven with recent observations of surface mass balance and glacier elevation. We calibrated OGGM ice thickness by optimizing the ice creep parameter in the model. The optimized OGGM yields an ice volume for Columbia Basin of 122.5 ± 22.4 km3for the year 2000, which is 23% greater than the range of previous estimates. At current rates of glacier mass loss for this region, glaciers would disappear from the basin in about 65–80 years. Disappearance of these glaciers will negatively affect the basin's surface hydrology, freshwater availability and aquatic ecosystems.

The combined use of unmanned aerial vehicles (UAVs) and structure-from-motion (SfM) is rapidly growing as a cost-effective alternative to airborne laser scanning (lidar) for reconstructing glacier surfaces. Here we present a thorough analysis of the precision and accuracy of a photogrammetric point cloud (PPC) constructed through SfM from UAV-acquired imagery over the spring snow surface at Haig Glacier, Alberta, Canada, the first of its kind in a glaciological setting. An aerial lidar survey conducted concurrently with UAV surveys was used to examine spatial patterns in the PPC accuracy. We found a median error in the PPC of −0.046 ± 0.067 m, with a 95% quantile of 0.218 m. Mean precision of the PPC was 0.199 m, with large spatially clustered outliers. We found an association between high-error, low-precision, and high-surface roughness in the PPC, likely due to illumination characteristics of the snow surface. Glacier surface reconstructions are important for geodetic mass balance measurements, giving key insights into changing climate where in situ measurements are difficult to obtain. The PPC errors are small enough that they would have minimal effects on total mass balance, should the technique be applied across the glacier.

Some lateral moraines contain a rich record of Holocene glacial expansion. Previous workers have used such evidence to document glacial fluctuations in western Canada, Alaska, and the U.S. Pacific Northwest, but similar studies in Patagonia are uncommon. Here we report on the late Holocene behavior of Stoppani Glacier, a 75 km2 glacier sourced in the Cordillera Darwin, southernmost Patagonia. Based on radiocarbon-dated wood and organic material contained in the glacier's northeast lateral moraine, we infer that Stoppani Glacier advanced shortly before 3.8–3.6, at 3.2–2.8, 2.3–2.1, and 0.3–0.2, and possibly sometime before 1.4–1.3 and 0.8–0.7 cal ka BP. These advances culminated at 0.3–0.2 cal ka BP, when the glacier constructed a prominent end moraine, marking its greatest extent of the past 4000 years. Although the timing of several of the advances overlap with the age range of glacial expansion recognized elsewhere in Patagonia, some do not. Asynchronous behavior observed in the glacial record may arise from the type of evidence (e.g., lateral stratigraphy vs. end moraine) used to document glacial fluctuations or variations in climate or glacial response times. A significant difference between the Stoppani record and some other Patagonian records is that the former indicates general expansion of ice over the last 4000 years, whereas the latter indicate a net decrease in extent over that period.

Two catastrophic landslides occurred in quick succession on 13 and 16 May 2019, from the north face of Joffre Peak, Cerise Creek, southern Coast Mountains, British Columbia. With headscarps at 2560 m and 2690 m elevation, both began as rock avalanches, rapidly transforming into debris flows along middle Cerise Creek, and finally into debris floods affecting the fan. Beyond the fan margin, a flood surge on Cayoosh Creek reached bankfull and attenuated rapidly downstream; only fine sediment reached Duffey Lake. The toe of the main debris flow deposit reached 4 km from the headscarp, with a travel angle of 0.28, while the debris flood phase reached the fan margin 5.9 km downstream, with a travel angle of 0.22. Photogrammetry indicates the source volume of each event is 2–3 Mm3, with combined volume of 5 Mm3. Lidar differencing, used to assess deposit volume, yielded a similar total result, although error in the depth estimate introduced large volume error masking the expected increase due to dilation and entrainment. The average velocity of the rock avalanche-debris flow phases, from seismic analysis, was ~ 25–30 m/s, and the velocity of the 16 May debris flood on the upper fan, from super-elevation and boulder sizes, was 5–10 m/s. The volume of debris deposited on the fan was ~ 104 m3, 2 orders of magnitude less than the avalanche/debris flow phases. Progressive glacier retreat and permafrost degradation were likely the conditioning factors; precursor rockfall activity was noted at least ~6 months previous; thus, the mountain was primed to fail. The 13 May landslide was apparently triggered by rapid snowmelt, with debuttressing triggering the 16 May event.

Finely resolved geodetic data provide an opportunity to assess the extent and morphology of crevasses and their change over time. Crevasses have the potential to bias geodetic measurements of elevation and mass change unless they are properly accounted for. We developed a framework that automatically maps and extracts crevasse geometry and masks them where they interfere with surface mass-balance assessment. Our study examines airborne light detection and ranging digital elevation models (LiDAR DEMs) from Haig Glacier, which is experiencing a transient response in its crevassed upper regions as the glacier thins, using a self-organizing map algorithm. This method successfully extracts and characterizes ~1000 crevasses, with an overall accuracy of 94%. The resulting map provides insight into stress and flow conditions. The crevasse mask also enables refined geodetic estimates of summer mass balance. From differencing of September and April LiDAR DEMs, the raw LiDAR DEM gives a 9% overestimate in the magnitude of glacier thinning over the summer: −5.48 m compared with a mean elevation change of −5.02 m when crevasses are masked out. Without identification and removal of crevasses, the LiDAR-derived summer mass balance therefore has a negative bias relative to the glaciological surface mass balance.

Seasonal measurements of glacier mass balance provide insight into the relation between climate forcing and glacier change. To evaluate the feasibility of using remotely sensed methods to assess seasonal balance, we completed tandem airborne laser scanning (ALS) surveys and field-based glaciological measurements over a 4-year period for six alpine glaciers that lie in the Columbia and Rocky Mountains, near the headwaters of the Columbia River, British Columbia, Canada. We calculated annual geodetic balance using coregistered late summer digital elevation models (DEMs) and distributed estimates of density based on surface classification of ice, snow, and firn surfaces. Winter balance was derived using coregistered late summer and spring DEMs, as well as density measurements from regional snow survey observations and our glaciological measurements. Geodetic summer balance was calculated as the difference between winter and annual balance. Winter mass balance from our glaciological observations averaged 1.95±0.09 m w.e. (meter water equivalent), 4 % larger than those derived from geodetic surveys. Average glaciological summer and annual balance were 3 % smaller and 3 % larger, respectively, than our geodetic estimates. We find that distributing snow, firn, and ice density based on surface classification has a greater influence on geodetic annual mass change than the density values themselves. Our results demonstrate that accurate assessments of seasonal mass change can be produced using ALS over a series of glaciers spanning several mountain ranges. Such agreement over multiple seasons, years, and glaciers demonstrates the ability of high-resolution geodetic methods to increase the number of glaciers where seasonal mass balance can be reliably estimated.

The mass-balance—elevation relation for a given glacier is required for many numerical models of ice flow. Direct measurements of this relation using remotely-sensed methods are complicated by ice dynamics, so observations are currently limited to glaciers where surface mass-balance measurements are routinely made. We test the viability of using the continuity equation to estimate annual surface mass balance between flux-gates in the absence of in situ measurements, on five glaciers in the Columbia Mountains of British Columbia, Canada. Repeat airborne laser scanning surveys of glacier surface elevation, ice penetrating radar surveys and publicly available maps of ice thickness are used to estimate changes in surface elevation and ice flux. We evaluate this approach by comparing modeled to observed mass balance. Modeled mass-balance gradients well-approximate those obtained from direct measurement of surface mass balance, with a mean difference of +6.6 ± 4%. Regressing modeled mass balance, equilibrium line altitudes are on average 15 m higher than satellite-observations of the transient snow line. Estimates of mass balance over flux bins compare less favorably than the gradients. Average mean error (+0.03 ± 0.07 m w.e.) between observed and modeled mass balance over flux bins is relatively small, yet mean absolute error (0.55 ± 0.18 m w.e.) and average modeled mass-balance uncertainty (0.57 m w.e.) are large. Mass conservation, assessed with glaciological data, is respected (when estimates are within 1σ uncertainties) for 84% of flux bins representing 86% of total glacier area. Uncertainty on ice velocity, especially for areas where surface velocity is low (−1) contributes the greatest error in estimating ice flux. We find that using modeled ice thicknesses yields comparable modeled mass-balance gradients relative to using observations of ice thickness, but we caution against over-interpreting individual flux-bin mass balances due to large mass-balance residuals. Given the performance of modeled ice thickness and the increasing availability of ice velocity and surface topography data, we suggest that similar efforts to produce mass-balance gradients using modern high-resolution datasets are feasible on larger scales.

Knowledge about the long-term response of High Mountain Asian glaciers to climatic variations is paramount because of their important role in sustaining Asian river flow. Here, a satellite-based time series of glacier mass balance for seven climatically different regions across High Mountain Asia since the 1960s shows that glacier mass loss rates have persistently increased at most sites. Regional glacier mass budgets ranged from −0.40 ± 0.07 m w.e.a−1 in Central and Northern Tien Shan to −0.06 ± 0.07 m w.e.a−1 in Eastern Pamir, with considerable temporal and spatial variability. Highest rates of mass loss occurred in Central Himalaya and Northern Tien Shan after 2015 and even in regions where glaciers were previously in balance with climate, such as Eastern Pamir, mass losses prevailed in recent years. An increase in summer temperature explains the long-term trend in mass loss and now appears to drive mass loss even in regions formerly sensitive to both temperature and precipitation., Multi-platform satellite observations document six decades of glacier mass balance variability across High Mountain Asia (HMA). Heterogeneous rates of ice loss reflect regional climatic differences, but ice loss is now pervasive across HMA even in regions formerly exhibiting slight mass gains.

Glaciers distinct from the Greenland and Antarctic ice sheets are shrinking rapidly, altering regional hydrology1, raising global sea level2 and elevating natural hazards3. Yet, owing to the scarcity of constrained mass loss observations, glacier evolution during the satellite era is known only partially, as a geographic and temporal patchwork4,5. Here we reveal the accelerated, albeit contrasting, patterns of glacier mass loss during the early twenty-first century. Using largely untapped satellite archives, we chart surface elevation changes at a high spatiotemporal resolution over all of Earth’s glaciers. We extensively validate our estimates against independent, high-precision measurements and present a globally complete and consistent estimate of glacier mass change. We show that during 2000–2019, glaciers lost a mass of 267 ± 16 gigatonnes per year, equivalent to 21 ± 3 per cent of the observed sea-level rise6. We identify a mass loss acceleration of 48 ± 16 gigatonnes per year per decade, explaining 6 to 19 per cent of the observed acceleration of sea-level rise. Particularly, thinning rates of glaciers outside ice sheet peripheries doubled over the past two decades. Glaciers currently lose more mass, and at similar or larger acceleration rates, than the Greenland or Antarctic ice sheets taken separately7–9. By uncovering the patterns of mass change in many regions, we find contrasting glacier fluctuations that agree with the decadal variability in precipitation and temperature. These include a North Atlantic anomaly of decelerated mass loss, a strongly accelerated loss from northwestern American glaciers, and the apparent end of the Karakoram anomaly of mass gain10. We anticipate our highly resolved estimates to advance the understanding of drivers that govern the distribution of glacier change, and to extend our capabilities of predicting these changes at all scales. Predictions robustly benchmarked against observations are critically needed to design adaptive policies for the local- and regional-scale management of water resources and cryospheric risks, as well as for the global-scale mitigation of sea-level rise. Analysis of satellite stereo imagery uncovers two decades of mass change for all of Earth’s glaciers, revealing accelerated glacier shrinkage and regionally contrasting changes consistent with decadal climate variability.

Glacier mass balance models are needed at sites with scarce long-term observations to reconstruct past glacier mass balance and assess its sensitivity to future climate change. In this study North American Regional Reanalysis (NARR) data are used to force a physically-based, distributed glacier mass balance model of Saskatchewan Glacier for the historical period 1979–2016 and assess it sensitivity to climate change. A two-year record (2014–2016) from an on-glacier automatic weather station (AWS) and a homogenized historical precipitation record from nearby permanent weather stations were used to downscale air temperature, relative humidity, wind speed, incoming solar radiation and precipitation from the nearest NARR gridpoint to the glacier AWS site. The model was run with fixed (1979, 2010) and time-varying (dynamic) geometry using a multi-temporal digital elevation model (DEM) dataset. The model showed a good performance against recent (2012–2016) direct glaciological mass balance observations as well as with cumulative geodetic mass balance estimates. The simulated mass balance showed a large sensitivity to the biases in NARR precipitation and solar radiation, as well as to the prescribed precipitation lapse rate and ice aerodynamic roughness lengths, showing the importance of constraining these parameters with ancillary data. The difference between the static (1979) and dynamic simulations showed small differences (mean = 0.06 m w.e. a−1 or 1.5 m w.e. over 37 yrs), indicating minor effects of elevation changes on the glacier specific mass balance. The static mass balance sensitivity to climate was assessed for prescribed changes in regional mean air temperature between 0 to 7 °C and precipitation between −20 to +20 %, which comprise the spread of ensemble IPCC representative concentration pathways climate scenarios for the mid (2041–2070) and late (2071–2100) 21st century. The climate sensitivity experiments showed that future changes in precipitation would have a small impact on glacier mass-balance, while the temperature sensitivity increases with warming, from −0.65 to −0.93 m w.e. °C−1. Increased melting accounted for 90 % of the temperature sensitivity while precipitation phase feedbacks accounted for only 10 %. Roughly half of the melt response to warming was driven by a positive albedo feedback, in which glacier albedo decreases as the snow cover on the glacier thins and recedes earlier in response to warming, increasing net solar radiation fluxes. About one quarter of the melt response to warming was driven by latent heat energy gains (positive humidity feedback). Our study underlines the key role of albedo and air humidity in modulating the response of winter-accumulation type mountain glaciers and upland icefield-outlet glacier settings to climate.

In 2015, the Olympic Mountains contain 255 glaciers and perennial snowfields totaling 25.34 ± 0.27 km2, half of the area in 1900, and about 0.75 ± 0.19 km

Mount Meager, a glacier-clad volcanic complex in British Columbia, Canada, is known for its large landslides, as well as a major eruption about 2360 years ago. In 2010, after decades of glacier retreat, the south flank of Mount Meager collapsed, generating a huge (53 Mm3) landslide. In 2016, fumaroles formed ice caves in one of the glaciers on the complex. This glacier is bordered by a large unstable slope presently moving about 3.5 cm per month. If this slope were to fail, a long-runout debris avalanche would reach the floor of the Lillooet River valley, with possible destructive effects on downstream infrastructure. The unloading of the volcanic edifice from an abrupt failure would also have unknown effects on the magmatic plumbing system. From geochemical, geophysical, and petrological data, we infer the presence of a magmatic chamber at 3–16 km depth. Based on numerical model simulations carried out to constrain the stress change, the failure would affect the stress field to depths of up to ~6 km, with changes in effective stress of up to ~4 MPa. The change in effective stress following such a landslide might destabilize the magmatic chamber and trigger an eruption. This result also suggests that a previously documented major flank collapse may have had a role in the 2360 cal yr BP eruption.

Glaciers provide cold, turbid runoff to many mountain streams in the late summer and buffer against years with low snowfall. The input of glacial meltwater to streams maintains unique habitats and support a diversity of stream flora and fauna. In western Canada, glaciers are anticipated to retreat by 60–80% by the end of the century, and this retreat will invoke widespread changes in mountain ecosystems. We used a space-for-time substitution along a gradient of glacierization in western Canada to develop insights into changes that may occur in glaciated regions over the coming decades. Here we report on observed changes in physical (temperature, turbidity), and chemical (dissolved and total nutrients) characteristics of mountain streams and the associated shifts in their diatom communities during de-glacierization. Shifts in habitat characteristics across gradients include changes in nutrient concentrations, light penetration, temperatures, and flow, all of which have led to distinct changes in diatom community composition. Importantly, glacial-fed rivers were 3–5 °C cooler than rivers without glacial contributions. Declines in glacial meltwater contribution to streams resulted in shifts in the timing of nutrient fluxes and lower concentrations of total phosphorus (TP), soluble reactive phosphorus (SRP), and higher dissolved inorganic nitrogen (DIN) and light penetration. The above set of conditions were linked to the overgrowth of the benthic diatom Didymosphenia geminata. These changes in stream condition and D. geminata colony development primarily occurred in streams with marginal (2–5%) to no glacier cover. Our data support a hypothesis that climate-induced changes in river hydrochemistry and physical condition lead to a phenological mismatch that favors D. geminata bloom development.

The world’s glaciers distinct from the Greenland and Antarctic ice sheets are shrinking rapidly, altering regional hydrology and raising global sea level. Yet, due to the scarcity of globally consistent observations, their recent evolution is only known as a heterogeneous temporal and geographic patchwork and future projections are thus not optimally constrained.Here, we present the first globally complete, consistent and resolved estimate of glacier mass change derived from more than half a million digital elevation models (DEMs) generated or extracted from multiple satellite archives including ASTER, ArcticDEM and REMA. Combining state-of-the-art numerical photogrammetry and novel statistical approaches, we reconstruct two decades of glacier surface elevation change at an unprecedented spatial and temporal resolution. We validate our results by comparing them to independent, high-precision elevation measurements from the ICESat and IceBridge campaigns, as well as to very high resolution DEM differences from LiDAR, Pléiades, and SPOT-6. The elevation time series are integrated to volume changes for every single glacier on Earth and, by assuming an average density, aggregated to regional and global mass changes. We compare our revised glacier mass changes to earlier estimates derived from altimetry, gravimetry, geodetic and field data. As an illustration, our integrated geodetic mass loss over all Icelandic glaciers yields -8.3 +- 1.1 Gt yr-1 over the period 2002-2016 in agreement with a recent gravimetry estimate of -8.3 +- 1.8 Gt yr-1 (Wouters et al., 2019), known to perform well in this region. Both estimates are more negative than -5.7 +- 1.2 Gt yr-1, compiled from glaciological observations and geodetic data (Zemp et al., 2019).Our global estimate of glacier mass change constitutes a new benchmark dataset that will help to: (i) assess present-day and future climate change impacts on glaciers; (ii) close the sea-level rise budget; (iii) assess the threat on water resources and (iv) facilitate research on natural hazards related to glaciers. Our results specifically provide a strong observational basis that holds a great potential to further our understanding of the multi-scale morphologic and climatic drivers of glacier mass change, essential to improve physically-based glaciological modelling and calibrate future projections.

The timing of Cordilleran Ice Sheet deglaciation along the central coast of British Columbia, Canada, informs climate forcing and early human migration. Thirty‐two 10Be exposure ages from glacial erratics and local bedrock on the western margin of the former ice sheet represent the earliest exposure ages for the Cordilleran Ice Sheet during the last deglaciation. These data show the western ice margin was retreating by 18.1 ± 0.2 ka, consistent with the global record of ice mass loss. In contrast, parts of the southern margin reached a maximum at ca. 17.6‐16.6 ka, and our data demonstrate a diachronous response of the Cordilleran Ice Sheet during deglaciation. We also show that a low altitude site was exposed by at least 17.7 ± 0.3 ka, implying that numerous ice‐free areas existed along the coastal margin by this time, providing a viable route for the first humans entering the Americas.

Combined use of radiocarbon-dated subfossil wood within lateral moraines and surface exposure ages on moraine boulders provides an approach to better constrain times of glacier advance and onset of retreat. We test this method at Gilbert Glacier in the southern Coast Mountains of British Columbia where units of sediments associated with glacier expansion date to 4.8–4.6, 4.5–4.3, 4.0–3.9, 3.8–3.6, 3.4, 3.2–2.9, 2.7, and 0.5–0.3 kilo calendar years BP (ka; 2-sigma age range). Surface exposure (10Be) ages reveal times of moraine stabilization at 1.83–1.78, 1.38–1.28, 0.85–0.76, and 0.13–0.06 ka (interquartile range). Analysis of both datasets, as well as previously published regional advance records, narrows the age range of four late Holocene advances to 2.0–1.8, ∼1.5–1.3, ∼0.9–0.8, and 0.4–0.1 ka. We advocate for widespread use of our tandem approach at other sites throughout Earth's high mountains to narrow the uncertainties associated with glacier expansion and better understand how glaciers respond to climate change.

Over the past decade, western North America glaciers experienced strong mass loss. Regional mass loss during the ablation season is influenced by air temperature, but the importance of other factors such as changes in surface albedo remains uncertain. We examine changes in surface albedo for 17 glaciated regions of western North America as documented in a 20-year record (2000 to 2019) of MODIS daily snow albedo (MOD10A1). Trend analysis reveals that albedo declined for 4% to 81% of the albedo grid cells, and the largest negative trends were situated south of 60°N and in the provinces of British Columbia and Alberta. Sen's slope estimates indicate that 15 of 17 regions showed a decline of which the majority of the largest declines occurred within 100 m of glacier median elevation, suggesting that these declines are driven by a rise of the transient snowline. For most regions, albedo correlates strongly to temperature, and albedo trend in the Chugach region of Alaska, the South Coast, Southern Interior and Central and Southern Rockies of Canada show a significant relationship to aerosols optical depth. Temperature is approximately 2–6 times more predictive of the variation in albedo than AOD for the majority of regions, except the Southern Interior and Southern Rockies where albedo shows a greater dependence on AOD. Investigation of broadband albedo (MCD43A3) for snow grid cells above glacier median elevation in the Central and Southern Rockies shows that declines in the visible and near infrared portions of the spectrum are linked to the presence of forest fire generated aerosols. The results of this study indicate that glacier surface mass balance experiences a regional dependence on forest fire generated light absorbing particles.

We use annually laminated lake sediments to estimate suspended sediment yield for a 3850 km2 montane catchment in the British Columbia Coast Mountains. Sediment yield over the past 369 years averages 213 ± 38 Mg·km−2·a−1. Sediment yield increases to 285 ± 50 Mg·km−2·a−1 during the first half of the twentieth century and declines thereafter. The frequency of high-yield events during the 369 year period is irregular: 11 of the 34 events occur in the early part of the twentieth century, a time when glaciers in the watershed underwent major retreat. We fitted a generalized extreme value (GEV) model to estimate quantiles of the sediment yield distribution, and we used epoch analysis to examine persistence in sediment yield following 34 of the largest events. Persistence is greatest for the most extreme events; it is more variable for events that recur, on average, every 10–25 years. Our results indicate that sediment yield is linked to long-term changes in sediment supply to the lake. The results of this study extend earlier sediment yield estimates and improve understanding of linkages to watershed geomorphology, recent glacier retreat, and landslides in the Lillooet River watershed.

As part of surface energy balance models used to simulate glacier melting, choosing parameterizations to adequately estimate turbulent heat fluxes is extremely challenging. This study aims to evaluate a set of four aerodynamic bulk methods (labeled as C methods), commonly used to estimate turbulent heat fluxes for a sloped glacier surface, and two less commonly used bulk methods developed from katabatic flow models. The C methods differ in their parameterizations of the bulk exchange coefficient that relates the fluxes to the near-surface measurements of mean wind speed, air temperature, and humidity. The methods' performance in simulating 30 min sensible- and latent-heat fluxes is evaluated against the measured fluxes from an open-path eddy-covariance (OPEC) method. The evaluation is performed at a point scale of a mountain glacier, using one-level meteorological and OPEC observations from multi-day periods in the 2010 and 2012 summer seasons. The analysis of the two independent seasons yielded the same key findings, which include the following: first, the bulk method, with or without the commonly used Monin–Obukhov (M–O) stability functions, overestimates the turbulent heat fluxes over the observational period, mainly due to a substantial overestimation of the friction velocity. This overestimation is most pronounced during the katabatic flow conditions, corroborating the previous findings that the M–O theory works poorly in the presence of a low wind speed maximum. Second, the method based on a katabatic flow model (labeled as the KInt method) outperforms any C method in simulating the friction velocity; however, the C methods outperform the KInt method in simulating the sensible-heat fluxes. Third, the best overall performance is given by a hybrid method, which combines the KInt approach with the C method; i.e., it parameterizes eddy viscosity differently than eddy diffusivity. An error analysis reveals that the uncertainties in the measured meteorological variables and the roughness lengths produce errors in the modeled fluxes that are smaller than the differences between the modeled and observed fluxes. This implies that further advances will require improvement to model theory rather than better measurements of input variables. Further data from different glaciers are needed to investigate any universality of these findings.

We analyzed a peat-bog sequence from the central region of Tierra del Fuego, southernmost Patagonia, to reconstruct the environmental changes over the past 15,000 years. Postglacial vegetation was mainly composed by grasslands and shrubby communities with sparsely distributed Nothofagus trees under dry conditions. The predominance of Glomus sp. and Gelasinospora sp. may also indicate less humid conditions related to the steppe communities. The presence of herbivorous grazers is suggested by the record of ascospores of coprophilous Sordaria-type and Sporormiella-type throughout the profile. Between 11,200 and 6500 cal yr BP, the Nothofagus woodland established under warmer and wetter climate than before. However, precipitation must have remained below present-day levels. Most fungal remains (Gaeumannomyces sp., Glomus sp., Gelasinospora sp., Microthyriaceae, Spegazzinia tessarthra, Alternaria sp.) are likely associated with the development of forest-steppe ecotone communities. Following 6500 cal yr BP, when climate conditions were more humid and colder than before, the closed-canopy forest dominated the landscape, and Microthyriaceae remains prevailed among fungi. During the last 1000 years, the record showed open forest communities accompanied by the presence of Glomus sp. These palaeoenvironmental changes observed along the sequence suggested variations in the amount of precipitation of westerly origin related to shifts in the position and/or strength of the southern margin of the westerlies. This study revealed that fungal remains complement pollen/spores analysis by providing important independent information on the palaeoenvironmental and palaeoclimate conditions that prevailed during the late Pleistocene-Holocene in the Lago Fagnano area. Fil: Musotto, Lorena Laura. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Bahía Blanca. Instituto Geológico del Sur. Universidad Nacional del Sur. Departamento de Geología. Instituto Geológico del Sur; Argentina Fil: Borromei, Ana Maria. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Bahía Blanca. Instituto Geológico del Sur. Universidad Nacional del Sur. Departamento de Geología. Instituto Geológico del Sur; Argentina Fil: Bianchinotti, Maria Virginia. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Bahía Blanca. Centro de Recursos Naturales Renovables de la Zona Semiárida. Universidad Nacional del Sur. Centro de Recursos Naturales Renovables de la Zona Semiárida; Argentina Fil: Coronato, Andrea Maria Josefa. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Austral de Investigaciones Científicas; Argentina Fil: Menounos, Brian. University of Northern British Columbia; Canadá Fil: Osborn, Gerald. University of Calgary; Canadá Fil: Marr, Robert. University of Calgary; Canadá

Alpine glaciers and perennial snow fields are important hydrologic elements in many mountain environments providing runoff during the late summer and during periods of drought. Because relatively long records of glacier mass-balance data are absent from many glacierized catchments, it remains unclear to what extent shrinking perennial snow and glaciers have affected runoff trends from these watersheds. Here we employ a hydrograph separation technique that uses a double mass curve in an attempt to isolate changes in runoff due to glacier retreat and disappearance of perennial snow. The method is tested using hydrometric data from 20 glacierized and 16 non-glacierized catchments in the Columbia Basin of Canada. The resulting estimates on cryosphere storage contribution to streamflow were well correlated to other regional estimates based on measurements as well as empirical and mechanistic models. Annual cryosphere runoff changed from +19 to -55% during the period 1975-2012, with an average decline of 26%. For August runoff, these changes ranged from +17 to -66%, with an average decrease of 24%. Reduction of cryosphere contributions to annual and late summer flows are expected to continue in coming decades as glaciers and the perennial snow patches shrink. Our method to isolate changes in late summer cryospheric storage contributions can be used as a first order estimate on changes in glacier contributions to flow and may help researchers and water managers target watersheds for further analysis.

The Yukon–Alaska Highway corridor in southern Yukon is subject to geohazards ranging from landslides to floods and earthquakes on faults in the St. Elias Mountains and Shakwak Valley. Here we discuss the late Holocene seismic history of the Denali fault, located at the eastern front of the St. Elias Mountains and one of only a few known seismically active terrestrial faults in Canada. Holocene faulting is indicated by scarps and mounds on late Pleistocene drift and by tectonically deformed Pleistocene and Holocene sediments. Previous work on trenches excavated against the fault scarp near the Duke River reveals paleoseismic sediment disturbance dated to ∼300–1,200, 1,200–1,900, and 3,000 years ago. Re-excavation of the trenches indicates a fourth event dated to 6,000 years ago. The trenches are interpreted to show a negative flower structure produced by extension of sediments by dextral strike-slip fault movement. Nearby Crescent Lake is ponded against the fault scarp. Sediment cores reveal four abrupt sediment and diatom changes reflecting seismic shaking at ∼1,200–1,900, 1,900–5,900, 5,900–6,200, and 6,500–6,800 years ago. At the Duke River, the fault offsets sediments, including two White River tephra layers (∼1,900 and 1,200 years old). Late Pleistocene outwash gravel and overlying Holocene aeolian sediments show in cross section a positive flower structure indicative of post-glacial contraction of the sediments by dextral strike-slip movement. Based on the number of events reflecting ∼M6, we estimate the average recurrence of large earthquakes on the Yukon part of the Denali fault to be about 1,300 years in the past 6,500–6,800 years.

Near the Pleistocene Termination, a glacier-dammed lake in central British Columbia suddenly drained to the south along the Fraser River valley. Floodwater travelled 330 km down the valley to Hope, British Columbia, and from there to the west into the Salish Sea near Vancouver. The flood was caused by the failure of an ice dam formed by the terminus of glaciers flowing from the central Coast Mountains across the British Columbia Interior Plateau. The ice dam impounded several hundred cubic kilometres of water to a maximum elevation of about 810 m asl (above sea level); at its maximum, the lake at the ice dam was over 250 m deep. Geomorphic and sedimentary evidence for the flood includes streamlined boulder-strewn bars, gravel dune fields, and terraces sloping up Fraser and lowermost Thompson valleys, opposite the present direction of river flow. The gravel bars and flood terraces are underlain by sheets of massive to poorly sorted gravel containing large boulders and rip-up clasts of silt and till. Shortly after the flood, a landslide near the northern margin of the former glacier dam impounded water to an elevation of about 550 m asl. This lake emptied due to overflow and incision of the landslide dam. The outburst flood from glacial Lake Fraser and the subsequent draining of the landslide-dammed lake deeply incised the older sediment fill in Fraser Valley and transported much of this sediment into the proto-Salish Sea west of Vancouver, British Columbia and Bellingham, Washington. TCN ages on flood-transported boulders at three localities along the flood path agree with radiocarbon ages on inferred flood layers in ODP cores collected from Saanich Inlet, a fiord on southern Vancouver Island, 80 km south-southwest of Vancouver.

International audience

We applied photogrammetric methods with aerial photography from 11 different years between 1946 and 2005 to assess changes in area and volume of 33 glaciers in the Cariboo Mountains of British Columbia for the latter half of the 20th century. These are used to identify changes in extent and elevation primarily for the periods 1952–1985, 1985–2005, and 1952–2005. All glaciers receded during the period 1952–2005; area retreat averaged −0.19 ± 0.05 % a−1. From 1952 to 1985, nine glaciers advanced; following 1985, retreat rates accelerated to −0.41 ± 0.12% a−1. Thinning rates of a subset of seven glaciers likewise accelerated, from −0.14 ± 0.04 m w.e. a−1 (1952–1985) to −0.50 ± 0.07 m w.e. a−1 for the period 1985–2005. Temperatures increased from the earlier to the latter period for the ablation (+0.38 °C) and accumulation (+0.87 °C) seasons, and average precipitation decreased, particularly in the accumulation season (−32 mm, −3.2%). Our comparison of surface area change with glacier morphometry corroborates previous studies that show primary relations between extent change and surface area. We also find that the strength and sign of these relations varied for different epochs. Our results also indicate that the 1985 glacier extent for the study area reported previously by other studies may be slightly overestimated due to errant mapping of late-lying snow cover.

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.

Glacier melt provides important contributions to streamflow in many mountainous regions. Hydrologic model calibration in glacier-fed catchments is difficult because errors in modelling snow accumulation can be offset by compensating errors in glacier melt. This problem is particularly severe in catchments with modest glacier cover, where goodness-of-fit statistics such as the Nash-Sutcliffe model efficiency may not be highly sensitive to the streamflow variance associated with glacier melt. While glacier mass balance measurements can be used to aid model calibration, they are absent for most catchments. We introduce the use of glacier volume change determined from repeated glacier mapping in a guided GLUE (generalized likelihood uncertainty estimation) procedure to calibrate a hydrologic model. This approach is applied to the Mica basin in the Canadian portion of the Columbia River Basin using the HBV-EC hydrologic model. Use of glacier volume change in the calibration procedure effectively reduced parameter uncertainty and helped to ensure that the model was accurately predicting glacier mass balance as well as streamflow. The seasonal and interannual variations in glacier melt contributions were assessed by running the calibrated model with historic glacier cover and also after converting all glacierized areas to alpine land cover in the model setup. Sensitivity of modelled streamflow to historic changes in glacier cover and to projected glacier changes for a climate warming scenario was assessed by comparing simulations using static glacier cover to simulations that accommodated dynamic changes in glacier area. Although glaciers in the Mica basin only cover 5% of the watershed, glacier ice melt contributes up to 25% and 35% of streamflow in August and September, respectively. The mean annual contribution of ice melt to total streamflow varied between 3 and 9% and averaged 6%. Glacier ice melt is particularly important during warm, dry summers following winters with low snow accumulation and early snowpack depletion. Although the sensitivity of streamflow to historic glacier area changes is small and within parameter uncertainties, our results suggest that glacier area changes have to be accounted for in future projections of late summer streamflow. Our approach provides an effective and widely applicable method to calibrate hydrologic models in glacier fed catchments, as well as to quantify the magnitude and timing of glacier melt contributions to streamflow.

Most models of glacier melt are forced by observed meteorological data in the vicinity of the glacier in question. In the absence of these observations, the forcing is commonly derived from statistical or dynamical downscaling of low resolution reanalysis datasets. Here we perform dynamical downscaling via the Weather Research and Forecasting (WRF) model in order to evaluate its performance against the observations from automatic weather stations (AWSs ) for three mountain glaciers in the interior of British Columbia, Canada over several summer seasons. The WRF model, nested within the ERA-Interim global reanalysis, produced output fields at 7.5 km and 2.5 km spatial resolution for all glaciers, as well as 1 km resolution for one of the glaciers. We find that the surface energy balance (SEB) model, forced by WRF at 2.5 km, adequately simulates the AWS-derived seasonal melt rates despite large biases in the individual SEB components. Overestimation of the number of clear sky days in WRF at 2.5 km explains the positive bias in the seasonal incoming shortwave radiation. This positive bias, however, is compensated by a negative bias in the seasonal incoming longwave radiation, and by an underestimation of sensible and latent heat fluxes. The underestimation of sensible heat fluxes down to −80 % of AWS-derived fluxes, as calculated by the bulk aerodynamic method, is due to the underestimated near-surface wind speeds. An increase of WRF spatial resolution from 7.5 to 1 km does not improve the simulation of downscaled variables other than near-surface air temperature. For relatively small glaciers (

Lillooet River in southwest British Columbia has produced damaging floods many times during the past century. The floods are recorded in Lillooet Lake, into which the river flows, as anomalously thick clastic varves. In order to determine whether an 825-year long varve record obtained from 12 percussion and vibracores can be used as flood proxy, we compare river discharge records dating back to 1914 to the thickness of the varves deposited during the same time period. Correlations between varve thickness and a variety of historical discharge measures are low to moderate for the periods 1914–2004 (r 2 = 0.37) and 1914–1945 (r 2 = 0.40), but higher for the period 1946–2004 (r 2 = 0.55). The best correlation (r 2 = 0.55) is between maximum fall discharge and varve thickness during the most recent period (1946–2004). Varve thickness for the earlier period of hydrometric data (1914–1946), which is a time of rapid glacier retreat and warmer temperatures in British Columbia, is best explained with a discharge proxy combining nival runoff, glacier runoff, and maximum fall discharge. Landslides, glacier fluctuations, river dyking, artificial lowering of Lillooet Lake, as well as lag effects of storms are responsible for the considerable unexplained variance in the relation between discharge measures and varve thickness over the historic period. The cores contain many anomalously thick varves, some of which we attribute to previously dated prehistoric landslides in the watershed or to local landslides into the lake. We conclude that many historic and prehistoric floods are faithfully recorded as anomalously thick clastic varves, but that other processes operating in the watershed preclude using this record as a reliable paleo-flood proxy.

Glaciers are major agents of erosion that increase sediment load to the downstream fluvial system. The Castle Creek Glacier, British Columbia, Canada, has retreated ~1.0 km in the past 70 years. Suspended sediment concentration (SSC) and streamflow (Q) were monitored independently at five sites within its pro-glacial zone over a 60 day period from July to September 2011, representing part of the ablation season. Meteorological data were collected from two automatic weather stations proximal to the glacier. The time-series were divided into hydrologic days and the shape and magnitude of the SSC response to hydro-meteorological conditions (‘cold and wet’, ‘hot and dry’, ‘warm and damp’, and ‘storm’) were categorized using principal component analysis (PCA) and cluster analysis (CA). Suspended sediment load (SSL) was computed and summarized for the categories. The distribution of monitoring sites and results of the multivariate statistical analyses describe the temporal and spatial variability of suspended sediment flux and the relative importance of glacial and para-glacial sediment sources in the pro-glacial zone. During the 2011 study period, ~ 60% of the total SSL was derived from the glacial stream and sediment deposits proximal to the terminus of the glacier; during ‘storm’ events, that contribution dropped to ~40% as the contribution from diffuse and point sources of sediment throughout the pro-glacial zone and within the meltwater channels increased. While ‘storm’ events accounted for just 3% of the study period, SSL was ~600% higher than the average over the monitoring period, and ~20% of the total SSL was generated in that time. Determining how hydro-meteorological conditions and sediment sources control sediment fluxes will assist attempts to predict how pro-glacial zones respond to future climate changes. Copyright © 2015 John Wiley & Sons, Ltd.

The glaciers in western Canada are experiencing rapid mass loss. Projections of their fate with a model that couples physics-based ice dynamics with a surface mass balance model suggest that glacier volume will shrink by 70% by 2100.

We present pollen and spore analysis from a peat-bog section from southwestern Isla de los Estados (Staaten Island), Tierra del Fuego to better understand late-Holocene environmental change in southernmost Patagonia. The island's position as the easternmost landmass at the southernmost of South America, places it squarely within the influence of the southern westerly winds (SWW) and the Antarctic Circumpolar Current, and so makes the island a unique location to document the palaeoecological response to climate change during the late Quaternary. We compare our Nothofagus pollen record from Bahía Franklin with other Nothofagus pollen records from Isla de los Estados and Tierra del Fuego. Nothofagus has been shown to respond to changes in mean annual temperature in the southern latitudes under cooler climate with little seasonality in precipitation. We also evaluate the correspondence between these changes in vegetation to other late-Holocene environmental proxies that include variations in sea surface temperature (SST) from the Beagle Channel and records of past glacier fluctuations. Our results show that low concentrations of Nothofagus pollen existed at Bahía Franklin at about 4000 cal yr BP. Climate ameliorated between ca. 3500 and 500 cal yr BP based on an increase in frequency and concentration values of Nothofagus and Drimys winteri. Between about 500 and 50 cal yr BP, the pollen record revealed a noticeable decline in the forest density. We interpret the decline in forest density in Isla de los Estados between 500 and ∼50 yr cal BP as a response to cold and windy conditions during the Little Ice Age (LIA). Our interpretation from our Nothofagus pollen record broadly agrees with regional cooling recorded in reconstructed SST of the Beagle Channel and with glacier expansion in the Fuegian Cordillera. Fil: Ponce, Juan Federico. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Austral de Investigaciones Científicas; Argentina. Universidad Nacional de Tierra del Fuego, Antártida e Islas del Atlántico Sur. Instituto de Ciencias Polares, Recursos Naturales y Ambiente. Departamento de Biología; Argentina Fil: Borromei, Ana Maria. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Bahía Blanca. Instituto Geológico del Sur. Universidad Nacional del Sur. Departamento de Geología. Instituto Geológico del Sur; Argentina Fil: Menounos, Brian. University of Northern British Columbia; Estados Unidos Fil: Rabassa, Jorge Oscar. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Austral de Investigaciones Científicas; Argentina. Universidad Nacional de Tierra del Fuego, Antártida e Islas del Atlántico Sur. Instituto de Ciencias Polares, Recursos Naturales y Ambiente. Departamento de Biología; Argentina

At least seven late Pliocene tills cap plateaus (mesetas) south of Lago Viedma, just east of the Andes in Argentine Patagonia. Chronologic constraints on the tills are provided by 40Ar/39Ar ages and magnetic polarities on associated basalt flows and sediments. The tills were deposited by piedmont glaciers that reached at least 80 km east of the crest of the Andes and flowed on a low-relief surface sloping gently downward in that direction. The oldest of the tills is about 3.6 Ma old. Glacial deposits dating to the Pliocene-Pleistocene transition are present at least 40 km beyond the east limit of the Pliocene tills at Lago Viedma, and tills of similar age occur at Condor Cliff in the Rio Santa Cruz valley to the southeast. A sequence of at least seven Early Pleistocene (2.1–1.1 Ma) tills is present between basalt flows in the Cerro del Fraile meseta south of Lago Argentino. The glaciers that deposited these Early Pleistocene tills reached far beyond the Marine Isotope Stage 2 limit in the Rio Santa Cruz valley. Based on positions, extents, and ages of the un-deformed, basalt-capped mesetas flanking Lago Viedma, we conclude that the topography in this area was profoundly changed during the Pleistocene – the low to moderate relief Pliocene surface was deeply incised by glaciers that became increasingly confined to, and flowed within, troughs. The valley floors today are up to 1350 m below the late Pliocene surface.

We estimate the glacier mass balance of a 9.5 km2 mountain glacier using three approaches for balance years 2009, 2010 and 2011. The photogrammetric, GPS and glaciological methods yielded sampling densities of 100, 5 and 2 points km-2, with measurement precisions of ± 0.40, ± 0.10 and ± 0.10 m w.e. respectively. Our glaciological measurements likely include a positive bias, due to omission of internal and basal mass balance, and uncertainty in determining the interface between snow and firn with a probe (±0.10 m w.e.). Measurements from our photogrammetric method include a negative bias introduced by the manual operator and our temperature index model used to correct for different dates of imaging (0.15 m w.e.), whereas GPS measurements avoid these biases. The photogrammetric and GPS methods are suitable for estimating glacier-wide annual mass balance, and thus provide a valuable measure that complements the glaciological method. These approaches, however, cannot be used to estimate mass balance at a point or mass-balance profiles without a detailed understanding of the vertical component of ice velocity.

We determined length, area, elevation and volume change of the Columbia Icefield using Interprovincial Boundary Commission Survey maps from 1919, eight sets of aerial photographs from 1948 to 1993, and satellite data from 1999 to 2009. Over the period 1919–2009, glaciers on average retreated 1150 ± 34 m and shrank by 2.4 ± 0.2 km2. Total area loss was 59.6 ± 1.2 km2 (23 ± 5%), and mean elevation change was −49 ± 25 m w.e., resulting in a total volume loss of 14.3 ± 2.0 km3 w.e. Large outlet glaciers experienced the greatest absolute ice loss, while small, detached glaciers lost the most relative length and area. Thinning rates of debris-covered ice were 30–60% lower than those for clean ice. All glacier changes were significantly correlated with each other (p < 0.01), with r values ranging from 0.54 to 0.82. Temperature is correlated with length and area change over periods lagged 1–5 years (p < 0.05), and with elevation and volume change over periods lagged 9–18 years (p < 0.05). Precipitation is correlated with glacier change over periods lagged 1–10 years (p < 0.05).

The Lillooet-Harrison watershed (7870 km2 in area, of which 459 km2 is glacierized) contains recurring morphological landscape elements that differ in terms of their level of connectivity with Lillooet River. First, ridge tops with low local relief are connected to Lillooet River primarily through dissolved solids moving along subsurface pathways. Second, small glaciers have little direct connection with Lillooet River because of their numerous proximal lakes, whereas larger valley glaciers are directly linked to Lillooet River. Third, hillslopes subject to infrequent, large mass movements are episodically connected, and fourth, debris flow channels are directly linked with valley-bottom rivers. Fifth, mid-valley and lower valley-bottom sites are sediment storage areas that are episodically coupled to the adjacent river channel. The Mount Meager Volcanic Complex dominates the supply of coarse clastic sediment to Lillooet River and is primarily responsible for the rapid rates of Lillooet delta progradation and aggradation. The 6 % of the watershed that is glacierized provides a high proportion of the fine clastic sediments to Lillooet and Harrison lakes, whereas solute production that has little influence on morphological elements is important in assessing the overall rate of landscape denudation in the Lillooet-Harrison watershed. More frequent extreme rainfall events and rapid glacier retreat together with intensified human occupancy in the valley can be expected to increase the coupling of each of the morphological elements to Lillooet River and expedite the movement of water, sediments, and solutes through the watershed.