Hostname: page-component-cd9895bd7-p9bg8 Total loading time: 0 Render date: 2024-12-14T13:42:01.643Z Has data issue: false hasContentIssue false

Frequent outburst floods from South Tahoma Glacier, Mount Rainier, U.S.A.: relation to debris flows, meteorological origin and implications for subglacial hydrology

Published online by Cambridge University Press:  20 January 2017

Joseph S. Walder
Affiliation:
U.S. Geological Survey, Cascades Volcano Observatory, 5400 MacArthur Boulevard, Vancouver, Washington 98661, U.S.A.
Carolyn L. Driedger
Affiliation:
U.S. Geological Survey, Cascades Volcano Observatory, 5400 MacArthur Boulevard, Vancouver, Washington 98661, U.S.A.
Rights & Permissions [Opens in a new window]

Abstract

Destructive debris flows occur frequently at glacierized Mount Rainier volcano, Washington, U.S.A. Twenty-three such flows have occurred in the Tahoma Creek valley since 1967. Hydrologic and geomorphic evidence indicate that all or nearly all of these flows began as outburst floods from South Tahoma Glacier. Flood waters are stored subglacially. The volume of stored water discharged during a typical outburst flood would form a layer several centimeters thick over the bed of the entire glacier, although it is more likely that large linked cavities account for most of the storage. Statistical analysis shows that outburst floods usually occur during periods of atypically hot or rainy weather in summer or early autumn, and that the probability of an outburst increases with temperature (a proxy measure of ablation rate) or rainfall rate. We suggest than outburst floods are triggered by rapid water input to the glacier bed. causing water-pressure transients that destabilize the linked-cavity) system, The correlation between outburst Hoods and meteorological factors casts doubt on an earlier hypothesis that melting around geothermal vents triggers outburst floods from South Tahoma Glacier.

Type
Research Article
Copyright
Copyright © International Glaciological Society 1995

Introduction

Destructive debris flows frequently move along stream valleys on glacierized Mount Rainier. Washington, the highest volcano in the coterminous United States. The smallest, but most frequent, of these debris Hows have generally been considered to originate as glacial outburst floods (Reference RichardsonRichardson, 1968; Reference CrandellCrandell, 1971: Reference Driedger and FountainDriedger and Fountain, 1989; Reference Scott, Pringle and VallanteScott and others, 1992). Transformation from water Hood to debris flow generally occurs in channels cut into stagnant ice and glacially derived sediment near glacier termini. In the last few decades, debris flows have been particularly frequent along Tahoma Creek, which drains South Tahoma Glaciers Fig. 1. Twenty-three debris flows moved along Tahoma Creek during the period 1967–92, including 15 in the years 1986–92, and repeatedly damaged a road and visitor facilities in Mount Rainier National Park.

Fig. 1. Index map of Mount Rainier showing location of South Tahoma Glacier and Tahoma Creek. Glaciers are drawn to scale, but streams, cultural features and weather stations (black dots) are shown only schematically. The Tahoma Creek picnic area was destroyed by debris flows in 1986 and 1987. In addition to the Tahoma Creek flows described in this paper, destructive debris flows triggered by outburst floods have also moved along Kautz Creek, Nisqually River and West Fork White River since the middle of the 20th century.

From 1988 to 1992. we served as consultants to the U.S. National Park Service, studying debris flows along Tahoma Creek, their relation to outburst floods, and associated geomorphic change in the Tahoma Creek drainage; assessing hazards presented by debris flows at Tahoma Creek; and examining whether those Hows might be predictable. A discussion of geomorphic evolution and of hazards is presented by Reference Walder and DriedgerWalder and Driedger (1994). In the present paper, we first summarize geomorphic and hydrologic data supporting an outburst-flood origin for most or all debris flows. We then present results of statistical analysis which indicate that debris flows (or outburst floods) are predictable in the sense that they are Strongly correlated with measured daily rainfall and maximum air temperature, data that may be considered proxies for the rate of water input to South Tahoma Glacier, Finally, we suggest a plausible hydrological model of the glacier than accounts for the link between weather and outburst floods.

Description of The Study Area

South Tahoma Glacier retreated and thinned substantially between the mid-19th century and about 1960 (Reference Sigafoos and HendricksSigafoos and Hendricks, 1972), resulting in a stagnant tongue of rock-rich ice. The glacier advanced slightly between about 1960 and 1975 (summarized in Driedger (1986)), with advancing ice overriding older, stagnant ice Fig. 2. Rapid retreat and thinning resumed about 1975, resulting in yet more rock-rich stagnant ice.

Fig. 2. Oblique aerial photograph (2 October 1984) showing South Tahoma Glacier prior to the outburst-flood phase than began in 1986. Stagnant-ice zones are indicated. Outburst floods in 1967 initiated incision of the reach of Tahoma Creek near center. (Photograph from U.S. Geological Survey archives in Tacoma. Washington.)

During the period 1988–92, outlet streams emerged from the active glacier ice and flowed over the upper stagnant-ice area (Fig. 2) before merging into a single stream. Downstream of that point, Tahoma Creek flows for 2 km at the bottom of a gorge deeply incised into stagnant ice, glacial till and fluvioglacial deposits. The nature and history of that incision, caused by outburst floods from South Tahoma Glacier, arc described in detail in Reference Walder and DriedgerWalder and Driedger (1994). Briefly, incision was initiated by a series of outburst floods in 1967 Fig. 2, and was increased by floods over the next decade. Outburst floods transform to debris flows as they pass through the incised reach. Debris-rich ice that formed the tongue of the active glacier as late as 1976 has itself been deeply incised since the renewal of outburst-flood activity in 1986.

Relation of Outburst Floods to Debris Flows

Most of the evidence for outburst floods from South Tahoma Glacier is circumstantial, owing largely to the glacier’s remoteness from the frequently visited parts of Mount Rainier National Park. A U.S. National Park Service employee observed one outburst flood in 1967 (cited by Reference CrandellCrandell (1971, p. 58)). There is also good evidence that, on at least a few occasions, flood waters emerged near an icefall and thence flowed over the glacier’s surface, eroding ice along the way (Reference CrandellCrandell, 1971; Reference Scott, Pringle and VallanteScott and others, 1992). Most of the evidence for outburst floods, however, has to do with their relation to debris flows. Because of this circumstance, it is necessary briefly to review pertinent data about debris flows at Tahoma Creek.

The debris flow of 26 July 1988 was serendipitously observed by U.S. Geological Survey hydrologist CH. Swift (written communication, 28 July 1988) in a reach where the flow was neither eroding nor depositing substantially and where it moved within a single, well-defined channel. For purposes of later discussion, we now summarize his observations. Swift visually estimated flow width, depth and velocity, as well as the relative fractions of water and sediment. A hydrograph based on his estimates is shown in Figure 3. The flow lasted about 80 min. Peak water discharge Q p was about 190 m3s−1, and a total water volume of approximately 3 × 105m3 passed before discharge returned to normal. Average water discharge during the debris flow was about 55 m3s−1 — about 15–20 times the base (low. These values are probably accurate to within a factor of 2.

Fig. 3. Hydrograph of the dry-weather debris flow of 26 July 1988 along Tahoma Creek. The estimated water flux within the debris flow is also shown. The abrupt rising limb is typical of a debris-flow hydrograph, although unlike that of a clear-water flood. Flow characteristics were visually estimated by C. H. Swift, U.S. Geological Survey, and are probably accurate to about a factor of 2.

Reference Walder and DriedgerWalder and Driedger (1994) considered whether sources of water other than outburst Hoods could account for debris flows along Tahoma Creek. Of the 23 debris flows since 1967, 15 (including the one seen by Swift) occurred during periods of sunny, dry weather, and it is plausible that at least these 15 flows began as outburst floods, because their estimated peak water discharge was much greater than typical stream (low in Tahoma Creek. Other potential water sources are (1) stream (low impounded by a slope failure and then released catastrophically; and (2) pore water held within stream-bank sediments that slumped and thence moved downstream as a debris flow. Reference Walder and DriedgerWalder and Driedger (1994) showed that neither of these sources would suffice. There is no geomorphic evidence of damming along Tahoma Creek, nor could a plausibly sized landslide dam impound the volume of water actually discharged in a debris flow of typical size at Tahoma Creek. There is also evidence (with one possible exception, discussed below) than slumping stream-bank sediments in the deeply incised reach of Tahoma Creek do not mobilize directly into debris flows, even during periods of heavy rainfall. Thus we conclude that dry-weather debris flows along Tahoma Creek originate as outburst floods.

Consider now the wet-weather debris flows. One of these — the debris flow of 16 October 1988 — evidently involved a flowslide of saturated morainal material Fig. 4. 16 October 1988 was one of the rainiest days on record in the period 1986–92 (U.S. National Oceanic and Atmospheric Administration, 1986–92), so it is conceivable that the flowslide was triggered simply by rainfall, without an outburst flood. However, between 28 September (the date of an earlier aerial photograph than Figure 4) and 28 October 1988, there had been substantial stream incision, bank failure, and stream avulsion upstream of the morainal slump, in a reach where the stream crossed stagnant ice. This observation at least suggests that an outburst flood was involved in the 16 October 1988 debris flow. Subsequently — once during each of the years 1989 through 1991 and twice in 1992 — debris flows moved along Tahoma Creek during rainy periods in late summer or autumn. There is no geomorphic evidence that these latter debris flows began as flowslides.

Fig. 4. Vertical aerial photograph taken on 28 October 1988 showing part of Tahoma Creek drainage. Upstream is toward the top. The circular feature near the center is the scar of a flowslide that transformed into the debris flow of 16 October 1988. The geometry of the scar is inconsistent with an origin of this feature as a rotational slump and strongly suggests that the sediment flowed in a liquified state into the gorge. Note that scars of other slumps along Tahoma Creek are quite different in shape. (Photograph by Washington State Department of Transportation.)

Wet-weather debris flows of 1989–92 were comparable to the flow observed by Swift Fig. 3 in their effects and in the extent of their deposits, so we assume, as an admittedly rough estimate, that Q p ≈ 200 m3 s−1 for each of these flows. For comparison, Reference RichardsonRichardson (1968), on the basis of flood-frequency relationships suggested by Reference Bodhaine and ThomasBodhaine and Thomas (1964), estimated Q p ≈ 11 m3 s−1 for the mean annual flood of meteorological origin in Tahoma Creek (at a point about 5 km downstream of Glacier Island); and Q p ≈ 26 m3 s−1 for a flood with a 2% annual exceedence probability the so-called 50 year flood). This order-of-magnitude discrepancy, along with the lack of evidence for a flowslide origin for these debris flows, leads us to conclude that the wet-weather debris flows of 1989–92 probably began as outburst floods, although part of the water in these debris flows must comprise pore water in the entrained sediment. In the rest of this paper, we therefore assume that all debris flows at Tahoma Creek, with the possible exception of the 16 October 1988 event, began as outburst floods.

Relation of Outburst Floods to Weather

All outburst floods from South Tahoma Glacier since 1967 have occurred during either summer or early autumn, as noted also by Reference Driedger and FountainDriedger and Fountain (1989) and Reference Scott, Pringle and VallanteScott and others (1992). We did statistical analysis of meteorological data collected near South Tahoma Glacier to test whether outburst Hoods have been correlated with atypical weather conditions, and whether meteorological data might be useful predictors of outburst-flood occurrence.

The U.S. National Weather Service maintains two stations near South Tahoma Glacier:one at the Paradise Ranger Station, about 9 km cast-southeast of Glacier Island, at an elevation of about 1650 m and near Nisqually Glacier; the other at Longmire, 9.5 km south-east of Glacier Island at an elevation of about 840 m Fig. 1. Statistical analysis shows that published values (U.S. National Oceanic and Atmospheric Administration, 1986–92) of daily rainfall and maximum daily temperature at these two stations are well correlated. For the period under consideration, we found from linear regression:

(1)

(2)

where T max is the maximum daily Celsius temperature, R is daily precipitation, and r 2 is the coefficient of determination,

During the summer and autumn of 1989 and 1990. we tried to measure temperature and rainfall at an elevation of about 2050 m on Glacier Island, near South Tahoma Glacier (Figs 1 and 2). but data collection was frequently interrupted by equipment failures and we had to abandon the expensive effort to maintain this remote site. For the 1989 data, linear regression gives the result:

(3)

We suspect [he temperature sensor at Glacier Island might have been inadequately protected from the effects of solar radiation, and that this problem bears on the relatively low r 2 value. The rain gauge malfunctioned, so we did not analyze precipitation data, The 1990 data were too discontinuous for meaningful analysis.

All in all, available data suggest that temperature and precipitation on the south side of Mount Rainier are spatially coherent. In what follows, we will assume that weather data from Paradise may be used as surrogates for conditions near South Tahoma Glacier; We restrict our analysis to weather data for the period 1986–92, because flood dates are not accurately known for earlier years. The analysis uses data for the months May through November, approximately encompassing the typical ablation season for glaciers on Mount Rainier. Actual outburst-flood dates ranged from 29 June to 9 November.

The statistical distributions of rainfall and maximum daily temperature at Paradise arc shown in Figures 5 and 6. (Days with measurable snowfall were excluded in compiling Figure 6.) The solid curve in each figure is the cumulative distribution function (henceforth CDF); the crosses connected by dashed lines show the distributions for days on which outburst floods from South Tahoma Glacier occurred. We now pose the question:What is the probability that the distributions of meteorological data for days on which outburst floods occurred are drawn aí random from the CDFs? We use the Kolmogorov–Smirnov test of goodness of lit Reference ConoverConover, 1971) to find an answer. The measured test statistic d m is the largest vertical distance between the CDF and the stepwise sample curve connecting values measured on days of outburst-flood occurrence (Figs 5 and 6). For a given number n of data points in the sample population, one may calculate the probability Ρ that for another sample of n points the largest vertical distance between the curves would exceed the measured value d m. Ρ may also be interpreted as the probability that the sample population (in our case, days with outburst floods) is drawn at random from the parent population (all snow-free days from May through November). The Kolmogorov-Smirnov test is preferred to the chi-square test for small sample sizes, as in the present case, because it docs not assume a Gaussian distribution of the test statistic (Reference ConoverConover, 1971).

Fig. 5. Cumulative distribution function of values of daily rainfall at the Paradise Ranger Station, for the months May–November in the years 1986–92. Solid curve:all snow-free days. Crosses connected by dashed line:all days with debris flows. The Kolmogorov–Smirnov statistic dm is indicated.

Fig. 6. Cumulative distribution function of values of maximum daily temperature at Paradise Ranger Station, for the months May–November in the years 1986–92. Solid curve:all days. Crosses connected by dashed line:all days with debris flows. The Kolmogorov-Smirnov statistic dm is indicated.

Results are shown in the top halves of Tables 1 and 2 for the cases illustrated by Figures 5 and 6, as well as for analogous tests using 2, 3, and 4 day moving averages of maximum daily temperature and rainfall. (For moving averages, the last day of the averaging window is always the day of the flood.) We conclude that values of rainfall and maximum daily temperature for days with outburst floods are unlikely to be drawn from the parent population, although this conclusion weakens as we consider progressively longer averaging periods.

We can make stronger statements once we note (Figs 5 and 6) that there are two distinct sub-populations in the data for days with outburst floods:warm, dry days (T max > 20 °C, R = 0) and cool rainy days (T max < 15 °C, R > 0). We therefore examined the statistical distribution of rainfall excluding dry days, and that of maximum temperature excluding rainy days. Results are presented in the bottom halves of Tables 1 and 2, for tests using daily-data and moving averages, as before. Values of Ρ are in some cases less than 0.001.

Table 1. Results of Kolmogorov-Smirnov (R S) test for rainfall data from Paradise Ranger Station. Meteorological data used in the analysis are for the months May-November in the years 1986-92. One outburst flood for which the exact date is uncertain was excluded from the analysis. The K—S statistic is defined in the text and illustrated in Figure 5

Table 2. Table 2. Results of Kolmogorov–Smirnov (Κ–S) test for maximum-temperature data from Paradise Ranger Station. Meteorological data used in the analysis are for the months May–November in the years 1986–92. One outburst flood for which the exact dale is uncertain was excluded from the analysis. The Κ–S statistic is defined in the text and illustrated in Figure 6

On the basis of these tests of goodness of fit, we conclude that outburst floods from South Tahoma Glacier arc highly correlated with atypical weather, and in particular, with relatively rainy or hot weather Reference Driedger and FountainDriedger and Fountain (1989) and Reference KambScott and others (1992) reached similar conclusions, but did not support their statements with statistical tests.

The conditional probability P c of an outburst flood from South Tahoma Glacier, as a function of R and T max during the years 1986–92, for the two sub-populations noted above, is shown in Figures 7 and 8. The conditional probability is the relative frequency of flood occurrence given a particular value of R or T max. For example, the probability of an outburst flood on any day with ≥40 mm of rain was about 0.22. There are fairly clear relations between P c and either R or T max measured on the day of a flood; for progressively longer averaging periods, the trends are not so obvious, owing at least partly to the fact that averaging “compresses” the range of R and T max. We caution against using Figures 7 and 8 to predict outburst-flood probabilities, however, because the glacier, including, presumably, its drainage system, is changing with time, making the statistical problem non-stationary.

Fig. 7. Conditional probability of an outburst flood from South Tahoma Glacier as a function of daily rainfall at Paradise Ranger Station, for days with measurable rainfall during the months May-November in the years 1986–92. The last day of the averaging period is always the day of the flood. The dashed Hue is an approximate fit to the data.

Fig. 8. Conditional probability of an outburst flood from South Tahoma Glacier as a function of maximum daily temperature at Paradise Ranger Station, for dry days during the months May–November in the years 1986–92. The last day of the averaging period is always the day of the flood. The dashed line is an approximate fit to the data.

If R and T max measured at Paradise are considered rough surrogates for rate of water input to South Tahoma Glacier, the conditional probabilities shown in Figures 7 and 8 are reasonably consistent with each other. To see this, consider that a plausible value for the ablation rate at South Tahoma Glacier on the hottest summer days (when T max ≈ 30 °C at Paradise) is perhaps 40 mm d−1 (cf. Reference FountainFountain, 1992). From Figures 7 and 8, we find (extrapolating from the dashed curves) Pc ≈ 0.13 for T max ≈ 30 °C, P c ≈ 0.22 for R = 40 mm. In light of the surrogate nature of the meteorological data, the two probabilities are reasonably close. We therefore suggest the following physical interpretation, elaborated in the next section: outburst floods from South Tahoma Glacier are meteorologically driven, with the probability of an outburst flood determined primarily by the rate of water input to the glacier.

Water Storage at the Glacier Bed and an Hypothesis for the Origin of Outburst Floods

Outburst floods from South Tahoma Glacier involve release of water stored either englacially or subglacially. For the 26 July 1988 event Fig. 3, which was probably typical, the total water volume (3 × 105 m3), if spread evenly over the glacier bed (area about 2 km2), would be equivalent to an average thickness of about 150 mm of water. Undoubtedly, part of the water in the debris flow was derived from sediments added to the flood waters as the) moved downstream, but even if such added water, accounted for 90% of the debris flow’s total water content, glacially derived water would still represent average basal storage of 15 mm. The most plausible way to store such a copious amount of water at the bed of South Tahoma Glacier is in basal cavities Fig. 9. Reference Driedger and FountainDriedger and Fountain (1989) previously reached the same conclusion by a similar line of reasoning. For comparison, the water volume released from subglacial storage during surge “slow-down” events at Variegated Glacier in 1983 represented average basal storage of about 50–100 mm (Reference Humphrey, Raymond and HarrisonHumphrey and others, 1986; Reference HumphreyHumphrey, 1987).

Fig. 9. Schematic drawing of the basal linked-cavity network. Areas of ice/bed contact are shaded, areas of ice/bed separation are blank. Ice separates from the bed along the hachured lines and re-contacts the bed along the plain lines. (After Reference KambR. Kamb (1987).)

South Tahoma Glacier rests on very steep slopes (averaging 20–30°) and probably (lows over large ledges eroded from the lavas and volcaniclastic deposits that constitute most of Mount Rainier. Indeed, one such ledge, several tens of meters high and formerly beneath an icefall, has been exhumed as the ice has thinned during the last decade. Supposing a subglacial cavity existed in the lee of such a ledge, with a cross-sectional area along the ice-flow direction of 50 m2 — corresponding to a cavity 10 m in height and length, with a triangular cross-section — a total cavity length of 6 km would be required to store 3 × 105m3 of water. Even if the stored-water volume were an order of magnitude less — corresponding to the case in which most debris-flow water comes from the sediments themselves a 600 m length of these enormous cavities would be required. Yet it is hard to imagine how the requisite volumes of water could be stored if giant cavities did not exist. If all cavities were decimetric in scale (such as those mapped by Reference Walder and HalletWalder and Hallet (1979) and Reference Hallet and AndersonHallet and Anderson (1982)), hundreds of kilometers of linked cavities would be required to hold the stored volume of water.

Except early in the ablation season, surficial meltwater or rainfall will enter the englacial drainage system rapidly and thence reach the glacier bed. Water at the base of an alpine glacier usually collects in one or more main channels incised into the ice (Reference ShreveShreve, 1972; Reference RöthlisbergerRöthlisberger, 1972), but at South Tahoma Glacier the route from englacial conduit to trunk channel is likely to involve passage through large linked cavities. Flow along the linked-cavity network Fig. 9 is regulated by the throttling effect of constrictions, or orifices, where the ice closely approaches the bed (Reference WalderWalder, 1986; Reference KambKamb, 1987). Based on these concepts and on results in the previous section, we propose the following hypothesis:outburst floods from South Tahoma Glacier occur when the rate of water input to the subglacial drainage system causes subglacial water pressure to rise rapidly, destabilizing the orifices in the linked-cavity network, in line with Kamb’s (1987) analysis. The linked-cavity drainage system then degenerates into a multiple-tunnel system that rapidly drains stored water.

In Kamb’s (1987) analysis, orifice stability is controlled by the dimensionless parameter Ξ, defined by

(4)

where α is the average hydraulic gradient, approximated as the ice-surface slope; is the length of orifices divided by the length of cavities; ω is the tortuosity of the water-flow path; D is a constant equal to 31 km; M is the Manning roughness of the orifices; η is the effective viscosity of the ice; ν is sliding velocity; σ is effective pressure (ice-overburden pressure minus water pressure); and h is the height of the orifice. For Ξ greater than about 1, orifices will enlarge unstably in response to transient water-pressure increases of any magnitude. For progressively smaller values of Ξ, the orifices are stable against progressively larger water-pressure transients.

To estimate a value of Ξ for the basal-cavity system of South Tahoma Glacier, we adopt the parameter values α = 0.3, = 10, ω = 4, M = 0.1 m−1/3s. η = 0.1 bar a, ν = 10 m a−1 and σ = 5 bar. The values for ∧, ω, M and η follow from Kamb, and should be reasonable for the South Tahoma Glacier system. The value of α. applies for the active ice and is based on a U.S. Geological Survey topographic map. We have no direct measurements of sliding velocity or effective pressure at the glacier bed. Our choice of ν is intuitive but plausible, based on the fact that South Tahoma Glacier has been thinning and retreating for nearly 2 decades, and our judgment that the glacier is sluggish. (For purposes of comparison, the average sliding speed of the nearby, larger Nisqually Glacier was about 50 m a −1 in 1969, while the glacier was advancing Reference HodgeHodge, 1974).) The value for σ is approximately equal to the estimated average ice-overburden pressure (based on glacier area-and-volume estimates by Reference DriedgerDriedger and Kennard (1986)) and thus a rough upper bound on effective pressure. With these parameter values, we find Ξ = 3.8h 7/6, when h is given in meters. We would then require h < 0.3 m for the orifices to have even marginal stability. If meter-scale cavities exist, as we have previously suggested, even the orifices may have decimetric dimensions, and thus be unstable in the event of transient increases in water pressure.

Indirect evidence possibly bearing on the size of orifices conies from a bizarre, macabre source. In December 1946, a U.S. Marine Corps airplane Crashed near the head of South Tahoma Glacier, with no survivors. The wreckage was quickly buried under snow and was never recovered. Deposits of the 1986 and 1987 debris flows contained many mangled bits of the aircraft, including propeller blades and pieces of sheet metal up to about 0.5 m on a side (personal communication from K.M. Scott, 1993). A plausible interpretation is that the airplane wreckage had been advected downward to the glacier bed. had melted out and accumulated in cavities, and had finally been Hushed out as cavities had drained through R channels during outburst floods. The size of the metal fragments would then set a lower bound on typical orifice dimensions. Other interpretations are possible, of course.

We conclude that Kamb’s model of orifice instability yields quantitatively reasonable results for the case in hand. Moreover, in the context of Kamb’s model, it is not surprising that the probability of an outburst flood increases as the rate of water delivery to the glacier bed increases (Figs 7 and 8). The greater the rate of water delivery, the larger must be the pressure transients, and thus the greater the likelihood that (for given Ξ) orifices in the cavity network will grow unstably.

The shortest period of time between two outburst floods from South Tahoma Glacier was about 56 h — the interval between the dry-weather floods of 29 and 31 August 1967 (Reference CrandellCrandell, 1971, p. 58). The dry-weather floods of 28 and 31 August 1987 were separated by only about 75 h. In the context of the cavity model discussed above, two (at least interpretations are possible. The linked-cavity system may rapidly reform itself after disruption leading to an outburst flood; alternatively, any individual outburst floods may represent break-down and drainage of just part of the cavity network. We favor the second interpretation, because the volume of water released in the second of each pair of floods could not plausibly have been generated by ablation in the period between floods, even if all surface-derived meltwater had been stored subglacially. The second interpretation is also in line with observations elsewhere (e.g. Reference FountainFountain, 1992). that alpine glaciers may contain several distinct internal drainage basins.

A puzzling aspect of outburst floods at South Tahoma Glacier is their probabilistic nature: few occurrences of hot weather or heavy rain do in fact coincide with outburst floods. One interpretation may lie that the glacier’s drainage system varies significantly with time. For example, the geometry of orifices in the cavity system might change with time, thereby affecting the value of Ξ. Another interpretation is that we simply have insufficient information about the rate of water input lo the drainage system through time, as well as about the geometry of the drainage system itself, to construct a deterministic model of outburst-flood release.

Discussion

Many other investigators have reported cases of glaciers that repeatedly spawn outburst floods, but these have involved drainage of either subaerial ice-dammed lakes (e.g. Reference ClarkeClarke, 1982) or subglacial lakes formed by geothermal activity (Reference BjörnssonBjörnsson. 1992). South Tahoma Glacier appears to be different. Wé are skeptical of Reference CrandellCrandell’s (1971) conclusion that geothermal activity causes outburst floods from South Tahoma Glacier. Crandell’s argument rested on anecdotal evidence, such as observations of supposed steam plumes, gathered from a variety of observers. We believe it would be a remarkable coincidence if the statistical relationship between outburst floods and meteorological variables existed independently of a causal relationship between those variables and outburst floods, even though we understand the pertinent physics inexactly. Furthermore, there is no glacier-surface relief indicative of melting over a geothermal area (cf. Reference Björnsson and EinarssonBjörnsson and Einarsson, 1991).

Other glaciers on Mount Rainier have released Outburst floods from subglacial storage, though none with nearly the frequency of South Tahoma Glacier (Reference Driedger and FountainDriedger and Fountain, 1989). It may be that the bed geometry of South Tahoma Glacier is particularly conducive to the formation of large subglacial cavities, bed geometry must be related to the structure and credibility of the bedrock, factors which on a strato-volcano like Mount Rainier are very heterogeneous. We note that the present phase of outburst-flood activity at South Tahoma Glacier coincides with a period of substantial ice retreat and thinning, whereas the previous phase of activity (1967 to the mid-1970s) was during a period of ice advance, with the glacier considerably thicker than at present. This suggests that the propensity for outburst floods is not related simply to ice thickness.

Conclusions

South Tahoma Glacier, located on Mount Rainier volcano, has released aí least 22 outburst floods since 1967, including at least 14 in the years 1986–92. The volume of water released during a typical outburst flood corresponds to a basal water layer several centimeters thick, a physically impossible situation; so water must be stored in large subglacial cavities. Geologic factors such as the volcano’s heterogeneous structure promote the existence of large basal cavities.

The likelihood of an outburst flood increases as either the temperature or rainfall rate increases. A reasonable interpretation is that flood probability depends on the rate of water input to the glacier bed. We therefore propose that the physical mechanism of flood release is the breakdown of the linked-cavity network over part of the glacier bed and rapid drainage through a multiple-tunnel system, much as envisaged by Reference KambKamb 1987 in bis discussion of glacier-surge termination.

Acknowledgements

The U.S. National Park Service provided partial funding for this study. L. Mastin introduced us to use of Kolmogorov-Smirnov statistics. We benefited from reviews of earlier versions of the paper by A. Fountain, N. Humphrey. R. Alley. R. LeB. Hooke and two anonymous referees.

References

Björnsson, H. 1992. Jökulhlaups in Iceland: prediction, characteristics and simulation. Ann. Glacial. 16. 95106.Google Scholar
Björnsson, H. and Einarsson, P. 1991. Volcanoes beneath Vainajökull, Iceland:evidence from radio echo-sounding, earthquakes and jökulhlaups. Jökull. 40. 1990, 147168.CrossRefGoogle Scholar
Bodhaine, G. L. and Thomas, D. M. 1964. Magnitude and frequency of floods in the United States. Part 12. Pacific Slope basins in Washington and upper Columbia River basin. U.S. Geol. Surv. Water-Supply Pap. 1687.Google Scholar
Clarke, G.K.C. 1982. Glacier outburst floods from “Hazard Lake”, Yukon Territory, and the problem of flood magnitude prediction, J. Glaciol., 28(98), 321.CrossRefGoogle Scholar
Conover, W.J. 1971. Practical nonparametric statistics. New York. John Wiley and Sons.Google Scholar
Crandell, D. R. 1971. Postglacial lahars from Mount Rainier volcano, Washington. U.S. Geol. Surv. Prof. Pap. 677.Google Scholar
Driedger, C. L. 1986. A visitor’s guide to Mount Rainier glaciers. Longmire, WA, Pacific Northwest National Parks and Forests Association.Google Scholar
Driedger, C. L. and Fountain, A. G. 1989. Glacier outburst floods at Mount Rainier. Washington State, U.S.A. Ann. Glaciol., 13, 5155.CrossRefGoogle Scholar
Driedger, G. I. and Kennard, P.M. 1986. Ire volumes on Cascade volcanoes:Mount Rainier. Mount flood. Three Sisters, and Mount Shasta. U.S. Geol. Surv. Prof. Pap. 1365.Google Scholar
Fountain, A.G. 1992. Subglacial water flow inferred from stream measurements at South Cascade Glacier. Washington. U.S.A. J. Glaciol., 38(128). 5164.Google Scholar
Hallet, B. and Anderson, R.S. 1982. Detailed glacial geomorphology of a proglacial bedrock area at Castleguard Glacier, Alberta. Canada. Z. Glescherkd. Glazialgeol., 16(2), 1980, 171184.Google Scholar
Hodge, S M. 1974. Variations in the sliding of a temperate glacier. J. Glaciol. 13(69), 349369.Google Scholar
Humphrey, N. F. 1987. Basal hydrology of a surge-type glacier: observations and theory relating to Variegated Glacier. (Ph.D. thesis. University of Washington.)Google Scholar
Humphrey, N. Raymond, C. and Harrison, W. 1986. Discharges of turbid water during mini-surges of Variegated Glacier, Alaska. U.S.A. J. Glacial., 32(111), 195207.Google Scholar
Kamb, B. 1987. Glacier surge mechanism based on linked cavity configuration of the basal water conduit system. J. Geophys. Res., 92(B9), 90839100.Google Scholar
Richardson, D. 1968. Glacier outburst floods in the Pacific Northwest. D.S. Geol. Surv. Prof. Pap. 600-D, D79D86.Google Scholar
Röthlisberger, H. 1972. Water pressure in intra- and subglacial channels. J. Glacial., 11(62). 177203.Google Scholar
Scott, K. M., Pringle, P.T. and Vallante, J.W. 1992. Sedimentology, behavior and hazards of debris flows at Mount Rainier. Washington. U.S. Geol. Surv. Open-file Rep. 90385.Google Scholar
Shreve, R.L. 1972. Movement of water in glaciers. J. Glaciol., 11(62). 205214.Google Scholar
Sigafoos, R.S. and Hendricks, E. L. 1972. Recent activity of glaciers of Mount Rainier, Washington. Botanical evidence of glacier activitv. U.S. Geol. Surv: Prof. Pap. 387-B.Google Scholar
U.S. National Oceanic and Atmospheric Administration. 1986-92. Climatological data — Washington. Asheville. NC, U.S. National Oceanographic and Atmospheric Administration, 90-96 (511) [in each volume].Google Scholar
Walder, J. S. 1986. Hydraulics of subglacial cavities. J. Glaciol., 32(112), 439445.Google Scholar
Walder, J. S. and Driedger, C. L. 1994. Geomorphic change caused by glacial outburst floods and debris flows at Mount Rainier. Washington, with emphasis on Tahoma Creek valley. U.S. Geol. Surv: Water-Resources Investigation Rep. 934093.Google Scholar
Walder, J. and Hallet, B. 1979. Geometry of former subglacial water channels and cavities. J. Glaciol., 23(89), 335346.Google Scholar
Figure 0

Fig. 1. Index map of Mount Rainier showing location of South Tahoma Glacier and Tahoma Creek. Glaciers are drawn to scale, but streams, cultural features and weather stations (black dots) are shown only schematically. The Tahoma Creek picnic area was destroyed by debris flows in 1986 and 1987. In addition to the Tahoma Creek flows described in this paper, destructive debris flows triggered by outburst floods have also moved along Kautz Creek, Nisqually River and West Fork White River since the middle of the 20th century.

Figure 1

Fig. 2. Oblique aerial photograph (2 October 1984) showing South Tahoma Glacier prior to the outburst-flood phase than began in 1986. Stagnant-ice zones are indicated. Outburst floods in 1967 initiated incision of the reach of Tahoma Creek near center. (Photograph from U.S. Geological Survey archives in Tacoma. Washington.)

Figure 2

Fig. 3. Hydrograph of the dry-weather debris flow of 26 July 1988 along Tahoma Creek. The estimated water flux within the debris flow is also shown. The abrupt rising limb is typical of a debris-flow hydrograph, although unlike that of a clear-water flood. Flow characteristics were visually estimated by C. H. Swift, U.S. Geological Survey, and are probably accurate to about a factor of 2.

Figure 3

Fig. 4. Vertical aerial photograph taken on 28 October 1988 showing part of Tahoma Creek drainage. Upstream is toward the top. The circular feature near the center is the scar of a flowslide that transformed into the debris flow of 16 October 1988. The geometry of the scar is inconsistent with an origin of this feature as a rotational slump and strongly suggests that the sediment flowed in a liquified state into the gorge. Note that scars of other slumps along Tahoma Creek are quite different in shape. (Photograph by Washington State Department of Transportation.)

Figure 4

Fig. 5. Cumulative distribution function of values of daily rainfall at the Paradise Ranger Station, for the months May–November in the years 1986–92. Solid curve:all snow-free days. Crosses connected by dashed line:all days with debris flows. The Kolmogorov–Smirnov statistic dm is indicated.

Figure 5

Fig. 6. Cumulative distribution function of values of maximum daily temperature at Paradise Ranger Station, for the months May–November in the years 1986–92. Solid curve:all days. Crosses connected by dashed line:all days with debris flows. The Kolmogorov-Smirnov statistic dm is indicated.

Figure 6

Table 1. Results of Kolmogorov-Smirnov (R S) test for rainfall data from Paradise Ranger Station. Meteorological data used in the analysis are for the months May-November in the years 1986-92. One outburst flood for which the exact date is uncertain was excluded from the analysis. The K—S statistic is defined in the text and illustrated in Figure 5

Figure 7

Table 2. Table 2. Results of Kolmogorov–Smirnov (Κ–S) test for maximum-temperature data from Paradise Ranger Station. Meteorological data used in the analysis are for the months May–November in the years 1986–92. One outburst flood for which the exact dale is uncertain was excluded from the analysis. The Κ–S statistic is defined in the text and illustrated in Figure 6

Figure 8

Fig. 7. Conditional probability of an outburst flood from South Tahoma Glacier as a function of daily rainfall at Paradise Ranger Station, for days with measurable rainfall during the months May-November in the years 1986–92. The last day of the averaging period is always the day of the flood. The dashed Hue is an approximate fit to the data.

Figure 9

Fig. 8. Conditional probability of an outburst flood from South Tahoma Glacier as a function of maximum daily temperature at Paradise Ranger Station, for dry days during the months May–November in the years 1986–92. The last day of the averaging period is always the day of the flood. The dashed line is an approximate fit to the data.

Figure 10

Fig. 9. Schematic drawing of the basal linked-cavity network. Areas of ice/bed contact are shaded, areas of ice/bed separation are blank. Ice separates from the bed along the hachured lines and re-contacts the bed along the plain lines. (After R. Kamb (1987).)