Hostname: page-component-7b7567c8c-2rsb6 Total loading time: 0 Render date: 2026-11-12T08:32:33.237Z Has data issue: false hasContentIssue false
  • English
  • Français

Application of the Gravity Flow Theory to the Percolation of Melt Water Through Firn

Published online by Cambridge University Press:  20 January 2017

W. Ambach,
M. Blumthaler and
P. Kirchlechner
W. Ambach
Affiliation:
Institut für Medizinische Physik, Universität Innsbruck, Müllerstrasse 44, A-6020 Innsbruck, Austria
M. Blumthaler
Affiliation:
Institut für Medizinische Physik, Universität Innsbruck, Müllerstrasse 44, A-6020 Innsbruck, Austria
P. Kirchlechner
Affiliation:
Institut für Medizinische Physik, Universität Innsbruck, Müllerstrasse 44, A-6020 Innsbruck, Austria
Rights & Permissions [Opens in a new window]

Abstract

Application of the gravity flow theory to the percolation of melt water through the firn in the accumulation area of a temperate glacier explains the occurrence of shock fronts in the melt-water flux. The time of propagation of a shock front moving from the surface through the entire firn was calculated under various assumptions. Various time input functions of melt-water flux at the surface with constant total input volumes yield only slight differences in the time of propagation of the shock front at greater depths. The dependence of the time of propagation of a shock front on the input volume, on snow parameters, and on the total thickness of the firn was calculated. An approximately linear relation was found to exist between the time of propagation of a shock front moving through the firn and the total thickness of the firn. The drainage of melt water from the firn after the summer ablation period is also quantitatively explained by the gravity flow theory. All results are in good agreement with experimental data.

Résumé

Résumé

L’application de la théorie des écoulements gravitaires à la percolation des eaux de fusion à travers le névé dans la zone d’accumulation d’un glacier tempéré, explique qu’il puisse se produire des à-coups dans le flux d’eau de fusion. Le temps de propagation d’un front d’écoulement depuis la surface jusqu’à la totalité du névéa été calculé sous différentes hypothèses. Différents types de modulation dans le temps des apports d’eau de fusion en surface pour un apport total constant ne provoquent que de faibles variations dans le temps de propagation du front d’écoulement aux plus grandes profondeurs. On a calculé comment le volume des apports, l’état de la neige et l’épaisseur totale du névé, influait sur le temps de propagation du front d’écoulement. On a trouvé une relation à peu près linéaire entre le temps de propagation du front d’écoulement à travers le névé, et l’épaisseur totale de celui-ci. La vidange de l’eau de fonte à partir du névé après la période d’ablation estivale est aussi expliquée quantitativement par la théorie de l’écoulement gravitaire. Tous les résultats sont en bon accord avec les données expérimentales.

Zusammenfassung

Zusammenfassung

Die Anwendung der gravity flow Theorie auf den Sickerfluss des Schmelzwassers durch den Firnkörper des Akkumulationsgebietes eines temperierten Gletschers erklärt das Auftreten von Schoekfronten im Schmelzwasserfluss. Die Laufzeit einer Schockfront von der Oberfläche durch den gesamten Firnkörper wurde unter verschiedenen Annahmen berechnet. Bei gleich grossem Schmelzwasservolumen als Input an der Oberfläche und verschiedenen zeitlichen Inputfunktionen ergeben sich in grösserer Tiefe nur geringe Unterschiede für die Laufzeit der Schockfront. Die Abhängigkeit der Laufzeit der Schockfront vom Inputvolumen, von Schneeparametern und von der Gesamtmächtigkeit des Firnkörpers wurde berechnet. Zwischen der Laufzeit der Schockfront durch den Firnkörper und der Gesamtmächtigkeit des Firns ergibt sich eine annähernd lineare Beziehung. Die Abnahme des Schmelzwasserflusses im Firnkörper am Ende der sommerlichen Ablationsperiode wurde ebenfalls durch die gravity flow Theorie quantitativ erklärt. Alle Ergebnisse sind mit experimentellen Daten in guter Überemstimmung.

Information

Type
Research Article
Information
Journal of Glaciology , Volume 27 , Issue 95 , 1981, pp. 67 - 75
Copyright
Copyright © International Glaciological Society 1981

List of symbols

Introduction

The gravity flow theory was applied by Reference Colbeck and DavidsonColbeck and Davidson (1973) and by Reference Denoth and DenothDenoth and others (1979[b]) for describing the percolation of melt water through homogeneous snow. In this theory it is assumed that percolation through snow is caused by gravity alone. Phenomena due to capillary forces are neglected (Reference ColbeckColbeck, 1971).

Gravity flow theory is based on Darcy’s law. It is furthermore assumed that the relative permeability of the water phase is related to the saturation by a power law.Footnote * The flux u then depends on the effective saturation S* as follows:

(1)

where α is a constant and k (m2) is the permeability at maximum saturation. The power n probably depends on the stage of metamorphosis (Reference Denoth, Denoth, Seidenbusch, Blumthaler, Kirchlechner, Colbeck and RayDenoth and others, 1979[a], Reference Denoth and Denoth[b]). Between the saturation S, the irreducible saturation Si, and the effective saturation S* , the relation is

This relation combines with the continuity equation to give :

(2)

where z is the depth, positive in direction of gravity, t the time, and Φ the porosity.

In the present paper, the formulation of gravity flow theory according to Reference Colbeck and DavidsonColbeck and Davidson (1973) and Reference Denoth and DenothDenoth and others (1979[b]), is applied to the percolation of melt water through firn in the accumulation area of a temperate glacier with non-uniform porosity. Thus a linear depth profile of porosity is assumed, inhomogeneities of firn are, however, neglected. Gravity flow theory is applied For the following analyses:

the propagation of the shock front for various histories of melt-water input at the surface, the influence of snow parameters on the propagation of the shock front, and the drainage of melt water from the firn after the summer ablation period.

2. Analysis of the shock front by means of the gravity flow theory

The formation of a shock front in the percolation of melt water may be explained qualitatively in the following simple manner: the percolation velocity of melt water increases with increasing flux. Larger fluxes that occur at a later time may thus catch up with smaller, earlier fluxes. The confluence may cause a sudden increase in the melt-water flux, which is called a shock front.

Figure 1 gives an example: two positive half-waves of a sinusoidal function each of 12 h duration (dashed curve) are assumed for the melt-water input at the surface. The resulting melt-water flux at a depth of 24 m is shown versus time (solid curve). After 77 h and 101 h, sudden increases of melt-water flux occur which represent the shock fronts. This variation of the melt-water flux, was calculated by means of the gravity flow theory assuming a linear depth profile of porosity. The analysis was made by means of the method of characteristics as applied to differential equation (2).

Fig. 1.

Calculated melt-water flux at 24 m firn depth (solid curve) and time function of the melt-water input flux at the surface (dashed curve). The calculation was based on the following assumptions: amplitude of the melt water input flux, umax = 1×10-6 m3/m2 s, mean values of snow parameters in accordance with Figure 2, curve B, linear porosity profile from Φ0 = 0.5 at the surface to Φ = 0.1 at the lower boundary of the firn and n = 3 for the power in Equation (1).

Fig. 2.

Dependence of the time of propagation of the shock front down to 24 m firn depth on the amplitude of the melt-water flux at the surface (positive part of a sine function with a duration of input of 12 h). The calculation was based on the following assumptions: Linear porosity profile from Φ0 = 0.5 at the surface to Φ = 0. 1 at the lower boundary of the firn and n = 3 for the power in Equation (I). Curve B holds for mean values of the snow parameters in the firn, curves A and c give possible limiting values.

2.1. Dependence of the propagation of the shock front on the input function

The propagation of the shock front has been calculated for the following variations with time of the melt-water input at the surface under otherwise constant conditions:

  • a sinusoidal input function (positive half-wave),

  • a semicircular input function, and

  • a triangular, non-symmetric input function.

It was furthermore assumed that the total daily volume of the melt water is constant irrespective of the input function. According to Figure 1, the shock front at a depth of 24 m is observed after about 77 h. The times at which the shock fronts occur with various input functions may differ by up to 3 h. The following values of snow parameters were used for the analysis:

The permeability at maximum saturation k was calculated following the method of Reference ShimizuShimizu (1970):

where ρi is the density of ice and ρw the density of water.

2.2. Dependence of the propagation of the shock front on the amount of melt-water percolation and on the snow structure

The time it takes to propagate a shock front through the entire firn depends on the amount of melt-water flux at the surface and on the snow parameters. First, the time of propagation was calculated for a total firn thickness of 24 m. This depth was chosen, because theoretical results should be directly comparable with experimental data.

Figure 2 shows the calculated relation between the time of propagation of the shock front and the assumed amplitude of the melt-water flux at the surface. The amplitude of the meltwater flux at the surface (umax) given on the vertical axis, corresponds to the amplitude of a positive half-wave of a sinusoidal input function having 12 h duration. The curves A and C are valid for possible limiting values of the snow parameters; curve B holds for mean values. For this analysis, the firn was assumed to have a total thickness of 24 m with a linear profile of porosity with depth. Snow parameters are the grain size d and irreducible saturation Si. The numerical data can be read directly from Figure 2, which shows, for example, that an amplitude of the melt-water input flux of umax = 1 X 10-6m3/m2s and the given limiting values of the snow parameters (curves A and C) may give between 55 and 93 h as the time of propagation for a shock front from the surface down to 24 m firn depth. The assumed meltwater flux corresponds to a net ablation rate of 2.4 g/cm2 d.

2.3. Dependence of the propagation of the shock front on the total firn thickness

Figure 3 gives the time of propagation of a shock front through the firn as a function of the total thickness of the firn. The curve expresses an approximately linear relation, which means that the mean velocity of propagation of the shock front is approximately independent of the total thickness of the firn. Figure 3, for example, shows that the propagation of a shock front from the surface down to the lower boundary of the firn, i.e. 20 m, takes 64 h. Hence the mean rate of propagation of a shock front in the firn is 0.31 m/h, and this value is approximately independent of the total thickness of the firn.

Fig. 3.

Propagation time of the shock front as dependent on the total thickness of the firn layer for mean values of snow parameters (Fig. 2, curve B) and for an amplitude of the melt-water flux at the surface umax = 1 x 10-6m3/m2 s. The porosity is assumed to decrease linearly from ϕ 0 = 0.5 at the surface to ϕ = 0.1 at the lower boundary of the firn. n = 3 was taken for the power in Equation (I).

The relation given in Figure 3 was calculated under the following assumptions:

Linear porosity profile having the values ϕ 0 = 0.5 at the surface and ϕ =0.1 at the lower boundary of the firn, independent of the total firn thickness.

Amplitude of the melt-water input flux at the surface umax = 1 X 10-6 m3 /m2 s having the shape of the positive half-wave of a sine function.

Mean snow parameters in accordance with curve B of Figure 2.

For the accumulation area of a temperate glacier with local differences in the total thickness of firn, the time lag between snow melt at the surface and inflow of the shock front into the water table at the lower boundary of the firn may be read for any value of total firn thickness from Figure 3.

3. Comparison between gravity flow theory and experimental data

In former experimental studies a time lag of 50–100 h was found between the maxima of snow melt at the surface resulting from weather conditions and the corresponding maxima of melt-water flux at the lower boundary of the firn at 24 m depth (Reference AmbachAmbach and others, 1978; H. Behrens and others, in press). These observed time lags, reviewed in Table I, are in good agreement with the time of propagation of a shock front through the firn calculated by means of the gravity flow theory. Figure 2 shows the time of propagation of the shock front to be 77 h for a total firn thickness of 24 m and for an amplitude of a melt-water flux at the surface of umax = 1 X 10-6 m3/m2 s, in accordance with a net ablation rate of 2.4 g/cm2 d, and for mean values of snow parameters (Fig, 2, curve B), A ±50% change in amplitude of the melt-water input at the surface causes the time of propagation to change from 60 to 120 h. These times of propagation correspond to mean velocities of propagation of the shock front between 0.4 and 0.2 m/h. The result compares with data given by other authors on the percolation velocities of melt water through the firn of a glacier: 0.1 m/h (Krimmel and others, 1973), 0.12 m/h (Reference Vallon, Vallon, Petit and FabreVallon and others, 1976), 0.22 m/h (Reference SharpSharp, [1952]) and 0.35 m/h (Reference AmbachAmbach and others, 1978). Furthermore, measurements in natural snow cover show similar results: 0.07–0.18 m/h (Reference FujinoFujino, 1968), 0.04–0.22 m/h (Reference KobayashiKobayashi, 1973), and 0.30–0.66 m/h (Reference WakahamaWakahama, 1968).Footnote *

Table I.Time lag Δt between the maximum snow melt at the surface and the maximum meltwater inflow qmax at 24 m firn depth in the accumulation area of Kesselwandferner in 1975 and 1977 (Ambach and others, 1978, fig. 2 and Behrens and others, in press, fig. 4). As the evaluation was based on daily mean values, the resolution for the time lag Δt is one day

During a period of strong ablation, diurnal fluctuations of the melt-water flux in 24 m firn depth have been observed in the form of shock fronts in the accumulation areas of Kesselwandferner (Reference AmbachAmbach and others, 1978, fig. 3d). These results are also in agreement with gravity flow theory. Figure 4 compares theoretical and experimental results. In these calculations, the amplitude of the melt-water input flux at the surface was again assumed to be umax =1x 10-6 m3/m2 s and again mean values of the snow parameters (Fig. 2, curve B) were used. Furthermore the melt-water input flux at the surface was assumed to be periodic over several days. The daily variations of the calculated melt-water flux and the occurrence of the shock front are in satisfactory agreement with the measured results. Here, however, it must be remembered that the maximum values of melt-water flux at a firn depth of 24 m correspond to the melt-water input at the surface with a time lag of 3 d. This displacement is due to the time of shock-front propagation.

Fig. 4.

Diurnal variations of the melt-water flux at 24 m depth calculated for mean values of the snow parameters (Fig. 2, curve B ) and for an amplitude of the melt-water flux at the surface of umax = 1 X 10 -6 m3/m2 s. The snow melt is assumed to be periodic over several days (positive half-wave of a sine function of 12 h duration ). Furthermore, a linear porosity profile was assumed from ϕ0 = 0.5 at th e surface to ϕ = O.1 at the lower boundary of th e firn and n = 3 for the power in Equation (1). The plotted dots are data obtained from continuous records (Reference AmbachAmbach and others, 1978, fig 3d).

Another feature confirming the applicability of the gravity flow theory to the calculation of melt-water flux in the accumulation area of an Alpine glacier is shown in Figure 5. It describes the drainage of melt water through the firn after the summer ablation period, being measured as inflow into the impermeable part of a 30 m deep pit (Reference AmbachAmbach and others, 1978, fig. 2, first half of October). The decrease in the melt-water inflow can be described by the relation:

(3)

The above relation follows from the gravity flow theory used so far (see Appendix). In Equation (3), q is the melt-water flow (expressed as volume flux in m3/s) at the lower boundary of the firn, t the time, and qo and to are constants. The parameter n (n > 1) is the power in the relation between the relative permeability and saturation (cf. Equation (1)).

Fig. 5.

Drainage of the melt-water flux from the firn of Kesselwandferner (Ötztal Alps) after the summer ablation Period 1975. The gravity flow theory (solid curve) was applied for the calculation. The data are daily means of the melt-water flux determined from continuous records (Reference AmbachAmbach and others, 1978, fig. 2, period: first half of October).

The parameters n and to were determined by a least-square calculation. The calculated melt-water inflow and the values measured are shown in Figure 5. The most important result is that by optimization the power n becomes n = 2.8. Experiments made by Reference Colbeck and DavidsonColbeck and Davidson ( 1973) on homogenized snow yielded a mean value of n = 3. Furthermore, extended experiments by Reference Denoth, Denoth, Seidenbusch, Blumthaler, Kirchlechner, Colbeck and RayDenoth and others (1979[a], Reference Denoth and Denoth[b]), on homogenized snow, showed n to vary between 1.4 and 4.6. It was found that this value probably depends on the snow structure, varying with the stage of metamorphosis.

Appendix Derivation of the relation for the drainage of melt water from firn after the ablation period

The conditions assumed are

  • (a) The differential equation for the melt-water flux in the snow is given by Equation (2)

  • (b) We assume a linear porosity profile in the firn

    with z = 0 at the surface and c the gradient of porosity.

  • (c) The permeability is calculated according to Reference ShimizuShimizu (1970)

For the assumed porosity profile the permeability also depends on z:

where ß = 6.0X 10-5d2(m2) and γ = 7.15.

For the characteristics of the differential equation we find:

After integration

where t = t0 for z = 0. Hence t0 is the starting time of a certain u value at the surface which is given by the input function. The starting time of the fluxes considered for drainage of the melt water from the firn after an ablation period lies immediately at the end of the last input. The starting time t0 is assumed to be constant for all u values, as the time of propagation of these u values are significantly longer than the differences between the start times. Thus the above equation is valid for all u values at constant 0 and for a given depth z:

The inflow into a reservoir is obtained by multiplying by the catchment area. Hence

where K and q0 are constants.

Footnotes

page note page 68 note * This relation is attributed in the hydrological literature to S. F. Aver’yanov in a paper in Russian entitled "Concerning the permeability of subsurface soils in the case of incomplete saturation", In zhenernyy Sbornik, Tom 7, 1950, but we have been unable to see the original paper. A briefer account appeared in Reference Aver’yanovAver’yanov (1949)

page note page 72 note * These measurements were carried out using tracers and the data give mean particle velocities v. From theoretical considerations, in general , where is the velocity of propagation of the shock front for a constant flux u.

References

Ambach, W., and others. 1978. Untersuchungen der Wassertafel am Kesselwandferner ( Ötztaler Alpen) an einem 30 Meter tiefen Firnschacht, von W. Ambach [and 9 others]. Zeitschrijl für Gletscherkunde und Glazialgeologie, Bd. 14, Ht. 1, p. 61–71 Google Scholar
Aver’yanov, S. F. 1949. Zavisimost’ vodopropitsayemosti pochvo-gruntov ot soderzhaniya v nikh vozdukha [Dependence of permeability of subsurface soils on their air content]. Doklady Akademii Nauk SSSR, Tom 69, No. 2, p. 141–44.Google Scholar
Behrens, H., and others. In press. Neue Ergebnisse zur Bewegung des Schmelzwassers im Firnkörper des Akkumulationsgebietes eines Alpengletschers (Kesselwandferner — Ötztaler Alpen), von H. Behrens [and 7 others], Zeitschrift für Gletscherkunde und Glazialgeologie.Google Scholar
Colbeck, S. C. 1971. One-dimensional water flow through snow. U.S. Cold Regions Research and Engineering Laboratory. Research Report 296.Google Scholar
Colbeck, S. C., and Davidson, G. 1973. Water percolation through homogeneous snow. (In [International Hydrological Decade.] The role of snow and ice in hydrology. Proceedings of the Banff symposia, September 1972. Paris, UNESCO; Geneva, WMO; Budapest, IAHS, Vol. 1, p. 242–57. ([Publication No. 107 de l’Association Internationale d’Hydrologie Scientifique.]))Google Scholar
Denoth, A., and others. 1979[a]. Some experimental data on water percolation through homogeneous snow, [by] Denoth, A., Seidenbusch, W., Blumthaler, M., Kirchlechner, P.. (In Colbeck, S. C., and Ray, M., ed. Proceedings of a meeting on modeling of snow cover runoff, 26-28 September 1978, [held in] Hanover, New Hampshire. Hanover, N.H., U.S. Cold Regions Research and Engineering Laboratory, p. 253–56.)Google Scholar
Denoth, A., and others. 1979[b]. Study of water drainage from columns of snow, [by] Denoth, A. [and 5 others]. U.S. Cold Regions Research and Engineering Laboratory. Report 79-1.Google Scholar
Fujino, K. 1968. Sekisetsu naibu de no yūsetsu-sui no ryūka sokudo no sokutei. I [Measurement of flowing down speed of melt water in a snow cover. I]. Teiott-kagaku: Low Temperature Science, Ser. A, [No.] 26, p. 87–100.Google Scholar
Kobayashi, D. 1973. Sekisetsuehū no yūsetsusui no himaku ryūka sokudo. I [Flow-down speed of melt water in snow cover. I]. Teion-kagaku: Low Temperature Science, Ser. A, [No.] 31, p. 135–42.Google Scholar
Krimmel, R. M., and others. 1973. Water flow through a temperate glacier, [by] Krimmel, R. M., Tangborn, W. V., and Meier, M. F.. (In [International Hydrological Decade.] The role of snow and ice in hydrology. Proceedings of the Banff symposia, September 1972. Paris, UNESCO; Geneva, WMO; Budapest, IAHS, p. 401–16. ([Publication No. 107 de 1’Association Internationale d’Hydrologie Scientifique.]))Google Scholar
Sharp, R. P. [1952.] Meltwater behavior in firn on upper Seward Glacier, St Elias Mountains, Canada. Union Géodésique et Géophysique Internationale. Association Internationale d’Hydrologic Scientifique. Assemblée générale de Bruxelles, 1951. Tom. 1, p. 246–53. (Publication No. 32 de l’Association Internationale d’Hydrologie.)Google Scholar
Shimizu, H. 1970. Air permeability of deposited snow. Contributions from the Institute of Low Temperature Science, Ser. A, No. 22.Google Scholar
Vallon, M., and others. 1976. Study of an ice core to the bedrock in the accumulation zone of an Alpine glacier, by Vallon, M., Petit, J.-R., and Fabre, B.. Journal of Glaciology, Vol. 17, No. 75, p. 13–27.Google Scholar
Wakahama, G. 1968. Sekisetsu-nai ni okeru yūsetsu-sui no idō. III. Suiro ryūka sokudo himaku ryūka sokudo [Infiltration of melt water into snow cover. III. Flowing down speed of melt water in a snow cover]. Teion- kagaku: Low Temperature Science, Ser. A, [No.] 26, p. 77–86.Google Scholar
Figure 0

Fig. 1.Calculated melt-water flux at 24 m firn depth (solid curve) and time function of the melt-water input flux at the surface (dashed curve). The calculation was based on the following assumptions: amplitude of the melt water input flux, umax = 1×10-6 m3/m2 s, mean values of snow parameters in accordance with Figure 2, curve B, linear porosity profile from Φ0 = 0.5 at the surface to Φ = 0.1 at the lower boundary of the firn and n = 3 for the power in Equation (1).

Figure 1

Fig. 2.Dependence of the time of propagation of the shock front down to 24 m firn depth on the amplitude of the melt-water flux at the surface (positive part of a sine function with a duration of input of 12 h). The calculation was based on the following assumptions: Linear porosity profile from Φ0 = 0.5 at the surface to Φ = 0. 1 at the lower boundary of the firn and n = 3 for the power in Equation (I). Curve B holds for mean values of the snow parameters in the firn, curves A and c give possible limiting values.

Figure 2

Fig. 3.Propagation time of the shock front as dependent on the total thickness of the firn layer for mean values of snow parameters (Fig. 2, curve B) and for an amplitude of the melt-water flux at the surface umax = 1 x 10-6m3/m2 s. The porosity is assumed to decrease linearly from ϕ 0 = 0.5 at the surface to ϕ = 0.1 at the lower boundary of the firn. n = 3 was taken for the power in Equation (I).

Figure 3

Table I.Time lag Δt between the maximum snow melt at the surface and the maximum meltwater inflow qmax at 24 m firn depth in the accumulation area of Kesselwandferner in 1975 and 1977 (Ambach and others, 1978, fig. 2 and Behrens and others, in press, fig. 4). As the evaluation was based on daily mean values, the resolution for the time lag Δt is one day

Figure 4

Fig. 4.Diurnal variations of the melt-water flux at 24 m depth calculated for mean values of the snow parameters (Fig. 2, curve B ) and for an amplitude of the melt-water flux at the surface of umax = 1 X 10 -6 m3/m2 s. The snow melt is assumed to be periodic over several days (positive half-wave of a sine function of 12 h duration ). Furthermore, a linear porosity profile was assumed from ϕ0 = 0.5 at th e surface to ϕ = O.1 at the lower boundary of th e firn and n = 3 for the power in Equation (1). The plotted dots are data obtained from continuous records (Ambach and others, 1978, fig 3d).

Figure 5

Fig. 5.Drainage of the melt-water flux from the firn of Kesselwandferner (Ötztal Alps) after the summer ablation Period 1975. The gravity flow theory (solid curve) was applied for the calculation. The data are daily means of the melt-water flux determined from continuous records (Ambach and others, 1978, fig. 2, period: first half of October).