We review new advances in volcano–ice interactions on Mars and focus additional attention on (1) recent analyses of the mechanisms of penetration of the cryosphere by dikes and sills; (2) documentation of the glacial origin of huge fan-shaped deposits on the northwest margins of the Tharis Montes and evidence for abundant volcano–ice interactions during the later Amazonian period of volcanic edifice construction and (3) the circumpolar Hesperian-aged Dorsa Argentea Formation, interpreted as an ice sheet and displaying marginal features (channels, lakes and eskers) indicative of significant melting and interior features interpreted to be due to volcano–ice interactions (e.g. subglacial volcanic edifices, pits, basins, channels and eskers). In this context, we describe and analyse several stages and types of volcano–ice interactions: (1) magmatic interactions with ice-rich parts of the cryosphere; (2) subglacial volcanism represented by intrusion under and into the ice and formation of dikes and moberg-like ridges, intrusion of sills at the glacier–volcano substrate interface and their evolution into subglacial lava flows, formation of subglacial edifices, marginal melting and channels; (3) synglacial (ice contact) volcanism represented by flows banking up against glacier margins, chilling and forming remnant ridges and (4) post-glacial volcanism and interactions with ice deposits.

The process of lava-dome emplacement through a glacier was observed for the first time after Mount St Helens reawakened in September 2004. The glacier that had grown in the crater since the cataclysmic 1980 eruption was split in two by the new lava dome. The two parts of the glacier were successively squeezed against the crater wall. Photography, photogrammetry and geodetic measurements document glacier deformation of an extreme variety, with strain rates of extraordinary magnitude as compared to normal alpine glaciers. Unlike normal temperate glaciers, the crater glacier shows no evidence of either speed-up at the beginning of the ablation season or diurnal speed fluctuations during the ablation season. Thus there is evidently no slip of the glacier over its bed. The most reasonable explanation for this anomaly is that meltwater penetrating the glacier is captured by a thick layer of coarse rubble at the bed and then enters the volcano’s groundwater system rather than flowing through a drainage network along the bed.

The debris-rich glacier that grew in the crater of Mount St Helens after the volcano’s cataclysmic 1980 eruption was split in two by a new lava dome in 2004. For nearly six months, the eastern part of the glacier was squeezed against the crater wall as the lava dome expanded. Glacier thickness nearly doubled locally and surface speed increased substantially. As squeezing slowed and then stopped, surface speed fell and ice was redistributed downglacier. This sequence of events, which amounts to a field-scale experiment on the deformation of debris-rich ice at high strain rates, was interpreted using a two-dimensional flowband model. The best match between modeled and observed glacier surface motion, both vertical and horizontal, requires ice that is about 5 times stiffer and 1.2 times denser than normal, temperate ice. Results also indicate that lateral squeezing, and by inference lava-dome growth adjacent to the glacier, likely slowed over a period of about 30 days rather than stopping abruptly. This finding is supported by geodetic data documenting dome growth.

The Gorshkov crater glacier at Ushkovsky volcano, Kamchatka, is characterized by a large aspect ratio and special thermodynamic conditions at the bedrock caused by a locally enhanced and spatially varying geothermal heat flux. Furthermore, large parts of this glacier consist of firn rather than pure ice, which alters the rheological properties (such as viscosity and compressibility) of the glacier. We present a newly developed, thermo-mechanically coupled, three-dimensional flow model based on the finite-element (FE) modeling software Elmer, and apply it to the Gorshkov crater glacier. By assuming steady-state conditions, the present-day velocity field, temperature field, basal melting rate and age distribution are simulated. We find that flow velocities are generally small (tens of centimeters per year). Horizontal and vertical velocities are of comparable magnitude, which shows that the shallow-ice approximation is not applicable. Owing to the spatially variable volcanic heat flux, the thermal regime at the ice base is cold in the deeper parts of the glacier and temperate in the shallower parts. The measured temperature profile and age horizons at the K2 borehole are reproduced quite well, and remaining discrepancies may be attributed to transient (non-steady-state) conditions. Firn compressibility is identified as a crucial element for the modeling approach.

Prestahnúkur is a 570m high rhyolite glaciovolcanic edifice in Iceland’s Western Rift Zone with a volume of 0.6 km3. Uniform whole rock, mineral and glass compositions suggest that Prestahnúkur was constructed during the eruption of one magma batch. Ar-Ar dating gives an age of 89± 24 ka, which implies eruption during the transition (Oxygen Isotope substages 5d to 5a) between the Eemian interglacial and the Weichselian glacial period. Prestahnu´kur is unique among published accounts of rhyolite tuyas because a base of magmatically-fragmented tephra appears to be absent. Instead, basal exposures consist of glassy lava lobes and coarse hyaloclastite, above which are single and multiple lava sheets with matrix-supported basal breccias and hyaloclastite upper carapaces. Steepening ramp structures at sheet termini are interpreted as ice-contact features. Interactions between erupting magma and water/ice have affected all lithologies. A preliminary model for the construction of Prestahnúkur involves an effusive subglacial eruption between 2–19 years duration which never became emergent, into an ice sheet over 700m thick. If 700m of ice had built up during this interglacial–glacial transition, this would corroborate models arguing for the swift accumulation of land-based ice in rapid response to global cooling.

Glaciological and related observations from 1961 to 2005 at the summit of Mt Wrangell (62.00° N, 144.02°W; 4317m a.s.l.), a massive glacier-covered shield volcano in south-central Alaska, show marked changes that appear to have been initiated by the Great Alaska Earthquake (Mw = 9.2) of 27 March 1964. The 4×6 km diameter, ice-filled Summit Caldera with several post-caldera craters on its rim, comprises the summit region where annual snow accumulation is 1–2m of water equivalent and the mean annual temperature, measured 10 m below the snow surface, is –20°C. Precision surveying, aerial photogrammetry and measurements of temperature and snow accumulation were used to measure the loss of glacier ice equivalent to about 0.03 km3 of water from the North Crater in a decade. Glacier calorimetry was used to calculate the associated heat flux, which varied within the range 20–140Wm–2; total heat flow was in the range 20–100MW. Seismicity data from the crater’s rim show two distinct responses to large earthquakes at time scales from minutes to months. Chemistry of water and gas from fumaroles indicates a shallow magma heat source and seismicity data are consistent with this interpretation.

Mt Veniaminof is a large active stratovolcano located on the Alaska Peninsula (56.2° N, 159° W). We present results of the first geophysical survey of the icefield that fills much of the 10 km×8 km caldera that was most recently modified during the last major eruption roughly 3700 BP. The subglacial topography and ice volume are derived from an 8MHz radio-echo sounding survey conducted in July 2005. Prominent internal reflectors are assumed to be isochronal ash/acid deposits related to local eruptions. Accumulation rates and basal melt rates are calculated using a Nye one-dimensional steady-state accumulation model applied at a location that approximates an ice divide and calibrated by matching internal reflectors with published records of recent local volcanic eruptions. The model yields order of magnitude estimates of the accumulation rate of 4ma–1 water equivalent and 2 ma–1 of basal melt. The subsequent geothermal flux of ∽19Wm–2 is similar to active hydrothermal vents in volcanic lakes. We suggest that these values represent an upper limit for the geothermal flux within the ice-covered regions of the main caldera. We also analyze likely subglacial meltwater flow paths to examine the implications of recent eruption activity at an active intra-caldera cinder cone. Two lava-producing eruptions from the cinder cone in 1983–84 and 1993–94 melted roughly 0.17km3 of ice. The lack of significant deformation of the internal stratigraphy to the south and east of the melt hole suggests that any subglacial drainage in those directions was entirely within subglacial deposits. We suggest that the more likely drainage route was northwest into a large outlet glacier.

The Katla caldera is located under the Mýrdalsjökull ice cap and is one of the most hazardous volcanoes in Iceland due to major jökulhlaups that accompany eruptions. Subglacial geothermal activity is manifested in several 10–50m deep depressions (ice cauldrons) within and at the caldera rim and the total geothermal heat output is of the order of a few hundred megawatts. A short-lived but powerful pulse in geothermal heat output took place in 1999, probably including a minor subglacial eruption, when new ice cauldrons formed in three places and an unexpected jökulhlaup occurred. Following these events, a comprehensive monitoring program was set up for Katla, including ice surface elevation profiling from aircraft, to monitor variations in geothermal heat and detect signs of subglacial water accumulation. A radar altimeter coupled with a kinematic GPS is used, achieving an absolute elevation accuracy of 3m and internal consistency of 1–2 m. Profiles across the caldera are flown twice a year. An annual accumulation-ablation cycle in surface elevation with amplitude of 5–10m is observed. By removing this cycle from the data, changes due to subglacial geothermal activity are obtained. After the events in 1999, a decline in geothermal activity was observed. In 2001–03 some ice cauldrons expanded and deepened by 10–15 m, indicating renewed increase in geothermal activity. This trend is also apparent for 2003–05. The increase in geothermal power amounts to a few tens of megawatts. It is likely that the increased thermal output is related to increased seismicity and caused by magma inflow.

Popocatépetl (19°02’ N, 98°62’W; 5424 m) is one of the largest active stratovolcanoes in the Transmexican Volcanic Belt. A glacier located on the north side has undergone severe ablation since the volcano reinitiated eruptive activity in December 1994. In our study, we calculate the extent of the glacier recession and the loss in glacial mass balance during the period of greatest laharic activity (1994–2002), using photogrammetric treatment of 20 pairs of aerial photographs. The results indicate that from November 1997 to December 2002, the glacier released approximately 3 967 000 m3 of water. A period of intense glacier melting occurred from 4 November 2000 to 15 March 2001 during which time 717 000 m3 of water was released. Much of the melting was attributed to the pyroclastic flow that took place on 22 January 2001 and produced a 14.2 km lahar with 68 000 m3 of water. Among the many types of volcanic events, pyroclastic flows were the most effective in causing sudden snowmelt, although small explosions were also effective since they deposited incandescent material on the glacier. The collapse of the plinian columns covered the glacier with pyroclasts and increased its volume. The existence of control points for georeferencing and a knowledge of the topography underlying the glacier previous to the eruption would have provided more accurate and useful results for hazard prevention.

The very high temperature contrast between magma/ lava and water ice commonly leads to the assumption that significant melting will take place immediately upon magma/ lava ice contact, yet observations of active flows show little evidence of voluminous melting upon contact. We use analytical thermal models to reassess the efficiency with which heat can be transferred from magma to ice in three situations: lava flows erupted on top of glacial ice, sill intrusions beneath glacial ice evolving into subglacial lava flows and dyke intrusions into the interiors of glaciers. We find that the maximum ratios of thickness of ice that can be melted to the thickness of magmatic heat source are likely to be ∽2–5 for subaerial lava flows encroaching onto glaciers, ∽6–7 for subglacial lava flows and ∽10 for dykes intruded into glacial ice. Rates of ice melt production are not linear functions of time and flow thickness, however, and this may account for the observations of minimal immediate water release from beneath advancing lava flows. Field observations during future eruptions should be directed at measuring the temperature of released water.

Simple analytical models of subglacial eruptions are presented, which simulate evolving subglacial cavities and volcanic edifices during rhyolitic eruptions beneath temperate glaciers. They show that the relative sizes of cavity and edifice may strongly influence the eruption mechanisms. Intrusive eruptions will occur if the edifice fills the cavity, with rising magma quenched within the edifice and slow melting of ice. Explosive magma–water interaction may occur if a water- or steam-filled gap develops above the edifice. Meltwater is assumed to drain away continuously, but any gap above the edifice will be filled by meltwater or steam. Ductile roof closure will occur if the glacier weight exceeds the cavity pressure and is modelled here using Nye’s law. The results show that the effusion rate is an important control on the eruption style, with explosive eruptions favoured by large effusion rates. The models are used to explain contrasting eruption mechanisms during various Quaternary subglacial rhyolite eruptions at Torfajökull, Iceland. Although the models are simplistic, they are first attempts to unravel the complex feedbacks between subglacial eruption mechanisms and glacier response that can lead to a variety of eruptive scenarios and associated hazards.

Despite the ubiquity of tunnel channels and valleys within formerly glaciated areas, their origin remains enigmatic. Few modern analogues exist for event-related subglacial erosion. This paper presents evidence of subglacial meltwater erosion and tunnel channel formation during the November 1996 jökulhlaup, Skeiðarárjökull, Iceland. The jökulhlaup reached a peak discharge of 45 000 to 50 000 m3 s–1, with flood outbursts emanating from multiple outlets across the entire 23 km wide glacier snout. Subsequent retreat of the southeast margin of Skeiðarárjökull has revealed a tunnel channel excavated into the surrounding moraine sediment and ascending 11.5m over a distance of 160 m from a larger trough to join the apex of an ice-contact fan formed in November 1996. The tunnel channel formed via hydro-mechanical erosion of 14 000m3 to 24 000 m3 of unconsolidated glacier substrate, evidenced by copious rip-up clasts within the ice-contact fan. Flow reconstruction provides peak discharge estimates of 680±140m3 s–1. The tunnel channel orientation, oblique to local ice flow direction and within a col, suggests that local jökulhlaup routing was controlled by (a) subglacial topography and (b) the presence of a nearby proglacial lake. We describe the first modern example of tunnel channel formation and illustrate the importance of pressurized subglacial jökulhlaup flow for tunnel channel formation.

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.

The catastrophic lahars triggered by the 13 November 1985 eruption of the ice-clad Nevado del Ruiz volcano, Colombia, demonstrate that the interaction of hot pyroclasts with snow and ice can release 30–50 millionm3 of meltwater in 30–90 minutes. The 1985 eruption caused a 16% loss in area and a 9% loss in volume of snow, firn and ice. Turbulent pyroclastic density currents mechanically mixed with snow and produced meltwater at a rate of 0.5–1.6mms–1. Laboratory experiments suggest that turbulent, fluidized pyroclastic density currents exert mechanical and thermal scour, thereby efficiently transferring heat from hot pyroclasts to snow. Ice cap loss at Nevado del Ruiz continued between 1985 and 2000, representing a ∽52% decline in area and a ∽30% fall in volume. Ice 60–190m thick caps the east and southeast summit plateau, whereas an ice field < 30m thick and devoid of snow is retreating on the north, northeast and west edges. This asymmetrical distribution of ice reflects combined long-term effects of the 1985 eruption and of the post-1985 ice cap retreat. Should volcanic activity resume, steep-sided glaciers can fail and pyroclastic flows and surges can sweep the snowpack and generate mixed avalanches and lahars. Although the potential source of meltwater has decreased since 1985, extensive debris at the ice cap margins can be incorporated to future lahars.

The Cordillera Central in Colombia hosts four important glacier-clad volcanoes, namely Nevado del Ruiz, Nevado de Santa Isabel, Nevado del Tolima and Nevado del Huila. Public and scientific attention has been focused on volcano–glacier hazards in Colombia and worldwide by the 1985 Nevado del Ruiz/Armero catastrophe, the world’s largest volcano–glacier disaster. Important volcanological and glaciological studies were undertaken after 1985. However, recent decades have brought strong changes in ice mass extent, volume and structure as a result of atmospheric warming. Population has grown and with it the sizes of numerous communities located around the volcanoes. This study reviews and reassesses the current conditions of and changes in the glaciers, the interaction processes between ice and volcanic activity and the resulting hazards. Results show a considerable hazard potential from Nevados del Ruiz, Tolima and Huila. Explosive activity within environments of snow and ice as well as non-eruption-related mass movements induced by unstable slopes, or steep and fractured glaciers, can produce avalanches that are likely to be transformed into highly mobile debris flows. Such events can have severe consequences for the downstream communities. Integrated monitoring strategies are therefore essential for early detection of emerging activity that may result in hazardous volcano–ice interaction. Corresponding efforts are currently being strengthened within the framework of international programmes.

From radio-echo sounding (RES) surveys and ice core data it can be seen that the ice sheet is melting at the base in a large area in Northern Greenland. The RES images reveal internal layers in the ice. The layers are former deposition surfaces and are thus isochrones. Undulations of the isochrones in regions where the base is smooth suggest that the basal melt rate changes over short distances. This indicates that the geothermal heat flux is very high and has large spatial variability in Northern Greenland. In this study, the basal melt rate at the NorthGRIP drill site in North-Central Greenland is calculated by inverse modelling. We use simple one- and two-dimensional flow models to simulate the ice flow along the NNW-trending ice ridge leading to NorthGRIP. The accumulation is calculated from a dynamical model. Several ice flow parameters are unknown and must be estimated along with the basal melt rate using a Monte Carlo method. The Monte Carlo inversion is constrained by the observed isochrones, dated from the timescale established for the NorthGRIP ice core. The estimates of the basal melt rates around NorthGRIP are obtained from both the one- and two-dimensional models. Combining the estimated basal melt rates with the observed borehole temperatures allows us to convert the basal melt rates to geothermal heat flow values. From the two-dimensional model we find the basal melt rate and geothermal heat flux at NorthGRIP to be 6.1 mma–1 and 129 mWm–2, respectively.

Processes and mechanisms of erosion, transport and deposition within high-magnitude outburst floods such as jökulhlaups and lahars are poorly understood and remain largely unquantified. This study therefore applies a two-dimensional or depth-averaged hydrodynamic model, with fully integrated sediment transport, to reconstruct a Holocene jökulhlaup to have occurred from Kverkfjöll volcano, Iceland. Results indicate simultaneous inundation of multiple channels, flow around islands, hydraulic jumps and multi-directional flow including backwater areas and hydraulic ponding. These flow characteristics are typical of outburst floods that are volcanically triggered, flow through steep volcanic terrain and contain high concentrations of volcaniclastic sediment. Kverkfjöll jökulhlaups had low frontal flow velocities but as stage increased, velocities reached 5–15ms–1. Peak stage was prolonged in zones of hydraulic ponding, but generally attenuated in magnitude and duration downstream. Suspended load transport persisted over the entire hydrograph but bed load transport was spatially discontinuous and comprised distinct pulses. A hierarchy of landforms is proposed, ranging from highest energy zones (erosional gorges, scoured bedrock, cataracts and spillways) to lowest (valley-fills, bars and slackwater deposits). Bedrock erosion was generally where flow exceeded ∽3m flow depth, ∽7ms–1 flow velocity, ∽1×102Nm–2 shear stress and 3 ×102Wm–2 stream power. Deposition occurred below ∽8m flow depth, 11 ms–1 flow velocity, 5 ×102Nm–2 shear stress and 3 ×103Wm–2 stream power. Hydraulic ranges associated with erosion and deposition have considerable overlap due to transitional flow phenomena, transitions in sediment concentration and the influence of upstream effects, such as hydraulic ponding behind topographic constrictions. These results are the first of coupled hydraulic and sediment transport phenomena in high-magnitude outburst floods with fluid rheology and high sediment content, such as jökulhlaups and dilute lahars. Modelled changes in sediment mass closely resembled field-mapped zones of erosion and deposition. This paper therefore introduces a capability to simulate rapid landscape change due to high-magnitude outburst flood.

A dam-break lahar resulting from the last eruption of Ruapehu is expected when Crater Lake reaches a critical high level, probably within the next 1–5 summers. A high level of public consultation, political decision-making and ongoing scientific input has occurred to address the risks. Decisions about managing lahar risk have taken into account the fact that lahars are common on the active Ruapehu Volcano (2797 m) and in valleys draining the mountain and that most lahars are generated by eruptions, about half of which have no useful precursors. Lahars threaten New Zealand’s largest ski area and nationally important infrastructure. Public safety has been the main consideration but the need for long-term risk mitigation and reducing impacts on Tongariro National Park World Heritage Area have also been important. Lahar mitigation on Ruapehu now includes six lahar warning systems, each with active response plans, and some infrastructure isolated from or hardened against lahars.

Ice cauldrons are depressions which form at the surface of ice sheets when an underlying subglacial lake empties, in particular when subglacial volcanic eruptions occur. Notable examples of such cauldrons occur on the surface of the Vatnajökull ice cap in Iceland. More generally, cauldrons will form when a subglacial lake empties during a jökulhlaup. The rate of subsidence of the ice surface is related to the rate at which the subglacial water empties from the lake. We use a viscous version of classical beam theory applied to the ice sheet to determine the relation between the subsidence rate and flood discharge. We use the results to make inferences concerning ring fracture spacings in cauldrons, the consequent effect on flood discharge dynamics and the likely nature of subsidence events in the Antarctic Ice Sheet.