Benjamin Brock

An energy-balance model for debris-covered glaciers including heat conduction through the debris layer

Extensive covers of supraglacial debris are often present in glacier ablation areas, and it is essential to assess exactly how the debris affects glacier melt rates. This paper presents a physically based energy-balance model for the surface of a debris-covered glacier. The model is driven by meteorological variables, and was developed using data collected at Miage glacier, Italy, during the ablation seasons of 2005, 2006 and 2007. The debris surface temperature is numerically estimated by considering the balance of heat fluxes at the air/debris interface, and heat conduction through the debris is calculated in order to estimate melt rates at the debris/ice interface. The predicted hourly debris surface temperatures and debris internal temperatures provide a good fit to temperatures measured on rock-covered Miage glacier (r 2 > 0.94) and the tephra-covered glacier on Villarrica volcano, Chile (r 2 > 0.82). The model can also be used to reproduce observed changes in melt rates below debris layers of varying types and thicknesses, an important consideration for the overall mass balance of debris-covered glaciers.

The surface energy balance of an active ice-covered volcano: Villarrica Volcano, Southern Chile

Abstract The energy balance of bare snow and tephra-covered ice near the glacier equilibrium line elevation on Villarrica Volcano, southern Chile, was investigated during 2004 and 2005, combining meteorological, surface temperature and ablation measurements with energy balance modelling. A tephra thermal conductivity of 0.35 Wm–1 K–1, and a critical tephra thickness of <5mm at which ablation is reduced compared to bare snow, were obtained from field data. These low values are attributable to the highly porous lapilli particles which make up most of the surface material. Modelled melt totals in the January to March period were 4.95 m and 3.96 m water equivalent (w.e.) in 2004 and 2005, respectively, compared with ∽0.5mw.e. melt for ice buried by >0.1m tephra. Windblown tephra impurities lowered snow albedo, but increased snowmelt by only an estimated 0.28mw.e. over the same period. The net mass balance impact of supraglacial tephra at Villarrica Volcano is therefore positive, as thick ash and lapilli mantle most of the glacier ablation zones, probably reducing annual ablation by several metres w.e. In the accumulation seasons, frequent melting events were recorded with modelled daily snowmelt rates of up to 50 mmw.e.

Ice volumetric changes on active volcanoes in southern Chile

Abstract Most of the glaciers in southern Chile have been retreating and shrinking during recent decades in response to atmospheric warming and decrease in precipitation. However, some glacier fluctuations are directly associated with the effusive and geothermal activity of ice-covered active volcanoes widely distributed in the region. The aim of this paper is to study the ice volumetric changes by comparing several topographic datasets. A maximum mean ice thinning rate of 0.81 ± 0.45 m a−1 was observed on the ash/debris-covered ablation area of Volcan Villarrica between 1961 and 2004, whilst on Volcan Mocho the signal-to-noise ratio was too small to yield any conclusion. An area reduction of 0.036 ±0.019 km2 a−1 since 1976 was obtained on Glaciar Mocho, while on Volcan Villarrica the area change was −0.090 ± 0.034 km2 a−1 between 1976 and 2005. Glaciers on active volcanoes are therefore shrinking, mainly in response to climatic driving factors. However, volcanic activity is affecting glaciers in two opposite ways: ash/debris advection is helping to reduce surface ablation at lower reaches by insulating the ice from solar radiation, while geothermal activity is probably enhancing melting and water production at the bedrock, resulting in negative ice-elevation changes.

Suitability of a constant air temperature lapse rate over an Alpine glacier: testing the Greuell and Böhm model as an alternative

Abstract Near-surface air temperature, typically measured at a height of 2 m, is the most important control on the energy exchange and the melt rate at a snow or ice surface. It is distributed in a simplistic manner in most glacier melt models by using constant linear lapse rates, which poorly represent the actual spatial and temporal variability of air temperature. In this paper, we test a simple thermodynamic model proposed by Greuell and Böhm in 1998 as an alternative, using a new dataset of air temperature measurements from along the flowline of Haut Glacier d’Arolla, Switzerland. The unmodified model performs little better than assuming a constant linear lapse rate. When modified to allow the ratio of the boundary layer height to the bulk heat transfer coefficient to vary along the flowline, the model matches measured air temperatures better, and a further reduction of the root-mean-square error is obtained, although there is still considerable scope for improvement. The modified model is shown to perform best under conditions favourable to the development of katabatic winds – few clouds, positive ambient air temperature, limited influence of synoptic or valley winds and a long fetch – but its performance is poor under cloudy conditions.

An enhanced temperature index model for debris-covered glaciers accounting for thickness effect

Highlights • We develop a melt model for debris-covered ice accounting for debris thickness.• The model can reproduce the melt reduction of the so-called Oestrem curve.• The model is transferable in time over seasons and in space to a second glacier.

A fifty year record of winter glacier melt events in southern Chile, 38??42?S

Little is known about the frequency and potential mass balance impact of winter glacier melt events. In this study, daily atmospheric temperature soundings from the Puerto Montt radiosonde (41.43°S) are used to reconstruct winter melting events at the glacier equilibrium line altitude in the 38°–42°S region of southern Chile, between 1960 and 2010. The representativeness of the radiosonde temperatures to near-surface glacier temperatures is demonstrated using meteorological records from close to the equilibrium line on two glaciers in the region over five winters. Using a degree-day model we estimate an average of 0.28 m of melt and 21 melt days in the 15 June–15 September period each year, with high inter-annual variability. The majority of melt events are associated with midlatitude migratory high pressure systems crossing Chile and northwesterly flows, that force adiabatic compression and warm advection, respectively. There are no trends in the frequency or magnitude of melt events over the period of record, but the annual frequency of winter melt days shows a significant, although rather weak and probably non-linear, relationship to late winter and early spring values of a multivariate El Nino Southern Oscillation Index (MEI).

Cryospheric ecosystems: a synthesis of snowpack and glacial research

The fourteen letters that contributed to this focus issue on cryospheric ecosytems provide an excellent basis for considering the state of the science following a marked increase in research attention since the new millennium. Research letters from the focus issue provide significant insights into the biogeochemical and biological processes associated with snow, glacier ice and glacial sediments. This has been achieved via a significant, empirical effort that has given particular emphasis to glacier surface habitats. However, far less is known about aerobiology, glacial snow covers, supraglacial lakes and sub-ice sedimentary habitats, whose access for sampling and in-situ monitoring remains a great challenge to scientists. Furthermore, the use of models to explore key fluxes, processes and impacts of a changing glacial cryosphere are conspicuous by their absence. As a result, a range of process investigations and modelling studies are required to address the increasing urgency and uncertainty that is associated with understanding the response of cryospheric ecosystems to global change.