Geomorphological impact and morphodynamic effects on flow conveyance of the 1999 jökulhlaup at sólheimajökull,Iceland

The 1999 jökulhlaup at Sólheimajökull was the first major flood to be routed through the proglacial system in over 600 years. This study reconstructed the flood using hydrodynamic, sediment transport and morphodynamic numerical modelling informed by field surveys, aerial photograph and digital elevation model analysis. Total modelled sediment transport was 469 800 m³ (+/‐ 20%). Maximum erosion of 8.2 m occurred along the ice margin. Modelled net landscape change was –86 400 m³ (+/‐ 40%) resulting from –275 400 m³ (+/‐ 20%) proglacial erosion and 194 400 m³ (+/‐ 20%) proglacial deposition. Peak erosion rate and peak deposition rate were 650 m³ s^‐1 (+/‐ 20%) and 595 m³ s^‐1 (+/‐ 20%), respectively, and coincided with peak discharge of water at 1.5 h after flood initiation. The pattern of bed elevation change during the rising limb suggested widespread activation of the bed, whereas more organisation, perhaps primitive bedform development, occurred during the falling limb. Contrary to simplistic conceptual models, deposition occurred on the rising stage and erosion occurred on the falling limb. Comparison of the morphodynamic results with a hydrodynamic simulation illustrated effects of sediment transport and bed elevation change on flow conveyance. The morphodynamic model advanced flood arrival and peak discharge timings by 100% and 19%, respectively. However, peak flow depth and peak flow velocity were not significantly affected. We suggest that morphodynamic processes not only increase flow mass and momentum but that they also introduce a feedback process whereby flood conveyance becomes more efficient via erosion of minor bed protrusions and deposition that infills or subdues minor bed hollows. A major implication of this study is that reconstructions of outburst floods that ignore sediment transport, such as those used in interpretation of long‐term hydrological record and flood risk assessments, may need considerable refinement. Copyright © 2015 John Wiley & Sons, Ltd.     

Glacier retreat and landform production on an overdeepened glacier foreland: the debris‐charged glacial landsystem at Kvíárjökull,Iceland

Glacier recession and landform development in a debris‐charged glacial landsystem characterized by an overdeepening is quantified using digital photogrammetry, digital elevation model (DEM) construction and mapping of the Icelandic glacier Kvíárjökull for the period 1945–2003. Melting of ice‐cores is recorded by surface lowering rates of 0·8 m yr^–1 (1945–1964), 0·3 m yr^–1 (1964–1980), 0·015 m yr^–1 (1980–1998) and 0·044 m yr^–1 (1998–2003). The distribution/preservation of pushed and stacked ice‐cored moraine complexes are determined by the location of the long‐term glacial drainage network in combination with retreat from the overdeepening, into which glacifluvial sediment is being directed and where debris‐rich ice masses are being reworked and replaced by esker networks produced in englacial meltwater pathways that bypassed the overdeepening and connected to outwash fans prograding over the snout. Recent accelerated retreat of Kvíárjökull, potentially due to increased mass balance sensitivity, has made the snout highly unstable, especially now that the overdeepening is being uncovered and the snout flooded by an expanding pro‐glacial, and partially supraglacial, lake. This case study indicates that thick sequences of debris‐charged basal ice/controlled moraine have a very low preservation potential but ice‐cored moraine complexes can develop into hummocky moraine belts in de‐glaciated terrains because they are related to the process of incremental stagnation, which at Kvíárjökull has involved periodic switches from transport‐dominant to ablation‐dominant conditions. Glacier recession is therefore recorded temporally and spatially by two suites of landforms relating to two phases of landform production which are likely typical for glaciers occupying overdeepenings: an early phase of active, temperate recession recorded by push moraines and lateral moraines and unconfined pro‐glacial meltwater drainage; and a later phase of incremental stagnation and pitted outwash head development initiated by the increasing topographic constraints of the latero‐frontal moraine arc and the increasing importance of the overdeepening as a depo‐centre. Copyright © 2012 John Wiley & Sons, Ltd.     

Sensitivity analysis of temperature‐index melt simulations to near‐surface lapse rates and degree‐day factors at Vestari‐Hagafellsjökull,Langjökull,Iceland

The performance of temperature‐index melt models is particularly affected by the choice of near‐surface lapse rate used to determine the sum of positive daily temperatures at different elevations, and by the choice of factor used to relate this sum to the rate of melting. Data from the Langjökull ice cap are used in this study to quantify the influence of lapse‐rate and degree‐day factor variation on temperature‐index melt simulations. The lapse rate was significantly lower during summer than in spring or autumn, as a result of diabatic cooling, reducing boundary‐layer sensitivity to free‐air temperature change. The summer lapse rate was also significantly lower than the saturated adiabatic lapse rate. A sensitivity of approximately 600 mm water equivalent (w.e.) cumulative June–August melt per 0.1 °C 100 m^–1 change in lapse rate was found across a 500‐m altitude range. The sensitivity to a 1‐mm w.e. °C^–1 day^–1 change in degree‐day factors varied more: from approximately 500 mm w.e. cumulative summer melt at low elevation to approximately 200 mm w.e. at high elevation, reflecting the decline in melt rates associated with the greater persistence of snow with increasing altitude. The determination of a degree‐day factor for snow is complicated by the densification of the ageing snowpack, but the application of a parameterization for near‐surface density on the basis of albedo helped account for the development of snow water equivalence. Lapse rate was parameterized as a function of standardized anomalies in 750 hPa reanalysis temperature and significantly improved the simulation of cumulative summer melt compared with models applying the saturated adiabatic lapse rate. Copyright © 2012 John Wiley & Sons, Ltd.     

Supra‐glacial deposition and flux of catastrophic rock–slope failure debris,south‐central Alaska

The ongoing debate over the effects of global environmental change on Earth's cryosphere calls for detailed knowledge about process rates and their variability in cold environments. In this context, appraisals of the coupling between glacier dynamics and para‐glacial erosion rates in tectonically active mountains remain rare. We contribute to filling this knowledge gap and present an unprecedented regional‐scale inventory of supra‐glacial sediment flux and hillslope erosion rates inferred from an analysis of 123 large (> 0·1 km²) catastrophic bedrock landslides that fell onto glaciers in the Chugach Mountains, Alaska, as documented by satellite images obtained between 1972 to 2008. Assuming these supra‐glacial landslide deposits to be passive strain markers we infer minimum decadal‐scale sediment yields of 190 to 7400 t km^–2 yr^–1 for a given glacier‐surface cross‐section impacted by episodic rock–slope failure. These rates compare to reported fluvial sediment yields in many mountain rivers, but are an order of magnitude below the extreme sediment yields measured at the snouts of Alaskan glaciers, indicating that the bulk of debris discharged derives from en‐glacial, sub‐glacial or ice‐proximal sources. We estimate an average minimum para‐glacial erosion rate by large, episodic rock–slope failures at 0·5–0·7 mm yr^–1 in the Chugach Mountains over a 50‐yr period, with earthquakes likely being responsible for up to 73% of this rate. Though ranking amongst the highest decadal landslide erosion rates for this size of study area worldwide, our inferred rates of hillslope erosion in the Chugach Mountains remain an order of magnitude below the pace of extremely rapid glacial sediment export and glacio‐isostatic surface uplift previously reported from the region. Copyright © 2012 John Wiley & Sons, Ltd.     

Englacial and subglacial water flow at Skálafellsjökull,Iceland derived from ground penetrating radar,in situ Glacsweb probe and borehole water level measurements

We reconstruct englacial and subglacial drainage at Skálafellsjökull, Iceland, using ground penetrating radar (GPR) common offset surveys, borehole studies and Glacsweb probe data. We find that englacial water is not stored within the glacier (water content ~0–0.3%). Instead, the glacier is mostly impermeable and meltwater is able to pass quickly through the main body of the glacier via crevasses and moulins. Once at the glacier bed, water is stored within a thin (1 m) layer of debris‐rich basal ice (2% water content) and the till. The hydraulic potential mapped across the survey area indicates that when water pressures are high (most of the year), water flows parallel to the margin, and emerges 3 km down glacier at an outlet tongue. GPR data indicates that these flow pathways may have formed a series of braided channels. We show that this glacier has a very low water‐storage capacity, but an efficient englacial drainage network for transferring water to the glacier bed and, therefore, it has the potential to respond rapidly to changes in melt‐water inputs. © 2015 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Pro‐glacial lake sedimentation from jökulhlaups (GLOF), Blåmannsisen,northern Norway

Outburst floods from glacier‐dammed lakes are major events associated with glacier thinning and volume reduction. This paper investigates jökulhlaups emanating from the glacier‐dammed lake Øvre Messingmalmvatn at Rundvassbreen, an outlet glacier of the Blåmannsisen ice cap in northern Norway. Since 2001, the lake has several times been observed to drain suddenly, causing jökulhlaup outbursts into the pro‐glacial lake Rundvatnet. Varve analysis and lead‐210 (²¹⁰Pb) dating were used to date sediment cores taken from Rundvatnet. It was found that sedimentation from jökulhlaups is recognizable in the lake as distinct sand layers embedded in the varved silt‐clay sequence which represents the normal lake sedimentation. Sand fractions were carried in suspension because of the extreme hydraulic conditions of jökulhlaups. The thickest sand layer was deposited during the 2001 jökulhlaup which lasted three days and had a total volume of 40 ×10⁶ m³. Jökulhlaups were also recorded in 2005, 2007, 2009, and 2010; they each resulted in a sand layer. Annual sediment accumulation in Rundvatnet increased up to 10‐fold during the years with jökulhlaup outburst floods, from a normal value of 1–2 mm yr^?1 to 8–10 mm yr^?1. Five other jökulhlaups were identified from the 1910–1930 sedimentation interval, in addition to those observed in 2001–2010; there appear to have been none for 70 years during 1931–2000. Each jökulhlaup was preceded by a period when the glacier thinned to a critical volume and could no longer withstand the hydrostatic pressure of Øvre Messingmalmvatn; consequently a tunnel developed beneath the glacier, leading to a jökulhlaup. Statistical analyses of the correlations between the pro‐glacial sedimentation rate and temperature and precipitation suggested that although climate conditions are expected to influence sedimentation in the pro‐glacial catchment, a host of other interacting factors moderate the availability and delivery of sediment to the pro‐glacial system, making the processes responsible for changes in pro‐glacial sedimentation to remain uncertain. Copyright © 2014 John Wiley & Sons, Ltd.     

Effects of helical flow in one‐dimensional modelling of sediment distribution at river bifurcations

Complex flow processes at river bifurcations and the influence of the layout of a bifurcation make it difficult to predict sediment distribution over the downstream branches in case bedload transport dominates. In one‐dimensional models we need a nodal point relationship that prescribes the distribution of sediment over the downstream branches. We have identified which factors need to be included in such a relationship for the division of bedload transport at bifurcations. Next, irrotational flow theory for idealized geometries has been used to derive a simple physics‐based nodal point relationship that accounts for the effects of helical flow in the situation that a channel takes off under an angle from a straight main channel. This first step towards a complete nodal point relationship is applicable to bedload transport situations if the flow is clearly curved and if there is no pronounced bed topography. The relationship has been tested against data from a unique set of laboratory measurements, numerical data and data from a scale model of the Rhine bifurcation at Pannerden in the Netherlands. We find that the derived model yields a reasonable prediction of the sediment division over the downstream branches, and yields better predictions than the Wang et al. model for the situation considered. Considering the relative complexity and limited accuracy of the nodal point relationship for the effect of helical flow alone, however, we conclude thatderiving a practical physics‐based 1‐D relationship including all relevant processes is not feasible. We therefore recommend 2‐D or 3‐D modelling for all cases in general where morphological evolution depends on the division of bedload transport at bifurcations. Copyright © 2012 John Wiley & Sons, Ltd.     

Hydro‐meteorological drivers and sources of suspended sediment flux in the pro‐glacial zone of the retreating Castle Creek Glacier,Cariboo Mountains,British Columbia,Canada

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.     

Post‐glacial range of variability in the Conejos River Valley,southern Colorado,USA: fluvial response to climate change and sediment supply

Effective river management strategies require an understanding of how fluvial processes vary both spatially and temporally. Here, we examine the natural range of variability in the Conejos River Valley, southern Colorado, through documentation of terrace morphostratigraphic and sedimentological characteristics as well as through investigation of sediment contributions from headwaters, hillslopes and tributary streams. Additionally, soil development and radiocarbon ages, together with local and regional paleoclimate reconstructions, were used to infer the range of processes acting in this system. Since de‐glaciation, the Conejos River has fluctuated between episodes of bedrock strath formation, aggradation and vertical incision. Morphostratigraphic relationships, soil development and radiocarbon ages enable us to propose a chronology for periods of alluvial deposition (around 8·9–7·6 ka, 5·5 ka and from 3·5 to 1·1 ka), separated by intervals of fluvial incision. We infer potential forcing mechanisms by utilizing multiple working hypotheses. Specifically, we discuss the potential for increases in sediment supply during periods of (1) para‐glacial adjustment, (2) climatic cooling, (3) increased frequency of climate change and (4) increased fire frequency or severity. We also consider the effects of changes in stream discharge and extreme storm occurrence. We conclude that combinations of these processes, operating at different times, have contributed to sediment mobilization since de‐glaciation. Stream and landform morphology also varies longitudinally due to the influence of remnant glacial topography. In particular, valley bottom overdeepening at tributary junctions has resulted in incision and strath formation into unlithified glacial deposits (i.e. fill‐cut terraces) rather than bedrock in some reaches. Overall, the Conejos fluvial system has varied significantly both temporally and spatially since de‐glaciation and appears to be sensitive to changes in sediment supply related to Holocene scale climate fluctuations. This natural range of variability must therefore be a key consideration in any future stream management policies. Copyright © 2012 John Wiley & Sons, Ltd.     

Controls on the location,morphology and evolution of complex esker systems at decadal timescales,Breiðamerkurjökull,southeast Iceland

This paper uses detailed mapping of eskers to address three questions which are important for reconstructing meltwater behaviour beneath contemporary and ancient ice masses: ‘What controls the morphology of simple and complex esker systems?’, ‘How do esker systems evolve through time?’ and ‘Are esker patterns compatible with groundwater controlled hydraulic spacing of esker tunnels?’. Esker crestlines and widths are mapped on the Breiðamerkurjökull foreland for eight time slices between 1945 and 2007, from high resolution (~50 cm) aerial photography, permitting their long‐term morphological evolution to be analysed in a high level of detail. We find that complex eskers develop where meltwater and sediment is abundant, such that sediment clogs channels, forming distributary eskers. Isolated eskers are simpler and smaller and reflect less abundant meltwater and sediment, which is unable to clog channels. Eskers may take several decades to emerge from outwash deposits containing buried ice and can increase or decrease in size when ice surrounding and underlying them melts out. It has been suggested that groundwater–channel coupling dictates the spacing between eskers at Breiðamerkurjökull. Our results do not dispute this, but suggest that the routing of sediment and meltwater through medial moraines is an additional important control on esker location and spacing. These results may be used to better understand the processes surrounding esker formation in a variety of geographical settings, enabling a more detailed understanding of the operation of meltwater drainage systems in sub‐marginal zones beneath glaciers and ice sheets. Copyright © 2015 John Wiley & Sons, Ltd.     

A multi‐temporal analysis of AMSR‐E data for flood and discharge monitoring during the 2008 flood in Iowa

The objective of this work is to demonstrate the potential of using passive microwave data to monitor flood and discharge conditions and to infer watershed hydraulic and hydrologic parameters. The case study is the major flood in Iowa in summer 2008. A new Polarisation Ratio Variation Index (PRVI) was developed based on a multi‐temporal analysis of 37 GHz satellite imagery from the Advanced Microwave Scanning Radiometer (AMSR‐E) to calculate and detect anomalies in soil moisture and/or inundated areas. The Robust Satellite Technique (RST) which is a change detection approach based on the analysis of historical satellite records was adopted. A rating curve has been developed to assess the relationship between PRVI values and discharge observations downstream. A time‐lag term has been introduced and adjusted to account for the changing delay between PRVI and streamflow. Moreover, the Kalman filter has been used to update the rating curve parameters in near real time. The temporal variability of the b exponent in the rating curve formula shows that it converges toward a constant value. A consistent 21‐day time lag, very close to an estimate of the time of concentration, was obtained. The agreement between observed discharge downstream and estimated discharge with and without parameters adjustment was 65 and 95%, respectively. This demonstrates the interesting role that passive microwave can play in monitoring flooding and wetness conditions and estimating key hydrologic parameters. Copyright © 2011 John Wiley & Sons, Ltd.     

Wavelet analysis of inter‐annual variability in the runoff regimes of glacial and nival stream catchments,Bow Lake,Alberta

Continuous wavelet analyses of hourly time series of air temperature, stream discharge, and precipitation are used to compare the seasonal and inter‐annual variability in hydrological regimes of the two principal streams feeding Bow Lake, Banff National Park, Alberta: the glacial stream draining the Wapta Icefields, and the snowmelt‐fed Bow River. The goal is to understand how water sources and flow routing differ between the two catchments. Wavelet spectra and cross‐wavelet spectra were determined for air temperature and discharge from the two streams for summers (June–September) 1997–2000, and for rainfall and discharge for the summers of 1999 and 2000. The diurnal signal of the glacial runoff was orders of magnitude higher in 1998 than in other years, indicating that significant ice exposure and the development of channelized glacial drainage occurred as a result of the 1997–98 El Niño conditions. Early retreat of the snowpack in 1997 and 1998 led to a significant summer‐long input of melt runoff from a small area of ice cover in the Bow River catchment; but such inputs were not apparent in 1999 and 2000, when snow cover was more extensive. Rainfall had a stronger influence on runoff and followed quicker flow paths in the Bow River catchment than in the glacial catchment. Snowpack thickness and catchment size were the primary controls on the phase relationship between temperature and discharge at diurnal time scales. Wavelet analysis is a fast and effective means to characterize runoff, temperature, and precipitation regimes and their interrelationships and inter‐annual variability. The technique is effective at identifying inter‐annual and seasonal changes in the relative contributions of different water sources to runoff, and changes in the time required for routing of diurnal meltwater pulses through a catchment. However, it is less effective at identifying changes/differences in the type of the flow routing (e.g. overland flow versus through flow) between or within catchments. Copyright © 2003 John Wiley & Sons, Ltd.     

Impacts of post‐glacial rebound on landslide spatial distribution at a regional scale in northern Iceland (Skagafjörður)

In de‐glaciated areas, para‐glaciation (i.e. the conditioning of landscapes by prior glaciation) has often been considered a major predisposing factor in landslide occurrence; its consequences have been particularly well identified at a fine scale (especially on bedrock jointing). Hitherto, the relative impacts of para‐glaciation on hillslope dynamics at a regional scale had nevertheless not been quantified statistically. We examine Skagafjörður area (northern Iceland) where landslides are widespread (at least 108 were mapped in an area of c. 3000 km²). We compare the role of para‐glaciation (debuttressing, influence of post‐glacial rebound) with that of classic factors (topography, lithology, etc.) in landslide occurrence and location, using a spatial analysis based on a chi‐square test. On the one hand, the results highlight that landslides are over‐represented in areas where post‐glacial rebound was at its maximum, with a stronger concentration of landslides in the northern part of the fjord. On the other hand, the distribution of landslides did not show any clear relationship with the pattern of glacial debuttressing. Tschuprow coefficient highlights that the influence of post‐glacial rebound on landslide location is higher than the combined influence of slope gradient, curvature or geological structure. This result is supported by our initial evidence for the timing of landslides in the area: most landslides occurred during the first half of the Holocene, and a period of hillslope instability was initiated when the post‐glacial uplift was at its maximum. Finally, the mechanisms that link post‐glacial rebound and landsliding as well as the geomorphic impacts of landslides, are discussed. Copyright © 2013 John Wiley & Sons, Ltd.     

Suspended sediment and total dissolved solid yield patterns at the headwaters of Urumqi River,northwestern China: a comparison between glacial and non‐glacial catchments

In order to understand the differences in the suspended sediment and total dissolved solid (TDS) yield patterns between the glacial and non‐glacial catchments at the headwaters of Urumqi River, northwestern China, water samples were collected from a glacier catchment and an empty cirque catchment within the region, during three melting seasons from 2006 to 2008. These samples were analyzed to estimate suspended sediment and TDS concentrations, fluxes and erosion rates in the two adjoining catchments. There were remarked differences in suspended sediment and TDS yield patterns between the two catchments. Suspended sediment concentrations were controlled mainly by the sediment source, whereas TDS concentrations were primarily related to the hydrologic interaction with soil minerals. Generally, the glacial catchment had much higher suspended sediment and TDS yields, together with higher denudation rates, than the non‐glacial catchment. Overall, glacial catchment was mainly dominated by physical denudation process, whereas the non‐glacial catchment was jointly influenced by physical and chemical denudation processes. The observed differences in material delivery patterns were mainly controlled by the runoff source and the glacial processes. The melting periods of glacier and snow were typically the most important time for the suspended sediment and TDS yields. Meanwhile, episodic precipitation events could generate disproportionately large yields. Subglacial hydrology dynamics, glaciers pluck and grind processes could affect erodibility, and the large quantities of dust stored on the glacier surface provided additional sources for suspended sediment transport in the glacial catchment. These mechanisms imply that, in response to climate change, the catchment behaviour will be modified significantly in this region, in terms of material flux. Copyright © 2013 John Wiley & Sons, Ltd.     

Subglacial to emergent basaltic volcanism at Hlöðufell, south-west Iceland: A history of ice-confinement

Hlöðufell is a familiar 1186 m high landmark, located about 80 km northeast of Reykjavík, and 9 km south of the Langkjökull ice-cap in south-west Iceland. This is the first detailed study of this well-exposed and easily accessible subglacial to emergent basaltic volcano. Eight coherent and eleven volcaniclastic lithofacies are described and interpreted, and its evolution subdivided into four growth stages (I–IV) on the basis of facies architecture. Vents for stages I, II, and IV lie along the same fissure zone, which trends parallel to the dominant NNE–SSW volcano-tectonic axis of the Western Volcanic Zone in this part of Iceland, but the stage III vent lies to the north, and is probably responsible for the present N–S elongation of the volcano. The basal stage (I) is dominated by subglacially erupted lava mounds and ridges, which are of 240 m maximum thickness, were fed from short fissures and locally display lava tubes. Some of the stage I lavas preserve laterally extensive flat to bulbous, steep, glassy surfaces that are interpreted to have formed by direct contact with surrounding ice, and are termed ice-contact lava confinement surfaces. These surfaces preserve several distinctive structures, such as lava shelves, pillows that have one flat surface and mini-pillow (< 10 cm across) breakouts, which are interpreted to have formed by the interplay of lava chilling and confinement against ice, ice melting and ice fracture. The ice-contact lava confinement surfaces are also associated with zones of distinctive open cavities in the lavas that range from about 1 m to several metres across. The cavities are interpreted as having arisen by lava engulfing blocks of ice, that had become trapped in a narrow zone of meltwater between the lava and the surrounding ice, and are termed ice-block meltout cavities. The same areas of the lavas also display included and sometimes clearly rotated blocks of massive to planar to cross-stratified hyaloclastite lapilli tuffs and tuff–breccias, termed hyaloclastite inclusions, which are interpreted as engulfed blocks of hyaloclastite/pillow breccia carapace and talus, or their equivalents reworked by meltwater. Some of the stage I lavas are mantled at the southern end of the mountain by up to 35 m thickness of well-bedded vitric lapilli tuffs (stage II), of phreatomagmatic origin, which were erupted from a now dissected cone, preserved in this area. The tephra was deposited dominantly by subaqueous sediment gravity flows (density currents) in an ice-bound lake (or less likely a sub-ice water vault), and was also transported to the south by sub-ice meltwater traction currents. This cone is onlapped by a subaerial pahoehoe lava-fed delta sequence, formed during stage III, and which was most likely fed from a now buried vent(s), located somewhere in the north-central part of the mountain. A 150 m rise in lake level submerged the capping lavas, and was associated with progradation of a new pahoehoe lava-fed delta sequence, produced during stage IV, and which was fed from the present summit cone vent. The water level rise and onset of stage IV eruptions were not associated with any obviously exposed phreatomagmatic deposits, but they are most likely buried beneath stage IV delta deposits. Stage IV lava-fed deltas display steep benches, which do not appear to be due to syn- or post-depositional mass wasting, but were probably generated during later erosion by ice. The possibility that they are due to shorter progradation distances than the underlying stage III deltas, due to ice-confinement or lower volumes of supplied lava is also considered.