Is the Iceland hot spot also wet? Evidence from the water contents of undegassed submarine and subglacial pillow basalts

Water contents have been measured in basaltic glasses from submarine and subglacial eruption sites along the Reykjanes Ridge and Iceland, respectively, in order to evaluate the hypothesis of Schilling et al. [Phil. Trans. R. Soc. London A 56 (1980) 147-178] that hot spots are also wet spots. Having erupted under pressure the water contents measured in these samples are potentially unaffected by degassing. After correcting these water contents for the effects of crystallisation (to give H₂O(8) values) they indicate that the concentration of water in the source regions increases from 165 ppm at the southern end of the Reykjanes Ridge to between 620 and 920 ppm beneath Iceland. This suggests that Iceland is a wet spot and the H₂O(8) values indicate that its influence on basalt compositions increases northwards along the Reykjanes Ridge from ∼61°N (650 km from the plume centre) towards Iceland. The existence of wetter Icelandic source regions have important implications for mantle melting, as enrichments of this magnitude depress the mantle solidus, increasing the degree of melting at a given temperature. Therefore the enhanced rates of volcanism on Iceland may be a result of wetter sources in addition to a thermal anomaly beneath Iceland.     

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.     

Large values of hydraulic roughness in subglacial conduits during conduit enlargement: implications for modeling conduit evolution

Hydraulic roughness accounts for energy dissipated as heat and should exert an important control on rates of subglacial conduit enlargement by melting. Few studies, however, have quantified how subglacial conduit roughness evolves over time or how that evolution affects models of conduit enlargement. To address this knowledge gap, we calculated values for two roughness parameters, the Darcy–Weisbach friction factor (f) and the Manning roughness coefficient (n), using dye tracing data from a mapped subglacial conduit at Rieperbreen, Svalbard. Values of f and n calculated from dye traces were compared with values of f and n calculated from commonly used relationships between surface roughness heights and conduit hydraulic diameters. Roughness values calculated from dye tracing ranged from 75–0.97 for f and from 0.68–0.09 s m^‐1/3 for n. Equations that calculate roughness parameters from surface roughness heights underpredicted values of f by as much as a factor of 326 and values of n by a factor of 17 relative to values obtained from the dye tracing study. We argue these large underpredictions occur because relative roughness in subglacial conduits during the early stages of conduit enlargement exceeds the 5% range of relative roughness that can be used to directly relate values of f and n to flow depth and surface roughness heights. Simple conduit hydrological models presented here show how parameterization of roughness impacts models of conduit discharge and enlargement rate. We used relationships between conduit relative roughness and values of f and n calculated from our dye tracing study to parameterize a model of conduit enlargement. Assuming a fixed hydraulic gradient of 0.01 and ignoring creep closure, it took conduits 9.25 days to enlarge from a diameter of 0.44 m to 3 m, which was 6–7‐fold longer than using common roughness parameterizations. Copyright © 2013 John Wiley & Sons, Ltd.     

Aluminium in glacial meltwater demonstrates an association with nutrient export (Werenskiöldbreen,Svalbard)

The aluminium (Al) cycle in glacierised basins has not received a great deal of attention in studies of biogeochemical cycles. As Al may be toxic for biota, it is important to investigate the processes leading to its release into the environment. It has not yet been ascertained whether filterable Al (passing through a pore size of 0.45 μm) is incorporated into biogeochemical cycles in glacierised basins. Our study aims to determine the relationship between the processes bringing filterable Al and glacier‐derived filterable nutrients (particularly Fe and Si) into glacierised basins. We investigated the Werenskiöldbreen basin (44.1 km², 60% glacierised) situated in SW Spitsbergen, Svalbard. In 2011, we collected meltwater from a subglacial portal at the glacier front and at a downstream hydrometric station throughout the ablation season. The Al concentration, unchanged between the subglacial system and proglacial zone, reveals that aluminosilicate weathering is a dominant source of filterable Al under subglacial conditions. By examining the Al:Fe ratio compared with pH and the sulphate mass fraction index, we found that the proton source for subglacial aluminosilicate weathering is mainly associated with sulphide oxidation and, to a lesser degree, with hydrolysis and carbonation. In subglacial outflows and in the glacial river, Al and Fe are primarily in the forms of Al(OH)₄^‐ and Fe(OH)₃. The annual filterable Al yield (2.7 mmol m^‐2) was of a magnitude similar to that of nutrients such as filterable Fe (3.0 mmol m^‐2) and lower than that of dissolved Si (18.5 mmol m^‐2). Our results show that filterable Al concentrations in meltwater are significantly correlated to filterable and dissolved glacier‐derived nutrients (Fe and Si, respectively) concentrations in glaciers worldwide. We conclude that a potential bioavailable Al pool derived from glacierised basins may be incorporated in biogeochemical cycles, as it is strongly related to the concentrations and yields of glacier‐derived nutrients.     

Initial results from a distributed,physically based model of glacier hydrology

This paper describes the development and testing of a distributed, physically based model of glacier hydrology. The model is used to investigate the behaviour of the hydrological system of Haut Glacier d'Arolla, Valais, Switzerland. The model has an hourly time-step and three main components: a surface energy balance submodel, a surface flow routing submodel and a subglacial hydrology submodel. The energy balance submodel is used to calculate meltwater production over the entire glacier surface. The surface routing submodel routes meltwater over the glacier surface from where it is produced to where it either enters the subglacial hydrological system via moulins or runs off the glacier surface. The subglacial hydrology submodel calculates water flow in a network of conduits, which can evolve over the course of a melt season simulation in response to changing meltwater inputs. The main model inputs are a digital elevation model of the glacier surface and its surrounding topography, start-of-season snow depth distribution data and meteorological data. Model performance is evaluated by comparing predictions with field measurements of proglacial stream discharge, subglacial water pressure (measured in a borehole drilled to the glacier bed) and water velocities inferred from dye tracer tests. The model performs best in comparison with the measured proglacial stream discharges, but some of the substantial features of the other two records are also reproduced. In particular, the model results show the high amplitude water pressure cycles observed in the borehole in the mid-melt season and the complex velocity/discharge hysteresis cycles observed in dye tracer tests. The results show that to model outflow hydrographs from glacierized catchments effectively, it is necessary to simulate spatial and temporal variations in surface melt rates, the delaying effect of the surface snowpack and the configuration of the subglacial drainage system itself. The model's ability to predict detailed spatial and temporal patterns of subglacial water pressures and velocities should make it a valuable tool for aiding the understanding of glacier dynamics and hydrochemistry. © 1998 John Wiley & Sons, Ltd.     

The tunnel through Torghatten,northern Norway: Landforms and discussion of formation

The tunnel through the mountain of Torghatten, in northern Norway, is generally regarded as a product of wave action. The tunnel is above the late Weichselian marine limit. Fresh looking polished bedrock that resembles subglacial ice-sculptured and meltwater forms, p-forms, occurs near the opening to the landward eastern side of the tunnel and inside. Most likely, the tunnel is a polygenetic formation. Storm action during deglaciations and also subglacial meltwater drainage and plastically sliding ice during glaciations have been active processes in the formation of the tunnel.     

Seasonal reorganization of subglacial drainage inferred from measurements in boreholes

The effect of the formation of a major subglacial drainage channel on the behaviour of the subglacial drainage system of Haut Glacier d'Arolla, Switzerland, was investigated using measurements of borehole water level and the electrical conductivity and turbidity of basal meltwaters. Electrical conductivity profiles were also measured within borehole water columns to identify the water sources driving water level changes, and to determine patterns of water circulation in boreholes. Prior to channel formation, boreholes showed idiosyncratic and poorly coordinated behaviour. Diurnal water level fluctuations were small and driven by supraglacial/englacial water inputs, even when boreholes were connected to a subglacial drainage system. This system appeared to consist of hydraulically impermeable patches interspersed with storage spaces, and transmitted a very low water flux. Drainage reorganization, which occurred around 31 July, 1993, in response to rapidly rising meltwater and rainfall inputs, seems to have involved the creation of a connection between an incipient channel and a well-established channelized system located further down-glacier. Once a major channel existed within the area of the borehole array, borehole water level fluctuations were forced by discharge-related changes in channel water pressure, although a diversity of responses was observed. These included (i) synchronous, (ii) damped and lagged, (iii) inverse, and (iv) alternating inverse/lagged responses. Synchronous responses occurred in boreholes connected directly to the channel, while damped and lagged responses occurred in boreholes connected to it by a more resistive drainage system. Pressure variations within the channel resulted in diurnal transfer of mechanical support for the ice overburden between connected and unconnected areas of the bed, producing inverse and alternating patterns of water level response. © 1998 John Wiley & Sons, Ltd.     

Glacial curvilineations found along the southern sector of the Laurentide Ice sheet and a hypothesis of formation involving subglacial slope failure in tunnel valleys and subglacial lakes

First discovered in Poland, glacial curvilineations (GCLs) are enigmatic landforms comprising parallel sets of sinuous ridges and troughs of metres amplitude and around 150 m wavelength, found within kilometres‐wide valleys interpreted as being produced by meltwater flowing subglacially. Their morphological and sedimentary characteristics and association with tunnel valleys has been described for some prominent Polish examples. From these observations the existing hypothesis is that they form as a consequence of erosion by longitudinal vortices that develop in subglacial floods. Here we report, for the first time, GCLs found along the southern sector of the Laurentide Ice Sheet in three northern states of the USA. Using mapping and topographic analysis from high resolution digital elevation models we report observations on their morphological properties and landform associations. We find aspects of their context and morphology difficult to explain using the existing hypothesis. We instead suggest that these glacial curvilineations are produced by subglacial bank and slope failures that locally widen tunnel valleys, or that occur near subglacial lake shorelines. Further investigation is required to test this hypothesis and to ascertain the mechanisms of proposed mass movements, which may have occurred by rotational or translational slope failure or by creep deformation. Our preferred mechanism is that such movements occurred where subglacial water was emplaced over previously perma‐frozen ground. Under such circumstances, sediment blocks thawed by the water may then easily glide over a frozen décollement at low slope angles; analogous to subaerial active‐layer glides in permafrost environments. Permafrost spring sapping may have provided lines of weakness for slope failure. If the requirement for permafrost is found to hold, then GCLs may become an important indicator of the palaeo‐distribution of permafrost. Copyright © 2017 John Wiley & Sons, Ltd.     

The emplacement of an obsidian dyke through thin ice: Hrafntinnuhryggur, Krafla Iceland

An eruption along a 2.5 km-long rhyolitic dyke at Krafla volcano, northern Iceland during the last glacial period formed a ridge of obsidian (Hrafntinnuhryggur). The ridge rises up to 80 m above the surrounding land and is composed of a number of small-volume lava bodies with minor fragmental material. The total volume is < 0.05 km³. The lava bodies are flow- or dome-like in morphology and many display columnar-jointed sides typical of magma–ice interaction, quench-fragmented lower margins indicative of interaction with meltwater and pumiceous upper surfaces typical of subaerial obsidian flows. The fragmental material compromises poorly-sorted perlitic quench hyaloclastites and poorly-exposed pumiceous tuffs. Lava bodies on the western ridge flanks are columnar jointed and extensively hydrothermally altered. At the southern end of the ridge the feeder dyke is exposed at an elevation  95 m beneath the ridge crest and flares upwards into a lava body.Using the distribution of lithofacies, we interpret that the eruption melted through ice only 35–55 m thick, which is likely to have been dominated by firn. Hrafntinnuhryggur is therefore the first documented example of a rhyolitic fissure eruption beneath thin ice/firn. The eruption breached the ice, leading to subaerial but ice/firn-contact lava effusion, and only minor explosive activity occurred. The ridge appears to have been well-drained during the eruption, aided by the high permeability of the thin ice/firn, which appears not to have greatly affected the eruption mechanisms. We estimate that the eruption lasted between 2 and 20 months and would not have generated a significant jökulhlaup (< 70 m³ s^− 1).     

The hydrochemistry of meltwaters draining a polythermal-based,high Arctic glacier,south Svalbard: I. The ablation season

Solute and runoff time-series at Finsterwalderbreen, Svalbard, provide evidence for considerable basal routing of water and the existence of at least two contrasting subglacial chemical weathering environments. The hydrochemistry of a subglacial upwelling provides evidence for a snowmelt-fed subglacial reservoir that dominates bulk runoff during recession flow. High concentrations of Cl⁻ and crustal ions, high pCO₂ and ratios of [*SO²⁻₄/(*SO²⁻₄+HCO⁻₃)] close to 0·5 indicate the passage of snowmelt through a subglacial weathering environment characterized by high rock:water ratios, prolonged residence times and restricted access to the atmosphere. At higher discharges, bulk runoff becomes dominated by icemelt from the lower part of the glacier that is conveyed through a chemical weathering environment characterized by low rock:water ratios, short residence times and free contact with atmospheric gases. These observations suggest that icemelt is routed via a hydrological system composed of basal/ice-marginal, englacial and supraglacial components and is directed to the glacier margins by the ice surface slope. Upwelling water flows relatively independently of icemelt to the terminus via a subglacial drainage system, possibly constituting flow through a sediment layer. Cold basal ice at the terminus forces it to take a subterranean routing in its latter stages. The existence of spatially discrete flow paths conveying icemelt and subglacial snowmelt to the terminus may be the norm for polythermal-based glaciers on Svalbard. Proglacial mixing of these components to form the bulk meltwaters gives rise to hydrochemical trends that resemble those of warm-based glaciers. These hydrochemical characteristics of bulk runoff have not been documented on any other glacier on Svalbard to date and have significance for understanding interactions between thermal regime and glacier hydrology. © 1998 John Wiley & Sons, Ltd.