Thesis (M. App. Sc.)--Carleton University, 2004. / Includes bibliographical references (p. 176-185). Also available in electronic format on the Internet.

Heavy rainfall on the evening of January 9 and morning of January 10, 1983 triggered debris avalanches and debris torrents at Smith Creek, western Whatcom County, Washington, USA. Nine debris avalanches are back analyzed in detail. Conclusions are drawn concerning, 1) climatic controls on debris avalanches and debris torrents; 2) debris avalanche characteristics; 3) hillslope hydrology; 4) slope stability. Rainfall data show that the January 9-10, 1983 storm had a 71-year recurrence interval in the 12-hour duration, with less than 6-year recurrence intervals in 1, 2, and 3-hour durations. In contrast, rainfall during a torrent event on January 29-30, 1971 had recurrence intervals of less than 2 years in all durations, but snowmelt was a contributing factor. The types of debris torrents produced by these contrasting storms are discussed. Four distinct failure geometries are defined, based on avalanche descriptions: 1) wedges; 2) drainage depressions; 3) logging roads; 4) discontinuity surfaces. Three scour zones are also distinguished, based on slope segment types observed. To model storm water table levels a one-dimensional, vertical, transient, saturated-unsaturated finite difference infiltration program is linked to a kinematic wave equation. Rainfall duration and intensity, initial conditions, soil hydraulic conductivity, and soil depth are factors controlling vertical soil discharge rates. January, 1983 discharges are clearly distinguishable from comparison storm discharges at all avalanches. Kinematic wave results help differentiate Coulomb shear and washout type failures, and provide pore pressures for stability analyses. The modified Mohr-Coulomb strength equation is used to outline factors controlling debris avalanche initiation. The factors are: 1) slope angle; 2) soil depth; 3) soil density; 4) vegetative cover; 5) bedrock surface characteristics; 6) snow. These factors are quantitatively assessed. Infinite slope analyses show limiting slope angles of 29.7° for Group I vegetation, and 24.6° for Group III vegetation. Vegetative cover and soil depth are the two controlling factors that change significantly over the short term. A root cohesion parameter, Cr, is used to assess the shear strength provided by vegetation. Four vegetative covers are distinguished, three of which were logged between 1918 and 1950: Group I - relatively weak understory vegetation (Cr range: 1.6 -2.0 kPa); Group II - understory plus stunted trees (Cr range: 2.3 - 2.6 kPa); Group III - understory plus mixed, regenerating forest (Cr range: 2.6 - 3.0 kPa); Group IV - old-growth forest of higher root strength. / Science, Faculty of / Earth, Ocean and Atmospheric Sciences, Department of / Graduate

Les glissements de terrain se produisent dans toutes les chaînes de montagnes où la résistance de massifs rocheux est insuffisante pour contrer l’action de la gravité. Les terrains volcaniques sont particulièrement susceptibles de s’effondrer car les édifices sont composés des lithologies diverses et variées qui peuvent être fortement fracturées. En plus, la croissance rapide des édifices volcaniques favorise leur instabilité et leur effondrement. L’activité magmatique est un facteur additionnel responsable de la déformation des édifices, tandis que l’activité hydrothermale réduit la résistance des roches volcaniques. Pour ces raisons, l’évaluation des aléas liés à l’effondrement des édifices et à la formation des avalanches des débris volcaniques mérite une attention particulière. Les caractéristiques physiques des composants des avalanches des débris ont une influence directe sur la dynamique de ce type d'écoulement. Les dépôts des avalanches de débris présentent une morphologie de surface composée des nombreuses collines (hummocks), qui montrent fréquemment les séquences volcaniques initiales, ce qui suggère un mécanisme de mise en place proche de celui des glissements de terrain. Cependant, d’autres dépôts présentent des crêtes allongées (rides) dont le mécanisme de formation est encore méconnu. Le volcan Tutucapa (sud du Pérou) a été affecté récemment par deux avalanches de débris. La plus ancienne, « Azufre », est d’âge Holocène et résulte de l’effondrement d’un complexe des dômes et d’une séquence volcanique altérée (hydrothermalisée) sous-jacente. La deuxième avalanche, « Paipatja », a eu lieu il y a seulement 200-230 ans BP et est associée à une grande éruption explosive du Tutupaca. Les dépôts de cette avalanche présentent notamment de nombreuses rides. Les deux dépôts d’avalanche montrent deux unités différentes : une unité inférieure, caractérisée par la présence des blocs altérés (hydrothermalisés) provenant de l’édifice basal, tandis que l’unité supérieure est constituée par des blocs du complexe de dômes actifs. Le travail de terrain montre que les rides de l’avalanche « Paipatja » présentent une forte variation de granulométrie entre leur partie centrale (enrichie en blocs grossiers) et leurs parties latérales, ce qui suggère un processus de ségrégation granulaire. Des expériences analogiques montrent que des écoulements de mélanges de particules des différentes tailles subissent un processus de ségrégation et de digitation granulaire qui engendre des rides par jonction de levées statiques qui délimitent un chenal d’écoulement. Le processus de formation des rides est facilité par de faibles différence de taille des particules dans des mélanges bidisperses. Ces résultats suggèrent que les rides observées au Tutupaca résultent d’un écoulement granulaire. Les principales caractéristiques morphologiques des structures formées lors de ces expériences de laboratoire ont été comparées qualitativement avec les structures observées dans les dépôts du Tutupaca. Les structures observées au Tutupaca montrent que deux mécanismes de mise en place peuvent coexister dans les avalanches de débris volcaniques : le glissement de blocs plus ou moins cohérents, et l’écoulement semblable à celui d’un matériau granulaire. Cela dépend probablement de la nature des différents matériaux à la source des avalanches. Cette information doit être prise en compte pour l’évaluation des aléas liés aux avalanches des débris car des mécanismes d’écoulement différents peuvent induire des fortes variations de la distance parcourue par ces avalanches. / Landslides occur in all mountainous terrain, where the rock strength is unable to support topographic loading. Volcanic rocks are particularly landslide prone, as they mix strong and weak lithologies and are highly pre-fractured. Also, volcanoes themselves, are peculiar mountains, as they grow, thus creating their own topographic instability. Magmatic activity also deforms the edifice, and hydrothermal activity reduces strength. For all these reasons, volcanoes need close consideration for hazards, especially for the landslide-derived rock avalanches. The characteristics and properties of different debris avalanche components influence their behavior during motion. Deposits are generally hummocky, preserving original layering, which indicates a slide-type emplacement. However, some deposits have ridged morphology for which the formation mechanisms are not well understood. Two recent debris avalanches occurred at the Tutupaca volcano (S Peru). The first one, “Azufre” is Holocene and involved the collapse of active domes and underlying older hydrothermally altered rocks. The second debris avalanche, “Paipatja” occurred 200-230 y BP and is associated with a large explosive event and this deposit is ridged. The excellent conservation state of the deposits and surface structures allows a comprehensive analysis of the ridges. Both deposits have two contrasting units: a lower basal edifice-derived hydrothermally-rich subunit and an upper dome-derived block-rich unit. Detailed fieldwork has shown that Paipatja ridges have coarser core material and are finer in troughs, suggesting grain size segregation. Using analog experiments, the process that allow ridge formation are explored. We find that the mixtures undergo granular segregation and differential flow that create fingering that forms ridges by junction of static léeves defining a channel flow. Granular segregation and fingering are favored by small particle size contrast during bi-dispersed flow. The results suggest that the ridges observed at Tutupaca are product of a granular flow We extract the morphological characteristics of the deposits of granular flows generated in the laboratory and make a qualitative comparison with the Tutupaca deposits. The description of the different landslide and debris avalanche features at Tutupaca shows that two types of debris avalanche motion can occur in volcanic debris avalanches: the sliding of blocks more or less coherent and a flow similar to a granular material. This probably depends on source materials and the conditions of different parts of the initial landslide. Such information should be taken into account when estimating hazards at other volcanic landslide sites, as the different behaviors may result in different run outs.

This study is an examination of the rate at which organic debris and clastic sediment accumulate in a gully after it is scoured by a debris torrent. Of particular interest is the effect that a change in land use from old-growth to clear-cut conditions may have on these rates. This change should result in a reduction in the delivery of large organic debris (LOD), which is a major factor in sediment storage in gullies. It is hypothesized that this change in land use, and the subsequent reduction in the LOD supply, should result in a significant difference in debris recharge rates between old-growth and clear-cut gullies. Twenty-nine gullies in both land-treatment groups were sampled on the west coast of the Queen Charlotte Islands. Sampling procedures involved the estimation of the volume of LOD and sediment in storage (normalized by the gully surface area) and the determination of the time elapsed since the last debris torrent. These data were then used to estimate recharge rates(3h1)am’year of LOD, sediment, and total debris. Recharge rates of each material were compared between land-treatment groups using the nonparametric Mann-Whitney test. This test revealed that LOD has been delivered to old-growth gullies at a significantly higher rate relative to clear-cut gullies. There was no significant difference in sediment and total debris recharge rates between gullies in the two groups, but this outcome was partially a result of the small samples and the different debris recharge times in each data set. Graphical representations of the data permitted the identification of possible temporal trends in sediment and debris accumulation, which may be strengthened with larger data sets. Debris recharge rates have several applications. The estimate of sediment volume stored in a gully can be used in the construction of local sediment budgets, as one component of a watershed sediment cascade is quantified. The calculation of debris recharge rates will provide insight into the transfer rate of sediment from hillslopes to low order channels and to the storage capacity of the channels. Finally, debris recharge rates can be used to improve knowledge of the frequency-magnitude characteristics of debris torrents in an area.

This study is an examination of the rate at which organic debris and clastic sediment accumulate in a gully after it is scoured by a debris torrent. Of particular interest is the effect that a change in land use from old-growth to clear-cut conditions may have on these rates. This change should result in a reduction in the delivery of large organic debris (LOD), which is a major factor in sediment storage in gullies. It is hypothesized that this change in land use, and the subsequent reduction in the LOD supply, should result in a significant difference in debris recharge rates between old-growth and clear-cut gullies. Twenty-nine gullies in both land-treatment groups were sampled on the west coast of the Queen Charlotte Islands. Sampling procedures involved the estimation of the volume of LOD and sediment in storage (normalized by the gully surface area) and the determination of the time elapsed since the last debris torrent. These data were then used to estimate recharge rates(3h1)am’year of LOD, sediment, and total debris. Recharge rates of each material were compared between land-treatment groups using the nonparametric Mann-Whitney test. This test revealed that LOD has been delivered to old-growth gullies at a significantly higher rate relative to clear-cut gullies. There was no significant difference in sediment and total debris recharge rates between gullies in the two groups, but this outcome was partially a result of the small samples and the different debris recharge times in each data set. Graphical representations of the data permitted the identification of possible temporal trends in sediment and debris accumulation, which may be strengthened with larger data sets. Debris recharge rates have several applications. The estimate of sediment volume stored in a gully can be used in the construction of local sediment budgets, as one component of a watershed sediment cascade is quantified. The calculation of debris recharge rates will provide insight into the transfer rate of sediment from hillslopes to low order channels and to the storage capacity of the channels. Finally, debris recharge rates can be used to improve knowledge of the frequency-magnitude characteristics of debris torrents in an area. / Arts, Faculty of / Geography, Department of / Graduate

Located in the Columbia River Gorge, The Red Bluff Landslide (18.8 km2) is one of four large landslides that make up the Cascade Landslide Complex. In its current form, the Red Bluff Landslide is a post-Missoula Flood feature made up of two components: an active upper lobe (8.6 km2) that is translational, creeping to the south at 25 cm/yr and spreading laterally to the east at 6 cm/yr over a semi-fixed portion (10.2 km2) of the Red Bluff Landslide area that has been "smoothed" by Missoula Floods. The upper active lobe is the landslide debris accumulated since Missoula Flood time (~15,000 yr. BP). Five separate collapse events have been identified and rock failures along the main scarp headwalls continue. Two rock avalanches on the Red Bluff Landslide were mapped. The Old Greenleaf Basin Rock Avalanche is estimated to have occurred 100 to 150 years ago, represents the fifth collapse event on the Red Bluff Landslide, and covers an area of 200,000 m2. It has a volume of 4.2 million m3; its length is 748 m and has a width of 215 m. On January 3, 2008, the Greenleaf Basin Rock Avalanche occurred, flowing over the Old Greenleaf Basin Rock Avalanche, covering an area of 100,000 m2 and deposited a volume of about 375,000 m3. Its length is 730 m with an average depth of 1.22 m. It contributed approximately 0.058% of the total volume and 0.01% of the surface area to the active upper lobe portion of the Red Bluff Landslide. The Greenleaf Basin Rock Avalanche was determined to be insignificant in the movement of the active part of the Red Bluff Landslide during the winter of 2007-2008. The original Cascade Landslide Complex map (Wise, 1961) included the Mosley Lakes Landslide which has now been removed because it lacked the characteristics of a landslide like a scarp. The original complex (35.5 km2) has been renamed the "Greater Cascade Landslide Complex" (43.0 km2), with the addition of the adjacent Stevenson Slide and the elimination of the Mosley Lakes Landslide.

The heavy precipitation event of November 3-8, 2006 dropped over 60 cm of rain onto the bare southern slopes of Mount St. Helens and generated debris flows in eight of the sixteen drainages outside the 1980 debris avalanche zone. Debris flows occurred on the upper catchments of the Muddy River, Shoestring Glacier, Pine Creek, June Lake, Butte Camp Dome, Blue Lake, Sheep Creek, and South Fork Toutle River. Debris flows were clustered on the west and south-east sides of the mountain. Of the eight debris flows, three were initiated by landslides, while five were initiated by headward or channel erosion. Six debris flows were initiated in deposits mapped as Holocene volcaniclastic deposits, while two were in 1980 pyroclastics on andesite flows. The largest (~975,000 m2) and longest (~8,900 m) debris flow was initiated by landslides in the upper South Fork Toutle River Drainage. The average debris flow initiation zone elevation was 1,750 m, with clusters around 1,700 m and 2,000 m elevation. The lower cluster is associated with basins that host modern or historic glaciers, while the upper is possibly associated with recent pyroclastic deposits. Upper drainages with debris flows averaged 41% slopes steeper than 33 degrees, while those without debris flows averaged 34%. The upper basins with debris flows averaged 6% snow and ice cover, 21% consolidated bedrock, and 74% unconsolidated deposits. Basins without debris flows averaged 3% snow and ice cover, 27% bedrock, and 67% unconsolidated deposits. Drainages with debris flows averaged an 89% loss of glacier area between 1998 and 2009, while those without debris flows lost 68%. Further comparing glacier coverage during that period found that only five of ten glaciers still existed in 2009. On average, the glaciers had reduced in area by 67%, decreased in length by 36%, and retreated by an average of 471 m during that period. Basin attributes were measured or calculated in order to construct a predictive debris flow model based on that of Pirot (2010) using multiple logistic regression. The most significant factors were the percentage of slopes steeper than 33 degrees, unconsolidated deposits in the upper basin, and average annual rainfall. These factors predicted the 2006 debris flows with an accuracy of 94% in a debris flow susceptibility map for Mount St. Helens.

Les avalanches de débris, qui résultent du démantèlement des flancs des édifices volcaniques et montagneux, sont des écoulements granulaires rapides et dangereux dont le monteur est la gravité et qui présentent des distances de transport extrêmement importantes. La dynamique de leur mise en place et leurs mécanismes de transport permettant cette très grande mobilité sont des phénomènes qui demeurent encore mal compris. De nombreux modèles existent pour expliquer la grande mobilité des avalanches de débris et incluent des processus basés sur la lubrification ou la fluidification de la masse granulaire mais également sur le phénomène de désintégration dynamique des éléments. Cependant la grande majorité des modèles proposés souffre du manque d’observations de terrain et de quantification de l’évolution des matériaux au cours de leur transport au sein de la masse granulaire. Afin d’identifier les principaux mécanismes de transport des avalanches de débris, nous proposons dans ce travail une étude de terrain détaillée de dépôts d’avalanches de débris volcaniques qui résultent du démantèlement d’un volcan bouclier océanique, le Piton des Neiges (île de La Réunion, océan Indien). L’approche est couplée à un examen morphométrique (dimension fractale et circularité), exoscopique et granulométrique des particules présentes dans les dépôts. Elle est complétée par l’examen de la fabrique des dépôts basée sur l’anisotropie de la susceptibilité magnétique (ASM). Les données obtenues nous permettent de mettre en évidence une évolution de la dynamique de transport et de mise en place des dépôts d’avalanches de débris depuis les zones sources jusqu’aux domaines de dépôt distaux. On montre également que la désintégration dynamique et le gonflement dispersif qui l’accompagne opèrent tout au long du transport et à toutes les échelles au-dessus d’une limite inférieure de broyage à 500 μm. En dessous de cette limite, la réduction granulométrique résulte uniquement de processus d’attrition par friction entre les particules. La grande mobilité des avalanches de débris pourrait ainsi être expliquée par l’effet combiné de la libération d’énergie élastique par la désintégration dynamique des particules > 500μm et par une réduction de la friction interne à la matrice liée aux interactions dispersives des particules fines (< 500 μm). L’ensemble des données permettent également de préciser les directions de transport et l’ampleur des avalanches de débris liées aux déstabilisations du massif du Piton des Neiges. / Debris avalanches, resulting from flank collapses that shape volcanic and mountainous edifices are rapidand dangerous gravity-driven granular flows that travel long run out distances. The dynamic and the transport mechanisms behind this high mobility remain poorly understood. The numerous models proposed to explain this high mobility include processes based on lubrication or fluidisation of the granular mass of the flow body, but also the dynamic disintegration of the transported particles. To date,all these proposed mechanisms lack observational support and quantification of the state of the particles of the granular mass during the transport. To identify the main transport mechanisms, we propose here detailed field studies of volcanic debris avalanches deposits resulting of flank-collapse events on an oceanic shield volcano, the Piton des Neiges (La Réunion Island, Indian Ocean). This study has been combined with morphometric (fractal dimension and circularity), exoscopic and grain-size analyses. Moreover, the fabric of the deposits has been investigated by with the characteristics of the anisotropy of the magnetic susceptibility (ASM). From these data we highlight a proximal to distal evolution of the debris avalanches transport and emplacement dynamics. We demonstrate that syn-transport dynamic disintegration continuously operates with the distance from the source down to a grinding limit of 500μm. Below this limit, the particle size reduction exclusively results from the attrition of the particles by frictional interactions. Thus, the exceptional mobility of debris avalanches may be explained by thecombined effect of elastic energy release during the dynamic disintegration of the larger clasts (> 500μm) and frictional reduction within the matrix due to the dispersive interactions between the finer particles (< 500 μm). All these data also allow to specify the transport direction and the approximate size of the debris avalanches related to the Piton des Neiges destabilisations.

Debris flows, common disturbances in many mountainous areas, initially scour or bury stream habitats; however, debris flows deliver vast amounts of wood, boulders, and gravel that may ultimately form complex stream habitat to potentially support a diverse salmonid assemblage. The materials deposited by debris flows would otherwise be inaccessible to streams, and thus deposits may play an important role in creating and maintaining complex salmonid habitat over time. Despite the potential of deposits for increasing habitat complexity, most fish studies have focused on the destructive effects of debris flows and short-term recovery and re-colonization in scour zones. Debris-flows that occurred during the record-setting winter storms of 1996 in western Oregon, USA, provide an opportunity to study intermediate-term effects of debris-flow deposits on abundances and habitat for juvenile salmonids. In this setting, I surveyed salmonid abundance and habitat in three Oregon Coast Range streams that contained several debris-flow deposits from the 1996 storms. I explained fish abundance using hierarchical models, accounting for heterogeneous detection probabilities with repeated counts from multiple-pass snorkeling. The "best" hierarchical model of detection probability and abundance was selected (QAIC) from pool and snorkel-pass characteristics separately for juvenile coho salmon (Oncorhynchus kisutch), age 0+ trout, and age 1+ trout (Oncorhynchus spp.) in each stream. Adding distance to the nearest 1996 debris-flow deposit (DDF) produced a significant drop-in-deviance for four of nine "best" models, including at least one in each stream and for each species/age-class. In these four models, salmonid abundance decreased with increasing distance from deposit. As a potential explanation, several pool habitat characteristics were correlated (Spearman's rank) with DDF. Results varied across streams, but generally, percent of substrate as bedrock was lower and boulder density and percent substrate as gravel were higher closer to deposits. Although repeat counts are increasingly used in hierarchical modeling of heterogeneous detection probabilities and abundance for other wildlife species, studies of fish often rely on uncalibrated, single-pass snorkel counts. When exploring the value of repeat counts, I found that juvenile salmonid abundance decreased with increasing distance from debris-flow deposits in more multiple-pass hierarchical models that accounted for heterogeneous detection probabilities than for single-pass models that did not. Thus, modeling heterogeneous detection probabilities with repeated snorkel counts may be beneficial in other situations, addressing limitations of uncalibrated indices without relying on methods such as electrofishing, which may be difficult or impossible for remote study areas, longer surveys, or sensitive species. My findings suggest that debris-flow deposits may influence salmonid abundances after 15 years, and support management of debris flow-prone hillslopes and low-order channels to deliver elements of stream habitat complexity. / Graduation date: 2013

Debris flows, which occur in mountain settings worldwide, have been particularly damaging in the glaciated basins flanking the stratovolcanoes in the Cascade Range of the northwestern United States. This thesis contains two manuscripts that respectively investigate the (1) initiation processes of debris flows in these glaciated catchments, and (2) debris flow occurrence and its effect on valley bottoms over the last thousand years. In a 2006 storm, seven debris flows initiated from proglacial gullies of separate basins on the flanks of Mount Rainier. Gully heads at glacier termini and distributed collapse of gully walls imply that clear water was transformed to debris flow through progressive addition of sediment along gully lengths. In the first study, we analyze gully changes, reconstruct runoff conditions, and assess spatial distributions of debris flows to infer the processes and conditions necessary for debris flow initiation in glaciated catchments. Gully measurements suggest that sediment bulking requires steep gradients, abundant unstable material, and sufficient gully length. Reconstruction of runoff generated during the storm suggests that glaciers are important for generating the runoff necessary for debris flow initiation, particularly because infiltration capacities on glacial till covered surfaces well exceed measured rainfall rates. Runoff generation from glaciers and abundant loose debris at their termini explain why all debris flows in the storm initiated from proglacial areas. Proglacial areas that produced debris flows have steeper drainage networks with significantly higher elevations and lower drainage areas, suggesting that debris flows are associated with high elevation glaciers with relatively steep proglacial areas. This correlation reflects positive slope-elevation trends for the Mount Rainier volcano. An indirect effect of glacier change is thus the change in the distribution of ice-free slopes, which influence a basin’s debris flow potential. These findings have implications for projections of debris flow activity in basins experiencing glacier change. The second study uses a variety of dating techniques to reconstruct a chronology of debris flows in the Kautz Creek valley on the southwest flank of Mount Rainier (Washington). Dendrochronologic dating of growth disturbances combined with lichenometric techniques constrained five debris flow ages from 1712 to 1915 AD. We also estimated ages of three debris flows ranging in age from ca. 970 to 1661. Run-out distances served as a proxy for debris flow magnitude, and indicate that at least 11, 2, and 1 debris flow(s) have traveled at least 1, 3, and 5 km from the valley head, respectively since ca. 1650. Valley form reflects the frequency-magnitude relationship indicated by the chronology. In the upper, relatively steep valley, discrete debris flow snouts and secondary channels are abundant, suggesting a process of debris flow conveyance, channel plugging, and channel avulsion. The lower valley is characterized by relatively smooth surfaces, an absence of bouldery debris flow snouts, few secondary channels, and relatively old surface ages inferred from the presence of tephra layers. We infer that the lower valley is deposited on by relatively infrequent, large magnitude, low-yield strength debris flows like an event in 1947, which deposited wide, tabular lobes of debris outside of the main channel. Debris flows during the Little Ice Age (LIA) predominantly traveled no further than the upper valley. Stratigraphic evidence suggests that the main Kautz Creek channel was filled during the LIA, enhancing debris flow deposition on the valley surface and perhaps reducing run-out lengths. Diminished areas and gradients in front of glaciers during the LIA also likely contributed to decreased run-out lengths. These findings suggest that changes in debris flow source and depositional zones resulting from temperature and glacier cycles influence the magnitude and run-out distances of debris flows, and the dynamics of deposition in valley bottoms. / Graduation date: 2013