Introduction
With the realisation that there may be a strong coupling between the hydrology of warm-based glaciers and their flow mechanics (Reference IkenIken, 1981), there has arisen a strong interest in determining the morphology of subglacial drainage systems and its evolution over time (Reference KambKamb and others, 1983; Reference Seaberg, Seaberg, Hooke and WibergSesberg and others, 1988; Reference Willis, Sharp and RichardsWillis and others, 1990; Reference FountainFountain, 1993; Reference Hock and HookeHock and Hooke, 1993; personal communication from P. Nienow, I. Willis, M. Sharp and K. Richards, 1994). Tracer techniques have been widely used for this purpose, and are perhaps the best method available for determining drainage-system characteristics on a glacier-wide scale (Reference HookeHooke, 1989).
Tracing involves the injection of a known quantity of tracer (salt, or more commonly a fluroescent dye such as rhodamine (B or WT) or fluorescein) at a known point (commonly a moulin or crevasse on the glacier surface, but occasionally the base of a borehole drilled to the glacier bed), and its detection in one or more meltwater streams draining from the glacier, Three forms of detection are widely practised:
-
1. Activated charcoal detectors are placed in the stream and absorb dye as it passes. After completion of a test, the detector is rinsed in an eluent which is then passed through a fluorometer. The presence or absence of dye is used to determine whether or not dye emerged in a particular stream (e.g. Reference Smart and BrownSmart and Brown, 1973).
-
2. Discrete samples of stream water are collected and run through a fluorometer to measure dye concentration (e.g. Reference Seaberg, Seaberg, Hooke and WibergSeaberg and others, 1988).
-
3. Water drawn from the stream is passed continuously through a fluorometer, and dye concentrations are recorded on a data logger or chart recorder (e.g. Reference FountainFountain, 1993).
Methods 2 and 3 allow the Construction of dye breakthrough curves (plots of dye concentration vs time).
Quantitative analyses of dye breakthrough curves form the basis of modern dye-tracing techniques (Reference SmartSmart, 1988). However, the success of such analyses is strongly dependent upon the accuracy with which the breakthrough curve is resolved (Reference FischerFischer, 1968; Reference Behrens, Bergmann, Moser, Ambach and JochumBehrens and others, 1975). This is a function of the accuracy of the fluorometric detection procedures employed and the frequency with which water samples are collected for analysis of dye concentration. In this paper we consider the influence of the sampling interval on the determination and interpretation of the properties of dye breakthrough curves. Our goal is to determine the minimum sampling frequency which allows consistent determination of these properties.
Methods Used in The Analysis of Dye Breakthrough Curves
In this section, we identify the parameters most commonly employed in quantitative analyses of dye breakthrough curves, in order to show how an accurate determination of the form of the breakthrough curve is essential to the accurate estimation of these parameters.
Establishing the time between dye injection and peak concentration in an outflow stream (tm) allows estimation of the flow velocity, u, through a glacial drainage system (distance travelled/tm ). This provides information concerning the hydraulic structure of the flow system. Rapid flow velocities (>0.2 ms−1) have been interpreted as indicating flow through hydraulically efficient channelised drainage systems (Reference StenborgStenborg, 1969; Reference Krimmel, Tangborn and MeierKrimmel and others, 1973; Reference Behrens, Bergmann, Moser, Ambach and JochumBehrens and others, 1975; Reference Moser and AmbachMoser and Ambach, 1977; Reference Lang, Leibundgut and FestelLang and others, 1979; Reference BurkimsherBurkimsher, 1983). Lower flow velocities have usually been interpreted as indicating flow through a hydraulically inefficient “distributed” hydrological system such as a linked-cavity configuration, a thin basal water film or a permeable basal sediment layer (Reference Theakstone and KnudsenTheakstone and Knudsen, 1981; Reference Iken and BindschadlerIken and Bindschadler, 1986; Reference KambKamb, 1987; Reference Willis, Sharp and RichardsWillis and others, 1990). In some recent work, velocities of less than 0.1 ms−1 have been interpreted as indicating flow through highly sinuous, hydraulically inefficient, anabranching channels (Reference FountainFountain, 1993; Reference Hock and HookeHock and Hooke, 1993).
The shape of a dye breakthrough Curve provides information about the flow conditions within the system traced, particularly with regard to the mechanisms which cause dispersion of the dye cloud. An increasingly dispersed dye cloud can result from two main effects: (i) an increase in the range of flow velocities within the drainage path (which will enhance mechanical dispersion); and (ii) an increase in the number of storage locations or immobile zones along the drainage path (which will enhance molecular diffusion). In general, dispersion resulting from variations in flow velocity is the dominant factor in channelised systems, whilst dispersion resulting from storage retardation becomes increasingly significant in distributed-flow systems. Two parameters are widely used to describe the degree of dispersion of a breakthrough curve: the dispersion coefficient, D, (Reference TaylorTaylor, 1954; Reference Levenspiel and BischoffLevenspiel and Bischoff, 1963) and the dispersivity, d (Reference FischerFischer, 1968; Reference Behrens, Bergmann, Moser, Ambach and JochumBehrens and others, 1975).
The dispersion coefficient, D, is used to describe the rate of dispersion of a dye cloud during its passage through a drainage system. It has units of m2s−1. There are numerous methods of obtaining D from a dye break-through curve. Within glaciology, a detailed analysis of the processes of dispersion was undertaken by Reference BrugmanBrugman (1986). One of the methods used by her (and subsequently by Reference Seaberg, Seaberg, Hooke and WibergSeaberg and others (1988) and Reference Willis, Sharp and RichardsWillis and others (1990)) to derive a value of D is used here:
Here, x is the straight-line distance between the points of dye injection and recovery, tm is the time to maximum dye concentration (but is not taken from field data, but instead is treated as a variable to be obtained in the solution), and ti represents t1 and t2 . which are the times when the dye concentration reaches half its peak value on the rise to and decline from the concentration peak. Equation (1) represents two equations (for i = 1, 2) and is solved iteratively for tm until a value is chosen such that the two equations are satisfied and D is obtained. Unfortunately, because D is dependent upon both the mean flow velocity and the range of flow velocities within the drainage system it does not, on its own. provide a reliable insight into drainage-system structure. It is therefore usual to look also at the relationship between D and the mean flow velocity determined from a tracer test.
In a given channel wiih a constant roughness and hydraulic radius, it has been found that D is approximately proportional to the average flow Velocity, u (Reference FischerFischer, 1968; Reference Behrens, Bergmann, Moser, Ambach and JochumBehrens and others, 1975; Reference Matoszewski, Hotzl and WernerMaloszewski, 1992), such that:
Here, d is referred to as the dispersivity and has units of m. Reference Hock and HookeHock and Hooke (1993, p. 544) state that “of the parameters derived from analysis of dye-return curves, the dispersivity is the one that most clearly illustrates evolution of the drainage system” at Storglaciären, Sweden, during the summer 1989 melt season.
The interpretation of dispersivity values derived from tracer tests in glacier drainage systems is discussed more fully elsewhere (personal communication from P. Nienow, I. Willis, M. Sharp and K. Richards, 1994). For our present purposes, it suffices to say that the dispersivity represents the rate of spreading of a dye cloud relative to the rate of advection of the dye during transit through a flow system, and that it provides a measure of a length scale characteristic of the drainage system (Reference FischerFischer, 1968; Reference Gillham, Sudicky, Cherry and FrindGillham and others, 1984). An accurate estimate of d can only be obtained if D and u are well known. To achieve this, it is necessary to have a well-defined breakthrough curve so that the values of tm and ti, (Equation (1)) obtained from the curve and used in the calculation of u and D are correct or have an acceptable error margin. The importance of this was acknowledged by Reference FischerFischer (1968) and Reference Behrens, Bergmann, Moser, Ambach and JochumBehrens and others (1975, p. 381), the latter noting that, in certain tracing experiments at the Hintereisferner, Austria, the sampling interval was “comparatively long as compared to the time of passage of a tracer cloud”, and “the amount of data was therefore too small for a dispersion analysis”. At the Hintereisferner, a 5 min sampling interval was used whilst the passage of a dye cloud often took no more than 15–20min. For detailed dispersion studies, a shorler sampling interval was required to generate a more accurate breakthrough curve and thus more accurate values of D and u. We now seek to determine how short this sampling interval must be.
Sampling-Rate Effects on The Determination of Values for The Flow Velocity, Dispersion Coefficient and Dispersivity
In this section, we evaluate the sensitivity of values of u, D and d to variations in the sampling interval used to define a dye breakthrough curve. To do so, we make use of breakthrough curves resulting from injections made into moulins on Haut Glacier d’Arolla, Valais, Switzerland during the summer of 1990. The tracer used was rhodamine B, and dye detection was by continuous flow fluorometry using a Turner Designs 10–005 field fluorometer. Full details of the methods employed are given by Reference SharpSharp and others (1993).
To illustrate the dependence of values of u, D and d on the sampling interval, we consider a series of breakthrough Curves which resulted from ten injections conducted at a moulin located 1.44 km from the glacier snout on 14 August (Fig. 1a). These curves were generated by 1 min averaging of dye concentrations measured every 10 s using a Campbell Scientific Limited 21 X data logger. The curves were sub-sampled to produce breakthrough curves based on dye concentrations taken at 5 and 10 min intervals (Fig. 1b and c). Values of time to peak dye concentration (tm), u, D and d obtained from each set of curves are given in Table 1.
The shape of the breakthrough curve generated by a given injection becomes increasingly angular as the sampling interval is increased, whilst the detail of the curve is lost (Fig. 1). The general shape of a breakthrough curve remains relatively consistent as the sampling interval is increased, but the peakedness of individual curves tends is decrease as the short-lived dye-concentration peak is missed. In the case of injection 3, the breakthrough curve resulting from the 10 min sampling interval is very different from that produced with shorter sampling intervals.
Although there are no systematic variations in the values of u obtained from breakthrough curves constructed with different sampling intervals, it is clear that the values obtained with 5 and 10 min sampling intervals rarely coincide with those obtained with a 1 min sampling interval (Fig. 2a; Table 1). The mean value of D obtained from 1 min sampling was 2.81 m2s−1 (standard deviation (S.D.) = 0.70). Such consistent results are to be expected from a series of closely spaced injections made into a single moulin over a period of time that is short enough to assume a “constant” drainage configuration. Sampling intervals of 5 and 10 min produced mean D values of 4.15 (S.D. = 2.04) and 9.54 (S.D. = 7.06) m2s−1, respectively. These larger and more variable values (Fig. 2b), which might be taken to suggest significant changes in the drainage-system morphology during the course of a single day, clearly result from inadequate definition of the breakthrough curve rather than from actual changes in the drainage configuration.
Since D tends to increase with the sampling interval, whilst u shows no systematic variation, it follows that d also tends to increase as the sampling interval increases. The dispersivity values determined from the 1, 5 and 10 min data sets by regressing D on u (forcing the regression to omit a constant) were, respectively, 5.28 (r2 = 0.57, p < 0.005), 7.93 (r2 = 0.03) and 17.92 (r2 = 0.04) (Fig. 2c). These results suggest that a 5 min sampling interval is too great to obtain an accurate estimate of d from the breakthrough curves resulting from the series of injections conduted on 14 August 1990.
In order to determine die precise sampling interval required to generate accurate estimates of D and thus d, the curves were further sub-sampled at 2, 3 and 4 min intervals. Figure 2d shows that certain values of D increase markedly from the 1 min sampling-interval estimates when the sampling interval is increased from 3 to 4 min. This fact is further illustrated by the standard deviation derived from the ten estimates of D which is 0.70 for 1 min sampling, and increases from 0.82 to 1.33 for the 3 and 4 min sampling intervals, respectively (Table 2). Similarly, the correlation coefficient obtained from regressing D on u to obtain d decreases from 0.43 (p < 0.025) to 0.23 between the 3 and 4 min sampling intervals (Table 2).
Thus, in the above example, in which the dye-return times were all less than 62 min and the passage of over 90% of the recovered tracer took no more than 20 min, a sampling interval of 3 min or less was required to prevent erroneous values of d being derived from the breakthrough curves. However, the sampling interval required to produce a breakthrough curve that is defined well enough to allow accurate estimation of d will increase as the time of passage of a tracer cloud and/or the the distance travelled by the tracer increases. This is a simple reflection of the rate of change of concentrations during a dye return, which tends to decrease as the tracer cloud becomes more dispersed.
To illustrate this, we consider how varying the sampling interval affects estimates of u, D and d derived from an additional six breakthrough curves resulting from tracer tests conducted during the 1990 melt season. The time to peak dye concentration in these experiments varied between 81 and 267 min (Table 3). Figure 3 shows how the value of D derived from the breakthrough curve varies with sampling interval for each of the six injections.
For each curve, the point at which D begins to increase with sampling interval is taken to indicate the maximum possible sampling interval Which will allow acceptable estimates of D and d to be made. The ratio of the spread of the return curve (characterised by the time between half-peak concentration on the rising and falling limbs of the curve) to maximum allowable sampling interval varies from 1.7 to 6.1 (Table 3) with a mean of 3.4 (S.D. = 1.7). The ratio of time to peak dye concentration to maximum allowable sampling interval varies from 7.3 to 16.0 (Table 3) with a mean of 11.2 (S.D. = 3.5). Given the variability of both sets of results, tm, is considered a more appropriate parameter for establishing a suitable sampling interval due to the ease with which it can be determined in the field. Thus, if a tracer experiment is conducted with the intention of analyzing the shape of the breakthrough curve, it is advisable to adopt a sampling interval no grealer than tm /16. Greater sampling intervals are likely to result in inconsistent and inaccurate estimates of D, and hence of d.
Conclusions
The principal impact of increasing the sampling interval used to construct a dye breakthrough curve is to increase estimated values of the dispersion coefficient, D. Estimates of through flow velocity, u, are changed by increasing the sampling interval, but not in any consistent way. Both D and u are used in calculate the dispersivity, d, so estimates of d also tend to increase with sampling interval. Since increases in d have commonly been taken to indicate an increase in the complexity of glacial drainage systems, and in the number of storage locations within them, it is important to be certain that such increases are not merely an artefact of the sampling procedure Used to construct dye breakthrough curves. An analysis of the variation of D with sampling interval for six breakthrough curves which produced widely varying times to peak dye concentration (tm ) suggests that accurate and consistent estimates of d will only be produced when the sampling interval used is <tm /l6. Where this criterion Cannol be met, caution should be applied when making inferences about the character of glacial drainage systems based on an analysis of the dispersion characteristics of dye breakthrough curves.
Acknowledgements
This work was supported by the U.K. Natural Environment Research Council through grant GR3/7004a, and by grants from Earthwatch, P. Nienow acknowledges receipt of an NERC studentship (GT4/89/AAPS/53) and Fellowship (GT5/93/AAPS/1). We thank W.H. Theakstone (University of Manchester) for the loan of his fluorometer, Grande Dixence SA for logistic support, and Y. Bams for all her help and support in Arolla. Field assistance was provided by K. Richards, C. Hill, B. Hubbard, N. Arnold, W. Lawson, J. Branson, J. Strike, J.-L. Tison, S. Lane, M. Skidmore, N. Spedding and 36 Earthwatch volunteers. We should also like to thank two anonymous reviewers for helpful comments on an earlier version of this paper.