Hostname: page-component-8448b6f56d-c4f8m Total loading time: 0 Render date: 2024-04-25T05:13:10.781Z Has data issue: false hasContentIssue false

A conceptual model of solute acquisition by Alpine glacial meltwaters

Published online by Cambridge University Press:  20 January 2017

Martyn Tranter
Affiliation:
Department of Geography, University of Bristol, Bristol BS8 1SS, England
Giles Brown
Affiliation:
Department of Geoography, University of Cambridge, Cambridge CB2 3EN, England
Robert Raiswell
Affiliation:
Department of Earth Sciences, University of Leeds, Leeds LS2 9JT, England
Martin Sharp
Affiliation:
Department of Geography, University of Cambridge, Cambridge CB2 3EN, England
Angela Gurnell
Affiliation:
Department of Geography, University of Southampton, Southampton SO9 5NH, England
Rights & Permissions [Opens in a new window]

Abstract

Solute acquisition by Alpine glacial meltwaters is the result of the coupling of different pairs of reactions, one of which usually involves dissolved gases. Hence, the availability of atmospheric gases to solution is an important control on the composition of glacial meltwaters. The chemical compositions of the two main components of the bulk meltwater, quick flow and delayed flow, are dominated by different geochemical processes. Delayed flow waters are solute-rich and exhibit high p(CO2) characteristics. The slow transit of these waters through a distributed drainage system and the predominance of relatively rapid reactions, such as sulphide oxidation and carbonate dissolution, in this environment maximize solute acquisition. Quick-flow waters are dilute, both because of their rapid transit through ice-walled conduits and open channels, and because the weathering reactions are fuelled by relatively slow gaseous diffusion of (CO2) into solution, despite solute acquisition being dominated by rapid surface exchange reactions. As a consequence, quick flow usually bears a low or open-system p(CO2) signature. Bulk meltwaters are more likely to exhibit low p(CO2) values when suspended-sediment concentrations are high, which promotes post-mixing reactions. This conceptual model suggests that the composition of both quick flow and delayed flow is likely to be temporally variable, since kinetic, rather than equilibrium, factors determine the composition.

Type
Research Article
Copyright
Copyright © International Glaciological Society 1993

Introduction

The chemistry of (bulk) meltwaters draining from Alpine glaciers may enable seasonal variations in water routing through the hydrological system to be inferred, provided that realistic estimates can be made of the compositions of waters which have followed different flowpaths through the glacier drainage system (Reference CollinsCollins, 1979a, Reference Collinsb, Reference Collina1981; Reference SharpSharp, 1991). Seasonal variations in water routing and in subglacial water pressure are thought to influence both glacier motion (Reference IkenIken, 1981; Reference BindschadlerBindschadler, 1983; Reference Iken, Röthlisberger, Flotron and HaeberliIken and others, 1983; Lliboutry, 1983; Reference Weertman and BirchfieldWeertman and Birchfield, 1983) and the tapping of suspended-sediment sources from subglacial environments (Reference Gurnell and FennGurnell and Fenn, 1984; Reference Gurnell, Gurnell and ClarkGurnell, 1987; Reference CollinsCollins, 1989). Earlier consensus regarding the chemical composition of the two or more components which constitute the bulk meltwaters (Reference CollinsCollins, 1978; Reference Oerter, Behrens, Hibsch, Rauert and StichlerOerter and others, 1980; Reference Gurnell and FennGurnell and Fenn, 1984) is that they may be assumed to be relatively constant throughout the ablation season (Reference CollinsCollins, 1978, Reference Collins1979b). However, recent work challenges this assumption (Reference Tranter, Davies, Tranter and JonesTranter and Raiswell, 1991; paper in preparation by M. Sharp and others). Hence, it seems timely to explain how chemical-weathering reactions in dilute solutions give rise to meltwater chemistry which is more variable than hitherto supposed.

The notion that the composition of each component may be constant arises because the existing glaciological literature tends to focus on individual, rapid weathering reactions, such as so-called “ion exchange” (Reference Lemmens and RogerLemmens and Roger, 1978), in isolation from the reactions which furbish protons, or H+ (aq). A supply of protons is usually required for the chemical weathering of glacial flour and may be generated by the relatively slow diffusion of gaseous (CO2) into solution (Reference RaiswellRaiswell, 1984). We show below that it is often the coupling of such reactions, rather than the exchange reaction alone, which results in solute acquisition by glacial meltwaters.

Our purpose in this paper is to outline a conceptual model of solute acquisition by Alpine glacial meltwaters, which identifies the likely sources of the solute, details the key chemical-weathering reactions and outlines the likely linkages between these reactions. Finally, we identify the relative rates of the reactions and show how meltwater chemistry bears the fingerprint of the differing solute-acquisition rates. We illustrate the text with examples of Alpine meltwater chemistry, mainly from Haut Glacier d’Arolla, Switzerland. Details of sample collection and treatment may be found in Reference BrownBrown (1991) and Reference Brown and TranterBrown and Tranter (1990).

The Components of Bulk Meltwater – Quick Flow and Delayed Flow – and Subglacial Drainage Routes

Bulk meltwaters are believed to be largely composed of two components, termed the “subglacial” and “en-glacial” components by (Reference CollinsCollins 1978, Reference Collins1979b). Both components are ultimately derived from surface-snow and ice melt. Less significant amounts of water can be derived from rainfall, ground water, melting due to viscous dissipation in deforming ice, and geothermal melting (Reference Collins and YoungCollins and Young, 1981). A two-component mixing model clearly simplifies the nature of water routing through Alpine glaciers. However, there is strong evidence to suggest that this type of model does indeed describe many features of the subglacial drainage system (paper in preparation by M. Sharp and others).

We prefer to use the terminology, “quick flow” and “delayed flow”, to describe the two major components. This better identifies the different residence times of the components in the subglacial drainage sytem, which are the result of different hydraulic conditions in the different parts of the system. Further, it allows both components to follow separate flowpaths prior to mixing and transit through the major conduits which terminate at the portal (paper in preparation by M. Sharp and others).

Quick flow pertains to waters in rapid transit through the hydroglacial system, predominantly through ice-walled conduits (Reference RöthlisbergerRöthlisberger, 1972) and basal channels. Delayed flow pertains to waters in slow transit through a distributed drainage system at the glacier bed (paper in preparation by M. Sharp and others), such as a linked-cavity system (Reference WalderWalder, 1986). Unlike temperate catchments, where “old” water or delayed flow may dominate peak run-off (Reference TurnerTurner and others, 1990), dye-tracing studies at Haut Glacier d’Arolla demonstrate that quick flow makes a significant contribution to the diurnal discharge maximum (paper in preparation by P. W. Nienow and others). These studies also demonstrate the existence of drainage systems with characteristically short (< 2h) and long (> 12h) residence times. Hence quick flow and delayed flow are pertinent descriptions of the two components.

Alpine glaciers may have an entirely distributed drainage sytem or a combined distributed-channel system, but they are unlikely to have only a channel system (Reference FowlerFowler, 1987). In the case of the combined system, the distributed system drains into the channel system, particularly during the height of the ablation season, and there is great potential for post-mixing chemical evolution of the bulk meltwaters via interactions with elevated suspended-sediment concentrations (Reference SharpSharp, 1991). The seasonal growth of the channel system and the implications for suspended-sediment concentrations in bulk meltwaters have been reviewed by Reference SharpSharp (1991), and an important observation is that suspended-sediment concentrations appear to be low when the quick-flow component is low. In general, delayed flow dominates bulk meltwater discharge early in the ablation season, when suspended-sediment concentrations are low, and the importance of quick flow increases as the ablation season progresses.

Quick flow and delayed flow are both derived mainly from surface-snow and ice melt, and the original solute content of these waters will be generally low (see below). Hence, the different composition of the components is derived from the differing types and the duration of chemical-weathering reactions that the waters undergo during transit through the hydroglacial system.

A Conceptual Model for the Acquisition of Solute by Alpine Glacial Meltwaters

Solute is defined as all ions, including base cations (e.g. Ca2 +, Mg2 +, K +, Na+), aqueous protons (H+ (aq)) and anions (HCO3 -, SO4 2-, NO3 -, Cl-), and neutral species (e.g. O2, N2, CO2, H4SiO4) dissolved in water (Reference Stumm and MorganStumm and Morgan, 1981). Solute acquisition (or chemical weathering) is therefore concerned with reactions which affect species in all phases: solid, liquid and gas. To date, the glaciological literature has given scant attention to the importance of dissolved gases and their impact on the rate and extent of chemical weathering.

Solute may be dervied from atmospheric deposition, either contained in snow or rain, or from dry deposition on to the catchment surface (Reference Barrie, Davies, Tranter and JonesBarrie, 1991; Reference Cadle, Davies, Tranter and JonesCadle, 1991). The first-derived snowmelt or rain-snowmelt mixture is more concentrated than that derived from leached snow cover, since solute from snow readily fractionates into snowmelt and rainfall (Reference Johannessen and HenriksenJohannessen and Henriksen, 1978; Reference Jones and DaviesJones and others, 1989; Reference Tranter, Davies, Tranter and JonesTranter, 1991). Thereafter, snow-and ice-melt is relatively pure, as confirmed by studies which have sampled supraglacial meltwaters in ice-walled channels near the glacier snout (Reference CollinsCollins, 1978; Reference Raiswell and ThomasRaiswell and Thomas, 1984; Reference Tranter and RaiswellTranter and Raiswell, 1991). The snow line retreats up-glacier as the ablation season proceeeds, and therefore snowpack leaching occurs later at higher altitudes. However, much of the solute found in the bulk meltwaters emerging at the portal is acquired by chemical weathering in subglacial environments (Reference Souchez, Lemmens, Gurnell and ClarkeSouchez and Lemmens, 1987). Exceptions to this rule are species such as NO3 and O2 (aq), which have a predominantly atmospheric source.

The potential for meltwaters to acquire solute in subglacial locales is greatest if there is slow transit of waters through an environment which promotes relatively rapid rates of chemical weathering. This type of subglacial environment is most likely to be found in distributed drainage sytems, such as the linked cavity system, where (i) there is a copious supply of freshly comminuted rock flour, (ii) the rock flour contains reactive minerals, such as carbonates and sulphides, (iii) the residence time of the meltwater is relatively long, and (iv) relatively rapid sulphide oxidation provides a proton supply. The ubiquitous minerals, calcite (a carbonate) and pyrite (a sulphide), are important sources of solute to meltwaters, despite being usually a minor component of the bedrock, since their dissolution kinetics are orders of magnitude more rapid than those of most other rock-forming minerals (Reference TranterTranter, 1982).

By contrast, chemical weathering in open channels is more likely to be driven by slower gaseous diffusion of CO2 into solution, since the trace amounts of sulphide and carbonate present in the glacial flour will have been largely dissolved in the distributed system. As will become evident, chemical weathering in full conduits may be limited because of the restricted access to a source of gaseous CO2.

Sources of Base Cations

Meltwaters acquire base cations via two major types of reaction, simple dissolution and acid hydrolysis. Simple dissolution reactions do not involve protons. Rather, they involve solids such as sea salt, derived from dry deposition on to the catchment surface or aerosol scavenged by snowfall, or evaporitic minerals in the bedrock, such as halite and gypsum (Reference Raiswell, Brimblecombe, Dent and LissRaiswell and others, 1980). Usually, these reactions supply minor quantities of base cations to Alpine glacial meltwaters because of the continental location of the Alps and the lithology of the bedrock. However, if sources of evaporitic minerals are found in glacial catchments, they are likely to have a marked impact on the composition of meltwater.

Instead, most base cations usually result from the broad class of reactions known as acid hydrolysis (Reference Raiswell, Brimblecombe, Dent and LissRaiswell and others, 1980). These reactions require a source of aqueous protons, which are exchanged for base cations from the rock-flour surface (see Equations (1) and (2)).

(1)
. (2)

The aqueous protons required for acid hydrolysis are derived from two main sources. The first is the dissociation of dissolved CO2, derived via diffusion, either directly from the atmosphere or indirectly from gas bubbles within melting ice:

(3)

The second source of protons is the oxidation of sulphides, such as pyrite (FeS2), which are commonly found in the bedrock of Alpine glacier basins in trace amounts:

. (4)

Oxidation of sulphides normally requires the presence of O2, but under certain circumstances, e.g. in water-filled subglacial channels where anaerobic conditions may develop, other oxidizing agents, such as NO3 , may be used (paper in preparation by M. Sharp and others). More detailed accounts of the mechanism and rate of pyrite oxidation may be found in Reference Singer and StummSinger and Stumm (1970) and Reference Toran and HarrisToran and Harris (1989).

Protons may also be derived from the dissolution of acids, either scavenged from the atmosphere by snow and rain or dry-deposited on to the glacier surface (Reference Barrie, Davies, Tranter and JonesBarrie, 1991; Reference Cadle, Davies, Tranter and JonesCadle, 1991). This is particularly the case early in the melt season, when snowmelt is usually most acidic (Reference Tranter, Davies, Tranter and JonesTranter, 1991).

(5)
. (6)

The major base cation in glacial meltwaters is Ca2 +(Reference RaiswellRaiswell, 1984), which partially reflects the relatively rapid rate of dissolution of Ca2 + from silicate and aluminosilicate minerals relative to other base cations (Reference LermanLerman, 1979), and the fact that carbonate dissolution is a major source of dissolved ions to delayed flow (see below).

The Source of Anions

The common acid anions, Cl, HCO3 , NO3 and SO4 2-, can be derived from a number of sources, as illustrated above, but it is normally safe to assume that Cl and NO3 are derived mainly from the dissolution of sea-salt aerosol (Reference HollandHolland, 1978; Soudiez and Lemmens, 1987) and acidic nitrate aerosol (see Equations (1) and (6)), respectively. The latter is true because significant numbers of N2-fixing bacteria are unlikely within the glacial system. By contrast, most SO4 2- is likely to be derived from the oxidation of sulphides in comminuted bedrock (Reference Tranter and RaiswellTranter and Raiswell, 1991). Smaller amounts are derived from the dissolution of atmospheric aerosol (Reference BrownBrown, 1991), as represented in Equation (5). Bicarbonate may be derived from two main sources, from the dissociation of atmospheric CO2 (see Equation (3)) and the dissolution of carbonates, such as calcite (see Equation (2)).

The two major anions in Alpine glacial meltwaters are HCO3 and SO4 2- (Reference Souchez, Lemmens, Gurnell and ClarkeSouchez and Lemmens, 1987), because the dissociation of atmospheric CO2 and the oxidation of sulphides provide the bulk of the protons used to weather chemically carbonate, silicate and aluminosilicate minerals in glacial flour. Hence, the supply of atmospheric CO2 and O2 to meltwaters is a fundamental control on the rate and the extent of chemical weathering in glacial meltwaters.

Surface Exchange and Carbonation – Dominant Reactions in Quick Flow and The Post-Mlxing Chemical Evolution of Bulk Meltwaters

Existing studies of chemical weathering in glacial environments have focused on so-called ion-exchange processes (Reference Lorrain and SouchezLorrain and Souchez, 1972; Reference Lemmens and RogerLemmens and Roger, 1978; Reference Souchez, Lemmens, Gurnell and ClarkeSouchez and Lemmens, 1987), which are relatively rapid reactions. We prefer to adopt the following nomenclature, since it gives a better descrip-tion of the reaction mechanism. Ion exchange is the rapid, reversible exchange of ions between solution and a solid surface (Reference Stumm and MorganStumm and Morgan, 1981). This type of exchange is said to be non-specific, because no chemical bonds are formed between the adsorbed ion and the mineral surface (Reference ParksParks, 1990). Cation exchange is the rapid, reversible exchange of positively charged ions. Ion and cation exchange, by definition, cannot add solute to meltwaters (cf. Reference Lemmens and RogerLemmens and Roger, 1978; Reference DrewryDrewry, 1986). These reactions simply modify the relative proportion of ions present in a fixed pool of dissolved ions. The initial reaction between water and freshly comminuted aluminosilicates or silicates usually involves a surface exchange of base cations for hydrogen ions (or protons) from solution (Reference WollastWollast, 1967; Reference Luce, Barlett and ParksLuce and others, 1972; Reference Busenberg and ClemencyBusenberg and Clemency, 1976; Reference LermanLerman, 1979). Surface, rather than ion, exchange is the term used for this reaction because although the reaction is rapid (typically occurring within 3 min; Reference Busenberg and ClemencyBusenberg and Clemency, 1976), it is largely irreversible (Reference Garrels and HowardGarrels and Howard, 1956). The surface of the mineral may be physically altered during surface exchange (Reference Casey and BunkerCasey and Bunker, 1990) rendering it difficult to replace (or exchange) the adsorbed protons for dissolved base cations (Reference Garrels and HowardGarrels and Howard, 1956). This type of exchange is said to be specific, since a chemical bond is formed between the adsorbed ion and the mineral surface (Reference ParksParks, 1990). Full and detailed accounts of adsorption phenomena can be found in Reference Hochella and WhiteHochella and White (1990).

Surface exchange has a marked effect on the ability of the solution to dissolve C02. The loss of protons from solution increases the pH and, as a consequence, markedly increases the quantity of dissolved CO2, HCO3 and CO3 2 that the solution can hold for a fixed partial pressure of carbon dioxide, pCO2 (see Fig. 1a). When isolated from the gas phase, the amount of dissolved CO2 in the solution, CO2 (aq), falls following surface exchange, since it is converted into HCO3 (see Equation (3)), as shown in Figure 1b. Hence, CO2 from the gas phase, which may be the atmosphere or bubbles in the ice walls for example diffuses into solution until the pCO2 in the gas phase, and CO2 (aq) are in equilibrium (see Appendix). As CO2 diffuses into the solution, a large proportion dissociates into HCO3 and H+ (aq). These protons are usually rapidly removed from solution by surface exchange and the slower dissolution of the mineral lattice, and are largely replaced by base cations from the rock-flour surface. Alumina and silica are also released into solution by the dissolution of the mineral lattice (Reference Chou and WollastChou and Wollast, 1984; Reference LasagaLasaga, 1984; Reference Casey and BunkerCasey and Bunker, 1990). This type of “carbonation” reaction (Reference Reynolds and JohnsonReynolds and Johnson, 1972) is most likely to occur in open channels or partially filled cavitities, where there is ample access to rock flour and gaseous C02 (cf. full conduits and water-filled cavitites, where access to gaseous CO2 is limited). Hence, carbonation will largely be responsible for solute acquisition by quick flow, because there are no other significant sources of protons and sulphide oxidation is largely confined to the distributed system (see below).

Fig. 1. a. The relationship between the sum of dissolved carbon dioxide species (CO2(aq) + HCO3 + CO3 2-) and pH. The ability of the solution to dissolve CO2 increases with increasing pH because CO2 (aq) dissociates into HCO3 and CO3 2-. The diagram was constructed by assuming that the solution is in equilibrium with a gas phase p(CO2) of 10−3.5 atm, that the Henry’s law constant, KH, is 10−3.5 mol atm 1−l and that the first and second dissociation constants for dissolved CO2 are 10−6.58 and 10−10.62 mol 1−1, respectively. These constants are for T= 0° C (Reference Garrels and ChristGarrels and Christ, 1965). b. The effect of pH on the speciation of dissolved CO2. The sum of the dissolved species remains constant at 24 μ,eq 1−1 i.e. no CO2 is allowed to diffuse into or out of solution, c. As the pH increases, CO2 (aq) is converted into HCO3 and the p(CO2) in the gas phase with which the solution is in equilibrium decreases (see Appendix). Equilibrium constants as in Figure 1a.

Carbonation will also be largely responsible for the solute acquisition which occurs after the quick flow and delayed flow mix in the basal channel system (Reference Raiswell and ThomasRaiswell and Thomas, 1984; Reference Tranter, Raiswell, Mills and MilesTranter and others, 1989) as confirmed by the non-conservative mixing properties of base cations and HCO3 relative to SO4 2- (paper in preparation by M. Sharp and others). Carbonation reactions therefore dominate the composition of bulk meltwaters when the proportion of quick flow is high and when post-mixing reactions add significant quantities of solute to solution.

Solute Acquisition as a Result of Sulphide Oxidation – A Dominant Control on The Composition of Delayed Flow

Protons may also be derived from the oxidation of suphides, such as pyrite (see Equation 6). Sulphide oxidation is likely to liberate protons and SO4 2- (i.e. sulphuric acid) at or near the site of rock-flour generation, when the rock flour makes contact with waters containing suitable oxidizing agents (e.g. O2). Freshly comminuted rock flour is most likely to contain unaltered sulphide minerals, since sulphide oxidation is relatively rapid and goes to completion in the presence of suitable oxidizing agents (Reference Moses, Nordstrom, Herman and MillsMoses and others, 1987). So, chemical weathering in the distributed system provides an additional proton source, which is largely limited by the oxygen content of the meltwater and the trace-sulphide content of the rock flour.

The protons liberated by sulphide oxidation are likely to be rapidly utilized by surface-exchange reactions and/ or carbonate dissolution at or near the site of rock-flour generation. Hence, base cations, rather than hydrogen ions, are the counter ions to SO4 2- in delayed flow (Reference Tranter and RaiswellTranter and Raiswell, 1991). Whereas the product of sulphide oxidation, SO4 2-, can be leached from the comminuted rock flour by the delayed flow, the rock flour itself only seems to leave the distributed system in relatively high concentrations during active periods of channel growth (Reference CollinsCollins, 1989; Reference SharpSharp, 1991). Hence, we believe that reactive sulphides, either present as comminuted particles or close to the surface of the rock flour, will be exhausted within the distributed system, and are largely absent from rock flour transported from the distributed system into the channel system.

The assumption that SO4 2- is largely confined to the delayed flow (Reference Tranter and RaiswellTranter and Raiswell, 1991) is likely to be significantly in error only during the early ablation season, when relatively concentrated snowmelt may directly run off the glacier surface and mix with bulk meltwaters, consisting predominantly of delayed flow. Even if the suspended sediment contains weatherable sulphide, it is likely that this will be relatively slow to react, since the reactive sulphide would have been consumed within the distribution system and the remaining sulphide would have to be relatively inert to survive transport through the distributed system.

Open-and Closed-System Weathering Regimes – Fingerprints of The Different Chemical-Weathering Reactions in The Distributed and Conduit/Basal Channel Systems

A solution is said to be open when the p(CO2) in the gas phase is in equilibrium with the amount of dissolved CO2 in the liquid phase (i.e. p(CO2) ≈ K H (CO2) (aq) – see the Appendix for the significance of KH). The system is said to be closed when this is not the case (i.e. p(CO2) ≠ K H (CO2) (aq)). Open-and closed-system weathering regimes are kinetic phenomena and depend on the rate at which CO2 diffuses into or out of solution relative to the rate of other chemical-weathering reactions (Reference RaiswellRaiswell, 1984; Reference Raiswell and ThomasRaiswell and Thomas, 1984; Reference Thomas and RaiswellThomas and Raiswell, 1984).

If the rate of weathering is greater than the rate of proton supply, the resultant increase in pH induces gaseous CO2 diffusion into solution because the solution can hold additional (CO2) (aq). We define this solution as a “low-p(CO2) system”. Conversely, a high-p(CO2) system is one in which protons added to an open system are not rapidly consumed by weathering reactions, and CO2 (aq) attempts to diffuse out of solution since some HCO3 is converted to CO2 (aq) by the additional protons.

Low-p(CO2) systems are common when a large supply of reactive, comminuted rock is brought into contact with relatively pure water (see Fig. 2). The glacial analogue is when large quantities of dilute quick-flow waters mix with sediment-laden delayed flow at the height of the ablation season (Reference SharpSharp, 1991; paper in preparation by P. W. Nienow and others; paper in preparation by M. Sharp and others). Given these circumstances, the rate of proton consumption is greater than the rate of gaseous diffusion of CO2 into solution, and (CO2) (aq) remains below open-system values (see Fig. 3).

Fig. 2. The p(CO2) of a mixture of Na-feldspar and deionized water, as a function of time and rock: water ratio. These results were obtained by adding a known weight of crushed feldspar to 11 of deionized water in open contact with the atmosphere at Τ = 20° C (Reference TranterTranter, 1982). Mote that most values are below atmospheric p(CO2) (= l(T−3.5 atm). When p(CO2) is declining, the rate of weathering is greater than the rate of CO2 influx into the solution. The converse argument holds when p(CO2) is rising.

Fig. 3. The (a) discharge, (b) suspended-sediment concentration, (c) p(CO2) and (d) sulphate concentration of meltwaters draining Haut Glacier d’Arolla during the 1989 ablation season. There is some evidence for a decrease in p(CO2) as the discharge and suspended-sediment concentrations increase (Reference BrownBrown, 1991).

High-p(CO2) systems may arise from three types of scenario. First, an input of protons to the drainage system may occur, so as to give a short-lived reduction in pH, for example when acidic snowmelt mixes with a subglacial water body. We believe that this scenario is transient in most subglacial environments, since rock flour has a great capacity for adsorbing protons. Secondly, the freezing of a water body may result in an increase in (CO2) (aq), since CO2 is not readily incorporated into the ice lattice. However, it seems unlikely that this mechanism of high p(CO2) generadon operates on any other than a local scale in most Alpine glaciers. Finally, neutralization of acidity by carbonates may produce high (CO2) (aq), particularly in water-filled cavities or channels where there is only limited opportunity for the solution to degas.

Waters of the delayed flow contain both high concentrations of HCO3 and SO2 2- (see Table 1). The following argument suggests that the delayed flow appears to exhibit high p(CO2). Figure 4 demonstrates that the mass fraction of SO4 2 (defined as the ratio of SO4 2 to the sum of SO4 2 and HCO3 , where concentrations are in units of equivalents) in the bulk meltwaters draining Haut Glacier d’Arolla decreases as discharge increases. This suggests that delayed flow, which dominates at lower discharge, contains a higher mass fraction of SO4 2 than quick flow. Figure 5 suggests that the pCO2 of bulk meltwater increases as the mass fraction of SO4 2- increases. Since bulk meltwaters dominated by delayed flow exhibit high p(CO2), delayed flow is also likely to exhibit high p(CO2).

Fig. 4. The mass fraction of and where units are equivalents), as a function of discharge for bulk meltwaters draining Haut Glacier d’Arolla throughout the 1989 ablation season.

Fig. 5. The association of p(CO2) with the mass fraction of SO2 for bulk meltwaters of Haut Glacier d’Arolla throughout the 1989 ablation season.

Table 1. Estimates of the composition of delayed flow

The combination of high concentrations of Ca2+ – HCO3 – S04 2 in the delayed flow (see Table 1) in conjunction with elevated p(CO2) values tends to suggest that the coupling of sulphide oxidation and carbonate dissolution are dominant controls on the composition of delayed flow, as illustrated by Equation (7). The high p(CO2) is generated by the production of HCO3 on the righthand side of the equation.

(7)

The resulting high p(CO2) signature so generated may be diminished or lost as a consequence of post-mixing chemical-weathering reactions and degassing of bulk meltwaters in open channels (Reference RaiswellRaiswell, 1984). This helps to explain some of the scatter displayed in Figure 5.

Low p(CO2) Signatures in Bulk Melt-Waters – A Consequence of Post-Mixing Reactions

Figure 4 demonstrates that bulk meltwaters at high discharge have a low mass fraction of SO4 2-. Figure 5 shows that these waters have low-p(CO2) characteristics. It is likely that the low-p(CO2) signature is attributable to post-mixing weathering reactions, of the type illustrated in Figure 2 and described by Equations (1) and (2). Quick flow passing through ice-walled conduits has little opportunity to gain access to rock flour and hence will maintain an open-system p(CO2) or a slightly lower p(CO2), depending on the amount of rock-water contact that is possible. By contrast, bulk meltwaters in basal channels clearly gain access to and transport large quantities of rock flour, which is capable of undergoing chemical weathering (Reference Tranter, Raiswell, Mills and MilesTranter and others, 1989). Therefore, it is most likely that post-mixing reactions promote the low-p(CO2) signatures of bulk meltwaters at high discharge.

Are Delayed Flow and Quick-Flow Compositions Constant Throughout The Ablation Season?

The above conceptual model of solute acquisition by Alpine glacial meltwaters stresses that dissolution and gas-diffusion kinetics are important controls on the composition of the two major discharge components and their post-mixing chemical evolution. These reaction kinetics are dependent on the residence time of waters within the channel and distributed systems, the supended-sediment concentration and the access of an atmospheric gas supply, factors which clearly change during the ablation season as the glacial drainage system evolves (paper in preparation by M. Sharp and others). Hence, it is most probable that the composition of each component will likewise evolve during the ablation season. A similar line of reasoning also suggests that the extent of post-mixing chemical evolution of the meltwaters will vary throughout the season.

Conclusions

Solute acquisition in glacial meltwaters results from reactions which involve species in all phases: solid, aqueous and gaseous. Typically, chemical weathering involves two coupled reactions, such that dissolved gases affect the nature and rate of solute acquisition. Surface-exchange reactions in open channels and conduits promote net solute acquisition as a result of influx of gaseous CO2 – so-called carbonation reactions. Oxidation of sulphides and the dissolution of carbonate is a characteristic coupled reaction of the distributed system.

Waters of the delayed flow are relatively concentrated because they flow slowly through the distributed system, where they are in contact with large quantities of relatively reactive, recendy comminuted rock flour, and importantly, because there is a renewable supply of protons available to sustain prolonged chemical weathering. These protons may be derived from the input of acidic snowmelt early in the ablation season, but oxidation of sulphides is the more sustainable source throughout the season. The coupling of sulphide oxidation and carbonate dissolution may impart a high-p(CO2) signature to the delayed flow, as well as relatively high concentrations of Ca2 +, HCO3 and SO4 2-.

Quick flow waters are dilute because of their short residence time in the conduit system and because the supply of protons is limited to slow gaseous diffusion of CO2 into solution. Quick flow waters may have open-system or low-p(CO2) signatures, depending on the extent of chemical weathering they undergo and the degree to which the channelized systems is in contact with a gaseous phase.

Mixing of the components in the channelized system may allow post-mixing chemical weathering of rock flour, fuelled by protons supplied from influxing CO2 and acidic snowmelt (paper in preparation by M. Sharp and others). Ultimately, the p(CO2) exhibited by bulk meltwaters is a balance between the rate at which CO2 can diffuse in or out of solution and the rate of proton consumption by chemical weathering. Hence, when suspended-sediment concentrations are high, the p(CO2) will be low and, when suspended-sediment concentrations are low, bulk meltwaters take on higher p(CO2) values characteristic of the delayed flow, or approach open-system values if degassing of CO2 can occur.

Given these kinetic controls on the composition of each component, it is unlikely that the components maintain a constant composition over the ablation season.

Acknowledgements

This paper arises from work undertaken in association with the U.K. Natural Environment Research Council grant No. GR3.7004A and U.K. NERC Studentship No. GT4/88/AAB/56.

The accuracy of references in the text and in this list is the responsibility of the authors, to whom queries should be addressed.

Appendix

Equilibrium between CO2 in the gas phase and the dissolved phase is dependent on the Henry’s Law constant, K H as described by Equation (Al).

(A1)
(A2)

It is important to remember that this situation pertains to (GO2) (aq) being in equilibrium with p(CO2). Now imagine the situation in which CO2 has been added to the solution, so that the new quantity, (CO2) (aq), is in excess of (CO2)* (aq). It follows from Equation (A2) that the apparent gas-phase partial pressure of CO2 with which the solution is in equilibrium, p(CO2), is likewise greater that p(CO2). So, the solution degasses until a new equilibrium between gas phase and dissolved phase CO2 is reached.

The terminology, “high-p(CO2) system”, derives from the fact that the solution is in apparent equilibrium with a higher p(CO2) than actually exists. A similar line of reasoning gives rise to the converse case of the low-p(CO2) system.

References

Barrie, L. A. 1991. Snow formation and processes in the atmosphere that influence its chemical composition. In Davies, T. D., Tranter, M. and Jones, H.G. eds. Seasonal snowpacks; processes of compositional change. Berlin, Springer-Verlag, 120.Google Scholar
Bindschadler, R. 1983. The importance of pressurized subglacial water in separation and sliding at the glacier bed. J. Glacial., 29(101), 319.Google Scholar
Brown, G. H. 1991. Solute provenance and transport pathways in Alpine glaciers. (Ph.D. thesis, University of Southampton.)Google Scholar
Brown, G.H. and Tranter, M. 1990. Hydrograph and chemograph separation of bulk meltwaters draining the Upper Arolla glacier, Valais, Switzerland. International Association of Hydrological Sciences Publication 193 (Symposium at Lausanne 1990 – Hydrology in Mountainous Regions. I. Hydrological Measurements; the Water Cycle), 429437.Google Scholar
Busenberg, Ε. and Clemency, C.V. 1976. The dissolution kinetics of feldspars at 25°C and 1 alms CO2 partial pressure. Geochim. Cosmochim. Acta., 40, 4149.Google Scholar
Cadle, S. H. 1991. Dry deposition to snowpacks. In Davies, T. D., Tranter, M. and Jones, H.G. eds. Seasonal snowpacks; processes of compositional change. Berlin, Springer-Verlag, 2166.Google Scholar
Casey, W. H. and Bunker, B. 1990. Leaching of mineral and glass surfaces during dissolution. Reviews in Mineralogy., 23, 397426.Google Scholar
Chou, L. and Wollast, R. 1984. Study of weathering of albite at room temperature and pressure with a fluidized bed reactor. Geochim. Cosmochim. Acta., 48, 22052217.Google Scholar
Collins, D.N. 1978. Hydrology of an Alpine glacier as indicated by the chemical composition of meltwater. Z. Gletscherkd. Glazialgeol., 13(1/2), 1977, 219238.Google Scholar
Collins, D.N. 1979a. Hydrochemistry of meltwaters draining from an Alpine glacier. Arct. Alp. Res., 11(3), 307324.CrossRefGoogle Scholar
Collins, D.N. 1979b. Quantitative determination of the subglacial hydrology of two Alpine glaciers. J. Glacioi., 23(89), 347362.CrossRefGoogle Scholar
Collina, D.N. 1981. Seasonal variation of solute concentration in melt waters draining from an Alpine glacier. Ann. Glacial., 2, 1116.Google Scholar
Collins, D.N. 1989. Seasonal development of subglacial drainage and suspended sediment delivery to melt waters beneath an Alpine glacier. Ann. Glacioi., 13, 4550.Google Scholar
Collins, D.N. and Young, G.J. 1981. Meltwater hydrology and hydrochemistry in snow-and ice-covered mountain catchments. Nord. Hydrol., 12(4–5), 319334.Google Scholar
Drewry, D. 1986. Glacial geologic processes . London, Edward Arnold.Google Scholar
Fowler, A.C. 1987. Sliding with cavity formation. J. Glacioi., 33(115), 255267.CrossRefGoogle Scholar
Garrels, R. M. and Christ, C. L. 1965. Solutions, minerals and equilibria. San Francisco, Freeman Cooper.Google Scholar
Garrels, R. M. and Howard, P. 1956. Reactions of feldspar and mica with water at low temperature and pressure. Clays Clay Miner., 6, 6888.CrossRefGoogle Scholar
Gurnell, A. M. 1987. Suspended sediment. In Gurnell, A. M. and Clark, M.J. eds. Glacio-fluvial sediment transfer. New York, John Wiley, 305354.Google Scholar
Gurnell, A. M. and Fenn, CR. 1984. Flow separation, sediment source areas and suspended sediment transport in a pro-glacial stream. Catena Supplement 5, 109119.Google Scholar
Hochella, M. F. and Α. White, F. eds. 1990. Mineral–water interface geochemistry. Reviews in Mineralogy., 23.Google Scholar
Holland, H. D. 1978. The chemistry of atmospheres and oceans . New York, Wiley-Interscience.Google Scholar
Iken, A. 1981. The effect of the subglacial water pressure on the sliding velocity of a glacier in an idealized numerical model. J. Glacioi., 27(97), 407421.Google Scholar
Iken, A., Röthlisberger, H. Flotron, A. and Haeberli, W. 1983. The uplift of Unteraargletscher at the beginning of the melt season – a consequence of water storage at the bed?. J. Glacioi., 29(101), 2847.Google Scholar
Johannessen, M. and Henriksen, Α. 1978. Chemistry of snow meltwater: changes in concentration during melting. Water Resour. Res., 14(4), 615619.Google Scholar
Jones, H. G., M. Tranter and Davies, T. D. 1989. The leaching of strong acid anions from snow during rain-on-snow events: evidence for two component mixing. International Association of Hydrological Sciences Publication 179 (Symposium at Baltimore 1987 – Atmospheric Deposition), 239250.Google Scholar
Lasaga, A. C. 1984. Chemical kinetics of water-rock interactions. J. Geophys. Res., 89(B6), 40094025.CrossRefGoogle Scholar
Lemmens, M. M. and Roger, M. 1978. Influence of ion-exchange on dissolved load of Alpine meltwaters. Earth Surface Processes., 3, 179187.Google Scholar
Lerman, A. 1979. Geochemical processes: water and sediment environments. New York, Wiley-Interscience.Google Scholar
Lliboutry, L. 1983. Modifications to the theory of intraglacial waterways for the case of subglacial ones. J. Glacioi., 29(102), 216226.Google Scholar
Lorrain, R. D. and Souchez, R. A. 1972. Sorption as a factor in the transport of major cations by meltwaters from an Alpine glacier. Qyat. Res., 2(2), 253256.Google Scholar
Luce, R. W., Barlett, R. W. and Parks, G. A. 1972. Dissolution kinetics of magnesium silicates. Geochim. Cosmochim. Acta., 36, 3550.Google Scholar
Moses, CO., Nordstrom, D.K. Herman, J.S. and Mills, A.L. 1987. Aqueous pyrite oxidation by dissolved oxygen and by ferric ions. Geochim. Cosmochim. Acta., 51, 15611571.Google Scholar
Oerter, H., Behrens, H. Hibsch, G. Rauert, W. and Stichler, W. 1980. Combined environmental isotope and electrical conductivity investigations at the runoff of Vernagtferner (Oetztal Alps). Mater. Glyatsiol. Issled. Khron. Obsuzhdeniya 39, 157161.Google Scholar
Parks, G.A. 1990. Surface energy and adsorption at mineral/water interfaces: an introduction. Reviews in Mineralogy., 23, 133175.Google Scholar
Raiswell, R. 1984. Chemical models of solute acquisition in glacial melt waters. J. Glacioi., 30(104), 4957.Google Scholar
Raiswell, R. and Thomas, A. G. 1984. Solute acquisition in glacial melt waters. I. Fjallsjökull (south-east Iceland): bulk melt waters with closed-system characteristics. J. Glacioi., 30(104), 3543.Google Scholar
Raiswell, R., Brimblecombe, P. Dent, D. and Liss, P.S. 1980. Environmental chemistry. London, Edward Arnold.Google Scholar
Reynolds, R. C. and Johnson, N. M. 1972. Chemical weathering in the temperate glacial environment of the Northern Cascade Mountains. Geochim. Cosmochim. Acta., 36, 537554.Google Scholar
Röthlisberger, H. 1972. Water pressure in intra-and subglacial channels. J. Glacioi., 11(62), 177203.Google Scholar
Sharp, M. 1991. Hydrological influences from meltwater quality data: the unfulfilled potential. Proceedings of tht BHS 3rd National Hydrological Symposium, Southampton, 5.15.8.Google Scholar
Singer, P. C. and Stumm, W. 1970. Acidic mine drainage: the rate-determining step. Science., 167, 11211123.Google Scholar
Souchez, R.A. and Lemmens, M.M. 1987. Solutes. In Gurnell, A.M. and Clarke, H.J. eds. Glacifluvial sediment transfer. New York, John Wiley, 285303.Google Scholar
Stumm, W. and Morgan, J. J. 1981. Aquatic chemistry. Second edition. New York, Wiley-Interscience.Google Scholar
Thomas, A. G. and Raiswell, R. 1984. Solute acquisition in glacial melt waters. II. Argentière (French Alps): bulk melt waters with open-system characteristics. J. Glcial., 30(104), 4448.Google Scholar
Toran, L. and Harris, R. F. 1989. Interpretation of sulfur and oxygen isotopes in biological and abiological sulfide oxidation. Geochim. Cosmochim. Acta., 53, 23412348.Google Scholar
Tranter, M. 1982. Controls on the chemical composition of Alpine glacial meltwaters. (Ph.D. thesis, University of East Anglia.)Google Scholar
Tranter, M. 1991. Controls on the composition of snowmelt. In Davies, M., Tranter, T. D., and Jones, H. G. eds. Seasonal snowpacks; processes of compositional change. Berlin, Springer-Verlag, 241271.CrossRefGoogle Scholar
Tranter, M. and Raiswell, R. 1991. The composition of the englacial and subglacial component in bulk meltwaters draining the Gornergletscher. J. Glacioi., 37(125), 5966.Google Scholar
Tranter, M., Raiswell, R. and Mills, R.A. 1989. Chemical weathering reactions in Alpine glacial meltwaters. In Miles, D. L., ed. Proceedings of the 6th International Symposium on Water-Rock Interaction, Malvern, 1989. Rotterdam, Balkema, 687690.Google Scholar
Turner, R. S. and 7 others. 1990. Watershed and lake processes affecting surface water acid-base chemistry. Acidic Deposition: State of Science and Technology. National Acid Precipitation Assessment Program Report 10. Washington, DC, Government Printing Office.Google Scholar
Walder, J. S. 1986. Hydraulics of subglacial cavities. J. Glacioi., 32(112), 439445.Google Scholar
Weertman, J. and Birchfield, G. E. 1983. Basal water film, basal water pressure, and velocity of traveling waves on glaciers. J. Glacial., 29(101), 2027.Google Scholar
Wollast, R. 1967. Kinetics of alteration of K-feldspar in buffered solutions at low temperature. Geochim. Cosmochim. Acta., 31, 635648.Google Scholar
Figure 0

Fig. 1. a. The relationship between the sum of dissolved carbon dioxide species (CO2(aq) + HCO3 + CO32-) and pH. The ability of the solution to dissolve CO2 increases with increasing pH because CO2 (aq) dissociates into HCO3 and CO32-. The diagram was constructed by assuming that the solution is in equilibrium with a gas phase p(CO2) of 10−3.5 atm, that the Henry’s law constant, KH, is 10−3.5 mol atm 1−l and that the first and second dissociation constants for dissolved CO2 are 10−6.58 and 10−10.62 mol 1−1, respectively. These constants are for T= 0° C (Garrels and Christ, 1965). b. The effect of pH on the speciation of dissolved CO2. The sum of the dissolved species remains constant at 24 μ,eq 1−1 i.e. no CO2 is allowed to diffuse into or out of solution, c. As the pH increases, CO2 (aq) is converted into HCO3 and the p(CO2) in the gas phase with which the solution is in equilibrium decreases (see Appendix). Equilibrium constants as in Figure 1a.

Figure 1

Fig. 2. The p(CO2) of a mixture of Na-feldspar and deionized water, as a function of time and rock: water ratio. These results were obtained by adding a known weight of crushed feldspar to 11 of deionized water in open contact with the atmosphere at Τ = 20° C (Tranter, 1982). Mote that most values are below atmospheric p(CO2) (= l(T−3.5 atm). When p(CO2) is declining, the rate of weathering is greater than the rate of CO2 influx into the solution. The converse argument holds when p(CO2) is rising.

Figure 2

Fig. 3. The (a) discharge, (b) suspended-sediment concentration, (c) p(CO2) and (d) sulphate concentration of meltwaters draining Haut Glacier d’Arolla during the 1989 ablation season. There is some evidence for a decrease in p(CO2) as the discharge and suspended-sediment concentrations increase (Brown, 1991).

Figure 3

Fig. 4. The mass fraction of and where units are equivalents), as a function of discharge for bulk meltwaters draining Haut Glacier d’Arolla throughout the 1989 ablation season.

Figure 4

Fig. 5. The association of p(CO2) with the mass fraction of SO2 for bulk meltwaters of Haut Glacier d’Arolla throughout the 1989 ablation season.

Figure 5

Table 1. Estimates of the composition of delayed flow