Distributed Drainage System Research Articles

Channelized, distributed, and disconnected: spatial structure and temporal evolution of the subglacial drainage under a valley glacier in the Yukon

Abstract. The subglacial drainage system is one of the main controls on basal sliding but remains only partially understood. Here we expand the analysis of the 8-year dataset of borehole observations on a small, alpine polythermal valley glacier in the Yukon Territory. We presented this dataset in Rada and Schoof (2018), where we described the seasonal evolution of the drainage system and underlined the importance of hydraulic isolation at the glacier bed. These borehole observations constitute a unique dataset, both due to the length of the records and the density of the observations, with up to 157 simultaneously working pressure sensors. Now, to explore the spatial structure of the drainage system and its seasonal progression, we automatically cluster boreholes based on similarities in their water pressure records and follow their evolution through the melt season. Some of these borehole clusters show water pressure variations that suggest they are part of a drainage system connected to the surface meltwater supply, while others show features consistent with hydraulic isolation. The distribution of connected and isolated boreholes suggests that the distributed drainage system we observe comprises a network of small conduits with spacings smaller than the borehole bottom diameter (approximately 25–50 cm). Within these hydraulically connected areas, pressure phase lags, and amplitude attenuation rarely shows the behaviour expected in a diffusive system. This observation suggests that the diffusivity distribution in such areas presents a fine structure at scales smaller than our minimum borehole spacing of 15 m. However, at a glacier-wide scale, we observe that hydraulic connections are ubiquitous in some regions of the bed and permanently absent in others, suggesting large contrasts in diffusivity. Within disconnected areas, boreholes often show small-amplitude water pressure variations associated with horizontal normal stress transfers. Such stress transfers seem to play a more important role than previously considered for controlling the effective pressure distribution at the bed. Through the melt season, the evolution of borehole clusters suggests that the diurnal meltwater supply promotes the growth of the low-efficiency drainage systems found early in the season while stimulating the shrinkage and fragmentation of the more efficient drainage systems that appear later in the season. Therefore, an increase in drainage efficiency is associated with the growth of disconnected areas. Our observations support the traditional view of a distributed drainage system early in the melt season that gradually evolves into a progressively more channelized system. However, the most notable difference is the highly heterogeneous distribution of diffusivity that our results suggest and the robust support for disconnected areas. The extent of disconnected areas could be an essential control of basal speed variations. It is possible that even relatively small disconnected areas could have a disproportionate effect on basal speed.

Moulin Density Controls the Timing of Peak Pressurization Within the Greenland Ice Sheet's Subglacial Drainage System

AbstractLinks between hydrology and sliding of the Greenland Ice Sheet (GrIS) are poorly understood. Here, we monitored meltwater's propagation through the glacial hydrologic system for catchments at different elevations by quantifying the lag cascade as daily meltwater pulses traveled through the supraglacial, englacial, and subglacial drainage systems. We found that meltwater's residence time within supraglacial catchments—depending upon area, snow cover, and degree of channelization—controls the timing of peak moulin head, resulting in the 2 hr later peak observed at higher elevations. Unlike at lower elevations where peak moulin head and peak sliding coincided, at higher elevations peak sliding lagged peak moulin head by ∼2.8 hr. This delay was likely caused by the area's lower moulin density, which required diurnal pressure oscillations to migrate further into the distributed drainage system to elicit the observed velocity response. These observations highlight the supraglacial drainage system's control on coupling GrIS subglacial hydrology and sliding.

Assessing Future Rainfall Intensity–Duration–Frequency Characteristics across Taiwan Using the k-Nearest Neighbor Method

This study analyzes the changes in rainfall intensities across Taiwan using the k-Nearest Neighbor method. Biases are corrected according to the identified discrepancy between the probability distribution of the model run and that of the observed data in the historical period. The projections of 21 weather stations in Taiwan under 10 (2RCP × 5GCM) scenarios for the near-(2021–2040) and far-future (2081–2100) are derived. The frequently occurred short-duration storm events in some regions decrease, but they are still vulnerable to flood considering the existing drainage capacities. More specifically, the land-subsidence region in the central, the landslide-sensitive mountainous region in the north and central, the pluvial- and fluvial-flood prone region in the north, and the eastern regions with vulnerable infrastructures should be especially aware of long-duration extreme events. Associations of the rainfall intensity with the different return period as well as the duration are further analyzed. The short-duration extreme events will become stronger, especially for 1-h events in the northern region and 1 or 2-h events in both the southern and eastern regions. In addition, places without experiences of long-lasting events may experience rainfall amounts exceeding 500 mm should be alert. Adaptation measures such as establishing distributed drainage system or renewing hydrological infrastructures in the eastern region are suggested considering the near future projection, and in the central and the southern regions for far future as well. Our findings can assist adaptation-related decision-making for more detailed stormwater/water resource management.

An analysis of instabilities and limit cycles in glacier-dammed reservoirs

Abstract. Glacier lake outburst floods are common glacial hazards around the world. How big such floods can become (either in terms of peak discharge or in terms of total volume released) depends on how they are initiated: what causes the runaway enlargement of a subglacial or other conduit to start the flood, and how big can the lake get before that point is reached? Here we investigate how the spontaneous channelization of a linked-cavity drainage system can control the onset of floods. In agreement with previous work, we show that floods only occur in a band of water throughput rates in which steady reservoir drainage is unstable, and we identify stabilizing mechanisms that allow steady drainage of an ice-dammed reservoir. We also show how stable limit cycle solutions emerge from the instability and identify parameter regimes in which the resulting floods cause flotation of the ice dam. These floods are likely to be initiated by flotation rather than the unstable enlargement of a distributed drainage system.

A model for interaction between conduits and surrounding hydraulically connected distributed drainage based on geomorphological evidence from Keewatin, Canada

Abstract. We identify and map visible traces of subglacial meltwater drainage around the former Keewatin Ice Divide, Canada, from high-resolution Arctic Digital Elevation Model (ArcticDEM) data. We find similarities in the characteristics and spatial locations of landforms traditionally treated separately (i.e. meltwater channels, meltwater tracks and eskers) and propose that creating an integrated map of meltwater routes captures a more holistic picture of the large-scale drainage in this area. We propose the grouping of meltwater channels and meltwater tracks under the term meltwater corridor and suggest that these features in the order of 10s–100s m wide, commonly surrounding eskers and transitioning along flow between different types, represent the interaction between a central conduit (the esker) and surrounding hydraulically connected distributed drainage system (the meltwater corridor). Our proposed model is based on contemporary observations and modelling which suggest that connections between conduits and the surrounding distributed drainage system within the ablation zone occur as a result of overpressurisation of the conduit. The widespread aerial coverage of meltwater corridors (5 %–36 % of the bed) provides constraints on the extent of basal uncoupling induced by basal water pressure fluctuations. Geomorphic work resulting from repeated connection to the surrounding hydraulically connected distributed drainage system suggests that basal sediment can be widely accessed and evacuated by meltwater.

Switch in chemical weathering caused by the mass balance variability in a Himalayan glacierized basin: a case of Chhota Shigri Glacier

ABSTRACTLong-term data (2003–2015) on meltwater chemistry, mass balance and discharge of a benchmark glacier (Chhota Shigri Glacier, India) were studied to determine any association between these variables. To infer the factors governing the alteration of chemical weathering processes in glacierized basins, multi-annual records of the hydrochemical indices (Ca2++Mg2+/Na++K+) and the C-ratio were also examined. A succession of negative mass balance years has resulted in a decline in solute concentrations in the runoff, as discharge has increased. The (Ca2++Mg2+/Na++K+) and C-ratio are highest during periods of negative annual mass balance, when the spatial extent of the channelized drainage system increases. Conversely, these ratios are lowest in positive mass balance years, when the spatial extent of the channelized drainage system decreases, and chemical weathering in the distributed drainage system becomes more dominant. This paper is the first to show the inter-annual linkages between meltwater chemistry, mass balance and discharge for a valley glacier.

Utility of 222Rn as a passive tracer of subglacial distributed system drainage

Water flow beneath the Greenland Ice Sheet (GrIS) has been shown to include slow-inefficient (distributed) and fast-efficient (channelized) drainage systems, in response to meltwater delivery to the bed via both moulins and surface lake drainage. This partitioning between channelized and distributed drainage systems is difficult to quantify yet it plays an important role in bulk meltwater chemistry and glacial velocity, and thus subglacial erosion. Radon-222, which is continuously produced via the decay of 226Ra, accumulates in meltwater that has interacted with rock and sediment. Hence, elevated concentrations of 222Rn should be indicative of meltwater that has flowed through a distributed drainage system network. In the spring and summer of 2011 and 2012, we made hourly 222Rn measurements in the proglacial river of a large outlet glacier of the GrIS (Leverett Glacier, SW Greenland). Radon-222 activities were highest in the early melt season (10–15 dpm L−1), decreasing by a factor of 2–5 (3–5 dpm L−1) following the onset of widespread surface melt. Using a 222Rn mass balance model, we estimate that, on average, greater than 90% of the river 222Rn was sourced from distributed system meltwater. The distributed system 222Rn flux varied on diurnal, weekly, and seasonal time scales with highest fluxes generally occurring on the falling limb of the hydrograph and during expansion of the channelized drainage system. Using laboratory based estimates of distributed system 222Rn, the distributed system water flux generally ranged between 1–5% of the total proglacial river discharge for both seasons. This study provides a promising new method for hydrograph separation in glacial watersheds and for estimating the timing and magnitude of distributed system fluxes expelled at ice sheet margins.

Ice-sheet dynamics through the Quaternary on the mid-Norwegian continental margin inferred from 3D seismic data

Reconstructing the evolution of ice sheets is critical to our understanding of the global environmental system, but most detailed palaeo-glaciological reconstructions have hitherto focused on the very recent history of ice sheets. Here, we present a three-dimensional (3D) reconstruction of the changing nature of ice-sheet derived sedimentary architecture through the Quaternary Ice Age of almost 3 Ma. An extensive geophysical record documents a marine-terminating, calving Fennoscandian Ice Sheet (FIS) margin present periodically on the mid-Norwegian shelf since the beginning of the Quaternary. Spatial and temporal variability of the FIS is illustrated by the gradual development of fast-flowing ice streams and associated intensification of focused glacial erosion and sedimentation since that time. Buried subglacial landforms reveal a complex and dynamic ice sheet, with converging palaeo-ice streams and several flow-switching events that may reflect major changes in topography and basal thermal regime. Lack of major subglacial meltwater channels suggests a largely distributed drainage system beneath the marine-terminating part of the FIS. This palaeo-environmental examination of the FIS provides a useful framework for ice-sheet modelling and shows that fragmentary preservation of buried surfaces and variability of ice-sheet dynamics should be taken into account when reconstructing glacial history from spatially limited datasets.

Sliding velocity fluctuations and subglacial hydrology over the last two decades on Argentière glacier, Mont Blanc area

ABSTRACTThe subglacial observatory beneath the Argentière glacier provides a rare opportunity to study the interactions between glacier sliding velocity and subglacial runoff. The sliding velocity has been monitored in this cavity almost continuously since 1997 and the resulting data indicate a decrease in annual sliding velocities over the last two decades. We found close relationships between annual surface velocity, sliding velocity and ice thickness. These relationships indicate that the ice-flow velocity changes do not depend on subglacial water runoff changes at the annual timescale. The seasonal magnitudes of sliding also show a decrease over the last two decades. At the seasonal timescale, sliding velocity increases before or simultaneously with the large runoff increase in May, indicating a distributed drainage system. Conversely, at the end of the melt season, sliding velocity continues to decrease after the runoff returns to low winter values. The simultaneous increases of runoff and sliding velocity occur mainly before the spring transition. Later, sliding velocity generally appears not to be related to water inputs coming from the surface, except for some large accelerations after midAugust that are always associated with periods of rapidly increasing water inputs to the subglacial drainage system.

A model for subglacial flooding through a preexisting hydrological network during the rapid drainage of supraglacial lakes

AbstractIncreasingly large numbers of supraglacial lakes form and drain every summer on the Greenland Ice Sheet. Some of these lakes drain rapidly within the timescale of a few hours, and the vertical discharge of water during these events may find a preexisting film of water potentially within a distributed drainage system of linked cavities. Here we present a model for subglacial flooding applied specifically to such circumstances. One of many interesting results we find is that water flows in from the far field prior to the arrival of flooding. We systematically evaluate the effect of initial ice/bed opening on the degree of perturbation to the subglacial system. Of particular importance, we find that floods propagate much faster and vertical displacements are much greater for larger openings. For 10 and 1 cm of initial opening, for example, floods travel about 68% and 50% farther than in the fully coupled ice/bed scenario after 2 h of drainage, respectively. For the same choices of initial opening, the elastostatic displacement at the injection point is about 1.39 and 1.26 times that of the fully coupled scenario. Using the framework with a preexisting water film results in avoiding the pressure singularity that is inherent to classical hydrofracture models, thus opening an avenue for integrating the likes of our model within continuum subglacial hydrological models. Furthermore, we foresee that the theory presented can be used to potentially infer subglacial hydrological conditions from surface observables.

Feedbacks between coupled subglacial hydrology and glacier dynamics

On most glaciers and ice sheet outlets the majority of motion is due to basal slip, a combination of basal sliding and bed deformation. The importance of basal water in controlling sliding is well established, with increased sliding generally related to high basal water pressure, but the details of the interactions between the ice and water systems has not received much study when there is coupling between the systems. Here we use coupled subglacial hydrology and ice dynamics models within the Community Ice Sheet Model to investigate feedbacks between the ice and water systems. The dominant feedback we find is negative: sliding over bedrock bumps opens additional cavity space, which lowers water pressure and, in turn, sliding. We also find two small positive feedbacks: basal melt increases through frictional heat during sliding, which raises water pressure, and strain softening of basal ice during localized speedup causes cavities to close more quickly and maintain higher water pressures. Our coupled modeling demonstrates that a sustained input of surface water to a distributed drainage system can lead to a speedup event that decays even in the absence of channelization, due to increased capacity of the system through opening of cavities, which is enhanced through the sliding-opening feedback. We find that the negative feedback resulting from sliding-opening is robust across a wide range of parameter values. However, our modeling also argues that subglacial channelization is required to terminate speedup events over timescales that are commensurate with observations of late summer slowdown on mountain glaciers.

Seasonal changes in ice sheet motion due to melt water lubrication

A numerical model is used to calculate how the motion of an idealized ice-sheet margin is affected by the subglacial drainage of melt water from its surface. The model describes the evolution of the drainage system and its coupling with ice flow through a sliding law that depends on the effective pressure. The results predict ice acceleration during early summer when the inefficient drainage system is temporarily overwhelmed. The growth of a more efficient drainage system leads to a subsequent slowdown of the ice very close to the margin, but high water pressure and ice velocity are maintained through much of the summer further inland. Annual mean ice velocity increases with the total quantity of melt water except close to the margin, where it is almost insensitive to the amount of melting. Short-term variability of melt water input leads to rapid changes in ice velocity that result in a slight increase in the mean velocity relative to a smoother input. Linked-cavity and poroelastic models for the distributed drainage system are compared, and their relative merits discussed. Two different sliding laws are considered, and the need for a holistic description of hydraulically controlled sliding is highlighted.

Reduced glacier sliding caused by persistent drainage from a subglacial lake

Abstract. We present velocity observations of a glacier outlet in Vatnajökull, Iceland, deduced from interferometric SAR (InSAR) data obtained during the ERS1/2 tandem mission in 1995–2000. More than a 50% decrease in glacier velocity was observed subsequent to a large jökulhlaup from the subglacial lake Grímsvötn in 1996. The glacier had not reached its former flow rate in 2000. The jökulhlaup damaged the lake's ice-dam causing persistent drainage from the lake. InSAR based studies of water accumulation within Grímsvötn suggest that a leakage of >3 m3 s−1 prevailed throughout our study period. We suggest that the lake leakage kept open a tunnel at low water pressure underneath the whole length of the glacier. The tunnel flow drained water from its surroundings, hence lowering the water pressure of a distributed drainage system, underneath the upper and centre parts of the glacier, which prior to the jökulhlaup sustained significant basal sliding. This is in accordance with theoretical prediction that tunnel flow in a steady state may cause slow-down in glacier velocity by reducing the subglacial water pressure. The width of the affected areas was ~5 km on the upper part of the glacier and ~8 km on the centre part of the glacier. This indicates that the water pressure reduction propagates laterally from the tunnel over a distance of a few km.

Subglacial drainage system structure and morphology of Brewster Glacier, New Zealand

AbstractA global positioning system and ground penetrating radar surveys is used to produce digital elevation models of the surface and bed of Brewster Glacier. These are used to derive maps of subglacial hydraulic potential and drainage system structure using three different assumptions about the subglacial water pressure (Pw): (i) Pw = ice overburden; (ii) Pw = half ice overburden; (iii) Pw = atmospheric. Additionally, 16 dye‐tracing experiments at 12 locations were performed through a summer melt season. Dye return curve shape, together with calculations of transit velocity, dispersivity and storage, are used to infer the likely morphology of the subglacial drainage system. Taken together, the data indicate that the glacier is underlain by a channelised but hydraulically inefficient drainage system in the early summer in which water pressures are close to ice overburden. By mid‐summer, water pressures are closer to half‐ice overburden and the channelised drainage system is more hydraulically efficient. Surface streams that enter the glacier close to the location of major subglacial drainage pathways are routed quickly to the channels and then to the glacier snout. Streams that enter the glacier further away from the drainage pathways are routed slowly to the channels and then to the snout because they first flow through a distributed drainage system. Copyright © 2008 John Wiley & Sons, Ltd.

Thermal conditions at the bed of the Laurentide ice sheet in Maine during deglaciation: implications for esker formation

The University of Maine Ice Sheet Model was used to study basal conditions during retreat of the Laurentide ice sheet in Maine. Within 150 km of the margin, basal melt rates average ∼5 mm a−1 during retreat. They decline over the next 100 km, so areas of frozen bed develop in northern Maine during retreat. By integrating the melt rate over the drainage area typically subtended by an esker, we obtained a discharge at the margin of ~1.2 m3 s-1. While such a discharge could have moved the material in the Katahdin esker, it was likely too low to build the esker in the time available. Additional water from the glacier surface was required. Temperature gradients in the basal ice increase rapidly with distance from the margin. By conducting upward into the ice all of the additional viscous heat produced by any perturbation that increases the depth of flow in a flat conduit in a distributed drainage system, these gradients inhibit the formation of sharply arched conduits in which an esker can form. This may explain why eskers commonly seem to form near the margin and are typically segmented, with later segments overlapping onto earlier ones.

Peat land oxidation enhances subsidence in the Venice watershed

The southernmost part of the Venice Lagoon catchment was progressively reclaimed from marshland starting from the end of the 19th century and finishing in the late 1930s (Figure 1). As a major result, the area was turned into a fertile farmland. At present, the area is kept dry by a distributed drainage system that collects the water from a capillary network of ditches, and pumps it into the lagoon or the sea. By its very origin, this area lies below sea level and progressively sinks mainly because of bio‐oxidation of the histosols (soils with high organic content) that represent a large fraction of the outcropping soil in the area. The bio‐oxidation process occurs in close connection with the agricultural practices and is currently responsible for a subsidence rate of between 1.5 and 2 cm/yr.

Integrated hydrologic and hydrochemical observations of Hidden Creek Lake jökulhlaups, Kennicott Glacier, Alaska

[1] Hidden Creek Lake (HCL), an ice-marginal lake impounded by Kennicott Glacier, Wrangell Mountains, Alaska, fills annually to ∼20 to 30 × 106 m3 and then drains subglacially within 2 to 3 days. During the 1999 and 2000 jokulhlaups, we carried out a series of planned observations around the lake and in the Kennicott River, which drains the glacier. Approximately 20% of the lake volume was contained within a subglacial water “wedge” beneath the ice dam. The entire volume of the lake drains through the wedge; hydraulic head loss through this constriction may be responsible for the fairly symmetrical shape of the HCL outflow hydrographs, deduced from lake level records, basin hypsometry, and collapse of the ice dam. The flood hydrographs in the Kennicott River are similar in shape to the outflow hydrographs, and within error, lake volume matched the river flood volume in both years. Up to 12 × 106 m3 of water was temporarily stored within the glacier during the 2000 jokulhlaup. During the 2000 jokulhlaup the background flow in the Kennicott River shifted to a dilute chemical composition. As the HCL jokulhlaup progressed, Donoho Falls Lake filled with water whose chemistry was closer to that of the background flow in Kennicott River than to HCL water. Comparison of these chemical signals with typical summer variations in Kennicott River chemistry suggests that the jokulhlaup created high subglacial water pressure that impeded normal drainage of solute-rich water from a distributed drainage system into a conduit system at the glacier bed and even caused flow direction locally to reverse.

Influence of subglacial drainage system evolution on glacier surface motion: Haut Glacier d'Arolla, Switzerland

The relationship between the evolution of subglacial drainage system morphology and spatial patterns of glacier surface velocity was investigated using dye tracing experiments and ground surveying throughout the 1995 melt season at Haut Glacier d'Arolla, Switzerland. With the onset of high and variable melt season discharges, subglacial drainage changed from a predominantly distributed system to a predominantly channelized system. The change occurred later farther up glacier. During the period of drainage evolution the glacier was subjected to three periods of rapidly rising meltwater discharge. The magnitude and spatial pattern of the glacier's velocity response differed between these periods and can be explained in terms of the impact of the evolving drainage system morphology on the amplitude and spatial distribution of basal hydrological forcing. Increasing discharge through a distributed drainage system caused widespread basal forcing and high glacier velocity. Increasing discharge through incipient channels below moulins, not yet connected to the main channel, caused more localized basal forcing and slightly increased glacier velocity. Increasing discharge through a fully channelized drainage system caused no significant basal forcing and glacier velocity was not significantly different from the annual deformation flow pattern. Empirical orthogonal function analysis of flow patterns defined two distinct spatial modes of surface velocity which corresponded closely with the drainage system morphologies inferred to be present during each event. The relative importance of these modes changed through the melt season, suggesting a temporal change in the spatial pattern of hydrologically induced basal forcing.

Geochemical weathering at the bed of Haut Glacier d'Arolla, Switzerland—a new model

AbstractWaters were sampled from 17 boreholes at Haut Glacier d'Arolla during the 1993 and 1994 ablation seasons. Three types of concentrated subglacial water were identified, based on the relative proportions of Ca2+, HCO3− and SO42− to Si. Type A waters are the most solute rich and have the lowest relative proportion of Si. They are believed to form in hydrologically inefficient areas of a distributed drainage system. Most solute is obtained from coupled sulphide oxidation and carbonate dissolution (SO–CD). It is possible that there is a subglacial source of O2, perhaps from gas bubbles released during regelation, because the high SO42− levels found (up to 1200 µeq/L) are greater than could be achieved if sulphides are oxidized by oxygen in saturated water at 0 °C (c.414 µeq/L). A more likely alternative is that sulphide is oxidized by Fe3+ in anoxic environments. If this is the case, exchange reactions involving FeIII and FeII from silicates are possible. These have the potential to generate relatively high concentrations of HCO3− with respect to SO42−. Formation of secondary weathering products, such as clays, may explain the low Si concentrations of Type A waters. Type B waters were the most frequently sampled subglacial water. They are believed to be representative of waters flowing in more efficient parts of a distributed drainage system. Residence time and reaction kinetics help determine the solute composition of these waters. The initial water–rock reactions are carbonate and silicate hydrolysis, and there is exchange of divalent cations from solution for monovalent cations held on surface exchange sites. Hydrolysis is followed by SO–CD. The SO42− concentrations usually are <414 µeq/L, although some range up to 580 µeq/L, which suggests that elements of the distributed drainage system may become anoxic. Type C waters were the most dilute, yet they were very turbid. Their chemical composition is characterized by low SO42− : HCO3− ratios and high pH. Type C waters were usually artefacts of the borehole chemical weathering environment. True Type C waters are believed to flow through sulphide‐poor basal debris, particularly in the channel marginal zone. The composition of bulk runoff was most similar to diluted Type B waters at high discharge, and was similar to a mixture of Type B and C waters at lower discharge. These observations suggest that some supraglacial meltwaters input to the bed are stored temporarily in the channel marginal zone during rising discharge and are released during declining flow.Little of the subglacial chemical weathering we infer is associated with the sequestration of atmospheric CO2. The progression of reactions is from carbonate and silicate hydrolysis, through sulphide oxidation by first oxygen and then FeIII, which drives further carbonate and silicate weathering. A crude estimate of the ratio of carbonate to silicate weathering following hydrolysis is 4 : 1. We speculate that microbial oxidation of organic carbon also may occur. Both sulphide oxidation and microbial oxidation of organic carbon are likely to drive the bed towards suboxic conditions. Hence, we believe that subglacial chemical weathering does not sequester significant quantities of atmospheric CO2 and that one of the key controls on the rate and magnitude of solute acquisition is microbial activity, which catalyses the reduction of FeIII and the oxidation of FeS2. Copyright © 2002 John Wiley & Sons, Ltd.

Character of the subglacial drainage system in the ablation area of Dokriani glacier, India, as revealed by dye-tracer studies

The character of the subglacial drainage system is best determined by the dye-tracer technique. Ten dye-tracer experiments were conducted at three moulins at 930, 1700 and 2300 m, respectively, from the measurement site on the Dokriani glacier, Garhwal Himalaya, India. The injections were made during July, August and September 2000. In July, the subglacial drainage appeared to be via efficient trunk channels, which resulted in dye-concentration curves with velocities between 0.37 and 0.47 m s −1. In August and September, drainage was active mainly in the inefficiently distributed system, as velocity declined and varied between 0.10 and 0.39 m s −1. Probably, the distributed system during August expanded to accommodate the increased runoff by the seasonal high discharges. Limited observations indicate that there was temporal switching between channelised and distributed drainage system.