Search

Spaceborne microwave remote sensing (300 MHz–100 GHz) provides a valuable method for characterizing environmental changes, especially in Arctic–boreal regions (ABRs) where ground observations are generally spatially and temporally scarce. Although direct measurements of carbon fluxes are not feasible, spaceborne microwave radiometers and radar can monitor various important surface and near-surface variables that affect terrestrial carbon cycle processes such as respiratory carbon dioxide (CO2) fluxes; photosynthetic CO2 uptake; and processes related to net methane (CH4) exchange including CH4 production, transport and consumption. Examples of such controls include soil moisture and temperature, surface freeze–thaw cycles, vegetation water storage, snowpack properties and land cover. Microwave remote sensing also provides a means for independent aboveground biomass estimates that can be used to estimate aboveground carbon stocks. The microwave data record spans multiple decades going back to the 1970s with frequent (daily to weekly) global coverage independent of atmospheric conditions and solar illumination. Collectively, these advantages hold substantial untapped potential to monitor and better understand carbon cycle processes across ABRs. Given rapid climate warming across ABRs and the associated carbon cycle feedbacks to the global climate system, this review argues for the importance of rapid integration of microwave information into ABR terrestrial carbon cycle science.

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 were used to force a physically based, distributed glacier mass balance model of Saskatchewan Glacier for the historical period 1979–2016 and assess its sensitivity to climate change. A 2-year record (2014–2016) from an on-glacier automatic weather station (AWS) and historical precipitation records from nearby permanent weather stations were used to downscale air temperature, relative humidity, wind speed, incoming solar radiation and precipitation from the NARR to the station sites. The model was run with fixed (1979, 2010) and time-varying (dynamic) geometry using a multitemporal digital elevation model 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 was not very sensitive to the NARR spatial interpolation method, as long as station data were used for bias correction. The simulated mass balance was however sensitive to the biases in NARR precipitation and air temperature, as well as to the prescribed precipitation lapse rate and ice aerodynamic roughness lengths, showing the importance of constraining these two parameters with ancillary data. The glacier-wide simulated energy balance regime showed a large contribution (57 %) of turbulent (sensible and latent) heat fluxes to melting in summer, higher than typical mid-latitude glaciers in continental climates, which reflects the local humid “icefield weather” of the Columbia Icefield. The static mass balance sensitivity to climate was assessed for prescribed changes in regional mean air temperature between 0 and 7 ∘C and precipitation between −20 % and +20 %, which comprise the spread of ensemble Representative Concentration Pathway (RCP) 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. a−1 ∘C−1. The mass balance response to warming was driven by a positive albedo feedback (44 %), followed by direct atmospheric warming impacts (24 %), a positive air humidity feedback (22 %) and a positive precipitation phase feedback (10 %). 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.

Accurate knowledge of snow depth distributions in forested regions is crucial for applications in hydrology and ecology. Understanding and assessing the effect of vegetation and topographic conditions on the snow depth variability is useful for the accurate prediction of snow depths. In this study, the spatial distribution of snow depth in two agro-forested sites and one coniferous site in eastern Canada was analyzed for topographic and vegetation effects on snow accumulation. Spatially distributed snow depths were derived by Unmanned Aerial Vehicle Light Detection and Ranging (UAV-lidar) surveys conducted in 2019 and 2020. Distinct patterns of snow accumulation and erosion in open areas (fields) versus adjacent forested areas were observed in lidar-derived snow depth maps at all sites. Omnidirectional semi-variogram analysis of snow depths showed the existence of a scale break distance less than 10 m in the forested area at all three sites, whereas open areas showed scale invariance or comparatively large scale break distances (i.e., 18 m). The effect of vegetation and topographic variables on the spatial variability of snow depths at each site was investigated with random forest models. Results show that including wind-related forest edge proximity effects improved the model accuracy by more than 50 % in agro-forested sites, whereas incorporating canopy characteristics improved the model accuracy by more than 60 % in the coniferous site. Hence the underlying topography and the wind-redistribution of snow along forest edges govern the snow depth variability at agro-forested sites, while forest structure variability dominates snow depth variability in the coniferous environment. These results highlight the importance of including and better representing these processes in process-based models for accurate estimates of snowpack dynamics. This study also demonstrates the usefulness of UAV-lidar to resolve and understand high-resolution snow depth heterogeneity in agro-forested environments and boreal forests.

In the Atlas Mountains range, streamflow is largely generated from meltwater supplied by the snowpack during spring and early summer. In this manner, snow is considered an important factor which determining water availability in semi-arid and arid mountains. This substantial part of freshwater stored in the form of snow contributes significantly to mountainous runoff. However, the contribution of snow and rain to the annual and multi-annual water balance remain largely unknown. Hydrological modeling is needed to support water resource assessment and management in the Atlas range. As meteorological data is often scarce, the models must be able to simulate the spatiotemporal heterogeneity of forcing variables while maintaining a low data input requirement.In this study, the performance of the snowmelt runoff model (SRM) is assessed to simulate and forecast daily runoff essentially from snowmelt and rainfall at the Rheraya watershed in the Moroccan High Atlas range over the 2010 - 2016 period. The SRM runoff simulation is tested under two forcing inputs: (i) four snowmelt rates previously estimated by a classical temperature-index model (TI) and three enhanced temperature index models that respectively include the potential clear-sky direct radiation (HTI), the incoming solar radiation (ETI-A), and net solar radiation (ETI-B); (ii) calculated snowmelt from the snow cover area (SCA) products of Moderate-Resolution Imaging Spectroradiometer (MODIS).All SRM simulated runoff were subjected to calibration and validation through the measured runoff in the Tahanaout weather station. The sensibility of recession coefficients was also evaluated. The SRM simulations results over the validation period show an acceptable performance.Keywords: Runoff, SRM, snowmelt, SCA, temperature index model; enhanced degree-day models, MODIS, semi-arid climate, Rheraya, High Atlas, Morocco.

Oceans and cryosphere are immediately affected by human-induced climate change. Our Editorial Board members present viewpoints on the most pressing and fruitful avenues of research on frozen and liquid water, and the transition from one to the other.

Cold region hydrology is conditioned by distinct cryospheric and hydrological processes. While snowmelt is the main contributor to both surface and subsurface flows, seasonally frozen soil also influences the partition of meltwater and rain between these flows. Cold regions of the Northern Hemisphere midlatitudes have been shown to be sensitive to climate change. Assessing the impacts of climate change on the hydrology of this region is therefore crucial, as it supports a significant amount of population relying on hydrological services and subjected to changing hydrological risks. We present an exhaustive review of the literature on historical and projected future changes on cold region hydrology in response to climate change. Changes in snow, soil, and streamflow key metrics were investigated and summarized at the hemispheric scale, down to the basin scale. We found substantial evidence of both historical and projected changes in the reviewed hydrological metrics. These metrics were shown to display different sensitivities to climate change, depending on the cold season temperature regime of a given region. Given the historical and projected future warming during the 21st century, the most drastic changes were found to be occurring over regions with near-freezing air temperatures. Colder regions, on the other hand, were found to be comparatively less sensitive to climate change. The complex interactions between the snow and soil metrics resulted in either colder or warmer soils, which led to increasing or decreasing frost depths, influencing the partitioning rates between the surface and subsurface flows. The most consistent and salient hydrological responses to both historical and projected climate change were an earlier occurrence of snowmelt floods, an overall increase in water availability and streamflow during winter, and a decrease in water availability and streamflow during the warm season, which calls for renewed assessments of existing water supply and flood risk management strategies.

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.

Excluding Antarctica and Greenland, 3.8% of the world&rsquo, s glacier area is concentrated in Chile. The country has been strongly affected by the mega drought, which affects the south-central area and has produced an increase in dependence on water resources from snow and glacier melting in dry periods. Recent climate change has led to an elevation of the zero-degree isotherm, a decrease in solid-state precipitation amounts and an accelerated loss of glacier and snow storage in the Chilean Andes. This situation calls for a better understanding of future water discharge in Andean headwater catchments in order to improve water resources management in glacier-fed populated areas. The present study uses hydrological modeling to characterize the hydrological processes occurring in a glacio-nival watershed of the central Andes and to examine the impact of different climate change scenarios on discharge. The study site is the upper sub-watershed of the Tinguiririca River (area: 141 km2), of which nearly 20% is covered by Universidad Glacier. The semi-distributed Snowmelt Runoff Model + Glacier (SRM+G) was forced with local meteorological data to simulate catchment runoff. The model was calibrated on even years and validated on odd years during the 2008&ndash, 2014 period and found to correctly reproduce daily runoff. The model was then forced with downscaled ensemble projected precipitation and temperature series under the RCP 4.5 and RCP 8.5 scenarios, and the glacier adjusted using a volume-area scaling relationship. The results obtained for 2050 indicate a decrease in mean annual discharge (MAD) of 18.1% for the lowest emission scenario and 43.3% for the most pessimistic emission scenario, while for 2100 the MAD decreases by 31.4 and 54.2%, respectively, for each emission scenario. Results show that decreasing precipitation lead to reduced rainfall and snowmelt contributions to discharge. Glacier melt thus partly buffers the drying climate trend, but our results show that the peak water occurs near 2040, after which glacier depletion leads to reducing discharge, threatening the long-term water resource availability in this region.

Four models of increasing complexity were tested and compared to simulate snow water equivalent at the local scale in the Moroccan High Atlas range. A classical temperature index model (TI) and thr...

Black carbon aerosol (BC), which is emitted from natural and anthropogenic sources (e.g., wildfires, coal burning), can contribute to magnify climate warming at high latitudes by darkening snow- and ice-covered surfaces, and subsequently lowering their albedo. Therefore, modeling the atmospheric transport and deposition of BC to the Arctic is important, and historical archives of BC accumulation in polar ice can help to validate such modeling efforts. Here we present a > 250-year ice-core record of refractory BC (rBC) deposition on Devon ice cap, Canada, spanning the years from 1735 to 1992. This is the first such record ever developed from the Canadian Arctic. The estimated mean deposition flux of rBC on Devon ice cap for 1963–1990 is 0.2 mg m−2 a−1, which is at the low end of estimates from Greenland ice cores obtained using the same analytical method ( ∼ 0.1–4 mg m−2 a−1). The Devon ice cap rBC record also differs from the Greenland records in that it shows only a modest increase in rBC deposition during the 20th century. In the Greenland records a pronounced rise in rBC is observed from the 1880s to the 1910s, which is largely attributed to midlatitude coal burning emissions. The deposition of contaminants such as sulfate and lead increased on Devon ice cap in the 20th century but no concomitant rise in rBC is recorded in the ice. Part of the difference with Greenland could be due to local factors such as melt–freeze cycles on Devon ice cap that may limit the detection sensitivity of rBC analyses in melt-impacted core samples, and wind scouring of winter snow at the coring site. Air back-trajectory analyses also suggest that Devon ice cap receives BC from more distant North American and Eurasian sources than Greenland, and aerosol mixing and removal during long-range transport over the Arctic Ocean likely masks some of the specific BC source–receptor relationships. Findings from this study suggest that there could be a large variability in BC aerosol deposition across the Arctic region arising from different transport patterns. This variability needs to be accounted for when estimating the large-scale albedo lowering effect of BC deposition on Arctic snow/ice.

In this work we propose a method of estimating the age of calving events of a tidewater glacier based on 3D point clouds. Terrestrial LiDAR scans of the calving front of Fuerza Aerea Glacier (Greenwich Island, South Shetland Islands, Antarctica) were analyzed to determine the timing of the ice face exposure by calving events. Median reflection intensity of laser beam within a calving event footprint (I) was shown to be the best proxy measure of the ice exposure time (tLC). Point clouds were used in concert with information about calving location to create a mathematical model of the intensity-time relationship. An exponential function in the form of I = a ⋅ exp(b/tLC) + c has been shown to describe the relationship between median reflection intensity and time elapsed since the last calving event in a given part of the glacier. The model explains over 60% of the variability in reflection intensity between calving events. The shape of the model curve limits the prediction capability to events occurring up to 3 days before a scan. Development of weathered or superimposed ice on the calving scars is proposed as the explanation for increase of ice reflectivity with time.

This study explores the potential impacts of climate change on soil erosion in an agricultural catchment in eastern Canada. The Modified Universal Soil Loss Equation (MUSLE) was used to calculate the sediment yields from the Acadie River Catchment for the historical 1996–2019 period. The runoff variables of the MUSLE were obtained from a physically based hydrological model previously built and validated for the catchment. Then, the hydrological model was perturbed using climate change projections and used to assess the climate sensitivity of the sediment yield. Two runoff types representing possible modes of soil erosion were considered. While type A represents a baseline case in which soil erosion occurs due to surface runoff only, type B is more realistic since it assumed that tile drains also contribute to sediment export, but with a varying efficiency throughout the year. The calibration and validation of the tile efficiency factors against measurements in 2009–2015 for type B suggest that tile drains export the sediments with an efficiency of 20% and 50% in freezing and non-freezing conditions, respectively. Results indicate that tile drains account for 39% of the total annual sediment yield in the present climate. The timing of highest soil erosion shifts from spring to winter in response to warming and wetting, which can be explained by increasing winter runoff caused by shifting snowmelt timing towards winter, a greater number of mid-winter melt events as well as increasing rainfall fractions. The large uncertainties in precipitation projections cascade down to the erosion uncertainties in the more realistic type B, with annual sediment yield increasing or decreasing according to the precipitation uncertainty in a given climate change scenario. This study demonstrates the benefit of conservation and no-till pratices, which could reduce the annual sediment yields by 20% and 60%, respectively, under any given climate change scenario.

Three glacier–rock glacier transitional landforms in the central Andes of Chile are investigated over the last decades in order to highlight and question the significance of their landscape and flow dynamics. Historical (1955–2000) aerial photos and contemporary (> 2000) Geoeye satellite images were used together with common processing operations, including imagery orthorectification, digital elevation model generation, and image feature tracking. At each site, the rock glacier morphology area, thermokarst area, elevation changes, and horizontal surface displacements were mapped. The evolution of the landforms over the study period is remarkable, with rapid landscape changes, particularly an expansion of rock glacier morphology areas. Elevation changes were heterogeneous, especially in debris-covered glacier areas with large heaving or lowering up to more than ±1 m yr−1. The use of image feature tracking highlighted spatially coherent flow vector patterns over rock glacier areas and, at two of the three sites, their expansion over the studied period; debris-covered glacier areas are characterized by a lack of movement detection and/or chaotic displacement patterns reflecting thermokarst degradation; mean landform displacement speeds ranged between 0.50 and 1.10 m yr−1 and exhibited a decreasing trend over the studied period. One important highlight of this study is that, especially in persisting cold conditions, rock glaciers can develop upward at the expense of debris-covered glaciers. Two of the studied landforms initially (prior to the study period) developed from an alternation between glacial advances and rock glacier development phases. The other landform is a small debris-covered glacier having evolved into a rock glacier over the last half-century. Based on these results it is proposed that morphological and dynamical interactions between glaciers and permafrost and their resulting hybrid landscapes may enhance the resilience of the mountain cryosphere against climate change.

The Las Liebres rock glacier is a large (~2.2 km long) Andean rock glacier whose internal composition and kinematics are known from previous studies. We investigate its development by posing and testing the following null hypothesis: the rock glacier has developed from a constant supply of debris and ground ice in periglacial conditions and resulting creep of the ice-rock mixture. A rheological model was formulated based on recent advances in the study of ice-rock mixture rheology, and calibrated on the known surface velocities and internal composition of the rock glacier. We show that the rock glacier viscosity is inversely related to both water and debris fractions, in agreement with recent field and theoretical studies of ice-rock mixture rheology. Taking into account the possible variations in water fraction, the model was used to estimate the time spans of development (0.91–7.11 ka), rates of rock wall retreat (0.44–4.18 mm/a), and rates of ground ice formation (0.004–0.026 m/a) for the rock glacier. These results support the null hypothesis of a periglacial origin of the Las Liebres rock glacier. Copyright © 2016 John Wiley & Sons, Ltd.

Using a coastal ice core collected from Prince of Wales (POW) Icefield on Ellesmere Island, we investigate source regions of sea ice-modulated chemical species (methanesulfonic acid (MSA) and chloride (Cl−)) to POW Icefield and the influence of large-scale atmospheric variability on the transport of these marine aerosols (1979–2001). Our key findings are (1) MSA in the POW Icefield core is derived primarily from productivity in the sea ice zone of Baffin Bay and the Labrador Sea, with influence from waters within the North Water (NOW) polynya, (2) sea ice formation processes within the NOW polynya may be a significant source of sea-salt aerosols to the POW core site, in addition to offshore open water source regions primarily in Hudson Bay, and (3) the tropical Pacific influences the source and transport of marine aerosols to POW Icefield through its remote control on regional winds and sea ice variability. Regression analyses during times of MSA deposition reveal sea level pressure (SLP) anomalies favorable for opening of the NOW polynya and subsequent oceanic dimethyl sulfide production. Regression analyses during times of Cl− deposition reveal SLP anomalies that indicate a broader oceanic region of sea-salt sources to the core site. These results are supported by Scanning Multichannel Microwave Radiometer- and Special Sensor Microwave/Imager-based sea ice reconstructions and air mass transport density analyses and suggest that the marine biogenic record may capture local polynya variability, while sea-salt transport to the site from larger offshore source regions in Baffin Bay is likely. Regression analyses show a link to tropical dynamics via an atmospheric Rossby wave.

A short-term series of quantitative observations of calving activity of Fuerza Aérea Glacier (Greenwich Island, the South Shetland Islands, Antarctica) was conducted in order to test new methods of monitoring calving. The volume of single calving events was quantified by combining terrestrial laser scanning (TLS) surveys with continuous video recording of the ice front. An empirical formula for area/volume scaling of the calved ice block was proposed based on the TLS measured calved ice volume and the calved ice front area obtained by manual delineation on the images acquired with the video camera. This combination of methods proves to be a valuable tool for glacier monitoring, providing both high-temporal resolution and precise quantitative measurements of the calving volume. The size distribution of calving events is best approximated by a power law and within the short period of observations (14 d) calving was found to be an intrinsic process not dependent on environmental forcings. Over the period of 21 January–04 February 2013 the ice flow velocity at the terminus of Fuerza Aérea Glacier was 0.26 ± 0.07 m d−1and the calving rate was 0.41 ± 0.07 m d−1.

Having a realistic estimation of snow cover by conceptual hydrological models continues to challenge hydrologists. The calibration of the free model parameters is an unavoidable step and the uncertainties resulting from the use of this optimal set remains a source of concern, especially in forecasting applications and climate changes impact assessments. This study seeks to improve the calibration of the conceptual hydrological model GR4J coupled with the Cemaneige snow model, in order to obtain a more realistic simulation of the snow water equivalent (SWE) and to reduce the uncertainty of the free parameters. The performance of the two models was tested over twelve snow-dominated basins in southern Quebec, Canada. Four calibration strategies were adopted and compared. In the first two strategies, the parameters were calibrated against observed streamflow alone using a local and a global algorithm. In the third and fourth strategies the calibration of snow and hydrological parameters was performed against observed streamflow and snow water equivalent (SWE) measured at snow course transects, first separately, and then with a multiobjective approach. An ensemble of equifinal parameters was used to compare the capacity of the global and multiobjective algorithms to improve the parameters identifiability and to assess the impact of parameter equifinality on the temperature sensitivity of spring peak streamflow. The large number of equifinal parameters found during calibration underscores the importance of structural non-identifiability of the coupled GR4J-Cemaneige model. The inclusion of snow observations within a multiobjective calibration improved the simulation of SWE, the identifiability of the parameters and their correlation with basins characteristics. Parameter equifinality caused a small but non negligible uncertainty in the simulated response of spring peak flow to warming temperatures. Parameter equifinality should be considered in climate impact studies in snow-dominated basins where poorly constrained snow parameters can affect the temperature sensitivity of streamflow.

Meteorological and surface change measurements collected during a 2.5 yr period are used to calculate surface mass and energy balances at 5324 m a.s.l. on Guanaco Glacier, a cold-based glacier in the semi-arid Andes of Chile. Meteorological conditions are marked by extremely low vapour pressures (annual mean of 1.1 hPa), strong winds (annual mean of 10 m s−1), shortwave radiation receipt persistently close to the theoretical site maximum during cloud-free days (mean annual 295 W m−2; summer hourly maximum 1354 W m−2) and low precipitation rates (mean annual 45 mm w.e.). Snowfall occurs sporadically throughout the year and is related to frontal events in the winter and convective storms during the summer months. Net shortwave radiation provides the greatest source of energy to the glacier surface, and net longwave radiation dominates energy losses. The turbulent latent heat flux is always negative, which means that the surface is always losing mass via sublimation, which is the main form of ablation at the site. Sublimation rates are most strongly correlated with net shortwave radiation, incoming shortwave radiation, albedo and vapour pressure. Low glacier surface temperatures restrict melting for much of the period, however episodic melting occurs during the austral summer, when warm, humid, calm and high pressure conditions restrict sublimation and make more energy available for melting. Low accumulation (131 mm w.e. over the period) and relatively high ablation (1435 mm w.e.) means that mass change over the period was negative (−1304 mm w.e.), which continued the negative trend recorded in the region over the last few decades.

Both satellite and ground-based broadband albedo measurements over rough and complex terrain show several limitations concerning feasibility and representativeness. To assess these limitations and understand the effect of surface roughness on albedo, firstly, an intrasurface radiative transfer (ISRT) model is combined with albedo measurements over different penitente surfaces on Glaciar Tapado in the semi-arid Andes of northern Chile. Results of the ISRT model show effective albedo reductions over the penitentes up to 0.4 when comparing the rough surface albedo relative to the albedo of the flat surface. The magnitude of these reductions primarily depends on the opening angles of the penitentes, but the shape of the penitentes and spatial variability of the material albedo also play a major role. Secondly, the ISRT model is used to reveal the effect of using albedo measurements at a specific location (i.e., apparent albedo) to infer the true albedo of a penitente field (i.e., effective albedo). This effect is especially strong for narrow penitentes, resulting in sampling biases of up to ±0.05. The sampling biases are more pronounced when the sensor is low above the surface, but remain relatively constant throughout the day. Consequently, it is important to use a large number of samples at various places and/or to locate the sensor sufficiently high in order to avoid this sampling bias of surface albedo over rough surfaces. Thirdly, the temporal evolution of broadband albedo over a penitente-covered surface is analyzed to place the experiments and their uncertainty into a longer temporal context. Time series of albedo measurements at an automated weather station over two ablation seasons reveal that albedo decreases early in the ablation season. These decreases stabilize from February onwards with variations being caused by fresh snowfall events. The 2009/2010 and 2011/2012 seasons differ notably, where the latter shows lower albedo values caused by larger penitentes. Finally, a comparison of the ground-based albedo observations with Landsat and MODIS (Moderate Resolution Imaging Spectroradiometer)-derived albedo showed that both satellite albedo products capture the albedo evolution with root mean square errors of 0.08 and 0.15, respectively, but also illustrate their shortcomings related to temporal resolution and spatial heterogeneity over small mountain glaciers. ispartof: Cryosphere vol:8 issue:3 pages:1069-1086 status: published

We investigated the Sachette rock glacier, Vanoise Massif, Northern French Alps, using former equilibrium line altitude reconstruction from glacial deposits, aerial photograph analysis, and ground-penetrating radar (GPR). The rock glacier is a young (probably −1. The constant-offset GPR data was processed and analysed in order to reconstruct the stratigraphy and model the radar wave velocity in two dimensions. The integration of the morphology, the velocity models, and the stratigraphy emphasized, in the upper half of the rock glacier, the good correspondence between high radar wave velocities (>0.15–0.16 m ns−1) and reflectors having a dipping-syncline structure, typical of true glaciers. Consequently, the rock glacier structure is described as being constituted of a glacial massive ice core embedded into diamictons. Our study of the Sachette rock glacier highlights possible significance of rock glaciers and interactions between glacier and permafrost in alpine environments.

The glacier to rock glacier transformation problem is revisited from a previously unseen angle. A striking case in the Juncal Massif (located in the upper Aconcagua Valley, Chilean central Andes) is documented. There, the Presenteseracae debris-covered glacier has advanced several tens of metres and has developed a rock glacier morphology in its lower part over the last 60 years. The conditions for a theoretically valid glacier to rock glacier transformation are discussed and tested. Permafrost probability in the area of the studied feature is highlighted by regional-scale spatial modelling together with on-site shallow ground temperature records. Two different methods are used to estimate the mean surface temperature during the summer of 2014, and the sub-debris ice ablation rates are calculated as ranging between 0.05 and 0.19 cm d− 1, i.e., 0.04 and 0.17 m over the summer. These low ablation rates are consistent with the development of a coherent surface morphology over the last 60 years. Furthermore, the rates of rock wall retreat required for covering the former glacier at Presenteseracae lie within the common 0.1–2 mm y− 1 range, assuming an average debris thickness and a range of debris-covering time intervals. The integration of the geomorphological observations with the numerical results confirms that the studied debris-covered glacier is evolving into a rock glacier.

Using ground-penetrating radar (GPR), we studied an entire 2.2km long rock glacier (3780–4350m asl) in the dry Andes of Chile with the aim of inferring its composition. In the high-quality, unmigrated data, we identified the active layer base and the rock glacier floor. In between, hyperbolae generated by diffracting boulders were inventoried; the ones along the rock glacier floor (n=51) allowed determination of the average electromagnetic (EM) velocity in the rock glacier, the latter being further used for migration. Within the rock glacier (16–39m thick), the EM velocity varies between 0.076 and 0.167m.ns -1 ; the main stratigraphic features observed are upward-dipping reflectors. The low EM velocities (